Skip to content

Bohr–Sommerfeld Quantization

Bohr–Sommerfeld quantization is the WKB rule for bound-state energies in a one-dimensional smooth potential well. It says that the classical action over one closed orbit is quantized in units of 2πℏ2\pi\hbar, with a turning-point correction. For a complete application with an analytic action and a converged spectral benchmark, see WKB Bound States in a Smooth Potential.

Required background. Turning-point connection formulas supply the two boundary phases.

For a bound orbit between two isolated simple smooth turning points x1x_1 and x2x_2,

p(x)=2m(E−V(x)),V(x1)=V(x2)=E,p(x) = \sqrt{2m\left(E-V(x)\right)}, \qquad V(x_1)=V(x_2)=E,

the WKB quantization condition is

∫x1x2p(x) dx=πℏ(n+12),n=0,1,2,….\int_{x_1}^{x_2}p(x)\,dx = \pi\hbar \left( n+\frac12 \right), \qquad n=0,1,2,\ldots.

Equivalently,

∮p dx=2πℏ(n+12).\oint p\,dx = 2\pi\hbar \left( n+\frac12 \right).

The closed integral is over a full classical period in phase space.

Inside the well, the WKB wavefunction is oscillatory:

ψ(x)≈1p(x)[C+e(i/ℏ)∫xp(x′) dx′+C−e−(i/ℏ)∫xp(x′) dx′].\psi(x) \approx \frac{1}{\sqrt{p(x)}} \left[ C_+e^{(i/\hbar)\int^x p(x')\,dx'} + C_-e^{-(i/\hbar)\int^x p(x')\,dx'} \right].

Outside the well, normalizability selects decaying forbidden-region tails. Each tail connects through a turning point to a sine-like allowed-region wave with a π/4\pi/4 phase shift. Compatibility of the left and right matched forms requires the phase accumulated across the well to satisfy

1ℏ∫x1x2p(x) dx+π2=(n+1)π.\frac{1}{\hbar} \int_{x_1}^{x_2}p(x)\,dx + \frac{\pi}{2} = (n+1)\pi.

This is equivalent to

∫x1x2p(x) dx=πℏ(n+12).\int_{x_1}^{x_2}p(x)\,dx = \pi\hbar \left( n+\frac12 \right).

The 1/21/2 shift is not put in by hand. It comes from the two smooth turning points. Each contributes a phase of π/4\pi/4, so together they shift the standing-wave condition by π/2\pi/2.

This is why the WKB condition differs from the old rule

∮p dx=nh.\oint p\,dx=nh.

For smooth one-dimensional wells, the modern WKB rule includes the Maslov correction associated with the turning points.

For

V(x)=12mω2x2,V(x) = \frac12m\omega^2x^2,

the turning points are

x±=±2Emω2.x_\pm = \pm \sqrt{\frac{2E}{m\omega^2}}.

The closed classical action is

∮p dx=2πEω.\oint p\,dx = \frac{2\pi E}{\omega}.

Using the WKB rule gives

2πEω=2πℏ(n+12),\frac{2\pi E}{\omega} = 2\pi\hbar \left( n+\frac12 \right),

so

En=ℏω(n+12).E_n = \hbar\omega \left( n+\frac12 \right).

For the harmonic oscillator, WKB gives the exact energy spectrum. This exactness is special; most potentials receive higher-order corrections.

WKB is usually most accurate for highly excited states because the action is large compared with ℏ\hbar. The rule estimates the zero-based level index associated with energy EE as

n(E)≈1πℏ∫x1(E)x2(E)p(x;E) dx−12.n(E) \approx \frac{1}{\pi\hbar} \int_{x_1(E)}^{x_2(E)}p(x;E)\,dx - \frac12.

The number of levels below EE is a staircase, not this index. Away from approximate eigenvalues and assuming a single sequence indexed from zero, leading WKB predicts

N<(E)≈⌊1πℏ∫x1x2p dx+12⌋.N_{<}(E)\approx\left\lfloor\frac{1}{\pi\hbar} \int_{x_1}^{x_2}p\,dx+\frac12\right\rfloor.

At a threshold, distinguish strict from non-strict counting; WKB errors can move a step. A smooth leading density discards this order-one staircase ambiguity.

The spacing of high-lying levels is controlled by the classical period. Since

ddE∮p dx=τ(E),\frac{d}{dE} \oint p\,dx = \tau(E),

where τ(E)\tau(E) is the classical period, the energy spacing is approximately

ΔE≈2πℏτ(E).\Delta E \approx \frac{2\pi\hbar}{\tau(E)}.

This connects the density of quantum levels to classical motion.

For the undergraduate finite-well counting version of this idea, see Bound-State Counting.

The formula above assumes two smooth turning points. It should not be blindly applied to hard-wall boxes, singular potentials, or discontinuous potentials. The infinite square well, for example, is better understood through exact boundary conditions; its walls are not smooth WKB turning points.

When endpoints are singular or when motion occurs on a multidimensional invariant torus, the phase correction must be modified. Those generalizations keep the same action-quantization spirit but require more geometry than the one-dimensional smooth-well rule.

  1. A particle moves in V(x)=FxV(x)=Fx for x>0x>0 and obeys ψ(0)=0\psi(0)=0, with F>0F>0. Find its leading quantization rule and energies.
Solution

There is one smooth turning point at E/FE/F and one hard wall. The decaying right tail gives sin⁡[ℏ−1∫xE/Fp dx+π/4]\sin[\hbar^{-1}\int_x^{E/F}p\,dx+\pi/4]. At the wall it must vanish, so

∫0E/Fp dx=πℏ(n+3/4),n=0,1,….\int_0^{E/F}p\,dx=\pi\hbar(n+3/4),\qquad n=0,1,\ldots.

The integral is 22mE3/2/(3F)2\sqrt{2m}E^{3/2}/(3F), hence

EnWKB=[3πℏF22m(n+3/4)]2/3.E_n^{\mathrm{WKB}}= \left[\frac{3\pi\hbar F}{2\sqrt{2m}}(n+3/4)\right]^{2/3}.

This is the large-index asymptotic spectrum of the negative zeros of Ai⁡\operatorname{Ai}, not the two-smooth-turning-point n+1/2n+1/2 rule.

  1. Derive dJ/dE=τdJ/dE=\tau for J(E)=2∫x1(E)x2(E)p dxJ(E)=2\int_{x_1(E)}^{x_2(E)}p\,dx.
Solution

The endpoint terms vanish because p(xi)=0p(x_i)=0. Since ∂Ep=m/p=1/v\partial_Ep=m/p=1/v, differentiation gives dJ/dE=2∫x1x2dx/v=τdJ/dE=2\int_{x_1}^{x_2}dx/v=\tau. For simple turning points the integrable endpoint singularities cause no divergence. Near a separatrix the period may diverge and the usual smooth level-spacing expansion needs extra care.

  • M. V. Berry and K. E. Mount, “Semiclassical approximations in wave mechanics”, Reports on Progress in Physics 35, 315-397, 1972.
  • L. D. Landau and E. M. Lifshitz, Quantum Mechanics: Non-Relativistic Theory, 3rd ed., Butterworth-Heinemann, 1981.
  • A. Messiah, Quantum Mechanics, Dover, 1999.
  • R. Shankar, Principles of Quantum Mechanics, 2nd ed., Springer, 1994.