Arxiv:Cond-Mat/0408683V1 [Cond-Mat.Mtrl-Sci] 31 Aug 2004 Onaisi Ocueteapaec Fanw”Lowered” New a of ¯ at Appearence Loss the Cause to Is Boundaries Mass

arXiv:cond-mat/0408683v1 [cond-mat.mtrl-sci] 31 Aug 2004 onaisi ocueteapaec fanw”lowered” new a of ¯ at appearence loss the cause to is boundaries mass. huddces togywt eraigfi thickness. foil decreasing loss with the Ritchie plasma find strongly reached for not decrease and probability did should the plasma he that bounded result, is conclusion a the field erroneus As in electric modes the surface. that surface the assumed at but zero foils, always thin in tions h aec lcrndniyand density electron valence the fteepam siltos h ai nto nrybeing energy plasmon: of excitation the unit the basic termed the to elec- oscillations, due plasma fast are these by metals of through experienced passing losses in energy trons plasma discrete of the Bohm form that and the Pines in itself oscillations. collec- manifesting yield to behaviour, known is tive metals in electrons valence tween h ufepamno ufc nriswsaddressed, was energies surface on point plasmon zero field, surfce shifted plasmon the the surface sur- in the metal originated of a be energy potential to and shown charge image was ex- classical classical face of point variety the a large between example, acting a For of of interpretation properties the periments. fundamental in of and understanding solids the in role the in and modes theoretical collective both interface. vacuum-solid in of advance there investigations then, significant experimental Since a Mg. been from and Al has scattered of electrons layers evaporated by newly experienced losses inelastic mnscridotb oeladSwan and Powell exper- by of series out a in carried confirmed iments was metals in losses energy etv siltos h unao hc tr n Ferrell plasmon. and surface Stern the which called of quanta the oscillations, lective Gabor h ogrnentr fteCuobitrcinbe- interaction Coulomb the of nature long-range The h ocp ftesraepamnhspae key a played has plasmon surface the of concept The ici’ rdcino ufc oaiaincuiglow- causing polarization surface of prediction Ritchie’s 3 hω 5 4 a h rtt n htteeeto h film the of effect the that find to first the was s netgtdteectto fpam oscilla- plasma of excitation the investigated ¯ = ASnmes 14.m 32.t 84.p 78.68.+m 78.47.+p, 73.20.At, 71.45.Gm, numbers: PACS low-en exist. new to a predicted introduce recently and been con structures, has basic the metallic review fun in we of Here plasmons areas applications. of photonic variety surface- of microscopy, range a surface-plasmon in to role dynamics significant surface a played have mons hω ic h ocp fasraecletv xiainwsfirs was excitation collective surface a of concept the Since 1 .INTRODUCTION I. p aei odnstae iiaSia inz Fakultatea, Zientzi Saila, Fisika Kondentsatuaren Materia / 2 √ ootaItrainlPyisCne DP)adCnr Mi Centro and (DIPC) Center Physics International Donostia .M Pitarke, M. J. u oteectto fsraecol- surface of excitation the to due 2 2 auld adzblPsaeu,E208Dnsi,Basque Donostia, E-20018 Pasealekua, Lardizabal de Manuel ¯ hω 3 aeilnFsk al,Kmk autta uklHerrik Euskal Fakultatea, Kimika Saila, Fisika Materialen p ¯ = 6 02Psaktaia -08 oota aqeCuty Sp Country, Basque Donostia, E-20080 kutxatila, Posta 1072 h 1 4 ot uxtl,E400Blo aqeCuty Spain Country, Basque Bilbo, E-48080 kutxatila, Posta 644 (4 eetefis osuggest to first the were ufc lsosi ealcstructures metallic in plasmons Surface πne m 1 , 8,9,10,11 e 2 2 .M Silkin, M. V. , /m stefe electron free the is 7 e h observed who , ) 1 / h matof impact the 2 where , Dtd etme ,2018) 5, September (Dated: 2 n .V Chulkov, V. E. 12 is ea ufc a hw ob u oteexcitation the to due a be outside plasmons, to moving surface shown particles of was charged surface of metal loss energy the ea ufc a eosrtdt edcae ythe by dictated plasmons. surface be a of to on dependence demonstrated acting wave-vector was potentials surface external by metal induced density tron utpl ufc lsosta a enpeitdby so-called predicted been the had of that Benett. existence plasmons the surface for multipole evidence showed also nlyo ipemdl,sc stehydrodynamic, the as such models, specular-reflection, simple on orig- relied inally dispersion surface-plasmon the for expressions esl,tesraepamncno ea yeitn a emitting by decay cannot photon. con- plasmon and, surface plasmons surface the on versely, excite incident cannot light surface that ideal so frequency, an same surface the and is of light plasmons there between curve mismatch Hence, momentum dispersion a radiation. always the electromagnetic crosses free-space never of curve dispersion ing in sarsl frtrain ufc lsoscouple plasmons than surface time larger retardation, interac- wavelengths the of retarded at result account true a the into As of take tion. solid propagation to inter- the with for necessary the light needed is or describe electrons it to in- relativistic surfaces, is is either one interaction of if Coulomb action However, the that stantaneous. approxima- assuming nonrelativistic by the in tion, Ritchie by derived was zdsrae fteakl metals. alkali the character- well of clean, surfaces on ized out carried been have that ments density-functional eigenval- of (DFT), hamiltonian and theory Kohn-Sham eigenfunctions the the of of performed ues knowledge been the also have from calculations numerical curate ag rmzr oa smttcvledtrie ythe by determined frequency value asymptotic surface-plasmon an to zero from frequency range entire the over surface- exists a polariton surface-plasmon called now polariton. is what plasmon yield and field tromagnetic h long-wavelength The lso eoac ehooy n wide a and technology, resonance plasmon et neligteeitneo surface of existence the underlying cepts 22 2 rysraecletv xiainthat excitation collective surface ergy nrdcdb ici,sraeplas- surface Ritchie, by introduced t , aetladapidrsac,from research, applied and damental 3 n .M Echenique M. P. and , uklHrioUnibertsitatea, Herriko Euskal 19 hwn ieareetwt h experi- the with agreement nice showing t CSIC-UPV/EHU, xto Unibertsitatea, o 17 24 onr,Spain Country, 13,14 n infinite-barrier and tteelrewvlnts the wavelengths, large these At 23 ain n h etodo h elec- the of centroid the and ω ufc-lso nry¯ energy surface-plasmon ∼ s oee,tecorrespond- the However, . 2 πc/ω 20,21 2 , 3 s hs experiments These ihtefe elec- free the with 18 oes Ac- models. 15 Explicit hω 16 s 2

Teng and Stern25 were the first to point out that any surface-plasmon excitation plays in the interaction of fast surface roughness permits the surface to impart some charged particles with solid surfaces, since it is precisely additional momentum to the surface-plasma oscillation, the investigation of electron energy loss in thin foils which so that it can couple to electromagnetic radiation. Al- brought Ritchie to the realization that surface collective 5 ternatively, prism coupling can be used to enhance the excitations exist at the lowered frequency ωs. The detec- momentum of incident light, as demonstrated by Otto26 tion of surface plasmons and their dispersion is discussed and by Kretchmann and Raether.27 Hence, surface plas- in Section III, in the framework of angle-resolved inelastic mons have been employed in a wide spectrum of stud- electron scattering experiments. Localized surface plas- ies from electrochemistry, catalysis, wetting, thin or- mons and the use of sum rules that provide insight of ganic condensates, and biosensing,28 to scanning tunnel- surface-plasmon energies in metallic structures of arbi- ing microscopy,29 the ejection of ions from surfaces,30 trary geometry are introduced in Section IV. Section V surface dynamics,31 surface-plasmon microscopy,32 and is devoted to the retarded region, where surface plasmons surface-plasmon resonance technology.33 Moreover, re- couple with the free electromagnetic field. Acoustic sur- cent advances that allow metals to be structured and face plasmons are introduced in Section VI. characterized on the nanometer scale have rekindled Unless stated otherwise, atomic units are used 2 the long-standing interest in surface plasmons, one of throughout, i.e., e =¯h = me = 1. the most attractive aspects of these collective excitations now being their use to concentrate light in sub- wavelength structures34 and to enhance transmission II. PLASMA LOSSES BY FAST CHARGED through periodic arrays of subwavelengths holes in op- PARTICLES IN SOLIDS tically thick metallic films,35,36 as well as the possible fabrication of nanoscale photonic circuits operating at Let us consider a recoiless fast point particle of charge 37 optical frequencies. Z1 moving in an arbitrary inhomogeneous many-electron Since the typical energy of bulk and surface plasmons system with nonrelativistic velocity v, for which retar- is of a few electronvolts, thermal excitation of these dation effects and radiation losses can be neglected.47 collective oscillations is improbable, so that the elec- The charge density of the probe particle is simply a delta tronic properties near the Fermi level cannot be directly function of the form influenced by these excitations. Much more effective ρext(r,t)= Z δ(r r v t), (1) than ordinary bulk and surface plasmons in mediating, 1 − 0 − e.g., superconductivity would be the so-called acoustic 38 and the energy that this classical particle loses per unit plasmons with sound-like long-wavelength dispersion, time due to electronic excitations in the medium can be which have spurred over the years a remarkable inter- written as48 est and research activity.39 Acoustic plasma oscillations were observed in two-dimensionally confined and spa- dE ∂V ind(r,t) = dr ρext(r,t) , (2) tially separated multicomponent structures such as quan- − dt − ∂t tum wells and heterojunctions,40,41 and were then pro- Z ind posed as possible candidates to mediate the attractive where V (r,t) is the potential induced by the probe interaction leading to the formation of Cooper pairs in particle at position r and time t. 42,43 high-Tc superconductors. To first order in the external perturbation, time- Recently, it has been shown that metal surfaces where dependent perturbation theory yields a partially occupied quasi two-dimensional (2D) surface- +∞ +∞ dω ′ state band coexists in the same region of space with the V ind(r,t) = dr′ dt′ e−iω(t−t ) 2π underlying three-dimensional (3D) continuum support a Z Z−∞ Z−∞ well-defined acoustic surface plasmon.44 This new low- energy collective excitation exhibits linear dispersion at W˜ (r, r′; ω) ρext(r′,t′), (3) low wave vectors, and might therefore affect electron-hole × (e-h) and phonon dynamics near the Fermi level.45 It has where been demonstrated that it is a combination of the nonlo- W˜ (r, r′; ω)= W (r, r′; ω) v(r, r′), (4) cality of the 3D dynamical screening and the spill out of − the 3D electron density into the vacuum which allows the v(r, r′) being the bare Coulomb interaction and formation of 2D electron-density acoustic oscillations at W (r, r′; ω) being the so-called screened interaction, metal surfaces, since these oscillations would otherwise be which is usually expresed in terms of the density-response 46 completely screened by the surrounding 3D substrate. function of the many-electron system.49 In this paper, we first present an overview with the ba- Eq. (3) is a general expression for the energy loss of a sic concepts underlying the existence of surface collective classical particle moving in an arbitrary inhomogeneous excitations in metallic structures, and we then introduce electron system that is characterized by the screened in- the new concept of acoustic surface plasmons. We be- teraction W (r, r′; ω). Here we consider a solid target con- gin in Section II with a brief discussion of the role that sisting of a fixed uniform positive background (jellium) 3 plus a neutralizing cloud of interacting valence electrons, assuming that v >>vF , Eq. (5) yields which will be described by either an infinite or a plane- bounded electron gas. dE ω2 k 2v = Z2 p ln c + ln , (9) − dx 1 v2 ω /v k p c A. Infinite electron gas or, equivalently, In the case of an infinite homogeneous electron gas that is translationally invariant in all directions, Eqs. (1)-(3) ω2 2 yield the following expression for the so-called stopping dE 2 p 2v = Z1 2 ln . (10) power, i.e., the energy that the probe particle loses per − dx v ωp unit path length: The first term of Eq. (9), which can also be obtained by using the Drude dielectric function of Eq. (8), represents dE d3k the contribution to the stopping power from losses to col- = Z2 k v ImW˜ (k, k v), (5) − dx − 1 (2π)3 · · lective excitations at wave vectors of magnitude smaller Z than kc. The second term of Eq. (9) represents the contri- W˜ (k,ω) being the 3D Fourier transform of W˜ (r, r′; ω), bution from losses to single-particle excitations at wave which is typically expressed in the form vectors above kc. 4π W˜ (k,ω)= ǫ−1(k,ω) 1 , (6) k2 − where ǫ−1(k,ω) is the so-called inverse dielectric function B. Plane-bounded electron gas of the electron gas. At high projectile velocities (v << v, v being the F F In the case of a plane-bounded electron gas that is Fermi velocity), the zero-point motion of the electron gas translationally invariant in two directions, which we take can be neglected and it can be considered, therefore, as to be normal to the z axis, Eqs. (1)-(3) yield the following if it were at rest. In this approximation, the dielectric expression for the energy that the probe particle loses per function ǫ(k,ω) of a homogeneous electron gas takes the unit time: form50

2 2 2 +∞ ∞ ωp dE Z d q ǫ(k,ω)=1 , (7) = i 1 dt′ dω ω − ω(ω + iη) k4/4 − dt π (2π)2 − Z Z−∞ Z0 which at long wavelengths (k 0) yields the classical ′ → e−i(ω−q·vk)(t−t ) W˜ [z(t),z(t′); q,ω], (11) Drude dielectric function × ω2 ǫ(ω)=1 p , (8) where q is a 2D wave vector in the plane of the surface, − ω(ω + iη) vk represents the component of the velocity that is parallel to the surface, z(t) represents the position of the ωp being the bulk plasma frequency and η a positive in- ˜ ′ finitesimal. projectile relative to the surface, and W (z,z ; q,ω) is the ˜ ′ The dielectric function of Eq. (7) describes both col- 2D Fourier transform of W (r, r ; ω). lective and single-particle excitations. At wave vec- In the simplest possible model of a bounded semi- tors k smaller than a cut-off wave vector of magnitude infinite electron gas, one characterizes the electron gas kc ωp/vF , the many-electron system can be expected at z 0 by a local dielectric function which jumps dis- to behave∼ collectively, losses being dominated by the ex- continuosly≥ at the surface from unity outside (z < 0) to citation of plasma oscillations.2 At wave vectors of mag- ǫ(ω) inside (z > 0). By imposing the ordinary boundary nitude larger than kc, losses are dominated by the ex- conditions of continuity of the potential and the normal citation of electron-hole (e-h) pairs. Using Eq. (7) and component of the displacement vector at z = 0, one finds

−q(|z|+|z′|) 2π g(ω)e z < 0 W˜ (z,z′; q,ω)= − (12) q  ′ ′ ǫ−1(ω) 1 e−q|z−z | + ǫ−1(ω) g(ω)e−q(|z|+|z |), z> 0,  −  where is the long-wavelength (q 0) limit of the so-called → ǫ(ω) 1 g(ω)= − (13) ǫ(ω)+1 4 surface-response function.51 In this limit, the dielectric the surface (z < 0) and enters the solid at z = t = 0. The function ǫ(ω) takes the Drude form dictated by Eq. (8). position of the projectile relative to the surface is then In the following, we shall explicitly consider particle z(t)= vt. Assuming that the electron gas at z 0 can be trajectories that are normal and parallel to the surface. described by the Drude dielectric function of Eq.≥ (8) and introducing Eq. (12) into Eq. (11), one ﬁnds the following expression for the energy that the probe particle loses per 1. Normal trajectory unit path length:

Let us consider a situation in which the probe particle moves along a normal trajectory from the vacuum side of

2 2 ωs f(2ωs z /v), z< 0 dE Z1 | | = 2 (14) − dz v  ω2 [ln(k v/ω ) h(ω z/v)] + ω2 h(ω z/v), z> 0,  p c p − p s s  where h(α) = 2cos(α) f(α) f(2α), (15) − with f(α) being given by the following expression: ∞ x e−αx f(α)= . (16) 1+ x2 Z0 The Drude dielectric function of Eq. (8), which as- sumes that infinitely long-lived plasmons at a single- frequency ωp are the only possible bulk excitations, is only sustainable at wave vectors below a cut-off k . c FIG. 1: Particle of charge Z1 passing through a finite foil of Hence, this cut-off has been introduced into the bulk (log- thickness a. arithmic) term of Eq. (14), which yields a contribution to the energy loss that coincides with the plasmon contribution of Eq. (9). Contributions to the energy loss that This is the result first derived by Ritchie in a differ- are due to the excitation of e-h pairs are not included in ent way,5 which brought him to the realization that sur- Eq. (14). face collective excitations exist at the lowered frequency In the absence of the boundary at z = 0, the position- ωs. The first term of Eq. (17), which is proportional to dependent f(α) and h(α) terms entering Eq. (14) would the thickness of the film represents the bulk contribu- not be present, and the energy loss would be that of tion, which would also be present in the absence of the charged particles moving in an infinite plasma. When boundaries. The second and third terms, which are both the probe particle is moving outside the solid, the effect due to the presence of the boundaries and become more of the boundary is to cause energy loss at the lowered important as the foil thickness decreases, represent the plasma frequency ωs. When the probe particle is moving decrease in the energy loss at the plasma frequency ωp inside the solid, the effect of the boundary is to cause and the energy loss at the lowered frequency ωs, respec- both a decrease in loss at the bulk plasma frequency ωp tively. Eq. (17) also shows that the net boundary effect is and an additional loss at the lowered plasma frequency an increase in the total energy loss above the value which 5 ωs, as predicted by Ritchie. would exist in its absence, as noted by Ritchie.5 Now we consider the real situation in which a fast Ritchie also considered the coupling that exists be- charged particle passes through a finite foil of thickness tween the two surfaces for finite values of the film thick- a (see Fig. 1). Assuming that the foil is thick enough ness a. He found the following dispersion relation be- for the effect of each boundary to be the same as in the tween the frequency of surface-plasma oscillations and case of a semi-infinite medium, and integrating along the the wave number q: whole trajectory from minus to plus infinity, one finds 1/2 the total energy that the probe particle loses to collec- ω = ω 1 e−aq , (18) tive excitations: s ± 2 Z kcv π the exponential factor being a consequence of the inter- ∆E = 1 aω2 ln ω + πω . (17) − v2 p ω − 2 p s action between the two surfaces. This equation has two p 5

which were later discussed by Stern53 and observed in ar- tiﬁcially structured semiconductors54 and more recently in a metallic surface-state band on a silicon surface.55

2. Parallel trajectory

Now we restrict our attention to the case of charged particles moving with constant velocity at a ﬁxed dis- FIG. 2: Particle of charge Z1 moving at a ﬁxed distance z from the surface of a plane-bounded electron gas. tance z from the surface (see Fig. 2). Eq. (11) then yields

52 dE 2 d2q ∞ limiting cases, as discussed by Ferrell. At short wave- = Z2 dω ω lengths (qa >> 1), the surface waves become decoupled − dx −v 1 (2π)2 Z Z0 and each surface sustains independent oscillations at the ImW˜ (z,z; q,ω) δ(ω q v). (20) reduced frequency ωs characteristic of a semi-inﬁnite elec- × − · tron gas with a single plane boundary. At long wavelengths (qa << 1), there are ”normal” oscillations at ωp and ”tangential” 2D oscillations at Assuming that the electron gas at z 0 can be described by the Drude dielectric function≥ of Eq. (8) and 1/2 ω2D = (2πnaq) , (19) introducing Eq. (12) into Eq. (20), one ﬁnds

2 2 ωs K0(2ωs z /v), z< 0 dE Z1 | | = 2 (21) − dx v  ω2 [ln(k v/ω ) K (2ω z/v)] + ω2 K (2ω z/v), z> 0,  p c p − 0 p s 0 s 

56 where K0(α) is the zero-order modified Bessel function. lifetime of this mode are determined by the correspond- For particle trajectories outside the solid (z < 0), ing energy-loss peak in the spectra, and the momentum Eq. (21) reproduces the classical expression of Echenique q parallel to the surface is obtained from the measured 13 and Pendry. For particle trajectories inside the solid angles θi and θf . (z > 0), Eq. (21) reproduces the result first obtained by Nuñez et al..57 As in the case of a normal trajectory, when the particle moves inside the solid the effect of the The inelastic-scattering cross section corresponding to boundary is to cause a decrease in loss at the bulk plasma a process in which an electronic excitation of energy ω and parallel momentum q is created at the surface of a frequency ωp and an additional loss at the lowered plasma frequency ω . semi-infinite electron gas is dictated by the imaginary s part of the surface-response function g(q,ω),51 which in the long-wavelength limit is given by Eq. (13). In a free- III. INELASTIC ELECTRON SCATTERING electron gas (jellium), the dielectric function ǫ(ω) entering Eq. (13) is the Drude dielectric function of Eq. (8). Hence, in a gas of free electrons Img(q,ω) becomes a delta The most commonly used experimental arrangement function peaked at the surface-plasmon energy ω . for detecting surface plasmons by the fields of moving s charged particles is based on angle-resolved inelastic electron scattering.58,59 Fig. 3 shows a schematic drawing of The classical picture leading to Eqs. (8) and (13), the scattering geometry. A monochromatic beam of elec- which is correct only in the long-wavelength limit, ig- trons, incident on a flat surface at an angle θi, is back nores both the nonlocality of the electronic response of scattered and detected by an angle-resolved energy an- the system and the microscopic spatial distribution of alyzer positioned at an angle θf . Inelastic events can the electron density near the surface. Feibelman showed occur, either before or after the elastic event, by excit- that up to first order in an expansion in powers of q, ing a surface mode of frequency ω(q). The energy and the surface-response function of a jellium surface can be 6

TABLE I: Relative widths ∆ωs/ωs of surface plasmons, as derived from the imaginary part of the surface-response function of Eq. (13) with measured values of the bulk dielectric function ǫ(ω) (theory)51 and from the surface-loss measurements at q = 0 (experiment).60 Ag K Al Mg Hg Li Theory 0.027 0.035 0.035 0.16 0.18 0.33 Experiment 0.027 0.1 0.24 0.16 0.16 0.35

FIG. 3: Schematic drawing of the scattering geometry in angle-resolved inelastic electron scattering experiments. By exciting a surface mode of frequency ω(q), the energy of detected electrons becomes εf = εi ω(q). The momentum q is − determined from q = √2 √εi sin θi √εf sin θf −

FIG. 5: Complementary systems in which the regions of plasma and vacuum are interchanged. The top panel represents the general situation. The bottom panel represents a FIG. 4: Schematic drawing of the bulk and surface plasmon half-space ﬁlled with metal and interfaced with vacuum. The energy dispersions in typical metal surfaces. While in the case surface-mode frequencies ωs1 and ωs2 of these systems fulﬁll of bulk plasmons the initial slope is usually (but not always) the sum rule of Eq. (23). positive, the surface-plasmon energy dispersion is negative for typical metal surfaces. the dielectric function ǫ(ω), the experimental surface-

15 plasmon widths ∆ωs at q = 0 should be approximately written as described by using in Eq. (13) the measured bulk di- ǫ(ω) 1 ǫ(ω) electric function ǫ(ω) instead of its Drude counterpart. g(q,ω)= − 1+2qd (ω) + O(q2), Table I exhibits the relative widths ∆ω /ω of surface ǫ(ω)+1 ⊥ ǫ(ω)+1 s s plasmons derived in this way,51 together with surface- (22) loss measurements at q = 0.60 Although in the case of where d (ω) represents the centroid of the induced elec- ⊥ Ag, Li, Hg, and Mg the surface-plasmon width is well de- tron density. Self-consistent density-functional calcula- scribed by introducing the measured bulk dielectric functions of g(q,ω) and d (ω) have demonstrated that in the ⊥ tion into Eq. (13), the surface-plasmon widths of K and case of jellium surfaces in the range of typical bulk den- Al are considerably larger than predicted in this sim- sities the surface-plasmon energy ω(q) [where Img(q,ω) ple way. This shows that an understanding of surface- is maximum] is at nonvanishing but small q wave vectors plasmon broadening mechanisms requires a careful anal- lower than ω , in agreement with experiment.20,21 This s ysis of the actual band structure of the solid. Approxi- is illustrated in Fig. 4, where the surface-plasmon energy mate treatments of the impact of the band structure on dispersion is drawn schematically, showing that the ini- the surface-plasmon energy dispersion have been devel- tial slope is negative. For an interpretation of negative oped by several authors,61,62,63,64,65 but a first-principles dispersion and comparison to experiment see Ref. 51. description of surface energy-loss measurements has not At jellium surfaces, the actual surface-response func- been carried out yet. tion reduces in the long-wavelength limit to Eq. (13) with the Drude ǫ(ω) of Eq. (8), and surface plasmons are therefore expected to be infinitely long-lived excitations. However, energy-loss measurements at simple IV. SUM RULES metal surfaces indicate that surface plasmons exhibit a finite width even at q = 0.60 Since surface plasmons are Sum rules have played a key role in providing insight dictated in this q = 0 limit by bulk properties through in the investigation of a variety of physical situations. A 7 useful sum rule for the surface modes in complementary media with arbitrary geometry was introduced by Apell et al.,66 which in the special case of a metal/vacuum interface implies that67 2 2 2 ωs1 + ωs2 = ωp, (23) where ωs1 is the surface-mode frequency of our system, and ωs2 represents the surface mode of a second complementary system in which the regions of plasma and vacuum are interchanged (see Fig. 5). For example, a half-space filled with a metal of bulk plasma frequency ωp and interfaced with vacuum maps FIG. 6: Fullerene molecule of inner and outer radii r1 and r2. into itself (see Fig. 5), and therefore Eq. (23) yields A Drude dielectric function is assigned to every point between the inner and outer surfaces. ωs1 = ωs2 = ωp/√2, (24) which is the frequency of plasma oscillations at a function g(q,ω),51 which in the long-wavelength limit is metal/vacuum planar interface. given by Eq. (13). This equation yields the classical non- Other examples are a metal sphere in vacuum, which dispersive surface-plasmon frequency ω given by sustains localized Mie plasmons at frequencies s ǫ(ωs)+1=0, (28) l ωl = ωp , (25) which in the case of a semi-infinite Drude metal [see r2l +1 Eq. (8)] is ωs = ωp/√2. with l = 1, 2,..., and a spherical void in a metal, which The long-wavelength surface-plasmon condition of shows Mie plasmons at frequencies Eq. (28) has been derived in the nonrelativistic approxi- l +1 mation, by neglecting retardation of the Coulomb inter- ωl = ωp . (26) action. Hence, Eq. (28) yields a good representation of r2l +1 surface plasma oscillations only at wavelengths that are The squared surface-mode frequencies of the sphere large compared to the Fermi wavelength (q << kF [Eq. (25)] and the void [(Eq. (26)] add up to ω2 for all l, −1 ∼ p 1 A˚ , kF being the Fermi momentum) but small com- as required by Eq. (23). pared to the wavelength of light at optical frequencies Now we consider a situation in which there are two in- −1 (q >>ω /c 0.005 A˚ ). In a typical inelastic electron terfaces, as occurs in the case of a thin film and approx- s scattering experiment,∼ however, the finite angular accep- imately occurs in the case of multishell fullerenes68 and tance of the energy-loss spectrometer garantees that the carbon nanotubes.69 Apell et al.66 have proved a second momentum transfer q be larger than ω /c, so that the sum rule, which relates the surface modes corresponding s retarded region of the surface plasmon dispersion is not to the in-phase and out-of-phase linear combinations of observed.71 the screening charge densities at the interfaces. In the Considering the full set of Maxwell equations and still case of metal/vacuum interfaces this sum rule takes the assuming that the wavelength is long enough for a classi- form of Eq. (23), but now ω 1 and ω 2 being in-phase and s s cal description of the metal/vacuum interface to be jus- out-of-phase modes of the same system. tified (q << k ), one finds that due to retardation the For a Drude metal film with equal and abrupt planar F surface-plasmon condition of Eq. (28) must be replaced surfaces, the actual values of ω 1 and ω 2 are those given s s by71 by Eq. (18). For a spherical fullerene molecule described by assigning a Drude dielectric function to every point ǫ(ωs) 1 + ′ =0, (29) between the inner and outer surfaces of radii r1 and r2 κ(ωs) κ (ωs) (see Fig. 6), one finds the following frequencies for the where in-phase and out-of-phase surface modes:70 ω2 ω2 1 κ(ω)= q2 ǫ(ω) (30) ω2 = p 1 1+4l(l + 1)(r /r )2l+1 , (27) − c2 s 2 ± 2l +1 1 2 r q and which fulfill the sum rule dictated by Eq. (23). 2 ′ 2 ω κ (ω)= q 2 . (31) r − c V. SURFACE-PLASMON POLARITONS In the case of a semi-infinite Drude metal, Eqs. (29)-(31) yield the surface-plasmon dispersion Planar surface plasmons are known to be traced to ω2(q)= ω2/2+ c2q2 ω4/4+ c4q4, (32) the peaks of the imaginary part of the surface response p − p q 8

10 (eV) ω

0 0 0.01 0.02 -1 q (Å )

FIG. 7: The solid line represents the dispersion of surface- plasmon polaritons on a semi-inﬁnite Drude metal with ωp = 15 eV, as obtained from Eq. (32). In the retarded region (q <ωp/c), the surface-plasmon polariton dispersion curve approaches the light line ω = cq (dotted line). At short wavelengths (q >> ω/c), the dispersion curve approaches asymptotically the non-retarded surface-plasmon frequency FIG. 8: Schematic representation of the surface band struc- ωs = ωp/√2 (dashed line). ture on Cu(111) near the Γ¯ point. The shaded region represents the projection of the bulk bands. or, equivalently,

2 2 ω ω ωp − 2D q(ω)= 2 2 , (33) TABLE II: Binding energies (ε ) of surface states at the Γ¯ c s2ω ωp F − point of Be(0001) and the (111) surfaces of the noble metals 2D 2D which we have represented in Fig. 7 by a solid line, to- Cu, Ag, and Au. vF and m represent the corresponding 2D 2D Fermi velocity and effective mass, respectively. vF is gether with the non-retarded surface-plasmon frequency 2 expressed in units of the Bohr velocity v0 = e /¯h. ωs (dashed line) and the light line ω = cq (dotted line). In the retarded region, where q<ω /c, surface plasmons 2D 2D 2D s εF (eV) vF /v0 m couple with the free electromagnetic field, thereby be- Be(0001) 2.8 0.41 1.18 coming what is called a surface-plasmon polariton. In Cu(111) 0.44 0.28 0.42 the non-retarded limit (q >> ωs/c), one finds the non- Ag(111) 0.065 0.11 0.44 dispersive surface-plasmon frequency ωs. Au(111) 0.48 0.35 0.28 Significant deviations from the classical surface- plasmon dispersion of Eq. (32) (like the negative dispersion drawn in Fig. 4), which are typically observed in electron scattering experiments, are only present at wave In the absence of the 3D substrate, Shockley surface vectors larger than those considered in Fig. 7. At q wave states would support a 2D collective oscillation, the en- vectors in the range ω /c<

which is the result first obtained by Chaplik in his study of charge-carrier crystallization in low-density inversion layers.73 If the 2D sheet is located inside the 3D substrate (z 0), however, the long-wavelength limit of the effective≤ 2D dielectric function of Eq. (34) has no zeros at low energies (ω<ωs), due to a complete screening at these energies of electron-electron interactions within the 2D sheet. This result has suggested over the years that acoustic plasmons should only exist in the case of spatially separated plasmas, as pointed out by Das Sarma FIG. 9: Surface-state electrons comprise a 2D sheet of in- and Madhukar.40 teracting free electrons at z = zd. All other states of the Nevertheless, Silkin et al.44 have shown that metal sur- semi-infinite metal comprise a plane-bounded 3D electron gas faces where a partially occupied quasi-2D surface-state at z 0. The metal surface is located at z = 0. ≤ band coexists in the same region of space with the underlying 3D continuum do support a well-defined acoustic surface plasmon, which could not be explained within A. A simple model the local model described above. Furthermore, it has been demonstrated that it is a combination of the nonlo- First of all, we consider a simplified model in which cality of the 3D dynamical screening and the spill out of surface-state electrons comprise a 2D electron gas at z = the 3D electron density into the vacuum which allows the zd (see Fig. 9), while all other states of the semi-infinite formation of 2D electron-density acoustic oscillations at metal comprise a 3D substrate at z 0 represented by ≤ metal surfaces, since these oscilations would otherwise be the Drude dielectric function of Eq. (8). Within this completely screened by the surrounding 3D substrate.46 model, one finds that both e-h and collective excitations occurring within the 2D gas can be described with the use of an effective 2D dielectric function, which in the B. Full calculation random-phase approximation (RPA) takes the form46 In order to achieve a full description of the dynami- 2π ǫ2D (q,ω)=1 + W˜ (z ,z ; q,ω) χ0 (q,ω), (34) cal response of real metal surfaces, we first consider a eff − q d d 2D one-dimensional single-particle potential that describes the main features of the surface band structure.74,75 We W˜ (z,z′; q,ω) being the screened interaction of Eq. (12), 0 then calculate the eigenfunctions and eigenvalues of the and χ2D(q,ω) being the noninteracting density-response 53 corresponding hamiltonian, and we evaluate the dynam- function of a 2D electron gas. ical density-response function χ0(z,z′; q,ω). Finally, we In the absence of the 3D substrate, W˜ (z,z′; q,ω) is 2D solve an integral equation to obtain the RPA interacting simply zero and ǫeff (q,ω) coincides, therefore, with the density-response function χ(z,z′; q,ω). From the knowl- RPA dielectric function of a 2D electron gas, which in edge of this function, which describes bulk and surface the long-wavelength (q 0) limit has one single zero → states on the same footing, one can obtain within linear- corresponding to collective excitations at ω = ω2D. response theory the electron density induced by an ex- In the presence of a 3D substrate that is spatially sep- ternal perturbation φext(z; q,ω): arated from the 2D sheet (zd > 0), the long-wavelength 2D limit of ǫeff (q,ω) has two zeros. One zero corresponds ′ ′ ext ′ 2 2 2 δn(z; q,ω)= dz χ(z,z ; q,ω)φ (z ; q,ω), (37) to a high-frequency oscillation of energy ω = ωs + ω2D in which 2D and 3D electrons oscillate in phase with one Z another. The other zero corresponds to a low-frequency and the collective excitations created by an external po- acoustic oscillation in which both 2D and 3D electrons tential of the form oscillate out of phase. The energy of this low-frequency φext(z; q,ω)= (2π/q)eqz (38) mode is found to be of the form46 −

2D can then be traced to the peaks of the imaginary part of ω = α vF q, (35) the surface-response function51 with α> 1. For small values of the z coordinate (z << d d qz 1), α 1 and the sound velocity approaches, therefore, Im [g(q,ω)] = dz e Im [δn(z; q,ω)] . (39) → the Fermi velocity of the 2D sheet. For zd >> 1, the non- Z interacting 2D density-response function takes the Drude In the bottom panel of Fig. 10 we show the result that 2D 2 form (1/2π)(vF q/ω) , and one ﬁnds we have obtained for the unperturbed electron density of the (0001) surface of Be. We see bulk states, whose total α = √2zd, (36) density extends to the interior of the solid at z < 0, and 10

2 of valence (2s ) electrons in Be (¯hωs = 12.8 eV). The second feature occurs at ω = 0.6eV, corresponding to a new low-energy acoustic collective oscilation, which had been overlooked over the years. In Fig. 11 we show by black lines our full RPA calcula- −1 tion of Im [δn(z; q,ω)] at q =0.05 a0 and the frequencies ω = 12.8 eV (top panel) and ω = 0.6 eV (bottom panel) at which surface collective oscillations occur. At the con- ventional surface-plasmon frequency ω = 12.8 eV we have also carried out a jellium density-functional calculation for a semi-infinite free-electron gas, which we have represented by a grey line in the top panel of Fig. 11. A comparison of our band-structure and jellium calculations (black and grey lines) indicates that the conven- tional surface plasmon is reasonably well described within a jellium model of the surface, although Friedel oscillations in the interior of the solid are considerably more damped in the presence of the actual band structure of the solid. At the acoustic surface-plasmon energy, which for q = −1 0.05 a0 is ω = 0.6eV, a quasi-2D surface-state band in the presence of a 3D substrate yields the new feature displayed in the bottom panel of Fig. 11. Also shown in this figure (shaded area) is the probability density of the partially occupied Shockley surface state, clearly indicating that the low-energy collective excitation at ω = 0.6 eV originates from this 2D surface-state band. Such a 2D electron gas alone would only support a plasmon that −1 for q = 0.05 a0 has energy ω2D = 2.7 eV, well above the low-energy excitation that is visible in the bottom FIG. 10: Unperturbed electron density of the (0001) surface of panel of Fig. 11, and it is only the combination of the Be (dark grey shaded areas), as obtained from a self-consistent strongly localized 2D surface-state band with 3D bulk jellium density-functional calculation (top panel) and from states which allows the formation of this new mode. the use of a one-dimensional potential that describes the main Fig. 12 shows the imaginary part of the surface- features of the surface band structure (bottom panel). The response function g(q,ω) of Be(0001), as obtained from light grey shaded areas represent the neutralizing uniform Eq. (39) for increasing values of q. This figure clearly positive background. The vertical lines and black area in the shows that the excitation spectra is dominated at low en- bottom panel, which are absent in the top panel, represent ergies by a well-defined acoustic peak with linear disper- the atomic positions of the solid and the unperturbed den- sion, the sound velocity being at long wavelengths very sity of occupied surface states, respectively. The crystal edge 2D (z = 0) is chosen to be located half a lattice spacing beyond close to the 2D Fermi velocity vF (see Table II). the last atomic layer, and z < 0 corresponds to the interior of In the bottom panel of Fig. 13 we show the energy of the solid. the acosutic surface plasmon in Be(0001) versus the wave number q (thick solid line), as derived from the max- ima of our calculated surface-loss function Im [g(q,ω)] surface states which are largely localized near the surface. of Fig. 12. The plasmon energy of electrons in an iso- For comparison, we have also carried out a self-consistent lated 2D electron gas would exhibit at long wavelengths jellium density-functional calculation of the electron den- a square-root dispersion with the wave number q [see sity (top panel of Fig. 10), by simply replacing the ions of Eq. (19)]. However, the combination of a 2D surface- the solid by a fixed uniform positive background. Surface state band with the underlying 3D system yields a new states are absent in this model. distinct mode whose energy lies just above the upper edge 2D 2D 2 RPA calculations of the imaginary part of the electron ωu = vF q + q /2 of the continuum of 2D e-h pair ex- density δn(z; q,ω) induced in a Be(0001) surface by the citations (shaded area), where momentum and energy external potential of Eq. (38) were reported in Ref. 44 conservation allows e-h pairs to be created within the 2D −1 ˚ for q = 0.05 a0 (a0 is the Bohr radius, a0 = 0.529 A) electron gas. and a wide range of frequencies. It was demonstrated For a well-defined acoustic surface plasmon to occur, that this quantity exhibits two distinct special features, it must exist for wave vectors q and energies ω where where Im [δn(q,ω)] is maximum near the surface. The decay cannot occur by exciting e-h pairs in the medium. first feature occurs near the surface-plasmon frequency Electron-hole pairs can be excited either within the 2D 11

FIG. 12: Energy-loss function Im [g(q,ω)] /ω of Be(0001) versus the energy ω, as obtained from Eq. (39) for various values −1 of the wave number q, in units of the inverse Bohr radius a0 . The peaks are dictated by the corresponding poles of the surface response function g(q,ω). In the long-wavelenth (q 0) limit, g(q,ω) is simply the total electron density induced→ by the potential of Eq. (38).

surface-state band, or within the 3D continuum of bulk states, or by promoting an electron from an occupied bulk state to an unoccupied 2D surface state. These three mechanisms for the production of e-h pairs at low energies are illustrated in Fig. 13. (i) Since the energy of acoustic surface plasmons always lies just above the up- 2D per edge ωu of the 2D e-h pair continuum (shaded area in the bottom panel of Fig. 13), they cannot possibly decay by creating e-h pairs within the 2D band. (ii) The 3D 3D Fermi velocity vF is typically larger than the Fermi 2D velocity vF of the 2D surface-state band. This means that acoustic surface plasmons can decay by exciting e-h pairs within the 3D continuum of bulk states. However, at the low energies involved the probability for this process to occur is small.46 (iii) An inspection of the upper panel of Fig. 13 shows that, due to the presence of the band gap, for optical (q = 0) transitions to occur from an occupied 3D bulk state to an unoccupied 2D surface state a minimum energy is required, which decreases as the momentum transfer q increases.

FIG. 11: RPA calculation of the imaginary part of the elec- The region of momentum space where transitions from tron density induced in the (0001) surface of Be (black lines), 3D to 2D states cannot occur corresponds to the area be- as obtained from Eqs. (37) and (38) for a wave vector of mag- low the thin solid line in the bottom panel of Fig. 13. This −1 −1 nitude q = 0.05 a0 , as a function of the z coordinate nor- figure shows that at long wavelengths with q < 0.06 a0 mal to the surface. The frequency has been chosen to be acoustic surface plasmons can decay by neither exciting ω = 0.6 eV (bottom panel) and ω = 12.8 eV (top panel). 2D e-h pairs nor exciting 3D-2D e-h pairs, which results The grey line in the top panel represents the corresponding in a very well-defined collective excitation (see Fig. 12). jellium density-functional calculation for a semi-infinite free- −1 At shorter wavelengths (q > 0.06 a0 ), the promotion electron gas. The probability density of the partially occupied of electrons from occupied 3D bulk states to unoccupied Shockley surface state is represented in the bottom panel by a 2D surface states becomes possible and the corresponding shaded area. As in Fig. 10, the crystal edge (z = 0) is chosen to be located half a lattice spacing beyond the last atomic plasmon peak broadenes considerably. layer, and z < 0 corresponds to the interior of the solid. 12

FIG. 14: Top panel: Schematic (out of scale) representation of typical energy dispersions of acoustic surface plasmons, acoustic phonons, and free-space electromagnetic radiation. In the actual scale, the dispersion lines of free light and acoustic phonons should be very close to the vertical and horizontal axes, respectively. Bottom panel: Periodic grating of constant L. The grating periodic structure can provide the impigning free electromagnetic radiation with additional momentum 2π/L.

C. Excitation of acoustic surface plasmons

We close this paper by discussing whether acoustic surface plasmons can be observed. As in the case of conven- tional surface plasmons, acoustic surface plasmons should be expected to be excited by either electrons or light. Here we focus on a possible mechanism that would lead to the excitation of acoustic surface plasmons by light in, e.g., vicinal surfaces with high indices. The top panel of Fig. 14 exhibits the energy dispersion of acoustic surface plasmons at low wave vectors. Also shown in this figure are the light line ω = cq and, for com- FIG. 13: Top panel: Schematic representation of the projec- parison, a typical energy dispersion of acoustic phonons. tion of the bulk bands onto a solid surface that supports a par- The sound velocity of acoustic surface plasmons, which is tially occupied Shockley surface-state band. The energy of oc- close to the Fermi velocity of the 2D surface-state band, is cupied and unoccupied states is displayed as a function of the typically a few orders of magnitude larger than the sound momentum parallel to the surface. The solid line represents a velocity of acoustic phonons in metals but still about two Shockley surface state. Dark and light shaded areas represent orders of magnitude smaller than the velocity of light. occupied and unoccupied bulk states, respectively. The white The acoustic surface plasmon dispersion curve is well area represents the band gap. Electron-hole pair excitations below the dispersion curve of free-space electromagnetic are represented by arrows, depending on whether they correspond to transitions within the surface-state band (2D), tran- radiation. Hence, there is, in principle, no way that in- sitions within the bulk bands (3D), or transitions from bulk cident light can provide an ideal surface with the correct to surface states (2D-3D). Bottom panel: Dispersion of the amount of momentum and energy for the excitation of acoustic surface collective excitation of Be(0001), as derived an acoustic surface plasmon. As in the case of conven- from the peaks that are visible in Fig. 12 (thick solid line). tional surface plasmons, however, a periodic corrugation 2D 2D 2 This curve stays just over the upper edge ωu = vF q + q /2 or grating in the metal surface should be able to provide of the 2D e-h pair continuum (shaded area). Momentum and the missing momentum. energy conservation prevent 2D e-h pairs from being produced 2D If light hits a surface with a periodic corrugation, the for energies above ωu . The area below the thin solid line represents the region of momentum space where transitions from 3D to 2D states cannot occur. Well-defined acoustic plasmons are expected to occur at wave vectors with magnitude smaller −1 than q 0.06 a0 . ∼ 13 grating (see bottom panel of Fig. 14) can provide the mons of energy ω 0.6 eV could be excited in this way. impigning free electromagnetic waves with additional mo- Although a grating∼ period of a few nanometers sounds mentum arising from the grating periodic structure. If unrealistic with present technology, the possible control free electromagnetic radiation hits the grating at an angle of vicinal surfaces with high indices could provide appro- θ, its wave vector along the grating surface has magni- priate grating periods in the near future. Alternatively, tude acoustic surface plasmons could be observed with the use of high-resolution electron energy-loss spectroscopy ω 2π q = sin θ n, (40) (EELS) under grazing incidence. c ± L where L represents the grating constant, and n =1, 2,.... Hence, the linear (nearly vertical) dispersion relation of Acknowledgments free light changes into a set of parallel straight lines, which can match the acoustic-plasmon dispersion relation. Partial support by the University of the Basque Coun- For a well-defined acoustic surface plasmon in Be(0001) try, the Basque Unibertsitate eta Ikerketa Saila, the to be observed, the wave number q needs to be smaller Spanish Ministerio de Educación y Cultura, and the −1 −1 than q 0.06 a0 . For q =0.05 a0 , Eq. (40) with n =1 Max Planck Research Award Funds is gratefully acknowl- yields a∼ grating constant L = 66 A.˚ Acoustic surface plas- edged.

1 D. Pines and D. Bohm, Phys. Rev. 85, 338 (1952). 26 A. Otto, Z. Phys. 216, 398 (1968). 2 D. Pines, Rev. Mod. Phys. 28, 184 (1956). 27 E. Kretschmann and H. Raether, Z. Naturforsch. A 23, 3 The valence electron density n is usually characterized by 2135 (1968). 3 28 the density parameter rs = (3/4πn) /a0 (a0 is the Bohr J. R. Sambles, G. W. Bradbery, and F. Z. Yang, Contemp. radius, a0 = 0.529 A).˚ In metals (2 < rs < 6), plasmons Phys. 32, 173 (1991). 29 have energies in the range 5eV <ωp < 20eV and frequen- R. Berndt, J. K. Gimzewski, and P. Johansson, Phys. Rev. cies that lie, therefore, in the optical regime. 67, 3796 (1991). 4 D. Gabor, Phil. Mag. 1, 1 (1956). 30 M. J. Shea and R. N. Compton, Phys. Rev. B 47, 9967 5 R. H. Ritchie, Phys. Rev. 106, 874 (1957). (1993). 6 E. A. Stern and R. A. Ferrell, Phys. Rev. 120, 130 (1960). 31 M. Rocca, Surf. Sci. Rep. 22, 1 (1995). 7 C. J. Powell and J. B. Swan, Phys. Rev. 115, 869 (1959); 32 B. Rothenhausler and W. Knoll, Nature 332, 615 (1988). 116, 81 (1959). 33 S. C. Schuster, R. V. Swanson, L. A. Alex, R. B. Bourret, 8 P. J. Feibelman, Surf. Sci. 27, 438 (1971). and M. I. Simon, Nature 365, 343 (1993). 9 R. H. Ritchie, Phys. Lett. A 38, 189 (1972). 34 J. B. Pendry, Science 265, 1687 (1999). 10 R. Ray amd G. D. Mahan, Phys. Lett. A 42, 306 (1972). 35 T. W. Ebbesen, H. J. Lezec, H. F. Ghaemi, T. Thio, and 11 A. Sunjic, G. Toulouse, and A. A. Lucas, Solid State Com- P. A. Wolﬀ, Nature 391, 667 (1998). mun. 11, 1629 (1972). 36 H. J. Lezec, A. Degiron, E. Devaux, R. A. Linke, L. Martin- 12 J. Schmit and A. A. Lucas, Solid State Commun. 11, 415 Moreno, F. J. Garc´ıa-Vidal, and T. W. Ebbesen, Science (1972). 297, 820 (2002). 13 P. M. Echenique and J. B. Pendry, J. Phys. C: Solid State 37 W. L. Barnes, A. Dereux, and T. W. Ebbesen, Nature 424, Phys. 8 (1975). 824 (2003). 14 P. M. Echenique, R. H. Ritchie, N. Barberan, and J. Ink- 38 D. Pines, Can. J. Phys. 34, 1379 (1956). son, Phys. Rev. B 23, 6486 (1981). 39 N. H. March and M. P. Tosi, Adv. Phys. 44, 299 (1995). 15 P. J. Feibelman, Prog. Surf. Sci. 12, 287 (1982). 40 S. Das Sarma and A. Madhukar, Phys. Rev. B 23, 805 16 R. H. Ritchie, Prog. Theor. Phys. 29, 607 (1963). (1981). 17 R. H. Ritchie and A. L. Marusak, Surf. Sci. 4, 234 (1966). 41 D. Olego, A. Pinczuk, A. C. Gossard, and W. Wiegmann, 18 D. E. Beck, Phys. Rev. B 4, 1555 (1971). Phys. Phys. Rev. B 32, 1921 (1985). 19 P. Hohenber and W. Kohn, Phys. Rev. 136 B864 (1964); 42 J. Ruvalds, Phys. Rev. B 35, 8869 (1987); Nature 328, W. Kohn and L. J. Sham, Phys. Rev. A1133 (1965). 299 (1987). 20 K.-D. Tsuei, E. W. Plummer, A. Liebsch, K. Kempa, and 43 A. Bill, H. Morawitz, and V. Z. Kresin, Phys. Rev. B 66, P. Bakshi, Phys. Rev. Lett. 64, 44 (1989). 100501 (2002). 21 W. Plummer, K.-D. Tsuei, and B.-O. Kim, Nucl. Instrum. 44 V. M. Silkin, A. Garc´ıa-Lekue, J. M. Pitarke, E. V. Methods B 96, 448 (1995). Chulkov, E. Zaremba, and P. M. Echenique, Eurphys. Lett. 22 A. J. Benett, Phys. Rev. B 1, 203 (1970). 66, 260 (2004). 23 Throughout this paper, long wavelengths are associated 45 The sound velocity of this acoustic mode is, however, close with low wave vectors. to the Fermi velocity of the 2D surface-state band, which is 24 R. H. Ritchie and H. B. Eldridge, Phys. Rev. 126, 1935 typically a few orders of magnitude larger than the sound (1962). velocity of acoustic phonons in metals but still about two 25 Y.-Y. Teng and E. A. Stern, Phys. Rev. Lett. 19, 511 orders of magnitude smaller than the velocity of light. (1967). 46 J. M. Pitarke, V. U. Nazarov, V. M. Silkin, E. V. Chulkov, 14

E. Zaremba, and P. M. Echenique (unpublished). Kempa, and P. Bakshi, Surf. Sci. 247, 302 (1991). 47 This approximation is valid for heavy charged particles, 60 P. T. Sprunger, G. M. Watson, and E. W. Plummer, Surf. e.g., ions, and also for swift electrons moving with nonrel- Sci. 269, 551 (1992). ativistic velocities that are large compared to the velocity 61 E. Zaremba and W. Kohn, Phys. Rev. B 13, 2270 (1976). 62 of target electrons. In the case of electrons, Z1 = 1. P. Apell and C. Holmberg, Solid State Commun. 49, 693 48 F. Flores and F. Garc´ıa Moliner, Introduction to the− Theory (1984). of Solid Surfaces (Cambridge University Press, Cambridge, 63 J. Tarriba and W. L. Mochan, Phys. Rev. 46, 12902 (1992). 1979). 64 P. J. Feibelman, Surf. Sci. 282, 129 (1993). 49 A. L. Fetter and J. D. Wallecka, Quantum Theory of Many- 65 A. Liebsch, Phys. Rev. Lett. 71, 145 (1993). Particle Systems (McGraw-Hill, New York, 1964). 66 S. P. Apell, P. M. Echenique, and R. H. Ritchie, Ultrami- 50 J. Lindhard, K. Dan. Vidensk. Selsk. Mat. Fys. Medd. 28, croscopy 65, 53 (1995). 8 (1954). 67 A. Ronveaux, A. Moussiaux, and A. A. Lucas, Can. J. 51 See, e.g., A. Liebsch, Electronic Excitations at Metal Sur- Phys. 55, 1407 (1977). faces (Plenum Press, New York, 1997). 68 H. W. Kroto, J. R. Heath, S. C. O’Brien, R. F. Curl, and 52 R. A. Ferrell, Phys. Rev. 111, 1214 (1958). R. E. Smalley, Nature 318, 162 (1985). 53 F. Stern, Phys. Rev. Lett. 18, 546 (1967); T. Ando, A. B. 69 S. Ijima, Nature 354, 56 (1991). Fowler, and F. Stern, Rev. Mod. Phys. 54, 437 (1982). 70 D. Ostling,¨ P. Apell, and A. Rosén, Europhys. Lett. 21, 54 S. J. Allen jr., D. C. Tsui, and R. A. Logan, Phys. Rev. 539 (1993). Lett. 38, 980 (1977). 71 H. Raether, Surface Plasmons, Springer Tracks in Modern 55 T. Nagao, T. Heldebrandt, M. Henzler, and S. Hasegava, Physics, 111 (New York, Springer, 1988). Phys. Rev. Lett. 86, 5747 (2001). 72 J. E. Inglesfield, Rep. Prog. Phys. 45, 223 (1982). 56 M. Abramowitz and I. A. Stegun, Handbokk of Mathemat- 73 A. V. Chaplik, Zh. Eksp. Teor. Diz. 62, 746 (1972) (Sov. ical Functions (Dover, New York, 1965). Phys. JETP, 35, 395 (1972)). 57 R. Nuñez, P. M. Echenique, and R. H. Ritchie, J. Phys. C: 74 E. V. Chulkov, V. M. Silkin, and P. M. Echenique, Surf. Solid State Phys. 13, 4229 (1980). Sci. 437, 330 (1999). 58 H. Ibach and D. L. Mills, Electron Energy Loss Spec- 75 P. M. Echenique, J. M. Pitarke, E. V. Chulkov, and A. troscopy and Surface Vibrations (Academic, New York, Rubio, Chem. Phys. 251, 1 (2000). 1982). 59 K.-D. Tsuei, E. W. Plummer, A. Liebsch, E. Pehlke, K.