Skip to content

WKB Barrier Tunneling Worked Example

This worked example estimates tunneling through a one-dimensional rectangular barrier using the WKB exponent, then compares the result with the exact opaque-barrier limit.

Use

V(x)={0,x<0,V0,0<x<a,0,x>a,V(x) = \begin{cases} 0, & x<0,\\ V_0, & 0<x<a,\\ 0, & x>a, \end{cases}

with

0<E<V0.0<E<V_0.

The exact rectangular barrier is not smooth at the edges, so it is not the ideal setting for WKB connection formulas. It is still a useful benchmark because the under-barrier exponent is simple and the exact result is known.

The WKB estimate should reproduce the leading exponential dependence when the barrier is opaque:

αa≫1,\alpha a\gg1,

where

α=2m(V0−E)ℏ\alpha = \frac{\sqrt{2m(V_0-E)}}{\hbar}

is the decay constant inside the barrier.

The forbidden region is 0<x<a0<x<a. The WKB barrier action in momentum units is

K(E)=∫0a2m(V0−E) dx.K(E) = \int_0^a \sqrt{2m(V_0-E)} \,dx.

Since the integrand is constant,

K(E)=a2m(V0−E)=ℏαa.K(E) = a\sqrt{2m(V_0-E)} = \hbar\alpha a.

The leading WKB transmission probability is

TWKB≈exp⁡[−2K(E)ℏ],T_{\mathrm{WKB}} \approx \exp \left[ - \frac{2K(E)}{\hbar} \right],

so

TWKB≈e−2αa.T_{\mathrm{WKB}} \approx e^{-2\alpha a}.

The exact transmission coefficient for the rectangular barrier is

Texact=[1+V02sinh⁡2(αa)4E(V0−E)]−1.T_{\mathrm{exact}} = \left[ 1+ \frac{V_0^2\sinh^2(\alpha a)} {4E(V_0-E)} \right]^{-1}.

In the opaque limit,

sinh⁡2(αa)≈14e2αa.\sinh^2(\alpha a) \approx \frac14e^{2\alpha a}.

Therefore

Texact≈16E(V0−E)V02e−2αa.T_{\mathrm{exact}} \approx \frac{16E(V_0-E)}{V_0^2} e^{-2\alpha a}.

The WKB estimate captures the leading exponential:

Texact∝TWKBT_{\mathrm{exact}} \propto T_{\mathrm{WKB}}

with an energy-dependent prefactor.

Take units where

ℏ=1,2m=1,V0=1,E=14,a=6.\hbar=1, \qquad 2m=1, \qquad V_0=1, \qquad E=\frac14, \qquad a=6.

Then

α=V0−E=32,\alpha = \sqrt{V_0-E} = \frac{\sqrt3}{2},

and

2αa=63≈10.39.2\alpha a = 6\sqrt3 \approx 10.39.

The WKB exponent gives

TWKB≈e−10.39≈3.1×10−5.T_{\mathrm{WKB}} \approx e^{-10.39} \approx 3.1\times10^{-5}.

The opaque-limit exact prefactor is

16E(V0−E)V02=3.\frac{16E(V_0-E)}{V_0^2} = 3.

Thus

Texact≈9.2×10−5T_{\mathrm{exact}} \approx 9.2\times10^{-5}

in the opaque approximation. The two estimates differ by a prefactor but agree on the exponential scale.

The leading WKB exponent is useful when

Kℏ≫1.\frac{K}{\hbar}\gg1.

For the rectangular barrier, edge matching is abrupt, so the prefactor is not expected to be captured by the simplest smooth-barrier WKB expression. The comparison is still valuable because it isolates the most robust part of tunneling: exponential suppression by the forbidden-region action.

The approximation becomes poor near the top of the barrier, E≈V0E\approx V_0, because α\alpha becomes small and the barrier is no longer opaque.

The exponent 2αa2\alpha a is dimensionless.

Increasing aa suppresses TT exponentially.

Increasing the mass mm suppresses TT because α\alpha grows like m\sqrt m.

Increasing EE toward V0V_0 reduces α\alpha and increases TT.

  • Using e−αae^{-\alpha a} for the probability instead of e−2αae^{-2\alpha a}.
  • Comparing the bare WKB exponent to the exact result and expecting the prefactor to match.
  • Calling the rectangular barrier a smooth WKB problem.
  • Forgetting that the left and right velocities matter if the asymptotic potentials differ.
  • Trusting the opaque approximation when αa\alpha a is not large.
  • L. D. Landau and E. M. Lifshitz, Quantum Mechanics: Non-Relativistic Theory, 3rd ed., Pergamon, 1977.
  • C. M. Bender and S. A. Orszag, Advanced Mathematical Methods for Scientists and Engineers, Springer, 1999.
  • D. J. Griffiths and D. F. Schroeter, Introduction to Quantum Mechanics, 3rd ed., Cambridge University Press, 2018.
  • R. Shankar, Principles of Quantum Mechanics, 2nd ed., Springer, 1994.
  1. In the same units as the numerical example, double the barrier width from a=6a=6 to a=12a=12. By what factor does the WKB transmission change?
Solution

The WKB probability scales as

TWKB∝e−2αa.T_{\mathrm{WKB}} \propto e^{-2\alpha a}.

Doubling aa multiplies the transmission by

e−2αa=e−63≈3.1×10−5e^{-2\alpha a} = e^{-6\sqrt3} \approx 3.1\times10^{-5}

relative to the original value.

  1. Explain why the exact opaque-limit prefactor does not invalidate the WKB exponent.
Solution

When K/ℏ≫1K/\hbar\gg1, the exponential factor can change by many orders of magnitude as parameters vary. A prefactor of order one, or even a modest numerical factor, is usually much less important than the exponent. The exact prefactor matters for precision, but the WKB exponent captures the leading asymptotic dependence.