Cross sections for low-energy (1-100 eV) electron elastic and inelastic scattering in amorphous ice

A model calculation of coherence effects in the elastic backscattering of very low energy electrons (1-20 eV) from amorphous ice

PURPOSE

Backscattering of very low energy electrons in thin layers of amorphous ice is known to provide experimental data for the elastic and inelastic cross sections and indicates values to be expected in liquid water. The extraction of cross sections was based on a transport analysis consistent with Monte Carlo simulation of electron trajectories. However, at electron energies below 20 eV, quantum coherence effects may be important and trajectory-based methods may be in significant error. This possibility is here investigated by calculating quantum multiple elastic scattering of electrons in a simple model of a very small, thin foil of amorphous ice.

METHOD

The average quantum multiple elastic scattering of electrons is calculated for a large number of simulated foils, using a point-scatterer model for the water molecule and taking inelastic absorption into account. The calculation is compared with a corresponding trajectory simulation.

RESULTS

The difference between average quantum scattering and trajectory simulation at energies below about 20 eV is large, in particular in the forward scattering direction, and is found to be almost entirely due to coherence effects associated with the short-range order in the amorphous ice. For electrons backscattered at the experimental detection angle (45° relative to the surface normal) the difference is however small except at electron energies below about 10 eV.

CONCLUSION

Although coherence effects are in general found to be strong, the mean free path values derived by trajectory-based analysis may actually be in fair agreement with the result of an analysis based on quantum scattering, at least for electron energies larger than about 10 eV.

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.

Damage induced to DNA by low-energy (0-30 eV) electrons under vacuum and atmospheric conditions

In this study, we show that it is possible to obtain data on DNA damage induced by low-energy (0-30 eV) electrons under atmospheric conditions. Five monolayer films of plasmid DNA (3197 base pairs) deposited on glass and gold substrates are irradiated with 1.5 keV X-rays in ultrahigh vacuum and under atmospheric conditions. The total damage is analyzed by agarose gel electrophoresis. The damage produced on the glass substrate is attributed to energy absorption from X-rays, whereas that produced on the gold substrate arises from energy absorption from both the X-ray beam and secondary electrons emitted from the gold surface. By analysis of the energy of these secondary electrons, 96% are found to have energies below 30 eV with a distribution peaking at 1.4 eV. The differences in damage yields recorded with the gold and glass substrates is therefore essentially attributed to the interaction of low-energy electrons with DNA under vacuum and hydrated conditions. From these results, the G values for low-energy electrons are determined to be four and six strand breaks per 100 eV, respectively.

Comparisons of calculations with PARTRAC and NOREC: transport of electrons in liquid water

Monte Carlo computer models that simulate the detailed, event-by-event transport of electrons in liquid water are valuable for the interpretation and understanding of findings in radiation chemistry and radiation biology. Because of the paucity of experimental data, such efforts must rely on theoretical principles and considerable judgment in their development. Experimental verification of numerical input is possible to only a limited extent. Indirect support for model validity can be gained from a comparison of details between two independently developed computer codes as well as the observable results calculated with them. In this study, we compare the transport properties of electrons in liquid water using two such models, PARTRAC and NOREC. Both use interaction cross sections based on plane-wave Born approximations and a numerical parameterization of the complex dielectric response function for the liquid. The models are described and compared, and their similarities and differences are highlighted. Recent developments in the field are discussed and taken into account. The calculated stopping powers, W values, and slab penetration characteristics are in good agreement with one another and with other independent sources.

Electron Emission from Foils and Biological Materials after Proton Impact

Electron emission spectra from thin metal foils with thin layers of water frozen on them (amorphous solid water) after fast proton impact have been measured and have been simulated in liquid water using the event-by-event track structure code PARTRAC. The electron transport model of PARTRAC has been extended to simulate electron transport down to 1 eV by including low-energy phonon, vibrational and electronic excitations as measured by Michaud et al. (Radiat. Res. 159, 3-22, 2003) for amorphous ice. Simulated liquid water yields follow in general the amorphous solid water measurements at higher energies, but overestimate them significantly at energies below 50 eV.

Analysis of band broadening in vibrational high-resolution electron-energy-loss spectra of condensed methane

High-resolution vibrational electron-energy-loss spectra of multilayer condensed films of methane recorded at 20 K show a strong tailing of the vibrational bands that clearly exceeds the instrumental resolution. At low incident electron energy, this tailing is remarkably less important for the dipole-allowed bending vibration (nu(4)) than for other bands. Also, the tailing becomes less pronounced with increasing size of the molecule as demonstrated by spectra of ethane and heptane recorded under the same conditions. Dipole coupling, rotational broadening, and multiple inelastic scattering have been considered as origins of this band broadening. While the first two effects can be excluded, multiple scattering involving a low-frequency phonon band provides a reasonable explanation as demonstrated by simulations of the spectrum of methane using a classical two-stream model. A lower phonon frequency in the cases of the larger molecules is held responsible for the better resolved vibrational signals in the spectra of ethane and heptane.

Accurate electron inelastic cross sections and stopping powers for liquid water over the 0.1-10 keV range based on an improved dielectric description of the Bethe surface

Electron inelastic cross sections and stopping powers for liquid water over the 0.1-10 keV range are presented based on a recently developed dielectric response model for liquid water (D. Emfietzoglou, F. Cucinotta and H. Nikjoo, Radiat. Res. 164, 202-211, 2005) that is consistent with the experimental data over the whole energy-momentum plane. Both exchange and second-order Born corrections are included in a material-specific way using the dielectric functions of liquid water. The numerical results are fitted by simple analytic functions to facilitate their further use. Compared to previous studies, differential cross sections are shifted toward smaller energy losses resulting in smaller inelastic and stopping cross sections with differences reaching, on average, the approximately 20% and approximately 50% level, respectively. Contrary to higher energies, it is shown that the dispersion model for the momentum dependence of the dielectric functions (Bethe ridge) is as important as the optical model used. Within the accuracy of the experimental data (a few percent) upon which our dielectric model is based, the calculations are "exact" to first order, while the uncertainty of the results beyond first order is estimated at the 5-10% level. The present work overcomes the limitations of Bethe's theory at low energies by a self-consistent account of inner-shell effects and may serve to extend the ICRU electron stopping power database for liquid water down to 100 eV with a level of uncertainty similar to that for the higher-energy values.

Excitation-energy dependence of the mechanism for two-photon ionization of liquid H2O and D2O from 8.3to12.4eV

Transient absorption measurements monitor the geminate recombination kinetics of solvated electrons following two-photon ionization of liquid water at several excitation energies in the range from 8.3 to 12.4 eV. Modeling the kinetics of the electron reveals its average ejection length from the hydronium ion and hydroxyl radical counterparts and thus provides insight into the ionization mechanism. The electron ejection length increases monotonically from roughly 0.9 nm at 8.3 eV to nearly 4 nm at 12.4 eV, with the increase taking place most rapidly above 9.5 eV. We connect our results with recent advances in the understanding of the electronic structure of liquid water and discuss the nature of the ionization mechanism as a function of excitation energy. The isotope dependence of the electron ejection length provides additional information about the ionization mechanism. The electron ejection length has a similar energy dependence for two-photon ionization of liquid D(2)O, but is consistently shorter than in H(2)O by about 0.3 nm across the wide range of excitation energies studied.

Monte Carlo simulation of water radiolysis for low-energy charged particles

The paper describes the development of chemical modules simulating the prechemical and chemical stages of charged particle tracks in pure liquid water. These calculations are based on our physical track structure codes for electrons and ions (KURBUC, LEPHIST and LEAHIST) which provide the initial spatial distribution of H2O+, H2O* and subexcitation electrons at approximately 10(-15) s. We considered 11 species and 26 chemical reactions. A step-by-step Monte Carlo approach was adopted for the chemical stage between 10(-12) s and 10(-6) s. The chemistry codes enabled to simulate the non-homogeneous chemistry that pertains to electron, proton and alpha-particle tracks of various linear energy transfers (LET). Time-dependent yields of chemical species produced by electrons and ions of different energies were calculated. The calculated primary yields (G values at 10(-6) s) of 2.80 for OH and 2.59 for e(aq)- for 1 MeV electrons are in good agreement with the published values. The calculated G values at 10(-6) s for a wide range LETs from of 0.2 to 235 keV microm(-1) were obtained. The calculations show the LET dependence for OH and H2O2. The electron penetration ranges were calculated in order to discuss the role of low energy electrons.

A complete dielectric response model for liquid water: a solution of the Bethe ridge problem

We present a complete yet computationally simple model for the dielectric response function of liquid water over the energy-momentum plane, which, in contrast to earlier models, is consistent with the recent inelastic X-ray scattering spectroscopy data at both zero and finite momentum transfer values. The model follows Ritchie's extended-Drude algorithm and is particularly effective at the region of the Bethe ridge, substantially improving previous models. The present development allows for a more accurate simulation of the inelastic scattering and energy deposition process of low-energy electrons in liquid water and other biomaterials. As an example, we calculate the stopping power of liquid water for electrons over the 0.1-10 keV range where direct experimental measurements are still impractical and the Bethe stopping formula is inaccurate. The new stopping power values are up to 30-40% lower than previous calculations. Within the range of validity of the first Born approximation, the new values are accurate to within the experimental uncertainties (a few percent). At the low end, the introduction of Born corrections raises the uncertainty to perhaps approximately 10%. Thus the present model helps extend the ICRU electron stopping power database for liquid water down to about two orders of magnitude with a comparable level of uncertainty.

High-LET Radiolysis of Liquid Water with 1H+, 4He2+, 12C6+, and 20Ne9+ Ions: Effects of Multiple Ionization

Monte Carlo simulations are used to investigate the effects of multiple ionization of water molecules on the yields of formation of free radical and molecular species, including molecular oxygen, in the radiolysis of pure, deaerated liquid water by using different types of radiation (1H+, 4He2+, 12C6+, and 20Ne9+ ions) up to approximately 900 keV/microm, at neutral pH and 25 degrees C. Taking into account the double, triple, and quadruple ionizations of water, the primary (or "escape") yields (at 10(-6) s) of the various radiolytic species (G(e(aq)-), G(H*), G(H2), G(*OH), G(HO2*/O2*-), and G(H2O2) are calculated as a function of the linear energy transfer (LET) of the radiation. Our results quantitatively reproduce the large increase observed in G(HO2*/O2*-) at high LET. Under the conditions of this study, the mechanisms of triple and quadruple ionizations contribute only weakly to the production of HO2*/O2*-. With the exception of protons, our calculations also simultaneously predict a maximum in G(H2O2) corresponding to the LET of approximately 4.5-MeV helium ions (approximately 100 keV/microm) and approximately 110-MeV carbon ions (approximately 180 keV/microm). This maximum occurs where G(HO2*/O2*-) begins to rise sharply, suggesting, in agreement with previous experimental data, that the yields of HO2*/O2*- and H2O2 are closely linked. Moreover, our results show a steep increase in the initial and primary yields of molecular oxygen with increasing LET, giving support to the "oxygen in heavy-ion tracks" hypothesis. By contrast, it is found that, in the whole LET range considered, the incorporation of multiple ionization in the simulations has only little effect on the variation of our computed G(e(aq)-), G(H*), G(H2), and G(*OH) values as a function of LET. As expected, G(e(aq)-) and G(*OH) decrease continuously with increasing LET. G(H*) at first increases and then decreases at high LET. Finally, G(H2) monotonically rises with increasing LET. Our calculated yield values compare generally very well with experiment.

The effect of model approximations on single-collision distributions of low-energy electrons in liquid water

The development of cross sections for the inelastic interaction of low-energy electrons with condensed tissue-like media is best accomplished within the framework of the dielectric theory. In this work we investigate the degree to which various model approximations, used in the above methodology, influence electron single-collision distributions. These distributions are of major importance to Monte Carlo track structure codes, namely, the energy-loss spectrum, the inelastic inverse mean free path, and the ionization efficiency. In particular, we make quantitative assessment of the influence of (1) the optical data set, (2) the dispersion algorithm, and (3) the perturbation and exchange Born corrections. It is shown that, although the shape and position of the energy-loss spectrum remains almost fixed, its peak height may vary by up to a factor of 1.5. Discrepancies in the calculated inelastic inverse mean free path are largely within 20-30% above 100 eV; they increase drastically, though, at lower energies. Exchange and perturbation Born corrections increase gradually below 1 keV leading to a approximately 30 to 40% reduction of the inverse mean free path at 100 eV. The perturbation effect contributes more than the exchange effect to this reduction. Similar to the dispersion situation, the effect of Born corrections at lower energies is also unclear since the models examined disagree strongly below 100 eV. In comparison, the vapor data are higher than the liquid calculations by 20 to 50% as the energy decreases from 1 to 0.1 keV, respectively. The excitation contribution is the main cause of this difference, since the ionization efficiency in the liquid levels off at approximately 90%, whereas the plateau value for the vapor is approximately 70%. It is concluded that electron inelastic distributions for liquid water, although in some respects distinctively different from the vapor phase, have associated uncertainties that are comparable in magnitude to the phase differences. The situation below 100 eV is uncertain.

Absolute vibrational and electronic cross sections for low-energy electron (2–12 eV) scattering from condensed pyrimidine

Low-energy vibrational and electronic electron-energy-loss (EEL) spectra of pyrimidine condensed on a thin film of solid argon held at 18 K are reported for the incident-energy range of 2-12 eV. Sensitivity to symmetry and spin forbidden transitions as well as correlations to the triplet states of benzene make it possible to ascribe the main features, below 7 eV in the electronic part of the EEL spectrum, to triplet transitions. The lowest EEL feature with an energy onset at 3.5 eV is attributed to a transition to the (3)B(1)(n-->pi(*)) valence electronic state and the next triplet n-->pi(*) transition to a (3)A(2) state located around 4.5 eV. The remaining EEL features at 4.3, 5.2, 5.8, and 6.5 eV are all assigned to pi-->pi(*) transitions to states of symmetry (3)B(2), (3)A(1), (3)B(2), and (3)B(2)+(3)A(1), respectively. The most intense maximum at 7.6 eV is found to correspond to both (1)B(2) and (1)A(1) transitions, as in the vacuum ultraviolet spectra. Absolute inelastic cross sections per scatterer are derived from a single collision treatment described herein. Their values are found to lie within the 10(-17) cm(2) range for both the electronic and the vibrational excitations. Features in the energy dependence of the cross sections are discussed, whenever possible, by comparison with data and mechanisms found in the gas phase. A maximum over the 4-5 eV range is attributed to a B (2)B(1) shape resonance and another one observed in the 6-7 eV range is ascribed to either or both sigma(*) shape resonances of (2)A(1) and (2)B(2) symmetries.

Intrinsic and extrinsic factors in anion electron-stimulated desorption: D− from deuterated hydrocarbons condensed on Kr and water ice films

The results of D(-) ion desorption induced by 3-20 eV electrons incident on condensed CD(4), C(2)D(6), C(3)D(8), C(2)D(4), and C(2)D(2) are presented. These compounds were deposited in submonolayer amounts on the surfaces of multilayer solid films of Kr and nonporous and porous amorphous ice. While desorption of the D(-) anions proceeds via well-known processes, i.e., dissociative electron attachment (DEA) and dipolar dissociation, significant perturbations of these processes due to presence of the different film substrates are observed. We have shown that it is possible to distinguish between the character and nature of these perturbations. The presence of the nonporous ice perturbs the D(-) desorption intensity by affecting the intrinsic properties of the intermediate anion states through which dissociation proceeds. On the other hand, the presence of the porous ice introduces extrinsic effects, which can affect electron energy losses prior to their interaction with the hydrocarbon molecule and/or the energies and intensities of the fragment species after dissociation. Simple mechanisms responsible for the observed variations in the intensities of desorbed anionic signals are proposed and discussed. Electron transfer from transient anion states to electron states of the substrate film or nearby hydrocarbon molecules appear as the most efficient mechanism to reduce the magnitude of the DEA process.

Analysis of Low-Energy Electron Track Structure in Liquid Water

An implementation is presented of interaction cross sections for non-relativistic electron track structure simulations. The model, incorporating liquid-phase cross sections for inelastic interactions and improved algorithms for elastic scattering, is applied to Monte Carlo simulation of the track structure of low-energy electrons. Benchmark distributions and mean values are presented for several measures of penetration distances that characterize the general physical extent of the track structure. The results indicate that, except for the last approximately 500 eV of energy loss, electron tracks have a quasi-linear character; this suggests that a major part of an electron track may be reasonably described by a lineal-energy-like characterization.

Vibrational and electronic excitations of H[sub 2]O on thymine films induced by low-energy electrons

We investigated vibrational and electronic excitations of 0.1-layer up to 2.4-layer film of H(2)O deposited on a 1.4-layer film of thymine condensed on Ar at a temperature of 18 K using high-resolution electron-energy loss (EEL) spectroscopy at the incident energy of 12 eV. The spectral contribution originating essentially from the H(2)O overlayer is obtained by separating the measured contribution from the underlying film of thymine, considering the electron beam attenuation in the H(2)O overlayer. The vibrational EEL spectrum of submonolayer amount of H(2)O on thymine, which excepts for small energy shift of the vibrational bands, is found to compare in intensity to that of the same amount of H(2)O deposited directly on the argon. The electronic energy-loss intensity near 8.6 eV, which is attributed to the excitation of (3,1)B(1) states of H(2)O in condensed phase, is observed to decrease by a factor of about 3 by the presence of the underlying film of thymine. This indicates that the corresponding cross section for excitation the (3,1)B(1) states of H(2)O by the electron impact is reduced significantly by the close proximity of the thymine molecules.

Kinetics of electron-induced decomposition of CF[sub 2]Cl[sub 2] coadsorbed with water (ice): A comparison with CCl[sub 4]

The kinetics of decomposition and subsequent chemistry of adsorbed CF(2)Cl(2), activated by low-energy electron irradiation, have been examined and compared with CCl(4). These molecules have been adsorbed alone and coadsorbed with water ice films of different thicknesses on metal surfaces (Ru; Au) at low temperatures (25 K; 100 K). The studies have been performed with temperature programmed desorption (TPD), reflection absorption infrared spectroscopy (RAIRS), and x-ray photoelectron spectroscopy (XPS). TPD data reveal the efficient decomposition of both halocarbon molecules under electron bombardment, which proceeds via dissociative electron attachment (DEA) of low-energy secondary electrons. The rates of CF(2)Cl(2) and CCl(4) dissociation increase in an H(2)O (D(2)O) environment (2-3x), but the increase is smaller than that reported in recent literature. The highest initial cross sections for halocarbon decomposition coadsorbed with H(2)O, using 180 eV incident electrons, are measured (using TPD) to be 1.0+/-0.2 x 10(-15) cm(2) for CF(2)Cl(2) and 2.5+/-0.2 x 10(-15) cm(2) for CCl(4). RAIRS and XPS studies confirm the decomposition of halocarbon molecules codeposited with water molecules, and provide insights into the irradiation products. Electron-induced generation of Cl(-) and F(-) anions in the halocarbon/water films and production of H(3)O(+), CO(2), and intermediate compounds COF(2) (for CF(2)Cl(2)) and COCl(2), C(2)Cl(4) (for CCl(4)) under electron irradiation have been detected using XPS, TPD, and RAIRS. The products and the decomposition kinetics are similar to those observed in our recent experiments involving x-ray photons as the source of ionizing irradiation.