Evidence for a spin acoustic surface plasmon from inelastic atom scattering

Closed-shell atoms scattered from a metal surface exchange energy and momentum with surface phonons mostly via the interposed surface valence electrons, i.e., via the creation of virtual electron-hole pairs. The latter can then decay into surface phonons via electron-phonon interaction, as well as into acoustic surface plasmons (ASPs). While the first channel is the basis of the current inelastic atom scattering (IAS) surface-phonon spectroscopy, no attempt to observe ASPs with IAS has been made so far. In this study we provide evidence of ASP in Ni(111) with both Ne atom scattering and He atom scattering. While the former measurements confirm and extend so far unexplained data, the latter illustrate the coupling of ASP with phonons inside the surface-projected phonon continuum, leading to a substantial reduction of the ASP velocity and possibly to avoided crossing with the optical surface phonon branches. The analysis is substantiated by a self-consistent calculation of the surface response function to atom collisions and of the first-principle surface-phonon dynamics of Ni(111). It is shown that in Ni(111) ASP originate from the majority-spin Shockley surface state and are therefore collective oscillation of surface electrons with the same spin, i.e. it represents a new kind of collective quasiparticle: a Spin Acoustic Surface Plasmon (SASP).


Observation of a spin acoustic surface plasmon with Ne atom scattering
Ni(111) has actually been the very first metal surface investigated with IAS time-of-flight (TOF) spectroscopy from a metal surface. In these measurements, performed by Feuerbacher and Willis in 1981 with Ne-atom scattering (NeAS) 20 , evidence was obtained, besides the observation of the Rayleigh wave (RW) dispersion curve up to 2/3 of the surface Brillouin zone (SBZ) in the [112] direction, of additional excitations at very small wavevectors and energies extending beyond the maximum of the phonon spectrum. These authors underlined the basic difference between IAS from insulator and metal surfaces, the latter involving surface conduction electrons in all energy and momentum exchanges with the solid, whereas IAS from insulator surfaces was well described by simple two-body interatomic potentials 21 . However, a possible association of the additional branch to singleparticle electron-hole (e-h) excitations was, at about the same time, proved theoretically to be very unlikely due to the very weak cross-section, well below the sensitivity of conventional IAS spectroscopy [21][22][23][24] . Only in the spectral region where e-h excitations are strongly hybridized with surface phonons the corresponding sharp non-adiabatic anomalies could be observed with He-atom scattering (HAS) 25,26 and explained in terms of single e-h surface excitations 27 . In the latter paper it was however anticipated that also low-energy surface collective electron excitations like two-dimensional (2D) plasmons and ASP would be within present IAS sensitivity.
The experimental proof of the existence of SASP mentioned in the introduction suggested that the unexplained branch observed with NeAS could be ascribed to these excitations. To prove this assignment, new NeAS and HAS studies of the Ni(111) surface together with calculations of the self-consitent surface-electron response function and of the first-principles surface phonon dispersion curves have been undertaken in the present work. This allows for a clear assignment of the inelastic features observed with both atom probes within the continuum of the surface-projected phonon to SASP. A density-functional perturbation theory (DFPT) calculation of the Ni(111) surface phonon dispersion curves, besides reproducing the experimental surface phonon branches, permits to assess the effects of SASP-phonon hybridization.
A series of inelastic NeAS TOF spectra has been measured with an incident energy of 66 meV along the ŴM direction of the Ni(111) surface at a temperature of 400 K. The angle between the incident and scattered beams-the source-to-detector angle in planar geometry-is fixed at a total angle θ SD = θ i + θ f = 105.4 • , while the incident angle varies by steps of 1 • from 52.7 o + 9 • (incident parallel wavevector K i = 22.23 Å −1 ) to 52.7 • − 8 ( K i = 17.76 Å −1 ). Hereafter, upper (lower)-case symbols shall be used for 2D (3D) vectors.
The resulting TOF spectra are plotted as functions of the energy transfer in Fig. 1a for the corresponding values of K i . The parallel wavevector transfer K corresponding to the maxima (broken red curve) of the inelastic peaks, are indicated aside the r.h. ordinate scale, and show that the corresponding excitations occur at comparatively small wavevectors. The attribution of these peaks to multi-phonon processes meets the following difficulty: in the limit of a Debye-Waller exponent 2W >> 1 , as is the case in these experiments with a large incident energy and a large surface temperature, the full width at half maximum (FWHM) of the multi-phonon peak can be approximated by (8 ln 2ǫ max ǫ 0 ) 1/2 , where ǫ max is the peak maximum energy and ǫ 0 is an average phonon energy 28 , and this, with, e.g., ǫ max ∼ ǫ 0 ∼ 20 meV, gives FWHM = 47 meV, which is clearly much larger than observed. Indeed, as appears from the following analysis, the observed sequence of peaks can be safely attributed to SASP, with possible avoided crossings with surface optical phonon modes.
These new NeAS data from Ni(111) are here integrated by an analysis of the extended set of NeAS data (Fig. 1b), privately communicated by Feuerbacher to V. Bortolani and reproduced in their theoretical study 29 . This set actually shows a few additional data points at wavevector transfers close to G-vector multiples and energies up to 4 times the maximum phonon energy ǫ M ≈ 38 meV (Fig. 1b). The latter fact should exclude the possible association of the additional peaks to multiphonon processes 30 . The multiphonon background is clearly present beneath the one-phonon spectral region but does not seem to extend beyond ±2ǫ M (cfr. Fig. 1b,  www.nature.com/scientificreports/ On the other hand, the observation of inelastic Umklapp processes with various G vectors is somewhat surprising in NeAS, since diffraction studies from metal surfaces normally do not reveal any appreciable surface corrugation. Actually the excitation of phonons by atom collisions on metal surfaces occurs through the creation of a virtual e-h pair, eventually recombining into a phonon 28 . The same holds for the creation of a SASP, and the matrix element converting an e-h pair into a SASP is governed by a potential reflecting lattice periodicity, and therefore allowing for Umklapp processes. This supports the assignment of the additional and long mysterious steep dispersion curve observed in Ni(111) to SASP. Also subsequent HAS data from a 10 ML Ni(111) film grown on diamond 31 showed additional modes (see, e.g., the TOF spectrum in Fig. 2a) forming a very steep dispersion curve, which was tentatively ascribed to the diamond substrate. However, the steepness of that additional branch, much larger than that of the diamond RW 32 provides a further argument in favour of Ni(111) SASP. It is interesting to note that in similar NeAS measurements on Ir(111) 33 , where the surface electronic structure at the Fermi level does not support any ASP, the observed peak dispersion and FWHM can actually be fairly well assigned to multiphonon processes.

Surface phonons and SASP of Ni(111) from He atom scattering
A comparison between two of the present TOF spectra (Fig. 2b) and a HAS spectrum for the 10ML-Ni(111) film measured by Braun et al. 31 , reproduced in Fig. 2a, (all converted to the energy-transfer scale) shows, aside the diffuse elastic peak (E) and the dominant RW peak, smaller features at wavevector Q ∼ 0, ascribed to SASPs. Note that the SASP localization within the 10 ML quantum well can justify the comparably larger SASP intensity in the film than in the semi-infinite crystal. Panels (c, d) of Fig. 2 show other examples of SASP features seen at higher energies, but still within the phonon spectrum.
The full sets of present and past HAS (full and open circles) and NeAS data (five-stars: present; crosses: Refs. 20,34 ) measured along [112 ] are plotted in Fig. 3 after folding into the positive quadrant of the first SBZ, together with the available HREELS data (lozenges) by Menezes et al. 35 . In Fig. 3a the experimental data are superimposed for comparison to ab initio calculation based on DFPT 36 . The surface dynamics of Ni(111) has been extensively investigated theoretically on the basis of force-constant models 29,35,[37][38][39][40][41][42][43][44][45][46] of effective-medium theory 47 , of a first-principle Green's function method 48 , and of molecular dynamics with the Finnis-Sinclair www.nature.com/scientificreports/ potential 49 . To our knowledge, the present calculation displayed in Fig. 3 is the first one performed with DFPT. The good agreement with the calculation permits a clear assignment of the observed phonons to the RW, the anomalous resonance S3, the sagittal Lucas resonance S4 and the optical Wallis S2 branches 50 . In particular, the DFPT calculation reproduces very well the EELS data 35 at the zone-boundary for the RW (S1), the S3 resonance and the optical gap mode S2.
SASP-phonon coupling. The plot in Fig. 3a, b of the full sets of HAS and NeAS data, whether assigned to phonons or SASP, allows to appreciate the effects of the SASP-phonon coupling. The broken-line eye-guides interpolating the data points suggest, in Fig. 3a,possible avoided crossings between the steep linear dispersion  31 ) and (b-d) for a clean Ni(111) surface (present work). The SASP peak is more intense for the 10 ML film due to quantum-well confinement. Besides the Rayleigh wave (RW), narrow features corresponding to the SASP are clearly visible in the spectra recorded from Ni(111). The rather broad SASP feature seen at 21 meV reflects the avoided crossing with the S4 phonon (see Fig. 3a). Other smaller features, falling where the optical phonon modes S2 and S3 are expected, are visible above the background noise. www.nature.com/scientificreports/ of SASP and the phonon branches, and in Fig. 3b on a compressed energy scale, the continuation of the SASP branch above the phonon spectrum (shown as a grey area). While at high energy the slope of the SASP branch is close to the calculated Fermi velocity of the surface majority-spin electronic band (broken line in Fig. 3b), the SASP branch running inside the phonon spectrum (Fig. 3a) exhibits a considerable slower velocity due to the strong e-p coupling which yields a larger carrier effective mass. In other words, the collective carrier oscillations drag the atoms due to the coincidence of the SASP energy and momentum with those of some bulk phonon of the spectrum. The same coupling leads to the avoided crossing with the surface-localized phonons, and a gradual change of an SASP into a phonon and viceversa. These mixed modes have been termed plasmarons in the original works on semiconductor quantum wells 8 . A similar hybridization is expected to occur for SASP dispersion curves and the surface optical phonon branches of other metal surfaces (notably with the shear-vertical (S2) and longitudinal Lucas (S4) optical modes).

Theory of atom scattering from spin-acoustic surface plasmons
Distorted-wave Born approximation. The way an incident closed-shell atoms can exchange energy and momentum with a metal surface can be understood within the Esbjerg-Nørskov (EN) approximation for the atom-surface potential 16 . In this approximation only the steep repulsive part V (r, t) of the atom-surface potential energy is considered and written in the form where n(r, t) is the surface electron density and A n = 364 eV/a.u. is the Esbjerg-Nørskov constant as currently used for He atoms on metal surfaces 51 . A different value of A n is expected for Ne, due to its larger attractive dispersion forces and larger size with respect to He. If one adopts the simple EN potential, these contributions are www.nature.com/scientificreports/ accounted for by an effective EN constant, which is therefore dependent on the probe atom. There is however a large discrepancy between the theoretical value of A n obtained for Ne by Puska et al. (which IAS quite larger than that of He) 52 , and the value, much smaller than that of He, derived by Baumberger et al. from experiments 53 . This discrepancy may be attributed to the difficulty of earlier density functional methods to account for dispersion forces. It is noted, however, that it concerns the constant prefactor in Eq. (3) and is therefore irrelevant for the present discussion. The time dependence of n(r, t) due to the surface elementary excitations determines the energy exchange, as well as the dynamical space modulation determines the exchange of parallel momentum with the surface. The differential reflection coefficient in the distorted-wave Born approximation (DWBA), for an IAS casting the atom of incident energy E i and normal wavevector k iz into the final solid angle f with final energy E f and wavevector k f , is proportional to the Fourier transform of the potential-potential autocorrelation function 28 : where V fi (t) is the time-dependent matrix element of the potential between the final and the incident atom wavefunctions as distorted by the static surface potential. Let ω = E i − E f and Q = K i − K f − G be the energy and parallel wavevector of the surface elementary excitation created in the IAS process, where K i and K f are the parallel components of the incident and final atom wavevectors, respectively, and G is a surface reciprocal lattice vector. In the EN approximation, Eq. (1), one actually needs the autocorrelation function of the charge density matrix elements n fi (t) . For a static flat surface, as generally found for close-packed metal surfaces such as Ni (111) and for small incident energies E i , the static surface charge density only depends on the coordinate z normal to the surface: n 0 (r) = n 0 (z t )e −κz , where the He-surface classical turning point z t = z t (E i ) is taken as the origin of z, and κ is the repulsive parameter. The matrix element A n n 0fi (r) between final and incident distorted waves of the corresponding static exponential potential can be calculated analytically in the DWBA 54 , which allows to express the IAS differential reflection coefficient as where is the Fourier transform of the dynamic density-density autocorrelation function, providing the response function of the surface electron density, and is the Jackson-Mott kinematic factor. For a hard-wall potential ( κ → ∞ ) S fi = 1.
As far as it concerns the use of the above equations for NeAS, it is noted that Jackson-Mott theory provides an analytical form of the scattering matrix element for a negligible surface corrugation and an exponential form like the EN potential. These approximations allow for the factorization of Eq. (3), thus separating what is relevant in the present analysis, the density autocorrelation function. Since also Ne diffraction peaks for fcc(111) metal surfaces are found to be much smaller than the specular peak, the assumption of a flat surface is reasonable also for NeAS, so that the EN exponential potential has been adopted in the past also for NeAS with the appropriate constant 53 . Although the J-M kinematic factor does not account for some of the features mentioned above of real surface potentials, like dispersion forces and corrugation effects, the approximations made in deriving Eq. (3) should provide a sufficient basis for a first theoretical analysis of the SASP dispersion relations. As regards the kinematic factor, it is noted, however, that the hard-wall potential model (where S fi = 1 ) is inappropriate for metal surfaces, especially for Ne, and therefore the prefactor S 2 fi needs to be included due to its dependence on the initial and final momenta of the probe atom.
Ni bulk and surface response functions. In this work, we explore two routes to determine a possible collective excitation in the electronic system responsible for the experimental observations. One, based on the first-principle calculations, consists in looking for the possibility of such low-energy excitations in the electronic system of bulk nickel. The second step consists in the evaluation of the surface response function of Ni(111) associated to the IAS. This calculation shows that the shallow majority-spin Shockley surface state occurring around the SBZ center supports the occurrence of a SASP branch.
The calculated normalized bulk loss function N(q, ω) ≡ L(q, ω)/ω at the momentum transfers q along three symmetry directions are displayed in Fig. 4. One can find several peaks with an acoustic dispersion in N(q, ω) in this momentum-energy region. Their existence can be explained by the presence of several energy bands crossing the Fermi level with different velocities. Detailed analysis showed that the prominent peaks, highlighted by dotted lines, can be considered as a signature of a bulk acoustic plasmon. Eventually, upon approaching the q = 0 limit its energy vanishes. However, its group velocities in all three directions are relatively large, substantially exceeding that of the mode interacting with phonons (Fig. 3) observed in the experiment. www.nature.com/scientificreports/ In order to address possible electronic collective excitations at the surface of the solid probed by atom scattering one should evaluate the response function g(Q, ω) , defined in Eq. (4), and here re-written as the Fourier transform of the electronic susceptibility: Its imaginary part, Im[g(Q, ω)] , determines the differential cross section for the scattering with energy ω and in-plane momentum transfer Q and coincides with the surface loss function 55,56 at Q = 0.
In the surface loss function, Im[g(Q, ω)] , reported in Fig. 5b one can observe a peak labeled as SASP which we ascribe to the spin acoustic surface plasmon. We attribute its origin to the partly occupied majority-spin surface  www.nature.com/scientificreports/ state, Fig. 5a. The Fermi velocity in this surface state band is significantly lower in comparison with that in the sp bulk energy bands. The resulting incomplete dynamical screening, involving slow carriers in the surface state band and the faster ones in the bulk electron system, is the origin of this mode 57 .
In Fig. 5b one can see that the variation in temperature T introduces notable modifications in the calculated surface loss function. First, with increasing T the SASP peak shifts to larger energies. Second, the broadening of the SASP peak increases as well. Nevertheless, even at room temperature, the SASP peak can be resolved in the calculated surface loss function.
In the upper panel of Fig. 5c we show the calculated imaginary part of the induced density, Im[n ind ], corresponding to the SASP at Q = 0.001 a.u. and ω=2.6 meV. One can see, that the peak positions in Im[n ind ] correspond to the charge density distribution observed in the majority-spin Shockley surface state reported in the bottom panel of Fig. 5c. As such, the positive parts in Im[n ind ] of Fig. 5c have mostly majority-spin orientation. This is schematically shown by the upward arrow. On the contrary, the negative parts of Im[n ind ] corresponding to screening by faster bulk electrons contain contributions from both the majority-and minority-spin electrons. Therefore, in a first approximation, this part in Im[n ind ] can be considered as paramagnetic. Since, the total charge in Im[n ind ] is almost zero due to the quadrupole nature of the SASP, the resulting spin of the oscillating charge related to the SASP should have excess of the majority-spin type. Moreover, since in the vacuum side the majorityspin surface state dominates, the oscillating charge related to the SASP should be mostly majority-spin resolved.
In Fig. 6 we report the two-dimensional plot of the calculated surface loss function as a function of Q and ω . The area in Im[g(Q, ω)] denoted as 2D and corresponding mainly to intra-band excitations inside the majorityspin surface state clearly has larger amplitude in comparison to the regions where only excitations inside the bulk system are allowed. The Ne scattering data substantially fall into the 2D area of the electronic surface loss function. Due to the comparatively large experimental uncertainty it cannot be decided how far is the predicted SASP dispersion curve reproduced by experiment. On the other hand, the early conclusion by Schönhammer and Gunnarsson 21,22 about the very small IAS cross-section of single e-h excitations, supports the assignment of the observed inelastic processes at small wavevector transfers and energy rapidly increasing well beyond the phonon maximum energy to electron collective excitations, notably to SASP. Thus, the change of slope from inside to outside the phonon spectrum, as well as the apparent gap between 30 and 35 meV associated to the avoided crossing between the SASP and the S2 phonon branch, can be attributed to e-p interaction. This is compatible with the comparatively large value of the mass-enhancement factor in nickel ( = 0.56) 58 derived from the temperature dependence of the HAS Debye-Waller exponent 59 . www.nature.com/scientificreports/

Conclusions
The analysis of inelastic atom scattering data from the (111) surface of nickel, based on calculations of the electronic response function and the surface phonon dispersion curves proves IAS as a valuable tool to investigate ASPs in the THz region, that is hardly accessible to other spectroscopies. The experimental findings were explained by a self-consistent calculation of the bulk and surface response functions and of the first-principle surface-phonon dynamics of Ni(111). We have shown that in Ni(111), the ASP originates from the majority-spin Shockley surface state and is therefore a collective oscillation of surface electrons with the same spin, i.e. it represents a new kind of collective quasiparticle: a spin acoustic surface plasmon (SASP). Present monochromatic atom scattering spectrometers, based on supersonic beams with a very high speed ratio, definitely attain a far better resolution than available in the early 80s, and are therefore suitable to investigate surface plasmon excitations in the THz frequency domain. 3 He spin-echo (3HeSE) spectrometry, in particular, with its unprecedented resolution in the neV range, can in principle measure the ASP dispersion in the GHz domain. Recent HAS experiments have provided the dispersion curves of other THz electronic collective excitations such as phasons and amplitons occurring in charge density wave systems 60,61 and, with 3 HeSE, in the long-period reconstructed Au(111)23 × √ 3 surface 62 . As shown in the case of Bi 2 Se 3 (111), bound-state resonance enhancement can be exploited in HAS experiments to amplify the response of low-energy electronic collective excitations 61 . On the other hand, fast neutral atom beams in grazing-incidence scattering experiments in the keV energy domain have also proved to have sufficient resolution to detect surface excitons in insulators 63 and plasmon coherent excitations at metal surfaces 64 : nothing prevents in principle their use at smaller incident energy and higher monochromaticity for ASP at metal surfaces.

Methods
Experimental details. The present HAS measurements have been performed with a high-resolution HAS TOF spectrometer described in detail elsewhere 65 . A pulsed atom beam is produced by free jet expansion of the high pressure gas and a rotating disk chopper. The atoms scattered from the sample travel through three differential pumping stages along the 1.7-m-long drift tube before reaching the detector. The angle between the incident and scattered beams (the source-to-detector angle in planar geometry) is fixed at a total angle θ SD = θ i + θ f = 105.4 • . The low intensity of the diffuse elastic peak indicates the high quality and cleanness of the crystal surface. The TOF spectra were integrated for typically 60-70 min at a surface temperature T = 460 K to avoid surface contamination with H 2 . We discussed in more detail the TOF spectra and phonon modes of the Ni(111) single crystal in an earlier work 66 .
DFPT calculations. The DFPT calculations were carried out using the Quantum-ESPRESSO package 67 . An ultrasoft pseudopotential 68 with ten valence electrons and the Perdew-Burke-Ernzerhof (PBE) approximation 69 for the exchange-correlation energy functional were used. The Kohn-Sham orbitals were expanded in plane waves up to an energy cutoff of 45 Ry while an 800 Ry cutoff was used for the charge density. A Gaussian smearing of 0.02 Ry was introduced in the occupation of states to deal with the metallic character of Ni while collinear spin-polarization was used to describe its ferromagnetic ground state. We optimized the bulk geometry by integrating the Brillouin zone (BZ) over a 12 × 12 × 12 Monkhorst-Pack (MP) mesh 70 obtaining an equilibrium lattice parameters of 3.518 Å. The surface phonon dispersion curves were first calculated for a 12-layer slab separated by its periodic replica by 25 Å of vacuum. The hexagonal SBZ was sampled over a 12 × 12 × 1 MP grid. Atomic positions were relaxed until forces became lower than 0.1 mRy/bohr and phonons were computed on a 4 × 4 × 1 grid. In order to obtain denser bands, thicker 38 layers slab was constructed through the insertion of a 26 layers bulk slab into the middle of the original 12-layer slab. The surface 6-layer slabs are sufficient to make the effect of the inserted bulk layers on the surface phonon negligible.
Bulk and surface response functions calculations. A way to obtain information on the excitations in a bulk electron system is established in the framework of time-dependent density functional theory 71,72 . The well-defined peaks in the loss function L(q, ω) ≡ −Im[ǫ −1 gg (q ′ , ω)] , directly probed in the energy-loss experiments, correspond to the plasmonic excitations. Here q ≡ g + q ′ with g being the bulk reciprocal lattice vectors. Vector q ′ is chosen to be inside the BZ. The inverse dielectric-function matrix ǫ −1 gg (q ′ , ω) is related to the densityresponse function of interacting electrons χ as where V g ′ (q ′ ) = 4π/|q ′ + g ′ | 2 is the Fourier transform of the bare Coulomb potential. χ gg ′ (q ′ , ω) is obtained from the matrix equation The kernel K xc accounts for the exchange-correlation effects, which are set to zero within the random-phase approximation (RPA) employed here. In general, impact of these effects beyond the RPA in the dielectric properties of solids is small at small momentum transfers q which are of interest in this work. In the above equation, χ o (r, r ′ , ω) is the response function of the non-interacting Kohn-Sham electrons, which in the reciprocal space reads www.nature.com/scientificreports/ Here is the normalization volume, f nsk is the Fermi occupation number, and summation over spin s is explicitly included. The sampling over k in the BZ was performed on a 300 × 300 × 300 mesh. The one-particle energies ε nsk and wave functions ψ nsk of Ni were obtained in the electron density functional theory formalism implemented in the Vienna Ab-initio Simulation Package (VASP) [73][74][75] . Since we are interested in the plasmons determined by the intra-band transitions 5,57 , in Eq. (9) we include only such transitions. The corresponding matrix elements were set to unity. The occupied (occ) and unoccupied (unocc) states in the energy bands crossing the Fermi level were taken into account. In the Fourier expansion of χ o , χ , and ǫ −1 we employed the g = 0 vector only, i.e. the local field effects were not considered, which in the low-q and low-ω domain is a good approximation. Other calculation details can be found, e.g. in Ref. 76 . Nowadays to perform a full ab initio calculation of g(Q, ω) 77 for a Ni surface with a meV energy resolution using a sufficiently thick slab represents a formidable task. Therefore we employed significantly less timedemanding approach based on a model potential 78 capturing essential features of the surface electronic structure of Ni(111). We constructed a spin-resolved potential for the Ni(111) surface following the recipe of Ref. 78 . In particular, the bulk energy gap positions in both spin channels and the sp Shockley surface state in the majorityspin case were taken into account. The resulting model band structure of the Ni(111) surface around the Fermi level obtained with this potential is presented in Fig. 5a. The data for the fitting of the model potential parameters were taken from the density-functional theory (DFT) calculations 15,[79][80][81] . Since the presence of the minority-spin sp surface state in the Fermi level vicinity is not confirmed by the DFT calculation, it was not taken into account in the model potential construction.
Notice that there is a large spread in the energy positions of the Shockley surface states in different experiments. For instance, the absolute binding energies of the majority-spin surface state at the band bottom range from 225 meV in the scanning tunneling microscopy measurements 14 to almost zero in the photoemission experiment 82 . It was pointed out in Ref. 83 that the energy position of this surface state is strongly temperature dependent. In turn, the Fermi velocity v F of the surface state changes as well. In consequence, the SASP dispersion, that is proportional 57 to v F , may be sensitive to the factors determining the surface state energy position. Our fitted Ni(111) electronic structure is close to that measured in the spin-resolved photoemission experiment carried out at room temperature 83 .
With this model potential we calculated the electronic structure of a slab consisting of 51 atomic layers. The resulting one-particle energies and wave functions where employed in the evaluation of the g function along the lines described in Ref. 84 . An important novel issue of the present calculations is the inclusion of the electron temperature in the evaluation of the occupation factors f nsk entering Eq. (9). Since the temperature of the experiment is comparable with the binding energy of the surface state at the bottom, it may have an impact in the surface excitation spectra.