Laser fields at flat interfaces: II. Plasmon resonances in aluminium photoelectron spectra

Using the model derived in paper I [G. Raşeev, Eur. Phys. J. D 66, 167 (2012)], this work presents calculations of the photoelectron spectrum (PES) of low index aluminium surfaces in the 10–30 eV region. The laser is p or transverse magnetic linearly polarized incident on a flat structureless surface and its fields are modeled in I using the vector potential in the temporal gauge. This model uses a tensor and non-local isotropic (TNLI) susceptibility and solves the classical Ampère-Maxwell equation through the use of the vector potential from the electron density-coupled integro-differential equations (VPED-CIDE). The PE cross sections are the squares of the PE transition moments calculated using the VPED-CIDE vector potential function of the penetration coordinate. The PES is obtained in a one step model using either the Fermi golden rule or the Weisskopf-Wigner (WW) expressions. The WW cross section PES compares favorably with the experimental angle and energy resolved photoelectron yield (AERPY) spectrum of Levinson et al. [Phys. Rev. Lett. 43, 952 (1979)], Levinson and Plummer [Phys. Rev. B 24, 628 (1981)] for Al(001) and of Barman et al. [Phys. Rev. B 58, R4285 (1998)], Barman [Curr. Sci. 88, 54 (2005)] for Al(111) surfaces. As in the experiment, our theoretical AERPY displays the multipole surface plasmon resonance at 11.32/12.75 eV for Al(001)/Al(111), mainly due to the surface contribution |⟨ψf|p·A|ψi⟩|2, the bulk plasmon minimum at 15 eV and the two single particle excitation resonances at about 16 and 22 eV. The nature of the plasmon resonances of the PES is analyzed using the reflectance, the electron density induced by the laser and Feibelman’s parameter d⊥ all introduced in paper I