ORIGINAL RESEARCH article

Front. Mater., 06 December 2023

Sec. Computational Materials Science

Volume 10 - 2023 | https://doi.org/10.3389/fmats.2023.1249517

Electronic cross section, stopping power and energy-loss straggling of metals for swift protons, alpha particles and electrons

  • 1. Centro de Investigación en Óptica y Nanofísica, Departamento de Física, Universidad de Murcia, Murcia, Spain

  • 2. Departament de Física Aplicada, Universitat d’Alacant, Alacant, Spain

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(EEth,i) is a switching function which can be either a simple Heaviside step Θ(EEth,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 σij are the cross sections for electron loss or capture, changing from the charge state i to the j, and ΔTij is the energy loss associated to these processes. When charge equilibrium is reached, ϕq(T) σqq+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 TEF, 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 = EBa, 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

TargetAlFeCuMoPtAu
Zt132629427879
At26.9855.8563.5595.94195.09196.97
ρt (g/cm3)2.707.878.9610.221.419.3
EF (eV)11.7 ()11.1 ()7 ()5.83 (Nguyen-Truong, 2017)8.78 ()5.53 ()
Epl (eV) ()15.023.019.325.222.624.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

TargetAlFeCuMoPtAu
Outer levels3s2 3p14s2 3d63d10 4s14d5 5s15d9 6s15d10 6s1
2s2 2p63s2 3p63s2 3p64s2 4p65s2 5p6 4f145s2 5p6 4f14
3s2 3p6 3d104s2 4p64s2 4p6 4d10
Inner shells1s22s2 2p62s2 2p62s2 2p63s2 3p6 3d103s2 3p6 3d10
1s21s21s22s2 2p62s2 2p6
1s21s2

Outer and inner excitation levels of the targets described by the MELF-GOS model.

TABLE 3

TargetAlFeCuMoPtAu
Palik and Ghosh’s compilation (Palik and Ghosh, 1999)165.1329.7347.5469.2738.6761.6
Werner et al. (Werner et al., 2009)283.3309.9443.4739.0775.4
Xu et al. (Xu et al., 2017)294.0
Pauly et al. (Pauly et al., 2020)503.9
ICRU Report 37 ()166 ± 2286 ± 9322 ± 10424 ± 15790 ± 30790 ± 30
ICRU Report 49 (ICRU, 1993)158.3291.1326.3424.0751.1741.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

) are also included. (C) and (D): effective number of electrons in Al as a function of the excitation energy E, as calculated from both datasets.

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)
Leveli-th MELFNexpectNeff
plasmon1 112.76 3.29 11.41
3p+3s230.53
2p+2s388.12
1sK-shell GOS21.59
Total1313.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 (). This is a general feature for all materials, which makes the number of electrons in the outer shells always slightly larger than expected. Overall, each assumed transition contains nearly the expected number of electrons, and the total f-sum rule is fulfilled within a negligible relative error (ɛrel = 0.07%). The error in the KK-sum rule is larger (ɛrel = 6.47%), but still within a tolerable level. Furthermore, for the mean excitation energy we obtain a value of I = 165.1 eV, which is in reasonable agreement with the values appearing in ICRU Reports 37 () and 49 (ICRU, 1993) (see Table 3). All the previous validates the experimental optical ELF, as well as the partial excitation assignment in terms of the MELF-GOS model. Please note that the ELF data obtained in this work go up to excitation energies of 30 keV. We have checked that the integrals over the ELF, at this upper energy range, increase around 10–5% per energy step, while at the bottom of the range these increments can be as high as 100% or more, so at 30 keV all integrals are converged, as can be seen in Figure 1D. The same applies to the rest of metals analysed in this work.

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; ; Moussa et al., 2015; Trzaska et al., 2018), where it can be seen that our calculations are within error bars of the experiments above keV/u.

FIGURE 2

; Moussa et al., 2015; Trzaska et al., 2018)), ions straggling (Kido and Hioki, 1983; Kido, 1986; Kido, 1987; ; Yang et al., 1991; Yang et al., 1991; ; ; ; Møller et al., 2008), electrons ICS (Kanter, 1970b; Kanter, 1970a; Tracy, 1974; ; Powell et al., 1977; Lesiak et al., 1996; ) and electrons stopping power (Kalil et al., 1959; Garber et al., 1971; Fitting et al., 1974; Ishigure et al., 1978; ; Hovington et al., 1996; Joy, 1995), as well as with other theoretical models (ICRU, 1993; Grande and Schiwietz, 2002; Schiwietz and Grande, 2012; Ziegler, 2013; Schleife et al., 2015; ; Montanari and Miraglia, 2017). For the other lines, please see details on the text.

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; ; Yang et al., 1991; Yang et al., 1991; ; ; ; Møller et al., 2008), shown by symbols. It can be seen how the calculations follow very well the shape of the experimental data and, in general, also their absolute values, approaching for energies around 104 keV/u. The dispersion of the experimental data is large, but lines go within or close to the point clouds within the whole energy range.

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; ; Powell et al., 1977; Lesiak et al., 1996; ) in the panel E, where also a comparison with the results of ab initio calculations is included (). We obtain an excellent agreement with experiments down to energies as low as 5 eV above the Fermi energy. The contributions of the plasmon and electronic transitions to the ICS are also presented by thin solid lines. Clearly, the plasmon excitation at around 15 eV makes the main contribution to the ICS over the entire energy range, with the contribution from the semi-core 2s and 2p shells slightly increasing the cross section for energies larger than 300 eV.

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; ) is satisfactory, especially at energies keV. At lower energies, we include results by Hovington et al. (Hovington et al., 1996) taken from Joy’s database (Joy, 1995) (assuming that the energies quoted there correspond to TEF), which compare well with our calculations for energies larger than 100 eV and below 30 eV. It should be noted that these latter data, although claimed to be experimental by their authors, were obtained by calculations within the dielectric formalism performed using the excitation spectrum measured by electron energy-loss spectroscopy (Hovington et al., 1996). Therefore, their similitude with our calculations is not surprising.

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 (; ; Nguyen-Truong, 2013; ). While in a previous work of our group we considered the latter to be given by E+T/2 due to indistinguishability () (see Section 4.2), Bourke and Chantler considered that excitations in metals are dominated by collective excitations, thus making this limit too restrictive (). They then assumed that it has to be determined exclusively by Pauly’s exclusion principle as E+ = TEF. In Figure 2E, the solid lines show the calculations assuming the limit E+ = TEF for the plasmon excitation (which dominates the spectrum), while the dashed lines use the limit E+ = T/2. As can be seen for the ICS, calculations reproduce very well the experimental data over the entire energy range when the maximum energy loss for the plasmon excitation is solely determined by Pauly’s exclusion principle, confirming Bourke and Chantler’s hypothesis () for the case of Al. Results for other metals in the following subsections, where interband transitions coexist with the plasmon, will serve to further elucidate this issue. Regarding the stopping power, it should be noted in panel F of Figure 2 that when the maximum energy loss is (incorrectly) considered to be solely determined by exchange, our calculations converge with those by Hovington et al. (Hovington et al., 1996; Joy, 1995), who probably used the same upper limit (E+T/2).

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. () reported theoretical calculations, which are depicted in the figure by a gray dash-dotted line. As can be seen, both curves are rather close, which indicates that these probably provide reasonable estimates of the energy-loss straggling for electrons in Al. It would be desirable to perform experimental measurements of this quantity.

FIGURE 3

).

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) () are also included, which are only valid at large excitations energies.

FIGURE 4

) (gray stars) are included. (C) and (D): Effective number of electrons, Neff, in Au as a function of the excitation energy E, as obtained from both datasets (Henke et al., 1993; ).

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. (), which coincides with the main second peak in the spectrum. The rest of main features of the optical ELF are attributed following the order of growing binding energies for the different atomic shells (). Features below 20 eV are attributed to the excitation of the outermost 6s electrons to the conduction band, those immediately above the plasmon energy to the 5d excitations (as also done for Au by Verkhovtsev et al. (Verkhovtsev et al., 2015)), then the 5p electrons above 45 eV and the 5s above 80 eV. The excitation of the semi-core levels 4s, 4p, 4d and 4f is contained in the broad feature above 100 eV. Please note that different assignements of the peaks in the spectrum of gold have been proposed in other references [see e.g. Ref. (Taioli et al., 2023)], and that the current assignment may be subject to revision in future work. The fitting parameters are given in Supplementary Tables S2, S3 in the Supplementary Material.

TABLE 5

Target Au (Palik and Ghosh, 1999)Target Au (Werner et al., 2009)
Leveli-th MELFNexpectNeffNeff
6s1–31 111.16 2.16 11.630.87 1.78 11.11
plasmon41.000.91
5d5–7109.479.34
5p866.096.09
5s921.742.06
4s+4p+4d+4f103242.2542.25
3dGOS M106.876.87
3pM-shell GOS64.684.68
3sM-shell GOS21.291.29
2pL-shell GOS62.672.67
2sL-shell GOS21.011.01
1sK-shell GOS20.940.94
Total7979.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) () and ICRU49 (I = 741.9 eV) (ICRU, 1993).

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; ; Primetzhofer, 2012) and alpha particles (Trzaska et al., 2002; Hsu et al., 2004; Heredia-Avalos et al., 2007; Zhang et al., 2005; Primetzhofer, 2012; Kumar et al., 2018; Trzaska et al., 2018), which are in a fair agreement with our calculations. The results of other widely used codes are also displayed: the recommended values of ICRU49 (ICRU, 1993) are shown by a cyan dashed line, the semiempirical results from SRIM (Ziegler, 2013) by a dotted green line, the theoretical model from Montanari et al. (Montanari et al., 2009) by an orange dash-dot-dotted line, as well as the Casp 5.2 calculations by Grande and Schiwietz (Schiwietz and Grande, 2012) by a magenta dash-dotted line. As for Al, for proton energies similar and larger than the maximum stopping power, our results agree well with those by Montanari et al. (Montanari et al., 2009), although this time with smaller deviations at lower energies as well as in better agreement with the experiments.

FIGURE 5

; Zhang et al., 2005; Primetzhofer, 2012; Primetzhofer, 2012; Kumar et al., 2018; Trzaska et al., 2018)), ions straggling (; ; ; ; Hsu et al., 2004; ; Møller et al., 2008), electrons ICS (Sze et al., 1963; Kanter, 1970b; Lindau et al., 1976; Gergely et al., 2004; Tanuma et al., 2005; ) and electrons stopping power (; Luo et al., 1991), as well as with other theoretical models (Montanari et al., 2009; Schiwietz and Grande, 2012; Ziegler, 2013; ). For the other lines, please see details on the text.

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 (; ; ; ; Hsu et al., 2004; ; Møller et al., 2008), The lines reproduce very well the experimental points, very often within error bars. The experimental datasets are less scattered in this case, supporting the reliability of our calculations.

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; ) and the cyan shaded region depicts results from ab initio calculations (). The dashed red and blue lines in this figure represent our calculations as described in Section 4.2, using Palik and Ghosh’s (Palik and Ghosh, 1999) and Werner et al.’s (Werner et al., 2009) optical ELFs, respectively; thin lines depict the contributions from the different excitation channels. In this case, despite the plasmon has a noticeable contribution, the ICS is dominated by the excitation of the 6s electrons at low energies and by the 5d electrons around the maximum and larger energies. These calculations (which are not very dependent on the optical ELF used for Au) agree well with the experimental data for energies eV and eV, but underestimate the maximum around 100 eV. In Ref. (), it was argued that the excitation of the 5d electrons in transition metals may be of collective character. This was justified by the claim by Verkhovtsev et al. that the 5d excitations in Au are a result of atomic giant resonances (Verkhovtsev et al., 2015). Other authors also comment about the presence of multiple plasmons in transition metals due to the excitation of the d-block electrons (; Iakoubovskii et al., 2008). If the 5d excitations are considered as collective, the maximum allowed energy loss grows (as discussed in Sections 4.2 and 5.1) and their contribution (depicted by a thin dashed purple line in Figure 5E if considered as interband transitions) grows to the solid thin purple line. The total ICS are now represented by the thick solid red and blue lines in the figure, which match much better the experimental data around the maximum, including its characteristic shape. Therefore, these results seem to support the interpretation of collective character of the 5d electron excitations. Moreover, the consideration of plasmon and intra- and interband transitions (taking into account appropriately their collective or individual character) provides an almost perfect reproduction of the experimental ICS of electrons in Au, a for which there exist empirical data over almost the entire energy range.

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 (; Luo et al., 1991). As in the case of the ICS, red and blue lines correspond to calculations using Palik and Ghosh’s (Palik and Ghosh, 1999) or Werner et al.’s (Werner et al., 2009) optical ELFs, considering (solid) and not considering (dashed) the 5d-electron excitations as collective. The green line corresponds to the calculation without including exchange and low-energy corrections. For the stopping power, these low-energy corrections are not as crucial as for the ICS (as stopping power arises mainly from the excitation of higher-energy shells), but the individual or collective character of the 5d-electron excitations changes significantly the stopping power for energies ≤500–1,000 eV. The agreement with the experimental data by Al-Ahmad et al. () at high energies ( keV) is very good. Again, the data from Luo’s database (Luo et al., 1991) are not to be really considered experimental, but rather dielectric formalism calculations based on experimentally-determined EELS data. Our results agree well with these data for T ≥ 100–200 eV when the 5d-electron excitations are not considered as collective, but overestimate Luo’s data when they are.

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) () calculations based on atomic data are also included, which are valid at large excitation energies.

FIGURE 6

) are included. (C) and (D): Effective number of electrons in Cu as a function of the excitation energy E, obtained from the MELF-GOS model fitted to Palik and Ghosh (Palik and Ghosh, 1999) (red solid line), to Werner et al. (Werner et al., 2009) (blue dotted line) and to the effective ELF from Zhang (Zhang et al., 2004) (dash-dotted green line). (E) Effective optical ELF (k = 0) of Cu as a function of the excitation energy E derived from experimental REELS measurements (solid dark cyan line) (Zhang et al., 2004), and surface ELF (dashed orange line) (Zhang et al., 2004). The MELF-GOS model fitted to the effective ELF is presented by a dotted green line and other symbols are as in panels (A) and (B).

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 () (second major feature in the ELF). Excitation of the 4s electrons is at lower energies, and of 3d, 3p, and 3s electrons at larger energies in order of ascending binding energies. Inner-shell electrons 1s, 2s, and 2p are described by K and L GOS. The fitting parameters can be found in Supplementary Tables S3, S4 in the Supplementary Material.

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 () and ICRU49 (ICRU, 1993), see Table 3.

TABLE 6

Target Cu (Palik and Ghosh, 1999)Target Cu (Werner et al., 2009)
Leveli-th MELFNexpectNeffNeff
4s11 110.49 1.06 11.260.27 1.18 11.97
plasmon20.570.91
3d3–51010.2010.79
3p666.435.79
3s722.422.41
2pL-shell GOS66.066.06
2sL-shell GOS21.541.54
1sK-shell GOS21.291.29
Total2929.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; ; Khodyrev et al., 1984; Sirotinin et al., 1984; Semrad et al., 1986) and He (; Kuldeep and Jain, 1985; Harith et al., 1987; Majackij and Pucherov, 1988). In general, the calculations based on optical data seem to coincide somewhat better than those based on REELS data (which seems to go in line with the somewhat better fulfillment of the sum rules), although the differences are not very significant in any case. We also compare with the recommended values of ICRU49 (ICRU, 1993) (cyan dash line), the semiempirical SRIM2013 code (Ziegler, 2013) (green short-dotted line) and the theoretical model by Montanari and Miraglia (Montanari and Miraglia, 2013).

FIGURE 7

; ; Khodyrev et al., 1984; Sirotinin et al., 1984; Kuldeep and Jain, 1985; Semrad et al., 1986; Harith et al., 1987; Majackij and Pucherov, 1988)), ions straggling (Hoffman and Powers, 1976; Nomura et al., 1976; Friedland and Kotze, 1981; Kido, 1987; Kawano and Kido, 1988; Kido and Koshikawa, 1991; ), electrons ICS (Sze et al., 1963; Kanter, 1970b; Lindau et al., 1976; Gergely et al., 2004; Tanuma et al., 2005; ) and electrons stopping power (; Luo et al., 1991), as well as with other theoretical models (ICRU, 1993; Montanari and Miraglia, 2013; Ziegler, 2013; ). For the other lines, please see details on the text.

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; ). The lines agree well with the experimental data for protons, and reasonably reproduce the few data for helium.

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; ), and for lower energies there is a discrepancy between the classical experiments (Knapp et al., 1979; Ogawa et al., 1997; ) and the most recent ones obtained from XAFS measurements (), which differ in more than one order of magnitude. Ab initio calculations are also shown by a cyan shaded area () which, although very noisy, seem to correlate better with the classical experimental determinations. Dashed lines treat the 3d-electron excitations as individual, while solid lines as collective. Both sets of optical ELF give ICS which agree very well with the experimental data for energies eV when the 3d electrons are treated as collective. Particularly, calculations using Palik and Ghosh’s ELF (Palik and Ghosh, 1999) perfectly agree with the experiments by Tanuma et al. (Tanuma et al., 2005). At low energies, the agreement with the classical experimental data (Knapp et al., 1979; Ogawa et al., 1997; ) is remarkable down to 2 eV above the Fermi energy (notably, the calculations using the REELS ELF (Werner et al., 2009) match perfectly with the data by Knapp et al. (Knapp et al., 1979)).

Our calculations (as well as the ab initio results from Ref ()) seem to support the classical experimental data for the ICS of electrons in bulk Cu. However, the question arises of why there is such a big discrepancy with the most recent XAFS measurements by Bourke and Chantler (). In Ref (). we suggested that these most recent experiments may be affected in some way by surface excitations. In Figure 6E we plot, together with the two already discussed sets of optical ELF for Cu (Palik and Ghosh, 1999; Werner et al., 2009), the surface excitation ELF and the effective ELF (corresponding to 500 eV primary electrons) determined by Zhang et al. from REELS experiments (Zhang et al., 2004). The concept of effective ELF was developed as an empirical approximation to reproduce by simulations the REELS spectra obtained experimentally in conditions where surface excitations affect the energy spectrum of reflected electrons. As can be seen, the effective ELF (which depends on the energy of the primary electrons) mixes the surface and bulk contributions to the ELF in a way that enhances the surface excitations, while also retaining the bulk ones. In order to test the surface excitation hypothesis, we have fitted this effective ELF (corresponding to 500 eV, the lowest energy reported by Zhang et al. (Zhang et al., 2004)) by means of the MELF-GOS methodology and then performed new calculations of the ICS. The results are shown by a green thick solid line in Figure 7E. Strikingly, the ICS obtained using the effective ELF, which weights bulk and surface excitations, matches almost perfectly Bourke and Chantler’s data at low energies (), within error bars down to 5 eV above the Fermi energy. The calculation converges with the high energy data (Tanuma et al., 2005; ), for which surface excitations are deemed to be less relevant. Therefore, these results seem to point out to the contribution of these surface excitations in the most recent XAFS data (), which would explain the discrepancy with the classical experiments (Knapp et al., 1979; Ogawa et al., 1997; ).

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 (), which agree very well with our calculations. Red lines represent results using the optical ELF (Palik and Ghosh, 1999) and blue lines using the REELS ELF (Werner et al., 2009). Solid lines are calculations treating the 3d-electron excitations as collective while dashed lines as individual. The latter coincide well with the results reported by (Luo et al., 1991) and Hovington et al. (Hovington et al., 1996), although it should be remembered that these are rather calculations than experiments, as discussed previously. Surface excitations (accounted by means of the effective ELF (Zhang et al., 2004) and represented by a green solid line) tend to diminish the stopping power around the maximum at eV, while they increase it at low energies eV above the Fermi energy. More experimental measurements of the stopping power of low electrons in Cu would be needed to check these calculations.

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 () that will serve us to check the assumptions made for the case of copper. For this material, there are three sets of experimental data available for the optical ELF, one relying on measured optical data (Palik and Ghosh, 1999) and two based on REELS experiments (Werner et al., 2009; Pauly et al., 2020). These sets of data are depicted in Figures 8A,B: magenta circles denote the experiments compiled by Palik and Ghosh (Palik and Ghosh, 1999), which are very different, above 30 eV, from those from Werner et al. (blue squares) (Werner et al., 2009) and Pauly et al. (red triangles) (Pauly et al., 2020) (these being rather similar to the results reported by Windt et al. (Windt et al., 1988)). In fact, the MELF-GOS fitting to Palik and Ghosh’s compiled data (Palik and Ghosh, 1999) did not fulfill the total f-sum rule. It should be mentioned that the measured data from Weaver (Weaver et al., 1974) goes up to 30 eV, and at larger excitation energies it is just an extrapolation. These experiments may possibly be affected by oxidation, impurities or roughness in the surface of the target, as mentioned in Palik and Ghosh’s book (Palik and Ghosh, 1999). Therefore, above 30 eV we better rely on the atomic information by Henke et al. (black crosses) (Henke et al., 1993), which, together with NIST’s () (gray starts), are also represented, and are more relevant for Figure 8B where a log-log scale is used at large excitation energies.

FIGURE 8

) are included. (C) and (D): Effective number of electrons in Mo as a function of the excitation energy E, obtained from the MELF-GOS model fitted to Palik and Ghosh (Palik and Ghosh, 1999) (red solid line), to Werner et al. (Werner et al., 2009) (blue dotted line), to Pauly et al. (Pauly et al., 2020) and to the effective ELF from Zhang () (dashed green line). (E) Effective optical ELF (k = 0) of Mo as a function of the excitation energy E derived from experimental REELS measurements (solid dark cyan line) (). The MELF-GOS model fitted to the effective ELF is presented by a dotted green line and other symbols are as in panels (A) and (B).

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 () and 424 eV (ICRU, 1993).

TABLE 7

Target Mo (Palik and Ghosh, 1999)Target Mo (Pauly et al., 2020)Target Mo (Werner et al., 2009)
Leveli-th MELFNexpectNeffNeffNeff
5s1–22 60.32 4.430.50 5.210.15 4.18
Plasmon33.443.863.78
4d440.670.840.25
4p5–668.385.528.59
4s722.092.082.08
3d81011.2712.4711.28
3p966.837.596.83
3s1022.182.182.18
2pL-shell GOS64.394.394.39
2sL-shell GOS21.351.351.35
1sK-shell GOS21.171.171.17
Total4242.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. (; Lin et al., 1973; Leminen and Fontell, 1974). For H, all calculations underestimate the scarce experimental data, although those based on Palik and Ghosh’s (Palik and Ghosh, 1999) and Werner et al.’s (Werner et al., 2009) data coincide very well with experiments above 200 keV/u. For the case of He, calculations based on Palik and Ghosh’s and, especially, on Werner et al.’s ELF coincide very well with the data data by Ref. (Leminen and Fontell, 1974). Calculations based on Pauly et al.’s data (Pauly et al., 2020) significantly underestimate all experimental data. However, further measurements will be needed to confirm such trends. A comparison with the semiempirical code SRIM (Ziegler, 2013) is also included (green dotted line), together with the theoretical calculations by Montanari and Miraglia (Montanari and Miraglia, 2019).

FIGURE 9

; Lin et al., 1973; Leminen et al., 1974; Izmailow et al., 1980; Sirotinin et al., 1984), electrons ICS (Tanuma et al., 2005; ; ; Jablonski, 2016), as well as with other theoretical models (Ziegler, 2013; Montanari and Miraglia, 2019). For the other lines, please see details on the text.

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 (), which is shown by a solid line in Figure 8E. The MELF-GOS fitting is shown by a green dotted line, which is used to calculate the energy-loss quantities.

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. () (dotted green line). For the case of Werner et al. ELF, we present contributions from the excitation of each band, including calculations treating the excitation of the 4d electrons as individual transitions (thin blue dashed line) and collective excitations (thin blue solid line). As in the case of Cu, treating the 4d electrons as collective causes an increment of the ICS value around the maximum, although in this case the increase is considerably smaller, as Mo has less d electrons than Cu. Symbols in the figure correspond to experimental data from Tanuma et al. (Tanuma et al., 2005), Chantler and Bourke (), Bauer et al. (), and Jablonski (Jablonski, 2016). The comparison between the calculated and experimental ICS is excellent at energies larger than 100 eV and very good even around the maximum for all three optical ELFs. At very-low electron energies all calculations agree reasonably well with the measurements by Bauer et al. (), but are much lower than those by Chantler and Bourke (). Notably, the latter high accuracy XAFS measurements (), as well as the data by Jablonski around the maximum (Jablonski, 2016), perfectly agree with the results obtained from the effective ELF (), again indicating that they might be affected by surface excitations.

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 () (dotted green line). In all cases the values at high and low electron energies are rather similar, and they only differ around the maximum. Unfortunately, to our knowledge there are no experimental data or other calculations to compare with. Pauly et al.‘s ELF (Pauly et al., 2020) leads to a somewhat lower stopping power around the maximum as compared to the other experimental ELFs, which is similar to that obtained using the effective ELF (). However, the latter ELF gives rise to slightly larger stopping powers at energies lower than 50 eV.

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 (; Xu et al., 2017)) seems to give, overall, a somewhat lower error in the individual sum rule for the electrons described by Mermin functions.

FIGURE 10

) are included. (C) and (D): Effective number of electrons in Fe as a function of the excitation energy E, obtained from the MELF-GOS model fitted to different optical ELF. (E) Effective optical ELF (k = 0) of Fe as a function of the excitation energy (green line), estimated as explained in the text. Other symbols are as in panels (A) and (B).

TABLE 8

Target Fe (Palik and Ghosh, 1999)Target Fe (Werner et al., 2009)Target Fe (Xu et al., 2017)
Leveli-th MELFNexpectNeffNeffNeff
4s1–22 80.42 1.94 7.950.37 2.93 7.740.43 2.68 8.41
plasmon31.522.562.25
3d4–5 (+9 for (Xu et al., 2017))66.014.815.73
3p6–7 (+9 for (Werner et al., 2009))66.466.776.61
3s822.442.321.80
2pL-shell GOS66.066.066.06
2sL-shell GOS21.541.541.54
1sK-shell GOS21.321.321.32
Total2626.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

) as well as with other theoretical models (ICRU, 1993; Montanari and Miraglia, 2013; Ziegler, 2013; ). For the other lines, please see details on the text.

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; ). For that purpose, we determined the ratio of the effective ELF for Cu with respect to its measured surface ELF (Zhang et al., 2004), and applied this ratio to the surface ELF measured by Werner et al. for Fe (Werner et al., 2009). The result is depicted by a solid line in Figure 10F.

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; ). Thick dashed lines are calculations treating the excitation of the 3d electrons as individual transitions, and solid thick lines treat them as collective. As the number of d electrons is not very large in this case, both calculations are not so different and they agree rather well with Tanuma’s experimental data (Tanuma et al., 2005) above 50 eV. The dispersion of the data around the maximum avoids drawing conclusions about the quality of the optical ELFs used, although both Werner et al.’s (Werner et al., 2009) and Xu et al.’s ELFs (Xu et al., 2017) agree better with the experimental ICS above 100 eV. None of the calculations agree with the experimental data at lower energies (Zhukov et al., 2006; Zdyb et al., 2013; ; Pappas et al., 1991), even though for other metals the agreement with the experiments by Bauer () was typically good. The use of the estimated effective ELF somewhat improves the agreement, bringing the calculation (green line) closer to the data by (Pappas et al., 1991; Zdyb et al., 2013; Zhukov et al., 2006), but still underestimating Bauer’s data (). This calls for the measurement of the effective ELF for electrons in Fe, which might help to close this gap.

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 () (grey stars) are also included which are valid at large excitations energies. Lines represent the fit by the MELF-GOS model to Palik and Ghosh’s (Palik and Ghosh, 1999) (red solid line) and to Werner et al.’s ELFs (Werner et al., 2009) (blue dotted line), which parameters are given in Supplementary Tables S12, S13 in the Supplementary Material.

FIGURE 12

) calculations are included. (C) and (D): Effective number of electrons, Neff, in Pt as a function of the excitation energy E.

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 () and by ICRU49 741.9 eV (ICRU, 1993).

TABLE 9

Target Pt (Palik and Ghosh, 1999)Target Pt (Werner et al., 2009)
Leveli-th MELFNexpectNeffNeff
6s1–22 100.67 2.18 10.440.35 1.45 8.89
plasmon31.511.10
5d4–588.267.44
5p6–7 (+13 for (Werner et al., 2009))66.767.56
4f81417.6318.22
5s922.402.40
4d101015.1315.13
4p1165.765.76
4s1222.352.35
3dM-shell GOS106.936.93
3pM-shell GOS64.734.73
3sM-shell GOS21.301.30
2pL-shell GOS62.702.70
2sL-shell GOS21.011.01
1sK-shell GOS20.860.86
Total7878.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

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+ = TEF 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 (; ; Nguyen-Truong, 2013; ), a point which is also further demonstrated for other metals.

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+ = TEF), 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 (; Iakoubovskii et al., 2008). However, here we limit ourselves to remark the almost perfect agreement with experiments when this is considered. Also, it was demonstrated with the case of electrons in gold that, when exchange and low-energy correction terms are not included in the calculations, results become unreliable below ∼ 100 eV, highlighting the importance of these improvements.

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 (), which are off from calculations and previous experiments by more than an order of magnitude for energies below 50 eV.

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 (). This fact seems to point out to a possible influence of surface excitations in these experimental measurements, a point which is further analysed for other targets. The use of the effective ELF also prompts a reduction in the stopping power around the maximum and an increase at energies below 50 eV, which could have relevant implications for situations in which low energy electrons are ejected in the vicinity of surfaces, such as in nanoparticle radiosensitisation.

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 which again there are recent XAFS measurements for its low-energy ICS available (). Full calculations give good results with the high energy ICS measured by Tanuma et al. (Tanuma et al., 2005), as well as with very low energy determinations compiled by Bauer et al. (), while calculations using the effective ELF provide again an excellent comparison with the XAFS measurements () for energies around 100 eV and below, and with the experiments by Jablonski (Jablonski, 2016) around 100 eV and above, remarking the potential impact of surface excitations on the determination of low energy electron ICS.

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; ) disagree with the calculations. However, although there is no effective ELF measured for this material, we estimated it based on the data for Cu (Zhang et al., 2004), again finding a very significant improvement in the agreement with the cited experimental data at low energies.

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 () (see Table 3), in general the optical ELF derived from REELS by Werner et al. (Werner et al., 2009) seem to give a somewhat better agreement for all the metals analysed. Total and partial sum rules are also advisable to perform in order to detect possible inconsistencies. The clearer examples shown in this work are those of the optical ELFs of Mo derived from REELS experiments by Pauly et al. (Pauly et al., 2020) and of Fe obtained from Palik and Ghosh’s compilation (Palik and Ghosh, 1999); the former seems to underestimate the number of electrons expected in the excitation of the 4p-band, while the latter underestimates all outer-shell electrons. This turns out to be important for correctly assessing the stopping power for light ions. By judging the available optical ELFs in the light of the calculation results, for Au, Mo and Pt, both Palik and Ghosh’s (Palik and Ghosh, 1999) and Werner et al.’s (Werner et al., 2009) ELFs give correct and similar energy-loss quantities for both ions and electrons. The differences are small for Cu, with Palik and Ghosh’s data providing somewhat closer agreement with the experimental stopping powers for H and He ions. For Fe, both Xu et al.’s and Werner et al’s ELF give correct stopping powers for H and He around the maximum, the latter slightly overestimating the experiments above 200 keV/u. Both datasets provide a better agreement with experimental electron ICS for energies larger than 100 eV.

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, 22572267. 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, 231238. 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, 3437. 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, 217225. 10.1016/s0168-583x(02)00692-4

  • 5

    AristaN. R. (1978). Erergy loss of correlated charges in an electron gas. Phys. Rev. B18, 18. 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, 497501. 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, 319376. 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, 18. 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, 10891095. 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, 115. 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, 127134.

  • 17

    BichselH. (1988). Straggling in thin silicon detectors. Rev. Mod. Phys.60, 663699. 10.1103/revmodphys.60.663

  • 18

    BordoloiA. K.AuluckS. (1983). Electronic structure of platinum. J. Phys. F - Metal Phys.13, 21012105. 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, 32023205. 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, 1319. 10.1016/0029-554x(82)90482-7

  • 22

    BrandtW.SizmannR. (1975). Atomic collisions in solids. New York: Plenum Press, 305306.

  • 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, 957966. 10.1021/acsnano.5b06199

  • 24

    CallcottT. A.ArakawaE. T. (1975). Volume and surface photoemission processes from plasmon resonance fields. Phys. Rev. B11 (8), 27502758. 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, 370377. 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, 24222427. 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, 553558. 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, 131140. 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, 439437. 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, 50795095. 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, 20752083. 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, 14811487. 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, 14981504. 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, 408411. 10.1119/1.13626

  • 43

    DingfelderM.HantkeD.InokutiM.ParetzkeH. G. (1998). Electron inelastic-scattering cross sections in liquid water. Radiat. Phys. Chem.53, 118. 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, 2037.

  • 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, 7884. 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, 9397. 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, 373385. 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, 410. 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, 98111. 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, 3338. 10.1016/0168-583x(92)95735-a

  • 54

    FanoU. (1963). Penetration of protons, alpha particles, and mesons. Annu. Rev. Nucl. Sci.13, 166. 10.1146/annurev.ns.13.120163.000245

  • 55

    FermiE. (1940). The Ionization Loss of Energy in Gases and in Condensed Materials. Phys. Rev.57, 485493. 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, 35933610. 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, 525535. 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, 41394147. 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, 490494. 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, 11491158. 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, 64756493. 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, 10981101. 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, 6677. 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, 5563. 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, 219224. 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, 181342. 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, 1519. 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, 20422048. 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, 653666.

  • 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, 251255. 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, 155159. 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, 928. 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, 297347. 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, 973978. 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, 8184. 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, 2445. 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, 555563. 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, 56005618. 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, 60616069. 10.1021/ac300904c

  • 90

    JoyD. C. (1995). A database on electron-solid interactions. Scanning17, 270275. 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), 23572358. 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), 522536. 10.1103/physrevb.1.522

  • 94

    KawanoA.KidoY. (1988). Effects of energy straggling on surface analysis with fast ion beams. J. Appl. Phys.63, 7579. 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, 2137. 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, 7377. 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, 347352. 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), 26672673. 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, 17591767. 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, 49524964. 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, 259262. 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, 970977. 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, 202210. 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, 3944. 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, 619622.

  • 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), 18811888. 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, 487491. 10.1016/0368-2048(76)80036-9

  • 110

    LindhardJ. (1954). On the Properties of a Gas of Charged Particles. Danske Matematisk-fysiske Meddelelser28 (8), 157.

  • 111

    LoefflerJ. S.DuranteM. (2013). Charged particle therapy—Optimization, challenges and future directions. Clin. Oncol.10, 411424. 10.1038/nrclinonc.2013.79

  • 112

    LuoS.ZhangX.JoyD. C. (1991). Experimental determinations of electron stopping power at low energies. Radiat. Eff. Defects Solids117, 235242. 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, 44064412. 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), 31313138. 10.1103/physreva.54.3131

  • 116

    MathewsonA. G.MyersH. P. (1971). Absolute values of the optical constants of some pure metals. Phys. Scr.4, 291292. 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, 23622363. 10.1103/physrevb.1.2362

  • 118

    MoiseiwitschB. L.SmithS. J. (1968). Electron impact excitation of atoms. Rev. Mod. Phys.40, 238353. 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, 8992. 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, 259269. 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, 2729. 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, 259267. 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, 4447. 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, 24952496. 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, 11751178.

  • 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, 1086910877. 10.1103/physrevb.55.10869

  • 137

    OppenheimerJ. R. (1928). On the quantum theory of electronic impacts. Phys. Rev.32 (3), 361376. 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, 504507. 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, 45274532. 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, 1917319187. 10.1039/d1cp01810h

  • 142

    PenalbaM.ArnauA.EcheniqueP. M.FloresF.RhR. (1992). Stopping power for protons in aluminum. Europhys. Lett.19, 4550. 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), 13701379. 10.1103/physrevb.16.1370

  • 145

    PrasadS. S. (1965). Electron exchange in impact ionization of atomic hydrogen. Proc. Phys. Soc.85, 5759. 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, 644654. 10.1103/physrev.114.644

  • 149

    RitchieR. H. (1957). Plasma losses by fast electrons in thin films. Phys. Rev.106, 874881. 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, 441490. 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, 763772. 10.1088/0370-1328/86/4/311

  • 153

    RudgeM. R. (1968). Theory of the ionization of atoms by electron impact. Rev. Mod. Phys.40, 564590. 10.1103/revmodphys.40.564

  • 154

    SchardtD.ElsässerT.Schulz-ErtnerD. (2010). Heavy-ion tumor therapy: physical and radiobiological benefits. Rev. Mod. Phys.82, 383425. 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, 125131. 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, 15. 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, 811816. 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, 7982. 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, 8690. 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), 16121628. 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, 123132. 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, 337345. 10.1016/0168-583x(84)90577-9

  • 167

    SlaterJ. C. (1930). Atomic shielding constants. Phys. Rev.36, 5764. 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, 46894694. 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, 487493. 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, 978988. 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, 95109. 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, 147165. 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, 112. 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, 12761197. 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, 1100011013. 10.1021/jp511419n

  • 184

    VriensL. (1966). Binary-encounter electron-atom collision theory. Phys. Rev.141, 8892. 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, 501516. 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, 10131092. 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, 246278. 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, 149155. 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, 17051712. 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, 479484. 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

Copyright

*Correspondence: Pablo de Vera,

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.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics