Skip to content

Free-Particle Propagator

The one-dimensional free-particle propagator is

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

for tf≠tit_f\ne t_i, with the square-root phase chosen by the usual real-time evolution prescription. It is the coordinate-space kernel of

U0(tf,ti)=exp⁡[−iℏp22m(tf−ti)].U_0(t_f,t_i) = \exp\left[ -\frac{i}{\hbar}\frac{p^2}{2m}(t_f-t_i) \right].

This page derives the kernel, checks its normalization and composition law, and explains why it is the basic short-time building block for path integrals.

Let

T=tf−ti,Δx=xf−xi.T=t_f-t_i, \qquad \Delta x=x_f-x_i.

Insert the momentum resolution of identity into

K0(xf,tf;xi,ti)=⟨xf∣U0(T)∣xi⟩.K_0(x_f,t_f;x_i,t_i) = \langle x_f|U_0(T)|x_i\rangle.

With

⟨x∣p⟩=12πℏeipx/ℏ,\langle x|p\rangle = \frac{1}{\sqrt{2\pi\hbar}}e^{ipx/\hbar},

one obtains

K0(xf,tf;xi,ti)=∫−∞∞dp2πℏ exp⁡[ipΔxℏ−ip2T2mℏ].\begin{aligned} K_0(x_f,t_f;x_i,t_i) &= \int_{-\infty}^{\infty} \frac{dp}{2\pi\hbar}\, \exp\left[ \frac{ip\Delta x}{\hbar} - \frac{ip^2T}{2m\hbar} \right]. \end{aligned}

This is a Fourier transform of the free-particle energy phase e−ip2T/(2mℏ)e^{-ip^2T/(2m\hbar)}.

Complete the square:

pΔxℏ−p2T2mℏ=−T2mℏ(p−mΔxT)2+m(Δx)22ℏT.\frac{p\Delta x}{\hbar} - \frac{p^2T}{2m\hbar} = - \frac{T}{2m\hbar} \left( p-\frac{m\Delta x}{T} \right)^2 + \frac{m(\Delta x)^2}{2\hbar T}.

Including the factor of ii in the exponent,

K0=eim(Δx)2/(2ℏT)∫−∞∞dp2πℏ exp⁡[−iT2mℏ(p−mΔxT)2].K_0 = e^{im(\Delta x)^2/(2\hbar T)} \int_{-\infty}^{\infty} \frac{dp}{2\pi\hbar}\, \exp\left[ -\frac{iT}{2m\hbar} \left( p-\frac{m\Delta x}{T} \right)^2 \right].

The remaining Fresnel integral is defined by a convergence prescription, for example T→T−i0+T\to T-i0^+ for T>0T\gt0. The result is

K0=(m2πiℏT)1/2exp⁡[im(Δx)22ℏT].K_0 = \left( \frac{m}{2\pi i\hbar T} \right)^{1/2} \exp\left[ \frac{im(\Delta x)^2}{2\hbar T} \right].

The prefactor is not optional: it enforces the delta-function initial condition and unitarity.

At equal times, the propagator must become the identity kernel:

lim⁡T→0K0(xf,ti+T;xi,ti)=δ(xf−xi)\lim_{T\to0} K_0(x_f,t_i+T;x_i,t_i) = \delta(x_f-x_i)

in the distributional sense. This does not mean the expression has an ordinary pointwise limit. Instead, for a sufficiently nice test wavefunction,

lim⁡T→0∫−∞∞dxi K0(xf,ti+T;xi,ti)ψ(xi)=ψ(xf).\lim_{T\to0} \int_{-\infty}^{\infty}dx_i\, K_0(x_f,t_i+T;x_i,t_i)\psi(x_i) = \psi(x_f).

The rapidly oscillating phase away from xi=xfx_i=x_f is what produces the delta distribution.

The kernel satisfies the free Schrödinger equation in the final coordinate:

iℏ∂K0∂tf=−ℏ22m∂2K0∂xf2.i\hbar\frac{\partial K_0}{\partial t_f} = -\frac{\hbar^2}{2m} \frac{\partial^2K_0}{\partial x_f^2}.

Together with the delta-function initial condition, this characterizes the free-particle kernel. The same expression also satisfies the adjoint equation in the initial variables, as required by unitarity.

For ti<t<tft_i\lt t\lt t_f, the kernel composes as

K0(xf,tf;xi,ti)=∫−∞∞dx K0(xf,tf;x,t)K0(x,t;xi,ti).K_0(x_f,t_f;x_i,t_i) = \int_{-\infty}^{\infty}dx\, K_0(x_f,t_f;x,t) K_0(x,t;x_i,t_i).

The cleanest verification uses the momentum representation. Each short kernel is an integral over a momentum phase; the intermediate position integral gives a delta distribution equating the two momenta. The result is the single momentum integral for the total time tf−tit_f-t_i.

This composition law is the starting point for the time-sliced construction in From Propagators to Path Integrals.

The classical path from (xi,ti)(x_i,t_i) to (xf,tf)(x_f,t_f) is a straight line with velocity

vcl=xf−xitf−ti.v_{\rm cl} = \frac{x_f-x_i}{t_f-t_i}.

The classical action is

Scl=∫titf12mvcl2 dt=m(xf−xi)22(tf−ti).S_{\rm cl} = \int_{t_i}^{t_f} \frac12m v_{\rm cl}^2\,dt = \frac{m(x_f-x_i)^2}{2(t_f-t_i)}.

Therefore

K0(xf,tf;xi,ti)=[12πiℏ(−∂2Scl∂xf ∂xi)]1/2eiScl/ℏ.K_0(x_f,t_f;x_i,t_i) = \left[ \frac{1}{2\pi i\hbar} \left( -\frac{\partial^2 S_{\rm cl}} {\partial x_f\,\partial x_i} \right) \right]^{1/2} e^{iS_{\rm cl}/\hbar}.

For the free particle this semiclassical-looking expression is exact because the action is quadratic.

The propagator evolves an initial wavefunction by

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

For an initial Gaussian packet

ψ(x,0)=1(2πσ02)1/4exp⁡[−(x−x0)24σ02+ip0(x−x0)ℏ],\psi(x,0) = \frac{1}{(2\pi\sigma_0^2)^{1/4}} \exp\left[ -\frac{(x-x_0)^2}{4\sigma_0^2} + \frac{ip_0(x-x_0)}{\hbar} \right],

the probability density remains Gaussian. Its center follows the classical trajectory

xc(t)=x0+p0mt,x_c(t)=x_0+\frac{p_0}{m}t,

while its width becomes

σt=σ0[1+(ℏt2mσ02)2]1/2.\sigma_t = \sigma_0 \left[ 1+ \left( \frac{\hbar t}{2m\sigma_0^2} \right)^2 \right]^{1/2}.

The packet spreads because the free-particle dispersion relation ω=ℏk2/(2m)\omega=\hbar k^2/(2m) makes different momentum components accumulate different phases. For the canonical wave-packet discussion, see Gaussian Wave Packets and Wave-Packet Spreading.

For a short time step ϵ\epsilon, the free-particle kernel is

K0(xj+1,tj+ϵ;xj,tj)=(m2πiℏϵ)1/2exp⁡[im(xj+1−xj)22ℏϵ].K_0(x_{j+1},t_j+\epsilon;x_j,t_j) = \left( \frac{m}{2\pi i\hbar\epsilon} \right)^{1/2} \exp\left[ \frac{im(x_{j+1}-x_j)^2}{2\hbar\epsilon} \right].

When a potential is present, a short-time kinetic-potential splitting gives the approximate kernel

K(xj+1,tj+ϵ;xj,tj)≈K0(xj+1,tj+ϵ;xj,tj)e−iϵV(xj)/ℏ.K(x_{j+1},t_j+\epsilon;x_j,t_j) \approx K_0(x_{j+1},t_j+\epsilon;x_j,t_j) e^{-i\epsilon V(x_j)/\hbar}.

Multiplying many such short-time kernels and integrating over intermediate positions is the formal route to the real-time path integral.

  • Treating K0K_0 as a probability density. It is an amplitude and is not normalizable as a function of xfx_f for fixed xix_i.
  • Dropping the phase of the square-root prefactor.
  • Forgetting that the equal-time limit is a distribution, not an ordinary function.
  • Reusing the free kernel in a box, on a half-line, or on a ring without enforcing the correct boundary conditions.
  • Confusing the free-particle propagator with a relativistic or QFT propagator.
  • Ignoring Fourier-transform conventions when deriving the prefactor.
  • 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.
  • R. Shankar, Principles of Quantum Mechanics, 2nd ed., Springer, 1994.
  • J. J. Sakurai and J. Napolitano, Modern Quantum Mechanics, 3rd ed., Cambridge University Press, 2020.
  • H. Kleinert, Path Integrals in Quantum Mechanics, Statistics, Polymer Physics, and Financial Markets, 5th ed., World Scientific, 2009.
  1. Starting from the momentum integral for K0K_0, complete the square and recover the stated prefactor.
Solution

Start with

K0=∫−∞∞dp2πℏexp⁡[ipΔxℏ−ip2T2mℏ].K_0 = \int_{-\infty}^{\infty} \frac{dp}{2\pi\hbar} \exp\left[ \frac{ip\Delta x}{\hbar} - \frac{ip^2T}{2m\hbar} \right].

Completing the square gives

ipΔxℏ−ip2T2mℏ=im(Δx)22ℏT−iT2mℏ(p−mΔxT)2.\frac{ip\Delta x}{\hbar} - \frac{ip^2T}{2m\hbar} = \frac{im(\Delta x)^2}{2\hbar T} - \frac{iT}{2m\hbar} \left( p-\frac{m\Delta x}{T} \right)^2.

The shifted Fresnel integral is

∫−∞∞dq2πℏexp⁡(−iTq22mℏ)=(m2πiℏT)1/2,\int_{-\infty}^{\infty} \frac{dq}{2\pi\hbar} \exp\left( -\frac{iTq^2}{2m\hbar} \right) = \left( \frac{m}{2\pi i\hbar T} \right)^{1/2},

with the standard convergence prescription. Combining factors gives the result.

  1. Verify the free Schrödinger equation for the kernel by differentiating with respect to TT and Δx\Delta x.
Solution

Write

K0=CT−1/2exp⁡[im(Δx)22ℏT].K_0=C T^{-1/2} \exp\left[ \frac{im(\Delta x)^2}{2\hbar T} \right].

Then

∂K0∂T=[−12T−im(Δx)22ℏT2]K0,\frac{\partial K_0}{\partial T} = \left[ -\frac{1}{2T} - \frac{im(\Delta x)^2}{2\hbar T^2} \right]K_0,

so

iℏ∂K0∂T=[−iℏ2T+m(Δx)22T2]K0.i\hbar\frac{\partial K_0}{\partial T} = \left[ -\frac{i\hbar}{2T} + \frac{m(\Delta x)^2}{2T^2} \right]K_0.

Also,

∂2K0∂xf2=[imℏT−m2(Δx)2ℏ2T2]K0.\frac{\partial^2K_0}{\partial x_f^2} = \left[ \frac{im}{\hbar T} - \frac{m^2(\Delta x)^2}{\hbar^2T^2} \right]K_0.

Multiplying by −ℏ2/(2m)-\hbar^2/(2m) gives the same expression as iℏ ∂TK0i\hbar\,\partial_TK_0.

  1. Use the momentum representation to verify the composition law for the free-particle propagator.
Solution

Write each kernel as a momentum integral:

K0(xf,tf;x,t)=∫dp2πℏeip(xf−x)/ℏe−ip2(tf−t)/(2mℏ)K_0(x_f,t_f;x,t) = \int\frac{dp}{2\pi\hbar} e^{ip(x_f-x)/\hbar} e^{-ip^2(t_f-t)/(2m\hbar)}

and similarly for K0(x,t;xi,ti)K_0(x,t;x_i,t_i). The intermediate integral over xx gives

∫dx ei(q−p)x/ℏ=2πℏ δ(q−p).\int dx\,e^{i(q-p)x/\hbar} = 2\pi\hbar\,\delta(q-p).

The delta distribution sets q=pq=p, leaving

∫dp2πℏeip(xf−xi)/ℏe−ip2(tf−ti)/(2mℏ),\int\frac{dp}{2\pi\hbar} e^{ip(x_f-x_i)/\hbar} e^{-ip^2(t_f-t_i)/(2m\hbar)},

which is the single free kernel for the total time.

  1. Why does the free-particle kernel not decay at large ∣xf−xi∣|x_f-x_i| for fixed nonzero TT?
Solution

The magnitude of the real-time free kernel is set by the prefactor:

∣K0∣=(m2πℏ∣T∣)1/2.\left|K_0\right| = \left( \frac{m}{2\pi\hbar |T|} \right)^{1/2}.

The dependence on xf−xix_f-x_i is purely in the phase. This is not a probability distribution over endpoints. Physical probabilities arise after convolving the kernel with a normalizable initial wavefunction and then taking the modulus squared.