Hydrated Electron Yields in the Heavy Ion Radiolysis of Water

Review of chemical models and applications in Geant4-DNA: Report from the ESA BioRad III Project

A chemistry module has been implemented in Geant4-DNA since Geant4 version 10.1 to simulate the radiolysis of water after irradiation. It has been used in a number of applications, including the calculation of G-values and early DNA damage, allowing the comparison with experimental data. Since the first version, numerous modifications have been made to the module to improve the computational efficiency and extend the simulation to homogeneous kinetics in bulk solution. With these new developments, new applications have been proposed and released as Geant4 examples, showing how to use chemical processes and models. This work reviews the models implemented and application developments for modeling water radiolysis in Geant4-DNA as reported in the ESA BioRad III Project.

Unraveling the Complexity of DNA Radiation Damage Using DNA Nanotechnology

Radiation cancer therapies use different ionizing radiation qualities that damage DNA molecules in tumor cells by a yet not completely understood plethora of mechanisms and processes. While the direct action of the radiation is significant, the byproducts of the water radiolysis, mainly secondary low-energy electrons (LEEs, <20 eV) and reactive oxygen species (ROS), can also efficiently cause DNA damage, in terms of DNA strand breakage or DNA interstrand cross-linking. As a result, these types of DNA damage evolve into mutations hindering DNA replication, leading to cancer cell death. Concomitant chemo-radiotherapy explores the addition of radiosensitizing therapeutics commonly targeting DNA, such as platinum derivatives and halogenated nucleosides, to enhance the harmful effects of ionizing radiation on the DNA molecule. Further complicating the landscape of DNA damage are secondary structures such as G-quadruplexes occurring in telomeric DNA. These structures protect DNA from radiation damage, rendering them as promising targets for new and more selective cancer radiation treatments, rather than targeting linear DNA. However, despite extensive research, there is no single paradigm approach to understanding the mysterious way in which ionizing radiation causes DNA damage. This is due to the multidisciplinary nature of the field of research, which deals with multiple levels of biological organization, from the molecular building blocks of life toward cells and organisms, as well as with complex multiscale radiation-induced effects. Also, intrinsic DNA features, such as DNA topology and specific oligonucleotide sequences, strongly influence its response to damage from ionizing radiation. In this Account, we present our studies focused on the absolute quantification of photon- and low-energy electron-induced DNA damage in strategically selected target DNA sequences. Our methodology involves using DNA origami nanostructures, specifically the Rothemund triangle, as a platform to expose DNA sequences to either low-energy electrons or vacuum-ultraviolet (VUV, <15 eV) photons and subsequent atomic force microscopy (AFM) analysis. Through this approach, the effects of the DNA sequence, incorporation of halogenated radiosensitizers, DNA topology, and the radiation quality on radiation-induced DNA strand breakage have been systematically assessed and correlated with fundamental photon- and electron-driven mechanisms underlying DNA radiation damage. At lower energies, these mechanisms include dissociative electron attachment (DEA), where electrons attach to DNA molecules causing strand breaks, and dissociative photoexcitation of DNA. Additionally, further dissociative processes such as photoionization and electron impact contribute to the complex cascade of DNA damage events induced by ionizing radiation. We expect that emerging DNA origami-based approaches will lead to a paradigm shift in research fields associated with DNA damage and suggest future directions, which can foster the development of technological applications in nanomedicine, e.g., optimized cancer treatments or the molecular design of optimized radiosensitizing therapeutics.

GEANT4-DNA simulation of temperature-dependent and pH-dependent yields of chemical radiolytic species

Objective.GEANT4-DNA can simulate radiation chemical yield (G-value) for radiolytic species such as the hydrated electron (eaq-) with the independent reaction times (IRT) method, however, only at room temperature and neutral pH. This work aims to modify the GEANT4-DNA source code to enable the calculation ofG-values for radiolytic species at different temperatures and pH values.Approach.In the GEANT4-DNA source code, values of chemical parameters such as reaction rate constant, diffusion coefficient, Onsager radius, and water density were replaced by corresponding temperature-dependent polynomials. The initial concentration of hydrogen ion (H+)/hydronium ion (H3O+) was scaled for a desired pH using the relationship pH = -log10[H+]. To validate our modifications, two sets of simulations were performed. (A) A water cube with 1.0 km sides and a pH of 7 was irradiated with an isotropic electron source of 1 MeV. The end time was 1μs. The temperatures varied from 25 °C to 150 °C. (B) The same setup as (A) was used, however, the temperature was set to 25 °C while the pH varied from 5 to 9. The results were compared with published experimental and simulated work.Main results.The IRT method in GEANT4-DNA was successfully modified to simulateG-values for radiolytic species at different temperatures and pH values. Our temperature-dependent results agreed with experimental data within 0.64%-9.79%, and with simulated data within 3.52%-12.47%. The pH-dependent results agreed well with experimental data within 0.52% to 3.19% except at a pH of 5 (15.99%) and with simulated data within 4.40%-5.53%. The uncertainties were below ±0.20%. Overall our results agreed better with experimental than simulation data.Significance.Modifications in the GEANT4-DNA code enabled the calculation ofG-values for radiolytic species at different temperatures and pH values.

CPU-GPU coupling independent reaction times method in NASIC and application in water radiolysis by FLASH irradiation

The mechanism of the FLASH effect remains unclear and could be revealed by studying chemical reactions during irradiation. Monte Carlo simulation of the radiolytic species is an effective tool to analyze chemical reactions, but the simulation is limited by computing costs of the step-by-step simulation of radiolytic species, especially when considering beam with complex time structure. The complexity of the time structure of beams from accelerators in FLASH radiotherapy requires a high-performance Monte Carlo code. In this work, we develop a CPU-GPU coupling accelerating code with the independent reaction times (IRT) method to extend the chemical module of our nanodosimetry Monte Carlo code NASIC. Every chemical molecule in the microenvironment contains time information to consider the reactions from different tracks and simulate beams with complex time structures. Performance test shows that our code significantly improved the computing efficiency of the chemical module by four orders of magnitude. Then the code is used to study the oxygen depletion hypothesis in FLASH radiotherapy for different conditions by setting different parameters. The transient oxygen consumption rate values in the water are calculated when the pulses width ranges from 2 ps to 2 μs, the total dose ranges from 0.5 Gy to 100 Gy and the initial oxygen concentration ranges from 0.1% to 21%. The time evolution curves are simulated to study the effect of the time structure of an electron linear accelerator. Results show that the total dose in several microseconds is a better indicator reflecting the radiolytic oxygen consumption rate than the dose rate. The initial oxygen greatly affects the oxygen consumption rate because of the reaction competition. The diffusion of oxygen determined by the physiological parameters is the key factor affecting oxygen depletion during the radiation using electron linear accelerators. Our code provides an efficient tool for simulating water radiolysis in different conditions

Real-Time Electron Solvation Induced by Bursts of Laser-Accelerated Protons in Liquid Water

Understanding the mechanisms of proton energy deposition in matter and subsequent damage formation is fundamental to radiation science. Here we exploit the picosecond (10^{-12}  s) resolution of laser-driven accelerators to track ultrafast solvation dynamics for electrons due to proton radiolysis in liquid water (H_{2}O). Comparing these results with modeling that assumes initial conditions similar to those found in photolysis reveals that solvation time due to protons is extended by >20  ps. Supported by magnetohydrodynamic theory this indicates a highly dynamic phase in the immediate aftermath of the proton interaction that is not accounted for in current models.

A self-referenced in-situ arrival time monitor for X-ray free-electron lasers

We present a novel, highly versatile, and self-referenced arrival time monitor for measuring the femtosecond time delay between a hard X-ray pulse from a free-electron laser and an optical laser pulse, measured directly on the same sample used for pump-probe experiments. Two chirped and picosecond long optical supercontinuum pulses traverse the sample with a mutually fixed time delay of 970 fs, while a femtosecond X-ray pulse arrives at an instant in between both pulses. Behind the sample the supercontinuum pulses are temporally overlapped to yield near-perfect destructive interference in the absence of the X-ray pulse. Stimulation of the sample with an X-ray pulse delivers non-zero contributions at certain optical wavelengths, which serve as a measure of the relative arrival time of the X-ray pulse with an accuracy of better than 25 fs. We find an excellent agreement of our monitor with the existing timing diagnostics at the SACLA XFEL with a Pearson correlation value of 0.98. We demonstrate a high sensitivity to measure X-ray pulses with pulse energies as low as 30 \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\upmu $$\end{document}μJ. Using a free-flowing liquid jet as interaction sample ensures the full replacement of the sample volume for each X-ray/optical event, thus enabling its utility even at MHz repetition rate XFEL sources.

A review of simulation codes and approaches for radiation chemistry

Historically, the field of radiation chemistry began shortly after the discovery of radioactivity, and its development has been closely related to discoveries in other related fields such as radiation and nuclear physics. Radiolysis of water and radiation chemistry have been very important in elucidating how radiation affects living matter and how it induces DNA damage. Nowadays, we recognize the importance of chemistry to understanding the effects of radiation on cells; however, it took several decades to obtain this insight, and much is still unknown. The radiolysis of water and aqueous solutions have been the subject of much experimental and theoretical research for many decades. One important concept closely related to radiation chemistry is radiation track structure. Track structure results from early physical and physicochemical events that lead to a highly non-homogenous distribution of radiolytic species. Because ionizing radiation creates unstable species that are distributed non-homogenously, the use of conventional reaction kinetics methods does not describe this chemistry well. In recent years, several methods have been developed for simulating radiation chemistry. In this review, we give a brief history of the field and the development of the simulation codes. We review the current methods used to simulate radiolysis of water and radiation chemistry, and we describe several radiation chemistry codes and their applications.

Independent reaction times method in Geant4-DNA: Implementation and performance

PURPOSE

The 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.

METHODS

The 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 e aq - 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.

RESULTS

The 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 e aq - 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 e aq - , 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.

CONCLUSION

The 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.

LET-Dependent Intertrack Yields in Proton Irradiation at Ultra-High Dose Rates Relevant for FLASH Therapy

FLASH radiotherapy delivers a high dose (≥10 Gy) at a high rate (≥40 Gy/s). In this way, particles are delivered in pulses as short as a few nanoseconds. At that rate, intertrack reactions between chemical species produced within the same pulse may affect the heterogeneous chemistry stage of water radiolysis. This stochastic process suits the capabilities of the Monte Carlo method, which can model intertrack effects to aid in radiobiology research, including the design and interpretation of experiments. In this work, the TOPAS-nBio Monte Carlo track-structure code was expanded to allow simulations of intertrack effects in the chemical stage of water radiolysis. Simulation of the behavior of radiolytic yields over a long period of time (up to 50 s) was verified by simulating radiolysis in a Fricke dosimeter irradiated by 60Co γ rays. In addition, LET-dependent G values of protons delivered in single squared pulses of widths, 1 ns, 1 µs and 10 µs, were obtained and compared to simulations using no intertrack considerations. The Fricke simulation for the calculated G value of Fe3+ ion at 50 s was within 0.4% of the accepted value from ICRU Report 34. For LET-dependent G values at the end of the chemical stage, intertrack effects were significant at LET values below 2 keV/µm. Above 2 keV/µm the reaction kinetics remained limited locally within each track and thus, effects of intertrack reactions remained low. Therefore, when track structure simulations are used to investigate the biological damage of FLASH irradiation, these intertrack reactions should be considered. The TOPAS-nBio framework with the expansion to intertrack chemistry simulation provides a useful tool to assist in this task.

Chemical Analysis of Secondary Electron Emission from a Water Cathode at the Interface with a Nonthermal Plasma

When a nonthermal plasma and a liquid form part of the same circuit, the liquid may function as a cathode, in which case electrons are emitted from the liquid into the gas to sustain the plasma. As opposed to solid electrodes, the mechanism of this emission has not been established for a liquid, even though various theories have attempted to explain it via chemical processes in the liquid phase. In this work, we tested the effects of the interfacial chemistry on electron emission from water, including the role of pH as well as the hydroxyl radical, the hydrogen atom, the solvated electron, and the presolvated electron; it was found that none of these species are critical to sustain the plasma. We propose an emission mechanism where electrons, generated from ionized water molecules in the uppermost monolayers of solution, are emitted into the plasma directly from the conduction band of the water.

Dosimetry as a Catch in Radiobiology Experiments

Experimental radiobiological studies in which the effects of ionizing radiation on a biological model are examined often highlight the biological aspects while missing detailed descriptions of the geometry, sample and dosimetric methods used. Such omissions can hinder the reproducibility and comparability of the experimental data. An application based on the Geant4 simulation toolkit was developed to design experiments using a biological solution placed in a microtube. The application was used to demonstrate the influence of the type of microtube, sample volume and energy of a proton source on the dose distribution across the sample, and on the mean dose in the whole sample. The results shown here are for samples represented by liquid water in the 0.4-, 1.5- and 2.0-ml microtubes irradiated with 20, 30 and 100 MeV proton beams. The results of this work demonstrate that the mean dose and homogeneity of the dose distribution within the sample strongly depend on all three parameters. Furthermore, this work shows how the dose uncertainty propagates into the scored primary DNA damages in plasmid DNA studies using agarose gel electrophoresis. This application is provided freely to assist users in verifying their experimental setup prior to the experiment.

Inhibition of Radiolytic Molecular Hydrogen Formation by Quenching of Excited State Water

Comparison of experimental measurements of the yield of molecular hydrogen produced in the gamma radiolysis of water and aqueous nitrate solutions with predictions of a Monte Carlo track chemistry model shows that the nitrate anion scavenging of the hydrated electron, its precursor, and hydrogen atom cannot account for the observed decrease in the yield at high nitrate anion concentrations. Inclusion of the quenching of excited states of water (formed by either direct excitation or reaction of the water radical cation with the precursor to the hydrated electron) by the nitrate anion into the reaction scheme provides excellent agreement between the stochastic calculations and experiment demonstrating the existence of this short-lived species and its importance in water radiolysis. Energy transfer from the excited states of water to the nitrate anion producing an excited state provides an additional pathway for the production of nitrogen containing products not accounted for in traditional radiation chemistry scenarios. Such reactions are of central importance in predicting the behavior of liquors common in the reprocessing of spent nuclear fuel and the storage of highly radioactive liquid waste prior to vitrification.

Efficient Reductive Decomposition of Perfluorooctanesulfonate in a High Photon Flux UV/Sulfite System

Hydrated electron (e_aq^-) induced reduction techniques are promising for decomposing recalcitrant organic pollutants. However, its vigorous reactivity with copresent scavenging species and the difficulty in minimizing the competitive reactions make the proportion of e_aq^- participating in pollutant decomposition low, reflecting by slow decomposition kinetics. In this study, a high photon flux UV/sulfite system was employed to promote e_aq^- production. Its feasibility in enhancing a notorious recalcitrant pollutant, PFOS, decomposition was investigated. The effective photon flux utilized for producing e_aq^- was 9.93 × 10^-8 einstein/cm²·s. At initial solution pH 9.2, with DO about 5 mg/L, and at around 25 °C, 98% PFOS was decomposed within 30 min from its initial concentration of 32 μM. The k_obs of PFOS decomposition was 0.118 min^-1 (7.08 h^-1), and about 8-400 folds faster than those obtained in other reductive approaches. In this system, PFOS decomposition showed can tolerate copresent 7 mg N/L of NO₃^-. Suggested by molecular orbitals and thermodynamic analyses, the mechanisms responsible for PFOS decomposition involve defluorination, desulfonation, and centermost C-C bond scission. By demonstrating a more practical relevant treatment process, the outcomes of this study would be helpful for facilitating future applications of e_aq^- induced reduction techniques for efficient recalcitrant pollutants decomposition.

Track structures, DNA targets and radiation effects in the biophysical Monte Carlo simulation code PARTRAC

This review describes the PARTRAC suite of comprehensive Monte Carlo simulation tools for calculations of track structures of a variety of ionizing radiation qualities and their biological effects. A multi-scale target model characterizes essential structures of the whole genomic DNA within human fibroblasts and lymphocytes in atomic resolution. Calculation methods and essential results are recapitulated regarding the physical, physico-chemical and chemical stage of track structure development of radiation damage induction. Recent model extension towards DNA repair processes extends the time dimension by about 12 orders of magnitude and paves the way for superior predictions of radiation risks.

Characterization of condensed plasmid DNA models for studying the direct effect of ionizing radiation

We have examined the changes in physical properties of aqueous solutions of the plasmid pUC18 that take place on the addition of the cationic oligopeptide penta-arginine. An increase in sedimentation rate and static light scattering, and changes in the nucleic acid CD spectrum all suggest that this ligand acts to condense the plasmid. Dynamic light scattering suggests the hydrodynamic radii of the condensate particles are a few micrometers, ca. 50-fold larger than that of the monomeric plasmid. Condensation of the plasmid also produces a ca. 100-fold decrease in the strand break yield produced by gamma irradiation. This extensive protection against reactive intermediates in the bulk of the solution implies that condensed plasmid DNA may offer a model system with which to study the direct effect of ionizing radiation (ionization of the DNA itself). The use of peptide ligands as condensing agents in this application is attractive because the derivatives of several amino acids (particularly tryptophan and tyrosine) have been shown to modify the radiation chemistry of DNA extensively.

High-LET ion radiolysis of water: oxygen production in tracks

It is known that molecular oxygen is a product of the radiolysis of water with high-linear energy transfer (LET) radiation, a result that is of particular significance in radiobiology and of practical relevance in radiotherapy. In fact, it has been suggested that the radiolytic formation of an oxygenated microenvironment around the tracks of high-LET heavy ions is an important factor in their enhanced biological efficiency in the sense that this may be due to an "oxygen effect" by O(2) produced by these ions in situ. Using Monte Carlo track simulations of pure, deaerated water radiolysis by 4.8 MeV (4)He(2+) (LET approximately 94 keV/microm) and 24 MeV (12)C(6+) (LET approximately 490 keV/microm) ions, including the mechanism of multiple ionization of water, we have calculated the yields and concentrations of O(2) in the tracks of these irradiating ions as a function of time between approximately 10(-12) and 10(-5) s at 25 and 37 degrees C. The track oxygen concentrations obtained compare very well with O(2) concentrations estimated from the "effective" amounts of oxygen that are needed to produce the observed reduction in oxygen enhancement ratio (OER) with LET (assuming this decrease is attributable to the sole radiolytic formation of O(2) in the tracks). For example, for 24 MeV (12)C(6+) ions, the initial track concentration of O(2) is estimated to be more than three orders of magnitude higher than the oxygen levels present in normally oxygenated and hypoxic tumor regions as well as in normal human cells. Such results, which largely plead in favor of the "oxygen in the heavy-ion track" hypothesis, could explain at least in part the greater efficiency of high-LET radiation for cell inactivation (at equal radiation dose).

Aqueous solution of UCl6(2-) in O2 saturated acidic medium: an efficient system to scavenge all primary radicals in spurs produced by irradiation

Absorbance measurements find the yield of the oxidation of U(IV) to be (8.75 +/- 0.05) x 10(-7) mol J(-1) in the (60)Co gamma radiolysis of aqueous solutions containing 4.4 x 10(-3) mol L(-1) U(IV) in the presence of O(2) saturated 2 mol L(-1) Cl(-) at pH = 0. This high value of oxidation yield suggests that all primary radicals formed by water decomposition are scavenged in these solutions. Simulations using a nonhomogeneous stochastic kinetic track model agree with the experimental results and are used to explain the mechanism for scavenging radicals and oxidation of U(IV).

An apparatus for the study of high temperature water radiolysis in a nuclear reactor: calibration of dose in a mixed neutron/gamma radiation field

The cooling water of nuclear reactors undergoes radiolytic decomposition induced by gamma, fast electron, and neutron radiation in the core. To model the process, recombination reaction rates and radiolytic yields for the water radical fragments need to be measured at high temperature and pressure. Yields for the action of neutron radiation are particularly hard to determine independently because of the beta/gamma field also present in any reactor. In this paper we report the design of an apparatus intended to measure neutron radiolysis yields as a function of temperature and pressure. A new methodology for separation of neutron and beta/gamma radiolysis yields in a mixed radiation field is proposed and demonstrated.

Neutron and β/γ Radiolysis of Water up to Supercritical Conditions. 1. β/γ Yields for H2, H• Atom, and Hydrated Electron

Yields for H2, H(.) atom, and hydrated electron production in beta/gamma radiolysis of water have been measured from room temperature up to 400 degrees C on a 250 bar isobar, and also as a function of pressure (density) at 380 and 400 degrees C. Radiolysis was carried out using a beam of 2-3 MeV electrons from a van de Graaff accelerator, and detection was by mass spectrometer analysis of gases sparged from the irradiated water. N2O was used as a specific scavenger for hydrated electrons giving N2 as product. Ethanol-d(6) was used to scavenge H(.) atoms, giving HD as a stable product. It is found that the hydrated electron yield decreases and the H(.) atom yield increases dramatically at lower densities in supercritical water, and the overall escape yield increases. The yield of molecular H2 increases with temperature and does not tend toward zero at low density, indicating that it is formed promptly rather than in spur recombination. A minimum in both the radical and H2 yields is observed around 0.4 kg/dm(3) density in supercritical water.

Reactivity between Biphenyl and Precursor of Solvated Electrons in Tetrahydrofuran Measured by Picosecond Pulse Radiolysis in Near-Ultraviolet, Visible, and Infrared

The initial decrease of solvated electrons in tetrahydrofuran (THF) upon addition of biphenyl was investigated by picosecond pulse radiolysis. Transient absorption spectra derived from the biphenyl radical anion (centered at 408 and 655 nm) and solvated electrons of THF (infrared) were successfully measured in the wavelength region from 400 to 900 nm by the extension of a femtosecond continuum probe light to near-ultraviolet using a second harmonic generation of Ti:sapphire laser and a CaF2 plate. From the analysis of kinetic traces at 1300 nm considering the overlap of primary solvated electrons and partial biphenyl radical anion, C37, which is defined by the solute concentration to reduce the initial yield of solvated electrons to 1/e, was found to be 87 +/- 3 mM. The rate constant of solvated electrons with biphenyl was determined as 5.8 +/- 0.3 x 10(10) M(-1) s(-1). We demonstrate that the kinetic traces at both 408 nm mainly due to biphenyl radical anion and 1300 nm mainly due to solvated electrons are reproduced with high accuracy and consistency by a simple kinetic analysis. Much higher concentrations of biphenyl (up to 2 M) were examined, showing further increase of the initial yield of biphenyl radical anion accompanying a fast decay component. This observation is discussed in terms of geminate ion recombination, scavenging, delayed geminate ion recombination, and direct ionization of biphenyl at high concentration.