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

BACKGROUND: 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. PURPOSE: 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. MATERIAL AND METHODS: 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.2 Gy/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 HO • ${\mathrm{HO}}^ \bullet $ , e aq - ${\mathrm{e}}_{{\mathrm{aq}}}^ - $ , H 3 O + ${{\mathrm{H}}}_3{{\mathrm{O}}}^ + $ , H 2 O 2 ${{\mathrm{H}}}_2{{\mathrm{O}}}_2$ , H 2 ${{\mathrm{H}}}_2$ , HO 2 • ${\mathrm{HO}}_2^ \bullet $ , HO 2 - , O 2 • - ${\mathrm{HO}}_2^ - ,{\mathrm{\ O}}_2^{ \bullet - }$ and HO - ${\mathrm{HO}}^ - $ 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 H 2 O 2 ${{\mathrm{H}}}_2{{\mathrm{O}}}_2$ and Fe3+. RESULTS: GATE simulations reproduced, within 4%, the depth dose profile in liquid water. With Geant4-DNA, we were able to reproduce experimental H 2 O 2 ${{\mathrm{H}}}_2{{\mathrm{O}}}_2$ 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/100 eV which is 11% higher than the measured value (14.4 ± 04 molecules/100 eV). CONCLUSIONS: These results aim to be consolidated by new comparisons involving other radiolytic species, such as e aq - ${\mathrm{e}}_{{\mathrm{aq}}}^ - $ or , O 2 • - $,{\mathrm{\ O}}_2^{ \bullet - }$ to further study the mechanisms underlying the FLASH effect observed at ultra-high dose rates (UHDR).

Proton Irradiations at Ultra-High Dose Rate vs. Conventional Dose Rate: Strong Impact on Hydrogen Peroxide Yield.

During ultra-high dose rate (UHDR) external radiation therapy, healthy tissues appear to be spared while tumor control remains the same compared to conventional dose rate. However, the understanding of radiochemical and biological mechanisms involved are still to be discussed. This study shows how the hydrogen peroxide (H2O2) production, one of the reactive oxygen species (ROS), could be controlled by early heterogenous radiolysis processes in water during UHDR proton-beam irradiations. Pure water was irradiated in the plateau region (track-segment) with 68 MeV protons under conventional (0.2 Gy/s) and several UHDR conditions (40 Gy/s to 60 kGy/s) at the ARRONAX cyclotron. Production of H2O2 was then monitored using the Ghormley triiodide method. New values of GTS(H2O2) were added in conventional dose rate. A substantial decrease in H2O2 production was observed from 0.2 to 1.5 kGy/s with a more dramatic decrease below 100 Gy/ s. At higher dose rate, up to 60 kGy/s, the H2O2 production stayed stable with a mean decrease of 38% ± 4%. This finding, associated to the decrease in the production of hydroxyl radical (•OH) already observed in other studies in similar conditions can be explained by the well-known spur theory in radiation chemistry. Thus, a two-step FLASH-RT mechanism can be envisioned: an early step at the microsecond scale mainly controlled by heterogenous radiolysis, and a second, slower, dominated by O2 depletion and biochemical processes. To validate this hypothesis, more measurements of radiolytic species will soon be performed, including radicals and associated lifetimes.

Technical note: Proton beam dosimetry at ultra-high dose rates (FLASH): Evaluation of GAFchromic™ (EBT3, EBT-XD) and OrthoChromic (OC-1) film performances.

PURPOSE: The ARRONAX cyclotron facility offers the possibility to deliver proton beams from low to ultra-high dose rates (UHDR). As a good control of the dosimetry is a prerequisite of UHDR experimentations, we evaluated in different conditions the usability and the dose rate dependency of several radiochromic films commonly used for dosimetry in radiotherapy. METHODS: We compared the dose rate dependency of three types of radiochromic films: GAFchromic™ EBT3 and GAFchromic™ EBT-XD (Ashland Inc., Wayne, NJ, USA), and OrthoChromic OC-1 (OrthoChrome Inc., Hillsborough, NJ, USA), after proton irradiations at various mean dose rates (0.25, 40, 1500, and 7500 Gy/s) and for 10 doses (2-130 Gy). We also evaluated the dose rate dependency of each film considering beam structures, from single pulse to multiple pulses with various frequencies. RESULTS: EBT3 and EBT-XD films showed differences of response between conventional (0.25 Gy/s) and UHDR (7500 Gy/s) conditions, above 10 Gy. On the contrary, OC-1 films did not present overall difference of response for doses except below 3 Gy. We observed an increase of the netOD with the mean dose rate for EBT3 and EBT-XD films. OC-1 films did not show any impact of the mean dose rate up to 7500 Gy/s, above 3 Gy. No difference was found based on the beam structure, for all three types of films. CONCLUSIONS: EBT3 and EBT-XD radiochromic films should be used with caution for the dosimetry of UHDR proton beams over 10 Gy. Their overresponse, which increases with mean dose rate and dose, could lead to non-negligible overestimations of the absolute dose. OC-1 films are dose rate independent up to 7500 Gy/s in proton beams. Films response is not impacted by the beam structure. A broader investigation of the usability of OC-1 films in UHDR conditions should be conducted at intermediate and higher mean dose rates and other beam energies.

A Monte Carlo Determination of Dose and Range Uncertainties for Preclinical Studies with a Proton Beam.

Proton therapy (PRT) is an irradiation technique that aims at limiting normal tissue damage while maintaining the tumor response. To study its specificities, the ARRONAX cyclotron is currently developing a preclinical structure compatible with biological experiments. A prerequisite is to identify and control uncertainties on the ARRONAX beamline, which can lead to significant biases in the observed biological results and dose-response relationships, as for any facility. This paper summarizes and quantifies the impact of uncertainty on proton range, absorbed dose, and dose homogeneity in a preclinical context of cell or small animal irradiation on the Bragg curve, using Monte Carlo simulations. All possible sources of uncertainty were investigated and discussed independently. Those with a significant impact were identified, and protocols were established to reduce their consequences. Overall, the uncertainties evaluated were similar to those from clinical practice and are considered compatible with the performance of radiobiological experiments, as well as the study of dose-response relationships on this proton beam. Another conclusion of this study is that Monte Carlo simulations can be used to help build preclinical lines in other setups.

THE RADIOBIOLOGICAL PLATFORM AT ARRONAX.

The cyclotron ARRONAX can deliver different types of particles (protons, deuterons, alpha-particles) in an energy range up to 68 MeV. One of its six experimental halls is dedicated to studying the interactions of radiation with matter including living matter. A horizontal beamline for cell irradiation has been setup and characterized. The radiobiological characterization was done in terms of V79 cells survival after irradiation with 68 MeV protons. The results demonstrate that radiobiological studies can be successfully performed confirming the high potential of the facility.