|q⟩ Bad Qubits

beginner · Physics · Quantum Tunneling & Barriers

Exponential Decay in the Barrier

The classically-forbidden region

Imagine sending a classical particle toward a hill whose height in energy units exceeds the particle's total energy. The particle cannot cross — it bounces back. Quantum mechanics tells a different story. The time-independent Schrödinger equation,

−ℏ22md2ψdx2+V(x)ψ=Eψ,-\frac{\hbar^2}{2m}\frac{d^2\psi}{dx^2} + V(x)\psi = E\psi,

must still hold inside the barrier. Rearranging for the second derivative gives

d2ψdx2=2m(V0−E)ℏ2 ψ.\frac{d^2\psi}{dx^2} = \frac{2m(V_0 - E)}{\hbar^2}\,\psi.

When E<V0E < V_0 the coefficient on the right is positive. A positive coefficient means the solution is not oscillatory (as it is outside the barrier) but exponential. The two independent solutions are e+κxe^{+\kappa x} and e−κxe^{-\kappa x}, where we define the decay constant

κ=2m(V0−E)ℏ.\kappa = \frac{\sqrt{2m(V_0 - E)}}{\hbar}.

For a particle coming from the left, the growing exponential e+κxe^{+\kappa x} must be discarded because it would diverge — only the decaying solution e−κxe^{-\kappa x} survives inside a thick barrier. The wavefunction therefore falls off as

ψbarrier(x)∝e−κx,\psi_{\text{barrier}}(x) \propto e^{-\kappa x},

and the probability density decays as e−2κxe^{-2\kappa x}. The word "exponential" in the lesson title refers exactly to this spatial decay.

What κ\kappa tells you

The decay constant κ\kappa has units of inverse length. Its reciprocal 1/κ1/\kappa is the decay length — the distance over which the amplitude falls by a factor of e≈2.718e \approx 2.718.

For an electron with V0−E=1.0 eVV_0 - E = 1.0\,\text{eV}, inserting me=9.109×10−31 kgm_e = 9.109 \times 10^{-31}\,\text{kg} and ℏ=1.055×10−34 J⋅s\hbar = 1.055 \times 10^{-34}\,\text{J}\cdot\text{s} yields

κ=2×(9.109×10−31 kg)×(1.602×10−19 J)1.055×10−34 J⋅s≈5.12×109 m−1,\kappa = \frac{\sqrt{2 \times (9.109 \times 10^{-31}\,\text{kg}) \times (1.602 \times 10^{-19}\,\text{J})}}{1.055 \times 10^{-34}\,\text{J}\cdot\text{s}} \approx 5.12 \times 10^9\,\text{m}^{-1},

so the decay length is 1/κ≈0.20 nm1/\kappa \approx 0.20\,\text{nm} — roughly the diameter of two hydrogen atoms. Tunneling currents therefore vanish on the scale of a few atoms, which is precisely why scanning tunneling microscopes achieve atomic resolution.

Connecting κ\kappa to the transmission probability

For a rectangular barrier of width LL and when κL≫1\kappa L \gg 1, the transmission probability simplifies to the WKB approximation

T≈e−2κL.T \approx e^{-2\kappa L}.

This factor of 2κ2\kappa in the exponent explains why the tunneling probability is extraordinarily sensitive to barrier parameters — doubling LL squares TT. The exact expression from matching boundary conditions is more involved, but in the thick-barrier limit the exponential dominates and κ\kappa is the one number that governs tunneling.

Try it

This is a numerical exercise — your code should return a number, not a circuit. Compute κ\kappa in m−1\text{m}^{-1} for an electron facing a barrier with V0−E=1.0 eVV_0 - E = 1.0\,\text{eV}.

Run your code to see the quantum state.

Sign in on the full site to ask questions and join the discussion.