Skip to content

Propagator Table

This table collects common real-time propagator kernels for nonrelativistic quantum mechanics. Each kernel is a coordinate-space transition amplitude for the time-evolution operator:

K(qf,tf;qi,ti)=⟨qf∣U(tf,ti)∣qi⟩.K(q_f,t_f;q_i,t_i) = \langle q_f\rvert U(t_f,t_i)\lvert q_i\rangle.

It is an amplitude, not a probability density. The formulas below use

T=tf−ti.T=t_f-t_i.

Unless stated otherwise, take T>0T\gt0 and use the usual real-time convergence prescription, equivalently T→T−i0+T\to T-i0^+ in Gaussian kernels. For other configuration spaces, the integration measure and boundary conditions are part of the kernel.

KernelConfiguration spaceBest useCanonical page
Free particlelineone-dimensional derivation, wave-packet spreadingFree-Particle Propagator
Free particleRd\mathbb R^dproduct structure, multidimensional short-time kernelPropagators in Multiple Dimensions
Harmonic oscillatorlinequadratic dynamics, coherent states, Gaussian path integralsHarmonic-Oscillator Propagator
Infinite square wellinterval 0<x<L0\lt x\lt Lhard-wall spectral expansionInfinite Square Well
Particle on a ringcircle S1S^1winding sectors, angular motionParticle on a Ring
Constant forcelineuniform-field motion, Airy-energy contrastthis table
Semiclassical Van Vleckgeneral smooth configuration spacelarge-action approximationSemiclassical Propagator Preview

The kernel evolves wavefunctions by

ψ(qf,tf)=∫dqi K(qf,tf;qi,ti)ψ(qi,ti),\psi(q_f,t_f) = \int dq_i\, K(q_f,t_f;q_i,t_i)\psi(q_i,t_i),

where dqidq_i denotes the correct configuration-space measure.

The equal-time limit is the identity kernel:

K(qf,ti;qi,ti)=δ(qf−qi),K(q_f,t_i;q_i,t_i) = \delta(q_f-q_i),

with the delta distribution understood relative to the same measure.

The composition law is

K(qf,tf;qi,ti)=∫dq K(qf,tf;q,t)K(q,t;qi,ti).K(q_f,t_f;q_i,t_i) = \int dq\, K(q_f,t_f;q,t)K(q,t;q_i,t_i).

For a time-independent Hamiltonian with a complete orthonormal energy basis, the spectral form is

K(qf,T;qi,0)=∑nψn(qf)ψn∗(qi)e−iEnT/ℏ,K(q_f,T;q_i,0) = \sum_n \psi_n(q_f)\psi_n^*(q_i) e^{-iE_nT/\hbar},

with integrals added for continuous spectra.

For

H=p22mH=\frac{p^2}{2m}

on the line,

K0(xf,T;xi,0)=(m2πiℏT)1/2exp⁡[im(xf−xi)22ℏT].K_0(x_f,T;x_i,0) = \left( \frac{m}{2\pi i\hbar T} \right)^{1/2} \exp\left[ \frac{im(x_f-x_i)^2}{2\hbar T} \right].

Use this kernel on the full line. On an interval, half-line, ring, or box, the same differential expression requires different boundary conditions and therefore a different kernel; see Propagators and Boundary Conditions.

For a free particle on Rd\mathbb R^d,

K0(d)(rf,T;ri,0)=(m2πiℏT)d/2exp⁡[im∣rf−ri∣22ℏT].K_0^{(d)}(\mathbf r_f,T;\mathbf r_i,0) = \left( \frac{m}{2\pi i\hbar T} \right)^{d/2} \exp\left[ \frac{im\lvert\mathbf r_f-\mathbf r_i\rvert^2} {2\hbar T} \right].

The kernel factorizes into one-dimensional free kernels in Cartesian coordinates. Propagators in Multiple Dimensions derives the formula and develops its configuration-space measure and product structure.

Let

H=p22m−Fx,H = \frac{p^2}{2m}-Fx,

so the classical force is +F+F. The exact real-time kernel on the line is

KF(xf,T;xi,0)=(m2πiℏT)1/2exp⁡[iℏSF],K_F(x_f,T;x_i,0) = \left( \frac{m}{2\pi i\hbar T} \right)^{1/2} \exp\left[ \frac{i}{\hbar}S_F \right],

where

SF=m(xf−xi)22T+FT2(xf+xi)−F2T324m.S_F = \frac{m(x_f-x_i)^2}{2T} + \frac{FT}{2}(x_f+x_i) - \frac{F^2T^3}{24m}.

Setting F=0F=0 recovers the free-particle kernel. The energy eigenfunctions of a linear potential are Airy functions, but the time-domain propagator remains a simple classical-action Gaussian because the Lagrangian is quadratic plus a linear term.

For a particle on 0<x<L0\lt x\lt L with hard-wall boundary conditions,

ψ(0)=ψ(L)=0,\psi(0)=\psi(L)=0,

the normalized eigenfunctions are

ψn(x)=2Lsin⁡nπxL,n=1,2,…,\psi_n(x) = \sqrt{\frac{2}{L}} \sin\frac{n\pi x}{L}, \qquad n=1,2,\ldots,

and

En=n2π2ℏ22mL2.E_n = \frac{n^2\pi^2\hbar^2}{2mL^2}.

The kernel is the spectral sum

Kbox(xf,T;xi,0)=∑n=1∞2Lsin⁡nπxfLsin⁡nπxiLexp⁡[−iEnTℏ].K_{\rm box}(x_f,T;x_i,0) = \sum_{n=1}^{\infty} \frac{2}{L} \sin\frac{n\pi x_f}{L} \sin\frac{n\pi x_i}{L} \exp\left[ -\frac{iE_nT}{\hbar} \right].

This is not the free-particle kernel restricted to 0<x<L0\lt x\lt L. The boundary conditions are part of the Hamiltonian domain and change the kernel.

For

H=p22m+12mω2x2,H = \frac{p^2}{2m} + \frac12m\omega^2x^2,

away from caustic times sin⁡ωT=0\sin\omega T=0,

Kho(xf,T;xi,0)=[mω2πiℏsin⁡ωT]1/2×exp⁡{imω2ℏsin⁡ωT[(xf2+xi2)cos⁡ωT−2xfxi]}.\begin{aligned} K_{\rm ho}(x_f,T;x_i,0) &= \left[ \frac{m\omega}{2\pi i\hbar\sin\omega T} \right]^{1/2} \\ &\quad\times \exp\left\{ \frac{im\omega}{2\hbar\sin\omega T} \left[ (x_f^2+x_i^2)\cos\omega T -2x_fx_i \right] \right\}. \end{aligned}

At

T=nπω,T=\frac{n\pi}{\omega},

the closed form must be replaced by the corresponding delta-function limit with Maslov phase. The oscillator kernel is exact because the action is quadratic.

Let a particle move on a ring of radius RR with moment of inertia

I=mR2.I=mR^2.

Using the angular coordinate θ∼θ+2π\theta\sim\theta+2\pi, one spectral form is

KS1(θf,T;θi,0)=12π∑n∈Zexp⁡[in(θf−θi)−iℏn2T2I].K_{S^1}(\theta_f,T;\theta_i,0) = \frac{1}{2\pi} \sum_{n\in\mathbb Z} \exp\left[ in(\theta_f-\theta_i) - \frac{i\hbar n^2T}{2I} \right].

An equivalent winding-sector form is

KS1(θf,T;θi,0)=∑w∈Z(I2πiℏT)1/2exp⁡[iI(θf−θi+2πw)22ℏT],K_{S^1}(\theta_f,T;\theta_i,0) = \sum_{w\in\mathbb Z} \left( \frac{I}{2\pi i\hbar T} \right)^{1/2} \exp\left[ \frac{iI(\theta_f-\theta_i+2\pi w)^2} {2\hbar T} \right],

with the same real-time prescription as the free particle. The winding sum makes the topology visible; the spectral sum makes angular-momentum quantization visible.

For a smooth dd-dimensional configuration-space problem, the leading semiclassical propagator is

Ksc(qf,tf;qi,ti)=∑γ:i→f(12πiℏ)d/2∣Dγ∣1/2exp⁡[iℏSγ−iπ2νγ].K_{\rm sc}(q_f,t_f;q_i,t_i) = \sum_{\gamma:i\to f} \left( \frac{1}{2\pi i\hbar} \right)^{d/2} \lvert D_\gamma\rvert^{1/2} \exp\left[ \frac{i}{\hbar}S_\gamma - i\frac{\pi}{2}\nu_\gamma \right].

Here

Dγ=det⁡(−∂2Sγ∂qf ∂qi).D_\gamma = \det\left( - \frac{\partial^2S_\gamma} {\partial q_f\,\partial q_i} \right).

This is not an exact formula for arbitrary potentials. It is the leading isolated-trajectory stationary-phase approximation. It must be repaired near caustics and supplemented by uniform approximations when classical branches coalesce.

Propagators and Boundary Conditions explains why the same differential expression can define different Hamiltonians on different domains:

DomainBoundary dataKernel consequence
full linesquare-integrabilitycontinuum free kernel
interval with hard wallsψ=0\psi=0 at endpointssine spectral sum
ringperiodicity in anglewinding or angular-momentum sum
half-lineself-adjoint boundary condition at endpointimage or spectral kernel depends on boundary condition

Always check the measure. A kernel on a ring evolves wavefunctions with dθd\theta, while a radial kernel in three dimensions uses the radial measure or a transformed radial wavefunction convention.

  • Treating propagators as probabilities rather than amplitudes.
  • Reusing the full-line free kernel in a box or on a ring without imposing boundary conditions.
  • Dropping the square-root phase in real-time Gaussian kernels.
  • Forgetting that equal-time limits are distributions.
  • Confusing time-domain propagator kernels with energy-domain Green functions.
  • Using closed-form oscillator kernels at caustic times without the delta-function and Maslov phase.
  • Ignoring the measure associated with angular, radial, or constrained coordinates.
  • 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.
  • R. Shankar, Principles of Quantum Mechanics, 2nd ed., Springer, 1994.
  • M. C. Gutzwiller, Chaos in Classical and Quantum Mechanics, Springer, 1990.
  1. Recover the one-dimensional free kernel from the dd-dimensional formula.
Solution

Set d=1d=1 and rf−ri=xf−xi\mathbf r_f-\mathbf r_i=x_f-x_i. Then

K0(1)=(m2πiℏT)1/2exp⁡[im(xf−xi)22ℏT],K_0^{(1)} = \left( \frac{m}{2\pi i\hbar T} \right)^{1/2} \exp\left[ \frac{im(x_f-x_i)^2}{2\hbar T} \right],

which is the one-dimensional free-particle kernel.

  1. Check that the constant-force kernel reduces to the free kernel when F=0F=0.
Solution

The constant-force action is

SF=m(xf−xi)22T+FT2(xf+xi)−F2T324m.S_F = \frac{m(x_f-x_i)^2}{2T} + \frac{FT}{2}(x_f+x_i) - \frac{F^2T^3}{24m}.

Setting F=0F=0 gives

SF→m(xf−xi)22T,S_F \to \frac{m(x_f-x_i)^2}{2T},

and the prefactor is already the free-particle prefactor. Thus KF→K0K_F\to K_0.

  1. Explain why the infinite-square-well kernel is not the full-line free kernel restricted to 0<x<L0\lt x\lt L.
Solution

The infinite square well imposes the hard-wall boundary conditions

ψ(0)=ψ(L)=0.\psi(0)=\psi(L)=0.

These boundary conditions change the domain of the Hamiltonian and force a sine eigenbasis. The full-line free kernel does not vanish at the walls and does not preserve the hard-wall domain. The correct kernel must be built from the well eigenfunctions or an equivalent image construction.

  1. Show that the ring spectral kernel is periodic in θf\theta_f.
Solution

The spectral kernel is

KS1=12π∑n∈Zein(θf−θi)e−iℏn2T/(2I).K_{S^1} = \frac{1}{2\pi} \sum_{n\in\mathbb Z} e^{in(\theta_f-\theta_i)} e^{-i\hbar n^2T/(2I)}.

Replacing θf\theta_f by θf+2π\theta_f+2\pi multiplies each term by

ei2πn=1e^{i2\pi n}=1

because nn is an integer. Therefore the kernel is periodic in θf\theta_f.

  1. Why does the Van Vleck formula need a caustic warning?
Solution

The Van Vleck prefactor contains

∣Dγ∣1/2,Dγ=det⁡(−∂2Sγ∂qf ∂qi).\lvert D_\gamma\rvert^{1/2}, \qquad D_\gamma = \det\left( - \frac{\partial^2S_\gamma} {\partial q_f\,\partial q_i} \right).

At a caustic, the classical map between endpoint data becomes singular, and the isolated-trajectory prefactor can diverge. The exact quantum kernel is not normally divergent there. The divergence signals that nearby branches must be treated together using a uniform approximation and the correct Maslov phase.