Search

Ultra-high dose rate FLASH radiotherapy, a promising cancer treatment approach, offers the potential to reduce healthy tissue damage during radiotherapy. As the mechanisms underlying this process remain unknown, several hypotheses have been proposed, including the altered production of radio-induced species under ultra-high dose rate (UHDR) conditions. This study explores realistic irradiation scenarios with various dose-per-pulse and investigates the role of pulse temporal structure. Using the Geant4 toolkit and its Geant4-DNA extension, we modeled the Oriatron eRT6 linac, a FLASH-validated electron beam, and conducted simulations covering four distinct dose-per-pulse scenarios - 0.17 Gy, 1 Gy, 5 Gy, and 10 Gy - all featuring a 1.8 µs pulse duration. Results show close agreement between simulated and experimental dose profiles in water, validating the eRT6 model for Geant4-DNA simulations. We observed important changes in the temporal evolution of certain species, such as the earlier fall in hydroxyl radicals ([Formula: see text]) and reduced production and lifetime of superoxide ([Formula: see text]) with higher dose-per-pulse levels. The pulse temporal structure did not influence the long-term evolution of species. Our findings encourage further investigation into different irradiation types, such as multi-pulse configurations, and emphasize the need to add components in water to account for relevant cellular processes., (© 2024. The Author(s).)

Purpose: This paper presents the capabilities of the Geant4-DNA Monte Carlo toolkit to simulate water radiolysis with scavengers using the step-by-step (SBS) or the independent reaction times (IRT) methods. It features two examples of application areas: (1) computing the escape yield of H 2 O 2 following a 60 Co γ-irradiation and (2) computing the oxygen depletion in water irradiated with 1 MeV electrons., Methods: To ease the implementation of the chemical stage in Geant4-DNA, we developed a user interface that helps define the chemical reactions and set the concentration of scavengers. The first application area example required two computational steps to perform water radiolysis using NO 2 - and NO 3 - as scavengers and a 60 Co irradiation. The oxygen depletion computation technique for the second application area example consisted of simulating track segments of 1 MeV electrons and determining the radio-induced loss and gain of oxygen molecules., Results: The production of H 2 O 2 under variable scavenging levels is consistent with the literature; the mean relative difference between the SBS and IRT methods is 7.2 % ± 0.5 %. For the oxygen depletion 1 µs post-irradiation, the mean relative difference between both methods is equal to 9.8 % ± 0.3 %. The results in the microsecond scale depend on the initial partial pressure of oxygen in water. In addition, the computed oxygen depletions agree well with the literature., Conclusions: The Geant4-DNA toolkit makes it possible to simulate water radiolysis in the presence of scavengers. This feature offers perspectives in radiobiology, with the possibility of simulating cell-relevant scavenging mechanisms., Competing Interests: Declaration of Competing Interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2023 Associazione Italiana di Fisica Medica e Sanitaria. Published by Elsevier Ltd. All rights reserved.)

Purpose: We explored different technologies to minimize simulation time of the Monte-Carlo method for track generation following the Geant4-DNA processes for electrons in water., Methods: A GPU software tool is developed for electron track simulations. A similar CPU version is also developed using the same collision models. CPU simulations were carried out on a single user desktop computer and on the computing grid France Grilles using 10 and 100 computing nodes. Computing time results for CPU, GPU, and grid simulations are compared with those using Geant4-DNA processes., Results: The CPU simulations better performs when the number of electrons is less than 10 4 with 100 eV initial energy, this number decreases as the energy increases. The GPU simulations gives better results when the number of electrons is more than 10 4 with initial energy of 100 eV, this number decreases to 10 3 for electrons with 10KeV and increases back with higher energy. The use of the grid introduces an additional queuing time which slows down the overall simulation performance. Thus, the Grid gives better performance when the number of electrons is over 10 5 with initial energy of 10KeV, and this number decreases as the energy increases., Conclusions: The CPU is best suited for small numbers of primary incident electrons. The GPU is best suited when the number of primary incident particles occupies sufficient resources on GPU card in order to get an important computing power. The grid is best suited for simulations with high number of primary incident electrons with high initial energy., Competing Interests: Declaration of Competing Interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2022 Associazione Italiana di Fisica Medica e Sanitaria. Published by Elsevier Ltd. All rights reserved.)

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-dependent G values in pure liquid water for · OH, e - aq , and H 2 with published experimental data. For H 2 O 2 and 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 by 137 Cs γ -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-dependent G values 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., (© 2021 Institute of Physics and Engineering in Medicine.)

This paper demonstrates the impact of the pre-chemical stage, especially the dissociation scheme and the associated probabilities, on water radiolysis simulation using the Geant4-DNA Monte Carlo track structure simulation toolkit. The models and parameters provided by TRACs have been collected and implemented into Geant4-DNA. In order to evaluate their influence on water radiolysis simulation, the radiochemical yields (G-values) are evaluated as a function of time and LET using the "chem6" Geant4-DNA example, and they are compared with published experimental and calculated data. The new pre-chemical models lead to a better agreement with literature data than the default pre-chemical models of Geant4-DNA, especially for OH radicals and H 2 O 2 . The revised chemistry constructor "G4EmDNAChemistry_option3" is available in Geant4-DNA version 10.7., (Copyright © 2021 Associazione Italiana di Fisica Medica. Published by Elsevier Ltd. All rights reserved.)

In this work, we use the next sub-volume method (NSM) to investigate the possibility of using the compartment-based ("on-lattice") model to simulate water radiolysis. We, first, start with a brief description of the reaction-diffusion master equation (RDME) in a spatially discretized simulation volume ("mesh"), which is divided into sub-volumes (or "voxels"). We then discuss the choice of voxel size and merging technique of a given mesh, along with the evolution of the system using the hierarchical algorithm for the RDME ("hRDME"). Since the compartment-based model cannot describe high concentration species of early radiation-induced spurs, we propose a combination of the particle-based step-by-step ("SBS") Brownian dynamics model and the compartment-based model ("SBS-RDME model") for the simulation. We, finally, use the particle-based SBS Brownian dynamics model of Geant4-DNA as a reference to test the model implementation through several benchmarks. We find that the compartment-based model can efficiently simulate the system with a large number of species and for longer timescales, beyond the microsecond, with a reasonable computing time. Our aim in developing this model is to study the production and evolution of reactive oxygen species generated under irradiation with different dose rate conditions, such as in FLASH and conventional radiotherapy.

Purpose: Track structure simulation codes can accurately reproduce the stochastic nature of particle-matter interactions in order to evaluate quantitatively radiation damage in biological cells such as DNA strand breaks and base damage. Such simulations handle large numbers of secondary charged particles and molecular species created in the irradiated medium. Every particle and molecular species are tracked step-by-step using a Monte Carlo method to calculate energy loss patterns and spatial distributions of molecular species inside a cell nucleus with high spatial accuracy. The Geant4-DNA extension of the Geant4 general-purpose Monte Carlo simulation toolkit allows for such track structure simulations and can be run on CPUs. However, long execution times have been observed for the simulation of DNA damage in cells. We present in this work an improvement of the computing performance of such simulations using ultraparallel processing on a graphical processing unit (GPU)., Methods: A new Monte Carlo simulator named MPEXS-DNA, allowing high computing performance by using a GPU, has been developed for track structure and radiolysis simulations at the subcellular scale. MPEXS-DNA physics and chemical processes are based on Geant4-DNA processes available in Geant4 version 10.02 p03. We have reimplemented the Geant4-DNA process codes of the physics stage (electromagnetic processes of charged particles) and the chemical stage (diffusion and chemical reactions for molecular species) for microdosimetry simulation by using the CUDA language. MPEXS-DNA can calculate a distribution of energy loss in the irradiated medium caused by charged particles and also simulate production, diffusion, and chemical interactions of molecular species from water radiolysis to quantitatively assess initial damage to DNA. The validation of MPEXS-DNA physics and chemical simulations was performed by comparing various types of distributions, namely the radial dose distributions for the physics stage, and the G-value profiles for each chemical product and their linear energy transfer dependency for the chemical stage, to existing experimental data and simulation results obtained by other simulation codes, including PARTRAC., Results: For physics validation, radial dose distributions calculated by MPEXS-DNA are consistent with experimental data and numerical simulations. For chemistry validation, MPEXS-DNA can also reproduce G-value profiles for each molecular species with the same tendency as existing experimental data. MPEXS-DNA also agrees with simulations by PARTRAC reasonably well. However, we have confirmed that there are slight discrepancies in G-value profiles calculated by MPEXS-DNA for molecular species such as H 2 and H 2 O 2 when compared to experimental data and PARTRAC simulations. The differences in G-value profiles between MPEXS-DNA and PARTRAC are caused by the different chemical reactions considered. MPEXS-DNA can drastically boost the computing performance of track structure and radiolysis simulations. By using NVIDIA's GPU devices adopting the Volta architecture, MPEXS-DNA has achieved speedup factors up to 2900 against Geant4-DNA simulations with a single CPU core., Conclusion: The MPEXS-DNA Monte Carlo simulation achieves similar accuracy to Monte Carlo simulations performed using other codes such as Geant4-DNA and PARTRAC, and its predictions are consistent with experimental data. Notably, MPEXS-DNA allows calculations that are, at maximum, 2900 times faster than conventional simulations using a CPU., (© 2018 The Authors Medical Physics published by Wiley Periodicals, Inc. on behalf of American Association of Physicists in Medicine.)

Purpose: Indirect biological damage due to reactive species produced in water radiolysis reactions is responsible for the majority of biological effect for low linear energy transfer (LET) radiation. Modeling water radiolysis and the subsequent interactions of reactive species, as well as track structures, is essential to model radiobiology on the microscale. Recently, chemistry models have been developed for Geant4-DNA to be used in combination with the comprehensive existing physics models. In the current work, the first detailed, independent, in silico validation of all species yields with published experimental observations and comparison with other radiobiological simulations is presented. Additionally, the effect of LET of protons and heavier ions on reactive species yield in the model was examined, as well as the completeness of the chemical reactions following the radiolysis within the time after physical interactions simulated in the model., Methods: Yields over time of reactive species were simulated for water radiolysis by incident electrons, protons, alpha particles, and ions with various LETs using Geant4 and RITRACKS simulation tools. Water dissociation and recombination was simulated using Geant4 to determine the completeness of chemical reactions at the end of the simulation. Yield validation was performed by comparing yields simulated using Geant4 with experimental observations and other simulations. Validation was performed for all species for low LET radiation and the solvated electron and hydroxyl radical for high LET ions., Results: It was found that the Geant4-DNA chemistry yields were generally in good agreement with experimental observations and other simulations. However, the Geant4-DNA yields for the hydroxyl radical and hydrogen peroxide at the end of the chemistry stage were found to be respectively considerably higher and lower than the experimentally observed yields. Increasing the LET of incident hadrons increased the yield of secondary species and decreased the yield of primary species. The effect of LET on the yield of the hydroxyl radical at 100 ns simulated with Geant4 was in good agreement with experimental measurements. Additionally, by the end of the simulation only 40% of dissociated water molecules had been recombined and the rate of recombination was slowing., Conclusions: The yields simulated using Geant4 are within reasonable agreement with experimental observations. Higher LET radiation corresponds with increased yields of secondary species and decreased yields of primary species. These trends combined with the LET having similar effects on the 100 ns hydroxyl radical yield for Geant4 and experimental measurements indicate that Geant4 accurately models the effect of LET on radiolysis yields. The limited recombination within the modeled chemistry stage and the slowing rate of recombination at the end of the stage indicate potential long-range indirect biological damage., (© 2018 American Association of Physicists in Medicine.)

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)., (Copyright © 2019 Associazione Italiana di Fisica Medica. Published by Elsevier Ltd. All rights reserved.)

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., (Copyright © 2016 Associazione Italiana di Fisica Medica. Published by Elsevier Ltd. All rights reserved.)

Radiation damage to the central nervous system (CNS) has been an on-going challenge for the last decades primarily due to the issues of brain radiotherapy and radiation protection for astronauts during space travel. Although recent findings revealed a number of molecular mechanisms associated with radiation-induced impairments in behaviour and cognition, some uncertainties exist in the initial neuronal cell injury leading to the further development of CNS malfunction. The present study is focused on the investigation of early biological damage induced by ionizing radiations in a sample neural network by means of modelling physico-chemical processes occurring in the medium after exposure. For this purpose, the stochastic simulation of incident particle tracks and water radiation chemistry was performed in realistic neuron phantoms constructed using experimental data on cell morphology. The applied simulation technique is based on using Monte-Carlo processes of the Geant4-DNA toolkit. The calculations were made for proton, 12 C, and 56 Fe particles of different energy within a relatively wide range of linear energy transfer values from a few to hundreds of keV/μm. The results indicate that the neuron morphology is an important factor determining the accumulation of microscopic radiation dose and water radiolysis products in neurons. The estimation of the radiolytic yields in neuronal cells suggests that the observed enhancement in the levels of reactive oxygen species may potentially lead to oxidative damage to neuronal components disrupting the normal communication between cells of the neural network., (Copyright © 2016 Associazione Italiana di Fisica Medica. Published by Elsevier Ltd. All rights reserved.)

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., (Copyright © 2015 Associazione Italiana di Fisica Medica. Published by Elsevier Ltd. All rights reserved.)

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.

This work presents a Monte Carlo study of energy depositions due to protons, alpha particles and carbon ions of the same linear-energy-transfer (LET) in liquid water. The corresponding track structures were generated using the Geant4-DNA toolkit, and the energy deposition spatial distributions were analyzed using an adapted version of the DBSCAN clustering algorithm. Combining the Geant4 simulations and the clustering algorithm it was possible to compare the quality of the different radiation types. The ratios of clustered and single energy depositions are shown versus particle LET and frequency-mean lineal energies. The estimated effect of these types of radiation on biological tissues is then discussed by comparing the results obtained for different particles with the same LET.

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.

This paper presents the influence of electron elastic scattering models, electron thermalization models, and chemical parameters on Geant4-DNA simulations of liquid water radiolysis under mega-electron-volt electron irradiation. The radiochemical yields are simulated using a new Geant4-DNA example. In particular, the influence of the new elastic scattering model recently developed is presented as well as the influence of improved electron thermalization models. The influence of a new chemistry constructor using parameters of another Monte Carlo track structure code is also described. The results calculated using these different models are compared with each other and with experimental data. For sub-mega-electron-volt electron simulations, the combination of the "G4EmDNAPhysics_option2" physics constructor with the recently developed elastic scattering model, the Meesungnoen electron thermalization model, and the "G4EmDNAChemistry_option1" chemistry constructor is recommended. [ABSTRACT FROM AUTHOR]

Abstract: In the perspective of building an open source simulation platform dedicated to the modelling of early biological molecular damages due to ionising radiation at the DNA scale, the general-purpose Geant4 Monte Carlo simulation toolkit has been recently extended with specific very low energy electromagnetic physics processes for liquid water medium. These processes – also called “Geant4-DNA” processes – simulate the physical interactions induced by electrons, hydrogen and helium atoms of different charge states. The present work reports on the energy deposit distributions obtained for incident electrons, protons and alpha particles in nanometre-size volumes comparable to those present in the genetic material of mammalian cells. The frequency distributions of the energy deposition obtained for three typical geometries of nanometre-size cylindrical targets placed in a spherical phantom are found to be in reasonable agreement with prior works. Furthermore, we present a combination of the Geant4-DNA processes with a simplified geometrical model of a cellular nucleus allowing the evaluation of energy deposits in volumes of biological interest. [Copyright &y& Elsevier]

Accurate modeling of DNA damages resulting from ionizing radiation remains a challenge of today’s radiobiology research. An original set of physics processes has been recently developed for modeling the detailed transport of protons and neutral hydrogen atoms in liquid water and in DNA nucleobases using the Geant4-DNA extension of the open source Geant4 Monte Carlo simulation toolkit. The theoretical cross sections as well as the mean energy transfers during the different ionizing processes were taken from recent works based on classical as well as quantum mechanical predictions. Furthermore, in order to compare energy deposition patterns in liquid water and DNA material, we here propose a simplified cellular nucleus model made of spherical voxels, each containing randomly oriented nanometer-size cylindrical targets filled with either liquid water or DNA material (DNA nucleobases) both with a density of 1g/cm3. These cylindrical volumes have dimensions comparable to genetic material units of mammalian cells, namely, 25nm (diameter)×25nm (height) for chromatin fiber segments, 10nm (d)×5nm (h) for nucleosomes and 2nm (d)×2nm (h) for DNA segments. Frequencies of energy deposition in the cylindrical targets are presented and discussed. [ABSTRACT FROM AUTHOR]

Abstract: The present work provides an accurate description of the elastic scattering process for low-energy electrons (10eV–10keV) in liquid water by means of a free-parameter quantum-mechanical treatment. The calculations are performed in the partial-wave formalism by means of a total interaction potential taking into account a static contribution as well as fine effects like exchange and polarization contributions. The obtained results in terms of singly differential and total cross sections exhibit relatively good agreement with available experimental data (in gaseous water). They have been incorporated into the Geant4 toolkit, which has been recently extended with physics processes for microdosimetry applications in liquid water down to the electronvolt scale. They offer an improved alternative to the semi-empirical and to the screened Rutherford models already available in this very low-energy extension. [Copyright &y& Elsevier]