Electron Emission from Foils and Biological Materials after Proton Impact

Ionizing Radiation and Complex DNA Damage: Quantifying the Radiobiological Damage Using Monte Carlo Simulations

Ionizing radiation is a common tool in medical procedures. Monte Carlo (MC) techniques are widely used when dosimetry is the matter of investigation. The scientific community has invested, over the last 20 years, a lot of effort into improving the knowledge of radiation biology. The present article aims to summarize the understanding of the field of DNA damage response (DDR) to ionizing radiation by providing an overview on MC simulation studies that try to explain several aspects of radiation biology. The need for accurate techniques for the quantification of DNA damage is crucial, as it becomes a clinical need to evaluate the outcome of various applications including both low- and high-energy radiation medical procedures. Understanding DNA repair processes would improve radiation therapy procedures. Monte Carlo simulations are a promising tool in radiobiology studies, as there are clear prospects for more advanced tools that could be used in multidisciplinary studies, in the fields of physics, medicine, biology and chemistry. Still, lot of effort is needed to evolve MC simulation tools and apply them in multiscale studies starting from small DNA segments and reaching a population of cells.

Ionizing Radiation and Complex DNA Damage: From Prediction to Detection Challenges and Biological Significance

Biological responses to ionizing radiation (IR) have been studied for many years, generally showing the dependence of these responses on the quality of radiation, i.e., the radiation particle type and energy, types of DNA damage, dose and dose rate, type of cells, etc. There is accumulating evidence on the pivotal role of complex (clustered) DNA damage towards the determination of the final biological or even clinical outcome after exposure to IR. In this review, we provide literature evidence about the significant role of damage clustering and advancements that have been made through the years in its detection and prediction using Monte Carlo (MC) simulations. We conclude that in the future, emphasis should be given to a better understanding of the mechanistic links between the induction of complex DNA damage, its processing, and systemic effects at the organism level, like genomic instability and immune responses.

Backscattered electron emission after proton impact on gold nanoparticles with and without polymer shell coating

This work aims at measuring experimentally proton induced secondary electron energy spectra after interaction with gold nano particles (GNPs) and polymer-coated GNPs. Backscattered electron energy spectra were collected over a 0 to 1000 eV energy range using a retarding field analyzer (RFA). This paper presents the spectra obtained for proton beam energies of 0.5 and 2 MeV and diameter 2.5 and 3.8 nm GNPs. The spectra were also measured for 3.8 nm GNPs after 5 and 10 MeV proton irradiations. GNPs were deposited on a 100 nm carbon film. Each experimental spectrum was compared with dedicated simulations based on existing numerical models used in the TRAX and Geant4 Monte Carlo codes. For 100 nm carbon target, good agreement between experimental, TRAX and Geant4 simulation results can be observed. For 3.8 nm GNPs, the TRAX simulations reproduce with good agreement the electron energy spectra produced after 0.5, 2, 5 and 10 MeV proton irradiations, while Geant4 spectra display a lower secondary electron yield at low energy (<600 eV) for all the studied energies. This underestimation can mostly be explained by the 790 eV threshold applied in the condensed history model used by Geant4 which impacts the secondary electron energy distribution. Results obtained for carbon and gold targets highlight the impact of the secondary electron production threshold for proton ionization process considered in condensed history models. The experimental results demonstrate that the single interaction approach used in TRAX is adapted to reproduce secondary electron emission from GNPs. On the other hand, the standard electron generation threshold implement in G4BetheBlochModel and G4BraggModel condensed-history models used in Geant4 is not adapted to reproduce low energy electron emission in gold targets. Finally, the results highlight that the GNP coating leads to a decrease of the electron yield and mostly affects low energy electrons (<500 eV) emitted from GNPs.

Energy Spectra of Protons and Generated Secondary Electrons around the Bragg Peak in Materials of Interest in Proton Therapy

The number and energy of secondary electrons generated around the trajectories of swift protons interacting with biological materials are highly relevant in proton therapy, due to the prominent role of low-energy electrons in the production of biodamage. For a given material, electron energy distributions are determined by the proton energy; and it is imperative that the distribution of proton energy at depths around the Bragg peak region be described as accurately as possible. With this objective, we simulated the energy distributions of proton beams of clinically relevant energies (50-300 MeV) at depths around the Bragg peak in liquid water and the water-equivalent polymer poly(methyl methacrylate) (PMMA). By using a simple model, this simulation has been conveniently extended to account for nuclear fragmentation reactions, providing depth-dose curves in excellent agreement with available experimental data. Special care has been taken to describe the electronic excitation spectrum of the target, taking into account its condensed phase nature. A predictive formula has been obtained for the mean value and the width of the proton energy distribution at the Bragg peak depth, quantities which are found to grow linearly with the initial energy of the beam, in good agreement with available data. To accurately characterize (in number and energy) the electrons generated around the proton paths, the energy distributions of the latter at each depth have been convoluted with the energy-dependent ionization inverse mean free paths. This results in a number of low-energy electrons around the Bragg peak larger than when only the proton beam average energy at the given depths is considered. The convoluted ionization inverse mean free path closely resembles the Bragg curve shape. The average energy of the secondary electrons is nearly constant (∼55 eV for liquid water and ∼43 eV for PMMA) in the plateau of the Bragg curve, independent of the proton incident energy and suddenly decaying once the Bragg peak is reached. These findings highlight the importance of a precise calculation of the proton beam energy distribution as a function of the target depth to reliably characterize the secondary electrons generated around the Bragg peak region.

Cross sections for track structure codes: volume versus surface transport

Cross-section calculations and transport models for Monte Carlo track structure codes are discussed as well as the simulation of secondary electron emission yields from thin metal foils. Inelastic cross sections for volume (bulk) and surface transport of electrons in copper are presented and implemented into PARTRAC. Simulations for the volume and surface excitation model are presented and analysed.

Accurate Measurements of Dielectric and Optical Functions of Liquid Water and Liquid Benzene in the VUV Region (1-100 eV) Using Small-Angle Inelastic X-ray Scattering

Using a third-generation synchrotron source (the BL12XU beamline at SPring-8), inelastic X-ray scattering (IXS) spectra of liquid water and liquid benzene were measured at energy losses of 1-100 eV with 0.24 eV resolution for small momentum transfers (q) of 0.23 and 0.32 au with ±0.06 au uncertainty for q. For both liquids, the IXS profiles at these values of q converged well after we corrected for multiple scattering, and these results confirmed the dipole approximation for q ≤ ∼0.3 au. Several dielectric and optical functions [including the optical oscillator strength distribution (OOS), the optical energy-loss function (OLF), the complex dielectric function, the complex index of refraction, and the reflectance] in the vacuum ultraviolet region were derived and tabulated from these small-angle (small q) IXS spectra. These new data were compared with previously obtained results, and this comparison demonstrated the strong reproducibility and accuracy of IXS spectroscopy. For both water and benzene, there was a notable similarity between the OOSs of the liquids and amorphous solids, and there was no evidence of plasmon excitation in the OLF. The static structure factor [S(q)] for q ≤ ∼0.3 au was also deduced and suggests that molecular models that include electron correlation effects can serve as a good approximation for the liquid S(q) values over the full range of q.

Track-structure simulations for charged particles

Monte Carlo track-structure simulations provide a detailed and accurate picture of radiation transport of charged particles through condensed matter of biological interest. Liquid water serves as a surrogate for soft tissue and is used in most Monte Carlo track-structure codes. Basic theories of radiation transport and track-structure simulations are discussed and differences compared to condensed history codes highlighted. Interaction cross sections for electrons, protons, alpha particles, and light and heavy ions are required input data for track-structure simulations. Different calculation methods, including the plane-wave Born approximation, the dielectric theory, and semi-empirical approaches are presented using liquid water as a target. Low-energy electron transport and light ion transport are discussed as areas of special interest.

Track structure, radiation quality and initial radiobiological events: considerations based on the PARTRAC code experience

PURPOSE

The role of track structures for understanding the biological effects of radiation has been the subject of research activities for decades. The physics that describes such processes is the core Monte Carlo codes, such as the biophysical PARTRAC (PARticle TRACks) code described in this review, which follow the mechanisms of radiation-matter interaction from the early stage. In this paper a review of the track structure theory (and of its possible extension concerning non-DNA targets) is presented.

MATERIALS AND METHODS

The role of radiation quality and track structure is analyzed starting from the heavy ions results obtained with the biophysical Monte Carlo code PARTRAC (PARticles TRACks). PARTRAC calculates DNA damage in human cells based on the superposition of simulated track structures in liquid water to an 'atom-by-atom' model of human DNA.

RESULTS

Calculations for DNA fragmentation compared with experimental data for different radiation qualities are illustrated. As an example, the strong dependence of the complexity of DNA damage on radiation track structure, and the very large production of very small DNA fragments (lower than 1 kbp (kilo base pairs) usually not detected experimentally) after high LET (high-Linear Energy Transfer) irradiation is shown. Furthermore the possible importance of non-nuclear/non-DNA targets is discussed in the particular case of cellular membrane and mitochondria.

CONCLUSIONS

The importance of the track structure is underlined, in particular the dependence of a given late cellular effect on the spatial distribution of DNA double-strand breaks (DSB) along the radiation track. These results show that the relative biological effectiveness (RBE) for DSB production can be significantly larger than 1. Moreover the cluster properties of high LET radiation may determine specific initial targets and damage evolution.

Simulation of secondary electron yields from thin metal foils after fast proton impact

Differential and total inelastic cross sections are derived for the interaction between fast protons and Cu. The calculations are done under the non-relativistic plane-wave first-Born approximation and the dielectric theory. A semi-empirical optical oscillator strength density function and a simple linear-momentum dispersion algorithm are used to construct the energy loss function or Bethe surface of the medium. A transport model using these inelastic cross sections is implemented in the Monte Carlo code PARTRAC to simulate the spectra of secondary electron emissions from this homogeneous and isotropic thin copper foil target. Comparisons with experimental results show general agreement for impact energies >50 eV up to non-relativistic values. The model, however, overestimates the secondary electron yields at lower energies.

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.

Electron emission from amorphous solid water induced by passage of energetic protons and fluorine ions

Absolute doubly differential electron emission yields were measured from thin films of amorphous solid water (ASW) after the transmission of 6 MeV protons and 19 MeV (1 MeV/nucleon) fluorine ions. The ASW films were frozen on thin (1-microm) copper foils cooled to approximately 50 K. Electrons emitted from the films were detected as a function of angle in both the forward and backward direction and as a function of the film thickness. Electron energies were determined by measuring the ejected electron time of flight, a technique that optimizes the accuracy of measuring low-energy electron yields, where the effects of molecular environment on electron transport are expected to be most evident. Relative electron emission yields were normalized to an absolute scale by comparison of the integrated total yields for proton-induced electron emission from the copper substrate to values published previously. The absolute doubly differential yields from ASW are presented along with integrated values, providing single differential and total electron emission yields. These data may provide benchmark tests of Monte Carlo track structure codes commonly used for assessing the effects of radiation quality on biological effectiveness.