Abstract
Understanding and quantifying the electronic inelastic interactions of swift ions and electrons in metals is fundamental for many applications of charged particle beams. A common theoretical approach is moreover desirable for the case of both types of projectiles, as large numbers of secondary electrons arise as the result of ion interaction with metals. The electronic cross section, stopping power and energy-loss straggling resulting from the interaction of swift protons, alpha particles and electrons when moving through the metals aluminum, iron, copper, molybdenum, platinum and gold, are calculated theoretically for a wide energy range of the projectiles. The model is based on the dielectric formalism, which realistically accounts for the excitation spectrum of each metal through the Mermin Energy-Loss Function–Generalized Oscillator Strength (MELF-GOS) methodology. The impact of the complexity of the excitation spectrum of each metal (encompassing interband transitions and collective excitations), as well as the different sources of (sometimes conflicting) optical data is analysed in detail. Specific interactions are considered for each projectile, such as electron capture/loss and electron cloud polarisation for ions, and indistinguishability, exchange and low-energy corrections for electrons. An estimate of possible contributions of surface excitations to the interaction probabilities of low energy electrons is given. Comparison of our results with a large collection of available experimental data shows good agreement. As a practical and useful outcome of the work, we provide analytical expressions fitting all our calculated quantities, which can be applied for simulation or comparison purposes.
1 Introduction
The interaction of energetic charged particles (either electrons or ions) with condensed matter underpins numerous techniques for both analysing and tailoring materials properties (Nastasi et al., 1996; Sigmund, 2006). Among the former, electron microscopy is nowadays routinely used not only for microscopic imaging, but also for various spectroscopies, such as electron energy-loss (EELS) or energy-dispersive X-ray (EDX) spectroscopies (; Goldstein et al., 2018; ), giving important information on the electronic structure and composition of materials. Lately, focused ion beam (FIB) imaging is also gaining popularity (Utke et al., 2008; ), and Rutherford backscattering (RBS) is a well established analytical technique (Jeynes et al., 2012), among others.
Materials modification can be accomplished by focused beams of charged particles in techniques such as electron- or ion-beam lithography (Udalagama et al., 2009; Manfrinato et al., 2014; Jesse et al., 2016; Huth et al., 2018) or ion-induced sputtering (Mookerjee et al., 2008; Sigmund, 2013). For these purposes, ion beams are particularly useful, as they deposit larger amounts of energy per unit path length as compared to electrons, but they also do it in nanometric radial distances from their path (a quality referred to as a high linear energy transfer, LET). This feature is exploited both in nanofabrication techniques (Utke et al., 2008; Sigmund, 2013; Jesse et al., 2016) and in cancer treatment by hadrontherapy (Schardt et al., 2010; Loeffler and Durante, 2013; Solov’yov, 2017). The latter benefit from the complex damage patterns induced in the scales of sensitive molecules, such as lithography resists, writing precursor molecules or the DNA in living cells. However, it should be kept in mind that the high LET of ion beams originates from the generation and propagation of large numbers of (mainly low energy) secondary electrons, being the latter the responsible of the energy propagation around the ion’s path, so ion and electron interactions are intimately interlinked in all ion beam applications.
Particularly, the interaction of charged particles with metals is of great interest, both basic and applied. Metals are not only essential components for advanced nanodevices for electronics (Janes et al., 2000), catalysis (Heiz and Landman, 2014) or many other fields, but they are also finding new applications in nanomedicine. Transition metal (such as Au or Pt) nanoparticles are being explored as sensitisers for hadrontherapy (Kuncic and Lacombe, 2018; Kempson, 2021), as they can be coated with specific molecules that allow targeting cancer cells. Ion interaction with the metals prompts the production of a considerable number of electrons, which enhance the large nanoscopic doses along ion tracks, improving their relative biological effectiveness.
From the fundamental point of view, the interaction of ions and electrons with metals is also a complex problem. Charge exchange between the impinging ion and the metal results in a collection of possible charge states, whose distribution depends on both ion’s velocity and ion-target combination, strongly affecting the ion electronic energy loss. In the case of primary electrons, the interaction is even more intricate due to their indistinguishability with the target’s electrons (Mott, 1930; ). This brings about several complexities related to exchange, but also to the complex electronic structure of metals, as individual and collective excitations need to be treated differently (). Last, but not least, in most cases low energy electrons (≲ 100 eV) are the most relevant in practice. However, larger discrepancies are reported in the literature with regard to their interaction probabilities, both from the experimental (Knapp et al., 1979; Ogawa et al., 1997; ; ; ) and theoretical (Emfietzoglou and Nikjoo, 2005; ; ; Nguyen-Truong, 2017) sides.
In this context, the dielectric formalism (Lindhard, 1954; Ritchie, 1959) (a theory dating back to Fermi (Fermi, 1940)) is gaining popularity in recent years for obtaining the necessary cross sections for different materials (Nguyen-Truong, 2017; ; Flores-Mancera et al., 2020; Poignant et al., 2020; Gibaru et al., 2021), as it represents a very convenient theoretical approach for the study of this problem. First of all, its independent description of projectile’s and target’s properties facilitates the implementation of methods to both treat primary ions and electrons. The target’s electronic excitation spectrum is modelled through its complex dielectric function, which can be both obtained from optical or electron energy-loss experiments (Palik and Ghosh, 1999; Werner et al., 2009; Xu et al., 2017; Pauly et al., 2020) and phenomenological models (Garcia-Molina et al., 2012), or from advanced ab initio calculations (Nguyen-Truong, 2017; Pedrielli et al., 2021; Taioli et al., 2021), providing a good balance between easiness of implementation and theoretical rigour. Finally, several extensions of the dielectric formalism allow the consideration of the particularities arising from low energy electrons, including their increased probabilities for surface excitation (; Zhang et al., 2004; ; ) and the improvement of the method accuracy at low energies (Emfietzoglou and Nikjoo, 2005; ). Thus, the dielectric formalism represents a method of choice for the practical yet physically motivated modelling of ion and electron interaction with metals in a wide range of energies.
In the present article, we focus our study on the electronic interactions of protons, alpha particles and electrons beams with bulk Al, Fe, Cu, Mo, Pt, and Au targets by means of the dielectric formalism. The different models to deal with ion and electron projectiles will be discussed, including the improvements necessary to describe low energy electrons. The metals chosen include Al, a representative free electron-like metal with a spectrum dominated by an intense plasmon excitation, and a selection of transition metals along the d-block in the periodic table presenting complex excitation spectra, combining plasmon, intra- and interband excitations. These materials count on with a rather extensive compilation of experimental data on electronic interactions (both for electron and ion beams) with which the calculations can be compared. Particularly, for Cu and Mo there are recent X-ray absorption fine structure spectroscopy (XAFS) measurements of very low energy electrons inelastic mean free paths (; ), which are in conflict with previous experimental data (Knapp et al., 1979; Ogawa et al., 1997; ), allowing for the discussion of some theoretical considerations. The noble metals Au and Pt are commonly used as radioenhancing nanoparticles (Kuncic and Lacombe, 2018; Kempson, 2021), and accurate information on ion and electron interaction with them is necessary to advance in the understanding of the mechanisms underlying their application in therapy. For most of these metals there are several datasets for their optical properties available (sometimes in conflict among them) (Palik and Ghosh, 1999; Werner et al., 2009; Xu et al., 2017; Pauly et al., 2020), whose consistency can be analysed by means of theoretical tests and comparison to experimental observables obtained from ion and electron interactions.
2 Theory of the energy loss of swift charged particles in metals
In the quantum theory of electronic excitations of atoms or molecules by a charged particle impact (see, e.g., (Fano, 1963; Moiseiwitsch and Smith, 1968; Inokuti, 1971; )), the most fundamental quantity is the inelastic doubly differential cross section (DDCS), , giving the probability for the incident particle of kinetic energy T to lose some specific amount of energy E = ℏω and of momentum ( and being, respectively, the momenta of the incident and the scattered particle) as a result of the interaction. The macroscopic cross section (or electronic inverse mean free path) Λ = λ−1 is related to the microscopic electronic interaction cross section σ by the relation , where is the target’s atomic density. The DDCS depends on the direct f and the exchange g scattering amplitudes as (Prasad, 1965; Rudge, 1965):The last term in Eq. (1) accounts for the interference between the direct and exchange scattering amplitudes.
In the first Born approximation (FBA) only the direct scattering amplitude is taken into account, so the exchange and interference terms in Eq. (1) are neglected. Within the FBA, the dielectric formalism provides a compact expression for the DDCS for a charged particle of charge Z, mass M and kinetic energy T inelastically interacting with the electrons of a condensed-phase medium (see Ref (Lindhard, 1954) for the derivation from electrodynamics, Refs (Fano, 1963; Moiseiwitsch and Smith, 1968; Inokuti, 1971; ) for quantum-mechanical treatments, and Ref (Nikjoo et al., 2012) for a comprehensive review):Here e is the elementary charge and ρq(k) represents the Fourier transform of the electronic charge density corresponding to the incident particle, for a given charge state q. For electrons, which are point charges, and the electron mass m appears instead of M. The quantity is the so-called material’s energy-loss function (ELF), representing its electronic excitation spectrum, where ϵ(k, E) = ϵ1(k, E) + i ϵ2(k, E) is the complex dielectric function of the material. In this work we limit ourselves to non-relativistic charged particles, so the above expression only includes longitudinal excitations. For relativistic particles, transverse excitations must be accounted for (Fano, 1963; ). The treatment of the target properties is discussed in Section 3, while the particularities of the ion and electron projectiles are discussed in detail in Section 4.
For projectiles of sufficient velocity (i.e., when the energy loss is small with respect to the projectile’s kinetic energy), the FBA is sufficiently accurate, so Eq. (2) will be used in what follows; all equations without subscript will correspond to FBA results, except when otherwise stated. At low and intermediate energies, different projectile charge states q are possible for the case of ion beams, which have to be taken into account, as will be explained in Section 4.1. However, when the primary projectile is an electron, its indistinguishability and exchange with the secondary (target) electron must be considered; the influence of this contribution will be more sizeable when both electrons have similar velocities, i.e., for the lower primary electron energies (of the order of tens and hundreds of eV). The treatment of the exchange, together with low-energy corrections for electrons, are discussed in Section 4.2.
Appropriate integration of the DDCS over energy and momentum transfers yields integral energy-loss quantities that are necessary for modelling the propagation of charged particles in the condensed phase metals, and which can be experimentally measured. The inverse mean free path (IMFP) is given by:the stopping power is given by:and the energy-loss straggling is:While the mean free path λ = Λ−1 represents the average distance between inelastic collisions, Λ denotes the average number of inelastic interactions per unit path length of the charged particle. The stopping power S characterises the average energy lost by the projectile per unit path length, and Ω2 corresponds to the second moment of the energy-loss distribution, related to the statistical deviations from the average energy lost per unit path length. As it will be discussed in Section 4, the energy and momentum transfer limits depend on both the type of projectile and on the characteristics of the electronic excitation spectrum. Whereas Eq. (3) is used mainly for electrons, particular versions of Eqs. (3), (4) and (5) will be given later for ions and electrons.
3 Description of the target excitation spectrum
In the dielectric framework (Lindhard, 1954; Fermi, 1940; Ritchie, 1957) the electronic excitation spectrum of a material is connected to its energy-loss function (ELF), which can be divided into the excitation of the loosely-bound outer-shell electrons and of the atomic-like inner-shell electrons (Heredia-Avalos et al., 2005):The second term in the right hand side of the equation accounts for the excitation of the inner-shells (with quantum numbers nl) of the atoms j forming the material with stoichiometric weight νj. Their ELF can be obtained by means of atomic hydrogenic generalised oscillator strengths (GOS) for the K-, L- and M-shells (Heredia-Avalos et al., 2005), which analytically provide the dispersion relation over the whole (k, E)-space:where are the GOS of the hydrogenic wavefunctions that are obtained using an effective nuclear charge for each inner shell. is a Heaviside step function with the threshold at the ionization energy of the nl-shell of the atom j.
The ELF describing the excitation of the outer-shell electrons can be obtained either from optical experiments (k = 0) extended to finite momentum transfer by appropriate models (Garcia-Molina et al., 2012), or alternatively by means of time-dependent density functional theory (TDDFT) calculations (Nguyen-Truong, 2017; Pedrielli et al., 2021; Taioli et al., 2021; Taioli et al., 2023). While the latter approach allows a totally ab initio calculation of the excitation spectrum (which could be very convenient for targets that are difficult to study experimentally, such as nanostructured materials), the description of bulk matter has typically relied on the use of optical experimental data, which is usually rather reliable (Raether, 1965; ). There exist many techniques to study the optical properties of materials, and several sources for metals are available, such as the classical compilation of refractive indexes and extinction coefficients by Palik and Ghosh (Palik and Ghosh, 1999) or the most recent electron energy-loss measurements by Werner et al. (Werner et al., 2009), among others.
In the optical limit, the experimental ELF (k = 0, E) can be described by a linear combination of Drude-type functions with energy thresholds ():where Ai, Ei and γi are fitting parameters denoting the intensities, positions and widths of the peaks comprising the optical spectrum; notice that the form was sometimes used, instead of , but both provide the same result. Eth,i are threshold energies, and F(E − Eth,i) is a switching function which can be either a simple Heaviside step Θ(E − Eth,i) or a continuous function of the form:which provides a soft onset for the excitation of some outer shells. The additional fitting parameter Δi determines the steepness of the soft threshold.
Among the different available approaches to extend the optical outer-shells ELF to finite momentum transfers (ℏk ≠ 0) (Garcia-Molina et al., 2012), here we use Mermin energy-loss functions (Mermin, 1970), which have been successfully applied to many condensed-phase targets by means of the so-called Mermin Energy-Loss Function–Generalised Oscillator Strengths (MELF-GOS) methodology (Heredia-Avalos et al., 2005; Heredia-Avalos and Garcia-Molina, 2007; ; ; Garcia-Molina et al., 2011; ; ). The Mermin functions, which are identical to the Drude functions at the optical limit, automatically account for the dispersion to finite k-values including the damping of the excitations.
Traditionally, the sum of Drude-type functions in Eq. (8) is just meant to provide an appropriate fitting of the optical ELF, assessed by the fulfillment of physically motivated sum rules (Smith and Shiles, 1978; Tanuma et al., 1993; Smith, 1998). However, following detailed parameterisations of the excitation spectra (ϵ2 (k = 0, E)) of materials such as liquid water (; ) and some metals and oxides (Kwei et al., 1993), it was suggested in Ref. () that Drude ELFs can be also approximately identified with particular energy-loss channels, such as plasmon, intra- and interband excitations.
The quality of a model ELF can be checked in terms of sum-rules (Tanuma et al., 1993), of which it is particularly important the f-sum rule, which gives the effective number of electrons participating in the excitations with energy transfers lower or equal than E:This equation, when using the full ELF, Eq. (6), must converge to the total number Zt of electrons per atom (or molecule, in the case of compounds) in the target material, when E → ∞. However, it was also shown in Ref. () that individual contributions of the outer-shells ELF given by the corresponding terms in Eq. (8) should also approximately converge to the number of electrons expected in a particular transition. Therefore, the f-sum rule will be very useful for assessing the quality of a particular set of optical data, both in terms of total and partial numbers of electrons. Additionally, the Kramers-Kronig (KK) sum rule must be also fulfilled for the complete ELF, which for conductors becomes (Tanuma et al., 1993):
4 Description of the projectile
4.1 Ion projectiles
As commented in Section 2, exchange and interference terms in the DDCS can be ignored for ion projectiles, and FBA results are discussed in the following. Ions large mass M in comparison with electron mass m makes them more energetic than electrons for the same velocity, thus the dielectric formalism becomes valid for ions at lower velocities as compared to electron projectiles (Rudd et al., 1992). Later on, the range of validity of the dielectric approach will be checked in comparison to the available experimental data for ion projectiles.
However, ions can lose electrons and capture them from the target, dynamically changing their charge state q, which affects their energy-loss quantities through the term in Eq. (2). Besides, the ions can also lose energy during the electron exchange events (), what is accounted for through the term SC&L, to be discussed below. The displacement of the electronic cloud with respect to the ion nucleus, due to the induced electric field at the position of the projectile (i.e., the self-induced electric field), leads to a modification of Eq. (4) for the stopping power, which is now denoted by Spol,q because it explicitly accounts for the polarization of the projectile charge distribution (Heredia-Avalos and Garcia-Molina, 2002) and whose form will be detailed later. For a given kinetic energy T, an equilibrium between electron capture and loss processes is quickly reached, giving place to equilibrium charge fractions ϕq(T). The stopping power is evaluated in the following way:where the sum goes over all possible charge states q, typically from 0 to Z. Analogously, the energy-loss straggling is calculated as:The charge fractions ϕq(T) can be obtained from the parameterisation included in the CasP code (Schiwietz and Grande, 2001), based on a large compilation of experimental measurements for several projectiles and targets.
The model implemented in Ref. (Heredia-Avalos and Garcia-Molina, 2002) for Spol,q is used in this work for both protons and alpha particles:The previous equation can be written as Spol,q = (Sq + Se) + Sinterf, where the terms in parentheses correspond to the stopping due to the nucleus and to the electron cloud, while the last term accounts for their interference, in clear analogy with the result obtained by Arista () for the energy loss of a pair of charges in correlated motion. The present expression corrects some misprints appearing in Refs. (; ), which did not affect the results presented in those works. The energy-loss straggling including polarisation, , is calculated according to the same equation, but replacing E by E2 in the argument of the integrals.
As for the integration limits, they are obtained from energy and momentum conservation in an inelastic collision. Assuming scattering with a free electron at rest, the upper limit in the energy transfer is E+ = 4 (m/M)T (Rudd et al., 1992). For metals, which have no band gap, the lower limit for the energy transfer is E− = 0. The limits for the momentum transfer are:which correspond to complete forward (ℏk−) or backward (ℏk+) scattering.
The displacement of the center of the electron cloud from its nucleus is , where αq is the projectile’s polarisability. The projectile self-induced electric field is given by (Heredia-Avalos and Garcia-Molina, 2002):This expression corrects misprints appearing in Ref. (Heredia-Avalos and Garcia-Molina, 2002), which did not affect its results.
The electron capture and loss processes can also lead to additional energy losses. For both protons and alpha particles, an approximate model to account for this energy loss was implemented in Ref. (), based on the charge-exchange cross section from Ref. (). The electron capture and loss (C&L) contribution, if only single-electron exchange processes are considered, is given by ():where σi→j are the cross sections for electron loss or capture, changing from the charge state i to the j, and ΔTi→j is the energy loss associated to these processes. When charge equilibrium is reached, ϕq(T) σq→q+1(T) = ϕq+1(T) σq+1→q(T), so the previous expression reduces to:
The energy lost by the projectile in an electron loss and capture cycle can be estimated as , where v is the ion’s velocity, with Bp and Bt being, respectively, the first ionization potentials of the projectile and the target atoms. These values can be obtained from the literature (Huheey et al., 1993) or estimated as the energy of an electron in a hydrogenic Coulomb field with an effective Slater’s nuclear charge (Slater, 1930).
The loss cross section for K-shell electrons (for hydrogen and helium ions, in our case) (; ) is given, following Ref. (), by:where nK is the number of K-shell electrons in the projectile, v0 is the Bohr velocity, v is the projectile velocity, ϒ is a screening length (introduced in the next paragraph), and Zt is the target’s atomic number.
The Fourier transform of the electronic density ρq(k) is obtained here through the modified Brandt-Kitagawa approach (), which is suitable for ions with 1 or 2 bound electrons:Here, the screening length ϒ is given by ():with a0 being the Bohr radius.
4.2 Electron projectiles
One of the main aspects that influences the energy loss for primary electrons (i.e., projectiles) is their indistinguishability from the secondary electrons (i.e., those resulting from the ionisation of the target atoms). This affects their maximum energy loss E+ in two ways. First, the primary electron moves in the conduction band of the metal with a kinetic energy T that is measured from the bottom of that band, which is filled with target electrons up to the Fermi energy EF. Therefore, Pauli’s exclusion principle forbids that the primary electron losses an energy larger than T − EF, as otherwise it would fall into already occupied levels.
Apart from this, any energy loss E that results in the transition of a bound electron from a band with binding energy Ba to the conduction band, with final kinetic energy W = E − Ba, results into two indistinguishable electrons. By convention, the primary electron is always the most energetic one after the collision, so this situation imposes an additional restriction to the maximum energy loss E+,a = (T + Ba)/2. Since this maximum energy loss is now shell-dependent through Ba, Eqs. (3), (4) and (5) have to be rewritten for electron projectiles, within the FBA, as:and:Notice that we have recovered in these equations the subscript FBA because, later on, we will introduce corrections to the first Born approximation. For any intra- or interband transition, the energy integration limits are:where the subscript ‘a’ refers either to an outer-shell Mermin function i or to an inner-shell with quantum numbers nl of an atom j, with Ba = Eth,i for the former and for the latter. It should be noted that the above energy integration limit E+ corresponds to individual electronic transitions. However, as noted in Refs. (; ), when the primary electron excites a plasmon or collective excitation, the quasi-particle is distinguishable from the primary electron, and in such a case this upper limit is simply:The integration limits for the momentum transfer are given by Eq. (15).
The FBA ignores the electron exchange and interference terms in Eq. (1). However, for primary electrons with an energy closer to that of the target electrons ( eV), this approximation does not hold anymore, and Eq. (1), which includes terms beyond the FBA, should be written as:where d2ΛFBA(T, E, k)/dE dk ∝|f|2 and d2Λxc(T, E, k)/dE dk ∝|g|2 − Re[fg*].
Despite the existence of more rigorous and complex methodologies accounting for the exchange amplitude (Rudge, 1968; Emfietzoglou et al., 2017), Ochkur developed convenient analytical approximations (based on the Born-Oppenheimer perturbation treatment (Oppenheimer, 1928)), which are proportional to the FBA direct amplitude (Ochkur, 1965). The Born-Ochkur (BO) exchange amplitude for the ionisation process is given by (Ochkur, 1965; Prasad, 1965; Rudge, 1968):where fFBA denotes the direct scattering amplitude calculated within the FBA, , and is the momentum of the secondary electron (note that Ochkur proposed a different exchange amplitude for excitation to discrete levels (Rudge, 1968; )). Any individual transition, which promotes a bound electron to the conduction band, is considered in this context an ionization, so the exchange term of the DDCS is rewritten, according to Eq. (28), as:where:which reproduces the classical Mott exchange for binary collision theory (Rudd et al., 1996) at high electron energies and is similar to the expression for ionisation given by Fernández-Varea et al. (Fernández-Varea et al., 1993). Note, however, that plasmon or collective excitations are distinguishable from the primary electron, so in this case Fxc (T, k) = 0.
Additionally, for primary electrons with energies T ≤ 100 eV, primary and secondary electrons may have similar kinetic energies, so the FBA expressions, Eqs. (22), (23), and (24), require further corrections to increase their accuracy. Some theoretical approaches have been developed to try to account for these corrections (Emfietzoglou et al., 2017), although usually they result too complex to be easily generalised to arbitrary materials. As a consequence, other simpler phenomenological approaches are usually preferred (Emfietzoglou and Nikjoo, 2005). In this work, we implement a classical Coulomb-field correction, in which the primary electron gains an additional amount of potential energy due to the interaction with the field of the target atoms, then increasing its kinetic energy T to T ′ as (Emfietzoglou and Nikjoo, 2005):where Ba is the electron binding energy (Ba = Eth,i for the target outer shells and for the inner nl-shell of the j-th atom). No particular correction is applied for unbound plasmon excitations.
As a result of the above low-energy corrections to the FBA plus addition of the exchange term, the inverse mean free path for electrons is finally written as:Similar expressions arise for the stopping power S(T) = Scorr(T) + Sxc(T) and the energy-loss straggling , where the integrands are further multiplied by E and E2, respectively. The integration limits for the energy transfer are those given by Eqs. (25) or (26), depending on the individual or collective nature of the excitations, while the momentum transfer limits ℏk± are given by Eq. (15). The Born-corrected momentum limits are (Vriens, 1966):with Ba = Eth,i for the outer shells and for the inner nl shells of the j-th atom.
5 Results and discussion
In the following, the methodology described above will be used to obtain the stopping power and energy-loss straggling of protons and alpha particles, as well as the inelastic cross section (ICS), stopping power and energy-loss straggling of electrons in the bulk (i.e., neglecting surface effects) metals Al, Fe, Cu, Mo, Pt and Au. Rather than discussing the results for the metals in order of ascending atomic number, we rather proceed in a way that facilitates the discussion of the different physical considerations involved. The materials whose electron inelastic mean free paths were already analysed in Ref (). (Al, Au and Cu, which results for protons and alpha particles also prsented here) are discussed first, and then continue with the other targets Mo, Fe, and Pt. Aluminum will serve to discuss the nature of collective excitations for electrons, as well as the role of electron capture and loss, and polarisation processes for ions. With gold, the importance of low-energy corrections to the dielectric formalism for electrons will become evident. Copper, molybdenum, iron and platinum will serve to discuss the possible influence of surface excitations for low energy electrons (when required for some experimental situations), as well as to point out the influence of different and conflicting experimental datasets for the optical ELF of the materials.
The main physical properties of the targets analysed are gathered in Table 1, namely, atomic number Zt, atomic weight At, mass density ρt, Fermi energy EF, and bulk plasmon energy Epl. As outer- and inner-shell electron excitations are treated differently in our model (see section 3), we specify in Table 2 which shells are considered as outer or as inner for each material. An additional fundamental property of each metal can be obtained from the integration of its optical ELF over the entire excitation energy range, namely, the mean excitation energy I (Shiles et al., 1980):which is the main parameter determining the stopping power at high ion or electron energies through the Bethe equation (Inokuti, 1971). In Table 3 the values of I obtained in this work for each target, from different sources of the optical experimental ELF (Palik and Ghosh, 1999; Werner et al., 2009; Pauly et al., 2020; Xu et al., 2017), are presented and compared to the values compiled in the ICRU Reports 37 () and 49 (ICRU, 1993). These I -values correspond to substances in the condensed phase, and the uncertainties reported by ICRU attempt to take into account the uncertainties of the underlying measurements, the errors inherent to the data analysis and the dispersion in the I-values derived from various sources.
TABLE 1
| Target | Al | Fe | Cu | Mo | Pt | Au |
|---|---|---|---|---|---|---|
| Zt | 13 | 26 | 29 | 42 | 78 | 79 |
| At | 26.98 | 55.85 | 63.55 | 95.94 | 195.09 | 196.97 |
| ρt (g/cm3) | 2.70 | 7.87 | 8.96 | 10.2 | 21.4 | 19.3 |
| EF (eV) | 11.7 () | 11.1 () | 7 () | 5.83 (Nguyen-Truong, 2017) | 8.78 () | 5.53 () |
| Epl (eV) () | 15.0 | 23.0 | 19.3 | 25.2 | 22.6 | 24.8 |
Physical properties of the targets discussed in this work: atomic number (Zt), atomic mass (At), mass density (ρt), Fermi energy (EF) and plasmon energy (Epl).
TABLE 2
| Target | Al | Fe | Cu | Mo | Pt | Au |
|---|---|---|---|---|---|---|
| Outer levels | 3s2 3p1 | 4s2 3d6 | 3d10 4s1 | 4d5 5s1 | 5d9 6s1 | 5d10 6s1 |
| 2s2 2p6 | 3s2 3p6 | 3s2 3p6 | 4s2 4p6 | 5s2 5p6 4f14 | 5s2 5p6 4f14 | |
| 3s2 3p6 3d10 | 4s2 4p6 | 4s2 4p6 4d10 | ||||
| Inner shells | 1s2 | 2s2 2p6 | 2s2 2p6 | 2s2 2p6 | 3s2 3p6 3d10 | 3s2 3p6 3d10 |
| 1s2 | 1s2 | 1s2 | 2s2 2p6 | 2s2 2p6 | ||
| 1s2 | 1s2 |
Outer and inner excitation levels of the targets described by the MELF-GOS model.
TABLE 3
| Target | Al | Fe | Cu | Mo | Pt | Au |
|---|---|---|---|---|---|---|
| Palik and Ghosh’s compilation (Palik and Ghosh, 1999) | 165.1 | 329.7 | 347.5 | 469.2 | 738.6 | 761.6 |
| Werner et al. (Werner et al., 2009) | 283.3 | 309.9 | 443.4 | 739.0 | 775.4 | |
| Xu et al. (Xu et al., 2017) | 294.0 | |||||
| Pauly et al. (Pauly et al., 2020) | 503.9 | |||||
| ICRU Report 37 () | 166 ± 2 | 286 ± 9 | 322 ± 10 | 424 ± 15 | 790 ± 30 | 790 ± 30 |
| ICRU Report 49 (ICRU, 1993) | 158.3 | 291.1 | 326.3 | 424.0 | 751.1 | 741.9 |
Mean excitation energy I (in eV) of the metals discussed in this work, obtained by the MELF-GOS model from different sources of the optical ELF. Values appearing in ICRU Reports 37 () and 49 (ICRU, 1993) are also provided for comparison purposes.
5.1 Aluminum
Aluminum is a paradigmatic free electron-like metal, whose spectrum is dominated by an intense plasmon excitation at eV due of its outermost 3p and 3s electrons. Its optical properties have been intensively studied during decades, with its optical (k = 0) ELF being well known (Palik and Ghosh, 1999), which is shown in Figures 1A, B as a function of the excitation energy E. Magenta circles are experimental data from Palik and Ghosh’s compilation (Palik and Ghosh, 1999) obtained from optical constants and reflectivity of evaporated metallic Al films at room temperature in high vacuum from Refs. (; ; Mathewson and Myers, 1971; Shiles et al., 1980; Smith et al., 1983). Panel A shows in linear scale the ELF due to the outer-shell electron excitations, whereas the panel B depicts in a log-log scale the contribution due to the semi-core 2p and 2s shells, and the 1s inner-shell excitation. At large excitation energies, Henke’s (black crosses) (Henke et al., 1993) and NIST’s calculations (gray stars) (; ), based on atomic data, are also included. The red solid line in both panels represent the MELF-GOS fit to the previous data (as explained in section 3), whose parameters are given in Supplementary Table S1 in the Supplementary Material.
FIGURE 1
The effective number of electrons Neff in Al as a function of the excitation energy E, as calculated from the MELF-GOS fit, is depicted in Figures 1C, D. Dashed gray lines correspond to 3 electrons (outer shells), 11 electrons (where the K-shell excitations start) and 13 electrons (the total number of electrons). The MELF-GOS fitting allows to approximately assign Mermin functions to particular electronic excitations. In this case, we assumed that the first Mermin function accounts for the intense plasmon excitation, the second one for the interband transitions of the 3s and 3p electrons to the conduction band (which should be the ones appearing at lower excitation energies), while the third function represents the excitation of the semi-core 2p and 2s electrons to the conduction band. The excitation of the 1s inner-shell is modelled through the K-shell GOS. The quality of such an assignment can be assessed by means of their partial f-sum rules (obtained by means of Eq. (10) applied to each Mermin or GOS function), which are given in Table 4; these values are in excellent agreement with the results reported by Shiles et al. (Shiles et al., 1980). As can be seen, the first two Mermin functions nearly contain the 3 outermost 3p and 3s electrons, which mainly interact through the plasmon excitation. Similarly, the third Mermin function contains around the 8 expected 2p and 2s electrons.
TABLE 4
| Target Al (Palik and Ghosh, 1999) | ||||||
|---|---|---|---|---|---|---|
| Level | i-th MELF | Nexpect | Neff | |||
| plasmon | 1 | 11 | 2.76 | 3.29 | 11.41 | |
| 3p+3s | 2 | 3 | 0.53 | |||
| 2p+2s | 3 | 8 | 8.12 | |||
| 1s | K-shell GOS | 2 | 1.59 | |||
| Total | 13 | 13.01 | ɛrel = 0.07% | f-sum rule | ||
| 1.06 | ɛrel = 6.47% | KK-sum rule | ||||
Comparison of the expected number of electrons Nexpect in each excitation level of aluminum, and the corresponding effective number of electrons Neff, as obtained from Eq. (10) with the MELF-GOS model fit to the experimental optical ELF (Palik and Ghosh, 1999). The values of the f- and KK-sum rules (with their corresponding relative errors, compared to the expected values) appear in the last two rows.
In total, 11.41 outer-shell electrons are found. The slight excess of 0.41 electrons in the outer shells compensates for the electrons missing from the K-shell GOS. These differences arise from Pauli’s exclusion principle: as inner-shell electrons cannot transit to upper energy levels occupied by the outer-shell electrons, they lack some amount of oscillator strength; similarly, the outer electrons cannot decay to lower energy levels occupied by the inner-shell electrons, thus lacking some amount of negative oscillator strength (
The ELF of Al obtained by means of the MELF-GOS methodology over the whole (k, E) plane is used to calculate the stopping power of protons and alpha particles as explained in Section 4.1. In Figures 2A,B our results are shown by red solid lines, which are compared with available experimental data, depicted by letters, from the extensive database gathered by the late Helmut Paul (International Atomic Energy Agency Nuclear Data Services, 2021), now maintained by the International Atomic Energy Agency. The latest measurements by Moro et al. (Moro et al., 2016) are also included in the figure, although they are completely indistinguishable from the overlapped data from the database. Our theoretical results agree with the experiments at energies larger than around 50 keV/u including the region of the maximum stopping power. The underestimation at lower energies may be due to the possible inaccuracy of the dielectric formalism at this energy range, or perhaps to the existence of negatively charged ions due to electron capture and loss processes in aluminum (Penalba et al., 1992), a case that is not considered in our model. The insets of Figures 2A,B compare our calculations with the most recent experimental data around the maximum of the stopping power curve, for protons (Eppacher and Semrad, 1992; Martínez-Tamayo et al., 1996; Shiomi-Tsuda et al., 1999; Møller et al., 2002) and for helium ions (Bianconi et al., 2005; Hsu et al., 2005;
FIGURE 2

Energy-loss quantities for charged particles in Al as a function of energy: (A) electronic stopping power for H and (B) for He, (C) electronic energy-loss straggling for H and (D) for He, (E) inelastic cross section (ICS) for electrons and (F) stopping power for electrons. Red solid lines correspond to calculations from the dielectric formalism with the MELF-GOS model fit to the experimental optical ELF (Palik and Ghosh, 1999). A comparison with experimental data (symbols) is made for ions stopping power (International Atomic Energy Agency Nuclear Data Services, 2021) (the insets show the corresponding stopping power around its maximum value for the most recent experimental data (Eppacher and Semrad, 1992; Martínez-Tamayo et al., 1996; Shiomi-Tsuda et al., 1999; Møller et al., 2002; Bianconi et al., 2005; Hsu et al., 2005;
We further report results from the present model in Figures 2A,B, but without considering some of the terms appearing in Eq. (12), i.e., ignoring both the energy losses due to the projectile polarisation and the electron capture and loss processes (dotted blue line), only ignoring polarisation (dash-dotted red line) and only ignoring capture and loss (dashed dark green line). As can be seen, disregarding these processes makes the calculated stopping power noticeably deviate from the experimental data around and below the maximum (T ≤ 200 keV/u). For protons, both processes present similar contributions around the maximum, with capture and loss having an increasing influence at lower energies. For helium ions polarisation has an almost negligible contribution.
Results from other widely used codes are also included in Figures 2A,B. The cyan dashed line corresponds to the stopping values recommended by the International Commission on Radiation Units and Measurements (ICRU) (ICRU, 1993) based on the PSTAR and ASTAR programs, where the stopping power at the lower energies is based on a fit to experimental data, extrapolated to high energies using a corrected Bethe equation. The dotted green line depicts the results of the semiempirical code SRIM (Ziegler, 2013). Results of the CasP code within the Unitary Convolution Approximation (UCA) by Grande and Schiwietz (Grande and Schiwietz, 2002; Schiwietz and Grande, 2012) are depicted by a magenta dot-dashed line. The orange dot-dot-dashed line corresponds to dielectric calculations for protons by Montanari and Miraglia (Montanari and Miraglia, 2017) where, for low and intermediate energies, a nonlinear model based on a central screened potential for a projectile moving in a free-electron gas was used. We note the similarity between Montanari and Miraglia’s and our calculations at energies around the maximum stopping power and beyond, and our departure from them and from experiments for T ≤ 40 keV/u. This may point to the deficiencies of the dielectric formalism at this low energy range, due to its lack of nonlinearity. First-principles calculations based on real-time propagation within time-dependent density functional theory are shown by a short-dashed olive line (Schleife et al., 2015) for He in crystalline aluminum at energies lower than 40 keV/u. Overall, our methodology based on the dielectric formalism provides reliable stopping powers both for protons and alpha particles for energies larger than 40 keV/u, covering the maximum of the stopping curve and all larger energies within a single theoretical formalism, without any particular fitting to experimental data (apart from that to the ELF in the optical limit).
In Figures 2C,D we present results for the energy-loss straggling of H and He beams in Al as a function of the ions kinetic energy. High energy ions present an asymptotic value, known as Bohr straggling, , so present results are normalised to . Red lines depict our calculated results based on the experimental optical ELF of Al (Palik and Ghosh, 1999), which are compared to a large collection of experimental determinations (Kido and Hioki, 1983; Kido, 1986; Kido, 1987;
The results for the inelastic cross section (ICS) and stopping power for electrons in aluminum (based on the same model ELF, and obtained according to the procedures explained in section 4.2) are shown by red solid lines in Figures 2E,F. The theoretical ICS is compared with available experimental data (symbols) (Kanter, 1970b; Kanter, 1970a; Tracy, 1974;
The panel F of Figure 2 reports by a thick solid line the calculated stopping power of electrons in aluminum, along with the contribution from each type of excitation (thin lines). Again, at energies around and lower than the maximum stopping power the main contribution comes from the collective plasmon excitation, while interband transitions from the 2s and 2p bands only have an influence at high energies. The comparison with available experimental data (symbols) at energies higher than 200 eV (Kalil et al., 1959; Garber et al., 1971; Fitting, 1974; Ishigure et al., 1978;
The comparison of our calculations and the experimental data for the ICS for the case of electrons in aluminum allows us to discuss the role of collective electronic excitations, as this material is basically dominated by its bulk plasmon. Particularly, for a long time there has been a controversy around the maximum energy loss which is allowed for primary electrons in metals (
Experimental measurements for the stopping power of low energy electrons would be needed to further assess the current calculations around the maximum of the energy loss. However, experimental measurements over the entire energy range are available for the ICS, which coincide very well with our calculations, supporting the assumptions of our energy-loss model for the case of electron projectiles.
Figure 3 shows by lines the calculated energy-loss straggling for electrons in all the materials studied in this work. The black solid line corresponds to the case of Al, based on the experimental optical ELF. Unfortunately, it has been impossible for us to find experimental information on the energy-loss straggling for electron beams. However, Ashley et al. (
FIGURE 3

Electronic energy-loss straggling of electrons in all metals studied in this work, as a function of the incident energy. Solid lines correspond to calculations using the experimental optical ELF from Ref. (Palik and Ghosh, 1999), dashed lines are calculations based on the ELF obtained from REELS experiments (Werner et al., 2009), and the dotted lines for Mo and Fe are additional calculations using the optical ELFs from Refs. (Pauly et al., 2020; Xu et al., 2017), respectively. The grey dash-dotted line for Al depicts a previous calculation by Ashley et al. (
Once all the calculated results have been presented for H, He and electron beams in Al, we emphasize the fact that it is possible to provide empirical fitting formulas for all the quantities, i.e., stopping power, straggling, and inelastic cross section. Such a fitting is very convenient for the straightforward implementation of these energy-loss quantities in Monte Carlo codes. We provide in the Supplementary Material a description of all the empirical formulas used, and we provide in tabular form the parameters that can be used to fit all the calculated results. This has been done for aluminum, as well as for the rest of metals analysed in the rest of the manuscript.
5.2 Gold
Now we turn our attention to gold, which is a very important material, whose complex excitation spectrum contains both collective excitations as well as intra- and interband transitions. Two important experimental sources for the optical ELF of gold exist (Palik and Ghosh, 1999; Werner et al., 2009). Figures 4A, B depict by symbols the data from Palik and Ghosh’s compilation of refractive index and extinction coefficient (Palik and Ghosh, 1999) (magenta circles), together with the most recent data by Werner et al. (Werner et al., 2009) (blue squares) obtained from reflection electron energy-loss spectroscopy (REELS). Panel A shows in linear scale the complex excitation spectrum of the outer-shell electrons, whereas panel B depicts in log-log scale the excitation of the semi-core and inner shells. Calculations based on atomic data from Henke et al. (crosses) (Henke et al., 1993) and NIST (stars) (
FIGURE 4

(A) and (B): Optical energy-loss function, ELF (k = 0, E), of Au as a function of the energy transfer E. The experimental data come from Palik and Ghosh compilation (Palik and Ghosh, 1999) (magenta circles), and from Werner et al. (Werner et al., 2009) (blue squares); solid lines correspond to the fit by the MELF-GOS model. Calculations from Henke et al. (Henke et al., 1993) (black crosses) and NIST (
The thick solid red and dashed blue lines in Figure 4 represent the MELF-GOS fittings to these two sets of experimental data for Au, where the K, L and M electron shells are treated by GOS. The outer-shells are described by a sum of Mermin-type ELFs which, as in the case of Al, have been assigned to particular excitation channels. Thin lines in Figure 4A depict the Mermin functions attributed to each excitation (for the particular case of the fitting to Palik and Ghosh’s data (Palik and Ghosh, 1999)), and a summary of the assignment is given in Table 5. First of all, the plasmon energy of 24.8 eV has been taken from Ref. (
TABLE 5
| Target Au (Palik and Ghosh, 1999) | Target Au (Werner et al., 2009) | ||||||||
|---|---|---|---|---|---|---|---|---|---|
| Level | i-th MELF | Nexpect | Neff | Neff | |||||
| 6s | 1–3 | 1 | 11 | 1.16 | 2.16 | 11.63 | 0.87 | 1.78 | 11.11 |
| plasmon | 4 | 1.00 | 0.91 | ||||||
| 5d | 5–7 | 10 | 9.47 | 9.34 | |||||
| 5p | 8 | 6 | 6.09 | 6.09 | |||||
| 5s | 9 | 2 | 1.74 | 2.06 | |||||
| 4s+4p+4d+4f | 10 | 32 | 42.25 | 42.25 | |||||
| 3d | GOS M | 10 | 6.87 | 6.87 | |||||
| 3p | M-shell GOS | 6 | 4.68 | 4.68 | |||||
| 3s | M-shell GOS | 2 | 1.29 | 1.29 | |||||
| 2p | L-shell GOS | 6 | 2.67 | 2.67 | |||||
| 2s | L-shell GOS | 2 | 1.01 | 1.01 | |||||
| 1s | K-shell GOS | 2 | 0.94 | 0.94 | |||||
| Total | 79 | 79.16 | ɛrel = 0.20% | 78.96 | ɛrel = −0.04% | f-sum rule | |||
| 1.10 | ɛrel = 10.23% | 1.03 | ɛrel = 4.18% | KK-sum rule | |||||
Comparison of the expected number of electrons Nexpect in each excitation level of gold, and the corresponding effective number of electrons Neff, as obtained from Eq. (10) with the MELF-GOS model fit to the experimental optical ELF from Palik and Ghosh (Palik and Ghosh, 1999) and Werner et al. (Werner et al., 2009). The values of the f- and KK-sum rules appear in the last two rows.
In Figures 4C, D, the effective number of electrons Neff, Eq. (10), is plotted as a function of the transferred energy E. As can be seen, the fitting to both datasets for the optical ELF leads to the correct atomic number at high excitation energies. Thin grey dashed lines correspond to 51 electrons (including the outer shells as well as the semi-core electrons 4s, 4p, 4d, and 4f, all described by MELF), 69 (including the M-shell electrons), 77 (with the L-shell electrons already accounted for) and 79 electrons (the total number of electrons in Au). The results of the evaluation of individual f-sum rules for each excitation are gathered in Table 5, together with the total sum rules, for the optical ELFs by Palik and Ghosh (Palik and Ghosh, 1999) and Werner et al. (Werner et al., 2009). The total errors in the f-sum rules are very small for both sets of data, although the total KK-sum rule is better fulfilled for the case of Werner et al.’s data (Werner et al., 2009), which may indicate a better representation of the optical ELF at low excitation energies.
Let us take a closer look to the number of electrons contained in the outermost shells. Particularly, we should consider together excitation of the 6s and 5d electrons together with the plasmon, as both 6s and 5d electrons may contribute to the bulk plasmon. As can be seen in Table 5, for both Palik and Ghosh’s (Palik and Ghosh, 1999) and Werner et al.‘s ELF (Werner et al., 2009) the plasmon is built up mainly by 2s electrons, with a small contribution coming from 5d electrons. The third Mermin function sums up, in both cases, almost the 10 electrons expected in the 5d shell. In any case, for both experimental ELFs, the three first functions contain nearly the 11 6s+5d expected electrons. However, Werner et al.‘s data (Werner et al., 2009) matches slightly better the expected number of 6s+5d electrons, which could be related to the better fulfillment of the overall KK-sum rule. For the rest of shells, both sets of optical ELFs satisfy rather well the expected number of electrons, with a slight greater agreement in the REELS data for the 5s shell. Again, the electrons missing from the inner-shell GOS are redistributed among the outer-shell Mermin functions. The mean excitation energies obtained from optical (I = 761.6 eV) and REELS data (I = 775.4 eV) are both in between the values recommended by ICRU37 (I = 790 ± 30 eV) (
Once the two sources of optical ELFs for Au have been analysed, they can be extended to arbitrary momentum transfers by means of the MELF-GOS method in order to obtain the stopping power for protons and alpha particles. The results are depicted by solid red lines (using Palik and Ghosh’s ELF (Palik and Ghosh, 1999)) and by dotted blue lines (using Werner et al.‘s ELF (Werner et al., 2009)) in Figures 5A,B. As the two optical ELFs are not radically different, they provide very similar results, which are in rather good agreement with the large collection of experimental data (letters) (International Atomic Energy Agency Nuclear Data Services, 2021) almost in the entire energy range, both for protons and helium ions. The insets present the most recent experimental data around the maximum stopping power with the corresponding error bars, for protons (Semrad et al., 1990; Eppacher and Semrad, 1992; Martínez-Tamayo et al., 1996; Heredia-Avalos et al., 2007;
FIGURE 5

Energy-loss quantities for charged particles in Au as a function of energy: (A) electronic stopping power for H and (B) for He, (C) electronic energy-loss straggling for H and (D) for He, (E) inelastic cross section (ICS) for electrons and (F) stopping power for electrons. Red solid and blue dotted lines correspond to calculations from the dielectric formalism with the MELF-GOS model fit to the experimental optical ELF from Refs. (Palik and Ghosh, 1999; Werner et al., 2009), respectively. A comparison with experimental data (symbols) is made for ions stopping power (International Atomic Energy Agency Nuclear Data Services, 2021) (the insets show the corresponding stopping power around its maximum value for the most recent experimental data (Semrad et al., 1990; Eppacher and Semrad, 1992; Martínez-Tamayo et al., 1996; Trzaska et al., 2002; Hsu et al., 2004; Heredia-Avalos et al., 2007; Heredia-Avalos et al., 2007;
Figures 5C,D shows the energy-loss straggling of H and He beams in Au. Red solid lines depict our calculated results based on the experimental optical ELF by Palik and Ghosh (Palik and Ghosh, 1999), while dotted blue lines correspond to Werner et al.’s ELF (Werner et al., 2009). Symbols depict experimental data (
Regarding electron projectiles, Figure 5E shows the results of the inelastic cross section (ICS) in Au. The symbols correspond to experimental data (Sze et al., 1963; Kanter, 1970b; Lindau et al., 1976; Gergely et al., 2004; Tanuma et al., 2005;
Surprisingly, the calculated ICS for electrons in Au agrees rather well with experimental data down to energies as low as 2 eV above the Fermi energy, while the dielectric formalism is not expected to work as well in this energy range. At this point, it is reasonable to wonder if this is due to the use of exchange and low-energy corrections in our formalism (Section 4.2). Figure 5E shows by a green line our results when such exchange and low-energy corrections are not used. As can be seen, the agreement with experimental data (especially in the 5–10 eV range) is significantly worsened, what supports the algorithms implemented in this work for these low-energy corrections to the dielectric formalism.
Figure 5F presents the results for the stopping power of electrons in Au, where experimental data are shown by symbols (
Our calculated energy-loss straggling for electrons in Au is shown in Figure 3, for the two sets of ELF, as indicated in the figure caption.
5.3 Copper
In the case of copper, we now face a more noticeable discrepancy in its optical ELFs, depending whether it comes from optical parameters measurements (Palik and Ghosh, 1999) or from REELS experiments (Werner et al., 2009). These ELFs are depicted, respectively, by magenta circles and blue squares in Figures 6A, B, in linear scale for the outer shells (A) and log-log scale for the semi-core and inner shells (B). Henke and Ghosh’s (crosses) (Henke et al., 1993) and NIST’s (stars) (
FIGURE 6

(A) and (B): Optical energy-loss function, ELF (k = 0, E), of Cu as a function of the energy transfer E. The experimental data from Palik and Ghosh’s compilation (Palik and Ghosh, 1999) (magenta circles) and from Werner et al. (Werner et al., 2009) (square blue symbols) are presented as well as the MELF-GOS fit to these experimental data. Calculations from Henke et al. (Henke et al., 1993) and the NIST (
The MELF-GOS fitting to these data (solid red line for optical data, dotted blue line for REELS data in Figure 6) have been done assuming that the bulk plasmon is located at 19.3 eV (
As in the previous cases, the effective number of electrons as a function of the excitation energy are plotted in Figures 6C, D and the effective number of electrons per shell are those indicated in Table 6. For Cu, Palik and Ghosh’s optical ELF (Palik and Ghosh, 1999) provides a closer approximation to the 11 electrons expected in the excitation of the 4s and 3s shells (including the plasmon excitation) than Werner et al’s (Werner et al., 2009). Similarly, the errors both in the total f- and KK-sum rules are significantly smaller for the optical measurements dataset (Palik and Ghosh, 1999). Also, the mean excitation energy I can be obtained from the optical ELF from Eq. (34), getting the values of I = 347.5 eV for Palik and Ghosh’s ELF (Palik and Ghosh, 1999) and I = 309.9 eV for Werner et al.‘s ELF (Werner et al., 2009), which are in rather good agreement with the values recommended by ICRU37 (
TABLE 6
| Target Cu (Palik and Ghosh, 1999) | Target Cu (Werner et al., 2009) | ||||||||
|---|---|---|---|---|---|---|---|---|---|
| Level | i-th MELF | Nexpect | Neff | Neff | |||||
| 4s | 1 | 1 | 11 | 0.49 | 1.06 | 11.26 | 0.27 | 1.18 | 11.97 |
| plasmon | 2 | 0.57 | 0.91 | ||||||
| 3d | 3–5 | 10 | 10.20 | 10.79 | |||||
| 3p | 6 | 6 | 6.43 | 5.79 | |||||
| 3s | 7 | 2 | 2.42 | 2.41 | |||||
| 2p | L-shell GOS | 6 | 6.06 | 6.06 | |||||
| 2s | L-shell GOS | 2 | 1.54 | 1.54 | |||||
| 1s | K-shell GOS | 2 | 1.29 | 1.29 | |||||
| Total | 29 | 29.00 | ɛrel < −0.01% | 29.05 | ɛrel = 0.18% | f-sum rule | |||
| 1.02 | ɛrel = 2.05% | 1.04 | ɛrel = 4.00% | KK-sum rule | |||||
Comparison of the expected number of electrons Nexpect in each excitation level of copper, and the corresponding effective number of electrons Neff, as obtained from Eq. (10) with the MELF-GOS model fit to the experimental optical ELF by Palik and Ghosh (Palik and Ghosh, 1999) or Werner et al. (Werner et al., 2009). The values of the f- and KK-sum rules appear in the last two rows.
Let us analyse the impact of both optical ELFs, extended to arbitrary momentum transfer by means of the MELF-GOS methodology, in the calculation of the stopping power of H and He ions in Cu. Figures 7A,B show by red solid and blue dotted lines the calculations done with the optical and REELS ELFs, respectively, together with the large compilation of experimental data (letters) (International Atomic Energy Agency Nuclear Data Services, 2021). For Cu, the slightly more intense ELF by Werner et al. (Werner et al., 2009) gives slightly larger stopping powers around the maxima as compared to Palik and Ghosh’s ELF (Palik and Ghosh, 1999), although both calculations are close to the cloud of experimental points in the entire energy range. The insets focus around the maxima of the stopping power, showing the most recent sets of experimental determinations with their error bars for H (Kido and Hioki, 1983; Semrad et al., 1983;
FIGURE 7

Energy-loss quantities for charged particles in Cu as a function of energy: (A) electronic stopping power for H and (B) for He, (C) electronic energy-loss straggling for H and (D) for He, (E) inelastic cross section (ICS) for electrons and (F) stopping power for electrons. Red solid and blue dotted lines correspond to calculations from the dielectric formalism with the MELF-GOS model fit to the experimental optical ELF from Refs. (Palik and Ghosh, 1999; Werner et al., 2009), respectively. A comparison with experimental data (symbols) is made for ions stopping power (International Atomic Energy Agency Nuclear Data Services, 2021) (the insets show the corresponding stopping power around its maximum value for the most recent experimental data (Kido and Hioki, 1983; Semrad et al., 1983;
The energy-loss straggling of H and He beams in Cu appears in Figures 7C,D. Red solid lines depict our calculated results based on the experimental optical ELF by Palik and Ghosh (Palik and Ghosh, 1999), while dotted blue lines correspond to Werner et al.’s ELF (Werner et al., 2009). Experimental data are represented by symbols (Hoffman and Powers, 1976; Nomura et al., 1976; Friedland and Kotze, 1981; Kido, 1987; Kawano and Kido, 1988; Kido and Koshikawa, 1991;
Figure 7E shows the results for the ICS of electrons in Cu, calculated using the ELF from optical data (red lines) (Palik and Ghosh, 1999) and that from REELS (blue lines) (Werner et al., 2009), in comparison with a relatively wide collection of experimental data (symbols). In this case, there are experiments performed for energies larger than 50 eV (Tanuma et al., 2005;
Our calculations (as well as the ab initio results from Ref (
Results for the stopping power of electrons in Cu are given in Figure 7F. Here, pure experimental data (squares) are only available from Al-Ahmad et al. at high energies (
Finally, Figure 3 shows by lines our calculated energy-loss straggling for electrons in Cu, for the two sets of ELF, as indicated in the figure caption which, as for the other cases, cannot be compared with any other experiment or reference calculation.
5.4 Molybdenum
The next metal we analyse is molybdenum, for which, as in the case of copper, there exist very low-energy electron measurements of the ICS by XAFS (
FIGURE 8

(A) and (B): Optical energy-loss function, ELF (k = 0, E), of Mo as a function of the energy transfer E. The experimental data from Palik and Ghosh’s compilation (Palik and Ghosh, 1999) (magenta circles), from Werner et al. (Werner et al., 2009) (square blue symbols) and from Pauly et al. (Pauly et al., 2020) (red triangle symbols) are presented as well as the MELF-GOS fit to these experimental data. Calculations from Henke et al. (Henke et al., 1993) and the NIST (
The MELF-GOS fittings to these sets of data are shown in Figures 8A,B by lines, and the corresponding parameters are gathered in Supplementary Tables S6–S8 in the Supplementary Material. The effective number of electrons Neff in Mo (calculated from the MELF-GOS model) as a function of the excitation energy E is depicted in Figures 8C,D. No important differences appear in the evolution of Neff with the excitation energy, except for a slight underestimation of the number of outer-shell electrons below 100–300 eV for Pauly et al.’s data (Pauly et al., 2020). The grey dashed lines correspond to 42 (the total number of electrons), 40 (where the K-shell starts), 32 (where the L-shell starts) and 14 (where the M-shell starts). Table 7 indicates that the total f- and KK-sum rules are fulfilled for the three fittings when the excitation energy tends to infinity. The effective number of electrons for each electronic excitation is reasonably satisfied, with Pauly et al.’s data (Pauly et al., 2020) slightly better respecting the expected number of 5s and 4d electrons. Noticeably, for Mo a vast majority of the outermost electrons participate in the plasmon excitation slightly above 20 eV. Another noticeable difference is observed for the excitation of the 4p electrons above 35 eV, which is underestimated by Pauly et al.’s data (Pauly et al., 2020), hence giving the lower number of outer-shell electrons with respect to Palik and Ghosh’s (Palik and Ghosh, 1999) (combined with Henke’s (Henke et al., 1993)) and Werner et al.’s (Werner et al., 2009) datasets. Regarding the mean excitation energies I obtained from these optical ELFs, they have values of 443.4 eV when Werner et al.’s ELF (Werner et al., 2009) is used, 469.2 eV for Palik and Ghosh’s (Palik and Ghosh, 1999) and 503.9 eV for Pauly et al.’s (Pauly et al., 2020), while the recommended values by ICRU37 and ICRU49 are, respctively, 424 ± 15 eV (
TABLE 7
| Target Mo (Palik and Ghosh, 1999) | Target Mo (Pauly et al., 2020) | Target Mo (Werner et al., 2009) | |||||||
|---|---|---|---|---|---|---|---|---|---|
| Level | i-th MELF | Nexpect | Neff | Neff | Neff | ||||
| 5s | 1–2 | 2 | 6 | 0.32 | 4.43 | 0.50 | 5.21 | 0.15 | 4.18 |
| Plasmon | 3 | 3.44 | 3.86 | 3.78 | |||||
| 4d | 4 | 4 | 0.67 | 0.84 | 0.25 | ||||
| 4p | 5–6 | 6 | 8.38 | 5.52 | 8.59 | ||||
| 4s | 7 | 2 | 2.09 | 2.08 | 2.08 | ||||
| 3d | 8 | 10 | 11.27 | 12.47 | 11.28 | ||||
| 3p | 9 | 6 | 6.83 | 7.59 | 6.83 | ||||
| 3s | 10 | 2 | 2.18 | 2.18 | 2.18 | ||||
| 2p | L-shell GOS | 6 | 4.39 | 4.39 | 4.39 | ||||
| 2s | L-shell GOS | 2 | 1.35 | 1.35 | 1.35 | ||||
| 1s | K-shell GOS | 2 | 1.17 | 1.17 | 1.17 | ||||
| Total | 42 | 42.08 | ɛrel = 0.19% | 41.95 | ɛrel = −0.11% | 42.04 | ɛrel = 0.10% f-sum rule | ||
| 1.00 | ɛrel = −0.39% | 1.02 | ɛrel = 1.74% | 1.02 | ɛrel = 1.76% KK-sum rule |
Comparison of the expected number of electrons Nexpect in each excitation level of molybdenum, and the corresponding effective number of electrons Neff, as obtained from Eq. (10) with the MELF-GOS model fit to the experimental optical ELF from Palik and Ghosh (Palik and Ghosh, 1999) and Werner et al. (Werner et al., 2009). The values of the f- and KK-sum rules appear in the last two rows.
The stopping powers for H and He ions in Mo obtained by means of the MELF-GOS model from the different optical ELFs are represented by lines in Figures 9A,B respectively. Results obtained from Palik and Ghosh’s data (Palik and Ghosh, 1999) are shown by solid red lines, from Werner et al.’s data (Werner et al., 2009) by blue dotted lines and from Pauly et al.’s experiments (Pauly et al., 2020) by orange dashed lines. All calculations agree at high and low energies but are different around the maximum. The experimental data available from Paul’s database (International Atomic Energy Agency Nuclear Data Services, 2021), which are rather scarce for Mo, are depicted by letters. We note that the most recent measurements for H by Moro et al. (Moro et al., 2016) (yellow symbols) coincide almost perfectly with the calculated stopping powers based on Palik and Ghosh’s (Palik and Ghosh, 1999) and Werner et al.’s (Werner et al., 2009) ELFs, while the use of Pauly et al.’s data (Pauly et al., 2020) underestimates them. In the insets, the most recent measurements around the maximum are presented. In the case of H in Mo we include the data from Refs. (Izmailov et al., 1980; Sirotinin et al., 1984), whereas for He we represent measurements from Refs. (
FIGURE 9

Energy-loss quantities for charged particles in Mo as a function of energy: (A) electronic stopping power for H and (B) for He, (C) electronic energy-loss straggling for H and (D) for He, (E) inelastic cross section (ICS) for electrons and (F) stopping power for electrons. Red solid and blue dotted lines correspond to calculations from the dielectric formalism with the MELF-GOS model fit to the experimental optical ELF from Refs. (Palik and Ghosh, 1999; Werner et al., 2009), respectively. A comparison with experimental data (symbols) is made for ions stopping power (International Atomic Energy Agency Nuclear Data Services, 2021) (the insets show the corresponding stopping power around its maximum value for the most recent experimental data (
Figures 9C,D depict the energy-loss straggling of H and He beams in Mo. Red solid and blue dotted lines show our calculated results based either on the experimental optical ELF by Palik and Ghosh (Palik and Ghosh, 1999) or by Werner et al.’s ELF (Werner et al., 2009). In this case, there are no experimental data to compare with.
In what follows, we discuss the energy-loss quantities of electrons in Mo. As we are going to analyse the possible influence of surface excitation effects, as we did for Cu, we will use the effective ELF for Mo obtained from 500 eV-electron REELS measurements by Chen and coworkers (
The inelastic cross section (ICS) of electrons in Mo is presented in Figure 9E. Calculations based on Palik and Ghosh’s (red line) (Palik and Ghosh, 1999), Werner et al.‘s (blue lines) (Werner et al., 2009) and Pauly et al.’s (Pauly et al., 2020) (orange dashed line) optical ELFs are depicted, as well as those obtained from the effective ELF by Chen et al. (
The stopping power of electrons in Mo is shown in Figure 9F. Calculations are derived from the three different experimental optical ELFs (Werner et al., 2009): (blue lines), (Pauly et al., 2020), (dashed orange line), and the effective ELF (
Figure 3 shows by lines our calculated energy-loss straggling for electrons in Mo for the two sets of ELF, as indicated in the figure caption.
5.5 Iron
Iron is a metal with a reasonable amount of experimental data for the energy-loss quantities, both for ions and electrons, and for which three conflicting sources of optical ELF are available, from optical (Palik and Ghosh, 1999) and REELS (Werner et al., 2009; Xu et al., 2017) measurements. These are shown, respectively, by magenta circles, blue squares and orange triangles in Figure 10A. The MELF-GOS fittings are depicted by solid red, dotted blue and dashed orange lines, respectively, with parameters given in Supplementary Tables S9–S11 in the Supplementary Material. As can be seen, REELS-based measurements seem to provide much sharper and intense peaks, while the optical data are broader and in general less intense, with the maxima being around half of the source of data by Werner et al. (Werner et al., 2009). As a consequence, although all optical ELFs respect the total f-sum rule, as well as the KK-sum rule to a good extent, Palik and Ghosh’s data (Palik and Ghosh, 1999) give a somewhat lower number of electrons for excitation energies up to several hundreds of eV, as can be seen in Figure 10C. Total and individual sum-rules are given in Table 8. Individual effective number of electrons are close to expected numbers for all excitation channels for all optical ELFs, although the error in the total KK-sum rule is significantly larger for Palik and Ghosh’s (Palik and Ghosh, 1999) and Xu et al.’s data (Xu et al., 2017) as compared to Werner et al.’s (Werner et al., 2009). Nonetheless, Xu et al.’s ELF (the most recent one, based on a reverse Monte Carlo algorithm using REELS measurements (
FIGURE 10

(A) and (B): Optical energy-loss function, ELF (k = 0, E), of Fe as a function of the energy transfer E. The experimental data from Palik and Ghosh’s compilation (Palik and Ghosh, 1999) (magenta circles), from Werner et al. (Werner et al., 2009) (square blue symbols) and from Xu et al. (Xu et al., 2017) (orange triangles) are presented as well as the MELF-GOS fit to these experimental data. Calculations from Henke et al. (Henke et al., 1993) and the NIST (
TABLE 8
| Target Fe (Palik and Ghosh, 1999) | Target Fe (Werner et al., 2009) | Target Fe (Xu et al., 2017) | ||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|---|
| Level | i-th MELF | Nexpect | Neff | Neff | Neff | |||||||
| 4s | 1–2 | 2 | 8 | 0.42 | 1.94 | 7.95 | 0.37 | 2.93 | 7.74 | 0.43 | 2.68 | 8.41 |
| plasmon | 3 | 1.52 | 2.56 | 2.25 | ||||||||
| 3d | 4–5 (+9 for (Xu et al., 2017)) | 6 | 6.01 | 4.81 | 5.73 | |||||||
| 3p | 6–7 (+9 for (Werner et al., 2009)) | 6 | 6.46 | 6.77 | 6.61 | |||||||
| 3s | 8 | 2 | 2.44 | 2.32 | 1.80 | |||||||
| 2p | L-shell GOS | 6 | 6.06 | 6.06 | 6.06 | |||||||
| 2s | L-shell GOS | 2 | 1.54 | 1.54 | 1.54 | |||||||
| 1s | K-shell GOS | 2 | 1.32 | 1.32 | 1.32 | |||||||
| Total | 26 | 26.01 | ɛrel = 0.05% | 26.00 | ɛrel < −0.01% | 25.99 | ɛrel < −0.03% | f-sum rule | ||||
| 0.92 | ɛrel = −7.89% | 1.03 | ɛrel = 2.79% | 1.07 | ɛrel = 6.91% | KK-sum rule |
Comparison of the expected number of electrons Nexpect in each excitation level of iron, and the corresponding effective number of electrons Neff, as obtained from Eq. (10) with the MELF-GOS model fit to the experimental optical ELF from Palik and Ghosh (Palik and Ghosh, 1999), Werner et al. (Werner et al., 2009) and Xu et al. (Xu et al., 2017). The values of the f- and KK-sum rules appear in the last two rows.
The calculated stopping powers for H and He in Fe are presented by lines (solid red line for Palik and Ghosh’s ELF (Palik and Ghosh, 1999), blue dotted line for Werner et al.’s (Werner et al., 2009), orange dashed line for Xu et al.’s (Xu et al., 2017)) in Figures 11A,B, and compared to the experimental data (letters) compiled in Paul’s database (International Atomic Energy Agency Nuclear Data Services, 2021). Although the amount of experimental information is not very large for Fe, we can clearly see, both for H and for He ions, that Werner et al.’s (Werner et al., 2009) and Xu et al.’s (Xu et al., 2017) optical ELFs gives place to stopping powers in rather good agreement with experiments, while Palik and Ghosh’s ELF (Palik and Ghosh, 1999) leads to an important underestimation around the maximum. The inset of Figure 11A zooms, for protons, the most recent experimental data around the maximum, revealing that the agreement with calculations based on Xu et al.’s ELF coincide very well in the maximum and around, while those based on Werner et al.’s ELF slightly overestimate them for energies larger than 200 keV/u.
FIGURE 11

Energy-loss quantities for charged particles in Fe as a function of energy: (A) electronic stopping power for H and (B) for He, (C) electronic energy-loss straggling for H and (D) for He, (E) inelastic cross section (ICS) for electrons and (F) stopping power for electrons. Red solid, blue dotted and orange dashed lines correspond to calculations from the dielectric formalism with the MELF-GOS model fit to the experimental optical ELF from Refs. (Palik and Ghosh, 1999; Werner et al., 2009; Xu et al., 2017), respectively. A comparison with experimental data (symbols) is made for ions stopping power (International Atomic Energy Agency Nuclear Data Services, 2021) (the insets show the corresponding stopping power around its maximum value for the most recent experimental data) electrons ICS (Pappas et al., 1991; van Dijken et al., 2002; Tanuma et al., 2005; Zdyb et al., 2013;
Figures 11C,D show the energy-loss straggling of H and He beams in Fe. Red solid lines correspond to our calculations based on the optical ELF by Palik and Ghosh (Palik and Ghosh, 1999), while dotted blue lines correspond to Werner et al.’s ELF (Werner et al., 2009) and dashed orange lines to Xu et al.’s data (Xu et al., 2017). For this metal, there are no experimental data to compare with.
Let us now study electrons in Fe. For this material there is no effective ELF available. However, we estimated an effective ELF, equivalent to those measured for Cu and Mo for 500 eV electrons (Zhang et al., 2004;
The calculated ICS for electrons in Fe are shown by lines in Figure 11E, together with experimental data (symbols) by Refs. (Pappas et al., 1991; van Dijken et al., 2002; Tanuma et al., 2005; Zhukov et al., 2006; Zdyb et al., 2013;
The calculated stopping powers of electrons are shown by thick lines in Figure 11F, although in this case there are no experimental data to compare with. Thin curves showing the contribution of different excitations correspond to the values obtained from Palik and Ghosh’s optical data (Palik and Ghosh, 1999). The difference appearing by treating the excitation of the 3d electrons as individual or collective is more clearly seen here, with the slight increase of the stopping around the maximum when they are treated as collective. The stopping power obtained by using Palik and Ghosh’s optical ELF (Palik and Ghosh, 1999) is significantly lower than that using Werner et al.‘s (Werner et al., 2009) and Xu et al.‘s ELFs (Xu et al., 2017). The use of the effective ELF (based on Werner et al.‘s data) produces a maximum value of the stopping power similar to calculations based on Palik and Ghosh’s data, and it is slightly increased for low electron energies below 50 eV.
Figure 3 contains our calculated energy-loss straggling for electrons in Fe for the three sets of ELF (see figure caption).
5.6 Platinum
We finalise our analysis with platinum, a relevant material for catalysis and nanoparticle enhancement of radiotherapy, with again two sets of experimental data available for the optical ELF, coming from optical data (Palik and Ghosh, 1999) and REELS (Werner et al., 2009), which are shown in Figure 12 as a function of the excitation energy E. The optical ELF from Palik and Ghosh’s compilation (Palik and Ghosh, 1999), which is based on reflectivity measurements, is shown by magenta circles, and the most recent optical ELF based on REELS experiments developed by Werner et al. (Werner et al., 2009) is depicted by blue squares. Panel A shows data in linear scale in order to appreciate the complex excitation spectrum, whereas panel B appears in a log-log scale, which is more appropriate at high excitation energies, where inner shells excitation are presented. Atomic calculations by Henke et al. (Henke et al., 1993) (black crosses) and NIST (
FIGURE 12

(A) and (B): Optical ELF (k = 0, E) of Pt as a function of the excitation energy E. Experimental data from Palik and Ghosh (Palik and Ghosh, 1999) and Werner et al. (Werner et al., 2009) are depicted by symbols and the fitting by the MELF-GOS model by lines. Henke et al. (Henke et al., 1993) and the NIST (
Panels C and D in Figure 12 correspond to the effective number Neff of electrons in Pt, as a function of the excitation energy E. The dashed gray lines correspond to 78 (the total number of electrons), 76 (where the K-shell starts), 68 (where the L-shell starts) and 50 (where the M-shell starts). It is immediate to verify that the total f-sum rule is verified and that at high transferred energies the number of electrons approaches the atomic number of the target, 78 electrons, for both sources of the optical ELF. The partial numbers of effective electrons appear in Table 9, where it can be seen that, in general, both optical ELFs fulfill rather well the expected number of electrons. The total KK-sum rule is better satisfied by Werner et al.’s optical ELF (Werner et al., 2009). The I-values obtained are 738.6 eV for Palik and Ghosh’s and 739.0 eV for Werner et al.’s ELFs, while the value recommended by ICRU37 is 790 ± 30 eV (
TABLE 9
| Target Pt (Palik and Ghosh, 1999) | Target Pt (Werner et al., 2009) | ||||||||
|---|---|---|---|---|---|---|---|---|---|
| Level | i-th MELF | Nexpect | Neff | Neff | |||||
| 6s | 1–2 | 2 | 10 | 0.67 | 2.18 | 10.44 | 0.35 | 1.45 | 8.89 |
| plasmon | 3 | 1.51 | 1.10 | ||||||
| 5d | 4–5 | 8 | 8.26 | 7.44 | |||||
| 5p | 6–7 (+13 for (Werner et al., 2009)) | 6 | 6.76 | 7.56 | |||||
| 4f | 8 | 14 | 17.63 | 18.22 | |||||
| 5s | 9 | 2 | 2.40 | 2.40 | |||||
| 4d | 10 | 10 | 15.13 | 15.13 | |||||
| 4p | 11 | 6 | 5.76 | 5.76 | |||||
| 4s | 12 | 2 | 2.35 | 2.35 | |||||
| 3d | M-shell GOS | 10 | 6.93 | 6.93 | |||||
| 3p | M-shell GOS | 6 | 4.73 | 4.73 | |||||
| 3s | M-shell GOS | 2 | 1.30 | 1.30 | |||||
| 2p | L-shell GOS | 6 | 2.70 | 2.70 | |||||
| 2s | L-shell GOS | 2 | 1.01 | 1.01 | |||||
| 1s | K-shell GOS | 2 | 0.86 | 0.86 | |||||
| Total | 78 | 78.08 | ɛrel < 0.01% | 77.94 | ɛrel = −0.08% | f-sum rule | |||
| 1.11 | ɛrel = 10.78% | 1.03 | ɛrel = 3.28% | KK-sum rule | |||||
Effective number Neff of electrons in each excitation level for Pt obtained by the MELF-GOS model fit to the experimental optical ELF from Palik and Ghosh (Palik and Ghosh, 1999) and Werner et al. (Werner et al., 2009).
The electronic stopping powers of Pt for H and He ions are presented in Figures 13A,B as a function of the incident ion energy. The experimental data (letters) are from Paul’s database (International Atomic Energy Agency Nuclear Data Services, 2021). Our calculations based on the ELF from Palik and Ghosh (Palik and Ghosh, 1999) are depicted by a red solid line and the results based on Werner et al.’s ELF (Werner et al., 2009) are shown by a blue dotted line. Both results are very similar for the energy range analysed, with only slight differences below 200 keV/u, and are compatible with the collection of experimental data down to 30–40 keV/u. The insets in the figure compare the calculated stopping powers in Pt with the most recent experimental data around the maximum stopping, from Refs. (Selau et al., 2020; Moro et al., 2020) for H, and from Ref. (Moro et al., 2020) for He. As can be seen, the agreement is excellent, with the calculations within error bars for H and also for He above 100 keV/u. The recommended values of the stopping power from ICRU49 (ICRU, 1993) (cyan dashed line) and from the semiempirical SRIM code (Ziegler, 2013) (green dotted line) are also presented. It should be noted that the values predicted by the SRIM code (Ziegler, 2013) for H in Pt do not agree with the most recent experimental determinations.
FIGURE 13

Energy-loss quantities for charged particles in Pt as a function of energy: (A) electronic stopping power for H and (B) for He, (C) electronic energy-loss straggling for H and (D) for He, (E) inelastic cross section (ICS) for electrons and (F) stopping power for electrons. Red solid and blue dotted lines correspond to calculations from the dielectric formalism with the MELF-GOS model fit to the experimental optical ELF from Refs. (Palik and Ghosh, 1999; Werner et al., 2009), respectively. A comparison with experimental data (symbols) is made for ions stopping power (International Atomic Energy Agency Nuclear Data Services, 2021) [the insets show the corresponding stopping power around its maximum value for the most recent experimental data (Selau et al., 2020; Moro et al., 2020)], ions straggling (Kido and Koshikawa, 1991), electrons ICS (Tanuma et al., 2005; Jablonski, 2016) and electrons stopping power (Luo, 1991; Hovington et al., 1996), as well as with other theoretical models (ICRU, 1993; Ziegler, 2013). For the other lines, please see details on the text.
In Figures 13C,D we show the energy-loss straggling of H and He beams in Pt. Red solid lines and blue dotted lines represent our calculations based on the experimental optical ELF by Palik and Ghosh’s (Palik and Ghosh, 1999) and by Werner et al.’s ELF (Werner et al., 2009) respectively, which are close to the experimental data for protons (Kido and Koshikawa, 1991).
For electron interactions in Pt, experimental data for the ICS are very scarce (being only available at the maximum value and larger energies), while there are no experimental data for the stopping power. The experimental ICS by Tanuma et al. (Tanuma et al., 2005) and Jablonski (Jablonski, 2016) are shown by symbols in Figure 13E, along with our calculations based on Palik and Ghosh’s (red line) (Palik and Ghosh, 1999) and Werner et al.’s (blue line) (Werner et al., 2009) ELFs. In the theoretical calculation, we can treat the 5d-electron excitations as collective (solid thin line) or as individual transitions (dashed thin line). At high electron energies the ICS are rather similar independently of the optical ELF used and they match quite well with the experimental data (Tanuma et al., 2005; Jablonski, 2016). However at low electron energies the ICS calculated from Palik and Ghosh’s ELF (Palik and Ghosh, 1999) is somewhat larger than the values obtained from Werner et al.’s ELF (Werner et al., 2009), but the scattering of the experimental results and the absence of data for lower energies prevents drawing further conclusions. In any case, considering the d electrons as collective excitations again seems to improve agreement with experiments.
The calculated stopping powers are similarly shown by lines in Figure 13F, together with the data by Luo et al. (Luo et al., 1991) and Hovington et al. (Hovington et al., 1996) (symbols) that appear in Joy’s compilation (Joy, 1995), which are dielectric calculations based on EELS measured spectra. Our calculations are very similar to those by Hovington et al. (Hovington et al., 1996) at high energies, and also at lower energies around the maximum when the excitation of the 5d electrons is considered as collective. However, more measurements of the stopping power would be desirable to draw more conclusions.
Our calculated energy-loss straggling for electrons in Pt appears in Figure 3, for the two sets of ELF, as indicated in the figure caption.
6 Summary and conclusion
The dielectric formalism approach for charged particle interaction with condensed-phase metals has been reviewed, and models for the particular treatment of the electronic inelastic scattering of light ion (protons and alpha particles) and electron beams have been explained.
For ions, we considered the electron capture and loss processes with the target (resulting in the dependence of their charge fractions with energy, which affect their stopping), the energy loss due to these charge exchange processes themselves, as well as the effect of the ion’s electron cloud polarisation to the energy loss.
For electrons, especial attention has been payed to the indistinguishability between the primary and target’s electrons in terms of exchange (through the Born-Ochkur exchange term) and how it affects the maximum allowed energy loss in inelastic collisions during individual and collective electronic excitations. A simple low-energy modification to the dielectric formalism (based on a Coulomb-field correction) has been also implemented.
The application of the present methodology for both light ion and electron beams has been tested by analysing and comparing energy-loss quantities for these projectiles in the relevant metals Al, Fe, Mo, Cu, Pt, and Au. These targets are good representatives of metals with complex excitation spectra, featuring both collective plasmon excitations and individual interband transitions. They also count on with a relatively large collection of experimental data to compare with, as well as with different sources for their optical excitation spectra (the ELF at vanishing momentum transfer k = 0), derived either from light or electron scattering experiments. These provide, in some cases, conflicting results, and the same happens for some of the experimental energy-loss quantities, particularly for very low energy electrons. Such discrepancies have been analysed in terms of the theoretical calculations, which use the MELF-GOS methodology to extend the experimental optical ELF to finite momentum transfers, both for the excitation of the outer and the inner-shell electrons, through the Mermin dielectric function and hydrogenic generalised oscillator strengths, respectively. A possible interpretation of the ELF in terms of contributions from the excitation of different electronic bands is proposed, which is based on the approximate fulfillment of partial f-sum rules.
Aluminum is an example of a metal with an excitation spectrum dominated by a very intense plasmon. An excellent agreement of the calculated stopping powers of Al for H and He ions is obtained in comparison with the abundant experimental data compiled in the database by the late Helmut Paul (International Atomic Energy Agency Nuclear Data Services, 2021), down to energies around 40 keV/u, well below the stopping maximum. To achieve these results, it is necessary to include in the model both the projectile’s polarisation and the energy loss due to electron capture and loss processes as, otherwise, the agreement with experiments is only attained for energies larger than 200 keV/u (already above the maximum S). This highlights the role of energy loss due to projectile’s polarisation and electron capture and loss for ions at the energies around the maximum of the stopping, so relevant for practical applications.
For electrons in Al, an excellent agreement with the experimental inelastic cross sections (ICS) available in the whole energy range is obtained, down to energies as low as 5 eV above the Fermi energy. The best agreement is obtained when, for the collective plasmon excitation, the maximum energy loss E+ = T − EF is considered, instead of the limit E+ ∼ T/2 appropriate for individual transitions due to electron indistinghishability. This clarifies a controversy on this point arisen in a series of works by several authors (
In the case of gold, we find two different sets of its optical ELF, derived from Palik and Ghosh’s optical parameters’ compilations (Palik and Ghosh, 1999) and from REELS experiments by Werner et al. (Werner et al., 2009), as it will happen for all the rest of materials. In this case, however, the analysis of the sum rules does not reveal major differences between these two optical ELFs, what is confirmed by the calculation of the energy-loss quantities both for ions and electrons, which are very similar using both sources of data.
For H and He ions in Au, full calculations give excellent agreement with the most recent experimental stopping power data down to ∼ 40 keV/u (or even lower energies for He), despite of the scattering in the abundant experimental measurements.
The comparison with the rich collection of experimental ICS for electrons in gold permits to assess the good performance of the calculations in the entire energy range, down to energies as low as 2 eV above the Fermi energy, and to further test other assumptions in the model. An almost perfect agreement with experiments is got when all excitations (individual transitions) are treated with the maximum energy loss E+ = (T + Ba)/2 (where Ba is the binding energy of the involved band), except for the plasmon and the excitation of 5d electrons. When the latter are considered as collective excitations (and hence E+ = T − EF), an excellent agreement with experiments is found. To treat the excitation of the 5d electrons as collective may be justified on the claim from Verkhovtsev et al. (Verkhovtsev et al., 2015) that these are excited through giant atomic resonances, and from other authors who discuss the presence of multiple plasmons in transition metals due to the excitation of the d-electrons (
Similar arguments can be used for the case of copper. Calculations based on the optical ELF from Palik and Ghosh’s compilation (Palik and Ghosh, 1999) give very good results for H and He ions stopping powers in comparison with experimental data down to a few tens of keV/u. Calculations based on the REELS optical ELF (Werner et al., 2009) also give good results, although apparently slightly off from the most recent experimental data. This may be the result of a slight better fulfilling of the total and partial sum rules for Palik and Ghosh’s data, although both sets of optical experiments seem to be very reasonable in this case.
As for electrons in Cu, the better agreement with most of the experimental ICS is got again considering the excitation of the 3d electrons as collective. Calculations based on both optical ELFs give excellent results. However, a big discrepancy is found when comparing with the most recent XAFS experimental determinations by Bourke and Chantler (
This issue has been here analysed in the context of potential surface excitations. These can be qualitatively taken into account in the calculations by using the effective ELF obtained experimentally by Zhang et al. for 500 eV-electrons in copper (Zhang et al., 2004). The effective ELF is claimed to include a weighted sum of bulk and surface electronic excitations, which can be used to interpret REELS experiments. Notably, when this effective ELF is used for performing the calculation of the electrons ICS in Cu, results almost perfectly matching the XAFS determinations are obtained (
The case of molybdenum serves to further discuss some of the above aspects. In this case, there is the ELF from optical measurements (Palik and Ghosh, 1999) as well as two sets obtained from REELS (Werner et al., 2009; Pauly et al., 2020). While Palik and Ghosh’s (Palik and Ghosh, 1999) and Werner et al.’s (Werner et al., 2009) data give almost identical results for H and He stopping powers (which agree well with the most recent experimental data for H above 200 keV/u and also for He around the maximum), Pauly et al.’s optical ELF (Pauly et al., 2020) gives stopping powers significantly underestimated. This seems to be mainly related to the somewhat week excitation of the 4p electrons starting at ∼ 35 eV, which underestimates the partial number of expected electrons. This remarks the important role of semi-core levels in the stopping power.
The effective ELF for 500 eV-electron REELS has also been measured for Mo (
For iron, there are also conflicting data for the optical ELF coming from optical (Palik and Ghosh, 1999) and electron experiments (Werner et al., 2009; Xu et al., 2017) as, even though their shapes are similar, the ones based on REELS are much sharper and intense. In this case, the calculated stopping powers of Fe for H and He ions are only in reasonable agreement with the scarce experimental data when Werner et al.’s (Werner et al., 2009) or Xu et al.’s (Xu et al.., 2017) ELFs are used, whereas calculations based on Palik and Ghosh’s (Palik and Ghosh, 1999) significantly underestimate S around the maximum. Particularly, Xu et al.’s data (the most recent determination based on the reverse Monte Carlo method) gives results in very good agreement with the ion stopping power experiments. This case highlights the importance of the semi-core excitations and of counting with reliable methods for the ELF determination.
The differences between the optical ELF datasets does not seem to be so critical for calculating the electron ICS, although again Werner et al.’s and Xu et al.’s data (Werner et al., 2009; Xu et al., 2017) seem to give results closer to experiments at high energies. In this case, low-energy experiments (Pappas et al., 1991; Zhukov et al., 2006; Zdyb et al., 2013;
Finally, similar results are obtained for the relevant metal platinum, for which both sets of optical ELF based on optical measurements and REELS also give almost identical results for H and He stopping powers, in excellent agreement with very recent experimental determinations from Refs. (Selau et al., 2020; Moro et al., 2020), both within experimental error bars. Experimental data is very scarce for electron ICS in Pt, but calculations are in reasonable good agreement with them around and above 100 eV.
Overall, the proposed methodology is capable of yielding, within a common theoretical framework, accurate electronic energy-loss quantities for both H and He ions, as well as electrons, in comparison with experimental measurements, in a very wide energy range down to energies as low as ∼ 40 keV/u for light ions and ∼ 2 eV for electrons, sometimes even for lower energies. For that, it is important that energy-loss due to electron capture and loss and projectile electronic cloud polarisation is considered for ions, and that exchange and indistinguishability and low-energy corrections are appropriately accounted for electrons. Unfortunately, there is a lack of experimental stopping powers for electrons, especially around the maximum value, at energies T ∼ 100 eV, data which would be very useful to further check models at the low energy range.
In some cases, it is also important to appropriately chose the experimental optical ELF to be used for calculations. Judging by the comparison between calculated I-values and the ones recommended by ICRU (
For the assessment of optical ELF data, apart from relying on total and partial sum rules (which could give some guidance in their assessment, but are not always so conclusive), one could nowadays also count on with the help of ab initio time-dependent density functional theory calculations (TDDFT). It was already shown by Werner et al. (Werner et al., 2009) that DFT calculations on the random-phase approximation can provide reasonable estimates of the optical ELF. However, more recently it is becoming clear that already established TDDFT codes are able to provide prediction of the ELF that almost perfectly match the experimental data available for some materials (Pedrielli et al., 2021; Taioli et al., 2021; Taioli et al., 2023). These methods might then be very useful to check conflicting sets of optical ELF data which could affect the reliability of the dielectric formalism calculations. Both experimental optical data and ab initio calculations can be also combined to provide recommended optical ELFs in wide energy ranges, as it was done, for example, for Si by Bichsel (
Furthermore, the proposed theoretical method for calculating the electronic cross sections may further benefit from TDDFT calculations. It has been already made evident here how surface excitations can very importantly affect the interaction probabilities for very low energy electrons. This fact, apart from affecting the experimental determination of electron mean free paths, can have important practical applications in contexts in which low energy electrons are produced in the vicinity of metal surfaces, as it is the case of using metal nanoparticles for the enhancement of ion beam cancer therapy (Kuncic and Lacombe, 2018; Kempson, 2021). Ab initio TDDFT calculations could be exploited for understanding the excitation spectrum not only of bulk metals, but also of nanostructured ones, which could then extend the use of the dielectric formalism, already demonstrated very effective for bulk materials, to accurately describe the electronic interactions of ions and (especially low energy) electrons with metallic nanosystems, an endeavour which is currently underway.
We would like to highlight that analytical expressions have been provided to calculate (by fitting) all the energy-loss quantities obtained and discussed in this paper. These formulas and the corresponding parameters are gathered in the Supplementary Material, due to the interest for providing analytical expression to compare with or to be used in simulation codes.
Statements
Data availability statement
The raw data supporting the conclusion of this article will be made available by the authors, under reasonable request.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmats.2023.1249517/full#supplementary-material
References
1
Al-AhmadK. O.WattD. E. (1983). Stopping powers and extrapolated ranges for electrons (1-10 keV) in metals. J. Phys. D - Appl. Phys.16, 2257–2267. 10.1088/0022-3727/16/11/028
2
AlbertsH. W.MalherbeJ. B. (1983). Energy loss and straggling of p, d and alpha-particles in Au in the energy region 0.2 to 2.4 MeV. Radiat. Eff.69, 231–238. 10.1080/00337578308217825
3
AmadonS.LanfordW. A. (2006). He stopping power and straggle in Al, Ti, Co, Cu, Ag, Ta and Au from 1.5 to 4MeV. Nucl. Instrum. Methods Phys. Res. B - Beam Interact. Mater. Atoms249, 34–37. 10.1016/j.nimb.2006.03.017
4
AndersenH.CseteA.IchiokaT.KnudsenH.MøllerS.UggerhøjU. (2002). An apparatus to measure stopping powers for low-energy antiprotons and protons. Nucl. Instrum. Methods Phys. Res. B - Beam Interact. Mater. Atoms194, 217–225. 10.1016/s0168-583x(02)00692-4
5
AristaN. R. (1978). Erergy loss of correlated charges in an electron gas. Phys. Rev. B18, 1–8. 10.1103/physrevb.18.1
6
AshcroftN. W.MerminN. D. (1976). Solid state physics. Fort Worth, USA: Harcourt College Published.
7
AshleyJ. C.TungC. J.AndersonV. E.RitchieR. H. (1975). Inverse mean free path, stopping power, CSDA range, and straggling in aluminum and aluminum oxide for electrons of energy = or > 10 keV - Air Force Cambridge Research Laboratory Report No. AFCRL-TR-75-0583.
8
AzzoliniM.AngelucciM.CiminoR.LarcipreteR.PugnoN. M.TaioliS.et al (2019). Secondary electron emission and yield spectra of metals from Monte Carlo simulations and experiments. J. Phys. Condens. Matter31, 055901. 10.1088/1361-648x/aaf363
9
BauerG. H.AntolakA. J.PontauA. E.MorseD. H.HeikkinenD. W.ProctorI. D. (1989). Proton energy straggling measurements in aluminum, titanium, silver and tungsten foils. Nucl. Instrum. Methods Phys. Res. B - Beam Interact. Mater. Atoms43, 497–501. 10.1016/0168-583x(89)90396-0
10
BauerM.MarienfeldA.AeschlimannM. (2015). Hot electron lifetimes in metals probed by time-resolved two-photon photoemission. Prog. Surf. Sci.90, 319–376. 10.1016/j.progsurf.2015.05.001
11
BauerP.AumayrF.SemradD.ScherzerB. (1984). Measurement of the stopping cross section for protons in copper by backscattering using various methods for foil-thickness determination. Nucl. Instrum. Methods Phys. Res. B - Beam Interact. Mater. Atoms1, 1–8. 10.1016/0168-583x(84)90469-5
12
BennettH. E.SilverM.AshleyE. J. (1963). Infrared reflectance of aluminum evaporated in ultra-high vacuum. J. Opt. Soc. Am.53, 1089–1095. 10.1364/josa.53.001089
13
BergerM. J.InokutiM.AndersonH. H.BichselH.DennisJ. A.PowersD.et al (1984). ICRU Report 37: Stopping powers for electrons and positrons.
14
BergerM. J.SeltzerS. M. (1982). “Stopping powers and ranges of electrons and positrons,” in NBSIR 82-2550 Tech. (Washington, DC: National Bureau of Standards).
15
BesenbacherF.AndersenJ. U.BonderupE. (1980). Straggling in energy loss of energetic hydrogen and helium ions. Nucl. Instrum. Methods168, 1–15. 10.1016/0029-554x(80)91224-0
16
BianconiM.BarradasN. P.CorreraL. (2005). The stopping cross-section of aluminum for He ions. Nucl. Instrum. Methods Phys. Res. B - Beam Interact. Mater. Atoms239, 127–134.
17
BichselH. (1988). Straggling in thin silicon detectors. Rev. Mod. Phys.60, 663–699. 10.1103/revmodphys.60.663
18
BordoloiA. K.AuluckS. (1983). Electronic structure of platinum. J. Phys. F - Metal Phys.13, 2101–2105. 10.1088/0305-4608/13/10/019
19
BourkeJ. D.ChantlerC. T. (2012). Electron energy loss spectra and overestimation of inelastic mean free paths in many-pole models. J. Phys. Chem. A116, 3202–3205. 10.1021/jp210097v
20
BourkeJ. D.ChantlerC. T. (2010). Measurements of electron inelastic mean free paths in materials. Phys. Rev. Lett.104, 206601. 10.1103/physrevlett.104.206601
21
BrandtW. (1982). Effective charges of ions and the stopping power of dense media. Nucl. Instrum. Methods Phys. Res.194, 13–19. 10.1016/0029-554x(82)90482-7
22
BrandtW.SizmannR. (1975). Atomic collisions in solids. New York: Plenum Press, 305–306.
23
BrownA. M.SundararamanR.NarangP.GoddardW. A.AtwaterH. A. (2016). Nonradiative plasmon decay and hot carrier dynamics: effects of phonons, surfaces, and geometry. ACS Nano10, 957–966. 10.1021/acsnano.5b06199
24
CallcottT. A.ArakawaE. T. (1975). Volume and surface photoemission processes from plasmon resonance fields. Phys. Rev. B11 (8), 2750–2758. 10.1103/physrevb.11.2750
25
ChantlerC. T. (2003). Atomic form factors and photoelectric absorption cross‐sections near absorption edges in the soft X‐ray region. AIP Conf. Proc.652, 370–377. 10.1063/1.1536398
26
ChantlerC. T.BourkeJ. D. (2010). X-Ray spectroscopic measurement of photoelectron inelastic mean free paths in molybdenum. J. Phys. Chem. Lett.1, 2422–2427. 10.1021/jz100776h
27
ChantlerC. T.OlsenK.DragosetR. A.ChangJ.KishoreA. R.KotochigovaS. A.et al (2005). X-ray form factor, attenuation and scattering tables. Available at: http://physics.nist.gov/ffast.
28
Chateau-ThierryA.GladieuxA.DelaunayB. (1976). Experimental neutral charge fractions in proton beams emerging from solids. Nucl. Instrum. Methods132, 553–558. 10.1016/0029-554x(76)90791-6
29
ChenR. B.ChangC. P.LeeC. H.LinM. F. (2007). Electron energy loss spectra of finite carbon nanotubes. J. Appl. Phys.101, 114305. 10.1063/1.2737627
30
ChenY. F.KweiC. M. (1996). Electron differential inverse mean free path for surface electron spectroscopy. Surf. Sci.364, 131–140. 10.1016/0039-6028(96)00616-4
31
ChuW.ZieglerJ.MitchellI.MackintoshW. (1973). Energy‐loss measurements of 4He ions in heavy metals. Appl. Phys. Lett.22, 439–437. 10.1063/1.1654703
32
CórdobaR.OrúsP.StrohauerS.TorresT. E.De TeresaJ. M. (2019). Ultra-fast direct growth of metallic micro- and nano-structures by focused ion beam irradiation. Sci. Rep.9, 14076. 10.1038/s41598-019-50411-w
33
DaB.ShinotsukaH.YoshikawaH.DingZ. J.TanumaS. (2014). Extended Mermin method for calculating the electron inelastic mean free path. Phys. Rev. Lett.113, 063201. 10.1103/physrevlett.113.063201
34
DaB.SunY.MaoS. F.ZhangZ. M.JinH.YoshikawaH.et al (2013). A reverse Monte Carlo method for deriving optical constants of solids from reflection electron energy-loss spectroscopy spectra. J. Appl. Phys.113, 214303. 10.1063/1.4809544
35
DaporM. (2022). Aluminum electron energy loss spectra. A comparison between Monte Carlo and experimental data. Front. Mater.9, 1068196. 10.3389/fmats.2022.1068196
36
DaporM. (2020). Transport of energetic electrons in solids. Computer simulation with applications to materials analysis and characterization. Cham, Switzerland: Springer International Publishing AG.
37
de VeraP.AbrilI.Garcia-MolinaR. (2021). Excitation and ionisation cross-sections in condensed-phase biomaterials by electrons down to very low energy: application to liquid water and genetic building blocks. Phys. Chem. Chem. Phys.23, 5079–5095. 10.1039/d0cp04951d
38
de VeraP.Garcia-MolinaR. (2019). Electron inelastic mean free paths in condensed matter down to a few electronvolts. J. Phys. Chem. C123, 2075–2083. 10.1021/acs.jpcc.8b10832
39
de VeraP.SurdutovichE.AbrilI.Garcia-MolinaR.Solov’yovA. V. (2014). Analytical model of ionization and energy deposition by proton beams in subcellular compartments. Eur. Phys. J. D68, 96. 10.1140/epjd/e2014-50041-7
40
DentonC. D.AbrilI.Garcia-MolinaR.Moreno-MarínJ. C.Heredia-AvalosS. (2008a). Influence of the description of the target energy-loss function on the energy loss of swift projectiles. Surf. Interface Analysis40, 1481–1487. 10.1002/sia.2936
41
DentonC. D.AbrilI.Moreno-MarínJ. C.Heredia-AvalosS.Garcia-MolinaR. (2008b). Energy loss of swift H and He projectiles in Al, Si, Ni and Cu targets. Phys. Status Solidi B- Basic Sol. State Phys.245, 1498–1504. 10.1002/pssb.200743283
42
DesmaraisD.DugganJ. L. (1984). An undergraduate α-particle time-of-flight experiment for determining the mean excitation energy for electronic stopping power of Al, Cu, Ag, and Au. Am. J. Phys.52, 408–411. 10.1119/1.13626
43
DingfelderM.HantkeD.InokutiM.ParetzkeH. G. (1998). Electron inelastic-scattering cross sections in liquid water. Radiat. Phys. Chem.53, 1–18. 10.1016/s0969-806x(97)00317-4
44
DitchburnR. W.FreemanG. H. C. (1966). The optical constants of aluminium from 12 to 36 eV, in Proceedings of the Royal Society of London. Series A. Math. Phys. Sci.294, 20–37.
45
DiwanP.KumarS. (2015). d E/d x and range of α-radiations in Al, Ti and Ni metallic foils. Nucl. Instrum. Methods Phys. Res. B - Beam Interact. Mater. Atoms359, 78–84. 10.1016/j.nimb.2015.07.093
46
EckardtJ. C.LantschnerG. H. (2001). Experimental energy straggling of protons in thin solid foils. Nucl. Instrum. Methods Phys. Res. B - Beam Interact. Mater. Atoms175-177, 93–97. 10.1016/s0168-583x(00)00623-6
47
EgertonR. F. (2011). Electron energy-loss spectroscopy in the electron microscope. New York: Springer.
48
EgertonR. F. (2009). Electron energy-loss spectroscopy in the TEM. Rep. Prog. Phys.72, 016502. 10.1088/0034-4885/72/1/016502
49
EmfietzoglouD. (2003). Inelastic cross-sections for electron transport in liquid water: A comparison of dielectric models. Radiat. Phys. Chem.66, 373–385. 10.1016/s0969-806x(02)00504-2
50
EmfietzoglouD.KyriakouI.Garcia-MolinaR.AbrilI. (2013). The effect of static many-body local-field corrections to inelastic electron scattering in condensed media. J. Appl. Phys.114, 144907. 10.1063/1.4824541
51
EmfietzoglouD.KyriakouI.Garcia-MolinaR.AbrilI. (2017). Inelastic mean free path of low-energy electrons in condensed media: beyond the standard models. Surf. Interface Analysis49, 4–10. 10.1002/sia.5878
52
EmfietzoglouD.NikjooH. (2005). The effect of model approximations on single-collision distributions of low-energy electrons in liquid water. Radiat. Res.163, 98–111. 10.1667/rr3281
53
EppacherC.SemradD. (1992). Dependence of proton and helium energy loss in solids upon plasma properties. Nucl. Instrum. Methods Phys. Res. B - Beam Interact. Mater. Atoms69, 33–38. 10.1016/0168-583x(92)95735-a
54
FanoU. (1963). Penetration of protons, alpha particles, and mesons. Annu. Rev. Nucl. Sci.13, 1–66. 10.1146/annurev.ns.13.120163.000245
55
FermiE. (1940). The Ionization Loss of Energy in Gases and in Condensed Materials. Phys. Rev.57, 485–493. 10.1103/PhysRev.57.485
56
Fernández-VareaJ. M.MayolR.LiljequistD.SalvatF. (1993). Inelastic scattering of electrons in solids from a generalized oscillator strength model using optical and photoelectric data. J. Phys. Condens. Matter5, 3593–3610. 10.1088/0953-8984/5/22/011
57
FittingH. J. (1974). Transmission, energy distribution, and SE excitation of fast electrons in thin solid films. Phys. Status Solidi (a)26, 525–535. 10.1002/pssa.2210260216
58
Flores-ManceraM. A.VillarrubiaJ. S.Massillon-JLG. (2020). Electron inelastic mean free paths for LiF, CaF2, Al2O3, and liquid water from 433 keV down to the energy gap. ACS Omega5, 4139–4147. 10.1021/acsomega.9b03872
59
FriedlandE.KotzeC. P. (1981). Energy-loss straggling of protons, deuterons and α-particles in copper. Nucl. Instrum. Methods Phys. Res.191, 490–494. 10.1016/0029-554x(81)91050-8
60
GarberF. W.NakaiM. Y.HarterJ. A.BirkhoffR. D. (1971). Low-energy electron beam studies in thin aluminum foils. J. Appl. Phys.42, 1149–1158. 10.1063/1.1660159
61
Garcia-MolinaR.AbrilI.Heredia-AvalosS.KyriakouI.EmfietzoglouD. (2011). A combined molecular dynamics and Monte Carlo simulation of the spatial distribution of energy deposition by proton beams in liquid water. Phys. Med. Biol.56, 6475–6493. 10.1088/0031-9155/56/19/019
62
Garcia-MolinaR.AbrilI.KyriakouI.EmfietzoglouD. (2012). “Energy loss of swift protons in liquid water: Role of optical data input and extension algorithms,” in Radiation Damage to Biomolecular Systems. García Gómez-TejedorG.FussM. C. (Editors). Dordrecht: Springer Science+Business Media B.V.
63
GergelyG.MenyhardM.GurbanS.TothJ.VargaD. (2004). Experimental measurements of the surface excitation parameters of Cu, Au, Ni, Ag, Ge and Pd based on Si and other reference standard materials. Surf. Interface Analysis36, 1098–1101. 10.1002/sia.1849
64
GibaruQ.InguimbertC.CaronP.RaineM.LambertD.PuechJ. (2021). Geant4 physics processes for microdosimetry and secondary electron emission simulation: extension of MicroElec to very low energies and 11 materials (C, Al, Si, Ti, Ni, Cu, Ge, Ag, W, Kapton and SiO2). Nucl. Instrum. Methods Phys. Res. B - Beam Interact. Mater. Atoms487, 66–77. 10.1016/j.nimb.2020.11.016
65
GoldsteinJ. I.NewburyD. E.MichaelJ. R.RitchieN. W. M.ScottJ. H. J.JoyD. C. (2018). Scanning electron microscopy and X-ray microanalysis. New York: Springer.
66
GrandeP. L.SchiwietzG. (2002). The unitary convolution approximation for heavy ions. Nucl. Instrum. Methods Phys. Res. B - Beam Interact. Mater. Atoms195, 55–63. 10.1016/s0168-583x(01)01164-8
67
HarithM.OsmanW.GaafarN.El-NadiL. (1987). Stopping power measurements for megaelectronvolt 4He+ ions via thin films of metals and binary metal alloys. Thin Solid Films149, 219–224. 10.1016/0040-6090(87)90298-7
68
HeizU.LandmanU. (Editors) (2014). Nanocatalysis (Berlin Heidelberg: Springer-Verlag).
69
HenkeB.GulliksonE.DavisJ. (1993). X-ray interactions: photoabsorption, scattering, transmission, and reflection at E = 50-30,000 eV, Z = 1-92. Atomic Data Nucl. Data Tables54, 181–342. 10.1006/adnd.1993.1013
70
Heredia-AvalosS.AbrilI.DentonC. D.Moreno-MarínJ. C.Garcia-MolinaR. (2007). Target inner-shells contributions to the stopping power and straggling for H and He ions in gold. J. Phys. Condens. Matter19, 466205. 10.1088/0953-8984/19/46/466205
71
Heredia-AvalosS.Garcia-MolinaR.Fernández-VareaJ. M.AbrilI. (2005). Calculated energy loss of swift He, Li, B, and N ions in SiO2, Al2O3, and ZrO2. Phys. Rev. A - Atomic, Mol. Opt. Phys.72, 052902. 10.1103/PhysRevA.72.052902
72
Heredia-AvalosS.Garcia-MolinaR. (2002). Projectile polarization effects in the energy loss of swift ions in solids. Nucl. Instrum. Methods Phys. Res. B - Beam Interact. Mater. Atoms193, 15–19. 10.1016/s0168-583x(02)00720-6
73
Heredia-AvalosS.Garcia-MolinaR. (2007). Reduction of the energy loss of swift molecular ions in solids due to vicinage effects in the charge state. Phys. Rev. A - Atomic, Mol. Opt. Phys.76, 032902. 10.1103/physreva.76.032902
74
HoffmanG. E.PowersD. (1976). Energy straggling ofαparticles in solid materials. Phys. Rev. A - Atomic, Mol. Opt. Phys.13, 2042–2048. 10.1103/physreva.13.2042
75
HovingtonP.JoyD. C.GauvinR.EvansN. (1996). Scanning Microscopy. Unpublished data. For a summary see. Editors JoyD. C.LuoS.GauvinR.HovingtonP.EvansN.10, 653–666.
76
HsuJ. Y.YuY. C.LiangJ. H.ChenK. M.NiuH. (2004). Energy loss of He, Li and B isotopes with MeV energies in Au. Nucl. Instrum. Methods Phys. Res. B - Beam Interact. Mater. Atoms219-220, 251–255. 10.1016/j.nimb.2004.01.063
77
HsuJ.YuY.LiangJ.ChenK. (2005). Experimental stopping forces in aluminum and silver by 3He/4He, 6Li/7Li and 10B/11B ions. Nucl. Instrum. Methods Phys. Res. B - Beam Interact. Mater. Atoms241, 155–159. 10.1016/j.nimb.2005.07.077
78
HuheeyJ. E.KeiterE. A.KeiterR. L. (1993). Inorganic chemistry: Principles of structure and reactivity. New York: Harper Collins College Publishers.
79
HuthM.PorratiF.DobrovolskiyO. V. (2018). Focused electron beam induced deposition meets materials science. Microelectron. Eng.185–186, 9–28. 10.1016/j.mee.2017.10.012
80
IakoubovskiiK.MitsuishiK.NakayamaY.FuruyaK. (2008). Mean free path of inelastic electron scattering in elemental solids and oxides using transmission electron microscopy: atomic number dependent oscillatory behavior. Phys. Rev. B - Condens. Matter Mater. Phys.77, 104102. 10.1103/physrevb.77.104102
81
ICRU (1993). Report 49 - Stopping powers and ranges for protons and alpha particles. Bethesda, Maryland: International Commission on Radiation Units and Measurements.
82
InokutiM. (1971). Inelastic collisions of fast charged particles with atoms and molecules—the Bethe theory revisited. Rev. Mod. Phys.43, 297–347. 10.1103/revmodphys.43.297
83
International Atomic Energy Agency Nuclear Data Services (2021). Electronic stopping power of matter for ions. Graphs, Data, Comments and Programs. Available at: https://www-nds.iaea.org/stopping.
84
IshigureN.MoriC.WatanabeT. (1978). Electron stopping power in aluminum in the energy region from 2 to 10.9 keV. J. Phys. Soc. Jpn.44, 973–978. 10.1143/jpsj.44.973
85
IzmailovSh. Z.SirotininE. I.TulinovA. F. (1980). Energy loss of protons in Si, Ge and Mo. Nucl. Instrum. Methods Phys. Res. B - Beam Interact. Mater. Atoms168, 81–84. 10.1016/0029-554x(80)91235-5
86
JablonskiA. (2016). Analytical theory of elastic electron backscattering from elements, alloys and compounds: comparison with experimental data. J. Electron Spectrosc. Relat. Phenom.206, 24–45. 10.1016/j.elspec.2015.10.006
87
JanesD. B.BatistutaM.DattaS.MellochM. R.AndresR. P.LiuJ.et al (2000). Interface and contact structures for nanoelectronic devices using assemblies of metallic nanoclusters, conjugated organic molecules and chemically stable semiconductor layers. Superlattices Microstruct.27, 555–563. 10.1006/spmi.2000.0882
88
JesseS.BorisevichA. Y.FowlkesJ. D.LupiniA. R.RackP. D.UnocicR. R.et al (2016). Directing matter: toward atomic-scale 3D nanofabrication. ACS Nano10, 5600–5618. 10.1021/acsnano.6b02489
89
JeynesC.BarradasN. P.SzilágyiE. (2012). Accurate determination of quantity of material in thin films by Rutherford backscattering spectrometry. Anal. Chem.84, 6061–6069. 10.1021/ac300904c
90
JoyD. C. (1995). A database on electron-solid interactions. Scanning17, 270–275. 10.1002/sca.4950170501
91
KalilF.StoneW. G.HubellH. H.BirkhoffR. D. (1959). ORNL Report 2731.
92
KanterH. (1970a). Electron mean free path near 2 keV in aluminum. Phys. Rev. B - Condens. Matter Mater. Phys.1 (5), 2357–2358. 10.1103/physrevb.1.2357
93
KanterH. (1970b). Slow-electron mean free paths in aluminum, silver, and gold. Phys. Rev. B - Condens. Matter Mater. Phys.1 (2), 522–536. 10.1103/physrevb.1.522
94
KawanoA.KidoY. (1988). Effects of energy straggling on surface analysis with fast ion beams. J. Appl. Phys.63, 75–79. 10.1063/1.340465
95
KempsonI. (2021). Mechanisms of nanoparticle radiosensitization. Wiley Interdiscip. Rev. - Nanomedicine and Nanobiotechnology13, e1656. 10.1002/wnan.1656
96
KhodyrevV. A.MizgulinV. N.SirotininE. I.TulinovA. F. (1984). Stopping cross sections of 80- to 500-KeV protons in phosphorus compounds. Radiat. Eff.83, 21–37. 10.1080/00337578408215788
97
KidoY. (1986). Energy and Z2 dependences of energy straggling for fast proton beams passing through solids. Phys. Rev. B - Condens. Matter Mater. Phys.34, 73–77. 10.1103/physrevb.34.73
98
KidoY. (1987). Energy straggling for fast proton beams passing through solid materials. Nucl. Instrum. Methods Phys. Res. B - Beam Interact. Mater. Atoms24-25, 347–352. 10.1016/0168-583x(87)90658-6
99
KidoY.HiokiT. (1983). Measurements of energy loss and straggling for fast H+ in metals and their compounds by means of a nuclear resonant reaction. Phys. Rev. B - Condens. Matter Mater. Phys.27 (5), 2667–2673. 10.1103/physrevb.27.2667
100
KidoY.KoshikawaT. (1991). Energy straggling for medium-energy H+ beams penetrating Cu, Ag, and Pt. Phys. Rev. A - Atomic, Mol. Opt. Phys.44, 1759–1767. 10.1103/physreva.44.1759
101
KnappJ. A.HimpselF. J.EastmanD. E. (1979). Experimental energy band dispersions and lifetimes for valence and conduction bands of copper using angle-resolved photoemission. Phys. Rev. B - Condens. Matter Mater. Phys.19, 4952–4964. 10.1103/physrevb.19.4952
102
KuldeepJainA. K.1985Stopping cross sections of He+ ions in bismuth, Nucl. Instrum. Methods Phys. Res. B - Beam Interact. Mater. Atoms9, 259–262. 10.1016/0168-583x(85)90748-7
103
KumarS.DiwanK. P. (2018). Energy loss and range of α-particles in different metallic foils. Radiat. Eff. Defects Solids173, 970–977. 10.1080/10420150.2018.1513002
104
KuncicZ.LacombeS. (2018). Nanoparticle radio-enhancement: principles, progress and application to cancer treatment. Phys. Med. Biol.63, 02TR01. 10.1088/1361-6560/aa99ce
105
KweiC.ChenY. F.TungC.WangJ. (1993). Electron inelastic mean free paths for plasmon excitations and interband transitions. Surf. Sci.293, 202–210. 10.1016/0039-6028(93)90314-a
106
LeminenE.FontellA. (1974). Stopping power of Ti, Mo, Ag, Ta and W for 0.5-1.75 MeV 4He ions. Radiat. Eff.22, 39–44. 10.1080/00337577408232143
107
LesiakB.JablonskiA.ZommerL.KosinskiA. A. N. D.GergelyG. K. A.SulyokA.et al (1996). Proc. 6th Eur. Conf. Appl. Surf. Interface Analysis - ECASIA ’9511, 619–622.
108
LinW. K.OlsonH. G.PowersD. (1973). Alpha-particle stopping cross section of solids from 0.3 to 2.0 MeV. Phys. Rev. B - Condens. Matter Mater. Phys.8 (5), 1881–1888. 10.1103/physrevb.8.1881
109
LindauI.PianettaP.YuK. Y.SpicerW. E. (1976). Determination of the escape depth of photoemitted electrons in gold in the energy range 25-75 eV by use of synchrotron radiation. J. Electron Spectrosc. Relat. Phenom.8, 487–491. 10.1016/0368-2048(76)80036-9
110
LindhardJ. (1954). On the Properties of a Gas of Charged Particles. Danske Matematisk-fysiske Meddelelser28 (8), 1–57.
111
LoefflerJ. S.DuranteM. (2013). Charged particle therapy—Optimization, challenges and future directions. Clin. Oncol.10, 411–424. 10.1038/nrclinonc.2013.79
112
LuoS.ZhangX.JoyD. C. (1991). Experimental determinations of electron stopping power at low energies. Radiat. Eff. Defects Solids117, 235–242. 10.1080/10420159108220619
113
MajackijV. D.PucherovN. N. (1988). Ukr. Fiz. Zurnal.33, 1285.
114
ManfrinatoV. R.WenJ.ZhangL.YangY.HobbsR. G.BakerB.et al (2014). Determining the resolution limits of electron-beam lithography: direct measurement of the point-spread function. Nano Lett.14, 4406–4412. 10.1021/nl5013773
115
Martínez-TamayoG.EckardtJ. C.LantschnerG. H.AristaN. R. (1996). Energy loss of H+ and He+ in Al, Zn, and Au in the very low-to intermediate-energy range. Phys. Rev. A - Atomic, Mol. Opt. Phys.54 (4), 3131–3138. 10.1103/physreva.54.3131
116
MathewsonA. G.MyersH. P. (1971). Absolute values of the optical constants of some pure metals. Phys. Scr.4, 291–292. 10.1088/0031-8949/4/6/009
117
MerminN. D. (1970). Lindhard dielectric function in the relaxation-time approximation. Phys. Rev. B - Condens. Matter Mater. Phys.1, 2362–2363. 10.1103/physrevb.1.2362
118
MoiseiwitschB. L.SmithS. J. (1968). Electron impact excitation of atoms. Rev. Mod. Phys.40, 238–353. 10.1103/revmodphys.40.238
119
MøllerS. P.CseteA.IchiokaT.KnudsenH.KristiansenH. P.UggerhøjU. I.et al (2008). Antiproton and proton energy loss straggling at keV energies. Eur. Phys. J. D46, 89–92. 10.1140/epjd/e2007-00314-3
120
MøllerS. P.CseteA.IchiokaT.KnudsenH.UggerhøjU. I.AndersenH. H. (2002). Antiproton stopping at low energies: confirmation of velocity-proportional stopping power. Phys. Rev. Lett.88 (19), 193201. 10.1103/physrevlett.88.193201
121
MontanariC. C.ArchubiC. D.MitnikD. M.MiragliaJ. E. (2009). Energy loss of protons in Au, Pb, and Bi using relativistic wave functions. Phys. Rev. A - Atomic, Mol. Opt. Phys.79, 032903. 10.1103/physreva.79.032903
122
MontanariC. C.MiragliaJ. E. (2013). The energy loss straggling of low Z ions in solids and gases. AIP Conf. Proc.1525, 259–269. 10.1063/1.4802331
123
MontanariC. C.MiragliaJ. E. (2017). Low- and intermediate-energy stopping power of protons and antiprotons in solid targets. Phys. Rev. A - Atomic, Mol. Opt. Phys.96 (1), 012707. 10.1103/physreva.96.012707
124
MontanariC. C.MiragliaJ. E. (2019). Stopping power of protons in transition metals of the groups V and VI. Nucl. Instrum. Methods Phys. Res. B - Beam Interact. Mater. Atoms460, 27–29. 10.1016/j.nimb.2018.10.046
125
MookerjeeS.BeuveM.KhanS. A.ToulemondeM.RoyA. (2008). Sensitivity of ion-induced sputtering to the radial distribution of energy transfers: A molecular dynamics study. Phys. Rev. B - Condens. Matter Mater. Phys.78 (4), 045435. 10.1103/physrevb.78.045435
126
MoroM. V.BauerP.PrimetzhoferD. (2020). Experimental electronic stopping cross section of transition metals for light ions: systematics around the stopping maximum. Phys. Rev. A - Atomic, Mol. Opt. Phys.102, 022808. 10.1103/physreva.102.022808
127
MoroM. V.SilvaT. F.MangiarottiA.Guimarães FilhoZ. O.RizzuttoM. A. (2016). Traceable stopping cross sections of Al and Mo elemental targets for 0.9–3.6-MeV protons. Phys. Rev. A - Atomic, Mol. Opt. Phys.93, 022704. 10.1103/physreva.93.022704
128
MottN. F. (1930). The collision between two electrons, Proceedings of the Royal Society A: Mathematical. Phys. Eng. Sci.126, 259–267. 10.1098/rspa.1930.0006
129
MoussaD.DamacheS.OuichaouiS. (2015). Accurate stopping power measurements for (0.21–2.68)MeV/u 1H+ and 4He+ ions crossing thin Al foils; extraction of the (I, b) parameters. Nucl. Instrum. Methods Phys. Res. B - Beam Interact. Mater. Atoms343, 44–47. 10.1016/j.nimb.2014.11.007
130
NastasiM.MayerJ. W.HirvonenJ. K. (1996). Ion-solid interactions: Fundamentals and applications. Cambridge: Cambridge University Press.
131
Nguyen-TruongH. T. (2013). Determination of the maximum energy loss for electron stopping power calculations and its effect on backscattering electron yield in Monte-Carlo simulations applying continuous slowing-down approximation. J. Appl. Phys.114, 163513. 10.1063/1.4827843
132
Nguyen-TruongH. T. (2017). Electron inelastic mean free path at energies below 100 eV. J. Phys. Condens. Matter29, 215501. 10.1088/1361-648x/aa6b9d
133
NikjooH.UeharaS.EmfietzoglouD. (2012). Interaction of radiation with matter. Boca Raton, FL, USA: CRC Press.
134
NomuraA.MatsubaraF.KiyonoS. (1976). Monte Carlo simulation on energy distribution of low energy protons passing through copper films. J. Appl. Phys.15, 2495–2496. 10.1143/jjap.15.2495
135
OchkurV. I. (1965). Ionization of the hydrogen atom by electron impact with allowance for the exchange. Sov. Phys. - JETP20, 1175–1178.
136
OgawaS.NaganoH.PetekH. (1997). Hot-electron dynamics at Cu(100), Cu(110), and Cu(111) surfaces: Comparison of experiment with Fermi-liquid theory. Phys. Rev. B - Condens. Matter Mater. Phys.55, 10869–10877. 10.1103/physrevb.55.10869
137
OppenheimerJ. R. (1928). On the quantum theory of electronic impacts. Phys. Rev.32 (3), 361–376. 10.1103/physrev.32.361
138
PalikE. D.GhoshG. (1999). The electronic handbook of optical constants of solids. San Diego: Academic Press.
139
PappasD. P.KamperK. P.MillerB. P.HopsterH.FowlerD. E.BrundleC. R.et al (1991). Spin-dependent electron attenuation by transmission through thin ferromagnetic films. Phys. Rev. Lett.66, 504–507. 10.1103/physrevlett.66.504
140
PaulyN.YuberoF.TougaardS. (2020). Optical properties of molybdenum in the ultraviolet and extreme ultraviolet by reflection electron energy loss spectroscopy. Appl. Opt.59, 4527–4532. 10.1364/ao.391014
141
PedrielliA.de VeraP.TrevisanuttoP. E.PugnoN. M.Garcia-MolinaR.AbrilI.et al (2021). Electronic excitation spectra of cerium oxides: from ab initio dielectric response functions to Monte Carlo electron transport simulations. Phys. Chem. Chem. Phys.23, 19173–19187. 10.1039/d1cp01810h
142
PenalbaM.ArnauA.EcheniqueP. M.FloresF.RhR. (1992). Stopping power for protons in aluminum. Europhys. Lett.19, 45–50. 10.1209/0295-5075/19/1/008
143
PoignantF.IpatovA.ChakchirO.LartaudP. J.TestaÉ.GervaisB.et al (2020). Theoretical derivation and benchmarking of cross sections for low-energy electron transport in gold. Eur. Phys. J. Plus135, 358. 10.1140/epjp/s13360-020-00354-3
144
PowellC. J.SteinR. J.NeedhamP. B.DriscollT. J. (1977). Attenuation lengths of low-energy electrons in solids derived from the yield of proton-excited auger electrons: beryllium and aluminum. Phys. Rev. B - Condens. Matter Mater. Phys.16 (4), 1370–1379. 10.1103/physrevb.16.1370
145
PrasadS. S. (1965). Electron exchange in impact ionization of atomic hydrogen. Proc. Phys. Soc.85, 57–59. 10.1088/0370-1328/85/1/308
146
PrimetzhoferD. (2012). Inelastic energy loss of medium energy H and He ions in Au and Pt: Deviations from velocity proportionality. Phys. Rev. B - Condens. Matter Mater. Phys.86, 094102. 10.1103/PhysRevB.86.094102
147
RaetherH. (1965). Solid state excitations by electrons. Berlin, Heidelberg: Springer.
148
RitchieR. H. (1959). Interaction of charged particles with a degenerate Fermi-Dirac electron gas. Phys. Rev.114, 644–654. 10.1103/physrev.114.644
149
RitchieR. H. (1957). Plasma losses by fast electrons in thin films. Phys. Rev.106, 874–881. 10.1103/physrev.106.874
150
RuddM. E.KimY. K.MadisonD. H.GayT. J. (1992). Electron production in proton collisions with atoms and molecules: energy distributions. Rev. Mod. Phys.64, 441–490. 10.1103/revmodphys.64.441
151
RuddM.KimY. K.MärkT.SchouJ.StolterfohtN.ToburenL. (1996). ICRU Report 55 - Secondary electron spectra from charged particle interactions. Bethesda, Maryland: International Commission on Radiation Units and Measurements.
152
RudgeM. R. H. (1965). The scattering of electrons by hydrogen atoms. Proc. Phys. Soc.86, 763–772. 10.1088/0370-1328/86/4/311
153
RudgeM. R. (1968). Theory of the ionization of atoms by electron impact. Rev. Mod. Phys.40, 564–590. 10.1103/revmodphys.40.564
154
SchardtD.ElsässerT.Schulz-ErtnerD. (2010). Heavy-ion tumor therapy: physical and radiobiological benefits. Rev. Mod. Phys.82, 383–425. 10.1103/revmodphys.82.383
155
SchiwietzG.GrandeP. L. (2001). Improved charge-state formulas. Nucl. Instrum. Methods Phys. Res. B - Beam Interact. Mater. Atoms175–177, 125–131. 10.1016/S0168-583X(00)00583-8
156
SchiwietzG.GrandeP. L. (2012). Stopping of protons – improved accuracy of the UCA model. Nucl. Instrum. Methods Phys. Res. B - Beam Interact. Mater. Atoms273, 1–5. 10.1016/j.nimb.2011.07.023
157
SchleifeA.KanaiY.CorreaA. A. (2015). Accurate atomistic first-principles calculations of electronic stopping. Phys. Rev. B - Condens. Matter Mater. Phys.91 (1), 014306. 10.1103/physrevb.91.014306
158
SelauF. F.TrombiniH.MarmittG. G.De AndradeA. M.MoraisJ.GrandeP. L.et al (2020). Stopping and straggling of 60–250-keV backscattered protons on nanometric Pt films. Phys. Rev. A - Atomic, Mol. Opt. Phys.102, 032812. 10.1103/physreva.102.032812
159
SemradD.BauerP.AumayrF.HuberP.ObermannW. (1983). Search for an influence of the measuring method on stopping cross section data near the maximum. Nucl. Instrum. Methods Phys. Res.218, 811–816. 10.1016/0167-5087(83)91088-8
160
SemradD.EppacherE.ToberR. (1990). The stopping power of Ag and Au, with regard to higher-order Z1-effects. Nucl. Instrum. Methods Phys. Res. B - Beam Interact. Mater. Atoms48, 79–82. 10.1016/0168-583x(90)90078-9
161
SemradD.MertensP.BauerP. (1986). Reference proton stopping cross sections for five elements around the maximum. Nucl. Instrum. Methods Phys. Res. B - Beam Interact. Mater. Atoms15, 86–90. 10.1016/0168-583x(86)90259-4
162
ShilesE.SasakiT.InokutiM.SmithD. Y. (1980). Self-consistency and sum-rule tests in the Kramers-Kronig analysis of optical data: applications to aluminum. Phys. Rev. B - Condens. Matter Mater. Phys.22 (4), 1612–1628. 10.1103/physrevb.22.1612
163
Shiomi-TsudaN.SakamotoN.OgawaH.KitobaU. (1999). Stopping powers of Al and Mo for protons from 0.3 to 3.0 MeV. Nucl. Instrum. Methods Phys. Res. B - Beam Interact. Mater. Atoms159, 123–132. 10.1016/s0168-583x(99)00532-7
164
SigmundP. (2006). Particle penetration and radiation effects. General aspects and stopping of swift point charges. Berlin Heidelberg: Springer-Verlag.
165
SigmundP. (2013). “Elements of sputtering theory,” in Nanofabrication by Ion-Beam Sputtering: Fundamentals and applications. SomT.KanjilalD. (Editors). Pan Stanford Publishing.
166
SirotininE.TulinovA.KhodyrevV.MizgulinV. (1984). Proton energy loss in solids. Nucl. Instrum. Methods Phys. Res. B - Beam Interact. Mater. Atoms4, 337–345. 10.1016/0168-583x(84)90577-9
167
SlaterJ. C. (1930). Atomic shielding constants. Phys. Rev.36, 57–64. 10.1103/physrev.36.57
168
SmithD.ShilesE.InokutiM. (1983). “The optical properties and complex dielectric function of metallic aluminum from 0.04 to 104 eV,” Argonne National Laboratory Report ANL-83-24. Argonne, Illinois: Argonne National Laboratory Report.
169
SmithD. Y. (1998). Dispersion theory, sum rules, and their application to the analysis of optical data. San Diego, CA: Academic Press.
170
SmithD. Y.ShilesE. (1978). Finite-energy f-sum rules for valence electrons. Phys. Rev. B - Condens. Matter Mater. Phys.17, 4689–4694. 10.1103/physrevb.17.4689
171
Solov’yovA. V. (Editor) (2017). Nanoscale insights into ion-beam cancer therapy. Cham, Switzerland: Springer International Publishing AG.
172
SzeS. M.MollJ. L.SuganoT. (1963). Range-energy relation of hot electrons in Au. IEEE Trans. Electron Devices10, 336. 10.1109/t-ed.1963.15234
173
TaioliS.TrevisanuttoP. E.de VeraP.SimonucciS.AbrilI.Garcia-MolinaR.et al (2021). Relative role of physical mechanisms on complex biodamage induced by carbon irradiation. J. Phys. Chem. Lett.12, 487–493. 10.1021/acs.jpclett.0c03250
174
TaioliS.DaporM.DimiccoliF.FabiM.FerroniV.GrimaniC.et al (2023). The role of low-energy electrons in the charging process of LISA test masses,Class. Quantum Grav.40, 075001. 10.1088/1361-6382/acbadd
175
TanumaS.PowellC. J.PennD. R. (2005). Calculations of stopping powers of 100 eV to 30 keV electrons in 10 elemental solids. Surf. Interface Analysis37, 978–988. 10.1002/sia.2092
176
TanumaS.PowellC. J.PennD. R. (1993). Use of sum rules on the energy-loss function for the evaluation of experimental optical data. J. Electron Spectrosc. Relat. Phenom.62, 95–109. 10.1016/0368-2048(93)80008-a
177
TracyJ. C. (1974). Abstract: electron escape depths in aluminum. J. Vac. Sci. Technol.11, 280. 10.1116/1.1318597
178
TrzaskaW. H.LyapinV.AlankoT.MuttererM.RäisänenJ.TjurinG.et al (2002). New approach to energy loss measurements. Nucl. Instrum. Methods Phys. Res. B - Beam Interact. Mater. Atoms195, 147–165. 10.1016/s0168-583x(02)01255-7
179
TrzaskaW.KnyazhevaG.PerkowskiJ.AndrzejewskiJ.KhlebnikovS.KozulinE.et al (2018). New experimental stopping power data of 4He, 16O, 40Ar, 48Ca and 84Kr projectiles in different solid materials. Nucl. Instrum. Methods Phys. Res. B - Beam Interact. Mater. Atoms418, 1–12. 10.1016/j.nimb.2017.12.025
180
UdalagamaC.BettiolA. A.WattF. (2009). Stochastic spatial energy deposition profiles for MeV protons and keV electrons. Phys. Rev. B - Condens. Matter Mater. Phys.80, 224107. 10.1103/physrevb.80.224107
181
UtkeI.HoffmannP.MelngailisJ. (2008). Gas-assisted focused electron beam and ion beam processing and fabrication. J. Vac. Sci. Technol. B - Microelectron. Nanom. Struct.26, 1276–1197. 10.1116/1.2955728
182
van DijkenS.JiangX.ParkinS. S. P. (2002). Spin-dependent hot electron transport in Ni81Fe19 and Co84Fe16 films on GaAs(001). Phys. Rev. B - Condens. Matter Mater. Phys.66, 094417. 10.1103/PhysRevB.66.094417
183
VerkhovtsevA. V.KorolA. V.Solov’yovA. V. (2015). Electron production by sensitizing gold nanoparticles irradiated by fast ions. J. Phys. Chem. C119, 11000–11013. 10.1021/jp511419n
184
VriensL. (1966). Binary-encounter electron-atom collision theory. Phys. Rev.141, 88–92. 10.1103/physrev.141.88
185
WeaverJ. H.LynchD. W.OlsonC. G. (1974). Optical properties of V, Ta, and Mo from 0.1 to 35 eV. Phys. Rev. B - Condens. Matter Mater. Phys.10, 501–516. 10.1103/physrevb.10.501
186
WernerW. S.GlantschnigK.Ambrosch-DraxlC. (2009). Optical constants and inelastic electron-scattering data for 17 elemental metals. J. Phys. Chem. Reference Data38, 1013–1092. 10.1063/1.3243762
187
WindtD. L.CashW. C.ScottM.ArendtP.NewnamB.FisherR. F.et al (1988). Optical constants for thin films of Ti, Zr, Nb, Mo, Ru, Rh, Pd, Ag, Hf, Ta, W, Re, Ir, Os, Pt, and Au from 24 Å to 1216 Å. Appl. Opt.27, 246–278. 10.1364/ao.27.000246
188
XuH.DaB.TóthJ.TokésiK.DingZ. J. (2017). Absolute determination of optical constants by reflection electron energy loss spectroscopy. Phys. Rev. B - Condens. Matter Mater. Phys.95, 195417. 10.1103/physrevb.95.195417
189
YangQ.O’ConnorD. J.WangZ. (1991). Empirical formulae for energy loss straggling of ions in matter. Nucl. Instrum. Methods Phys. Res. B - Beam Interact. Mater. Atoms61, 149–155. 10.1016/0168-583x(91)95454-l
190
ZdybR.MenteşT. O.LocatelliA.NiñoM. A.BauerE. (2013). Inelastic mean free path from reflectivity of slow electrons. Phys. Rev. B - Condens. Matter Mater. Phys.87, 075436. 10.1103/physrevb.87.075436
191
ZhangY.WeberW. J. (2009). Response of materials to single ion events. Nucl. Instrum. Methods Phys. Res. B - Beam Interact. Mater. Atoms267, 1705–1712. 10.1016/j.nimb.2009.01.104
192
ZhangY.WeberW. J.WangC. M. (2004). Electronic stopping powers in silicon carbide. Phys. Rev. B - Condens. Matter Mater. Phys.69, 205201. 10.1103/physrevb.69.205201
193
ZhangY.WeberW. J.RazpetA.PossnertO. (2005). Electronic stopping powers for He, Be and F ions in Au. Nucl. Instrum. Methods Phys. Res. B - Beam Interact. Mater. Atoms227, 479–484. 10.1016/j.nimb.2004.10.082
194
ZhukovV. P.ChulkovE. V.EcheniqueP. M. (2006). Lifetimes and inelastic mean free path of low-energy excited electrons in Fe, Ni, Pt, and Au: Ab initio GW+T calculations. Phys. Rev. B - Condens. Matter Mater. Phys.73, 125105. 10.1103/PhysRevB.73.125105
195
ZieglerJ. F. (2013). SRIM - The Stopping and Range of Ions in Matter. Available at: http://www.srim.org/.
Summary
Keywords
ion beams, electron beams, stopping power, energy-loss straggling, electronic cross section, metals, dielectric formalism, energy-loss function
Citation
de Vera P, Abril I and Garcia-Molina R (2023) Electronic cross section, stopping power and energy-loss straggling of metals for swift protons, alpha particles and electrons. Front. Mater. 10:1249517. doi: 10.3389/fmats.2023.1249517
Received
28 June 2023
Accepted
28 August 2023
Published
06 December 2023
Volume
10 - 2023
Edited by
Shijun Zhao, City University of Hong Kong, Hong Kong SAR, China
Reviewed by
Karoly Tokesi, Institute for Nuclear Research (MTA), Hungary
Zine El Abidine Chaoui, University Ferhat Abbas of Setif, Algeria
Updates

Check for updates
Copyright
© 2023 de Vera, Abril and Garcia-Molina.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Pablo de Vera, pablo.vera@um.es
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.