Independent reaction times method in Geant4‐DNA: Implementation and performance
Langue
en
Article de revue
Ce document a été publié dans
Medical Physics. 2020-09-24, vol. 47, n° 11, p. 5919-5930
American Association of Physicists in Medicine
Résumé en anglais
PurposeThe 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 ...Lire la suite >
PurposeThe 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.MethodsThe 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 urn:x-wiley:00942405:media:mp14490:mp14490-math-0001 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.ResultsThe 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 urn:x-wiley:00942405:media:mp14490:mp14490-math-0002 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 urn:x-wiley:00942405:media:mp14490:mp14490-math-0003, 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.ConclusionThe 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.< Réduire
Origine
Importé de halUnités de recherche