Path Integral Formulation QFT

⚡ Key Takeaways

The path integral formulation expresses quantum amplitudes as a sum over all possible histories of a system, weighted by eiS/ℏe^{iS/\hbar} where SS is the classical action. In QFT, the fundamental objects become functional integrals over field configurations, ∫Dϕ eiS[ϕ]/ℏ\int \mathcal{D}\phi \, e^{iS[\phi]/\hbar}, generating all correlation functions and the S-matrix. It applies to all of relativistic quantum field theory, gauge theory, statistical mechanics, and quantum gravity. Conceived by Dirac (1933), developed fully by Feynman (1948), and made rigorous in Euclidean form by later workers.

💡 The Core Insight

Quantum amplitudes are obtained by superposing contributions from every conceivable history, not just the classical one. Each history ϕ(x)\phi(x) contributes a complex number of unit modulus, eiS[ϕ]/ℏe^{iS[\phi]/\hbar}. Constructive interference near the stationary point of SS reproduces classical physics; everywhere else, rapid phase oscillation produces cancellation. The transition from a single trajectory to a democratic sum over all of them is the conceptual heart of the formalism.

Formal Statement

DefinitionQuantum-Mechanical Path Integral (Single Particle)

For a particle with Lagrangian L(q,q˙)L(q,\dot q), the transition amplitude between position eigenstates is given by the configuration-space path integral:

1
⟨qf,tf∣qi,ti⟩=∫q(ti)=qiq(tf)=qfDq(t)  exp⁡ ⁣(iℏ∫titfL(q,q˙) dt)\quad \langle q_f, t_f \mid q_i, t_i\rangle = \int_{q(t_i)=q_i}^{q(t_f)=q_f} \mathcal{D}q(t)\; \exp\!\left(\frac{i}{\hbar}\int_{t_i}^{t_f} L(q,\dot q)\,dt\right)

where Dq(t)\mathcal{D}q(t) denotes integration over all paths with the stated endpoints, S[q]=∫L dtS[q]=\int L\,dt is the classical action, and ℏ\hbar is the reduced Planck constant.

Derivation hint: Insert N−1N-1 complete sets of position eigenstates between time-sliced evolution operators e−iH ϵ/ℏe^{-iH\,\epsilon/\hbar}, with ϵ=(tf−ti)/N\epsilon=(t_f-t_i)/N, evaluate each short-time kernel, and take N→∞N\to\infty.

DefinitionThe Time-Slicing Measure

The functional measure is the formal limit of a finite-dimensional integral over intermediate positions:

2
Dq(t)=lim⁡N→∞(m2πiℏϵ)N/2∏k=1N−1dqk\quad \mathcal{D}q(t) = \lim_{N\to\infty}\left(\frac{m}{2\pi i \hbar \epsilon}\right)^{N/2}\prod_{k=1}^{N-1} dq_k

where mm is the particle mass and the prefactor is the normalization of each Gaussian short-time kernel.

Derivation hint: The free-particle short-time kernel ⟨qk+1∣e−iHϵ/ℏ∣qk⟩\langle q_{k+1}|e^{-iH\epsilon/\hbar}|q_k\rangle is a Gaussian in (qk+1−qk)(q_{k+1}-q_k); its normalization (m/2πiℏϵ)1/2(m/2\pi i\hbar\epsilon)^{1/2} accumulates over NN slices.

DefinitionField-Theoretic Generating Functional

For a scalar field ϕ(x)\phi(x), all quantum information is packaged in the generating functional with external source J(x)J(x):

3
Z[J]=∫Dϕ  exp⁡ ⁣(iℏ∫d4x [L(ϕ,∂μϕ)+J(x)ϕ(x)])\quad Z[J] = \int \mathcal{D}\phi \; \exp\!\left(\frac{i}{\hbar}\int d^4x\,\big[\mathcal{L}(\phi,\partial_\mu\phi) + J(x)\phi(x)\big]\right)

where L\mathcal{L} is the Lagrangian density and J(x)J(x) is a classical source coupling linearly to the field.

Derivation hint: Promote the single coordinate q(t)q(t) to a field ϕ(x⃗,t)\phi(\vec x,t) indexed continuously by space, replacing ∫L dt\int L\,dt by ∫L d4x\int \mathcal{L}\,d^4x and the discrete product over qkq_k by a product over field values at every spacetime point.

TheoremCorrelation Functions from Functional Derivatives

Time-ordered vacuum correlation functions (Green's functions) are generated by repeated functional differentiation:

4
⟨0∣T ϕ(x1)⋯ϕ(xn)∣0⟩=1Z[0](−iℏ)n δnZ[J]δJ(x1)⋯δJ(xn)∣J=0\quad \langle 0|T\,\phi(x_1)\cdots\phi(x_n)|0\rangle = \frac{1}{Z[0]}\left.\frac{(-i\hbar)^n\,\delta^n Z[J]}{\delta J(x_1)\cdots\delta J(x_n)}\right|_{J=0}

where TT denotes time ordering and the division by Z[0]Z[0] removes vacuum bubble contributions.

Derivation hint: Each functional derivative δ/δJ(xk)\delta/\delta J(x_k) brings down a factor (i/ℏ)ϕ(xk)(i/\hbar)\phi(x_k) from the exponential; setting J=0J=0 leaves the nn-point integrand, and the path integral automatically enforces time ordering.

TheoremEuclidean Path Integral and Wick Rotation

Analytically continuing to imaginary time t→−iτt\to -i\tau converts the oscillatory weight into a real, positive, exponentially damped one:

5
ZE[J]=∫Dϕ  exp⁡ ⁣(−1ℏ SE[ϕ]+1ℏ∫d4xE Jϕ),SE=∫d4xE [12(∂ϕ)2+V(ϕ)]\quad Z_E[J] = \int \mathcal{D}\phi\; \exp\!\left(-\frac{1}{\hbar}\,S_E[\phi] + \frac{1}{\hbar}\int d^4x_E\, J\phi\right),\qquad S_E = \int d^4x_E\,\Big[\tfrac{1}{2}(\partial\phi)^2 + V(\phi)\Big]

where SE≥0S_E\ge 0 is the Euclidean action and d4xEd^4x_E is the Euclidean volume element.

Derivation hint: Substitute t=−iτt=-i\tau into the Minkowski action; the factor i dt=dτi\,dt = d\tau and the sign flip in (∂tϕ)2(\partial_t\phi)^2 turn iSiS into −SE-S_E, rendering the measure a genuine probability weight.

ProofDerivation of the Path Integral from Operator Mechanics

We derive Eq. [1] from canonical quantization in numbered steps. 1. Slice the evolution operator. Write ⟨qf∣e−iH(tf−ti)/ℏ∣qi⟩\langle q_f|e^{-iH(t_f-t_i)/\hbar}|q_i\rangle and split the interval into NN equal pieces of width ϵ\epsilon, so e−iHT/ℏ=(e−iHϵ/ℏ)Ne^{-iHT/\hbar}=\big(e^{-iH\epsilon/\hbar}\big)^N. 2. Insert completeness. Between adjacent factors insert ∫dqk ∣qk⟩⟨qk∣=1\int dq_k\,|q_k\rangle\langle q_k|=\mathbb{1}, yielding a product of short-time kernels ∏k⟨qk+1∣e−iHϵ/ℏ∣qk⟩\prod_k \langle q_{k+1}|e^{-iH\epsilon/\hbar}|q_k\rangle. 3. Evaluate the kernel. For H=p2/2m+V(q)H=p^2/2m+V(q), insert momentum completeness and perform the Gaussian pp-integral:

6
⟨qk+1∣e−iHϵ/ℏ∣qk⟩≈(m2πiℏϵ)1/2exp⁡ ⁣[iϵℏ(m2(qk+1−qkϵ)2−V(qk))].\quad \langle q_{k+1}|e^{-iH\epsilon/\hbar}|q_k\rangle \approx \left(\frac{m}{2\pi i\hbar\epsilon}\right)^{1/2}\exp\!\left[\frac{i\epsilon}{\hbar}\left(\frac{m}{2}\Big(\frac{q_{k+1}-q_k}{\epsilon}\Big)^2 - V(q_k)\right)\right].

Derivation hint: Complete the square in pp after writing e−iVϵ/ℏe^{-iV\epsilon/\hbar} to the right; the pp-integral is a shifted Gaussian. 4. Recognize the action. The exponent summed over slices is iϵℏ∑k[m2q˙k2−V(qk)]→iℏ∫L dt\frac{i\epsilon}{\hbar}\sum_k[\frac{m}{2}\dot q_k^2 - V(q_k)]\to \frac{i}{\hbar}\int L\,dt as ϵ→0\epsilon\to 0. 5. Take the continuum limit. The accumulated normalizations form the measure of Eq. [2], and the product of integrals becomes the functional integral of Eq. [1]. ■\blacksquare

∎
RP
Richard P. Feynman
1918–1988

Feynman, building on a 1933 paper by Dirac that noted the amplitude "corresponds to" eiL dt/ℏe^{iL\,dt/\hbar}, established the path integral as a full reformulation of quantum mechanics in his 1942 Princeton PhD thesis (supervised by Wheeler) and his landmark 1948 Reviews of Modern Physics paper. He used it to derive the Feynman rules and diagrams for QED, work for which he shared the 1965 Nobel Prize with Schwinger and Tomonaga. His view of physics as a "sum over histories" reshaped quantum field theory, statistical mechanics, and our intuition for quantum interference.

Building Intuition

🧠 Intuition

Picture the famous double-slit experiment, then push it to its logical extreme. Add more screens with more slits; the particle's amplitude is the sum over each path through the holes. Now remove the screens entirely and drill "all possible holes everywhere"—the particle effectively traverses every path through space simultaneously, each contributing a phase eiS/ℏe^{iS/\hbar}. Far from the classical trajectory, neighboring paths have wildly different actions, so their phases point in random directions and cancel. Near the classical path, where δS=0\delta S=0, neighboring paths share nearly the same phase and add coherently. Thus the classical trajectory is not "chosen" but emerges as the locus of constructive interference. In QFT, the "particle" becomes a field configuration over all spacetime, and the sum over histories becomes a sum over field shapes—a superposition of every way the field can wiggle between specified boundary data.

Figure
Stationary-phase: classical path dominates by constructive interference

```svg

Sum over paths q(t) q_i, t_i q_f, t_f classical: δS=0 others cancel by oscillating phase e^{iS/ℏ}

```

Worked Examples

ExampleFree-Particle Propagator

Given: A free particle of mass mm, with V=0V=0, propagating from qiq_i at tit_i to qfq_f at tft_f over duration T=tf−tiT=t_f-t_i.

Steps: 1. The classical action for the straight-line trajectory is Scl=m(qf−qi)22TS_{\rm cl}=\frac{m(q_f-q_i)^2}{2T}. 2. Because LL is quadratic, the path integral is exactly Gaussian; fluctuations contribute only a normalization A(T)A(T), giving K=A(T) eiScl/ℏK = A(T)\,e^{iS_{\rm cl}/\hbar}. 3. Determine A(T)A(T) from the Gaussian fluctuation determinant (or from the semigroup property ∫dq′ K(qf,T2;q′)K(q′,T1;qi)=K(qf,T1+T2;qi)\int dq'\,K(q_f,T_2;q')K(q',T_1;q_i)=K(q_f,T_1+T_2;q_i)), yielding A(T)=(m/2πiℏT)1/2A(T)=(m/2\pi i\hbar T)^{1/2}.

Answer:

 K(qf,tf;qi,ti)=(m2πiℏT)1/2exp⁡ ⁣[im(qf−qi)22ℏT] \boxed{\,K(q_f,t_f;q_i,t_i)=\left(\frac{m}{2\pi i\hbar T}\right)^{1/2}\exp\!\left[\frac{i m (q_f-q_i)^2}{2\hbar T}\right]\,}

Numerically, for an electron (m=9.11×10−31m=9.11\times10^{-31} kg), Δq=1\Delta q=1 nm, T=10−15T=10^{-15} s: the phase Scl/ℏ=mΔq2/(2ℏT)≈(9.11×10−31⋅10−18)/(2⋅1.055×10−34⋅10−15)≈4.3S_{\rm cl}/\hbar = m\Delta q^2/(2\hbar T)\approx (9.11\times10^{-31}\cdot10^{-18})/(2\cdot1.055\times10^{-34}\cdot10^{-15})\approx 4.3 radians (dimensionless), confirming strong quantum interference at this scale.

ExampleHarmonic Oscillator via Fluctuation Determinant

Given: A harmonic oscillator L=m2q˙2−mω22q2L=\frac{m}{2}\dot q^2-\frac{m\omega^2}{2}q^2, duration TT.

Steps: 1. Split q(t)=qcl(t)+η(t)q(t)=q_{\rm cl}(t)+\eta(t) with η\eta vanishing at endpoints; the action separates as S=Scl+m2∫(η˙2−ω2η2)dtS=S_{\rm cl}+\frac{m}{2}\int(\dot\eta^2-\omega^2\eta^2)dt. 2. The classical action is Scl=mω2sin⁡ωT[(qi2+qf2)cos⁡ωT−2qiqf]S_{\rm cl}=\frac{m\omega}{2\sin\omega T}\big[(q_i^2+q_f^2)\cos\omega T-2q_iq_f\big]. 3. The fluctuation integral gives det⁡−1/2\det^{-1/2} of the operator −∂t2−ω2-\partial_t^2-\omega^2; its ratio to the free case is ωT/sin⁡ωT\sqrt{\omega T/\sin\omega T}, fixing the prefactor.

Answer:

 K=(mω2πiℏsin⁡ωT)1/2exp⁡ ⁣[imω2ℏsin⁡ωT((qi2+qf2)cos⁡ωT−2qiqf)] \boxed{\,K=\left(\frac{m\omega}{2\pi i\hbar\sin\omega T}\right)^{1/2}\exp\!\left[\frac{i m\omega}{2\hbar\sin\omega T}\Big((q_i^2+q_f^2)\cos\omega T-2q_iq_f\Big)\right]\,}

For ωT→0\omega T\to 0 this reduces smoothly to the free-particle result; poles at ωT=nπ\omega T=n\pi signal the focusing of all classical paths (caustics).

Common Misconceptions

⚠️ Common mistake

The measure Dϕ\mathcal{D}\phi is not a genuine measure. In Minkowski signature the "weight" eiS/ℏe^{iS/\hbar} has unit modulus and the formal infinite-dimensional integral does not converge absolutely; it is defined only through limiting/regularization procedures (time-slicing, Euclidean continuation, lattice). Treating Dϕ eiS/ℏ\mathcal{D}\phi\,e^{iS/\hbar} as a literal Lebesgue-type measure is mathematically false—no translation-invariant countably-additive measure exists on infinite-dimensional spaces. Only the Euclidean version, Eq. [5], approaches a true (Wiener-like, Gaussian-based) measure, and even then interactions require renormalization.

Counterexample

Naïve path integration violates operator ordering. One might think the configuration-space measure of Eq. [2] always captures the Hamiltonian faithfully. It does not when the kinetic term is position-dependent—e.g. a particle on a curved manifold with H=12gij(q)pipjH=\frac{1}{2}g^{ij}(q)p_ip_j. The naïve discretization gij(qk)q˙iq˙jg_{ij}(q_k)\dot q^i\dot q^j misses an O(ℏ2)O(\hbar^2) curvature term (proportional to the Ricci scalar RR) that depends on the discretization point (midpoint vs. endpoint). The correct prescription comes from the phase-space path integral with careful ordering, demonstrating that the configuration-space form is not universally valid.

Connections

connection

The path integral sits at the crossroads of much of theoretical physics. It reduces to ordinary Canonical Quantization via time-slicing, and the two are equivalent for well-behaved Hamiltonians. The Euclidean version, Eq. [5], is mathematically identical to the partition function of Statistical Mechanics, with ℏ\hbar playing the role of temperature and the action playing the role of energy—this is the bedrock of Quantum Field Theory at Finite Temperature. It provides the natural language for Gauge Theory through the Faddeev–Popov procedure and ghost fields, and underlies the Renormalization Group when one integrates out high-momentum field modes shell by shell. The stationary-phase limit connects it directly to Classical Mechanics and the Hamilton-Jacobi Equation, while instanton solutions of Eq. [5] describe Quantum Tunneling non-perturbatively.

Real-World Applications

🔬 Where this appears

Historical Development

⏳ History
- 1923–1925: Wiener constructs the rigorous measure for Brownian motion—the Euclidean ancestor of the path integral.
- 1933: Dirac, in Physikalische Zeitschrift der Sowjetunion, observes that the quantum transition amplitude over a short time "corresponds to" exp⁡(iL dt/ℏ)\exp(iL\,dt/\hbar), planting the seed.
- 1942: Feynman's Princeton PhD thesis, "The Principle of Least Action in Quantum Mechanics," develops the sum-over-histories idea fully.
- 1948: Feynman publishes "Space-Time Approach to Non-Relativistic Quantum Mechanics" (Rev. Mod. Phys.), establishing the formalism; companion work derives QED Feynman rules.
- 1949: Dyson proves the equivalence of Feynman's, Schwinger's, and Tomonaga's formulations of QED.
- 1950s–60s: Matthews–Salam, and later Faddeev–Popov (1967), extend the path integral to non-abelian gauge theories, introducing gauge-fixing and ghosts.
- 1972: 't Hooft and Veltman use path-integral (dimensional) regularization to prove the renormalizability of gauge theories.
- 1974: Wilson formulates lattice gauge theory, making the Euclidean path integral a computational tool and explaining confinement.
- 1975–83: Coleman (instantons/false vacuum), Gibbons–Hawking and Hartle–Hawking apply Euclidean path integrals to tunneling, black-hole thermodynamics, and quantum cosmology.

Subtleties & Limits

open-problem

Despite its practical triumphs, the path integral resists full mathematical rigor. The Minkowski measure Dϕ eiS/ℏ\mathcal{D}\phi\,e^{iS/\hbar} has no known countably-additive realization; constructive quantum field theory (Glimm–Jaffe) succeeds only for super-renormalizable models in d≤3d\le 3 and remains open for four-dimensional Yang–Mills (a Clay Millennium Problem). The gravitational path integral faces the conformal-factor problem (the Euclidean Einstein–Hilbert action is unbounded below) and ambiguities in the contour of integration. Convergence and the sign problem plague Euclidean lattice simulations at finite chemical potential, where the weight becomes complex. Finally, the precise definition of the functional measure including its anomalous Jacobians (Fujikawa's derivation of chiral anomalies) reveals that the "measure" itself carries physical content not visible in the classical action—an aspect still being mapped in curved space and on the lattice.

See Also

Canonical Quantization · Feynman Diagrams · Gauge Theory · Renormalization Group · Statistical Mechanics · Wick Rotation · Effective Action · Instantons and Tunneling · Lattice Gauge Theory