Path Integral Formulation

⚡ Key Takeaways

The path integral formulation of quantum mechanics expresses transition amplitudes as a sum over all possible paths connecting initial and final configurations, each weighted by a phase $e^{iS/\hbar}$ where $S$ is the classical action. Formulated by Richard Feynman (1948), it is equivalent to canonical (operator) quantum mechanics but reorganizes the theory around the Lagrangian and trajectories rather than Hamiltonians and operators. It is the foundation of modern quantum field theory, statistical mechanics, and many areas of mathematical physics.

💡 The Core Insight

Classical mechanics selects the single path of stationary action via $\delta S = 0$ . Quantum mechanics keeps all paths but assigns each a complex phase of unit modulus, $e^{iS/\hbar}$ . Constructive interference near the classical trajectory (where $S$ is stationary, so phases align) reproduces classical behavior; destructive interference elsewhere suppresses non-classical paths. The amplitude is the coherent superposition of every conceivable history.

DefinitionPropagator / Kernel

The central object is the propagator (or kernel), the position-space matrix element of the time-evolution operator:

1

\quad K(q_b,t_b;q_a,t_a) = \langle q_b | e^{-i\hat{H}(t_b-t_a)/\hbar} | q_a \rangle

Derivation hint: Insert the time-evolution operator $\hat{U}(t_b,t_a)=e^{-i\hat H(t_b-t_a)/\hbar}$ (for time-independent $\hat H$ ) between position eigenstates $\hat q|q\rangle=q|q\rangle$ . Here $\hat H$ is the Hamiltonian operator, $|q\rangle$ are position eigenstates with $\langle q'|q\rangle = \delta(q'-q)$ , and $K$ is the amplitude that propagates the wavefunction: $\psi(q_b,t_b) = \int dq_a\, K(q_b,t_b;q_a,t_a)\psi(q_a,t_a)$ .

DefinitionFeynman Path Integral

The propagator is expressed as a functional integral over all paths $q(t)$ with $q(t_a)=q_a$ , $q(t_b)=q_b$ :

2

\quad K(q_b,t_b;q_a,t_a) = \int_{q(t_a)=q_a}^{q(t_b)=q_b} \mathcal{D}q(t)\; \exp\!\left(\frac{i}{\hbar} S[q]\right)

Derivation hint: Obtained by time-slicing the propagator into $N$ infinitesimal evolutions and taking $N\to\infty$ (see proof below). where the action functional is the time integral of the Lagrangian:

3

\quad S[q] = \int_{t_a}^{t_b} dt\; L\big(q(t),\dot q(t),t\big), \qquad L = \tfrac{1}{2}m\dot q^2 - V(q)

Derivation hint: $L$ is the classical Lagrangian; for a non-relativistic particle in potential $V$ , kinetic minus potential energy. The symbol $\mathcal{D}q(t)$ denotes the (formal) integration measure over path space, defined as a limit of finite-dimensional integrals.

TheoremTime-Slicing Definition of the Measure

The path-integral measure is rigorously defined (in the discretized sense) by

4

\quad \int \mathcal{D}q\; e^{iS/\hbar} = \lim_{N\to\infty} \left(\frac{m}{2\pi i\hbar\,\epsilon}\right)^{N/2} \int \prod_{k=1}^{N-1} dq_k \; \exp\!\left(\frac{i}{\hbar}\sum_{k=0}^{N-1} \left[\frac{m(q_{k+1}-q_k)^2}{2\epsilon} - \epsilon V(q_k)\right]\right)

Derivation hint: Split $[t_a,t_b]$ into $N$ slices of width $\epsilon=(t_b-t_a)/N$ ; insert $N-1$ completeness relations $\int dq_k|q_k\rangle\langle q_k|=\mathbb{1}$ ; evaluate each short-time kernel. with $q_0=q_a$ , $q_N=q_b$ , $\epsilon = (t_b-t_a)/N$ , and the normalization factor $(m/2\pi i\hbar\epsilon)^{N/2}$ absorbed into $\mathcal{D}q$ .

ProofTime-Slicing Derivation

We derive [2]–[4] from [1] for $\hat H = \hat p^2/2m + V(\hat q)$ .

Step 1 — Slice the evolution. Write $t_b - t_a = N\epsilon$ and split the exponential:

5

\quad e^{-i\hat H(t_b-t_a)/\hbar} = \left(e^{-i\hat H\epsilon/\hbar}\right)^N

Derivation hint: The operator exponential factorizes because $\hat H$ commutes with itself; $N$ identical factors.

Step 2 — Insert completeness. Between each factor insert $\int dq_k |q_k\rangle\langle q_k| = \mathbb{1}$ , giving

6

\quad K = \int \prod_{k=1}^{N-1} dq_k \prod_{k=0}^{N-1}\langle q_{k+1}|e^{-i\hat H\epsilon/\hbar}|q_k\rangle

Derivation hint: $N-1$ resolutions of identity produce $N-1$ integrals and a product of $N$ short-time matrix elements.

Step 3 — Trotter factorization. For small $\epsilon$ , the Baker–Campbell–Hausdorff error in splitting $e^{-i\epsilon(\hat T+\hat V)/\hbar}\approx e^{-i\epsilon \hat T/\hbar}e^{-i\epsilon \hat V/\hbar}$ is $O(\epsilon^2)$ :

7

\quad \langle q_{k+1}|e^{-i\hat H\epsilon/\hbar}|q_k\rangle \approx \langle q_{k+1}|e^{-i\epsilon \hat p^2/2m\hbar}|q_k\rangle\, e^{-i\epsilon V(q_k)/\hbar}

Derivation hint: $\hat V(\hat q)$ is diagonal in position; act it on $|q_k\rangle$ to get the c-number $V(q_k)$ . The Trotter product formula guarantees convergence as $N\to\infty$ .

Step 4 — Evaluate the kinetic factor. Insert momentum completeness $\int \frac{dp}{2\pi\hbar}|p\rangle\langle p|$ and use $\langle q|p\rangle=e^{ipq/\hbar}$ :

8

\quad \langle q_{k+1}|e^{-i\epsilon \hat p^2/2m\hbar}|q_k\rangle = \int\frac{dp_k}{2\pi\hbar}\exp\!\left[\frac{i}{\hbar}\left(p_k(q_{k+1}-q_k)-\frac{\epsilon p_k^2}{2m}\right)\right]

Derivation hint: A Gaussian integral in $p_k$ ; complete the square in the exponent.

Step 5 — Do the Gaussian momentum integral.

9

\quad \int\frac{dp_k}{2\pi\hbar}\, e^{\frac{i}{\hbar}(p_k\Delta q_k - \epsilon p_k^2/2m)} = \left(\frac{m}{2\pi i\hbar\epsilon}\right)^{1/2}\exp\!\left[\frac{im(\Delta q_k)^2}{2\hbar\epsilon}\right]

Derivation hint: Standard Fresnel/Gaussian integral $\int dp\, e^{-ap^2+bp}=\sqrt{\pi/a}\,e^{b^2/4a}$ with imaginary $a$ ; here $\Delta q_k=q_{k+1}-q_k$ .

Step 6 — Assemble. Combining the $N$ factors and recognizing the Riemann sum:

10

\quad \sum_{k=0}^{N-1}\left[\frac{m}{2}\left(\frac{\Delta q_k}{\epsilon}\right)^2 - V(q_k)\right]\epsilon \;\xrightarrow{N\to\infty}\; \int_{t_a}^{t_b} dt\, L = S[q]

Derivation hint: $\Delta q_k/\epsilon\to\dot q$ , $\sum_k \epsilon \to \int dt$ ; the exponent becomes $iS/\hbar$ , completing [2] and [4]. $\blacksquare$

∎

⚠️ Common mistake

The "limit" in [4] is not the limit of integrals against a genuine measure. The factor $(m/2\pi i\hbar\epsilon)^{N/2}$ diverges as $N\to\infty$ , and the oscillatory integrand is not absolutely convergent. The path integral exists as a limit of finite-dimensional oscillatory integrals (or via analytic continuation from the Euclidean theory), not as a Lebesgue integral over path space. There is no translation-invariant $\sigma$ -additive measure on the infinite-dimensional path space — see Wiener Measure for the rigorous Euclidean counterpart.

RP

Richard P. Feynman

1918–1988

Feynman developed the path integral in his 1942 Princeton PhD thesis under John Wheeler and published the definitive account in Reviews of Modern Physics (1948), building on a cryptic remark by Dirac (1933) that the short-time propagator "corresponds to" $e^{iL\epsilon/\hbar}$ . Feynman turned "corresponds to" into "equals" (up to normalization), reformulating all of quantum mechanics as a sum over histories. He later extended it to electrodynamics, producing the Feynman diagrams (1949) that organize perturbation theory. He shared the 1965 Nobel Prize with Schwinger and Tomonaga. Feynman's path integral remains his deepest single contribution to theoretical physics, prized for making quantum mechanics visualizable as interfering histories.

🧠 Intuition

Imagine sending the particle from $A$ to $B$ and asking it to "try" every conceivable route — straight lines, loops, paths that briefly run backward, infinitely jagged zigzags. Each route carries a stopwatch hand (a unit complex phasor) that rotates at rate $S/\hbar$ . To get the total amplitude, add all the little arrows tip-to-tail. Where $S$ changes rapidly between neighboring paths (far from the classical trajectory), the arrows spin wildly and cancel. Near the classical path, $\delta S = 0$ means neighboring paths have nearly equal $S$ , so their arrows point the same way and add coherently. This is why classical mechanics ( $\hbar\to 0$ ) emerges: the classical path is where the phasors stop canceling. The width of the "corridor" of contributing paths is set by $\hbar$ — vanishingly thin for macroscopic objects, broad for electrons.

Figure

Path interference: the classical trajectory is the stationary-phase corridor among all summed paths

ExampleFree Particle Propagator

Given: A free particle, $V=0$ , mass $m$ , propagating from $q_a$ at $t_a=0$ to $q_b$ at time $T$ .

Step 1 — Classical action. The classical path is the straight line $q_{cl}(t)=q_a+(q_b-q_a)t/T$ with constant velocity $v=(q_b-q_a)/T$ . Then

11

\quad S_{cl} = \int_0^T \tfrac{1}{2}m v^2\, dt = \frac{m(q_b-q_a)^2}{2T}

Derivation hint: Constant velocity gives constant Lagrangian; multiply by $T$ .

Step 2 — Fluctuation integral. Write $q=q_{cl}+\eta$ with $\eta(0)=\eta(T)=0$ . Since $S$ is quadratic, $S=S_{cl}+\frac{1}{2}\int \dot\eta^2 m\,dt$ and the fluctuation determinant is field-independent, yielding the prefactor $(m/2\pi i\hbar T)^{1/2}$ (the continuum limit of the slicing normalization).

Step 3 — Result.

12

\quad K(q_b,T;q_a,0) = \left(\frac{m}{2\pi i\hbar T}\right)^{1/2} \exp\!\left[\frac{im(q_b-q_a)^2}{2\hbar T}\right]

Derivation hint: Combine $S_{cl}$ from [11] with the Gaussian fluctuation prefactor; verify it satisfies the free Schrödinger equation and reduces to $\delta(q_b-q_a)$ as $T\to0$ .

Numerical check (units). For an electron ( $m=9.11\times10^{-31}$ kg), $q_b-q_a=1$ nm $=10^{-9}$ m, $T=10^{-15}$ s, $\hbar=1.055\times10^{-34}$ J·s: the phase $=\dfrac{m(q_b-q_a)^2}{2\hbar T}=\dfrac{9.11\times10^{-31}\cdot 10^{-18}}{2\cdot1.055\times10^{-34}\cdot10^{-15}}\approx 4.3\times10^{0}\ \text{rad}\approx \boxed{4.3\ \text{rad}}$ (dimensionless).

ExampleHarmonic Oscillator Propagator

Given: $L=\tfrac{1}{2}m\dot q^2-\tfrac{1}{2}m\omega^2 q^2$ , propagation over time $T$ (with $\omega T\neq n\pi$ ).

Step 1 — Classical action. Solving $\ddot q_{cl}+\omega^2 q_{cl}=0$ with boundary values and substituting gives the well-known result:

13

\quad S_{cl} = \frac{m\omega}{2\sin\omega T}\Big[(q_a^2+q_b^2)\cos\omega T - 2q_a q_b\Big]

Derivation hint: Use $q_{cl}(t)=A\cos\omega t+B\sin\omega t$ , integrate by parts so $S_{cl}=\tfrac{m}{2}[q\dot q]_0^T$ , evaluate boundary terms.

Step 2 — Fluctuation determinant. The quadratic fluctuation operator $-m(\partial_t^2+\omega^2)$ gives, via the Gelfand–Yaglom theorem, the prefactor $\big(m\omega/2\pi i\hbar\sin\omega T\big)^{1/2}$ .

Step 3 — Result.

14

\quad K = \left(\frac{m\omega}{2\pi i\hbar\sin\omega T}\right)^{1/2}\exp\!\left(\frac{i}{\hbar}S_{cl}\right)

Derivation hint: Reduces to the free result [12] as $\omega\to0$ ( $\sin\omega T\to\omega T$ ); poles at $\omega T=n\pi$ encode the energy spectrum via the Mehler/spectral expansion.

Numerical check. Take $\omega=10^{14}$ s⁻¹, $T=\pi/(4\omega)=7.85\times10^{-15}$ s, $q_a=q_b=0$ . Then $S_{cl}=0$ , and the amplitude magnitude is $|K|=\big(m\omega/2\pi\hbar\sin\omega T\big)^{1/2}$ . With $m=9.11\times10^{-31}$ kg: $m\omega/(2\pi\hbar\sin(\pi/4)) = (9.11\times10^{-17})/(2\pi\cdot1.055\times10^{-34}\cdot0.707)=\boxed{1.94\times10^{17}\ \text{m}^{-1}}$ , so $|K|\approx 4.4\times10^{8}\ \text{m}^{-1/2}$ .

⚠️ Common mistake

The "measure" is not a measure. Students often treat $\mathcal{D}q$ like a Lebesgue measure and manipulate the path integral with the same freedom as a finite integral. In Minkowski (real-time) signature, $e^{iS/\hbar}$ has unit modulus everywhere, so the integral does not converge absolutely and no genuine measure exists. Rigorous statements require either the discretized limit [4], analytic continuation to imaginary time (Euclidean Wiener Measure), or oscillatory-integral (Fresnel) machinery. Formal manipulations frequently give correct answers, but their validity is non-trivial.

⚠️ Common mistake

Operator ordering hides in the measure. The naive replacement $\hat H\to H(q,p)$ is ambiguous when $\hat q$ and $\hat p$ do not commute (e.g. for position-dependent mass or magnetic terms). Different discretizations (midpoint vs endpoint evaluation of $V$ or of velocity-dependent terms) correspond to different operator orderings (Weyl, normal, etc.). The midpoint prescription corresponds to Weyl ordering. A "Lagrangian" path integral silently encodes an ordering choice.

Counterexample

Naively integrating the kinetic phase gives the wrong normalization without the prefactor. Suppose one writes $K=\int\mathcal{D}q\,e^{iS/\hbar}$ and tries to evaluate the free case without tracking the $(m/2\pi i\hbar\epsilon)^{N/2}$ normalization. One then obtains a propagator missing the crucial $(m/2\pi i\hbar T)^{1/2}$ factor, so $K$ fails the composition law $\int dq_c\, K(b,c)K(c,a)=K(b,a)$ and does not reduce to $\delta(q_b-q_a)$ as $T\to0$ . The normalization is not optional decoration — it is fixed by unitarity and is part of the definition of $\mathcal{D}q$ .

connection

The path integral connects to the rest of physics along several axes. Via Wick Rotation $t\to -i\tau$ , the oscillatory factor $e^{iS/\hbar}$ becomes the convergent Boltzmann-like weight $e^{-S_E/\hbar}$ , linking quantum mechanics to Statistical Mechanics: the Euclidean partition function $Z=\int\mathcal{D}q\,e^{-S_E/\hbar}$ over periodic paths is the thermal partition function with $\beta\hbar=$ period. Generalizing the coordinate $q(t)$ to a field $\phi(x,t)$ yields the Quantum Field Theory path integral, the workhorse of particle physics. The stationary-phase ( $\hbar\to0$ ) limit reproduces Classical Mechanics through the Principle of Least Action and connects to WKB Approximation and semiclassical methods. The operator formulation is recovered through the Canonical Quantization correspondence, and for gauge theories the path integral requires Faddeev–Popov Gauge Fixing.

🔬 Where this appears

⏳ History

- 1933 — P. A. M. Dirac publishes "The Lagrangian in Quantum Mechanics," noting the short-time propagator corresponds to

\exp(iL\,dt/\hbar)

, planting the seed.

- 1942 — Richard Feynman, in his Princeton PhD thesis under John Wheeler, develops a space-time formulation of nonrelativistic quantum mechanics as a sum over paths.

- 1948 — Feynman publishes "Space-Time Approach to Non-Relativistic Quantum Mechanics" (Rev. Mod. Phys. 20, 367), the foundational paper of the path integral.

- 1948–49 — Mark Kac formalizes the Euclidean case via the Feynman–Kac formula, connecting the imaginary-time path integral to Brownian motion and Wiener measure.

- 1949 — Feynman introduces Feynman diagrams ("Space-Time Approach to Quantum Electrodynamics"), giving graphical rules for the perturbative path integral.

- 1950s — Extension to fields; Feynman–Vernon influence functional (1963) for open systems.

- 1965 — Feynman & Hibbs publish Quantum Mechanics and Path Integrals, the canonical textbook; Feynman shares the Nobel Prize.

- 1967 — Faddeev & Popov solve gauge-fixing for non-abelian path integrals, enabling the quantization of Yang–Mills theory.

- 1970s — Wilson, Polyakov, 't Hooft develop Euclidean lattice and instanton methods; the path integral becomes the central language of QFT and statistical field theory.

open-problem

Several foundational and technical issues remain active:

Mathematical rigor in Minkowski signature. There is still no fully satisfactory, universally applicable definition of the real-time oscillatory path integral as a genuine integration theory; constructions rely on analytic continuation from Euclidean space, which can fail for theories with sign problems or no Osterwalder–Schrader reflection positivity.
The fermionic sign / sign problem. Monte Carlo evaluation of Euclidean path integrals at finite density or with real-time dynamics is obstructed by oscillatory/negative weights — an NP-hard problem in general, blocking ab initio simulation of dense QCD and many strongly correlated systems.
Measure for gauge and gravity. Defining $\mathcal{D}g_{\mu\nu}$ for the gravitational path integral is unresolved; the Euclidean action is unbounded below (conformal factor problem), and the very meaning of "summing over geometries" (including topology change) is open.
Constructive QFT. Rigorously constructing interacting 4-dimensional path integrals (e.g. nontrivial $\phi^4_4$ , or Yang–Mills with a mass gap — a Clay Millennium Problem) remains unsolved.
Operator ordering & curved configuration spaces. The path integral on a curved manifold requires extra $\hbar^2$ counterterms whose universal form is still debated.

Formal Statement

Why It Holds / Derivation

Building Intuition

Worked Examples

Common Misconceptions

Connections

Real-World Applications

Historical Development

Subtleties & Limits

See Also