TY - JOUR
T1 - Modeling the effect of oxygen on the chemical stage of water radiolysis using GPU-based microscopic Monte Carlo simulations, with an application in FLASH radiotherapy
AU - Lai, Youfang
AU - Jia, Xun
AU - Chi, Yujie
N1 - Publisher Copyright:
© 2021 Institute of Physics and Engineering in Medicine.
PY - 2021/1/21
Y1 - 2021/1/21
N2 - Oxygen plays a critical role in determining the initial DNA damages induced by ionizing radiation. It is important to mechanistically model the oxygen effect in the water radiolysis process. However, due to the computational costs from the many body interaction problem, oxygen is often ignored or treated as a constant continuum radiolysis-scavenger background in the simulations using common microscopic Monte Carlo tools. In this work, we reported our recent progress on the modeling of the chemical stage of the water radiolysis with an explicit consideration of the oxygen effect, based upon our initial development of an open-source graphical processing unit (GPU)-based MC simulation tool, gMicroMC. The inclusion of oxygen mainly reduces the yields of eh and H. chemical radicals, turning them into highly toxic O2 .- and HO. 2 species. To demonstrate the practical value of gMicroMC in large scale simulation problems, we applied the oxygen-simulation-enabled gMicroMC to compute the yields of chemical radicals under a high instantaneous dose rate Di to study the oxygen depletion hypothesis in FLASH radiotherapy. A decreased oxygen consumption rate (OCR) was found associated with a reduced initial oxygen concentration level due to reduced probabilities of reactions. With respect to dose rate, for the oxygen concentration of 21% and electron energy of 4.5 keV, OCR remained approximately constant (∼0.22 μM Gy-1) for Di's of 106, 107 Gy s-1 and reduced to 0.19 μM Gy-1 at 108 Gy s-1, because the increased dose rate improved the mutual reaction frequencies among radicals, hence reducing their reactions with oxygen. We computed the time evolution of oxygen concentration under the FLASH irradiation setups. At the dose rate of 107 Gy s-1 and initial oxygen concentrations from 0.01% to 21%, the oxygen is unlikely to be fully depleted with an accumulative dose of 30 Gy, which is a typical dose used in FLASH experiments. The computational efficiency of gMicroMC when considering oxygen molecules in the chemical stage was evaluated through benchmark work to GEANT4-DNA with simulating an equivalent number of radicals. With an initial oxygen concentration of 3% (∼105 molecules), a speedup factor of 1228 was achieved for gMicroMC on a single GPU card when comparing with GEANT4-DNA on a single CPU.
AB - Oxygen plays a critical role in determining the initial DNA damages induced by ionizing radiation. It is important to mechanistically model the oxygen effect in the water radiolysis process. However, due to the computational costs from the many body interaction problem, oxygen is often ignored or treated as a constant continuum radiolysis-scavenger background in the simulations using common microscopic Monte Carlo tools. In this work, we reported our recent progress on the modeling of the chemical stage of the water radiolysis with an explicit consideration of the oxygen effect, based upon our initial development of an open-source graphical processing unit (GPU)-based MC simulation tool, gMicroMC. The inclusion of oxygen mainly reduces the yields of eh and H. chemical radicals, turning them into highly toxic O2 .- and HO. 2 species. To demonstrate the practical value of gMicroMC in large scale simulation problems, we applied the oxygen-simulation-enabled gMicroMC to compute the yields of chemical radicals under a high instantaneous dose rate Di to study the oxygen depletion hypothesis in FLASH radiotherapy. A decreased oxygen consumption rate (OCR) was found associated with a reduced initial oxygen concentration level due to reduced probabilities of reactions. With respect to dose rate, for the oxygen concentration of 21% and electron energy of 4.5 keV, OCR remained approximately constant (∼0.22 μM Gy-1) for Di's of 106, 107 Gy s-1 and reduced to 0.19 μM Gy-1 at 108 Gy s-1, because the increased dose rate improved the mutual reaction frequencies among radicals, hence reducing their reactions with oxygen. We computed the time evolution of oxygen concentration under the FLASH irradiation setups. At the dose rate of 107 Gy s-1 and initial oxygen concentrations from 0.01% to 21%, the oxygen is unlikely to be fully depleted with an accumulative dose of 30 Gy, which is a typical dose used in FLASH experiments. The computational efficiency of gMicroMC when considering oxygen molecules in the chemical stage was evaluated through benchmark work to GEANT4-DNA with simulating an equivalent number of radicals. With an initial oxygen concentration of 3% (∼105 molecules), a speedup factor of 1228 was achieved for gMicroMC on a single GPU card when comparing with GEANT4-DNA on a single CPU.
KW - FLASH radiotherapy
KW - GPU parallel computation
KW - Microscopic Monte Carlo simulation
KW - Oxygen depletion effect
UR - http://www.scopus.com/inward/record.url?scp=85100286854&partnerID=8YFLogxK
UR - http://www.scopus.com/inward/citedby.url?scp=85100286854&partnerID=8YFLogxK
U2 - 10.1088/1361-6560/abc93b
DO - 10.1088/1361-6560/abc93b
M3 - Article
C2 - 33171449
AN - SCOPUS:85100286854
SN - 0031-9155
VL - 66
JO - Physics in medicine and biology
JF - Physics in medicine and biology
IS - 2
M1 - 025004
ER -