Molecular basic data calculation for radiation transport in chromatin

Sub-keV corrections to binary encounter cross section models for electron ionization of liquid water with application to the Geant4-DNA Monte Carlo code

INTRODUCTION

The electron ionization cross section of water is one of the most important input in Monte Carlo studies of cellular radiobiological effects. Analytical cross section models of the binary-encounter type have the potential of reducing simulation time and facilitate application to a variety of biological materials (other than water). The Binary-Encounter-Bethe (BEB) and Binary-Encounter-Dipole (BED) models of NIST are perhaps the most popular of such models giving reliable results for atoms and molecules in the gas-phase over a wide energy range. However, the use of such models to sub-keV electron energies in liquid water raises concerns due to the neglect of condensed phase effects that leads to a significant overestimation when compared to medium-specific dielectric models.

PURPOSE

To modify the BEB and BED models towards better agreement with the recommended low-energy dielectric model of Geant4-DNA (Option 4). To implement the new modifications to the existing BEB model of the Option 6 physics constructor of Geant4-DNA and re-evaluate fundamental transport quantities for sub-keV electrons.

METHODS

In analogy to a Yukawa potential a simple, yet physically-motivated, modification of the Burgess correction term is proposed to account for the reduction of the Coulomb interaction due to the polarizability of the target. The magnitude of the correction is guided by the dielectric-based ionization cross section implemented in Option 4.

RESULTS

Differential, total and stopping ionization cross sections for low-energy electrons in liquid water are presented. When combined with the Vriens correction (which is not included in Option 6), the proposed modification to the BEB and BED models brings the ionization and stopping cross sections in much better agreement against those used in the Option 4 dielectric model of Geant4-DNA, with up to 30% and 10% deviation, respectively. Implementation of the new correction to the Option 6 constructor of Geant4-DNA and re-evaluation of fundamental transport quantities, such as electron penetration ranges and dose-point-kernels, reduced the discrepancies from Option 4 at sub-keV energies from 20 to 100% (or more) to well below 10% in most cases.

CONCLUSIONS

A simple modification to the BEB and BED analytic models was found to improve their performance for sub-keV electrons in liquid water medium. Implementation of the new modification to the Option 6 constructor of Geant4-DNA significantly improved the agreement with the recommended low-energy Option 4 constructor for a variety of fundamental quantities related to electron transport.

Electron Impact Ionization of Adenine: Partial Cross Sections

Electron ionization of a genetically important nucleobase, adenine, was investigated from threshold to 500 eV using crossed electron beam–effusive molecular beam geometry and time-of-flight mass spectrometry. We measured the complete set of absolute partial cross sections for adenine using the relative flow technique (RFT) up to an electron energy of 500 eV. Normalization to absolute values was performed using electron ionization cross sections for argon and the vapor pressure data of adenine. The total cross sections obtained by summing the partial cross sections were compared with the existing theoretical and experimental data. The appearance energies of various fragment ions were also measured and compared with the reported data. The prominence of ions with mass (HCN)n+ (n = 1 to 5) indicated a possible pathway to form adenine in the interstellar medium through aggregation of HCN units. Analysis of the partial cross sections for various groups of fragment ions as a function of electron energy was found to give insights into their composition.

Review of the Geant4-DNA Simulation Toolkit for Radiobiological Applications at the Cellular and DNA Level

Simple Summary

A brief description of the methodologies to simulate ionizing radiation transport in biologically relevant matter is presented. Emphasis is given to the physical, chemical, and biological models of Geant4-DNA that enable mechanistic radiobiological modeling at the cellular and DNA level, important to improve the efficacy of existing and novel radiotherapeutic modalities for the treatment of cancer.

Abstract

The Geant4-DNA low energy extension of the Geant4 Monte Carlo (MC) toolkit is a continuously evolving MC simulation code permitting mechanistic studies of cellular radiobiological effects. Geant4-DNA considers the physical, chemical, and biological stages of the action of ionizing radiation (in the form of x- and γ-ray photons, electrons and β^±-rays, hadrons, α-particles, and a set of heavier ions) in living cells towards a variety of applications ranging from predicting radiotherapy outcomes to radiation protection both on earth and in space. In this work, we provide a brief, yet concise, overview of the progress that has been achieved so far concerning the different physical, physicochemical, chemical, and biological models implemented into Geant4-DNA, highlighting the latest developments. Specifically, the “dnadamage1” and “molecularDNA” applications which enable, for the first time within an open-source platform, quantitative predictions of early DNA damage in terms of single-strand-breaks (SSBs), double-strand-breaks (DSBs), and more complex clustered lesions for different DNA structures ranging from the nucleotide level to the entire genome. These developments are critically presented and discussed along with key benchmarking results. The Geant4-DNA toolkit, through its different set of models and functionalities, offers unique capabilities for elucidating the problem of radiation quality or the relative biological effectiveness (RBE) of different ionizing radiations which underlines nearly the whole spectrum of radiotherapeutic modalities, from external high-energy hadron beams to internal low-energy gamma and beta emitters that are used in brachytherapy sources and radiopharmaceuticals, respectively.

Implementation of new physics models for low energy electrons in liquid water in Geant4-DNA

A new alternative set of elastic and inelastic cross sections has been added to the very low energy extension of the Geant4 Monte Carlo simulation toolkit, Geant4-DNA, for the simulation of electron interactions in liquid water. These cross sections have been obtained from the CPA100 Monte Carlo track structure code, which has been a reference in the microdosimetry community for many years. They are compared to the default Geant4-DNA cross sections and show better agreement with published data. In order to verify the correct implementation of the CPA100 cross section models in Geant4-DNA, simulations of the number of interactions and ranges were performed using Geant4-DNA with this new set of models, and the results were compared with corresponding results from the original CPA100 code. Good agreement is observed between the implementations, with relative differences lower than 1% regardless of the incident electron energy. Useful quantities related to the deposited energy at the scale of the cell or the organ of interest for internal dosimetry, like dose point kernels, are also calculated using these new physics models. They are compared with results obtained using the well-known Penelope Monte Carlo code.

Communication: Electron ionization of DNA bases

No reliable experimental data exist for the partial and total electron ionization cross sections for DNA bases, which are very crucial for modeling radiation damage in genetic material of living cell. We have measured a complete set of absolute partial electron ionization cross sections up to 500 eV for DNA bases for the first time by using the relative flow technique. These partial cross sections are summed to obtain total ion cross sections for all the four bases and are compared with the existing theoretical calculations and the only set of measured absolute cross sections. Our measurements clearly resolve the existing discrepancy between the theoretical and experimental results, thereby providing for the first time reliable numbers for partial and total ion cross sections for these molecules. The results on fragmentation analysis of adenine supports the theory of its formation in space.

Track structure modeling in liquid water: A review of the Geant4-DNA very low energy extension of the Geant4 Monte Carlo simulation toolkit

Understanding the fundamental mechanisms involved in the induction of biological damage by ionizing radiation remains a major challenge of today's radiobiology research. The Monte Carlo simulation of physical, physicochemical and chemical processes involved may provide a powerful tool for the simulation of early damage induction. The Geant4-DNA extension of the general purpose Monte Carlo Geant4 simulation toolkit aims to provide the scientific community with an open source access platform for the mechanistic simulation of such early damage. This paper presents the most recent review of the Geant4-DNA extension, as available to Geant4 users since June 2015 (release 10.2 Beta). In particular, the review includes the description of new physical models for the description of electron elastic and inelastic interactions in liquid water, as well as new examples dedicated to the simulation of physicochemical and chemical stages of water radiolysis. Several implementations of geometrical models of biological targets are presented as well, and the list of Geant4-DNA examples is described.

Absolute electron total ionization cross-sections: molecular analogues of DNA and RNA nucleobase and sugar constituents

Accurate ionization cross-sections for DNA and RNA constituents in the condensed or aqueous phase are important parameters for models simulating radiation damage to genetic material in living cells. In this work, absolute gas-phase electron total ionization cross-sections (TICSs) have been measured for a series of six aromatic and eight non-aromatic cyclic species that can be considered as prototype functional group analogues for the nucleobases and sugar backbone constituents of DNA and RNA. TICSs for water, hexane, and ethylacetamide (a peptide bond analogue) are also reported. The experimental apparatus utilizes a cylindrical ion collector that surrounds the ionization region, providing essentially unit detection efficiency. Two theoretical models, the polarizability-correlation method and binary-encounter Bethe theory, are able to reproduce the measured maximum TICS well for all species studied. An empirical energy-dependent correction is found to yield improvement in the agreement between experimental energy-dependent cross sections and the predictions of the BEB model. Having characterised and optimised the performance of both models, they are then used to predict TICSs for gas-phase DNA and RNA nucleobases and sugars. Direct experimental determinations of TICSs for these species are difficult because of their low volatility, which makes it difficult to prepare suitable gas-phase samples for measurement.

Semiempirical model for the ion impact ionization of complex biological media

We present a semiempirical model for calculating the electron emission from any organic compound after ion impact. With only the input of the density and composition of the target we are able to evaluate its ionization cross sections using plausible approximations. Results for protons impacting in the most representative biological targets (such as water or DNA components) show a very good comparison with experimental data. Because of its simplicity and great predictive effectiveness, the method can be immediately extended to any combination of biological target and charged particle of interest in ion beam cancer therapy.

Doses and risks from uranium are not increased significantly by interactions with natural background photon radiation

The impact of depleted uranium (DU) on human health has been the subject of much conjecture. Both the chemical and radiological aspects of its behaviour in the human body have previously been investigated in detail, with the radiological impact being assumed to be linked to the alpha decay of uranium. More recently, it has been proposed that the accumulation in tissue of high-Z materials, such as DU, may give rise to enhanced local energy deposition in the presence of natural background photon radiation due to the high photoelectric interaction cross sections of high-Z atoms. It is speculated that, in addition to producing short-range photoelectrons, these events will be followed by intense Auger and Coster-Kronig electron emission, thereby causing levels of cell damage that are unaccounted for in conventional models of radiological risk. In this study, the physical and biological bases of these claims are investigated. The potential magnitudes of any effect are evaluated and discussed, and compared with the risks from other radiological or chemical hazards. Monte Carlo calculations are performed to estimate likely energy depositions due to the presence of uranium in human tissues in photon fields: whole body doses, organ doses in anthropomorphic phantoms and nano-/micro-dosimetric scenarios are each considered. The proposal is shown generally to be based on sound physics, but overall the impact on human health is expected to be negligible.

Auger processes in the 21st century

PURPOSE

The extreme radiotoxicity of Auger electrons and their exquisite capacity to irradiate specific molecular sites has prompted scientists to extensively investigate their radiobiological effects. Their efforts have been punctuated by quadrennial international symposia that have focused on biophysical aspects of Auger processes. The latest meeting, the 6th International Symposium on Physical, Molecular, Cellular, and Medical Aspects of Auger Processes, was held 5-6 July 2007 at Harvard Medical School in Boston, Massachusetts, USA. This article provides a review of the research in this field that was published during the years 2004-2007, the period that has elapsed since the previous meeting.

CONCLUSION

The field has advanced considerably. A glimpse of the potential of this unique form of ionizing radiation to contribute to future progress in a variety of fields of study is proffered.

Computer simulation of strand break yields in plasmid pBR322: DNA damage following 125I decay

This paper presents results of (125)I effects on plasmid pBR322 in aqueous solution, simulating the complete transport of Auger and X rays up to the chemical phase. In addition to new sampling algorithms, new electronic cross sections are included. Simulations were carried out both with (125)I, bound to plasmid, or free, in its vicinity. The influence of the hydroxyl radical scavenger dimethyl sulfoxyde (DMSO) has also been tested, underlying that, in naked DNA, double strand breaks (caused by the decay of bound (125)I) are mainly due to direct hits. The calculated yields of relaxation events (RE) and linearization events (LE) show good agreement with experimental ones: when (125)I is bound to the plasmid pBR322, 0.16 RE and 0.83 LE per decay (without DMSO) are then observed. Then, when 2 mol DMSO is added, RE and LE probabilities become 0.22 and 0.76. The very light differences with those from literature could arise from experimental conditions.