MINAS TIRITH: a new tool for simulating radiation-induced DNA damage at the cell population level.

Objective. The mechanisms of radiation-induced DNA damage can be understood via the fundamental acquisition of knowledge through a combination of experiments and modeling. Currently, most biological experiments are performed by irradiating an entire cell population, whereas modeling of radiation-induced effects is usually performed via Monte Carlo simulations with track structure codes coupled to realistic DNA geometries of a single-cell nucleus. However, the difference in scale between the two methods hinders a direct comparison because the dose distribution in the cell population is not necessarily uniform owing to the stochastic nature of the energy deposition. Thus, this study proposed the MINAS TIRITH tool to model the distribution of radiation-induced DNA damage in a cell population.Approach. The proposed method is based on precomputed databases of microdosimetric parameters and DNA damage distributions generated using the Geant4-DNA Monte Carlo Toolkit. First, a specific energyzwas assigned to each cell of an irradiated population for a particular absorbed doseDabs,following microdosimetric formalism. Then, each cell was assigned a realistic number of DNA damage events according to the specific energyz,respecting the stochastic character of its occurrence.Main results. This study validated the MINAS TIRITH tool by comparing its results with those obtained using the Geant4-DNA track structure code and a Geant4-DNA based simulation chain for DNA damage calculation. The different elements of comparison indicated consistency between MINAS TIRITH and the Monte Carlo simulation in case of the dose distribution in the population and the calculation of the amount of DNA damage.Significance. MINAS TIRITH is a new approach for the calculation of radiation-induced DNA damage at the cell population level that facilitates reasonable simulation times compared to those obtained with track structure codes. Moreover, this tool enables a more direct comparison between modeling and biological experimentation.

Experimental validation in a neutron exposure frame of the MINAS TIRITH for cell damage simulation.

In the domains of medicine and space exploration, refining risk assessment models for protecting healthy tissue from ionizing radiation is crucial. Understanding radiation-induced effects requires biological experimentations at the cellular population level and the cellular scale modeling using Monte Carlo track structure codes. We present MINAS TIRITH, a tool using Geant4-DNA Monte Carlo-generated databases to study DNA damage distribution at the cell population scale. It introduces a DNA damage location module and proposes a method to convert double-strand breaks (DSB) into DNA Damage Response foci. We evaluate damage location precision and DSB-foci conversion parameters. MINAS TIRITH's accuracy is validated againstγ-H2AX foci distribution from cell population exposed to monoenergetic neutron beams (2.5 or 15.1 MeV) under different configurations, yielding mixed radiation fields. Strong agreement between simulation and experimental results was found demonstrating MINAS TIRITH's predictive precision in radiation-induced DNA damage topology. Additionally, modeling intercellular damage variability within a population subjected to a specific macroscopic dose identifies subpopulations, enhancing realistic fate models. This approach advances our understanding of radiation-induced effects on cellular systems for risk assessment improvement.

Calculation of electron interaction models in N2 and O2.

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.

The invariance of the total direct DNA strand break yield.

PURPOSE: The invariance of the total direct strand break yield when DNA is irradiated by different types of particles and energies has been reported by previous works. This study is intended to explain the physical causes of this behavior. METHODS: The GEANT4-DNA extension of the GEANT4 general purpose Monte Carlo simulation toolkit has been used to determine direct strand break yields induced by protons and alpha particles impacting on a B-DNA geometrical model, including five organization levels of the human genetic material. The linear energy transfer (LET) of such particles ranges from 4.8 keV/microm (10 MeV protons) to about 235 keV/microm (2 MeV alpha particles), at 5.225 pm depth (near the center of the region of interest). Direct total, single and double strand break probabilities have been determined in a liquid water homogeneous medium with a 1.06 g/cm3 density. The energetic spectra of single strand breaks (SSB), the number of energy deposition events, and the SSB/event ratio were determined. RESULTS: The target-hit probability was found to be independent of both the type and the energy of the incident particle, even if this latter is a secondary electron. This probability is determined by the geometrical properties of the system. The total strand break yield and the number of energy deposition events required to reach a certain absorbed dose were found nearly independent of the type and energy of the incident ion (proton or alpha). In contrast, the double strand break (DSB) yield was found strongly dependent on the LET of the incident radiation. CONCLUSIONS: The SSB generation process is homogeneous and independent of the LET of the particles involved, at least within the proton and alpha particle energy range here studied. The target-hit probability is only determined by the ratio between the total volume occupied by targets and that of the ROI where the radiation deposits its energy. The maximum separation distance between two adjacent SSBs to produce a DSB is the parameter that breaks the homogeneity of the target-hit process, making the DSB production process strongly heterogeneous.

Monte Carlo study of patient and medical staff radiation exposures during interventional cardiology.

The aim of this study is to assess the radiation exposure of the patient and the medical staff during interventional cardiology procedures. Realistic exposure scenarios were developed using the adult reference anthropomorphic phantoms adopted by the International Commission on Radiological Protection (ICRP110Male and ICRP110Female), and the radiation transport code Geant4 (version 10.3). The calculated equivalent and effective doses were normalised by the simulated Kerma-Area Product (KAP), resulting in two conversion coefficients HT/KAP and E/KAP. To properly evaluate the risk of exposure, several dose-dependent parameters have been investigated, namely: radiological parameters (tube kilovoltage peak (kVp), type of projection, field size (FOV)), and operator positions. Four projections (AP,PA,LAO25° and RAO25°) were simulated for three X-ray energy spectra (80,100 and 120 kVp) with four different values of FOV (15×15 cm2,20×20 cm2,25×25 cm2 and 30×30 cm2). The results showed that the conversion coefficients values increase with increasing tube voltage as well as the FOV size. Recommended projection during the interventional cardiology procedures, whenever possible, should be the PA projection rather than AP projection. The most critical projection for the patient and the main operator is the RAO25° projection and the LAO25° projection respectively. The comparison of our results with the literature data showed good agreement allowing their use in the dosimetric characterization of interventional cardiology procedures.

DNA damage modeled with Geant4-DNA: effects of plasmid DNA conformation and experimental conditions.

The chemical stage of the Monte Carlo track-structure (MCTS) code Geant4-DNA was extended for its use in DNA strand break (SB) simulations and compared against published experimental data. Geant4-DNA simulations were performed using pUC19 plasmids (2686 base pairs) in a buffered solution of DMSO irradiated by60Co or137Csγ-rays. A comprehensive evaluation of SSB yields was performed considering DMSO, DNA concentration, dose and plasmid supercoiling. The latter was measured using the super helix density value used in a Brownian dynamics plasmid generation algorithm. The Geant4-DNA implementation of the independent reaction times method (IRT), developed to simulate the reaction kinetics of radiochemical species, allowed to score the fraction of supercoiled, relaxed and linearized plasmid fractions as a function of the absorbed dose. The percentage of the number of SB after •OH + DNA and H• + DNA reactions, referred as SSB efficiency, obtained using MCTS were 13.77% and 0.74% respectively. This is in reasonable agreement with published values of 12% and 0.8%. The SSB yields as a function of DMSO concentration, DNA concentration and super helix density recreated the expected published experimental behaviors within 5%, one standard deviation. The dose response of SSB and DSB yields agreed with published measurements within 5%, one standard deviation. We demonstrated that the developed extension of IRT in Geant4-DNA, facilitated the reproduction of experimental conditions. Furthermore, its calculations were strongly in agreement with experimental data. These two facts will facilitate the use of this extension in future radiobiological applications, aiding the study of DNA damage mechanisms with a high level of detail.

TOPAS-nBio validation for simulating water radiolysis and DNA damage under low-LET irradiation.

The chemical stage of the Monte Carlo track-structure simulation code Geant4-DNA has been revised and validated. The root-mean-square (RMS) empirical parameter that dictates the displacement of water molecules after an ionization and excitation event in Geant4-DNA has been shortened to better fit experimental data. The pre-defined dissociation channels and branching ratios were not modified, but the reaction rate coefficients for simulating the chemical stage of water radiolysis were updated. The evaluation of Geant4-DNA was accomplished with TOPAS-nBio. For that, we compared predicted time-dependentGvalues in pure liquid water for·OH, e-aq, and H2with published experimental data. For H2O2and H·, simulation of added scavengers at different concentrations resulted in better agreement with measurements. In addition, DNA geometry information was integrated with chemistry simulation in TOPAS-nBio to realize reactions between radiolytic chemical species and DNA. This was used in the estimation of the yield of single-strand breaks (SSB) induced by137Csγ-ray radiolysis of supercoiled pUC18 plasmids dissolved in aerated solutions containing DMSO. The efficiency of SSB induction by reaction between radiolytic species and DNA used in the simulation was chosen to provide the best agreement with published measurements. An RMS displacement of 1.24 nm provided agreement with measured data within experimental uncertainties for time-dependentGvalues and under the presence of scavengers. SSB efficiencies of 24% and 0.5% for·OH and H·, respectively, led to an overall agreement of TOPAS-nBio results within experimental uncertainties. The efficiencies obtained agreed with values obtained with published non-homogeneous kinetic model and step-by-step Monte Carlo simulations but disagreed by 12% with published direct measurements. Improvement of the spatial resolution of the DNA damage model might mitigate such disagreement. In conclusion, with these improvements, Geant4-DNA/TOPAS-nBio provides a fast, accurate, and user-friendly tool for simulating DNA damage under low linear energy transfer irradiation.

Report on G4-Med, a Geant4 benchmarking system for medical physics applications developed by the Geant4 Medical Simulation Benchmarking Group.

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.

Comparison of GEANT4 very low energy cross section models with experimental data in water.

PURPOSE: The GEANT4 general-purpose Monte Carlo simulation toolkit is able to simulate physical interaction processes of electrons, hydrogen and helium atoms with charge states (H0, H+) and (He0, He+, He2+), respectively, in liquid water, the main component of biological systems, down to the electron volt regime and the submicrometer scale, providing GEANT4 users with the so-called "GEANT4-DNA" physics models suitable for microdosimetry simulation applications. The corresponding software has been recently re-engineered in order to provide GEANT4 users with a coherent and unique approach to the simulation of electromagnetic interactions within the GEANT4 toolkit framework (since GEANT4 version 9.3 beta). This work presents a quantitative comparison of these physics models with a collection of experimental data in water collected from the literature. METHODS: An evaluation of the closeness between the total and differential cross section models available in the GEANT4 toolkit for microdosimetry and experimental reference data is performed using a dedicated statistical toolkit that includes the Kolmogorov-Smirnov statistical test. The authors used experimental data acquired in water vapor as direct measurements in the liquid phase are not yet available in the literature. Comparisons with several recommendations are also presented. RESULTS: The authors have assessed the compatibility of experimental data with GEANT4 microdosimetry models by means of quantitative methods. The results show that microdosimetric measurements in liquid water are necessary to assess quantitatively the validity of the software implementation for the liquid water phase. Nevertheless, a comparison with existing experimental data in water vapor provides a qualitative appreciation of the plausibility of the simulation models. The existing reference data themselves should undergo a critical interpretation and selection, as some of the series exhibit significant deviations from each other. CONCLUSIONS: The GEANT4-DNA physics models available in the GEANT4 toolkit have been compared in this article to available experimental data in the water vapor phase as well as to several published recommendations on the mass stopping power. These models represent a first step in the extension of the GEANT4 Monte Carlo toolkit to the simulation of biological effects of ionizing radiation.

Monte Carlo dosimetry for targeted irradiation of individual cells using a microbeam facility.

Microbeam facilities provide a unique opportunity to investigate the effects of ionising radiation on living biological cells with a precise control of the delivered dose. This paper describes dosimetry calculations performed at the single-cell level in the microbeam irradiation facility available at the Centre d'Etudes Nucléaires de Bordeaux-Gradignan in France, using the object-oriented Geant4 Monte Carlo simulation toolkit. The cell geometry model is based on high-resolution three-dimensional voxelised phantoms of a human keratinocyte (HaCaT) cell line. Such phantoms are built from confocal microscopy imaging and from ion beam chemical elemental analysis. Results are presented for single-cell irradiation with 3 MeV incident alpha particles.

A simulation study of gold nanoparticles localisation effects on radiation enhancement at the mitochondrion scale.

This paper presents a Monte-Carlo study focusing on the effects of gold nanoparticles on the energy deposition patterns produced by incident photons in the close vicinity of the mitochondrial network modeled as a tube. Spherical shaped gold nanoparticles of 30 nm diameter were placed in a micrometric (10 × 10 × 10 µm3) water phantom containing a tube of 300 nm diameter and 5 µm length. The tube represented a mitochondrial fragment and nanoparticles were distributed in the water phantom outside the tube. Photons of 120 keV were simulated using the Geant4 Livermore processes and the Geant4-DNA electron processes to account for secondary electrons collisions. The Livermore processes took into account the Auger cascade inside the gold material. A data mining algorithm was then used to analyze the energy deposition clusters inside the water phantom and the tube. A comparison was made between the results obtained for a uniform distribution of nanoparticles and a vesicle distribution model. The results including energy deposition clusters are also compared to dose enhancement ratios.

Influence of chromatin compaction on simulated early radiation-induced DNA damage using Geant4-DNA.

PURPOSE: In this work, we present simulated double-strand breaks (DSBs) obtained for two human cell nucleus geometries. The first cell nucleus represents fibroblasts, filled with DNA molecules in different compaction forms: heterochromatin or euchromatin only. The second one represents an endothelial cell nucleus, either filled with heterochromatin only or with a uniform distribution of 48% of heterochromatin and 52% of euchromatin, obtained from measurements carried out at IRSN. Protons and alpha particles of different energies were used as projectiles. Each cell nucleus model includes a multi-scale description of the DNA target from the molecular level to the whole human genome representation. METHODS: The cell nucleus models were generated using an extended version of the DnaFabric software in which a new model of euchromatin was implemented in addition to the existing model of heterochromatin. Thus, each nucleus model contains the complete human genome (a total of 6 Gbp) in the G0/G1 phase of the cycle, filled with a continuous chromatin fiber per chromosome that can take into account the heterochromatin and the euchromatin compaction. These geometries were then exported to a simulation chain using the Monte Carlo toolkit Geant4-DNA to perform computations of the physical, physicochemical, and chemical stages, in order to evaluate the influence of chromatin compaction on DSB induction and the contribution of direct and indirect damage, as well as DSB complexity. RESULTS: More direct damage and less indirect damage were observed in the heterochromatin than in the euchromatin. Nevertheless, no difference in terms of DSB complexity was observed between those formed in the heterochromatin or the euchromatin models. Yields of DSB/Gy/Gbp show an increase when both heterochromatin and euchromatin models are taken into account, compared to when only heterochromatin is considered. CONCLUSIONS: The results presented indicate that the chromatin compaction decreases DNA damage generated by ionizing radiation and thus, DNA compaction should be considered for the simulation of DNA repair and other cellular outcomes.

TOPAS-nBio: An Extension to the TOPAS Simulation Toolkit for Cellular and Sub-cellular Radiobiology.

The TOPAS Monte Carlo (MC) system is used in radiation therapy and medical imaging research, having played a significant role in making Monte Carlo simulations widely available for proton therapy related research. While TOPAS provides detailed simulations of patient scale properties, the fundamental unit of the biological response to radiation is a cell. Thus, our goal was to develop TOPAS-nBio, an extension of TOPAS dedicated to advance understanding of radiobiological effects at the (sub-)cellular, (i.e., the cellular and sub-cellular) scale. TOPAS-nBio was designed as a set of open source classes that extends TOPAS to model radiobiological experiments. TOPAS-nBio is based on and extends Geant4-DNA, which extends the Geant4 toolkit, the basis of TOPAS, to include very low-energy interactions of particles down to vibrational energies, explicitly simulates every particle interaction (i.e., without using condensed histories) and propagates radiolysis products. To further facilitate the use of TOPAS-nBio, a graphical user interface was developed. TOPAS-nBio offers full track-structure Monte Carlo simulations, integration of chemical reactions within the first millisecond, an extensive catalogue of specialized cell geometries as well as sub-cellular structures such as DNA and mitochondria, and interfaces to mechanistic models of DNA repair kinetics. We compared TOPAS-nBio simulations to measured and published data of energy deposition patterns and chemical reaction rates (G values). Our simulations agreed well within the experimental uncertainties. Additionally, we expanded the chemical reactions and species provided in Geant4-DNA and developed a new method based on independent reaction times (IRT), including a total of 72 reactions classified into 6 types between neutral and charged species. Chemical stage simulations using IRT were a factor of 145 faster than with step-by-step tracking. Finally, we applied the geometric/chemical modeling to obtain initial yields of double-strand breaks (DSBs) in DNA fibers for proton irradiations of 3 and 50 MeV and compared the effect of including chemical reactions on the number and complexity of DSB induction. Over half of the DSBs were found to include chemical reactions with approximately 5% of DSBs caused only by chemical reactions. In conclusion, the TOPAS-nBio extension to the TOPAS MC application offers access to accurate and detailed multiscale simulations, from a macroscopic description of the radiation field to microscopic description of biological outcome for selected cells. TOPAS-nBio offers detailed physics and chemistry simulations of radiobiological experiments on cells simulating the initially induced damage and links to models of DNA repair kinetics.

Influence of track structure and condensed history physics models of Geant4 to nanoscale electron transport in liquid water.

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).

A New Standard DNA Damage (SDD) Data Format.

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.

Geant4-DNA example applications for track structure simulations in liquid water: A report from the Geant4-DNA Project.

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.

A comparison of cellular irradiation techniques with alpha particles using the Geant4 Monte Carlo simulation toolkit.

A comparison of three cellular irradiation techniques using the Monte Carlo simulation toolkit Geant4 is presented in this paper. They involve electrodeposited source of alpha particle-emitting radionuclides, random classical alpha beam irradiation and single cell targeted irradiation using a focused alpha microbeam line. The simulation allows the calculation of hit distributions among the cellular population as well as the absorbed dose for two typical cellular geometries.

Review of Geant4-DNA applications for micro and nanoscale simulations.

Emerging radiotherapy treatments including targeted particle therapy, hadron therapy or radiosensitisation of cells by high-Z nanoparticles demand the theoretical determination of radiation track structure at the nanoscale. This is essential in order to evaluate radiation damage at the cellular and DNA level. Since 2007, Geant4 offers physics models to describe particle interactions in liquid water at the nanometre level through the Geant4-DNA Package. This package currently provides a complete set of models describing the event-by-event electromagnetic interactions of particles with liquid water, as well as developments for the modelling of water radiolysis. Since its release, Geant4-DNA has been adopted as an investigational tool in kV and MV external beam radiotherapy, hadron therapies using protons and heavy ions, targeted therapies and radiobiology studies. It has been benchmarked with respect to other track structure Monte Carlo codes and, where available, against reference experimental measurements. While Geant4-DNA physics models and radiolysis modelling functionalities have already been described in detail in the literature, this review paper summarises and discusses a selection of representative papers with the aim of providing an overview of a) geometrical descriptions of biological targets down to the DNA size, and b) the full spectrum of current micro- and nano-scale applications of Geant4-DNA.

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.