Skip to content

The propagator kernel is the coordinate-space matrix element of the time-evolution operator:

K(xf,tf;xi,ti)=⟨xf∣U(tf,ti)∣xi⟩.K(x_f,t_f;x_i,t_i) = \langle x_f\rvert U(t_f,t_i)\lvert x_i\rangle.

It is an amplitude density between delta-normalized positions, with units of inverse length on the line. It becomes a physical position wavefunction only after integration against a normalizable initial state; see transition amplitudes.

Given an initial wavefunction ψ(xi,ti)\psi(x_i,t_i), the later wavefunction is

ψ(xf,tf)=∫−∞∞K(xf,tf;xi,ti)ψ(xi,ti) dxi.\psi(x_f,t_f) = \int_{-\infty}^{\infty} K(x_f,t_f;x_i,t_i)\psi(x_i,t_i)\,dx_i.

For other configuration spaces, the integration domain and measure must match the Hilbert-space inner product.

At equal times, the propagator reduces to a delta distribution:

K(xf,ti;xi,ti)=δ(xf−xi).K(x_f,t_i;x_i,t_i) =\delta(x_f-x_i).

With zero elapsed time, the identity kernel leaves every initial wavefunction unchanged.

As a function of the final variables, the kernel satisfies the time-dependent Schrödinger equation:

iℏ∂∂tfK(xf,tf;xi,ti)=HxfK(xf,tf;xi,ti),i\hbar\frac{\partial}{\partial t_f} K(x_f,t_f;x_i,t_i) = H_{x_f} K(x_f,t_f;x_i,t_i),

with the delta-function initial condition above. Here HxfH_{x_f} means that the coordinate-space Hamiltonian acts on the final coordinate.

Time-evolution operators compose:

U(tb,ta)=U(tb,tc)U(tc,ta).U(t_b,t_a) =U(t_b,t_c)U(t_c,t_a).

In coordinate representation this becomes

K(xb,tb;xa,ta)=∫dxc K(xb,tb;xc,tc)K(xc,tc;xa,ta).K(x_b,t_b;x_a,t_a) = \int dx_c\, K(x_b,t_b;x_c,t_c) K(x_c,t_c;x_a,t_a).

This identity is the Composition Law and is the direct route to path-integral time slicing.

For a time-independent Hamiltonian with discrete eigenstates,

K(xf,t;xi,0)=∑nψn(xf)ψn∗(xi)e−iEnt/ℏ.K(x_f,t;x_i,0) = \sum_n \psi_n(x_f)\psi_n^*(x_i) e^{-iE_nt/\hbar}.

Continuous spectra add integrals over continuum labels. Spectral Decomposition of the Propagator develops the discrete, continuous, and mixed cases and makes clear how bound states and scattering states contribute to propagation.

For H=p2/(2m)H=p^2/(2m) on the line and T>0T>0,

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

The square root has phase e−iπ/4e^{-i\pi/4}, fixed by the damped momentum integral and the delta-function limit. The exponent is iScl/ℏiS_{\mathrm{cl}}/\hbar, with Scl=m(xf−xi)2/(2T)S_{\mathrm{cl}}=m(x_f-x_i)^2/(2T). The free-particle propagator derives this expression and checks composition. Its constant magnitude as a function of xf−xix_f-x_i is another reason not to interpret ∣K0∣2|K_0|^2 as a normalized probability density.

The path integral begins by applying the composition law many times, inserting intermediate positions, and taking a formal continuum limit. The result is a sum over histories weighted by eiS/ℏe^{iS/\hbar}.

Propagators evolve states in time. Green functions often solve inhomogeneous differential equations or resolvent equations with specified boundary conditions. The two are related, but not identical; boundary conditions and time-ordering conventions matter.

  • Treating KK as a probability instead of an amplitude.
  • Forgetting the integration over the initial coordinate when evolving a wavefunction.
  • Ignoring the measure and boundary conditions of the configuration space.
  • Confusing the propagator kernel with an energy-domain Green function.
  • Forgetting that the equal-time limit is a distribution.
  • Quoting the free-particle square root without its convergence and branch prescription.
  • Using the infinite-line kernel when walls, periodicity, or another configuration-space boundary condition changes the kernel.
  • R. P. Feynman and A. R. Hibbs, Quantum Mechanics and Path Integrals, McGraw-Hill, 1965.
  • D. J. Griffiths and D. F. Schroeter, Introduction to Quantum Mechanics, 3rd ed., Cambridge University Press, 2018.
  • L. S. Schulman, Techniques and Applications of Path Integration, Wiley, 1981.
  • R. Shankar, Principles of Quantum Mechanics, 2nd ed., Springer, 1994.
  1. Show that the composition law for KK follows from inserting a position resolution of identity at the intermediate time.
Solution

Start from

K(xb,tb;xa,ta)=⟨xb∣U(tb,tc)U(tc,ta)∣xa⟩.K(x_b,t_b;x_a,t_a) = \langle x_b\rvert U(t_b,t_c)U(t_c,t_a)\lvert x_a\rangle.

Insert

I=∫dxc ∣xc⟩⟨xc∣.I=\int dx_c\,\lvert x_c\rangle\langle x_c\rvert.

This gives

K(xb,tb;xa,ta)=∫dxc K(xb,tb;xc,tc)K(xc,tc;xa,ta).K(x_b,t_b;x_a,t_a) = \int dx_c\, K(x_b,t_b;x_c,t_c) K(x_c,t_c;x_a,t_a).
  1. Use dimensional analysis to verify that the one-dimensional free kernel has units of inverse length.
Solution

Because ℏT/m\hbar T/m has units of length squared, the prefactor m/(ℏT)\sqrt{m/(\hbar T)} has units of inverse length. The exponential is dimensionless. This is also required by ψ(xf,T)=∫dxi Kψ(xi,0)\psi(x_f,T)=\int dx_i\,K\psi(x_i,0).

  1. For the free kernel, apply stationary phase to ∫dp ei[p(xf−xi)−p2T/(2m)]/ℏ\int dp\,e^{i[p(x_f-x_i)-p^2T/(2m)]/\hbar} and identify the stationary momentum.
Solution

The phase is Φ(p)=p(xf−xi)−p2T/(2m)\Phi(p)=p(x_f-x_i)-p^2T/(2m). Solving dΦ/dp=0d\Phi/dp=0 gives

pcl=m(xf−xi)T,p_{\mathrm{cl}}=\frac{m(x_f-x_i)}{T},

the momentum of the classical straight path joining the endpoints in time TT.

  1. Derive the three-dimensional free kernel by factorizing the Cartesian momentum integrals.
Solution

Since p2=px2+py2+pz2p^2=p_x^2+p_y^2+p_z^2 and the measure factorizes, the integral is the product of three one-dimensional kernels. Therefore

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