Skip to content

Propagator Kernel

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 a transition amplitude, not a probability.

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).

This says that no time evolution maps a position amplitude to itself by the identity kernel.

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 real line and T=tf−ti>0T=t_f-t_i>0, the momentum resolution of the identity gives

K0(xf,T;xi,0)=∫−∞∞dp2πℏexp⁡ ⁣[iℏp(xf−xi)−iTℏp22m].K_0(x_f,T;x_i,0) = \int_{-\infty}^{\infty}\frac{dp}{2\pi\hbar} \exp\!\left[ \frac{i}{\hbar}p(x_f-x_i)-\frac{iT}{\hbar}\frac{p^2}{2m} \right].

The oscillatory Gaussian is defined by the causal convergence prescription T↦T−i0+T\mapsto T-i0^+, or equivalently by first evaluating it with a small damping factor and then taking the limit. The result is

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 branch is fixed by continuity from the damped integral and by the distributional condition K0→δ(xf−xi)K_0\to\delta(x_f-x_i) as T→0+T\to0^+. Its phase is the classical action,

Scl=m(xf−xi)22T,S_{\mathrm{cl}} = \frac{m(x_f-x_i)^2}{2T},

so stationary phase selects the classical momentum pcl=m(xf−xi)/Tp_{\mathrm{cl}}=m(x_f-x_i)/T. In dd Cartesian dimensions the independent Gaussian integrals give

K0(xf,T;xi,0)=(m2πiℏT)d/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)^{d/2} \exp\!\left[ \frac{im\lvert\mathbf x_f-\mathbf x_i\rvert^2}{2\hbar T} \right].

Convolving this kernel with an initial packet displays free spreading directly. The Free-Particle Propagator: First Encounter keeps the coordinate-space interpretation close to the canonical system; this page owns the derivation and general kernel structure.

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. Shankar, Principles of Quantum Mechanics, 2nd ed., Springer, 1994.
  • L. S. Schulman, Techniques and Applications of Path Integration, Wiley, 1981.
  • 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.
  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. Apply stationary phase to the momentum integral for K0K_0 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].