Feynman Diagrams

⚡ Key Takeaways

Feynman diagrams are graphical representations of terms in the perturbative expansion of scattering amplitudes in quantum field theory. Each diagram encodes a specific integral over internal momenta via the Feynman rules, with vertices, propagators, and external legs. They apply wherever a weakly-coupled QFT or many-body theory admits a perturbative series — QED, QCD (perturbative regime), electroweak theory, and condensed matter. Introduced by Richard Feynman (1948–1949), with parallel formalisms by Schwinger and Tomonaga, and unified by Dyson.

💡 The Core Insight

The time-ordered exponential of the interaction Hamiltonian — the Dyson series — when evaluated between asymptotic in/out states using Wick's theorem, decomposes into a sum of contractions. Each distinct contraction pattern is a graph. Feynman's genius was to read the graph backwards: instead of computing the algebra and drawing the picture afterward, one draws all topologically distinct diagrams of a given order and reads off the amplitude directly from a small set of rules. Combinatorics becomes geometry.

DefinitionPerturbative S-matrix and the Dyson Series

The scattering operator (S-matrix) is the time-ordered exponential of the interaction Hamiltonian density $\mathcal{H}_I$ :

1

\quad S = T\exp\left(-i\int d^4x\, \mathcal{H}_I(x)\right) = \sum_{n=0}^{\infty} \frac{(-i)^n}{n!}\int d^4x_1\cdots d^4x_n\, T\{\mathcal{H}_I(x_1)\cdots\mathcal{H}_I(x_n)\}

Derivation hint: Solve the interaction-picture equation $i\,\partial_t U = H_I(t)U$ iteratively (Picard iteration); the time-ordering symbol $T$ absorbs the $1/n!$ from symmetrizing the nested integration domains. Symbols: $T$ = time-ordering operator (later times to the left); $\mathcal{H}_I$ = interaction Hamiltonian density; $x_i$ = spacetime four-vectors; $n$ = perturbative order.

DefinitionFeynman Rules in Momentum Space (φ⁴ theory)

For the scalar theory $\mathcal{L} = \tfrac12(\partial\phi)^2 - \tfrac12 m^2\phi^2 - \tfrac{\lambda}{4!}\phi^4$ , the amputated amplitude $i\mathcal{M}$ is assembled from:

2

\quad \text{Propagator:}\quad \frac{i}{p^2 - m^2 + i\epsilon}, \qquad \text{Vertex:}\quad -i\lambda\,(2\pi)^4\delta^4\!\Big(\textstyle\sum p_{\text{in}} - \sum p_{\text{out}}\Big)

Derivation hint: The propagator is the Fourier transform of the Feynman two-point function $\langle 0|T\phi(x)\phi(y)|0\rangle$ ; the vertex factor $-i\lambda$ comes from one power of the interaction term, with $4! = 24$ canceling the $1/4!$ for distinct leg assignments. Each independent loop carries an integral $\int \frac{d^4k}{(2\pi)^4}$ ; momentum is conserved at each vertex; divide by the symmetry factor of the diagram. Symbols: $p$ = four-momentum, $m$ = mass, $\lambda$ = coupling, $k$ = loop momentum, $\epsilon \to 0^+$ = Feynman pole prescription.

TheoremWick's Theorem

Any time-ordered product of free fields equals the sum over all possible full and partial contractions, normal-ordered:

3

\quad T\{\phi(x_1)\cdots\phi(x_n)\} = \sum_{\text{contractions}} :\!\phi(x_{i_1})\cdots\,\contraction{}{\phi}{(x_a)}{\phi}\phi(x_a)\phi(x_b)\cdots\!:

where each contraction $\contraction{}{\phi}{}{\phi}\phi(x_a)\phi(x_b) = D_F(x_a - x_b)$ is the Feynman propagator. Derivation hint: Decompose each field into positive/negative frequency parts, $\phi = \phi^+ + \phi^-$ , and induct on the number of fields, commuting annihilation operators rightward; each commutator generates a propagator.

TheoremLinked-Cluster / Connected Diagrams

The generating functional satisfies $Z[J] = e^{iW[J]}$ , where $W[J]$ generates only connected diagrams. Consequently disconnected vacuum bubbles exponentiate and cancel in normalized correlation functions:

4

\quad \langle\Omega|T\{\phi(x_1)\cdots\phi(x_n)\}|\Omega\rangle = \frac{\langle 0|T\{\phi\cdots\phi\,e^{-i\int\mathcal{H}_I}\}|0\rangle}{\langle 0|T\{e^{-i\int\mathcal{H}_I}\}|0\rangle}

Derivation hint: Factor every diagram into connected pieces; the sum over all bubble configurations factorizes into $\exp(\sum \text{connected bubbles})$ , which is precisely the denominator.

ProofFrom Dyson Series to Diagrams

Step 1. Write a generic S-matrix element $\langle f|S|i\rangle$ between asymptotic Fock states, inserting the expansion [1].

Step 2. At each order $n$ , the matrix element is a vacuum expectation value of a time-ordered product of fields (from $\mathcal{H}_I$ ) and the creation/annihilation operators defining $|i\rangle,|f\rangle$ .

Step 3. Apply Wick's theorem [3]. Fully-contracted terms survive between the vacuum; each contraction is a propagator $D_F$ , and contractions with external operators give external legs.

Step 4. Catalogue the contraction patterns topologically. Each distinct pattern ↔ one Feynman diagram. Vertices ↔ insertions of $\mathcal{H}_I$ ; internal lines ↔ internal contractions; external lines ↔ contractions with $|i\rangle,|f\rangle$ .

Step 5. Count combinatorics. The $1/n!$ from [1] is cancelled by the $n!$ ways of assigning vertices to a labeled diagram, except for residual symmetry: when a diagram is invariant under a permutation of its lines/vertices, divide by the symmetry factor $S$ .

Step 6. Fourier transform to momentum space. Each vertex integral $\int d^4x_i$ produces a momentum-conserving delta. One overall delta enforces total conservation; the rest fix internal momenta, leaving $L = I - V + 1$ unconstrained loop integrals. This yields rule [2]. $\blacksquare$

∎

RP

Richard P. Feynman

1918–1988

Feynman introduced his diagrams at the 1948 Pocono Conference, initially to skepticism from Bohr and others who thought he was reviving naive spacetime trajectories. His 1949 papers "Space-Time Approach to Quantum Electrodynamics" and "The Theory of Positrons" established the rules and the interpretation of antiparticles as negative-energy solutions propagating backward in time. He shared the 1965 Nobel Prize with Julian Schwinger and Sin-Itiro Tomonaga. Freeman Dyson proved the equivalence of Feynman's intuitive diagrams with the operator formalisms of Schwinger and Tomonaga, giving the method its rigorous footing.

🧠 Intuition

Think of a scattering amplitude as a sum over stories of how an initial state could become a final state. Each story is a diagram. The propagator is the amplitude for a disturbance to travel between two interaction events; its pole at $p^2 = m^2$ says real particles propagate freely, but off-shell ("virtual") excitations contribute too, suppressed by how far $p^2$ is from $m^2$ . Vertices are where the coupling acts — each one costs a factor of the coupling constant, so higher-order diagrams (more vertices) are smaller when the coupling is weak. The total amplitude is the coherent sum of all stories: quantum mechanics adds amplitudes, not probabilities, so diagrams interfere.

Figure

Tree-level and one-loop diagrams for 2→2 scattering in φ⁴ theory. The s-channel tree diagram (gold) dominates; the one-loop bubble is the leading correction.

```svg

```

ExampleTree-level 2→2 scattering in φ⁴ theory

Given: Compute $i\mathcal{M}$ for $\phi\phi \to \phi\phi$ at lowest order (single vertex), coupling $\lambda$ .

Step 1. The leading diagram has one vertex connecting all four external legs. There are no internal propagators and no loops.

Step 2. Apply rule [2]: the single vertex contributes $-i\lambda$ . The symmetry factor is $S=1$ (all four legs are distinguishable external states).

Step 3. Therefore:

5

\quad i\mathcal{M} = -i\lambda \quad\Longrightarrow\quad \mathcal{M} = -\lambda

Derivation hint: With only one vertex and no internal lines, the amplitude is just the vertex factor; no loop integrals occur.

Step 4. The differential cross section in the CM frame for identical bosons:

6

\quad \frac{d\sigma}{d\Omega} = \frac{|\mathcal{M}|^2}{64\pi^2 s} = \frac{\lambda^2}{64\pi^2 s}

Derivation hint: Use the standard $2\to2$ CM cross-section formula $d\sigma/d\Omega = |\mathcal{M}|^2/(64\pi^2 s)$ with $|\mathbf{p}_f|=|\mathbf{p}_i|$ for equal masses.

Boxed answer: $\boxed{\mathcal{M} = -\lambda,\quad \dfrac{d\sigma}{d\Omega} = \dfrac{\lambda^2}{64\pi^2 s}\ \text{[GeV}^{-2}\text{]}}$ where $s$ is the Mandelstam invariant in GeV².

ExampleOne-loop bubble and its divergence

Given: Evaluate the s-channel one-loop bubble correction to $\phi\phi\to\phi\phi$ . External momentum flowing through the loop is $p = p_1 + p_2$ , with $p^2 = s$ .

Step 1. Apply the rules: two vertices give $(-i\lambda)^2$ , the symmetry factor of the bubble is $S = 2$ , and there is one loop integral:

7

\quad i\mathcal{M}_{\text{loop}} = \frac{(-i\lambda)^2}{2}\int\frac{d^4k}{(2\pi)^4}\,\frac{i}{k^2-m^2+i\epsilon}\,\frac{i}{(k-p)^2-m^2+i\epsilon}

Derivation hint: Two propagators share the loop momentum $k$ and $k-p$ ; the factor $1/2$ is the symmetry factor for exchanging the two internal lines.

Step 2. Combine denominators with a Feynman parameter $x$ and shift $k\to k+xp$ :

8

\quad \mathcal{M}_{\text{loop}} = \frac{\lambda^2}{2}\int_0^1 dx\int\frac{d^4k}{(2\pi)^4}\,\frac{1}{[k^2 - \Delta + i\epsilon]^2},\qquad \Delta = m^2 - x(1-x)s

Derivation hint: Use $\frac{1}{AB}=\int_0^1\frac{dx}{[xA+(1-x)B]^2}$ , then complete the square in $k$ .

Step 3. In $d = 4-2\varepsilon$ dimensions (dimensional regularization), the integral gives:

9

\quad \mathcal{M}_{\text{loop}} = \frac{\lambda^2}{32\pi^2}\int_0^1 dx\left[\frac{1}{\varepsilon} - \gamma_E + \ln 4\pi - \ln\frac{\Delta}{\mu^2}\right]

Derivation hint: Wick-rotate to Euclidean space ( $k^0\to ik_E^0$ ) and use the master formula $\int\frac{d^dk}{(2\pi)^d}\frac{1}{(k^2+\Delta)^2}=\frac{\Gamma(2-d/2)}{(4\pi)^{d/2}}\Delta^{d/2-2}$ .

Boxed answer: $\boxed{\mathcal{M}_{\text{loop}} = \dfrac{\lambda^2}{32\pi^2}\left[\dfrac{1}{\varepsilon} + \text{finite}\right]}$ — the $1/\varepsilon$ pole signals a UV divergence absorbed by coupling renormalization.

⚠️ Common mistake

Feynman diagrams are not pictures of particle trajectories. The internal lines are propagators integrated over all internal momenta; "virtual particles" are off-shell ( $p^2 \neq m^2$ ) and are not directly observable. A diagram does not depict what happened — only the full coherent sum of all diagrams has physical meaning, and individual diagrams are gauge-dependent in non-abelian theories.

⚠️ Common mistake

The perturbation series is asymptotic, not convergent. For QED the coefficients eventually grow factorially; summing all orders does not give a finite answer. The series is useful because early terms shrink while $\alpha \approx 1/137$ is small, but it must be truncated. At strong coupling (e.g., low-energy QCD) the diagrammatic expansion fails entirely.

Counterexample

A naïve student might count a single labeled vacuum-polarization bubble and forget the symmetry factor $S=2$ , overcounting the amplitude by a factor of two. Likewise, drawing the "tadpole" diagram and ignoring that it can be cancelled by a counterterm (or vanishes by normal-ordering) leads to spurious infinities. Symmetry factors and the connected/amputated prescriptions are not optional bookkeeping — omitting them yields the wrong cross section.

connection

Feynman diagrams sit at the crossroads of several frameworks. They are the perturbative face of the Path Integral Formulation, where the generating functional $Z[J]$ produces them by functional differentiation. They implement Wick's Theorem graphically and rely on Renormalization to render loop integrals finite. In condensed matter they reappear as the diagrammatics of Many-Body Perturbation Theory (Matsubara/Green's-function techniques). Their UV divergences are organized by the Renormalization Group, and at high orders their factorial growth connects to Resurgence and Trans-series and non-perturbative instanton effects. Modern amplitude methods (BCFW Recursion, the amplituhedron) bypass diagrams entirely while computing the same quantities.

🔬 Where this appears

⏳ History

- 1927–1932 — Dirac formulates QED and predicts antiparticles; Heisenberg & Pauli develop field quantization, but calculations are plagued by infinities.

- 1947 — Lamb–Retherford measure the Lamb shift; Kusch & Foley measure the electron anomalous moment, demanding a finite computational scheme.

- 1947–1948 — Bethe computes a non-relativistic Lamb shift; Schwinger and Tomonaga develop covariant operator renormalization independently.

- 1948 (Pocono Conference) — Feynman presents his diagrammatic spacetime approach; the audience is largely unconvinced.

- 1949 — Feynman publishes "Space-Time Approach to QED" and "The Theory of Positrons," establishing the rules.

- 1949 — Dyson's paper "The Radiation Theories of Tomonaga, Schwinger, and Feynman" proves their equivalence and systematizes the rules to all orders.

- 1950s — Dyson and Salam establish the power-counting theory of renormalizability; the method spreads to all of particle physics.

- 1965 — Feynman, Schwinger, Tomonaga share the Nobel Prize.

- 1972 — 't Hooft and Veltman introduce dimensional regularization, taming non-abelian loop diagrams and enabling the renormalization of the electroweak theory.

open-problem

The perturbative series defined by summing all Feynman diagrams is generically divergent (asymptotic). Its large-order behavior is governed by factorially-growing contributions whose Borel transform has singularities — renormalons and instantons — encoding genuinely non-perturbative physics inaccessible to any finite set of diagrams. Whether and how the full non-perturbative amplitude can be reconstructed from the diagrammatic series (the program of resurgence) remains an active research frontier. Separately, for strongly-coupled theories (low-energy QCD, certain condensed-matter systems) no diagrammatic expansion converges, motivating lattice methods and the gauge/gravity (holographic) correspondence. Even at the level of bookkeeping, the number of diagrams grows factorially with loop order, making high-precision computations a major computational challenge addressed by recursion relations, unitarity cuts, and integration-by-parts identities.

Formal Statement

Why It Holds / Derivation

Building Intuition

Worked Examples

Common Misconceptions

Connections

Real-World Applications

Historical Development

Subtleties & Limits

See Also