research32 min studyLesson 985 in Physics

Scanning tunneling spectroscopy and local density of states

Interpret a local tunneling spectrum while tracking the tip, barrier, feedback, and thermal kernel.

On this page

What you’ll understand

  • Declare a bias convention and derive the low-temperature conductance approximation.
  • Distinguish a constant-current topograph from a map of atomic height.
  • Identify matrix-element, setpoint, and broadening ambiguities in local spectra and interference maps.

A scanning tunneling microscope places a conducting tip close enough to a surface for electrons to tunnel across the vacuum barrier. Moving the tip provides exceptional spatial discrimination, but the measured current is an electronic overlap rather than a direct ruler for atomic positions. Spectroscopy varies bias while controlling the tip height and asks how accessible electronic states change with energy. Its local character complements momentum-resolved photoemission, yet both techniques require assumptions about matrix elements and the prepared surface.

From a tunneling current to an energy window

Choose positive sample bias V to access unoccupied sample states, and define positive tunneling current for that electron-transfer direction. Measure sample energy E from its chemical potential, with the tip chemical potential higher by eV. In weak elastic tunneling, a golden-rule expression combines sample and tip densities of states, the squared tunneling matrix element, and the difference of their occupation factors. Energy conservation selects pairs of states at the same absolute energy, not simply at the same energy measured from their separate chemical potentials.

With a featureless tip and nearly constant matrix element over the relevant energy window, the zero-temperature positive-bias current integrates the sample's local density of states from zero to eV. Differentiation then selects its upper boundary. This useful limit is often summarized as conductance measuring local density of states, but the local quantity is evaluated through the tip wave function and barrier. Tip orbital structure, a superconducting or magnetic tip, inelastic channels, and bias-dependent transmission change the relation.

I(V)∝∫dE ρs(r,E) [f(E−eV)−f(E)],dIdV∝e∫dE ρs(r,E)[−f′(E−eV)]I(V)\propto\int dE\,\rho_s(\mathbf r,E)\,[f(E-eV)-f(E)],\qquad \frac{dI}{dV}\propto e\int dE\,\rho_s(\mathbf r,E)[-f'(E-eV)]
This reduced expression assumes a constant tip density of states and matrix element, weak elastic tunneling, and fixed tip position. rho_s has units of states per energy per chosen local normalization; proportionality factors include the junction coupling.↔ Scroll wide equations horizontally.

Feedback mixes topography with electronic structure

In constant-current imaging a feedback loop changes tip height to hold an integrated tunneling current at its setpoint. More accessible sample spectral weight can therefore make the tip retract even when the atomic positions are unchanged. The approximate exponential dependence on tip distance amplifies small geometric changes, but it also makes the current sensitive to the barrier and tip condition. A bright topographic feature can combine true height, energy-dependent electronic structure, and local work-function changes.

For spectroscopy the feedback is commonly stabilized at a chosen bias and current, then disabled while a bias trace is recorded at approximately fixed height. Different stabilization conditions can give different multiplicative spectral weights. Comparing maps therefore requires tracking the setpoint and testing reproducibility with several tips and heights. Normalized conductance can reduce some matrix-element variations, but division by current or voltage can introduce its own singular or misleading behavior near zero current and does not universally yield an absolute density of states.

Iset≃Ce−2κz∫0eVsetρs(r,E) dE,κ≃2meϕbℏI_{\rm set}\simeq C e^{-2\kappa z}\int_0^{eV_{\rm set}}\rho_s(\mathbf r,E)\,dE,\qquad \kappa\simeq\frac{\sqrt{2m_e\phi_b}}{\hbar}
z is a gap distance, kappa is an inverse length, and phi_b is an approximate positive barrier energy. This rectangular-barrier scaling is illustrative; real barriers and tip wave functions need not have one spatially uniform kappa.↔ Scroll wide equations horizontally.

A local gap and a scattering pattern need interpretation

At finite temperature the derivative of the Fermi function broadens a sharp feature. Lock-in modulation, voltage noise, instrumental response, and intrinsic lifetime broadening add further kernels. Their widths cannot always be subtracted by a Gaussian rule. A suppressed zero-bias conductance can indicate a gap in accessible local spectral weight, but a tunneling matrix-element minimum, junction instability, or charging effect can also suppress current. Spatial reproducibility and temperature or field dependence help distinguish these mechanisms.

Fourier transforming conductance modulations can reveal quasiparticle interference: defects connect states by scattering wave vectors. The map depends on the impurity potential, orbital and coherence factors, and the energy-dependent electronic propagators. It is not generally a direct image of the Fermi surface or a joint density of states without weighting. A missing scattering vector need not mean the corresponding states are absent. Spatially resolved spectra should be compared with diffraction and bulk measurements before surface inhomogeneity is assigned to an entire sample.

−f′(E)=14kBTsech⁡2 ⁣(E2kBT),ΔEFWHM=4arcosh⁡(2)kBT≃3.53kBT-f'(E)=\frac{1}{4k_BT}\operatorname{sech}^2\!\left(\frac{E}{2k_BT}\right),\qquad \Delta E_{\rm FWHM}=4\operatorname{arcosh}(\sqrt2)k_BT\simeq3.53k_BT
The FWHM belongs to the thermal derivative kernel alone. It is an energy width, not a universal total experimental resolution or a direct intrinsic gap threshold.↔ Scroll wide equations horizontally.

An electronic feature can appear as a height

In a hypothetical constant-current image the integrated accessible local spectral weight doubles at a patch, with unchanged atomic height and barrier. Take κ=1 A˚−1\kappa=1\,\text{\AA}^{-1}. Find the feedback height change needed to preserve the same current.

  1. Write the current ratio after a height change as Inew/Iold=2e−2κΔzI_{\rm new}/I_{\rm old}=2e^{-2\kappa\Delta z}. The factor two is an electronic-weight change, not a geometric displacement.

  2. Holding the ratio at one gives 2κΔz=ln⁡22\kappa\Delta z=\ln2. The positive solution means that the tip retracts rather than approaches the patch.

  3. With the stated inverse-length scale, Δz=(ln⁡2)/2 A˚≃0.347 A˚\Delta z=(\ln2)/2\,\text{\AA}\simeq0.347\,\text{\AA}. The same atoms can thus produce an apparent topographic height in this simplified junction model.

The result

The feedback retracts by about 0.347 angstrom even though the atomic height is unchanged.

Explore the model

Interactive lab

Interactive lab · function

A nonlinear signal need not be its derivative

For fresh paused runs Reset and set coefficient a and position x to (0.5,1)(0.5,1) and (0.5,2)(0.5,2). Predict y values 0.5 and 2 but tangent slopes 1 and 2. This distinguishes a signal from its derivative; the quadratic curve is not a tunneling spectrum, a feedback simulation, or a density-of-states reconstruction.

A nonlinear signal need not be its derivativeThis lab uses only y = ax². Playback traces the curve within the fixed integration interval from 0 to x; it does not estimate numerical error or change the chosen point. Current readings: y at chosen x: 0.608; Tangent slope · 2ax: 1.35; Signed area · ax³/3: 0.182. All parameters are adjustable with the labeled controls below.−2.4−1.201.22.4−9−4.504.59xy(x, ax²)
y = ax²TangentSigned integral
y at chosen x
0.608
Tangent slope · 2ax
1.35
Signed area · ax³/3
0.182
t = 0 s

This lab uses only y = ax². Playback traces the curve within the fixed integration interval from 0 to x; it does not estimate numerical error or change the chosen point. 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 current-like quantity and its differential slope carry different information.

Check your understanding

Why can a constant-current STM image show apparent height where atomic positions are unchanged?

Try it yourself

For a hypothetical low-temperature spectrum let kBT=0.20 meVk_BT=0.20\,\mathrm{meV}. Estimate the FWHM of the thermal conductance kernel. If the voltage modulation adds further broadening, is this number the full experimental energy resolution?

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

Give me a hint

Use the coefficient 3.53 for the derivative-of-Fermi kernel, keeping energy and voltage widths distinct.

Compare with the solution

The thermal energy width is about 3.53(0.20)=0.706 meV3.53(0.20)=0.706\,\mathrm{meV}. It corresponds to 0.706 millivolt for a singly charged electron, but it is only the thermal kernel's width. Modulation, noise, and other instrumental kernels further shape the measured spectrum, so 0.706 meV is not the complete resolution.

Take it with you

  • Conductance approximates local spectral weight only with a specified tip, barrier, and occupation model.
  • Constant-current topography mixes electronic weight with distance.
  • Thermal broadening and scattering selection rules limit gap and interference inferences.

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.