Stationary Phase and the Classical Limit

The formal real-time path integral

$$K(q_f,t_f;q_i,t_i) = \int_{q(t_i)=q_i}^{q(t_f)=q_f} \mathcal Dq\, e^{iS[q]/\hbar}$$

is an oscillatory integral over histories. In a semiclassical regime, stationary phase organizes it around paths satisfying

$$\delta S[q_{\mathrm{cl}}]=0.$$

Those stationary paths obey the classical equations of motion. Quantum mechanics has not been replaced by one classical trajectory: Gaussian fluctuations determine the leading prefactor, several classical paths can interfere, and higher variations generate corrections.

This page is the canonical path-integral application of stationary phase. The general asymptotic theorem belongs to Stationary Phase, the phase intuition belongs to Action and Phase, and the full Van Vleck formula belongs to Semiclassical Propagator.

Oscillatory Integrals

Begin with an ordinary integral

$$I(\hbar) = \int d^nx\, a(x) \exp\left[ \frac{i}{\hbar}S(x) \right].$$

The useful small parameter is not dimensionful $\hbar$ by itself. Choose a characteristic action $S_0$ and define

$$\epsilon_{\mathrm{sc}} = \frac{\hbar}{S_0}.$$

When $\epsilon_{\mathrm{sc}}\ll1$, the phase varies rapidly except near stationary points

$$\nabla S(x_\star)=0.$$

If the Hessian at an isolated saddle is nondegenerate, the local quadratic approximation gives a Gaussian oscillatory integral. Nonstationary regions are suppressed by phase cancellation under suitable smoothness, endpoint, and decay assumptions. The exact multidimensional coefficient and Hessian-signature phase are derived in the canonical Stationary Phase page.

A path integral acquires meaning through a regulator such as time slicing. With $N$ intervals and $N-1$ intermediate coordinates, it is an ordinary finite-dimensional integral:

$$K_N = C_N \int dq_1\cdots dq_{N-1} \exp\left[ \frac{i}{\hbar}S_N(q_1,\ldots,q_{N-1}) \right].$$

Stationary phase should first be understood at this regulated level. The functional notation is the formal continuum summary of those finite-dimensional saddle equations and Hessians.

Stationary-Phase Approximation

Let $q_{\gamma,N}$ denote isolated stationary points of the discrete action:

$$\frac{\partial S_N}{\partial q_j} \bigg|_{q_{\gamma,N}} =0, \qquad j=1,\ldots,N-1.$$

Near one saddle, write

$$q_j = q_{\gamma,j}+\eta_j.$$

Then

$$S_N[q] = S_N[q_\gamma] + \frac12 \eta^{\mathsf T}H_{\gamma,N}\eta + O(\eta^3),$$

where

$$(H_{\gamma,N})_{jk} = \left. \frac{\partial^2S_N} {\partial q_j\partial q_k} \right|_{q_{\gamma,N}}.$$

The regulated saddle contribution has the structure

$$K_{N,\gamma} \sim C_N e^{iS_N[q_\gamma]/\hbar} \frac{ (2\pi\hbar)^{(N-1)/2} e^{i\pi\sigma_{\gamma,N}/4} }{ \sqrt{\lvert\det H_{\gamma,N}\rvert} },$$

up to convention-dependent factors of $i$ already carried by $C_N$. Here $\sigma_{\gamma,N}$ is the Hessian signature. The normalization and determinant separately depend on the discretization; only their consistently regulated combination has physical meaning.

If several saddles satisfy the endpoint data, their amplitudes must be added before taking an absolute square. A saddle-point approximation is therefore still an interference calculation.

Classical Equations of Motion

For

$$L(q,\dot q) = \frac{m}{2}\dot q^2 - V(q),$$

a simple discrete action is

$$S_N = \sum_{j=1}^{N} \left[ \frac{m}{2\Delta t} (q_j-q_{j-1})^2 - \Delta t\,V(q_j) \right],$$

with fixed $q_0=q_i$ and $q_N=q_f$. An endpoint or midpoint potential rule changes finite-step details but not the elementary continuum equation for this regular Cartesian example.

Stationarity with respect to an interior $q_j$ gives

$$m \frac{ q_{j+1}-2q_j+q_{j-1} }{(\Delta t)^2} + V'(q_j) =0.$$

As $\Delta t\to0$, this becomes

$$m\ddot q_{\mathrm{cl}}(t) + V'\bigl(q_{\mathrm{cl}}(t)\bigr) =0,$$

the Euler–Lagrange equation for the classical path joining the fixed endpoints.

The result depends on the boundary-value problem. A propagator fixes initial and final coordinates. A trace uses periodic paths. Coherent-state and phase-space path integrals impose different endpoint data. Stationary equations and boundary terms must be derived from the regulated object being approximated, not imported from a superficially similar integral.

Fluctuations Around Classical Paths

Write

$$q(t) = q_{\mathrm{cl}}(t)+\eta(t), \qquad \eta(t_i)=\eta(t_f)=0.$$

For the Lagrangian above,

$$\begin{aligned} S[q] &= S[q_{\mathrm{cl}}] + \frac12 \int_{t_i}^{t_f}dt\, \left[ m\dot\eta^2 - V''(q_{\mathrm{cl}})\eta^2 \right]\\ &\quad+ O(\eta^3). \end{aligned}$$

After integrating by parts, the quadratic form is

$$\delta^2S = \int_{t_i}^{t_f}dt\, \eta(t)\, \mathcal J_\gamma\eta(t),$$

with fixed-endpoint Jacobi operator

$$\mathcal J_\gamma = -m\frac{d^2}{dt^2} - V''\bigl(q_{\mathrm{cl}}(t)\bigr).$$

The formal Gaussian integral over $\eta$ supplies

$$\bigl(\det\mathcal J_\gamma\bigr)^{-1/2},$$

together with a phase fixed by the regulator and the signs of the fluctuation modes. The equivalent endpoint-stability expression is the Van Vleck Determinant.

Why powers of ℏ appear

Rescale the fluctuation by

$$\eta=\sqrt{\hbar}\,\xi.$$

Then the quadratic term in $S/\hbar$ is order unity, while a cubic variation contributes schematically as $\sqrt{\hbar}$ and a quartic variation as $\hbar$, after dimensionful scales have been factored out. This is the local origin of a semiclassical expansion around a nondegenerate saddle.

Zero modes require separate collective coordinates rather than an ordinary determinant. A negative mode changes the contour and phase structure. At a caustic, an eigenvalue passes through zero and the isolated Gaussian approximation fails; Maslov Index and uniform approximations repair the description.

Semiclassical Propagator Preview

After the regulator is removed consistently, the leading structure is

$$K_{\mathrm{sc}}(q_f,t_f;q_i,t_i) \sim \sum_{\gamma:i\to f} A_\gamma \exp\left[ \frac{i}{\hbar}S_\gamma - i\frac{\pi}{2}\nu_\gamma \right].$$

Each symbol carries information:

Quantity	Meaning	Canonical treatment
$S_\gamma$	classical action for a path joining the endpoints	Action and Phase
$A_\gamma$	normalization and classical stability	Van Vleck Determinant
$\nu_\gamma$	caustic or conjugate-point phase count	Maslov Index
sum over $\gamma$	interference among classical branches	Semiclassical Propagator

For quadratic actions, including the free particle and harmonic oscillator away from caustic times, the Gaussian expansion is exact after normalization. For a generic potential, the formula is asymptotic and higher variations generate corrections.

Limits of the Argument

Stationary phase explains one important classical-limit mechanism, not every emergence of classical behavior.

• No universal small parameter: one must identify dimensionless ratios such as $\hbar/S_0$ for the observable and time scale of interest.
• Noncommuting limits: the large-action, long-time, continuum, and regulator-removal limits need not commute.
• Caustics and degenerate saddles: ordinary Gaussian stationary phase fails when the fluctuation determinant vanishes.
• Endpoints and boundaries: leading contributions may come from boundaries even without an interior saddle.
• Tunneling: classically forbidden amplitudes often require Euclidean or complex saddles rather than real classical paths.
• Symmetries and zero modes: continuous families of saddles require collective-coordinate treatment.
• Long-time instability: in chaotic systems, semiclassical errors can grow and limit the useful propagation time.
• Decoherence is separate: stationary phase organizes amplitudes; it does not by itself turn a coherent superposition into an improper mixture.
• Measure and ordering matter: curvilinear coordinates, constrained systems, spin, and gauge theories require more structure than the elementary Cartesian formula.

What Is the Classical Limit? compares stationary phase with wave-packet localization, phase-space evolution, coarse graining, and decoherence. The Semiclassical Limit page owns the dimensionless scaling framework.

Common Mistakes

• Applying a formal functional saddle directly without identifying a regulator.
• Saying $\hbar$ is small without specifying a comparison action.
• Treating a stationary path as a probability maximum.
• Dropping Gaussian fluctuations and keeping only $e^{iS_{\mathrm{cl}}/\hbar}$.
• Assuming the finite-dimensional Hessian determinant has a regulator-independent meaning by itself.
• Keeping one classical path when several saddles have comparable magnitude.
• Using the isolated-saddle formula at a caustic or in the presence of zero modes.
• Concluding that nonclassical paths are absent from the exact integral.
• Equating stationary phase with environmental decoherence.
• Ignoring endpoint conditions, boundary saddles, or complex saddles.

Cross-Links

• Time Slicing
• Action and Phase
• Stationary Phase
• Stationary Phase in Quantum Mechanics
• What Is the Classical Limit?
• Semiclassical Propagator Preview
• Action Principles
• Semiclassical Limit
• Semiclassical Propagator
• Van Vleck Determinant
• Maslov Index
• Instantons in Quantum Mechanics Preview

References

• R. P. Feynman and A. R. Hibbs, Quantum Mechanics and Path Integrals, McGraw-Hill, 1965.
• L. S. Schulman, Techniques and Applications of Path Integration, Wiley, 1981.
• H. Kleinert, Path Integrals in Quantum Mechanics, Statistics, Polymer Physics, and Financial Markets, 5th ed., World Scientific, 2009.
• M. C. Gutzwiller, Chaos in Classical and Quantum Mechanics, Springer, 1990.
• R. G. Littlejohn, “The semiclassical evolution of wave packets,” Physics Reports 138, 193–291 (1986).
• J. H. Van Vleck, “The correspondence principle in the statistical interpretation of quantum mechanics,” Proceedings of the National Academy of Sciences 14, 178–188 (1928).

Exercises

1. Derive the discrete stationary equation for the time-sliced action of a particle in $V(q)$.

Solution

The terms depending on an interior $q_j$ are

$$\frac{m}{2\Delta t} (q_j-q_{j-1})^2 + \frac{m}{2\Delta t} (q_{j+1}-q_j)^2 - \Delta t\,V(q_j).$$

Differentiating and setting the result to zero gives

$$\frac{m}{\Delta t} \left( 2q_j-q_{j-1}-q_{j+1} \right) - \Delta t\,V'(q_j) =0.$$

Multiplying by $-1/\Delta t$ yields

$$m \frac{q_{j+1}-2q_j+q_{j-1}} {(\Delta t)^2} + V'(q_j) =0.$$

The central difference converges to $\ddot q$, producing $m\ddot q+V'(q)=0$.

2. Derive the fixed-endpoint Jacobi operator for $L=m\dot q^2/2-V(q)$.

Solution

Set $q=q_{\mathrm{cl}}+\eta$. The quadratic part of the action is

$$\frac12 \int dt\, \left[ m\dot\eta^2 - V''(q_{\mathrm{cl}})\eta^2 \right].$$

Since $\eta$ vanishes at both endpoints,

$$\int dt\,m\dot\eta^2 = -\int dt\, m\eta\ddot\eta.$$

Therefore the quadratic action is

$$\frac12 \int dt\, \eta \left[ -m\frac{d^2}{dt^2} - V''(q_{\mathrm{cl}}) \right] \eta,$$

so

$$\mathcal J = -m\frac{d^2}{dt^2} - V''(q_{\mathrm{cl}}).$$

3. Use $\eta=\sqrt{\hbar}\,\xi$ to power count the fluctuation expansion.

Solution

Write schematically

$$S = S_{\mathrm{cl}} + \frac12S^{(2)}\eta^2 + \frac{1}{3!}S^{(3)}\eta^3 + \frac{1}{4!}S^{(4)}\eta^4 + \cdots.$$

After dividing by $\hbar$ and substituting $\eta=\sqrt{\hbar}\,\xi$,

$$\frac{S}{\hbar} = \frac{S_{\mathrm{cl}}}{\hbar} + \frac12S^{(2)}\xi^2 + \frac{\sqrt{\hbar}}{3!}S^{(3)}\xi^3 + \frac{\hbar}{4!}S^{(4)}\xi^4 + \cdots,$$

after the necessary dimensionful scales are absorbed into the coefficients. The quadratic integral is leading, cubic corrections begin at order $\sqrt{\hbar}$ locally, and quartic corrections begin at order $\hbar$. Symmetry or Gaussian averaging can make odd contributions vanish in a final observable.

4. Find the fluctuation eigenvalues for a free particle with fixed endpoints over a duration $T$ and explain why there is no zero mode.

Solution

For $V''=0$,

$$\mathcal J = -m\frac{d^2}{dt^2}$$

on functions satisfying $\eta(0)=\eta(T)=0$. The normalized eigenfunctions are proportional to

$$\sin\left( \frac{n\pi t}{T} \right), \qquad n=1,2,\ldots,$$

with eigenvalues

$$\lambda_n = m\left( \frac{n\pi}{T} \right)^2.$$

Every eigenvalue is positive. The constant function would be a zero mode of $-d^2/dt^2$, but it violates the fixed-endpoint conditions unless it vanishes identically.

5. Why does stationary phase not by itself explain the disappearance of macroscopic interference?

Solution

Stationary phase reorganizes a coherent amplitude into contributions from saddles and their fluctuations. If several saddles contribute, their amplitudes still interfere. Environmental decoherence instead entangles the selected system with unobserved degrees of freedom and suppresses off-diagonal coherence in a reduced description. Coarse graining can also average rapidly varying fringes. These mechanisms can work together, but none is identical to the stationary-phase approximation.