Angleresolved photoemission spectroscopy
Angleresolved 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 twodimensional 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 hightemperature superconductors and materials exhibiting charge density waves.
The main components of an ARPES system are a source to deliver a highfrequency 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 ultrahigh 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 twodimensional maps, the sample is rotated, or the electrons are manipulated.
Instrumentation
A typical instrument for angleresolved photoemission consists of a light source, a sample holder attached to a manipulator, and an electron spectrometer. These are all part of an ultrahigh 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, highintensity photon beam to the sample (~10^{12} photons/s with a few meV energy spread).[2] Light sources range from compact noblegas discharge UV lamps and radiofrequency 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 Xrays).
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 electronbeam bombardment for temperatures up to 2000 °C, and light beam focusing and calibration.
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.
The energy dispersion is carried out for a narrow range of energies around the socalled 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 passenergy and slitwidth 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, timeofflight analyzers are used. These, however, require pulsed photon sources and are most common in laserbased ARPES labs.[8]
Theory
Principle
Angleresolved 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 lowdimensional systems such as twodimensional materials, ultrathin films, and nanowires. When it is used for threedimensional materials, the perpendicular component of the wave vector is usually approximated, with the assumption of a parabolic, freeelectronlike final state with the bottom at energy . This gives:
Fermi surface mapping
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 twodimensional 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 energyangleangle maps can be further processed to give energyk_{x}k_{y} 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
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 twodimensional 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, k_{x}, k_{y}) around origin in twodimensional momentum planes.
 Analyzing electrons emitted into polar angle ϑ and α≤8° around.
 Polar angle ϑ, tilt τ, and α≤8°.
 Polar angle ϑ, tilt τ, and α≤8°. Azimuth is set to φ.
Theoretical derivation of intensity relationship
The theory of photoemission[1][10][12] is that of direct optical transitions between the states and of an Nelectron system. Light excitation is introduced as the magnetic vector potential through the minimal substitution in the kinetic part of the quantummechanical 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 timedependent 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 N1 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 oneelectron removal spectral function known from the manybody 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 E_{F} and E_{k} 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 FermiDirac distribution function in the form of a temperaturedependent sigmoidshaped drop of intensity in the vicinity of E_{F}. In the case of a twodimensional, oneband electronic system the intensity relation further reduces to .[1]
Selection rules
The electronic states in crystals are organized in energy bands, which have associated energyband dispersions that are energy eigenvalues for delocalized electrons according to Bloch's theorem. From the planewave 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 oneelectron 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 atomicorbital origin of the initial and final states.[14]
Manybody effects
The oneelectron 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 energymomentum space. In the case of increased electron correlations, the spectral function broadens and starts developing richer features that reflect the interactions in the underlying manybody 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 highresolution ARPES spectra under a few reasonable assumptions. Namely, one can assume that the part of the spectrum is nearly constant along highsymmetry 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
This function is known from ARPES as a scan along a chosen direction in momentum space and is a twodimensional map of the form . When cut at a constant energy , a Lorentzianlike 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 selfconsistent way by enforcing the KramersKronig 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 KramersKronig 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 selfenergy, one can judge on the strength and shape of electronelectron correlations, electronphonon (more generally, electronboson) 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 electronphonon coupling factor λ can be determined from the linear dependence of the peak widths on temperature.[20]
Uses
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] oneatomthin 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 hightemperature 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, twophoton excitation in pumpprobe setups (2PPE) is used. There, the first photon of lowenough 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 lowenergy 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]
References
 Damascelli, Andrea; Shen, ZhiXun; Hussain, Zahid (April 17, 2003). "Angleresolved photoemission spectroscopy of the cuprate superconductors". Reviews of Modern Physics. 75 (2): 473–541. arXiv:condmat/0208504. doi:10.1103/RevModPhys.75.473. ISSN 00346861. S2CID 118433150.
 Hüfner, Stefan, ed. (2007). Very High Resolution Photoelectron Spectroscopy. Lecture Notes in Physics. 715. Berlin, Heidelberg: Springer Berlin Heidelberg. doi:10.1007/3540681337. ISBN 9783540681304. (subscription required)
 "MBScientific electron analysers and UV sources".
 "ARPES Lab". Scienta Omicron. 2020. Retrieved August 29, 2020.
 "Lab ARPES System with PHOIBOS Analyzer". SPECS. Retrieved August 29, 2020.
 "Products". Lumeras LLC. 2013. Retrieved August 29, 2020.
 "Light sources of the world".
 Zhou, Xingjiang; He, Shaolong; Liu, Guodong; Zhao, Lin; Yu, Li; Zhang, Wentao (June 1, 2018). "New Developments in LaserBased 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/13616633/aab0cc. ISSN 00344885. PMID 29460857. S2CID 3440746.
 Soper, Davison E. "Electromagnetic radiation is made of photons". Retrieved September 3, 2020.
 Hüfner, Stefan. (2003). "Introduction and Basic Principles". Photoelectron Spectroscopy: Principles and Applications (Third rev. and enlarged ed.). Berlin, Heidelberg: Springer Berlin Heidelberg. ISBN 9783662092804. OCLC 851391282.
 Damascelli, Andrea; Shen, ZhiXun; Hussain, Zahid (April 17, 2003). "Angleresolved photoemission spectroscopy of the cuprate superconductors". Reviews of Modern Physics. 75 (2): 473–541. arXiv:condmat/0208504. doi:10.1103/RevModPhys.75.473. ISSN 00346861. S2CID 118433150.
 Damascelli, Andrea (2004). "Probing the LowEnergy Electronic Structure of Complex Systems by ARPES". Physica Scripta. T109: 61. arXiv:condmat/0307085. doi:10.1238/Physica.Topical.109a00061. ISSN 00318949. S2CID 21730523.
 Wacker, Andreas. "Fermi's golden rule" (PDF). Teaching notes (Lund University).
 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). "InPlane 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 17452473.
 Pletikosić, Ivo; Kralj, Marko; Milun, Milorad; Pervan, Petar (April 24, 2012). "Finding the bare band: Electron coupling to two phonon modes in potassiumdoped graphene on Ir(111)". Physical Review B. 85 (15): 155447. arXiv:1201.0777. Bibcode:2012PhRvB..85o5447P. doi:10.1103/PhysRevB.85.155447. ISSN 10980121. S2CID 119170154.
 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:condmat/0405696. doi:10.1103/PhysRevB.71.214513. ISSN 10980121. S2CID 67784336.
 Norman, M. R.; Ding, H.; Fretwell, H.; Randeria, M.; Campuzano, J. C. (September 1, 1999). "Extraction of the Electron SelfEnergy from Angle Resolved Photoemission Data: Application to Bi2212". Physical Review B. 60 (10): 7585–7590. arXiv:condmat/9806262. doi:10.1103/PhysRevB.60.7585. ISSN 01631829. S2CID 4691468.
 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 01631829. (subscription required)
 Valla, T.; Fedorov, A. V.; Johnson, P. D.; Hulbert, S. L. (September 6, 1999). "ManyBody Effects in AngleResolved Photoemission: Quasiparticle Energy and Lifetime of a Mo(110) Surface State". Physical Review Letters. 83 (10): 2085–2088. arXiv:condmat/9904449. Bibcode:1999PhRvL..83.2085V. doi:10.1103/PhysRevLett.83.2085. ISSN 00319007. S2CID 55072153.
 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/13672630/11/12/125005. ISSN 13672630.
 Veenstra, C. N.; Goodvin, G. L.; Berciu, M.; Damascelli, A. (July 16, 2010). "Elusive electronphonon 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 10980121. S2CID 56044826.
 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/S01675729(00)000066. ISSN 01675729. (subscription required)
 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:condmat/0608069. Bibcode:2006NatPh...2..595Z. doi:10.1038/nphys393. ISSN 17452473. S2CID 119505122.
 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 00280836. PMID 18432240. S2CID 4402113.
 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 Threedimensional Topological Dirac Semimetal, Na3Bi". Science. 343 (6173): 864–867. arXiv:1310.0391. Bibcode:2014Sci...343..864L. doi:10.1126/science.1245085. ISSN 00368075. PMID 24436183. S2CID 206552029.
 Kordyuk, A. A. (May 2, 2014). "ARPES experiment in fermiology of quasi2D metals (Review Article)". Low Temperature Physics. 40 (4): 286–296. arXiv:1406.2948. Bibcode:2014LTP....40..286K. doi:10.1063/1.4871745. ISSN 1063777X. S2CID 119228462.
 Lu, Donghui; Vishik, Inna M.; Yi, Ming; Chen, Yulin; Moore, Rob G.; Shen, ZhiXun (January 3, 2012). "AngleResolved Photoemission Studies of Quantum Materials". Annual Review of Condensed Matter Physics. 3 (1): 129–167. doi:10.1146/annurevconmatphys020911125027. ISSN 19475454. OSTI 1642351. (subscription required)
 Weinelt, Martin (November 4, 2002). "Timeresolved twophoton photoemission from metal surfaces". Journal of Physics: Condensed Matter. 14 (43): R1099–R1141. doi:10.1088/09538984/14/43/202. ISSN 09538984. (subscription required)
 Ueba, H.; Gumhalter, B. (January 1, 2007). "Theory of twophoton photoemission spectroscopy of surfaces". Progress in Surface Science. 82 (4–6): 193–223. doi:10.1016/j.progsurf.2007.03.002. (subscription required)