Angle-resolved photoemission spectroscopy

From Wikipedia, the free encyclopedia
Jump to navigation Jump to search

ARPES spectrum of a two-dimensional electronic state localized at the (111) surface of copper. The energy has free-electron-like momentum dependence, p2/2m, where m=0.46me. Color scale represents electron counts per kinetic energy and emission angle channel. When 21.22 eV photons are used, the Fermi level is imaged at 16.64 eV. T=300K.

Angle-resolved photoemission spectroscopy (ARPES) is a powerful technique used in condensed matter physics to probe the structure of the electrons in a material, usually a crystalline solid. The technique is best suited for use in one- or two-dimensional materials. It is based on the photoelectric effect, in which an incoming photon of sufficient frequency dislodges an electron from the surface of a material. By directly measuring the kinetic energy and momentum distributions of the emitted photoelectrons, the technique can be used to map the electronic band structure, provide elemental information, and map Fermi surfaces. ARPES has been used by physicists to investigate high-temperature superconductors and materials exhibiting charge density waves.

The main components of an ARPES system are a source to deliver a high-frequency monochromatic beam of photons, a sample holder connected to a manipulator used to position and manipulate the material, and an electron spectrometer. The equipment is contained within an ultra-high vacuum (UHV) environment, which protects the sample and prevents the emitted electrons from being scattered. After being dispersed, the electrons are directed to a microchannel plate detector, which is linked to a camera. Energy dispersion is carried out for a narrow range of energies around the pass energy, which enables the electrons to reach the detector.

Some ARPES systems have an electron extraction tube alongside the detector, which measures the electrons’ spin polarization. Systems that use a slit can only make angular maps in one direction. For two-dimensional maps, the sample is rotated, or the electrons are manipulated.


Typical laboratory setup of an ARPES experiment: Helium discharge lamp as an ultraviolet light source, sample holder that attaches to a vacuum manipulator, and hemispherical electron energy analyzer.

A typical instrument for angle-resolved photoemission consists of a light source, a sample holder attached to a manipulator, and an electron spectrometer. These are all part of an ultra-high vacuum system that provides the necessary protection from adsorbates for the sample surface and eliminates scattering of the electrons on their way to the analyzer.[1][2]

The light source delivers a monochromatic, usually polarized, focused, high-intensity photon beam to the sample (~1012 photons/s with a few meV energy spread).[2] Light sources range from compact noble-gas discharge UV lamps and radio-frequency plasma sources (10–⁠40 eV),[3][4][5] ultraviolet lasers (5–⁠11 eV)[6] to synchrotron[7] insertion devices that are optimized for different parts of the electromagnetic spectrum (from 10 eV in the ultraviolet to 1000 eV X-rays).

The sample holder accommodates samples of crystalline materials, the electronic properties of which are to be investigated, and facilitates their insertion into the vacuum, cleavage to expose clean surfaces, precise manipulation as the extension of a manipulator (for translations along three axes, and rotations to adjust the sample's polar, azimuth and tilt angles), precise temperature measurement and control, cooling to temperatures as low as 1 kelvin with the help of cryogenic liquefied gases, cryocoolers, and dilution refrigerators, heating by resistive heaters to a few hundred °C or by backside electron-beam bombardment for temperatures up to 2000 °C, and light beam focusing and calibration.

Electron trajectories in an ARPES spectrometer shown in the plane of angular dispersion. The instrument shows a certain degree of focusing on the same detection channel of the electrons leaving the crystal at the same angle but originating from two separate spots on the sample. Here, the simulated separation is 0.5 mm.

The electron spectrometer disperses along with two spatial directions the electrons reaching its entrance concerning their kinetic energy and their emission angle when exiting the sample. In the type most commonly used, the hemispherical electron energy analyzer, the electrons first pass through an electrostatic lens that picks electrons emitted from its own small focal spot on the sample (conveniently located some 40 mm from the entrance to the lens), enhances the angular spread of the electron plume, and serves it to the narrow entrance slit of the energy dispersing element with adjusted energy.

Angle- and energy-resolving electron spectrometer for ARPES

The energy dispersion is carried out for a narrow range of energies around the so-called pass energy in the direction perpendicular to the slit, typically 25 mm long and >0.1 mm wide. The angular dispersion of the cylindrical lens is only preserved along the slit, and depending on the lens model and the desired angular resolution can amount to ±3°, ±7° or ±15°.[3][4][5] The hemispheres of the energy analyzer are kept at constant voltages so that the central trajectory is followed by electrons that have the kinetic energy equal to the set pass energy; those with higher or lower energies end up closer to the outer or the inner hemisphere at the other end of the analyzer. This is where an electron detector is mounted, usually in the form of a 40 mm microchannel plate paired with a fluorescent screen. Electron detection events are recorded using an outside camera and are counted in hundreds of thousands of separate angle vs. kinetic energy channels. Some instruments are additionally equipped with an electron extraction tube at one side of the detector to enable the measurement of electrons spin polarization.

Modern analyzers are capable of resolving the electron emission angles as low as nearly 0.1°. Energy resolution is pass-energy and slit-width dependent so the operator chooses between measurements with ultrahigh resolution and low intensity (<1 meV at 1 eV pass energy) or poorer energy resolutions of 10 or more meV at higher pass energies and with wider slits resulting in higher signal intensity. The instrument's resolution shows up as artificial broadening of the spectral features: a Fermi energy cutoff wider than expected from the sample's temperature, and the theoretical electron's spectral function convolved with the instrument's resolution function in both energy and momentum/angle.[3][4][5]

Sometimes, instead of hemispherical analyzers, time-of-flight analyzers are used. These, however, require pulsed photon sources and are most common in laser-based ARPES labs.[8]

Left: Analyzer angle - Energy map I0(α,Ek) around vertical emission. Right: Analyzer angle - Energy maps Iθ(α,Ek) at several polar angles away from vertical emission
Left: Constant energy map near EF in analyzer angle - polar angle units (polar motion perpendicular to analyzer slit). Right: Constant energy map near EF in crystal momentum units (transformed from the analyzer angle - polar angle map)



Angle-resolved photoemission spectroscopy is a potent refinement of ordinary photoemission spectroscopy. Photons with a frequency have an energy , defined by the equation:

where is Planck's constant.[9]

A photon is used to stimulate the transition of an electron from an occupied to unoccupied electronic state of the solid. If the photon's energy is greater than the electron's binding energy , the electron will eventually be emitted with a characteristic kinetic energy and angle relative to the surface normal. The kinetic energy is given by:


Electron emission intensity maps can be produced from these results. The maps represent the intrinsic distribution of electrons in the solid and are expressed in terms of and the Bloch wave is described by the wave vector , which is related to the electrons' crystal momentum and group velocity. In the process, the Bloch wave vector is linked to the measured electron's momentum , where the magnitude of the momentum, is given by the equation:


Only the component that is parallel to the surface is preserved. The component of the wave vector parallel to the direction of the crystal lattice is related to the parallel component of the momentum and , the reduced Planck constant, by the expression:

This component is known, and its magnitude is given by:


Because of this,[vague] and its pronounced surface sensitivity, ARPES is best suited to the complete characterization of the band structure in ordered low-dimensional systems such as two-dimensional materials, ultrathin films, and nanowires. When it is used for three-dimensional materials, the perpendicular component of the wave vector is usually approximated, with the assumption of a parabolic, free-electron-like final state with the bottom at energy . This gives:


Fermi surface mapping[edit]

Electron analyzers that need a slit to prevent the mixing of momentum and energy channels are only capable of taking angular maps along one direction. To take maps over energy and two-dimensional momentum space, either the sample is rotated in the proper direction so that the slit receives electrons from adjacent emission angles, or the electron plume is steered inside the electrostatic lens with the sample fixed. The slit width will determine the step size of the angular scans: if a 30 mm long slit is served with a 30° plume, this will, in the narrower (say 0.5 mm) direction of the slit average signal over a 0.5mm by 30°/30mm, that is, 0.5° span, which will be the maximal resolution of the scan in that other direction. Coarser steps will lead to missing data, and a finer step to overlaps. The energy-angle-angle maps can be further processed to give energy-kx-ky maps, and sliced in such a way to display constant energy surfaces in the band structure and most importantly the Fermi surface map when cutting near the Fermi level.

Emission angle to momentum conversion[edit]

Geometry of an ARPES experiment. In this position, ϑ=0° & τ=0°, the analyzer is accepting electrons emitted vertically from the surface and α≤8° around.

ARPES spectrometer measures angular dispersion in a slice α along its slit. Modern analyzers record these angles simultaneously, in their reference frame, typically in the range of ±15°.[3][4][5] To map the band structure over a two-dimensional momentum space, the sample is rotated while keeping the light spot on the surface fixed. The most common choice is to change the polar angle ϑ around the axis that is parallel to the slit and adjust the tilt τ or azimuth φ so emission from a particular region of the Brillouin zone can be reached. The measured electrons have these momentum components in the reference frame of the analyzer , where . The reference frame of the sample is rotated around the y axis by ( there has components ), then tilted around x by τ, resulting in . Here, are appropriate rotation matrices. The components of the electron's crystal momentum known from ARPES in this mapping geometry are thus

choose sign at depending on whether is proportional to or

If high symmetry axes of the sample are known and need to be aligned, a correction by azimuth φ can be applied by rotating around z, or by rotating the map I(E, kx, ky) around origin in two-dimensional momentum planes.

Theoretical derivation of intensity relationship[edit]

The theory of photoemission[1][10][12] is that of direct optical transitions between the states and of an N-electron system. Light excitation is introduced as the magnetic vector potential through the minimal substitution in the kinetic part of the quantum-mechanical Hamiltonian for the electrons in the crystal. The perturbation part of the Hamiltonian comes out to be:


In this treatment, the electron's spin coupling to the electromagnetic field is neglected. The scalar potential set to zero either by imposing the Weyl gauge [1] or by working in the Coulomb gauge in which becomes negligibly small far from the sources. Either way, the commutator is taken to be zero. Specifically, in Weyl gauge because the period of for ultraviolet light is about two orders of magnitude larger than the period of the electron's wave function. In both gauges it is assumed the electrons at the surface had little time to respond to the incoming perturbation and add nothing to either of the two potentials. It is for most practical uses safe to neglect the quadratic term. Hence, .

The transition probability is calculated in time-dependent perturbation theory and is given by the Fermi's golden rule:


The delta distribution above says that energy is conserved when a photon of energy is absorbed .

If the electric field of an electromagnetic wave is written as , where , the vector potential holds its polarization and equals to . The transition probability is then given in terms of the electric field as[13]


In the sudden approximation, which assumes an electron is instantaneously removed from the system of N electrons, the final and initial states of the system are taken as properly antisymmetrized products of the single particle states of the photoelectron , and the states representing the remaining N-1 electron systems.[1]

The photoemission current of electrons of energy and momentum is then expressed as the products of

  • , known as the dipole selection rules for optical transitions, and
  • , the one-electron removal spectral function known from the many-body theory of condensed matter physics

summed over all allowed initial and final states leading to the energy and momentum being observed.[1] Here, E is measured with respect to the Fermi level EF and Ek with respect to vacuum so where , the work function, is the energy difference between the two referent levels that is material, surface orientation, and surface condition dependent. Because the allowed initial states are only those that are occupied, the photoemission signal will reflect the Fermi-Dirac distribution function in the form of a temperature-dependent sigmoid-shaped drop of intensity in the vicinity of EF. In the case of a two-dimensional, one-band electronic system the intensity relation further reduces to .[1]

Selection rules[edit]

The electronic states in crystals are organized in energy bands, which have associated energy-band dispersions that are energy eigenvalues for delocalized electrons according to Bloch's theorem. From the plane-wave factor in Bloch's decomposition of the wave functions, it follows the only allowed transitions when no other particles are involved are between the states whose crystal momenta differ by the reciprocal lattice vectors , i.e. those states that are in the reduced zone scheme one above another (thus the name direct optical transitions).[12]

Another set of selection rules comes from (or ) when the photon polarization contained in (or ) and symmetries of the initial and final one-electron Bloch states and are taken into account. Those can lead to the suppression of the photoemission signal in certain parts of the reciprocal space or can tell about the specific atomic-orbital origin of the initial and final states.[14]

Many-body effects[edit]

ARPES spectrum of the renormalized π band of electron-doped graphene; p-polarized 40eV light, T=80K. Dotted line is the bare band. The kink at -0.2 eV is due to graphene's phonons.[15]

The one-electron spectral function that is directly measured in ARPES maps the probability the state of the system of N electrons from which one electron has been instantly removed is any of the ground states of the N−1 particle system:


If the electrons were independent of one another, the N electron state with the state removed would be exactly an eigenstate of the N−1 particle system and the spectral function would become an infinitely sharp delta function at the energy and momentum of the removed particle; it would trace the dispersion of the independent particles in energy-momentum space. In the case of increased electron correlations, the spectral function broadens and starts developing richer features that reflect the interactions in the underlying many-body system. These are customarily described by the complex correction to the single particle energy dispersion that is called the quasiparticle self energy, . It contains the full information about the renormalization of the electronic dispersion due to interactions and the lifetime of the hole created by the excitation. Both can be determined experimentally from the analysis of high-resolution ARPES spectra under a few reasonable assumptions. Namely, one can assume that the part of the spectrum is nearly constant along high-symmetry directions in momentum space and that the only variable part comes from the spectral function, which in terms of , where the two components of are usually taken to be only dependent on , reads

Constant energy cuts of the spectral function are approximately Lorentzians whose width at half maximum is determined by the imaginary part of the self energy, while their deviation from the bare band is given by its real part.

This function is known from ARPES as a scan along a chosen direction in momentum space and is a two-dimensional map of the form . When cut at a constant energy , a Lorentzian-like curve in is obtained whose renormalized peak position is given by and whose width at half maximum is determined by , as follows:[16][15]

The only remaining unknown in the analysis is the bare band . The bare band can be found in a self-consistent way by enforcing the Kramers-Kronig relation between the two components of the complex function that is obtained from the previous two equations. The algorithm is as follows: start with an ansatz bare band, calculate by eq. (2), transform it into using the Kramers-Kronig relation, then use this function to calculate the bare band dispersion on a discrete set of points by eq. (1), and feed to the algorithm its fit to a suitable curve as a new ansatz bare band; convergence is usually achieved in a few quick iterations.[15]

From this obtained self-energy, one can judge on the strength and shape of electron-electron correlations, electron-phonon (more generally, electron-boson) interaction, active phonon energies, and quasiparticle lifetimes.[17][18][19][20][21]

In simple cases of band flattening near the Fermi level because of the interaction with Debye phonons, the band mass is enhanced by (1+λ) and the electron-phonon coupling factor λ can be determined from the linear dependence of the peak widths on temperature.[20]


ARPES has been used to map the occupied band structure of many metals and semiconductors, states appearing in the projected band gaps at their surfaces,[10] quantum well states that arise in systems with reduced dimensionality,[22] one-atom-thin materials like graphene[23] transition metal dichalcogenides, and many flavors of topological materials.[24][25] It has also been used to map the underlying band structure, gaps, and quasiparticle dynamics in highly correlated materials like high-temperature superconductors and materials exhibiting charge density waves.[1][26][27][8]

When the electron dynamics in the bound states just above the Fermi level need to be studied, two-photon excitation in pump-probe setups (2PPE) is used. There, the first photon of low-enough energy is used to excite electrons into unoccupied bands that are still below the energy necessary for photoemission (i.e. between the Fermi and vacuum levels). The second photon is used to kick these electrons out of the solid so they can be measured with ARPES. By precisely timing the second photon, usually by using frequency multiplication of the low-energy pulsed laser and delay between the pulses by changing their optical paths, the electron lifetime can be determined on the scale below picoseconds.[28][29]

External links[edit]


  1. ^ a b c d e f g Damascelli, Andrea; Shen, Zhi-Xun; Hussain, Zahid (April 17, 2003). "Angle-resolved photoemission spectroscopy of the cuprate superconductors". Reviews of Modern Physics. 75 (2): 473–541. arXiv:cond-mat/0208504. doi:10.1103/RevModPhys.75.473. ISSN 0034-6861. S2CID 118433150.
  2. ^ a b Hüfner, Stefan, ed. (2007). Very High Resolution Photoelectron Spectroscopy. Lecture Notes in Physics. 715. Berlin, Heidelberg: Springer Berlin Heidelberg. doi:10.1007/3-540-68133-7. ISBN 978-3-540-68130-4. (subscription required)
  3. ^ a b c d "MBScientific electron analysers and UV sources".
  4. ^ a b c d "ARPES Lab". Scienta Omicron. 2020. Retrieved August 29, 2020.
  5. ^ a b c d "Lab ARPES System with PHOIBOS Analyzer". SPECS. Retrieved August 29, 2020.
  6. ^ "Products". Lumeras LLC. 2013. Retrieved August 29, 2020.
  7. ^ "Light sources of the world".
  8. ^ a b Zhou, Xingjiang; He, Shaolong; Liu, Guodong; Zhao, Lin; Yu, Li; Zhang, Wentao (June 1, 2018). "New Developments in Laser-Based Photoemission Spectroscopy and its Scientific Applications: a Key Issues Review". Reports on Progress in Physics. 81 (6): 062101. arXiv:1804.04473. Bibcode:2018RPPh...81f2101Z. doi:10.1088/1361-6633/aab0cc. ISSN 0034-4885. PMID 29460857. S2CID 3440746.
  9. ^ Soper, Davison E. "Electromagnetic radiation is made of photons". Retrieved September 3, 2020.
  10. ^ a b c Hüfner, Stefan. (2003). "Introduction and Basic Principles". Photoelectron Spectroscopy: Principles and Applications (Third rev. and enlarged ed.). Berlin, Heidelberg: Springer Berlin Heidelberg. ISBN 978-3-662-09280-4. OCLC 851391282.
  11. ^ Damascelli, Andrea; Shen, Zhi-Xun; Hussain, Zahid (April 17, 2003). "Angle-resolved photoemission spectroscopy of the cuprate superconductors". Reviews of Modern Physics. 75 (2): 473–541. arXiv:cond-mat/0208504. doi:10.1103/RevModPhys.75.473. ISSN 0034-6861. S2CID 118433150.
  12. ^ a b Damascelli, Andrea (2004). "Probing the Low-Energy Electronic Structure of Complex Systems by ARPES". Physica Scripta. T109: 61. arXiv:cond-mat/0307085. doi:10.1238/Physica.Topical.109a00061. ISSN 0031-8949. S2CID 21730523.
  13. ^ Wacker, Andreas. "Fermi's golden rule" (PDF). Teaching notes (Lund University).
  14. ^ Cao, Yue; Waugh, J. A.; Zhang, X.-W.; Luo, J.-W.; Wang, Q.; Reber, T. J.; Mo, S. K.; Xu, Z.; Yang, A.; Schneeloch, J.; Gu, G. (July 21, 2013). "In-Plane Orbital Texture Switch at the Dirac Point in the Topological Insulator Bi2Se3". Nature Physics. 9 (8): 499–504. arXiv:1209.1016. doi:10.1038/nphys2685. ISSN 1745-2473.
  15. ^ a b c Pletikosić, Ivo; Kralj, Marko; Milun, Milorad; Pervan, Petar (April 24, 2012). "Finding the bare band: Electron coupling to two phonon modes in potassium-doped graphene on Ir(111)". Physical Review B. 85 (15): 155447. arXiv:1201.0777. Bibcode:2012PhRvB..85o5447P. doi:10.1103/PhysRevB.85.155447. ISSN 1098-0121. S2CID 119170154.
  16. ^ Kordyuk, A. A.; Borisenko, S. V.; Koitzsch, A.; Fink, J.; Knupfer, M.; Berger, H. (June 9, 2005). "Bare electron dispersion from photoemission experiments". Physical Review B. 71 (21): 214513. arXiv:cond-mat/0405696. doi:10.1103/PhysRevB.71.214513. ISSN 1098-0121. S2CID 67784336.
  17. ^ Norman, M. R.; Ding, H.; Fretwell, H.; Randeria, M.; Campuzano, J. C. (September 1, 1999). "Extraction of the Electron Self-Energy from Angle Resolved Photoemission Data: Application to Bi2212". Physical Review B. 60 (10): 7585–7590. arXiv:cond-mat/9806262. doi:10.1103/PhysRevB.60.7585. ISSN 0163-1829. S2CID 4691468.
  18. ^ LaShell, S.; Jensen, E.; Balasubramanian, T. (January 15, 2000). "Nonquasiparticle structure in the photoemission spectra from the Be(0001) surface and determination of the electron self energy". Physical Review B. 61 (3): 2371–2374. Bibcode:2000PhRvB..61.2371L. doi:10.1103/PhysRevB.61.2371. ISSN 0163-1829. (subscription required)
  19. ^ Valla, T.; Fedorov, A. V.; Johnson, P. D.; Hulbert, S. L. (September 6, 1999). "Many-Body Effects in Angle-Resolved Photoemission: Quasiparticle Energy and Lifetime of a Mo(110) Surface State". Physical Review Letters. 83 (10): 2085–2088. arXiv:cond-mat/9904449. Bibcode:1999PhRvL..83.2085V. doi:10.1103/PhysRevLett.83.2085. ISSN 0031-9007. S2CID 55072153.
  20. ^ a b Hofmann, Ph; Sklyadneva, I Yu; Rienks, E D L; Chulkov, E V (December 11, 2009). "Electron–phonon coupling at surfaces and interfaces". New Journal of Physics. 11 (12): 125005. Bibcode:2009NJPh...11l5005H. doi:10.1088/1367-2630/11/12/125005. ISSN 1367-2630.
  21. ^ Veenstra, C. N.; Goodvin, G. L.; Berciu, M.; Damascelli, A. (July 16, 2010). "Elusive electron-phonon coupling in quantitative analyses of the spectral function". Physical Review B. 82 (1): 012504. arXiv:1003.0141. Bibcode:2010PhRvB..82a2504V. doi:10.1103/PhysRevB.82.012504. ISSN 1098-0121. S2CID 56044826.
  22. ^ Chiang, T. -C (September 1, 2000). "Photoemission studies of quantum well states in thin films". Surface Science Reports. 39 (7): 181–235. Bibcode:2000SurSR..39..181C. doi:10.1016/S0167-5729(00)00006-6. ISSN 0167-5729. (subscription required)
  23. ^ Zhou, S. Y.; Gweon, G.-H.; Graf, J.; Fedorov, A. V.; Spataru, C. D.; Diehl, R. D.; Kopelevich, Y.; Lee, D.-H.; Louie, Steven G.; Lanzara, A. (August 27, 2006). "First direct observation of Dirac fermions in graphite". Nature Physics. 2 (9): 595–599. arXiv:cond-mat/0608069. Bibcode:2006NatPh...2..595Z. doi:10.1038/nphys393. ISSN 1745-2473. S2CID 119505122.
  24. ^ Hsieh, D.; Qian, D.; Wray, L.; Xia, Y.; Hor, Y. S.; Cava, R. J.; Hasan, M. Z. (April 24, 2008). "A topological Dirac insulator in a quantum spin Hall phase : Experimental observation of first strong topological insulator". Nature. 452 (7190): 970–974. arXiv:0902.1356. doi:10.1038/nature06843. ISSN 0028-0836. PMID 18432240. S2CID 4402113.
  25. ^ Liu, Z. K.; Zhou, B.; Wang, Z. J.; Weng, H. M.; Prabhakaran, D.; Mo, S.-K.; Zhang, Y.; Shen, Z. X.; Fang, Z.; Dai, X.; Hussain, Z. (February 21, 2014). "Discovery of a Three-dimensional Topological Dirac Semimetal, Na3Bi". Science. 343 (6173): 864–867. arXiv:1310.0391. Bibcode:2014Sci...343..864L. doi:10.1126/science.1245085. ISSN 0036-8075. PMID 24436183. S2CID 206552029.
  26. ^ Kordyuk, A. A. (May 2, 2014). "ARPES experiment in fermiology of quasi-2D metals (Review Article)". Low Temperature Physics. 40 (4): 286–296. arXiv:1406.2948. Bibcode:2014LTP....40..286K. doi:10.1063/1.4871745. ISSN 1063-777X. S2CID 119228462.
  27. ^ Lu, Donghui; Vishik, Inna M.; Yi, Ming; Chen, Yulin; Moore, Rob G.; Shen, Zhi-Xun (January 3, 2012). "Angle-Resolved Photoemission Studies of Quantum Materials". Annual Review of Condensed Matter Physics. 3 (1): 129–167. doi:10.1146/annurev-conmatphys-020911-125027. ISSN 1947-5454. OSTI 1642351. (subscription required)
  28. ^ Weinelt, Martin (November 4, 2002). "Time-resolved two-photon photoemission from metal surfaces". Journal of Physics: Condensed Matter. 14 (43): R1099–R1141. doi:10.1088/0953-8984/14/43/202. ISSN 0953-8984. (subscription required)
  29. ^ Ueba, H.; Gumhalter, B. (January 1, 2007). "Theory of two-photon photoemission spectroscopy of surfaces". Progress in Surface Science. 82 (4–6): 193–223. doi:10.1016/j.progsurf.2007.03.002. (subscription required)