Skip to content

Theory of Operation

Published: at 12:39 AM (2 min read)

Tunneling Current Derivation

The scanning tunneling microscope (STM) relies on quantum tunneling of electrons across the vacuum gap between a sharp metal tip and a conductive sample. In a one-dimensional (1D) model, an electron with energy EE approaching a potential barrier of height U0> ⁣EU_0>\!E over a width dd satisfies the time-independent Schrödinger equation 22md2ψdz2+U(z)ψ=Eψ-\frac{\hbar^2}{2m}\frac{d^2\psi}{dz^2} +U(z)\psi=E\psi. Using the WKB approximation, the decaying wavefunction under the barrier implies a transmission probability

T    exp ⁣(20dκ(z)dz),κ(z)=2m[U(z)E]2.T \;\approx\; \exp\!\Bigl(-2\int_0^d \kappa(z)\,dz\Bigr), \quad \kappa(z)=\sqrt{\frac{2m[U(z)-E]}{\hbar^2}}.

For a rectangular vacuum barrier of height U0=ϕU_0=\phi (on the order of the work function) and assuming EEFϕE\approx E_F\ll\phi, κ\kappa is approximately constant. Thus Texp(2κd)T\approx \exp(-2\kappa d) with κ=2mϕ/\kappa=\sqrt{2m\phi}/\hbar. The resulting tunneling current II is proportional to the transmission probability times the number of available electrons. To leading order at low temperature and small bias VV, Bardeen’s formalism (as simplified by Tersoff–Hamann) yields

I(V)    e2κd0eVρs(E)dE,I(V)\;\propto\; e^{-2\kappa d}\int_0^{eV}\rho_s(E)\,dE,

where ρs(E)\rho_s(E) is the sample’s electronic density of states (DOS) and the tip DOS is assumed constant. In practice one often writes II0e2κdVI\approx I_0\,e^{-2\kappa d}V, showing the exponential decay with distance. Here κ=2mϕˉ/\kappa=\sqrt{2m\bar\phi}/\hbar uses the average work function ϕˉ\bar\phi of tip and sample. This formalism shows that the tunneling current depends both on tip-sample separation and on the integral of the sample’s DOS from the Fermi level to eVeV.

Exponential Dependence on Tip–Sample Separation

From the above derivation, I(d)e2κdI(d)\propto e^{-2\kappa d}, so the tunneling current decays exponentially with increasing gap dd. Equivalently,

I(d)  =  I(d0)exp[2κ(dd0)].I(d)\;=\;I(d_0)\,\exp\bigl[-2\kappa\,(d-d_0)\bigr].

Typical metal work functions ϕ4\phi\sim466\,eV give κ2mϕ/1010m1\kappa\approx\sqrt{2m\phi}/\hbar\sim10^{10}\text{m}^{-1} (i.e.\ 10 nm1\sim10\text{ nm}^{-1}). Thus increasing dd by 0.10.1\,nm (1) multiplies II by exp(2κ0.1nm)\exp(-2\kappa\cdot0.1\,{\rm nm}), on the order of e20.14e^{-2}\sim0.14 or roughly a factor of ten reduction. More precisely, experiments find roughly an order-of-magnitude change per 0.1 nm increment of gap. This extreme sensitivity underpins the STM’s ability to resolve atomic-scale height variations: a small 0.01\sim0.01\,nm change in dd can alter II by a significant percentage. For analysis, one often linearizes by writing lnI2κd+const\ln I\approx -2\kappa d + \text{const}, so that lnI/d=2κ\partial \ln I/\partial d = -2\kappa. This shows that in constant-current mode, the tip height dd responds logarithmically to changes in surface height.