ABSTRACT
Cosmic rays have the potential to significantly affect the atmospheric composition by increasing the rate and changing the types of chemical reactions through ion production. The amount and states of ionization, and the spatial distribution of ions produced are still open questions for atmospheric models. To precisely estimate these quantities, it is necessary to simulate particle-molecule interactions, down to very low energies. Models enabling such simulations require interaction probabilities over a broad energy range and for all energetically allowed scattering processes. In this paper, we focus on electron interaction with the two most abundant molecules in the atmosphere, i.e., N2 and O2, as an initial step. A set of elastic and inelastic cross section models for electron transportation in oxygen and nitrogen molecules valid in the energy range 10 eV - 1 MeV, is presented. Comparison is made with available theoretical and experimental data and a reasonable good agreement is observed. Stopping power is calculated and compared with published data to assess the general consistency and reliability of our results. Good overall agreement is observed, with relative differences lower than 6% with the ESTAR database.
Subject(s)
Electrons , Oxygen , Reproducibility of Results , Monte Carlo Method , Physical Phenomena , Ions , Water/chemistryABSTRACT
BACKGROUND: Geant4 is a Monte Carlo code extensively used in medical physics for a wide range of applications, such as dosimetry, micro- and nanodosimetry, imaging, radiation protection, and nuclear medicine. Geant4 is continuously evolving, so it is crucial to have a system that benchmarks this Monte Carlo code for medical physics against reference data and to perform regression testing. AIMS: To respond to these needs, we developed G4-Med, a benchmarking and regression testing system of Geant4 for medical physics. MATERIALS AND METHODS: G4-Med currently includes 18 tests. They range from the benchmarking of fundamental physics quantities to the testing of Monte Carlo simulation setups typical of medical physics applications. Both electromagnetic and hadronic physics processes and models within the prebuilt Geant4 physics lists are tested. The tests included in G4-Med are executed on the CERN computing infrastructure via the use of the geant-val web application, developed at CERN for Geant4 testing. The physical observables can be compared to reference data for benchmarking and to results of previous Geant4 versions for regression testing purposes. RESULTS: This paper describes the tests included in G4-Med and shows the results derived from the benchmarking of Geant4 10.5 against reference data. DISCUSSION: Our results indicate that the Geant4 electromagnetic physics constructor G4EmStandardPhysics_option4 gives a good agreement with the reference data for all the tests. The QGSP_BIC_HP physics list provided an overall adequate description of the physics involved in hadron therapy, including proton and carbon ion therapy. New tests should be included in the next stage of the project to extend the benchmarking to other physical quantities and application scenarios of interest for medical physics. CONCLUSION: The results presented and discussed in this paper will aid users in tailoring physics lists to their particular application.
Subject(s)
Benchmarking , Physics , Radiometry , Computer Simulation , Monte Carlo MethodABSTRACT
The Geant4 toolkit offers a range of electromagnetic (EM) models for simulating the transport of charged particles down to sub-keV energies. They can be divided to condensed-history (CH) models (like the Livermore and Penelope models) and the track-structure (TS) models included in the Geant4-DNA low-energy extension of Geant4. Although TS models are considered the state-of-the-art for nanoscale electron transport, they are difficult to develop, computationally intensive, and commonly tailored to a single medium (e.g., water) which prohibits their use in a wide range of applications. Thus, the use of CH models down to sub-keV energies is particularly intriguing in the context of general-purpose Monte Carlo codes. The aim of the present work is to compare the performance of the CH models of Geant4 against the recently implemented TS models of Geant4-DNA for nanoscale electron transport. Calculations are presented for two fundamental quantities, the dose-point-kernel and the microdosimetric lineal energy. The influence of user-defined simulation parameters (tracking and production cuts, and maximum step size) on the above calculations is also examined. It is shown that Livermore offers the best performance among the CH models of Geant4 for nanoscale electron transport. However, even under optimally-chosen simulation parameters, the differences between the CH and TS models examined may be sizeable for low energy electrons (<1 keV) and/or nanometer size targets (<100â¯nm).
Subject(s)
Monte Carlo Method , Water/chemistry , Electron Transport , RadiobiologyABSTRACT
Our understanding of radiation-induced cellular damage has greatly improved over the past few decades. Despite this progress, there are still many obstacles to fully understand how radiation interacts with biologically relevant cellular components, such as DNA, to cause observable end points such as cell killing. Damage in DNA is identified as a major route of cell killing. One hurdle when modeling biological effects is the difficulty in directly comparing results generated by members of different research groups. Multiple Monte Carlo codes have been developed to simulate damage induction at the DNA scale, while at the same time various groups have developed models that describe DNA repair processes with varying levels of detail. These repair models are intrinsically linked to the damage model employed in their development, making it difficult to disentangle systematic effects in either part of the modeling chain. These modeling chains typically consist of track-structure Monte Carlo simulations of the physical interactions creating direct damages to DNA, followed by simulations of the production and initial reactions of chemical species causing so-called "indirect" damages. After the induction of DNA damage, DNA repair models combine the simulated damage patterns with biological models to determine the biological consequences of the damage. To date, the effect of the environment, such as molecular oxygen (normoxic vs. hypoxic), has been poorly considered. We propose a new standard DNA damage (SDD) data format to unify the interface between the simulation of damage induction in DNA and the biological modeling of DNA repair processes, and introduce the effect of the environment (molecular oxygen or other compounds) as a flexible parameter. Such a standard greatly facilitates inter-model comparisons, providing an ideal environment to tease out model assumptions and identify persistent, underlying mechanisms. Through inter-model comparisons, this unified standard has the potential to greatly advance our understanding of the underlying mechanisms of radiation-induced DNA damage and the resulting observable biological effects when radiation parameters and/or environmental conditions change.
Subject(s)
DNA Damage , Computer Simulation , DNA Repair , Linear Energy Transfer , Models, Theoretical , Monte Carlo MethodABSTRACT
This Special Report presents a description of Geant4-DNA user applications dedicated to the simulation of track structures (TS) in liquid water and associated physical quantities (e.g., range, stopping power, mean free path ). These example applications are included in the Geant4 Monte Carlo toolkit and are available in open access. Each application is described and comparisons to recent international recommendations are shown (e.g., ICRU, MIRD), when available. The influence of physics models available in Geant4-DNA for the simulation of electron interactions in liquid water is discussed. Thanks to these applications, the authors show that the most recent sets of physics models available in Geant4-DNA (the so-called "option4" and "option 6" sets) enable more accurate simulation of stopping powers, dose point kernels, and W-values in liquid water, than the default set of models ("option 2") initially provided in Geant4-DNA. They also serve as reference applications for Geant4-DNA users interested in TS simulations.
ABSTRACT
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.
Subject(s)
DNA/chemistry , Models, Molecular , Monte Carlo Method , Water/chemistry , Chemical Phenomena , HumansABSTRACT
PURPOSE: The geant4-DNA physics models are upgraded by a more accurate set of electron cross sections for ionization and excitation in liquid water. The impact of the new developments on low-energy electron transport simulations by the geant4 Monte Carlo toolkit is examined for improving its performance in dosimetry applications at the subcellular and nanometer level. METHODS: The authors provide an algorithm for an improved implementation of the Emfietzoglou model dielectric response function of liquid water used in the geant4-DNA existing model. The algorithm redistributes the imaginary part of the dielectric function to ensure a physically motivated behavior at the binding energies, while retaining all the advantages of the original formulation, e.g., the analytic properties and the fulfillment of the f-sum-rule. In addition, refinements in the exchange and perturbation corrections to the Born approximation used in the geant4-DNA existing model are also made. RESULTS: The new ionization and excitation cross sections are significantly different from those of the geant4-DNA existing model. In particular, excitations are strongly enhanced relative to ionizations, resulting in higher W-values and less diffusive dose-point-kernels at sub-keV electron energies. CONCLUSIONS: An improved energy-loss model for the excitation and ionization of liquid water by low-energy electrons has been implemented in geant4-DNA. The suspiciously low W-values and the unphysical long tail in the dose-point-kernel have been corrected owing to a different partitioning of the dielectric function.
Subject(s)
Algorithms , Electron Transport , Models, Chemical , Water/chemistry , Computer Simulation , Electrons , Ions/chemistry , Monte Carlo MethodABSTRACT
PURPOSE: We provide a short overview of optical-data models for the description of inelastic scattering of low-energy electrons (10-10,000 eV) in liquid water. The effect on the inelastic scattering cross section due to different optical data and extension algorithms is examined with emphasis on some recent developments. MATERIALS AND METHODS: The optical-data method whereby experimental optical data and theoretical extension algorithms are used to describe the dependence of the dielectric response function on energy- and momentum-transfer and obtain the Bethe surface of the material, currently represents the most used method for computing the inelastic scattering of low-energy electrons in condensed media. Two sets of experimental optical data for liquid water obtained from reflectance and inelastic X-ray scattering spectroscopy, respectively, and the extension algorithms of Ritchie, Penn, and Ashley are examined. Recent developments are discussed along with the role of corrections to the random phase approximation (RPA) of electron gas theory. RESULTS: The inelastic scattering cross section in the energy range 200-10,000 eV was found to be rather insensitive (to within 10%) to the choice of optical data or the extension algorithm. In contrast, differences between model calculations increase rapidly below 200 eV with the influence of the extension algorithm being dominant. CONCLUSION: The choice of the extension algorithm used to extrapolate optical data to finite momentum transfer and obtain the Bethe surface is crucial in modelling the inelastic scattering of electrons with energies below 200 eV. A new set of measurements on the dielectric response function of liquid water beyond the optical limit and the development of extension algorithms that will go beyond RPA by considering the effect of (short-range) electron exchange and correlation should be of some priority.
Subject(s)
Elasticity , Electrons , Models, Theoretical , Optical Phenomena , Water/chemistry , Surface PropertiesABSTRACT
The electronic stopping power of liquid water for protons over the 50 keV to 10 MeV energy range is studied using an improved dielectric response model which is in good agreement with the best available experimental data. The mean excitation energy (I) of stopping power theory is calculated to be 77.8 eV. Shell corrections are accounted for in a self-consistent manner through analytic dispersion relations for the momentum dependence of the dielectric function. It is shown that widely used dispersion schemes based on the random-phase approximation (RPA) can result in sizeable errors due to the neglect of damping and local field effects that lead to a momentum broadening and shifting of the energy-loss function. Low-energy Born corrections for the Barkas, Bloch and charge-state effects practically cancel out down to 100 keV proton energies. Differences with ICRU Report 49 stopping power values and earlier calculations are found to be at the approximately 20% level in the region of the stopping maximum. The present work overcomes the limitations of the Bethe formula below 1 MeV and improves the accuracy of previous calculations through a more consistent account of the dielectric response properties of liquid water.