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.

Picosecond Pulse Radiolysis of Highly Concentrated Carbonate Solutions.

Highly concentrated potassium carbonate aqueous solutions are studied by picosecond pulse radiolysis with the purpose of exploring the formation processes of carbonate radical CO3(•-). The transient absorption band of solvated electron produced by ionizing is markedly shifted from 715 to 600 nm when the solute concentration of K2CO3 is 5 mol L(-1). This spectral shift is even more important than that observed for the solvated electron in 10 mol L(-1) KOH solutions. The broad absorption band of solvated electron in K2CO3 solutions overlaps with that of carbonate radical CO3(•-) formed at ultrashort time. Nitrate ion is used to scavenge the solvated electron and to observe the contribution of carbonate radical CO3(•-). The analysis of the amplitude and the kinetics of carbonate radical formation in highly concentrated solutions shows that CO3(•-) is formed within the electron pulse (7 ps) by two parallel mechanisms: a direct effect on the solute and the oxidation of the solute by water radical hole H2O(•+). These two mechanisms are followed by an additional one, by reaction between the solute and OH(•) radical especially in lower concentration. The radiolytic yield of each process is discussed.

Discrepancies in thorium oxide solubility values: study of attachment/detachment processes at the solid/solution interface.

The solubility of thorium under oxide and/or hydroxide forms has been extensively studied for many years. Nevertheless, a large discrepancy in the solubility values is noticed in the literature. We study Th atom exchange between thorium oxide surfaces and various aqueous solutions (0.01 mol·L(-1) NaCl for 0.0 < pH < 5.2) to address this issue. By solid-state characterization [X-ray photoelectron spectroscopy (XPS), scanning electron microscopy, and atomic force microscopy], we determined that 80% of the XPS accessible near the surface region of sintered thorium oxide is represented by the less reactive ThO(2)(cr) grains. The remaining 20% corresponds to ThO(x)(OH)(y)(H(2)O)(z), which is largely associated with grain boundaries. Only the latter fraction is involved in solid/solution exchange mechanisms. Local conditions (thorium concentrations, pH values, etc.) in grain boundaries lead to an adjustment of the "local solubility constraints" and explain the thorium concentration measured in our experiments. For pH <5.2, the thorium concentration and pH gradient between the bulk solution and grain-boundary regions imply that the solubility values mainly depend on the availability and accessibility of ThO(x)(OH)(y)(H(2)O)(z). We have performed two solubility experiments with a (232)ThO(2)(cr) solid in a 0.01 mol·L(-1) NaCl solution for 300 days. In a first experiment, we measured (232)Th concentrations in dissolution experiments in order to determine the detachment rates of Th atoms from the solid surface. In a subsequent step, we added (229)Th to the solution in order to measure the surface attachment rate for dissolved Th atoms. This allowed an assessment of the net balance of Th atom exchange at the solid/solution interface. The empirical solubility data do not correspond to the thermodynamic bulk phase/solution equilibrium because measured solution concentrations are controlled by site-specific exchange mechanisms at the solid/solution interface. Therefore, for sparingly soluble solids, one needs to quantify site-specific surface attachment and detachment rates if one wants to assess solubility constraints.

Interaction mechanisms between uranium(VI) and rutile titanium dioxide: from single crystal to powder.

This paper is devoted to the study of the mechanisms of interaction between uranyl ion and rutile TiO2. Among the radionuclides of interest, U(VI) can be considered as a model of the radionuclides oxo-cations. The substrate under study here is the rutile titanium dioxide (TiO2) which is an interesting candidate as a methodological solid since it can be easily found as powder and as manufactured single crystals. This material presents also a wide domain of stability as a function of pH. Then, it allows the study of the retention processes on well-defined crystallographic planes, which can lead to a better understanding of the surface reaction mechanisms. Moreover, it is well-established that the (110) crystallographic orientation is dominating the surface chemistry of the rutile powder. Therefore, the spectroscopic results obtained for the U(VI)/rutile (110) system and other relevant crystallographic orientations were used to have some insight on the nature of the uranium surface complexes formed on rutile powder. This goal was achieved by using time-resolved laser-induced fluorescence spectroscopy (TRLFS) which allows the investigation, at a molecular scale, of the nature of the reactive surface sites as well as the surface species. For rutile surfaces, oxygen atoms can be 3-fold, 2-fold (bridging oxygens), or single-fold (top oxygens) coordinated to titanium atoms. However, among these three types of surface oxygen atoms, the 3-fold coordinated ones are not reactive toward water molecules or aqueous metallic cations. This study led to conclude on the presence of two uranium(VI) surface complexes: the first one corresponds to the sorption of aquo UO22+ ion sorbed on two bridging oxygen atoms, while the second one, which is favored at higher surface coverages, corresponds to the retention of UO22+ by one bridging and one top oxygen atom. Thus, the approach presented in this paper allows the establishment of experimental constraints that have to be taken into account in the modeling of the sorption mechanisms.

Sorption of uranyl cations on a rutile (001) single crystal monitored by surface second-harmonic generation.

The rotational anisotropy of second-harmonic generation at the surface of a (001) single-crystal rutile is obtained in the presence of uranyl cations sorbed at the surface from acidic solutions at various concentrations. Surface second-harmonic generation appears to be sensitive to the presence of uranyl cations on the rutile samples. Evolution of the anisotropy pattern with initial uranyl concentration is analyzed through a phenomenological model. The elements obtained for the nonlinear susceptibility tensor Chi(2) for each sample significantly constrain the geometry of the possible sorption complexes between uranyl cations and rutile and lead to the proposition of two sorption sites involving different oxygen atoms of the rutile surface.