research32 min studyLesson 984 in Physics

Angle-resolved photoemission and spectral functions

Relate photoelectron kinematics to occupied spectral weight without treating intensity as an unfiltered band picture.

On this page

What you’ll understand

  • Convert photon energy, kinetic energy, and emission angle into binding energy and parallel momentum.
  • State the assumptions connecting ARPES intensity with an occupied single-particle spectral function.
  • Separate intrinsic widths from resolution and matrix-element suppression.

Angle-resolved photoemission removes an electron and records the outgoing electron's energy and direction. It is sensitive to the many-body cost of electron removal, not just to eigenvalues of an isolated one-electron Hamiltonian. The resulting image can reveal dispersions and gaps, but brightness also depends on occupation, photon polarization, escape depth, and instrumental resolution. Interpreting a dark region as an absent electronic state is therefore a hypothesis that needs experimental controls, not a consequence of the image alone.

The kinematic map has a surface boundary

Let binding energy be positive below the sample chemical potential. Energy conservation relates it to photon energy, measured kinetic energy, and the calibrated work function. For a clean crystalline surface, parallel momentum conservation includes the incident photon momentum, modulo a surface reciprocal vector. The displayed mapping neglects that photon contribution and assumes nonrelativistic vacuum-electron kinematics, appropriate approximations for the low-energy example below. More generally, initial parallel crystal momentum is the outgoing parallel electron momentum minus incident parallel photon momentum, modulo a surface reciprocal vector. Analyzer energy calibration and the reference Fermi edge should be established before comparing small shifts or gaps across temperature and photon-energy changes.

Momentum normal to the surface is less direct because translation symmetry is broken and the electron has a finite escape depth. Estimating it often introduces a final-state dispersion and inner potential; varying photon energy can help test a three-dimensional assignment. Surface reconstruction, charging, and cleavage-dependent termination may change the electronic structure itself. A surface-sensitive spectrum and a bulk oscillation experiment can therefore disagree for physical reasons even if both instruments and analyses are functioning correctly.

EB=hν−Φ−Ekin,k∥≃2meEkinℏsin⁡θE_B=h\nu-\Phi-E_{\rm kin},\qquad k_\parallel\simeq\frac{\sqrt{2m_eE_{\rm kin}}}{\hbar}\sin\theta
EB, photon energy, work function Phi, and kinetic energy must use consistent energy units. theta is measured from the surface normal. The initial parallel crystal momentum mapping neglects incident parallel photon momentum and uses nonrelativistic vacuum-electron kinematics; surface reciprocal vectors relate zone representatives.↔ Scroll wide equations horizontally.

The spectral function is multiplied before it is blurred

For a resolved band in equilibrium, the sudden approximation and a suitable photoemission model motivate intensity proportional to a matrix-element factor times the Fermi occupation times the single-particle spectral function. Use energy epsilon relative to the chemical potential, so occupied states lie mostly at negative epsilon. This epsilon is an energy, not an angular frequency. The spectral function includes both electron addition and removal, while ordinary direct photoemission predominantly accesses the occupied removal part.

For a canonical fermionic orbital, the complete spectral function integrates to one over all energies. Raw ARPES intensity generally does not: its matrix element varies with momentum and photon configuration, the positive-energy part is suppressed by occupation, and the experiment has a finite window. A polarization selection rule can suppress a band that still exists. Dividing by a small Fermi factor near or above the chemical potential can amplify noise dramatically and does not reconstruct missing information without assumptions.

A(k,ϵ)=−1πIm⁡GR(k,ϵ),∫−∞∞A dϵ=1,Iideal=C∣M∣2f(ϵ)AA(\mathbf k,\epsilon)=-\frac{1}{\pi}\operatorname{Im}G^R(\mathbf k,\epsilon),\qquad \int_{-\infty}^{\infty}A\,d\epsilon=1,\qquad I_{\rm ideal}=C|M|^2f(\epsilon)A
A has units of inverse energy for one normalized canonical orbital. M is the photoemission matrix element, f is the Fermi function, and C contains apparatus factors. Multiband measurements require the corresponding orbital and matrix-element structure.↔ Scroll wide equations horizontally.

Line shapes require an instrument model

A retarded self-energy changes both peak position and width. For an isolated simple pole with nearly constant negative imaginary self-energy, the spectral peak is Lorentzian; its half width is the positive broadening parameter. A lifetime interpretation also requires the appropriate quasiparticle residue and a clear definition of amplitude versus population decay. Momentum-distribution widths can be converted into energy scales only when dispersion, self-energy dependence, and the measurement cut justify that conversion.

Measured spectra convolve the occupied, matrix-weighted signal with energy and momentum resolution and add background. Gaussian widths add in quadrature only for genuinely Gaussian convolutions; a Lorentzian intrinsic peak convolved with a Gaussian becomes a Voigt profile. An unresolved pair can mimic one broad band, and a leading-edge shift can reflect thermal or resolution changes as well as a gap. Varying polarization, photon energy, temperature, and resolution tests these alternatives before a detailed many-body interpretation is assigned.

Imeas=RE,k∗(C∣M∣2fA)+B,Apole=1πΓ(ϵ−ϵ0)2+Γ2I_{\rm meas}=\mathcal R_{E,k}*\bigl(C|M|^2fA\bigr)+B,\qquad A_{\rm pole}=\frac{1}{\pi}\frac{\Gamma}{(\epsilon-\epsilon_0)^2+\Gamma^2}
The normalized resolution kernel acts on the whole ideal signal, not on A alone. B is background; the displayed Lorentzian assumes a unit-weight pole with energy-independent positive half width Gamma.↔ Scroll wide equations horizontally.

Kinematics before spectral interpretation

For a hypothetical spectrum take photon energy 25 eV25\,\mathrm{eV}, work function 4.5 eV4.5\,\mathrm{eV}, kinetic energy 16 eV16\,\mathrm{eV}, and emission angle 30∘30^\circ. Use 2meE/ℏ=0.512E/eV A˚−1\sqrt{2m_eE}/\hbar=0.512\sqrt{E/\mathrm{eV}}\,\text{\AA}^{-1}. Find binding energy and parallel momentum.

  1. Subtract work function and kinetic energy from photon energy: EB=25−4.5−16=4.5 eVE_B=25-4.5-16=4.5\,\mathrm{eV}. Relative to the chemical potential, the electron energy is ϵ=−4.5 eV\epsilon=-4.5\,\mathrm{eV}.

  2. The outgoing wave-vector magnitude is 0.51216=2.048 A˚−10.512\sqrt{16}=2.048\,\text{\AA}^{-1}. This is the vacuum value, not a direct determination of the initial perpendicular crystal momentum.

  3. Project parallel to the surface using sine of thirty degrees: k∥=1.024 A˚−1k_\parallel=1.024\,\text{\AA}^{-1}, before choosing a surface reciprocal-zone representative.

The result

The feature is 4.5 electronvolts below the chemical potential and has parallel momentum 1.024 inverse angstroms under the stated calibration.

Explore the model

Interactive lab

Interactive lab · vectors

An angular projection

For separate fresh runs Reset and set angle and magnitude to (30,2)(30,2) and (90,2)(90,2). Predict y projections 1 and 2. The projection illustrates an angular factor only; the vector magnitude has no electron-energy calibration, and this lab does not model a surface, a spectral function, photoemission matrix elements, or resolution.

An angular projectionPlayback traces a fixed vector and its projections. Angle and magnitude set the geometry; the moving marker does not change them. Current readings: x component: 1.393; y component: 0.975; Magnitude: 1.7. All parameters are adjustable with the labeled controls below.−7.46−3.7303.737.46−3.5−1.7501.753.5xyvvₓ = 1.39
Vectorx / y projectionsTrace marker
x component
1.393
y component
0.975
Magnitude
1.7
t = 0 s

Playback traces a fixed vector and its projections. Angle and magnitude set the geometry; the moving marker does not change them. Parameter changes pause playback and recompute the model at the current elapsed time. Reset restores both parameters and returns time to zero. Playback starts paused. Adjust either parameter, then choose Play or Step.

Look forA measured angle supplies a component, while the magnitude needs a separate energy calibration.

Check your understanding

What can cause a band to become dark in an ARPES image while its spectral states remain present?

Try it yourself

An illustrative intrinsic Gaussian feature has energy standard deviation 5 meV5\,\mathrm{meV} and is measured with independent Gaussian resolution 12 meV12\,\mathrm{meV}. Find the measured standard deviation. Would the same quadrature rule exactly deconvolve a Lorentzian intrinsic line?

Take a moment to make a prediction or work it out on paper. Open a hint when you need one.

Give me a hint

Add variances for two Gaussian convolution kernels, then distinguish a different line-shape family.

Compare with the solution

The measured variance is 52+122=169 meV25^2+12^2=169\,\mathrm{meV}^2, so its standard deviation is 13 meV13\,\mathrm{meV}. A Lorentzian convolved with a Gaussian is a Voigt profile, so subtracting squared Gaussian widths is not an exact Lorentzian deconvolution or an intrinsic lifetime measurement.

Take it with you

  • ARPES kinematics directly constrains parallel momentum and calibrated removal energy.
  • Intensity contains matrix elements, occupation, resolution, and background as well as spectral weight.
  • A width or a missing feature needs an instrument and line-shape model before a lifetime or gap claim.

Read further

Original sources and deeper explanations to follow at your own pace.

Optional progress, saved in this browser when available. Revisit any lesson at any time.

Lesson not marked complete.