QFT: Path Integrals & Methods
Path Integrals & Methods
Quantum Field Theory » Path Integrals & Methods
Knowing what a quantum field is tells you almost nothing about how to compute with one. This page is the workshop: the path integral that packages all of QFT into a single weighted sum over field configurations, the generating functionals that turn that sum into a machine for spitting out correlation functions, the perturbative expansion of that machine into Feynman diagrams, and the effective-field-theory mindset that lets you calculate without knowing the full theory. These are the tools that produce the numbers — the electron g−2 to twelve digits, cross sections at the LHC, the running of the strong coupling.
- Sum over histories. Every field configuration contributes an amplitude $e^{iS/\hbar}$; the classical path is just where the phase is stationary.
- Diagrams are the expansion. Each Feynman diagram is one term in a power series in the coupling — propagators for lines, factors for vertices.
- One functional, all correlators. Differentiate the generating functional $Z[J]$ and every Green’s function falls out automatically.
- Calculate at the right scale. Effective field theory lets you ignore physics you cannot reach and organize the rest by powers of energy.
What You’ll Find on This Page
| Section | What it covers |
|---|---|
| Path-Integral Formulation | Sum over histories, the action, Euclidean rotation |
| Generating Functionals | $Z[J]$, $W[J]$, the effective action, Gaussian integrals |
| Perturbation Theory | Expanding the interacting theory, Wick’s theorem |
| Feynman Diagrams | Reading and building diagrams; the Feynman rules |
| A Worked Amplitude | Tree-level scattering from start to finish |
| Loops & Functional Methods | Loop integrals, Ward identities, Schwinger–Dyson, BRST |
| Effective Field Theory | EFT as a practical organizing principle |
How the Pieces Fit Together
graph TD
PI["Path integral: sum over field configs"] --> Z["Generating functional Z[J]"]
Z --> W["Connected functional W[J] = -i ln Z"]
W --> G["Effective action Γ[φc] (Legendre transform)"]
Z --> PT["Perturbative expansion"]
PT --> WICK["Wick's theorem"]
WICK --> FD["Feynman diagrams"]
FD --> AMP["Scattering amplitudes / cross sections"]
EFT["Effective field theory"] --> PT
style PI fill:#11998e,color:#fff
style AMP fill:#38ef7d,color:#222
style EFT fill:#ccf,color:#222
The Path-Integral Formulation
The path integral, due to Feynman (building on a remark of Dirac), provides an alternative to canonical quantization that is fully equivalent but often far more powerful. Instead of promoting fields to operators and imposing commutation relations, it assigns a complex amplitude to every possible field configuration and sums them all up.
Sum Over Histories
In ordinary quantum mechanics, the amplitude for a particle to travel from $x_i$ at time $t_i$ to $x_f$ at time $t_f$ is obtained by summing $e^{iS/\hbar}$ over all paths connecting the endpoints — not just the classical one:
\[\langle x_f, t_f | x_i, t_i \rangle = \int_{x(t_i)=x_i}^{x(t_f)=x_f} \mathcal{D}x \; e^{iS[x]/\hbar}\]The classical trajectory is the one where the action is stationary ($\delta S = 0$), so nearby paths interfere constructively; far from it, the rapidly oscillating phase causes cancellation. This is how the classical limit ($\hbar \to 0$) emerges: only the stationary-action path survives.
The Field-Theory Path Integral
Promoting the single coordinate $x(t)$ to a field $\phi(x)$ defined at every spacetime point, the transition amplitude between field configurations becomes a functional integral:
\[\langle \phi_f, t_f | \phi_i, t_i \rangle = \int_{\phi(t_i)=\phi_i}^{\phi(t_f)=\phi_f} \mathcal{D}\phi \; e^{iS[\phi]/\hbar}\]where the action is the spacetime integral of the Lagrangian density:
\[S[\phi] = \int_{t_i}^{t_f} dt \int d^3x \; \mathcal{L}\big[\phi(x,t),\, \partial_\mu\phi(x,t)\big]\]The symbol $\mathcal{D}\phi$ denotes integration over every possible field configuration — a continuous infinity of ordinary integrals, one for the value of $\phi$ at each spacetime point. This object is the master quantity from which everything else is derived.
Euclidean Formulation
The oscillatory factor $e^{iS}$ makes the Minkowski path integral only conditionally convergent. The standard cure is Wick rotation: analytically continue to imaginary time, $t \to -i\tau$. The action picks up factors of $i$ that convert it into the Euclidean action $S_E$, and the integrand becomes a real, exponentially damped weight:
\[Z_E = \int \mathcal{D}\phi \; e^{-S_E[\phi]/\hbar}\]This has two enormous payoffs:
- Convergence. The integrand is a genuine probability-like weight; field configurations with large action are exponentially suppressed, making the integral well-defined and amenable to numerical evaluation (this is the basis of lattice QFT).
- Statistical-mechanics connection. $Z_E$ is mathematically identical to the partition function of a classical statistical system, with $S_E/\hbar$ playing the role of $\beta H$. Correlation lengths map to inverse masses, phase transitions to critical phenomena, and the renormalization group is shared between the two fields. (See Statistical Mechanics.)
Generating Functionals
The path integral becomes a calculational tool the moment we couple the field to an external source $J(x)$. Differentiating with respect to that source pulls down factors of the field, so a single functional encodes every correlation function at once.
The Generating Functional Z[J]
Add a linear source term $J(x)\phi(x)$ to the action:
\[Z[J] = \int \mathcal{D}\phi \; e^{\,i\left(S[\phi] + \int d^4x \, J(x)\phi(x)\right)}\]Functional derivatives with respect to $J$ bring down fields, and setting $J=0$ leaves the time-ordered vacuum correlation functions (Green’s functions):
\[\langle 0|T[\phi(x_1)\cdots\phi(x_n)]|0\rangle = \frac{1}{Z[0]} \frac{(-i)^n \, \delta^n Z[J]}{\delta J(x_1)\cdots\delta J(x_n)}\bigg|_{J=0}\]This is the central identity of the functional approach: all the physical content of the theory is packed into $Z[J]$, and any correlator is extracted by mechanical differentiation.
Connected Functional W[J]
The full Green’s functions contain redundant, disconnected pieces (processes happening independently in different regions). Taking the logarithm strips these away, leaving only connected correlators:
\[W[J] = -i \ln Z[J]\] \[\langle 0|T[\phi(x_1)\cdots\phi(x_n)]|0\rangle_c = (-i)^{n-1} \frac{\delta^n W[J]}{\delta J(x_1)\cdots\delta J(x_n)}\bigg|_{J=0}\]Physically, $W[J]$ generates exactly the diagrams that cannot be split into independent pieces — the ones that carry real interaction information.
The Effective Action
One Legendre transform further isolates the one-particle-irreducible (1PI) content — diagrams that stay connected when any single internal line is cut. Define the classical field $\phi_c = \delta W/\delta J$ and form:
\[\Gamma[\phi_c] = W[J] - \int d^4x \; J(x)\,\phi_c(x)\]$\Gamma[\phi_c]$ is the effective action. Its great virtues:
- Its derivatives are the 1PI vertex functions — the irreducible building blocks from which all diagrams are assembled.
- It includes all quantum corrections to the classical action, so extremizing $\Gamma$ (rather than $S$) gives the quantum equations of motion.
- Its minimum locates the true vacuum, making it the natural tool for analyzing spontaneous symmetry breaking: the effective potential $V_{\text{eff}}(\phi_c)$ is the part of $\Gamma$ with no derivatives, and one-loop corrections to it (the Coleman–Weinberg potential) can shift the vacuum away from the classical minimum.
Gaussian Integration: the Free Theory
For a free field the action is quadratic, and the path integral reduces to an infinite-dimensional Gaussian — which we can do exactly. Writing the quadratic form as $\phi K \phi$:
\[Z_0 = \int \mathcal{D}\phi \, \exp\left[\frac{i}{2} \int d^4x\, d^4y \; \phi(x)\,K(x,y)\,\phi(y)\right] = (\det K)^{-1/2}\]The two-point function is the inverse of the kinetic operator — which is precisely the Feynman propagator:
\[\langle 0|T[\phi(x)\phi(y)]|0\rangle_0 = K^{-1}(x,y) = D_F(x-y)\]For the Klein–Gordon field, $K = -(\Box + m^2)$ and its inverse in momentum space is the familiar
\[\tilde{D}_F(k) = \frac{i}{k^2 - m^2 + i\varepsilon}\]This is the cornerstone of perturbation theory: the free theory is solved in closed form, and interactions are added as a controlled expansion around it.
Perturbation Theory
Realistic theories are interacting, and almost none are exactly solvable. The strategy is to split the Lagrangian into a free (quadratic) part we can integrate exactly and an interaction part we expand in powers of the small coupling.
Expanding Around the Free Theory
Write $\mathcal{L} = \mathcal{L}0 + \mathcal{L}{\text{int}}$. The interacting generating functional can be written by pulling the interaction outside the integral as a differential operator acting on the free functional:
\[Z[J] = \exp\left[\,i\int d^4x \; \mathcal{L}_{\text{int}}\!\left(\frac{1}{i}\frac{\delta}{\delta J(x)}\right)\right] Z_0[J]\]Expanding the exponential generates an infinite series, each term containing some number of interaction vertices acting on the free propagators inside $Z_0[J]$. This series is the perturbative expansion, and each term corresponds to a Feynman diagram.
Wick’s Theorem
To evaluate a given term we must reduce a time-ordered product of many fields into propagators. Wick’s theorem does exactly this: it rewrites a time-ordered product as a sum over all possible pairwise contractions, where each contraction is a Feynman propagator:
\[T[\phi(x_1)\cdots\phi(x_n)] = \;:\!\phi(x_1)\cdots\phi(x_n)\!: \;+\; \text{(all contractions)}\]Each full contraction of $n$ fields (with $n$ even) into $n/2$ propagators contributes a product like $D_F(x_1-x_2)\,D_F(x_3-x_4)\cdots$. The normal-ordered terms vanish between vacuum states, so only the fully contracted pieces survive in vacuum correlators. Wick’s theorem is the bookkeeping device that turns the abstract functional expansion into a concrete, finite set of propagator products at each order — precisely the lines of a Feynman diagram.
Organizing the Expansion
Two complementary expansions are in play:
- Coupling expansion. Powers of the coupling $g$ (or $\lambda$, $e$) count vertices; weak coupling means few vertices dominate.
- Loop expansion. Powers of $\hbar$ count independent loop momenta. The number of loops is $L = I - V + 1$ for $I$ internal lines and $V$ vertices, so tree diagrams ($L=0$) are the classical/leading result and loops are successive quantum corrections.
For QED the relevant small parameter is the fine-structure constant $\alpha = e^2/4\pi \approx 1/137$, which is why the perturbative series converges so well numerically — each additional loop costs roughly a factor of $\alpha$.
Feynman Diagrams
Feynman diagrams are not mere illustrations; they are a precise, one-to-one shorthand for the terms of the perturbation series. Each diagram translates, via the Feynman rules, into a specific mathematical expression contributing to an amplitude.
The Dictionary
graph LR
A["External line"] -->|"incoming/outgoing particle"| B["wavefunction factor"]
C["Internal line"] -->|"virtual particle"| D["propagator i/(p²−m²+iε)"]
E["Vertex"] -->|"interaction"| F["coupling factor, e.g. −ieγ^μ"]
G["Closed loop"] -->|"unfixed momentum"| H["integrate ∫d⁴k/(2π)⁴"]
style A fill:#11998e,color:#fff
style E fill:#38ef7d,color:#222
Reading a diagram:
- External lines represent the incoming and outgoing real particles; each carries a wavefunction factor (a spinor $u(p)$, $\bar u(p)$, $v(p)$, $\bar v(p)$ for fermions, or a polarization $\varepsilon^\mu$ for photons).
- Internal lines are virtual particles — off-shell intermediate states that need not satisfy $p^2 = m^2$. Each contributes a propagator.
- Vertices encode the interaction and carry the coupling. Momentum is conserved at every vertex.
- Loops carry an undetermined internal momentum that must be integrated over; these are the source of both quantum corrections and ultraviolet divergences.
Feynman Rules for QED
The interaction in quantum electrodynamics is the single term $\mathcal{L}{\text{int}} = -e\,\bar\psi\gamma^\mu\psi\,A\mu$, giving a remarkably compact set of rules.
Vertex factor (one photon, one incoming and one outgoing electron line):
\[-i e \gamma^\mu\]Electron (fermion) propagator:
\[S_F(p) = \frac{i(\not{p} + m)}{p^2 - m^2 + i\varepsilon} = \frac{i}{\not{p} - m + i\varepsilon}\]Photon propagator (Feynman gauge):
\[D^{\mu\nu}_F(k) = \frac{-i g^{\mu\nu}}{k^2 + i\varepsilon}\]Closed fermion loop: include a factor of $(-1)$ and a trace over the Dirac indices around the loop.
Loop integration: integrate each unfixed loop momentum with $\displaystyle\int \frac{d^4k}{(2\pi)^4}$.
Assembling these factors according to the topology of the diagram, then summing over all diagrams at the desired order, yields the invariant amplitude $i\mathcal{M}$.
From Amplitude to Observable
The amplitude $\mathcal{M}$ feeds directly into measurable quantities. For a $2 \to 2$ scattering in the centre-of-mass frame, the differential cross section is
\[\frac{d\sigma}{d\Omega} = \frac{1}{64\pi^2 s}\,\frac{|\mathbf{p}_f|}{|\mathbf{p}_i|}\,\overline{|\mathcal{M}|^2}\]| where $\sqrt{s}$ is the total centre-of-mass energy and $\overline{ | \mathcal{M} | ^2}$ denotes the spin-averaged (initial) and spin-summed (final) squared amplitude. For decays, $\mathcal{M}$ enters the decay rate $\Gamma$ through an analogous phase-space formula. This is the bridge from diagrams to numbers an experimentalist can check. |
A Worked Amplitude
To make the machinery concrete, consider electron–muon scattering, $e^- \mu^- \to e^- \mu^-$, at lowest (tree) order in QED. Because the electron and muon are distinct particles, there is a single tree diagram: the two fermion lines exchange one virtual photon.
Step 1 — Identify the diagram. Incoming electron $p_1$ and muon $p_2$; outgoing electron $p_3$ and muon $p_4$. A virtual photon of momentum $q = p_1 - p_3$ connects the two vertices. Momentum conservation: $q = p_1 - p_3 = p_4 - p_2$.
Step 2 — Apply the rules. Write down the factors:
- Electron line: $\bar u(p_3)\,(-ie\gamma^\mu)\,u(p_1)$
- Muon line: $\bar u(p_4)\,(-ie\gamma^\nu)\,u(p_2)$
- Photon propagator: $\dfrac{-ig_{\mu\nu}}{q^2}$
Step 3 — Assemble the amplitude.
\[i\mathcal{M} = \big[\bar u(p_3)(-ie\gamma^\mu)u(p_1)\big]\,\frac{-i g_{\mu\nu}}{q^2}\,\big[\bar u(p_4)(-ie\gamma^\nu)u(p_2)\big]\]which simplifies to
\[\mathcal{M} = \frac{e^2}{q^2}\,\big[\bar u(p_3)\gamma^\mu u(p_1)\big]\big[\bar u(p_4)\gamma_\mu u(p_2)\big]\]Step 4 — Square and average over spins. Using the spin sums $\sum_s u(p)\bar u(p) = \not{p} + m$, the spin-averaged squared amplitude becomes a product of Dirac traces:
\[\overline{|\mathcal{M}|^2} = \frac{e^4}{4q^4}\,\text{Tr}\!\big[(\not{p}_3 + m_e)\gamma^\mu(\not{p}_1 + m_e)\gamma^\nu\big]\,\text{Tr}\!\big[(\not{p}_4 + m_\mu)\gamma_\mu(\not{p}_2 + m_\mu)\gamma_\nu\big]\]Evaluating the traces (using $\text{Tr}[\gamma^\mu\gamma^\nu] = 4g^{\mu\nu}$ and the four-gamma trace identity) and expressing the result in Mandelstam variables $s,t,u$ gives, in the high-energy limit $m_e, m_\mu \to 0$:
\[\overline{|\mathcal{M}|^2} = \frac{2e^4\,(s^2 + u^2)}{t^2}\]with $t = q^2$. Step 5 — Insert into the cross-section formula above to obtain $d\sigma/d\Omega$. This is the entire pipeline — diagram to rule to trace to observable — that underlies every QED prediction; adding the muon’s internal structure or moving to identical particles (Bhabha, Møller scattering) only changes which diagrams appear.
Loops & Functional Methods
Beyond tree level, diagrams contain closed loops with unconstrained internal momenta. These loops carry the genuine quantum corrections — and the divergences that motivate renormalization.
Loop Integrals
A one-loop correction requires integrating over the loop momentum. The QED electron self-energy, for instance, is
\[\Sigma(p) = -i e^2 \int \frac{d^4k}{(2\pi)^4} \; \frac{\gamma^\mu(\not{p}-\not{k}+m)\gamma_\mu}{\big[(p-k)^2 - m^2 + i\varepsilon\big]\big[k^2 + i\varepsilon\big]}\]and the vertex correction is
\[\Lambda^\mu(p',p) = -i e^2 \int \frac{d^4k}{(2\pi)^4} \; \frac{\gamma^\nu(\not{p}'-\not{k}+m)\gamma^\mu(\not{p}-\not{k}+m)\gamma_\nu}{\big[(p'-k)^2 - m^2\big]\big[(p-k)^2 - m^2\big]\big[k^2\big]}\]These integrals typically diverge in the ultraviolet. Handling them — via dimensional regularization, counterterms, and the renormalization group — is the subject of the renormalization page. The vertex correction $\Lambda^\mu$ is also exactly what produces the famous anomalous magnetic moment, the electron $g-2$, agreeing with experiment to twelve digits.
Practical Techniques for Loops
A handful of standard tools turn loop integrals into tractable expressions:
- Feynman parameters. Combine denominators into a single quadratic via $\frac{1}{AB} = \int_0^1 dx\,[xA + (1-x)B]^{-2}$, then shift the loop momentum to complete the square.
- Wick rotation. Rotate $k^0 \to ik^0_E$ to turn the Minkowski integral into a convergent Euclidean one.
- Dimensional regularization. Compute in $d = 4 - \varepsilon$ dimensions, isolating divergences as $1/\varepsilon$ poles while preserving gauge invariance.
- Integration by parts (IBP). Use $\int d^d k\,\partial_\mu[k^\mu f(k)] = 0$ to reduce a large family of integrals to a small basis of master integrals — the workhorse of modern multi-loop calculations.
Functional Identities: Ward & Schwinger–Dyson
The functional formalism does more than reproduce diagrams; it yields exact, all-orders relations.
Schwinger–Dyson equations are the quantum equations of motion. They follow from the fact that a path integral is invariant under a shift of the integration variable, $\phi \to \phi + \delta\phi$:
\[\left\langle \frac{\delta S}{\delta \phi(x)} \right\rangle = \text{(contact terms)}\]They form an infinite, coupled hierarchy relating $n$-point to $(n+1)$-point functions — the non-perturbative skeleton of the theory.
Ward–Takahashi identities are the special case enforced by a symmetry of the action. In QED, gauge invariance ties the vertex function to the electron propagator:
\[q_\mu \Gamma^\mu(p', p) = S_F^{-1}(p') - S_F^{-1}(p)\]These identities guarantee that gauge invariance survives renormalization (e.g. relating the charge and field renormalization constants), and they protect the photon from acquiring a mass. They are a powerful consistency check on any calculation.
Gauge Fixing and BRST
Naively, the photon path integral diverges because gauge-equivalent configurations are summed redundantly. The Faddeev–Popov procedure factors out this redundancy by inserting a gauge-fixing condition, at the cost of introducing ghost fields (anticommuting scalars) in non-abelian theories. The whole gauge-fixed action then possesses a residual rigid fermionic symmetry — BRST symmetry — which encodes the original gauge invariance at the quantum level and guarantees unitarity of physical amplitudes. BRST is the modern, systematic way to quantize Yang–Mills theories within the path integral.
Effective Field Theory as a Calculational Tool
Effective field theory (EFT) is less a specific model than a way of thinking that makes QFT usable: you never need the complete theory of everything to compute a low-energy process — you only need the right degrees of freedom at the scale you care about.
The Core Idea
If a process occurs at energy $E$ well below some heavy scale $\Lambda$, the heavy physics cannot be produced directly. Its effects appear only indirectly, suppressed by powers of $E/\Lambda$. EFT makes this systematic: integrate out the heavy fields and write the most general local Lagrangian for the light fields consistent with the symmetries, organized as an expansion in $1/\Lambda$:
\[\mathcal{L}_{\text{eff}} = \mathcal{L}_{d \le 4} + \sum_{i} \frac{c_i}{\Lambda^{n_i - 4}}\,\mathcal{O}_i\]Each operator $\mathcal{O}_i$ of mass dimension $n_i$ comes with a dimensionless Wilson coefficient $c_i$. Higher-dimension operators are suppressed by more powers of $\Lambda$, so at any desired accuracy only finitely many operators matter. This is power counting, and it is what makes a non-renormalizable theory perfectly predictive at low energies.
Why It Works
- Decoupling. Heavy particles decouple from low-energy physics, leaving only their imprint in the Wilson coefficients. You can compute without ever resolving the heavy sector.
- Predictivity with a cutoff. A theory with infinitely many couplings is still predictive: truncating the $1/\Lambda$ expansion at finite order leaves a finite, controlled error of order $(E/\Lambda)^{n}$.
- Matching and running. Wilson coefficients are fixed by matching the EFT to the full theory at the scale $\Lambda$, then run down to the experimental scale using the renormalization group — automatically resumming large logarithms $\ln(\Lambda/E)$.
EFT in Practice
| Effective theory | Light fields | Heavy scale integrated out |
|---|---|---|
| Fermi theory of $\beta$-decay | leptons, nucleons | $W$ boson mass $m_W$ |
| Chiral perturbation theory | pions (Goldstone bosons) | QCD scale $\Lambda_{\text{QCD}}$ |
| Heavy quark effective theory | light quarks, gluons | heavy quark mass $m_Q$ |
| Standard Model EFT (SMEFT) | Standard Model fields | unknown new-physics scale |
| Euler–Heisenberg | photons | electron mass $m_e$ |
The classic example is Fermi’s four-fermion theory of beta decay: long before the $W$ boson was known, the weak interaction could be described by a contact interaction with coupling $G_F \sim g^2/m_W^2$. The full electroweak theory “matches onto” this EFT in the limit $q^2 \ll m_W^2$, with the heavy $W$ propagator $\sim 1/(q^2 - m_W^2) \to -1/m_W^2$ shrinking to a point. Crucially, the entire Standard Model is itself almost certainly an EFT — the leading terms of a more complete theory whose new physics lives at a scale we have not yet reached. (See Effective Field Theory as a Calculational Tool below.)
Key Takeaways
- The path integral is the master tool. Summing $e^{iS/\hbar}$ over all field configurations reproduces all of QFT and connects directly to statistical mechanics after Wick rotation.
- Generating functionals automate correlators. $Z[J]$ gives all Green’s functions, $W[J]$ the connected ones, and the effective action $\Gamma[\phi_c]$ the 1PI vertices and the quantum vacuum.
- Diagrams are the perturbation series. Wick’s theorem turns the expansion into propagator products; each Feynman diagram is one term, with rules for lines, vertices, and loops.
- Amplitudes become observables. Squaring $\mathcal{M}$, averaging over spins, and inserting phase space yields cross sections and decay rates an experiment can test.
- Functional identities are exact. Schwinger–Dyson and Ward–Takahashi relations hold to all orders, protecting gauge invariance and constraining every calculation.
- EFT lets you calculate at the right scale. Integrate out heavy physics, organize by powers of $E/\Lambda$, and compute predictively without knowing the complete theory.
See Also
- Quantum Field Theory — fields, gauge symmetry, the Standard Model, and renormalization.
- Quantum Mechanics — the non-relativistic foundation the path integral generalizes.
- Statistical Mechanics — the Euclidean path integral is a statistical partition function.
- Relativity — special relativity makes the action Lorentz-invariant.
- String Theory — worldsheet path integrals extend these methods to extended objects.
- Physics Hub — browse all physics topics.