Premium
Independent reaction times method in Geant4‐DNA: Implementation and performance
Author(s) -
RamosMéndez José,
Shin WookGeun,
Karamitros Mathieu,
DomínguezKondo Jorge,
Tran Ngoc Hoang,
Incerti Sebastien,
Villagrasa Carmen,
Perrot Yann,
Štěpán Václav,
Okada Shogo,
MorenoBarbosa Eduardo,
Faddegon Bruce
Publication year - 2020
Publication title -
medical physics
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 1.473
H-Index - 180
eISSN - 2473-4209
pISSN - 0094-2405
DOI - 10.1002/mp.14490
Subject(s) - monte carlo method , linear energy transfer , radiolysis , electron , ionizing radiation , dosimetry , range (aeronautics) , materials science , computational physics , chemistry , physics , irradiation , nuclear physics , mathematics , nuclear medicine , medicine , statistics , composite material
Purpose The simulation of individual particle tracks and the chemical stage following water radiolysis in biological tissue is an effective means of improving our knowledge of the physico‐chemical contribution to the biological effect of ionizing radiation. However, the step‐by‐step simulation of the reaction kinetics of radiolytic species is the most time‐consuming task in Monte Carlo track‐structure simulations, with long simulation times that are an impediment to research. In this work, we present the implementation of the independent reaction times (IRT) method in Geant4‐DNA Monte Carlo toolkit to improve the computational efficiency of calculating G‐values, defined as the number of chemical species created or lost per 100 eV of deposited energy. Methods The computational efficiency of IRT, as implemented, is compared to that from available Geant4‐DNA step‐by‐step simulations for electrons, protons and alpha particles covering a wide range of linear energy transfer (LET). The accuracy of both methods is verified using published measured data from fast electron irradiations for • OH and e aq ‐ for time‐dependent G‐values. For IRT, simulations in the presence of scavengers irradiated by cobalt‐60 γ‐ray and 2 MeV protons are compared with measured data for different scavenging capacities. In addition, a qualitative assessment comparing measured LET‐dependent G‐values with Geant4‐DNA calculations in pure liquid water is presented. Results The IRT improved the computational efficiency by three orders of magnitude relative to the step‐by‐step method while differences in G‐values by 3.9% at 1 μs were found. At 7 ps, • OH and e aq ‐ yields calculated with IRT differed from recent published measured data by 5% ± 4% and 2% ± 4%, respectively. At 1 μs, differences were 9% ± 5% and 6% ± 7% for • OH and e aq ‐ , respectively. Uncertainties are one standard deviation. Finally, G‐values at different scavenging capacities and LET‐dependent G‐values reproduced the behavior of measurements for all radiation qualities. Conclusion The comprehensive validation of the Geant4‐DNA capabilities to accurately simulate the chemistry following water radiolysis is an ongoing work. The implementation presented in this work is a necessary step to facilitate performing such a task.