WKB Quantization and the Bohr-Sommerfeld Condition

The WKB quantization condition is the rule that turns the semiclassical wavefunction into a spectrum. Take the WKB solution inside a potential well, demand that it decay into the forbidden regions beyond both turning points, and consistency forces the action accumulated over one classical period to be quantized: ∮p dx=(n+12)h\oint p\,\mathrm{d}x=(n+\tfrac12)h. This is the Bohr–Sommerfeld condition — the heart of the old quantum theory — now derived, and crucially corrected by the 12\tfrac12, from wave mechanics. The article owns that derivation. The WKB wavefunction, its p−1/2p^{-1/2} amplitude, and its turning-point breakdown are supplied by The WKB Approximation; the connection formulae that carry a π/4\pi/4 phase across each soft turning point are supplied by WKB Connection Formulae; the raw phase-integral rule and the Maslov ledger belong to Bohr–Sommerfeld and Semiclassical Quantization in 1D. What we derive here is how those pieces combine: matching the decaying tails through both turning points, adding the two π/4\pi/4 shifts, and reading off the half-integer. We prove where the 12\tfrac12 comes from — it is not Bohr's integer but the sum of two Maslov quarter-turns — and we make the counting explicit for wells with soft turning points, hard walls, and mixtures. The centerpiece is a structural theorem: WKB is exact for the harmonic oscillator. We evaluate the action integral ∫−aa2mE−m2ω2x2 dx=πE/ω\int_{-a}^{a}\sqrt{2mE-m^2\omega^2x^2}\,\mathrm{d}x=\pi E/\omega in closed form and land on En=(n+12)ℏωE_n=(n+\tfrac12)\hbar\omega with zero error for every nn — a genuine coincidence of the two-soft-turning-point structure, proved symbolically. Against it we place the honest contrast: for a hard-walled box the correct Maslov count (π/2\pi/2 per wall) reproduces the exact spectrum too, but the naive soft-turning-point rule misapplied to the box is only asymptotically right, off by 75%75\% at the ground state. The condition also carries a deep reading: each quantum state occupies a phase-space cell of area hh. That single sentence is the semiclassical density of states, the correspondence principle (ΔE→ℏωcl\Delta E\to\hbar\omega_{\mathrm{cl}}), and the bridge to statistical mechanics. Tunneling — the analytic continuation into the forbidden region — is deferred by name to WKB Tunneling and the Gamow Factor.

Key takeaways
  • The Bohr–Sommerfeld / WKB quantization condition is ∮p dx=(n+12)h\oint p\,\mathrm{d}x=(n+\tfrac12)h, equivalently ∫x1x2p dx=(n+12)πℏ\int_{x_1}^{x_2}p\,\mathrm{d}x=(n+\tfrac12)\pi\hbar between the two turning points.
  • It is derived by matching the WKB solution to decaying exponentials beyond both turning points, using the connection formulae; each soft turning point contributes a π/4\pi/4 phase, and the two sum to the 12\tfrac12 (the Maslov index).
  • The 12\tfrac12 is the wave-mechanical correction to Bohr's integer rule ∮p dx=nh\oint p\,\mathrm{d}x=nh — it is Maslov bookkeeping, not an arbitrary offset.
  • WKB is EXACT for the harmonic oscillator: the action integral evaluates to πE/ω\pi E/\omega, giving En=(n+12)ℏωE_n=(n+\tfrac12)\hbar\omega for all nn with zero error (proved symbolically).
  • The Maslov count is potential-dependent: two soft turning points give (n+12)(n+\tfrac12); two hard walls give the box rule (exact); a soft-plus-hard pair gives (n+34)(n+\tfrac34). Using 12\tfrac12 blindly for a box is off by 75%75\% at n=0n=0.
  • The condition says each state occupies a phase-space cell of area hh — the semiclassical density of states, and the correspondence principle ΔE→ℏωcl\Delta E\to\hbar\omega_{\mathrm{cl}}.
  • Key idea: match the WKB wave to decaying tails through both turning points; the two π/4\pi/4 connection shifts force ∮p dx=(n+12)h\oint p\,\mathrm{d}x=(n+\tfrac12)h, exact for the oscillator and asymptotic in general.
The core insight

Between two classical turning points the WKB wave must connect to decaying exponentials on both sides. Applying the connection formula at each end attaches a −π/4-\pi/4 phase, and demanding a single consistent oscillation in between forces the accumulated action ∫x1x2p dx\int_{x_1}^{x_2}p\,\mathrm{d}x to equal (n+12)πℏ(n+\tfrac12)\pi\hbar. The 12\tfrac12 is the sum of the two quarter-turns — the Maslov index — and it makes WKB exact for the harmonic oscillator.

Where this lives today

The half-integer quantization rule is the semiclassical spine of modern physics. Its immediate application is molecular spectroscopy: vibrational term values of diatomics come straight from the Bohr–Sommerfeld integral, and the RKR (Rydberg–Klein–Rees) method runs it backwards to reconstruct a potential curve from measured levels. In atomic physics the same integral with the Langer replacement ℓ(ℓ+1)→(ℓ+12)2\ell(\ell+1)\to(\ell+\tfrac12)^2 recovers the Bohr energies and explains quantum defects of penetrating states. The phase-space-cell reading — one state per area hh — is the semiclassical density of states behind the Thomas–Fermi atom and the classical limit of statistical mechanics. Generalized to several degrees of freedom it becomes EBK (Einstein–Brillouin–Keller) torus quantization for integrable systems and, for chaotic ones, Gutzwiller's periodic-orbit trace formula. The tunneling continuation of this same machinery is WKB Tunneling and the Gamow Factor; the raw phase-integral rule and its Maslov ledger are Bohr–Sommerfeld and Semiclassical Quantization in 1D, which this article's derivation consumes by name.

The physical picture

Picture a marble rolling in a smooth bowl. Where it moves fast (deep in the well) the de Broglie wavelength λ=h/p\lambda=h/p is short and the wave wiggles rapidly; near the turning points the marble slows, λ\lambda stretches, and the amplitude swells. The bound state is a standing wave: as the disturbance runs from one turning point to the other and back, it must return to itself in phase. The accumulated phase is ∮p dx/ℏ\oint p\,\mathrm{d}x/\hbar — but each reflection off a soft turning point is not a hard bounce; the connection formula shows it costs an extra π/4\pi/4. Two soft reflections per round trip cost π/2\pi/2, which is exactly the 12\tfrac12. The cleanest way to see the rule is in phase space. A classical orbit of energy EE traces a closed loop in the (x,p)(x,p) plane, and ∮p dx\oint p\,\mathrm{d}x is the area it encloses. The quantization condition says the allowed orbits enclose areas (n+12)h(n+\tfrac12)h — one Planck cell per state, offset by half a cell. For the harmonic oscillator these loops are ellipses of area π xmax⁡pmax⁡=2πE\pi\,x_{\max}p_{\max}=2\pi E, so the levels sit at equal area intervals of 2πℏ=h2\pi\hbar=h, and the 12\tfrac12-cell offset is the zero-point energy. The figure shows the nested ellipses and the staircase of enclosed areas.

Figure
Phase-space orbits enclose (n+½)h: nested oscillator ellipses and the area staircase
each orbit encloses ∮p dx = (n+½)h — one Planck cell per state x p n=0 n=1 n=2 areas: ½h, 3⁄2 h, 5⁄2 h — equal spacing h between levels ∮p dx /h n ½ 3⁄2 5⁄2 7⁄2 the ½-cell offset = zero-point energy

There is a subtlety the figure cannot show but the derivation makes precise: the 12\tfrac12 is not universal. It is the sum of Maslov contributions, one quarter per soft turning point. A well with two soft turning points (like the oscillator) gives 12\tfrac12; a half-well against a hard wall gives 34\tfrac34; a hard-walled box gives the integer-node count. The second figure contrasts the two headline cases — the oscillator, where WKB is exact, and the box, where the soft rule misapplied is only asymptotic — plotting the fractional error against quantum number.

Figure
Where WKB is exact and where it is only asymptotic: oscillator vs the box with the wrong Maslov count
fractional error of WKB energy vs quantum number n n error harmonic oscillator: 0% for every n (WKB exact) box with the SOFT (n+½) rule misapplied: 75% at n=0 → 9.8% at n=9 (the box is EXACT with the correct hard-wall count; the error above is from using the wrong Maslov index)

The core equations

We assemble the condition from pieces the neighbours own. The starting point is the WKB wavefunction in a classically allowed well, which oscillates with the local momentum and amplitude p−1/2p^{-1/2}. In an allowed region (E>VE>V) the general WKB solution is a two-term oscillation.

ψ(x)≈1p(x)[A eiℏ∫xp dx′+B e−iℏ∫xp dx′](1)\psi(x)\approx \frac{1}{\sqrt{p(x)}}\left[A\,e^{\frac{i}{\hbar}\int^x p\,\mathrm{d}x'}+B\,e^{-\frac{i}{\hbar}\int^x p\,\mathrm{d}x'}\right] \qquad (1)

where: A,BA,B — constants fixed by the connection conditions; p=2m(E−V)p=\sqrt{2m(E-V)}. Physical read: a standing wave whose local wavelength tracks the classical momentum. Beyond a turning point the region is forbidden and a bound state must keep only the decaying exponential.

ψ(x)≈D∣p(x)∣ e−1ℏ∫xtx∣p∣ dx′(2)\psi(x)\approx \frac{D}{\sqrt{|p(x)|}}\,e^{-\frac{1}{\hbar}\int_{x_t}^{x}|p|\,\mathrm{d}x'} \qquad (2)

where: xtx_t — a turning point; DD — amplitude of the decaying tail. Physical read: normalizability forbids the growing branch. Crossing a soft turning point relates the decaying tail to a cosine on the allowed side, with a definite π/4\pi/4 phase — the connection formula (derived in WKB Connection Formulae, consumed here).

D∣p∣ e−1ℏ∫xxt∣p∣ dx′ ⟷ 2Dpcos⁡ ⁣(1ℏ∫xtxp dx′−π4)(3)\frac{D}{\sqrt{|p|}}\,e^{-\frac{1}{\hbar}\int_{x}^{x_t}|p|\,\mathrm{d}x'}\ \longleftrightarrow\ \frac{2D}{\sqrt{p}}\cos\!\left(\frac{1}{\hbar}\int_{x_t}^{x}p\,\mathrm{d}x'-\frac{\pi}{4}\right) \qquad (3)

where: the factor 22 and the −π/4-\pi/4 come from the Airy asymptotics; this is the piece owned by the connection-formula article. Combining the two connections gives the quantization condition, written between the turning points or as a loop.

∫x1x2p(x) dx=(n+12)πℏ ⟺ ∮p dx=(n+12)h,n=0,1,2,…(4)\int_{x_1}^{x_2}p(x)\,\mathrm{d}x=\left(n+\tfrac12\right)\pi\hbar\ \Longleftrightarrow\ \oint p\,\mathrm{d}x=\left(n+\tfrac12\right)h,\quad n=0,1,2,\dots \qquad (4)

where: x1,x2x_1,x_2 — the two turning points of the well; nn — the number of interior nodes; ∮\oint — one full classical period. Physical read: allowed orbits enclose (n+12)(n+\tfrac12) Planck cells of phase-space area.

Derivation ABohr–Sommerfeld from the two connection formulae

The goal is to obtain the quantization condition by requiring a WKB bound state to decay beyond both turning points simultaneously. We consume the connection formula (owned next door) at each end, write the interior oscillation two ways — connecting inward from the right turning point and inward from the left — and demand the two descriptions agree. Consistency of a single physical wave forces the accumulated action to a half-integer multiple of πℏ\pi\hbar. The two π/4\pi/4 Maslov shifts, one per soft turning point, are what produce the 12\tfrac12.

Result. ∫x1x2p dx=(n+12)πℏ\int_{x_1}^{x_2}p\,\mathrm{d}x=(n+\tfrac12)\pi\hbar.

Step 1. Beyond the right turning point x2x_2 keep only the decaying tail, then connect inward: the interior wave is ψ∝2pcos⁡ ⁣(1ℏ∫xx2p dx′−π4)\psi\propto \tfrac{2}{\sqrt p}\cos\!\big(\tfrac1\hbar\int_x^{x_2}p\,\mathrm{d}x'-\tfrac{\pi}{4}\big). Why: normalizability kills the growing exponential on the right; the connection formula attaches a −π/4-\pi/4 as the wave enters the allowed region.

Step 2. Beyond the left turning point x1x_1 do the same, connecting inward from the left: ψ∝2pcos⁡ ⁣(1ℏ∫x1xp dx′−π4)\psi\propto \tfrac{2}{\sqrt p}\cos\!\big(\tfrac1\hbar\int_{x_1}^{x}p\,\mathrm{d}x'-\tfrac{\pi}{4}\big). Why: the identical logic on the left produces its own −π/4-\pi/4; both turning points are soft, so both contribute a quarter-turn.

Step 3. These two cosines describe the same interior wave, so their phases may differ only by a multiple of π\pi (up to an overall sign).

[1ℏ∫x1xp dx′−π4]+[1ℏ∫xx2p dx′−π4]=nπ\left[\frac{1}{\hbar}\int_{x_1}^{x}p\,\mathrm{d}x'-\frac{\pi}{4}\right]+\left[\frac{1}{\hbar}\int_{x}^{x_2}p\,\mathrm{d}x'-\frac{\pi}{4}\right]=n\pi

Why: a single-valued standing wave cannot have two inconsistent phase descriptions; the freedom is exactly nπn\pi, which is what indexes the states.

Step 4. The xx-dependent parts combine into the full action integral; the two quarter-turns add to π/2\pi/2.

1ℏ∫x1x2p dx−π2=nπ\frac{1}{\hbar}\int_{x_1}^{x_2}p\,\mathrm{d}x-\frac{\pi}{2}=n\pi

Why: ∫x1x+∫xx2=∫x1x2\int_{x_1}^{x}+\int_{x}^{x_2}=\int_{x_1}^{x_2} independent of the interior point xx; the two −π/4-\pi/4's are the Maslov sum.

Step 5. Solve.

∫x1x2p dx=(n+12)πℏ ⟺ ∮p dx=(n+12)h\int_{x_1}^{x_2}p\,\mathrm{d}x=\left(n+\tfrac12\right)\pi\hbar\ \Longleftrightarrow\ \oint p\,\mathrm{d}x=\left(n+\tfrac12\right)h

Why: the 12\tfrac12 is the sum of the two soft-turning-point Maslov quarter-turns — the wave-mechanical correction to Bohr's integer rule ∮p dx=nh\oint p\,\mathrm{d}x=nh. *Verified: verify_aq-qm-wkb-quantization-condition.py confirms the two π/4\pi/4 shifts sum to the rational 12\tfrac12, and that the next derivation reproduces the exact oscillator spectrum from this rule.*

Go deeper · Derivation B · 1 minWKB is exact for the harmonic oscillator
Derivation BWKB is exact for the harmonic oscillator

This is the structural theorem the article owns: not merely that WKB is good for the oscillator, but that it is exact for every level, and why. The engine is a single closed-form integral. We compute the action ∫−aap dx\int_{-a}^{a}p\,\mathrm{d}x for the parabolic potential symbolically, insert it into the quantization condition, and land on the exact ladder En=(n+12)ℏωE_n=(n+\tfrac12)\hbar\omega with identically zero error — a genuine coincidence of the two-soft-turning-point structure, not an approximation that happens to be close.

Result. En=(n+12)ℏωE_n=(n+\tfrac12)\hbar\omega, exact for all nn.

Step 1. For V=12mω2x2V=\tfrac12 m\omega^2 x^2 the turning points are x1,2=∓ax_{1,2}=\mp a with a=2E/(mω2)a=\sqrt{2E/(m\omega^2)}, and p=2mE−m2ω2x2p=\sqrt{2mE-m^2\omega^2x^2}. Why: E=VE=V fixes the classical amplitude aa; the momentum vanishes there and is real in between.

Step 2. Evaluate the action integral in closed form (substitute x=asin⁡θx=a\sin\theta).

∫−aa2mE−m2ω2x2 dx=πEω\int_{-a}^{a}\sqrt{2mE-m^2\omega^2x^2}\,\mathrm{d}x=\frac{\pi E}{\omega}

Why: the integral is the area of a half-ellipse in (x,p)(x,p); the mm's cancel and it collapses to πE/ω\pi E/\omega exactly — verified symbolically, not numerically.

Step 3. Impose the quantization condition (4).

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

Why: the oscillator has two soft turning points, so the (n+12)(n+\tfrac12) form applies without modification.

Step 4. Solve for the spectrum.

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

Why: this matches the exact ladder-operator result for every nn, zero-point energy included. The exactness is special to the quadratic potential — for a general smooth well the higher-order ℏ\hbar corrections to the phase do not vanish and WKB is only asymptotically accurate. *Verified: verify_aq-qm-wkb-quantization-condition.py evaluates ∫−aap dx\int_{-a}^{a}p\,\mathrm{d}x symbolically to πE/ω\pi E/\omega and solves to En=(n+12)ℏωE_n=(n+\tfrac12)\hbar\omega exactly; it also checks dE/dn=ℏω\mathrm{d}E/\mathrm{d}n=\hbar\omega (the correspondence-principle spacing).*

Limiting cases

The quantization shift depends on the boundary structure; the validity depends on how slowly the wavelength varies (owned by The WKB Approximation).

CaseResultMeaning
Two soft turning points∮p dx=(n+12)h\oint p\,\mathrm{d}x=(n+\tfrac12)hGeneric smooth well (oscillator, Morse). Maslov index μ=2\mu=2, shift μ/4=12\mu/4=\tfrac12.
One soft, one hard wall∮p dx=(n+34)h\oint p\,\mathrm{d}x=(n+\tfrac34)hBouncing ball, half-harmonic well. Soft gives 14\tfrac14, hard wall (node) gives 12\tfrac12.
Two hard walls (box)∮p dx=(n+1)h\oint p\,\mathrm{d}x=(n+1)hInfinite square well. WKB with this count reproduces Ek=k2π2ℏ2/2mL2E_k=k^2\pi^2\hbar^2/2mL^2 exactly, k=n+1k=n+1.
Large quantum numberΔEn→ℏωcl(En)=2πℏT(En)\Delta E_n\to\hbar\omega_{\mathrm{cl}}(E_n)=\tfrac{2\pi\hbar}{T(E_n)}Level spacing equals the classical orbital frequency times ℏ\hbar — the correspondence principle.
Harmonic oscillatorEn=(n+12)ℏωE_n=(n+\tfrac12)\hbar\omegaThe action integral is exactly πE/ω\pi E/\omega; WKB is exact for all nn.

Worked examples

Example 1 — the oscillator action integral, end to end. Compute ∫−aap dx\int_{-a}^{a}p\,\mathrm{d}x for V=12mω2x2V=\tfrac12 m\omega^2x^2 and read off the spectrum. Work in natural units m=ω=ℏ=1m=\omega=\hbar=1 and take the level n=3n=3.

Step 1. At E=E3E=E_3 the turning points are a=2Ea=\sqrt{2E}; the momentum is p=2E−x2p=\sqrt{2E-x^2}.

Step 2. Substitute x=asin⁡θx=a\sin\theta: ∫−aa2E−x2 dx=2E∫−π/2π/2cos⁡2θ dθ=2E⋅π2=πE\int_{-a}^{a}\sqrt{2E-x^2}\,\mathrm{d}x=2E\int_{-\pi/2}^{\pi/2}\cos^2\theta\,\mathrm{d}\theta=2E\cdot\tfrac{\pi}{2}=\pi E. (With ℏ=ω=1\hbar=\omega=1 this is the general πE/ω\pi E/\omega.)

Step 3. Impose πE=(n+12)π\pi E=(n+\tfrac12)\pi, i.e. En=n+12E_n=n+\tfrac12.

Step 4. For n=3n=3: E3=3.5 ℏωE_3=3.5\,\hbar\omega, identical to the exact ladder-operator value — no error, as the theorem promises. The action integral evaluated to πE/ω\pi E/\omega symbolically confirms exactness for every nn, not just n=3n=3.

Example 2 — the box, right way and wrong way. Compare the WKB spectrum of an infinite square well of width LL (natural units m=ℏ=1m=\hbar=1, L=1L=1) computed with the correct hard-wall count versus the naive soft rule, at n=0n=0 and n=9n=9.

Step 1. The exact spectrum is Ek=k2π2/2E_k=k^2\pi^2/2, k=1,2,…k=1,2,\dots. The momentum is constant, p=2Ep=\sqrt{2E}, so ∫0Lp dx=L2E\int_0^L p\,\mathrm{d}x=L\sqrt{2E}.

Step 2. Correct (two hard walls): the rule is ∫0Lp dx=(n+1)πℏ\int_0^L p\,\mathrm{d}x=(n+1)\pi\hbar, giving E=(n+1)2π2/2E=(n+1)^2\pi^2/2. With k=n+1k=n+1 this is exact for every nn — zero error at n=0n=0 and n=9n=9 alike.

Step 3. Wrong (soft (n+12)(n+\tfrac12) rule misapplied): E=(n+12)2π2/2E=(n+\tfrac12)^2\pi^2/2. At n=0n=0 this gives 14π2/2\tfrac14\pi^2/2 versus the exact π2/2\pi^2/2 — a 75%75\% shortfall.

Step 4. At n=9n=9 the wrong rule gives (9.5)2(9.5)^2 versus the exact 10210^2: an error of 9.75%9.75\%. The mistake shrinks as nn grows (asymptotic recovery) but never vanishes — the honest lesson that the Maslov count, not a universal 12\tfrac12, is what makes WKB exact for the box.

Physicist

Origins · 1 minBohr-Sommerfeld Quantization and how this was discovered
BQ
Bohr-Sommerfeld Quantization
1913–1926 · Copenhagen, Munich

The condition ∮p dq=nh\oint p\,\mathrm{d}q=nh was the organizing principle of the old quantum theory. Bohr's 1913 hydrogen model quantized angular momentum; Sommerfeld generalized it in 1915–16 to ∮p dq=nh\oint p\,\mathrm{d}q=nh for each separable coordinate, quantizing elliptical Kepler orbits and — with a relativistic correction — even producing the hydrogen fine-structure formula before Dirac. The rule was empirically astonishing and conceptually opaque: no one knew why phase-space areas should be integer multiples of hh. Wave mechanics answered that in 1926: the WKB analysis reproduces the rule, explains it as a standing-wave condition, and corrects the integer to n+12n+\tfrac12 through the Maslov phases. The old quantum theory was not so much overthrown as derived — with a half-quantum of interest.

History
1913
— Bohr quantizes angular momentum for the hydrogen atom
1915
— Sommerfeld generalizes to ∮ p dq = nh for separable coordinates
1916
— Sommerfeld's relativistic elliptical orbits reproduce the hydrogen fine-structure formula
1926
— WKB analysis rederives the rule within wave mechanics and adds the ½ Maslov correction
1958
— Maslov formalizes the index counting the phase loss at caustics and turning points

Going further

References. Griffiths & Schroeter, Introduction to Quantum Mechanics, 3rd ed., Ch. 9 · Sakurai & Napolitano, Modern Quantum Mechanics, 2nd ed., Sec. 2.5 · Shankar, Principles of Quantum Mechanics, 2nd ed., Ch. 16 · Cohen-Tannoudji, Diu & Laloë, Quantum Mechanics, Complement H_III / B_XI · Townsend, A Modern Approach to Quantum Mechanics, Ch. 10. This article owns the derivation of the quantization condition and the WKB-exact-for-the-oscillator theorem. The WKB wavefunction and validity condition are The WKB Approximation; the π/4\pi/4 connection formulae are WKB Connection Formulae, consumed here; the phase-integral rule and Maslov ledger (the ν=2/3/4\nu=2/3/4 toll and the quartic-well accuracy staircase 18.2%→1.3%→0.6%→0.3%18.2\%\to1.3\%\to0.6\%\to0.3\% and the neutron-bouncer levels) are Bohr–Sommerfeld and Semiclassical Quantization in 1D, recalled by name and never re-owned. Tunneling is WKB Tunneling and the Gamow Factor; the Langer-correction radial theory and higher-order quantization corrections are deferred to the graduate register.

**Provenance.*All equations are standard semiclassical results. A verifier should: (1) confirm the ℏ\hbar-expansion of −ℏ22mψ′′+Vψ=Eψ-\frac{\hbar^2}{2m}\psi''+V\psi=E\psi yields the eikonal (S0′)2=p2(S_0')^2=p^2 and transport S1=i2ln⁡pS_1=\frac{i}{2}\ln p equations (Griffiths 3e Ch. 9; Shankar Ch. 16); (2) recompute the harmonic-oscillator action integral ∫−aap dx=πE/ω\int_{-a}^{a}p\,\mathrm{d}x=\pi E/\omega giving En=ℏω(n+1/2)E_n=\hbar\omega(n+1/2) exactly; (3) verify the Airy-asymptotic connection formulas with the π/4\pi/4 shift (Sakurai 2e Sec. 2.5; Cohen-Tannoudji Complement BXI); (4) check the Maslov-index rule (1/4 per soft turning point, 1/2 per hard wall) against the infinite square well result En=(n+1)2π2ℏ2/2mL2E_n=(n+1)^2\pi^2\hbar^2/2mL^2; (5) confirm Bohr energies −13.6 eV/n2-13.6\,\mathrm{eV}/n^2 from the radial integral with the Langer correction. Numerical constant 13.6 eV13.6\,\mathrm{eV} is the standard Rydberg energy.

Builds on: Time-Independent Schrödinger Equation · Classical Action and Hamilton–Jacobi Theory · Airy Functions · Correspondence Principle Leads to: WKB Tunneling and the Gamow Factor · Fine Structure of Hydrogen · Non-Degenerate Time-Independent Perturbation Theory

Why is it called both "Bohr–Sommerfeld" and "WKB" quantization?

They are the same rule reached from two directions. Bohr (1913) and Sommerfeld (1915–16) postulated ∮p dq=nh\oint p\,\mathrm{d}q=nh as the quantization axiom of the old quantum theory. WKB (1926) derives it from wave mechanics as a standing-wave/connection-formula condition — and in doing so corrects the integer to n+12n+\tfrac12 via the Maslov phases. "Bohr–Sommerfeld" honours the historical postulate; "WKB quantization" honours the derivation that explained and corrected it.

Is the ½ always ½?

No — this is the most important caveat. The additive constant is μ/4\mu/4 where μ\mu is the Maslov index, counting 11 per soft turning point and 22 per hard wall (a node). Two soft turning points give 12\tfrac12; a soft-plus-hard pair gives 34\tfrac34; two hard walls give the integer-node box count. Blindly using 12\tfrac12 for a box is off by 75%75\% at the ground state.

If WKB is approximate, how can it be exact for the oscillator?

The WKB error lives in the turning-point regions, where the potential is linearized and matched to an Airy function. For the quadratic oscillator the action integral is exactly πE/ω\pi E/\omega and the two π/4\pi/4 phases are exactly right, so the quantization integral has no residual error. It is a genuine coincidence of the parabolic potential — the higher-order ℏ\hbar corrections to the phase happen to sum to zero. No other elementary potential shares this except through similar special structure.

Why does each state occupy phase-space area h?

Differentiating ∮p dx=(n+12)h\oint p\,\mathrm{d}x=(n+\tfrac12)h with respect to nn gives d(∮p dx)=h\mathrm{d}(\oint p\,\mathrm{d}x)=h per unit increase in nn: consecutive orbits enclose areas differing by exactly hh. So the number of states below energy EE is the enclosed phase-space area divided by hh. This is the semiclassical density of states and the origin of the "one state per Planck cell" rule that underlies statistical mechanics' classical limit.

How does WKB handle the hydrogen atom?

Apply the radial integral ∮pr dr=(nr+12)h\oint p_r\,\mathrm{d}r=(n_r+\tfrac12)h with the Langer replacement ℓ(ℓ+1)→(ℓ+12)2\ell(\ell+1)\to(\ell+\tfrac12)^2 in the centrifugal term. Without Langer the WKB energies are wrong; with it they reproduce the exact Bohr levels En=−13.6 eV/n2E_n=-13.6\,\mathrm{eV}/n^2. The Langer correction fixes the r→0r\to0 singularity where a naive linearization fails, and its full justification is a graduate-register topic.

Does the derivation assume both turning points are "soft"?

Yes — the two-connection-formula derivation attaches a π/4\pi/4 at each linear (soft) turning point, giving the (n+12)(n+\tfrac12) result. If a boundary is a hard wall the wavefunction must vanish there (a node), which is a different phase condition (π/2\pi/2 contribution), and the count changes accordingly. The derivation here is the two-soft-turning-point case; the limiting-cases table lists the others.

Before this

Where this leads