Skip to content

Semiclassical Limits and Correspondence

A semiclassical limit is a controlled comparison between a family of quantum predictions and a classical description. It is not the assertion that a quantum state literally turns into a point in phase space, nor is it obtained by deleting every occurrence of ℏ\hbar from an equation.

The word correspondence is therefore incomplete until one specifies what is being compared. Energy spacings may approach classical frequencies, smoothed eigenstate densities may approach classical time-spent measures, packet centers may follow classical trajectories for a finite interval, and oscillatory propagators may be organized by classical paths. These are distinct statements with distinct errors and failure modes. The suppression of Bose and Fermi exchange at low phase-space density is a separate limit developed in Classical Limit of Quantum Statistics.

Helpful background. The classical limit distinguishes the physical mechanisms; semiclassical asymptotics explains the small action ratio. Here we compare spectra, densities, and packets while stating which state, observable, and time scales remain fixed.

Different limits agree with different classical quantities:

Part Question to state explicitly Typical choices
Parameter family What sequence of quantum problems is considered? ϵ=ℏ/S0→0\epsilon=\hbar/S_0\to0, n→∞n\to\infty, large spin, short wavelength
Fixed data What remains unchanged along the sequence? Classical action, energy surface, geometry, scaled potential, observation time
State class Which preparations are admitted? WKB eigenstates, coherent packets, Lagrangian states, thermal mixtures
Observable class What is actually measured? Smooth functions, finite-resolution projectors, transition amplitudes, spectral windows
Time window For how long is the comparison asserted? Fixed time, one period, algebraic time, Ehrenfest time
Sense of approximation How is agreement quantified? Relative error, weak convergence, coarse-grained convergence, asymptotic series

For example, high-energy box densities approach the uniform distribution after averaging, although they keep oscillating pointwise. A localized packet can follow an orbit for each fixed time while losing localization on a time interval that grows with the small-parameter limit.

Because ℏ\hbar has dimensions of action, the notation ℏ→0\hbar\to0 is shorthand for a dimensionless asymptotic family. Consider

[−ℏ22m∇x2+V(x)]ψ(x)=Eψ(x).\left[ -\frac{\hbar^2}{2m}\nabla_x^2 + V(\mathbf x) \right]\psi(\mathbf x) = E\psi(\mathbf x).

Choose a length LL and an energy E0E_0, and define

y=xL,p0=mE0,S0=p0L.\mathbf y = \frac{\mathbf x}{L}, \qquad p_0 = \sqrt{mE_0}, \qquad S_0 = p_0L.

With

v(y)=V(Ly)E0,E=EE0,ϵ=ℏS0,\begin{aligned} v(\mathbf y) &= \frac{V(L\mathbf y)}{E_0}, \\ \mathcal E &= \frac{E}{E_0}, \\ \epsilon &= \frac{\hbar}{S_0}, \end{aligned}

the stationary Schrödinger equation becomes

[−ϵ22∇y2+v(y)]ψ(y)=Eψ(y).\left[ -\frac{\epsilon^2}{2}\nabla_y^2 + v(\mathbf y) \right]\psi(\mathbf y) = \mathcal E\psi(\mathbf y).

Now ϵ→0\epsilon\to0 is meaningful: the scaled potential, geometry, and scaled energy can be held fixed while the wavelength becomes short compared with LL. The same dimensionless family can be approached physically by increasing a mass, momentum, length, or collective action scale. Those realizations need not have identical experimental constraints, but they share the same leading asymptotic structure when the nondimensionalized problem agrees.

The limit is usually singular. The highest derivative is multiplied by ϵ2\epsilon^2, so setting ϵ=0\epsilon=0 removes the differential order and cannot reproduce boundary conditions, turning-point layers, tunneling tails, or interference. WKB keeps the rapidly varying phase before expanding:

ψ=exp⁡[iϵ∑k=0∞ϵkSk].\psi = \exp\left[ \frac{i}{\epsilon} \sum_{k=0}^{\infty} \epsilon^k S_k \right].

The classical Hamilton–Jacobi equation appears at leading order, but the quantum solution still requires transport amplitudes, branch sums, phase indices, and uniform repairs. This is why a semiclassical limit is an asymptotic construction rather than a substitution.

Large Quantum Numbers and Fixed Classical Action

Section titled “Large Quantum Numbers and Fixed Classical Action”

Large quantum number is often a symptom of semiclassical scaling, not its definition. For a one-dimensional periodic orbit, use the action convention

I(E)=12π∮p dq.I(E) = \frac{1}{2\pi} \oint p\,dq.

The leading EBK rule is

In=ℏ(n+μ4),I_n = \hbar \left( n+\frac{\mu}{4} \right),

where μ\mu is the Maslov index of the closed cycle. To compare with a fixed classical orbit, take the joint limit

ℏ→0,n→∞,In→I.\hbar\to0, \qquad n\to\infty, \qquad I_n\to I.

Thus nn grows because a fixed classical action contains more quantum cells as ℏ\hbar decreases. Sending n→∞n\to\infty at fixed ℏ\hbar instead moves through different classical actions and often through different energies. That can also be a useful asymptotic regime, but it is not the same family unless an additional rescaling identifies the classical problems.

Let the classical energy be E=H(I)E=H(I) and define the angular frequency

ω(I)=dHdI.\omega(I) = \frac{dH}{dI}.

For a fixed integer rr, the EBK rule gives

In+r−In=rℏ.I_{n+r}-I_n = r\hbar.

For a smooth scalar quantization whose energy remainder is uniformly O(ℏ2)O(\hbar^2) near the chosen torus, Taylor expansion at fixed II yields

En+r−En=H(In+rℏ)−H(In)+O(ℏ2)=rℏ ω(In)+O(ℏ2).\begin{aligned} E_{n+r}-E_n &= H(I_n+r\hbar)-H(I_n)+O(\hbar^2) \\ &= r\hbar\,\omega(I_n) + O(\hbar^2). \end{aligned}

Therefore

En+r−Enℏ⟶rω(I).\frac{E_{n+r}-E_n}{\hbar} \longrightarrow r\omega(I).

This is a precise spectral form of correspondence: nearby Bohr frequencies approach harmonics of the classical orbital frequency. It does not say that all absolute level spacings vanish. For the harmonic oscillator, ΔE=ℏω\Delta E=\hbar\omega is independent of nn at fixed ℏ\hbar; for a box, adjacent spacings grow with nn. What often becomes small is the spacing relative to the total energy or to the experimental energy scale.

Write a classical observable along an invariant orbit as

A(I,θ)=∑r∈ZAr(I)eirθ,A(I,\theta) = \sum_{r\in\mathbb Z} A_r(I)e^{ir\theta},

with

Ar(I)=12π∫02πA(I,θ)e−irθ dθ.A_r(I) = \frac{1}{2\pi} \int_0^{2\pi} A(I,\theta)e^{-ir\theta}\,d\theta.

Under regularity and quantization assumptions, semiclassical correspondence gives the leading relation

⟨n+r∣A^∣n⟩∼Ar(In)\langle n+r\vert\hat A\vert n\rangle \sim A_r(I_n)

for fixed rr in the joint large-nn, small-ℏ\hbar limit. Phase conventions for the eigenstates can alter the phase assigned to each matrix element, but transition strengths and the matching of harmonic content remain physical. This relation refines the slogan that “quantum transitions reproduce classical radiation frequencies”: the energy differences supply the frequencies, while the matrix elements supply the corresponding classical Fourier amplitudes.

For several integrable degrees of freedom, II and rr become vectors:

En+r−Enℏ⟶r⋅ω(I),\frac{E_{\mathbf n+\mathbf r}-E_{\mathbf n}}{\hbar} \longrightarrow \mathbf r\cdot\boldsymbol\omega(\mathbf I),

on regular tori with a consistent local action lattice and controlled spectral remainders. Resonant frequency ratios can make the leading combination vanish; they do not by themselves invalidate the relation. Separatrices and singular tori need separate analysis.

A highly excited stationary state usually does not resemble a localized classical particle. Its density contains wavelength-scale oscillations, nodes, and interference. The appropriate classical comparison is often a probability measure on an energy shell, tested at finite spatial resolution.

For one-dimensional bound motion between turning points x−x_- and x+x_+, the normalized WKB wavefunction in the allowed region has the form

ψE(x)≃Cp(x;E)cos⁡ΦE(x),ΦE(x)=1ℏ∫x−xp(x′;E) dx′−π4.\begin{aligned} \psi_E(x) &\simeq \frac{C}{\sqrt{p(x;E)}} \cos\Phi_E(x), \\ \Phi_E(x) &= \frac{1}{\hbar} \int_{x_-}^{x}p(x';E)\,dx' -\frac{\pi}{4} . \end{aligned}

where

p(x;E)=2m[E−V(x)].p(x;E) = \sqrt{2m[E-V(x)]}.

A detector that averages over many local wavelengths but over a scale short compared with the variation of pp replaces cos⁡2\cos^2 by 1/21/2. Normalization then gives

∣ψE(x)∣2‾≃2T(E)v(x;E),\overline{\lvert\psi_E(x)\rvert^2} \simeq \frac{2}{T(E)v(x;E)},

where

v(x;E)=p(x;E)mv(x;E) = \frac{p(x;E)}{m}

and T(E)T(E) is the full classical period. The right-hand side is exactly the classical time-spent density: during one period the orbit crosses an interior position twice, spending total time 2 dx/v2\,dx/v there.

The statement is local only away from turning points. The WKB expression diverges as v→0v\to0, while the exact wavefunction remains finite and is described by an Airy uniform approximation. The integrated probability near the turning point can nevertheless have a controlled semiclassical limit. Pointwise formulas, weak limits, and uniform approximations answer different questions.

A box makes the mode of convergence visible

Section titled “A box makes the mode of convergence visible”

For an infinite well on 0<x<L0\lt x\lt L, define kn=2nπ/Lk_n=2n\pi/L. Then

ρn(x)=2Lsin⁡2(nπxL)=1L[1−cos⁡(knx)].\begin{aligned} \rho_n(x) &= \frac{2}{L} \sin^2\left( \frac{n\pi x}{L} \right) \\ &= \frac{1}{L} \left[ 1-\cos(k_nx) \right]. \end{aligned}

At a fixed interior point, the cosine generally keeps oscillating as n→∞n\to\infty, so ρn(x)\rho_n(x) has no pointwise limit. For a smooth detector response g(x)g(x), however,

∫0Lg(x)ρn(x) dx=1L∫0Lg(x) dx−1L∫0Lg(x)cos⁡(knx) dx.\int_0^L g(x)\rho_n(x)\,dx = \frac{1}{L} \int_0^L g(x)\,dx {}- \frac{1}{L} \int_0^L g(x) \cos(k_nx)\,dx.

The oscillatory integral tends to zero under mild regularity assumptions. Hence

ρn⇀1L,\rho_n \rightharpoonup \frac{1}{L},

which is weak convergence to the uniform classical position distribution.

A rapidly oscillating high-n box probability density around the uniform classical density, with a finite-width detector bin averaging over several oscillations.

A high-nn density need not settle point by point. A detector bin satisfying L/n≪Δx≪LL/n\ll\Delta x\ll L resolves the classical spatial scale but not each quantum fringe, so its averaged probability approaches the uniform classical measure.

Weak convergence is not permission to discard all fine structure. An observable deliberately tuned to the wavelength-scale oscillations can retain an order-one quantum signal. The observable class must therefore be fixed before the limit is taken.

Stationary-state correspondence concerns distributions and frequencies. Particle-like motion requires a different state class, usually a packet localized in both position and momentum.

For

H^=p^22m+V(x^),\hat H = \frac{\hat p^2}{2m} + V(\hat x),

Ehrenfest’s equations are exact:

d⟨x^⟩dt=⟨p^⟩m,d⟨p^⟩dt=−⟨V′(x^)⟩.\frac{d\langle\hat x\rangle}{dt} = \frac{\langle\hat p\rangle}{m}, \qquad \frac{d\langle\hat p\rangle}{dt} = -\langle V'(\hat x)\rangle.

They become a closed classical system only when

⟨V′(x^)⟩≃V′(⟨x^⟩).\langle V'(\hat x)\rangle \simeq V'(\langle\hat x\rangle).

Let

xˉ=⟨x^⟩,σx2=⟨(x^−xˉ)2⟩.\bar x = \langle\hat x\rangle, \qquad \sigma_x^2 = \langle(\hat x-\bar x)^2\rangle.

Expanding the force about the packet center gives

⟨V′(x^)⟩=V′(xˉ)+12V′′′(xˉ)σx2+16V(4)(xˉ)⟨(x^−xˉ)3⟩+⋯ .\langle V'(\hat x)\rangle = V'(\bar x) + \frac{1}{2} V'''(\bar x)\sigma_x^2 + \frac{1}{6} V^{(4)}(\bar x) \langle(\hat x-\bar x)^3\rangle + \cdots.

The packet center follows Newton’s equation while the higher-moment corrections remain small on the force scale of interest. The closure is exact for potentials at most quadratic, but generic packets spread, shear, or split in nonlinear potentials.

Three conclusions should be kept separate:

  1. Center correspondence: xˉ(t)\bar x(t) and pˉ(t)\bar p(t) remain close to one classical trajectory.
  2. Shape correspondence: the packet remains narrow and approximately Gaussian.
  3. Observable correspondence: selected expectation values agree with a classical ensemble, even if one-packet localization has failed.

The first can fail while the third remains useful. A semiclassical propagator may accurately sum many classical trajectories after a single Gaussian packet has stretched beyond recognition. Ehrenfest Theorem Revisited derives the hierarchy of moments, while Wave Packets and Classical Trajectories explains the localization and spreading criteria. Coherent-State Semiclassics Preview compares one-packet variational propagation with trajectory-ensemble methods.

A semiclassical approximation proved or derived for each fixed time need not remain accurate when time grows as ϵ\epsilon decreases. Schematically, a statement such as

lim⁡ϵ→0sup⁡0≤t≤T∣Qϵ(t)−Qcl(t)∣=0\lim_{\epsilon\to0} \sup_{0\le t\le T} \lvert Q_\epsilon(t)-Q_{\mathrm{cl}}(t) \rvert = 0

for every fixed TT does not establish the same result for T=T(ϵ)→∞T=T(\epsilon)\to\infty.

This order-of-limits issue is central:

lim⁡ϵ→0lim⁡t→∞Qϵ(t)andlim⁡t→∞lim⁡ϵ→0Qϵ(t)\lim_{\epsilon\to0} \lim_{t\to\infty} Q_\epsilon(t) \qquad \text{and} \qquad \lim_{t\to\infty} \lim_{\epsilon\to0} Q_\epsilon(t)

need not agree. At long times, small frequency errors accumulate into phase errors; neighboring classical trajectories shear apart; caustics proliferate; discrete quantum recurrences can become visible; and exponentially small terms may become comparable to the nominal leading approximation.

For regular motion, packet deformation is often controlled by algebraic shear and by derivatives of the frequencies with respect to action. For chaotic motion, a linearized separation grows locally as

δz(t)∼δz0eλt,\delta z(t) \sim \delta z_0 e^{\lambda t},

where λ\lambda is a positive Lyapunov exponent in the relevant region. If LzL_z is a classical phase-space scale at which localization is lost, the estimate

tE∼1λlog⁡(Lzδz0)t_{\mathrm E} \sim \frac{1}{\lambda} \log\left( \frac{L_z}{\delta z_0} \right)

defines an Ehrenfest time. For a minimum-uncertainty packet, δz0\delta z_0 scales with a power of ℏ\hbar, so tEt_{\mathrm E} grows only logarithmically with the inverse semiclassical parameter. The coefficient depends on which width, observable, flow direction, and phase-space norm define breakdown; there is no universal context-free prefactor.

The conclusion is not that “quantum trajectories become chaotic.” Closed quantum evolution is linear and does not possess classical trajectory separation in Hilbert space in that naive sense. Quantum signatures of classical chaos instead appear through spectral statistics, eigenfunction structure, periodic-orbit contributions, transport, scrambling diagnostics, and the time dependence of selected observables. Quantum Chaos Preview introduces those diagnostics.

Closed-System Semiclassics and Decoherence

Section titled “Closed-System Semiclassics and Decoherence”

Semiclassical approximation and decoherence solve different problems.

  • Semiclassical approximation organizes closed-system amplitudes in powers or asymptotic scales of ℏ/S0\hbar/S_0, often around classical trajectories or invariant manifolds.
  • Coarse graining restricts which oscillations or phase-space structures an observable resolves.
  • Decoherence suppresses locally observable interference between alternatives by entangling them with unobserved environmental degrees of freedom.

A closed high-action system can remain in a coherent superposition whose interference is visible to a sufficiently fine measurement. Conversely, environmental decoherence can make records robust even when a simple WKB expansion is not the right computational tool. In realistic classical behavior the mechanisms may cooperate, but neither should be smuggled into the assumptions of the other. See Decoherence and the Classical-Limit Bridge for the open-system analysis.

The classical object and the appropriate notion of agreement depend on the quantum question.

Quantum question Classical object Appropriate comparison Main method or repair
Bound-state energies Actions and invariant tori Level spacings, counting functions, smoothed spectral density Bohr–Sommerfeld Quantization, EBK Quantization
Highly excited eigenstate density Time-spent or invariant measure Weak or coarse-grained convergence WKB Approximation plus turning-point repair
Localized packet motion Trajectory and stability matrix Moments or smooth observables over a stated time Ehrenfest analysis, Gaussian methods, coherent-state semiclassics
Propagator or transition amplitude Boundary-value trajectories Oscillatory asymptotics including phase and prefactor Stationary Phase, Semiclassical Propagator
Tunneling probability Complex or forbidden classical action Logarithmic asymptotics of exponentially small terms Forbidden-region WKB and connection formulas
Chaotic spectrum or transport Unstable periodic orbits and classical correlations Smoothed, statistical, or finite-time quantities Periodic-orbit and quantum-chaos methods

The table also shows why “the quantum answer approaches the classical answer” is too coarse. A wavefunction is not a classical trajectory, an amplitude is not a probability density, and a discrete spectrum is not an orbit. Correspondence relates selected structures after the relevant representation and resolution have been specified.

1. Estimate the first nonclassical force correction

Section titled “1. Estimate the first nonclassical force correction”

Let a packet have mean xˉ\bar x, variance σx2\sigma_x^2, and vanishing third central moment. Expand ⟨V′(x^)⟩\langle V'(\hat x)\rangle through second order in the packet displacement. Apply the result to

V(x)=12mω2x2+αx4.V(x) = \frac{1}{2}m\omega^2x^2 + \alpha x^4.
Solution

Write x^=xˉ+δx^\hat x=\bar x+\delta\hat x, with ⟨δx^⟩=0\langle\delta\hat x\rangle=0. Taylor expansion gives

⟨V′(x^)⟩=V′(xˉ)+12V′′′(xˉ)σx2+⋯ .\langle V'(\hat x)\rangle = V'(\bar x) + \frac{1}{2} V'''(\bar x)\sigma_x^2 + \cdots.

For the quartic oscillator,

V′(x)=mω2x+4αx3,V′′′(x)=24αx.\begin{aligned} V'(x) &= m\omega^2x + 4\alpha x^3, \\ V'''(x) &= 24\alpha x. \end{aligned}

Hence

⟨V′(x^)⟩≃mω2xˉ+4αxˉ3+12αxˉσx2.\langle V'(\hat x)\rangle \simeq m\omega^2\bar x + 4\alpha\bar x^3 + 12\alpha\bar x\sigma_x^2.

For this quartic potential and zero third central moment, the displayed expression is exact, provided the moments exist. The last term is the width correction; it vanishes for a purely quadratic potential and grows as the packet spreads. Near xˉ=0\bar x=0 this particular correction vanishes by symmetry, so one should compare the first nonzero moment correction rather than divide mechanically by the classical force.

2. Distinguish pointwise and weak convergence in a box

Section titled “2. Distinguish pointwise and weak convergence in a box”

For

ρn(x)=1L[1−cos⁡(2nπxL)],\rho_n(x) = \frac{1}{L} \left[ 1-\cos\left( \frac{2n\pi x}{L} \right) \right],

show that no pointwise limit exists for generic fixed xx, but that for every continuously differentiable gg,

∫0Lg(x)ρn(x) dx⟶1L∫0Lg(x) dx.\int_0^L g(x)\rho_n(x)\,dx \longrightarrow \frac{1}{L} \int_0^L g(x)\,dx.

Obtain an explicit O(1/n)O(1/n) bound on the oscillatory term.

Solution

At generic fixed xx, the phase 2nπx/L2n\pi x/L continues to wind as nn increases, so the cosine does not settle to a single value. Thus ρn(x)\rho_n(x) does not converge pointwise in general.

Let

kn=2nπL.k_n = \frac{2n\pi}{L}.

Let

In=∫0Lg(x)cos⁡(knx) dx.I_n = \int_0^L g(x)\cos(k_nx)\,dx.

Integration by parts gives

In=[g(x)sin⁡(knx)kn]0L−1kn∫0Lg′(x)sin⁡(knx) dx.I_n = \left[ \frac{g(x)\sin(k_nx)}{k_n} \right]_0^L -\frac{1}{k_n} \int_0^L g'(x)\sin(k_nx)\,dx.

Because sin⁡(knL)=sin⁡(2nπ)=0\sin(k_nL)=\sin(2n\pi)=0, the boundary term vanishes. Therefore

∣In∣≤1kn∫0L∣g′(x)∣ dx=O(1n).\begin{aligned} \lvert I_n\rvert &\le \frac{1}{k_n} \int_0^L\lvert g'(x)\rvert\,dx \\ &= O\left( \frac{1}{n} \right). \end{aligned}

The oscillatory contribution to the measured probability vanishes, leaving

lim⁡n→∞∫0Lg(x)ρn(x) dx=1L∫0Lg(x) dx.\lim_{n\to\infty} \int_0^L g(x)\rho_n(x)\,dx = \frac{1}{L} \int_0^L g(x)\,dx.

This proves weak convergence for the stated observable class without claiming pointwise convergence.

  • V. I. Arnold, Mathematical Methods of Classical Mechanics, 2nd ed. (Springer, 1989). Action-angle variables, invariant tori, and the classical structures used in EBK correspondence.
  • G. P. Berman and G. M. Zaslavsky, “Condition of stochasticity in quantum nonlinear systems”, Physica A 91, 450–460 (1978). Early analysis of logarithmic correspondence times in unstable dynamics.
  • M. V. Berry and K. E. Mount, “Semiclassical approximations in wave mechanics”, Reports on Progress in Physics 35, 315–397 (1972). Broad review of WKB, classical correspondence, turning points, and multidimensional semiclassics.
  • M. V. Berry, “Semi-Classical Mechanics in Phase Space: A Study of Wigner’s Function”, Philosophical Transactions of the Royal Society A 287, 237–271 (1977). Semiclassical phase-space measures, eigenstates, and classical invariant sets.
  • M. Brack and R. K. Bhaduri, Semiclassical Physics, revised ed. (Westview Press, 2003). Spectral densities, periodic orbits, and regular and chaotic semiclassical systems.
  • M. S. Child, Semiclassical Mechanics with Molecular Applications, 2nd ed. (Oxford University Press, 2014). Action quantization, correspondence of frequencies and matrix elements, and molecular applications.
  • M. C. Gutzwiller, Chaos in Classical and Quantum Mechanics (Springer, 1990). Classical instability, semiclassical propagation, periodic orbits, and quantum chaos.
  • L. D. Landau and E. M. Lifshitz, Quantum Mechanics: Non-Relativistic Theory, 3rd ed. (Butterworth-Heinemann, 1981). Standard WKB correspondence between highly excited bound states and classical motion.
  • R. G. Littlejohn, “The semiclassical evolution of wave packets”, Physics Reports 138, 193–291 (1986). Systematic treatment of localized states, classical transport, and long-time limitations.
  • W. H. Zurek, “Decoherence, einselection, and the quantum origins of the classical”, Reviews of Modern Physics 75, 715–775 (2003). Authoritative review separating environmental decoherence from closed-system classical limits.