Potential barrier and tunneling

Roll a ball at a hill it does not have enough energy to climb, and classical mechanics is unambiguous: the ball rolls back down. A quantum particle facing the same situation — a region where its potential energy V0V_0 exceeds its total energy EE — obeys a stranger rule. Its wavefunction does not stop dead at the wall; it leaks into the forbidden region and decays there, and if the wall is thin enough some of that leaked wave re-emerges as a genuine, freely moving particle on the far side. This is quantum tunneling, and it is not a rare curiosity — it underlies radioactive alpha decay, the scanning tunneling microscope, and the leakage currents that limit how thin a transistor's insulating layer can be. This article builds the order-of-magnitude tunneling law from the time-independent Schrödinger equation, applies it to an electron and a proton hitting the same barrier, and shows why tunneling is essentially an electron phenomenon.

The core insight

A quantum particle with E<V0E<V_0 still has a wavefunction that decays exponentially — not zero — inside a potential barrier, so the transmission probability falls off as T∼e−2κLT\sim e^{-2\kappa L} with κ=2m(V0−E)/ℏ\kappa=\sqrt{2m(V_0-E)}/\hbar, making tunneling significant only for thin barriers and light, fast particles.

The core idea

Classically forbidden does not mean wavefunction-forbidden. Take a particle of energy EE moving in one dimension and meeting a rectangular barrier of height V0V_0 and width LL, with E<V0E<V_0. Outside the barrier the particle is free and ϕ(x)\phi(x) oscillates like e±ikxe^{\pm ikx}, k=2mE/ℏk=\sqrt{2mE}/\hbar. Inside the barrier the time-independent Schrödinger equation reads

−ℏ22md2ϕdx2+V0ϕ=Eϕ⟹d2ϕdx2=κ2ϕ,κ=2m(V0−E)ℏ.-\frac{\hbar^2}{2m}\frac{d^2\phi}{dx^2} + V_0\phi = E\phi \quad\Longrightarrow\quad \frac{d^2\phi}{dx^2} = \kappa^2\phi,\qquad \kappa=\frac{\sqrt{2m(V_0-E)}}{\hbar}.

where: mm is the particle's mass, EE its total energy, V0V_0 the barrier height, ℏ\hbar the reduced Planck constant, and κ\kappa (units of inverse length) sets how fast the wavefunction decays inside the barrier.

Why the solution is a real exponential. Because E<V0E<V_0, the classical kinetic energy E−V0E-V_0 is negative, so the equation flips sign compared with the free-particle case and the solutions become real exponentials, ϕ(x)∼e±κx\phi(x)\sim e^{\pm\kappa x}, instead of oscillations. A classical particle could never occupy this region at all — but ϕ\phi and dϕ/dxd\phi/dx must stay continuous at each wall of the barrier, matching the decaying interior solution smoothly onto the oscillating exterior waves. That continuity requirement is what forces ϕ\phi to be nonzero, and decaying rather than vanishing, all the way across the barrier.

The transmission law. For a barrier wide enough that the decay is substantial (κL≫1\kappa L\gg1), the wavefunction amplitude that survives to the far wall — and hence the probability of the particle being found beyond the barrier — is governed by how much ϕ\phi has decayed:

T∼e−2κL,κ=2m(V0−E)ℏ.T \sim e^{-2\kappa L}, \qquad \kappa=\frac{\sqrt{2m(V_0-E)}}{\hbar}.

where: TT is the tunneling probability (dimensionless, between 00 and 11), LL is the barrier width, and κ\kappa is as defined above. This is the order-of-magnitude estimate expected at this level — not the exact matched-boundary-condition transmission coefficient, which needs more algebra than the syllabus asks for.

The step case, for contrast. A potential step (no second wall) shows the same wave-mechanical stubbornness in a different guise. For E>V0E>V_0 at a step, classical mechanics guarantees transmission, but quantum mechanics predicts partial reflection at the discontinuity in V(x)V(x), because ϕ\phi and ϕ′\phi' must still match continuously across it. For E<V0E<V_0 at a step, ϕ\phi decays exponentially into the forbidden region but, with no far wall to escape through, reflection is eventually total.

Figure
Wavefunction oscillating on both sides of a barrier and decaying exponentially inside it, with the transmitted amplitude on the right smaller than the incident amplitude on the left
Tunneling through a rectangular barrier barrier, width L, height V0 incident + reflected wave decaying real exponential transmitted wave, smaller amplitude The wave never goes to zero inside the barrier, so a smaller wave survives on the far side.
Derivation ATunneling probability across a rectangular barrier

The goal is to justify the order-of-magnitude tunneling law T∼e−2κLT\sim e^{-2\kappa L} starting only from the interior Schrödinger equation and the requirement that the wavefunction and its slope stay continuous at the two walls of the barrier, without carrying out the full multi-region algebra that fixes the exact prefactor in front of the exponential.

Step 1. Write the interior equation. For 0<x<L0<x<L with E<V0E<V_0, the time-independent Schrödinger equation is −ℏ22mϕ′′=(E−V0)ϕ-\frac{\hbar^2}{2m}\phi''=(E-V_0)\phi. Since E−V0<0E-V_0<0, rearrange as ϕ′′=κ2ϕ\phi''=\kappa^2\phi with κ=2m(V0−E)/ℏ\kappa=\sqrt{2m(V_0-E)}/\hbar real and positive. The general interior solution is ϕ(x)=Ae−κx+Beκx\phi(x)=Ae^{-\kappa x}+Be^{\kappa x}, a mixture of a decaying and a growing real exponential — no oscillation is possible because the effective "kinetic energy" E−V0E-V_0 is negative here.

Step 2. Match at the entrance wall, x=0x=0. Outside, for x<0x<0, the particle is free and its wavefunction is a mix of an incident wave eikxe^{ikx} and a reflected wave, with k=2mE/ℏk=\sqrt{2mE}/\hbar. Continuity of ϕ\phi and ϕ′\phi' at x=0x=0 fixes AA and BB in terms of the incident amplitude; the key point is that neither AA nor BB can be zero in general, so ϕ\phi is nonzero everywhere inside the barrier, not just near the edges.

Step 3. Track the amplitude to the exit wall, x=Lx=L. The decaying piece Ae−κxAe^{-\kappa x} shrinks by a factor e−κLe^{-\kappa L} across the barrier while the growing piece BeκxBe^{\kappa x} grows by eκLe^{\kappa L}; continuity again at x=Lx=L hands this combination over to an outgoing free wave eikxe^{ikx} for x>Lx>L. When the barrier is wide (κL≫1\kappa L\gg1), the dominant surviving contribution to the transmitted amplitude scales as e−κLe^{-\kappa L}, because the growing solution's coefficient is itself suppressed by matching conditions that keep the total wave bounded.

Step 4. Convert amplitude to probability. Transmission probability is proportional to the square of the transmitted wave's amplitude relative to the incident one, so a transmitted amplitude scaling as e−κLe^{-\kappa L} gives T∼(e−κL)2=e−2κLT\sim\left(e^{-\kappa L}\right)^2=e^{-2\kappa L}. This is the working formula: TT depends exponentially on both the barrier width LL and on V0−E\sqrt{V_0-E} through κ\kappa, so a barrier twice as wide does not halve TT — it roughly squares the already-small value, since e−2κ(2L)=(e−2κL)2e^{-2\kappa(2L)}=\left(e^{-2\kappa L}\right)^2.

Worked example

An electron of kinetic energy E=1.0 eVE=1.0\ \text{eV} approaches a rectangular barrier of height V0=2.0 eVV_0=2.0\ \text{eV} and width L=0.50 nmL=0.50\ \text{nm}. Find the tunneling probability, then compare it with a proton facing the same energy deficit and width.

Step 1. Convert the energy deficit to SI: V0−E=1.0 eV=1.0×1.6×10−19 J=1.6×10−19 JV_0-E=1.0\ \text{eV}=1.0\times1.6\times10^{-19}\ \text{J}=1.6\times10^{-19}\ \text{J}.

Step 2. Compute κ\kappa for the electron, using me=9.11×10−31 kgm_e=9.11\times10^{-31}\ \text{kg} and ℏ=1.055×10−34 J ⋅ s\hbar=1.055\times10^{-34}\ \text{J}\,\cdot\,\text{s}: the product 2me(V0−E)=2×9.11×10−31×1.6×10−19=2.92×10−492m_e(V_0-E)=2\times9.11\times10^{-31}\times1.6\times10^{-19}=2.92\times10^{-49}, whose square root is 5.40×10−255.40\times10^{-25}. Dividing by ℏ\hbar gives κ=5.40×10−25/1.055×10−34=5.12×109 m−1\kappa=5.40\times10^{-25}/1.055\times10^{-34}=5.12\times10^{9}\ \text{m}^{-1}.

Step 3. Compute 2κL=2×5.12×109×0.50×10−9=5.122\kappa L=2\times5.12\times10^{9}\times0.50\times10^{-9}=5.12, so T∼e−5.12≈6.0×10−3T\sim e^{-5.12}\approx6.0\times10^{-3} — about 0.6%0.6\% of electrons tunnel through, small but very much observable at nanometre scales.

Step 4. Scale up to a proton of the same energy deficit and width. Since κ∝m\kappa\propto\sqrt{m} and mp/me≈1836m_p/m_e\approx1836, κp=κe1836≈5.12×109×42.8≈2.19×1011 m−1\kappa_p=\kappa_e\sqrt{1836}\approx5.12\times10^{9}\times42.8\approx2.19\times10^{11}\ \text{m}^{-1}, giving 2κpL≈2192\kappa_pL\approx219.

Step 5. Estimate Tp∼e−219≈10−95T_p\sim e^{-219}\approx10^{-95} — utterly negligible. The electron tunnels with probability roughly 6×10−36\times10^{-3}, while a proton under identical conditions tunnels with probability of order 10−9510^{-95}: because κ∝m\kappa\propto\sqrt{m}, tunneling is dramatically suppressed for heavier particles, which is why tunneling phenomena in electronics and atomic physics are almost always electron effects.

Common traps

  • Using κ=2mE/ℏ\kappa=\sqrt{2mE}/\hbar instead of 2m(V0−E)/ℏ\sqrt{2m(V_0-E)}/\hbar. Inside the barrier the relevant quantity is the energy deficit V0−EV_0-E, not the particle's own energy EE.
  • Thinking the particle "loses energy" while tunneling. It does not — the particle that emerges beyond the barrier still has kinetic energy EE, unchanged. Tunneling is a probability-of-transmission effect, not an energy-loss effect.
  • Assuming "enough energy always means getting through." At a potential step with E>V0E>V_0, classical mechanics guarantees transmission, but quantum mechanics still gives partial reflection because ϕ\phi and ϕ′\phi' must match continuously across the discontinuity.
  • Treating TT as roughly linear in LL. Because TT is exponential in LL, doubling the barrier width does not halve TT — it approximately squares the already-small probability.

Ask the tutor: If you want to see the full matched-boundary-condition transmission coefficient behind this order-of-magnitude formula, or want a fresh barrier-tunneling number checked step by step, ask the AI tutor to walk through it with you.

Why does kappa depend on V0 minus E and not on V0 or E separately?

Inside the barrier the Schrodinger equation becomes ϕ′′=κ2ϕ\phi''=\kappa^2\phi with κ=2m(V0−E)/ℏ\kappa=\sqrt{2m(V_0-E)}/\hbar, and this combination arises directly from rearranging −ℏ22mϕ′′=(E−V0)ϕ-\frac{\hbar^2}{2m}\phi''=(E-V_0)\phi. The physically meaningful quantity is how far short of the barrier's height the particle's energy falls — the energy deficit — not either EE or V0V_0 on its own. A particle with E=1.0 eVE=1.0\ \text{eV} facing V0=2.0 eVV_0=2.0\ \text{eV} behaves identically, in terms of κ\kappa, to one with E=5.0 eVE=5.0\ \text{eV} facing V0=6.0 eVV_0=6.0\ \text{eV}, because both have the same deficit of 1.0 eV1.0\ \text{eV}.

Does the particle exist inside the barrier the way it exists outside?

Its wavefunction is nonzero and real (a decaying exponential) inside the barrier, so there is a nonzero probability density there, but the particle cannot be found there classically because its kinetic energy would have to be negative. This is exactly the tension quantum mechanics resolves: ϕ\phi can be nonzero in a classically forbidden region as long as it decays smoothly and matches continuity conditions at the boundaries. Nothing about this violates energy conservation — the particle that eventually emerges beyond the barrier still carries its original energy EE, unchanged.

Why is tunneling essentially never seen for protons or heavier particles at nanometre scales?

Because κ=2m(V0−E)/ℏ\kappa=\sqrt{2m(V_0-E)}/\hbar grows with the square root of the mass, and T∼e−2κLT\sim e^{-2\kappa L} is exponentially sensitive to κ\kappa. Even though a proton is only about 18361836 times heavier than an electron, κ\kappa scales as 1836≈42.8\sqrt{1836}\approx42.8, which multiplies the exponent by that same factor. For the worked example's barrier this turns a modest exponent of 5.125.12 for the electron into an exponent of about 219219 for the proton, crushing TT from roughly 0.6%0.6\% down to about 10−9510^{-95} — effectively zero.

Before this

Where this leads