Angular and energy distributions of electrons produced in arbitrary biomaterials by proton impact

Simulation of depth-dose curves and water equivalent ratios of energetic proton beams in cortical bone

We have determined the depth-dose curve, penetration range, and water equivalent ratio for proton beams of clinical energies in cortical bone by means of a detailed and accurate simulation that combines molecular dynamics and Monte Carlo techniques. The fundamental input quantities (stopping power and energy loss straggling) for the simulation were obtained from a reliable electronic excitation spectrum of the condensed-phase target, which takes into account the organic and mineral phases that form it. Our simulations with these inputs, which are in excellent agreement with the scarce data available for a cortical bone target, deviate from simulations performed using other stopping quantities, such as those provided by the International Commission on Radiation Units and Measurements in its widely used Report No. 49 [M. J. Berger et al., Stopping powers and ranges for protons and alpha particles, International Commission on Radiation Units and Measurements, Bethesda, Maryland, 1993]. The results of this paper emphasize the importance of an accurate determination of the stopping quantities of cortical bone to advance towards the millimetric precision for the proton penetration ranges and deposited dose needed in radiotherapy.

Energy Deposition around Swift Carbon-Ion Tracks in Liquid Water

Energetic carbon ions are promising projectiles used for cancer radiotherapy. A thorough knowledge of how the energy of these ions is deposited in biological media (mainly composed of liquid water) is required. This can be attained by means of detailed computer simulations, both macroscopically (relevant for appropriately delivering the dose) and at the nanoscale (important for determining the inflicted radiobiological damage). The energy lost per unit path length (i.e., the so-called stopping power) of carbon ions is here theoretically calculated within the dielectric formalism from the excitation spectrum of liquid water obtained from two complementary approaches (one relying on an optical-data model and the other exclusively on ab initio calculations). In addition, the energy carried at the nanometre scale by the generated secondary electrons around the ion's path is simulated by means of a detailed Monte Carlo code. For this purpose, we use the ion and electron cross sections calculated by means of state-of-the art approaches suited to take into account the condensed-phase nature of the liquid water target. As a result of these simulations, the radial dose around the ion's path is obtained, as well as the distributions of clustered events in nanometric volumes similar to the dimensions of DNA convolutions, contributing to the biological damage for carbon ions in a wide energy range, covering from the plateau to the maximum of the Bragg peak.

Relative Role of Physical Mechanisms on Complex Biodamage Induced by Carbon Irradiation

The effective use of swift ion beams in cancer treatment (known as hadrontherapy) as well as appropriate protection in manned space missions rely on the accurate understanding of the energy delivery to cells that damages their genetic information. The key ingredient characterizing the response of a medium to the perturbation induced by charged particles is its electronic excitation spectrum. By using linear-response time-dependent density functional theory, we obtained the energy and momentum transfer excitation spectrum (the energy-loss function, ELF) of liquid water (the main constituent of biological tissues), which was in excellent agreement with experimental data. The inelastic scattering cross sections obtained from this ELF, together with the elastic scattering cross sections derived by considering the condensed phase nature of the medium, were used to perform accurate Monte Carlo simulations of the energy deposited by swift carbon ions in liquid water and carried away by the generated secondary electrons, producing inelastic events such as ionization, excitation, and dissociative electron attachment (DEA). The latter are strongly correlated with cellular death, which is scored in sensitive volumes with the size of two DNA convolutions. The sizes of the clusters of damaging events for a wide range of carbon-ion energies, from those relevant to hadrontherapy up to those for cosmic radiation, predict with unprecedented statistical accuracy the nature and relative magnitude of the main inelastic processes contributing to radiation biodamage, confirming that ionization accounts for the vast majority of complex damage. DEA, typically regarded as a very relevant biodamage mechanism, surprisingly plays a minor role in carbon-ion induced clusters of harmful events.

Excitation and ionisation cross-sections in condensed-phase biomaterials by electrons down to very low energy: application to liquid water and genetic building blocks

A model is presented for computing electron-impact electronic excitation and ionisation cross-sections for arbitrary condensed-phase biomaterials in a wide energy range, showing a general good agreement with the available experimental data.

Multiscale simulation of the focused electron beam induced deposition process

Focused electron beam induced deposition (FEBID) is a powerful technique for 3D-printing of complex nanodevices. However, for resolutions below 10 nm, it struggles to control size, morphology and composition of the structures, due to a lack of molecular-level understanding of the underlying irradiation-driven chemistry (IDC). Computational modeling is a tool to comprehend and further optimize FEBID-related technologies. Here we utilize a novel multiscale methodology which couples Monte Carlo simulations for radiation transport with irradiation-driven molecular dynamics for simulating IDC with atomistic resolution. Through an in depth analysis of [Formula: see text] deposition on [Formula: see text] and its subsequent irradiation with electrons, we provide a comprehensive description of the FEBID process and its intrinsic operation. Our analysis reveals that simulations deliver unprecedented results in modeling the FEBID process, demonstrating an excellent agreement with available experimental data of the simulated nanomaterial composition, microstructure and growth rate as a function of the primary beam parameters. The generality of the methodology provides a powerful tool to study versatile problems where IDC and multiscale phenomena play an essential role.

Efficient ab initio calculation of electronic stopping in disordered systems via geometry pre-sampling: Application to liquid water

Knowledge of the electronic stopping curve for swift ions, S_e(v), particularly around the Bragg peak, is important for understanding radiation damage. Experimentally, however, the determination of such a feature for light ions is very challenging, especially in disordered systems such as liquid water and biological tissue. Recent developments in real-time time-dependent density functional theory (rt-TDDFT) have enabled the calculation of S_e(v) along nm-sized trajectories. However, it is still a challenge to obtain a meaningful statistically averaged S_e(v) that can be compared to observations. In this work, taking advantage of the correlation between the local electronic structure probed by the projectile and the distance from the projectile to the atoms in the target, we devise a trajectory pre-sampling scheme to select, geometrically, a small set of short trajectories to accelerate the convergence of the averaged S_e(v) computed via rt-TDDFT. For protons in liquid water, we first calculate the reference probability distribution function (PDF) for the distance from the proton to the closest oxygen atom, ϕ_R(r_p→O), for a trajectory of a length similar to those sampled experimentally. Then, short trajectories are sequentially selected so that the accumulated PDF reproduces ϕ_R(r_p→O) to increasingly high accuracy. Using these pre-sampled trajectories, we demonstrate that the averaged S_e(v_p) converges in the whole velocity range with less than eight trajectories, while other averaging methods using randomly and uniformly distributed trajectories require approximately ten times the computational effort. This allows us to compare the S_e(v_p) curve to experimental data and assess widely used empirical tables based on Bragg's rule.

Proton transport modeling in a realistic biological environment by using TILDA-V

Whether it is in radiobiology to identify DNA lesions or in medicine to adapt the radiotherapeutic protocols, a detailed understanding of the radiation-induced interactions in living matter is required. Monte Carlo track-structure codes have been successfully developed to describe these interactions and predict the radiation-induced energy deposits at the nanoscale level in the medium of interest. In this work, the quantum-mechanically based Monte Carlo track-structure code TILDA-V has been used to compute the slowing-down of protons in water and DNA. Stopping power and range are then reported and compared with existing data. Then, a first application of TILDA-V to cellular irradiations is also reported in order to highlight the absolute necessity of taking into account a realistic description of the cellular environment in microdosimetry.

Simple and universal model for electron-impact ionization of complex biomolecules

We present a simple and universal approach to calculate the total ionization cross section (TICS) for electron impact ionization in DNA bases and other biomaterials in the condensed phase. Evaluating the electron impact TICS plays a vital role in ion-beam radiobiology simulation at the cellular level, as secondary electrons are the main cause of DNA damage in particle cancer therapy. Our method is based on extending the dielectric formalism. The calculated results agree well with experimental data and show a good comparison with other theoretical calculations. This method only requires information of the chemical composition and density and an estimate of the mean binding energy to produce reasonably accurate TICS of complex biomolecules. Because of its simplicity and great predictive effectiveness, this method could be helpful in situations where the experimental TICS data are absent or scarce, such as in particle cancer therapy.

Differential electron emission from polycyclic aromatic hydrocarbon molecules under fast ion impact

Interaction between polycyclic aromatic hydrocarbon (PAH) molecule and energetic ion is a subject of interest in different areas of modern physics. Here, we present measurements of energy and angular distributions of absolute double differential electron emission cross section for coronene (C₂₄H₁₂) and fluorene (C₁₃H₁₀) molecules under fast bare oxygen ion impact. For coronene, the angular distributions of the low energy electrons are quite different from that of simpler targets like Ne or CH₄, which is not the case for fluorene. The behaviour of the higher electron energy distributions for both the targets are similar to that for simple targets. In case of coronene, a clear signature of plasmon resonance is observed in the analysis of forward-backward angular asymmetry of low energy electron emission. For fluorene, such signature is not identified probably due to lower oscillator strength of plasmon compared to the coronene. The theoretical calculation based on the first-order Born approximation with correct boundary conditions (CB1), in general, reproduced the experimental observations qualitatively, for both the molecules, except in the low energy region for coronene, which again indicates the role of collective excitation. Single differential and total cross sections are also deduced. An overall comparative study is presented.

Double differential cross sections for proton induced electron emission from molecular analogues of DNA constituents for energies in the Bragg peak region

For track structure simulations in the Bragg peak region, measured electron emission cross sections of DNA constituents are required as input for developing parameterized model functions representing the scattering probabilities. In the present work, double differential cross sections were measured for the electron emission from vapor-phase pyrimidine, tetrahydrofuran, and trimethyl phosphate that are structural analogues to the base, the sugar, and the phosphate residue of the DNA, respectively. The range of proton energies was from 75 keV to 135 keV, the angles ranged from 15° to 135°, and the electron energies were measured from 10 eV to 200 eV. Single differential and total electron emission cross sections are derived by integration over angle and electron energy and compared to the semi-empirical Hansen-Kocbach-Stolterfoht (HKS) model and a quantum mechanical calculation employing the first Born approximation with corrected boundary conditions (CB1). The CB1 provides the best prediction of double and single differential cross section, while total cross sections can be fitted with semi-empirical models. The cross sections of the three samples are proportional to their total number of valence electrons.

Water versus DNA: new insights into proton track-structure modelling in radiobiology and radiotherapy

Water is a common surrogate of DNA for modelling the charged particle-induced ionizing processes in living tissue exposed to radiations. The present study aims at scrutinizing the validity of this approximation and then revealing new insights into proton-induced energy transfers by a comparative analysis between water and realistic biological medium. In this context, a self-consistent quantum mechanical modelling of the ionization and electron capture processes is reported within the continuum distorted wave-eikonal initial state framework for both isolated water molecules and DNA components impacted by proton beams. Their respective probability of occurrence-expressed in terms of total cross sections-as well as their energetic signature (potential and kinetic) are assessed in order to clearly emphasize the differences existing between realistic building blocks of living matter and the controverted water-medium surrogate. Consequences in radiobiology and radiotherapy will be discussed in particular in view of treatment planning refinement aiming at better radiotherapy strategies.