Monte Carlo simulations of microdosimetry and radiolytic species production at long time post proton irradiation using GATE and Geant4-DNA.

Abstract

Radiobiological effectiveness of radiation in cancer treatment can be studied at different scales (molecular till organ scale) and different time post irradiation. The production of free radicals and reactive oxygen species during water radiolysis is particularly relevant to understand the fundamental mechanisms playing a role in observed biological outcomes. The development and validation of Monte Carlo tools integrating the simulation of physical, physico-chemical and chemical stages after radiation is very important to maintain with experiments. Therefore, in this study, we propose to validate a new Geant4-DNA chemistry module through the simulation of water radiolysis and Fricke dosimetry experiments on a proton preclinical beam line. In this study, we used the GATE Monte Carlo simulation platform (version 9.3) to simulate a 67.5 MeV proton beam produced with the ARRONAX isochronous cyclotron (IBA Cyclone 70XP) at conventional dose rate (0.2Gy/s) to simulate the irradiation of ultra-pure liquid water samples and Fricke dosimeter. We compared the depth dose profile with measurements performed with a plane parallel Advanced PTW 34045 Markus ionization chamber. Then, a new Geant4-DNA chemistry application proposed from Geant4 version 11.2 has been used to assess the evolution of , , , , , , and reactive species along time until 1-h post-irradiation. In particular, the effect of oxygen and pH has been investigated through comparisons with experimental measurements of radiolytic yields for and Fe3+. GATE simulations reproduced, within 4%, the depth dose profile in liquid water. With Geant4-DNA, we were able to reproduce experimental radiolytic yields 1-h post-irradiation in aerated and deaerated conditions, showing the impact of small changes in oxygen concentrations on species evolution along time. For the Fricke dosimeter, simulated G(Fe3+) is 15.97±0.2 molecules/100eV which is 11% higher than the measured value (14.4±04 molecules/100 eV). These results aim to be consolidated by new comparisons involving other radiolytic species, such as or to further study the mechanisms underlying the FLASH effect observed at ultra-high dose rates (UHDR).