The Action Principle as Bridge to Quantum Field Theory

One evening in 1941, at a tavern in Princeton, a visiting physicist named Herbert Jehle told a graduate student about an obscure 1933 paper in which Dirac had remarked that the quantum amplitude for a particle to move between two points "corresponds to" eiS/ℏe^{iS/\hbar} — the exponential of the classical action. The next morning the student, Richard Feynman, worked out what "corresponds to" meant at the blackboard while Jehle watched. It meant equals. Compose that exponential over infinitesimal time steps and the Schrödinger equation falls out. Every possible history contributes an amplitude; the action — the same functional this whole track has been built on — is the phase of each contribution. That discovery is the summit toward which classical mechanics has been climbing since the first page. The action S[ϕ]=∫d4x L(ϕ,∂μϕ)S[\phi] = \int\mathrm{d}^4x\,\mathcal{L}(\phi,\partial_\mu\phi), whose stationarity δS=0\delta S = 0 gave the Euler–Lagrange field equations, turns out to be the exponent of the quantum weight: Z=∫Dϕ eiS[ϕ]/ℏZ = \int\mathcal{D}\phi\,e^{iS[\phi]/\hbar}, a sum over every field history. Classical physics is the stationary-phase saddle of that integral — the one history whose neighbours agree in phase. Meanwhile the Hamiltonian road runs in parallel: the field Poisson brackets you built become equal-time commutators by Dirac's substitution, and the two quantizations agree. Nothing you learned is discarded. Noether's currents become Ward identities; the stress-energy tensor becomes the source of gravity and the generator of translations; a Wick rotation t→−iτt \to -i\tau turns eiSe^{iS} into a Boltzmann weight e−SEe^{-S_E} and quantum field theory into statistical mechanics. The classical action is not a limit of the quantum theory — it is its genetic code.

Key takeaways
  • The classical field action S[ϕ]=∫L d4xS[\phi] = \int\mathcal{L}\,\mathrm{d}^4x, which yields the Euler–Lagrange equations by extremization δS=0\delta S = 0, becomes the central weight of the quantum theory: every field configuration contributes an amplitude eiS/ℏe^{iS/\hbar}.
  • Feynman's path integral sums those amplitudes over all field histories, Z=∫Dϕ eiS[ϕ]/ℏZ = \int\mathcal{D}\phi\,e^{iS[\phi]/\hbar}; the classical equations of motion re-emerge as the stationary-phase saddle in the ℏ→0\hbar \to 0 limit.
  • Observables are time-ordered correlation functions, generated by coupling a source: Z[J]=∫Dϕ eiℏ(S+∫Jϕ)Z[J] = \int\mathcal{D}\phi\,e^{\frac{i}{\hbar}(S + \int J\phi)} and differentiating functionally with respect to JJ.
  • The Hamiltonian route gives the same theory: promote the equal-time field Poisson bracket to the commutator [ϕ^,π^]=iℏ δ3[\hat\phi, \hat\pi] = i\hbar\,\delta^3 — Dirac's prescription { ,}→[ ,]/iℏ\{\ ,\} \to [\ ,]/i\hbar.
  • Classical symmetries survive intact: Noether's conserved currents become Ward–Takahashi operator identities, provided the path-integral measure shares the symmetry (no anomaly).
  • The Wick rotation t→−iτt \to -i\tau converts oscillatory eiSe^{iS} into positive e−SEe^{-S_E}: the integral becomes mathematically definable, numerically computable (lattice QCD), and formally identical to a statistical partition function.
  • Key idea: the classical field action S[ϕ]S[\phi] is the exponent of the quantum amplitude: Z=∫Dϕ eiS[ϕ]/ℏZ = \int\mathcal{D}\phi\,e^{iS[\phi]/\hbar}. Classical dynamics (δS=0\delta S = 0) is the stationary-phase saddle point of this integral, so QFT is built directly on the variational principle of classical field theory.

Where this lives today

Every prediction checked at the Large Hadron Collider is a functional derivative of Z[J]Z[J]. The pipeline is exactly this page's equations run at industrial scale: a classical action (the Standard Model's), correlators by eq. (6), amplitudes by LSZ reduction, organized as Feynman diagrams — millions of them, generated and integrated by computer algebra for each measured process. When experimenters quote a Higgs production cross section to percent precision, they are comparing nature against high-order terms in the expansion of one path integral about its Gaussian saddle. Where the expansion fails — quarks binding into protons — the Euclidean form (8) takes over. Lattice QCD discretizes spacetime into a four-dimensional grid, samples field configurations with probability ∝e−SE\propto e^{-S_E} on exascale supercomputers, and computes hadron masses from first principles to sub-percent accuracy; the proton's mass emerges from a classical action and Monte Carlo statistics. The sharpest current example is the muon's anomalous magnetic moment, where lattice evaluations of the hadronic vacuum polarization now rival data-driven determinations in the comparison that tests the Standard Model at the part-per-billion level. And the statistical bridge runs both ways: the same Euclidean path integral is the native language of critical phenomena in condensed matter, and stationary-phase-with-fluctuations — the mathematics of this page — reappears wherever a system must choose among histories, from semiclassical gravity to the sampling algorithms of machine learning.

The physical picture

Think of every possible history of the field as a contestant. The wild ones — the field leaping to a mountain and back in a femtosecond — are not forbidden; each carries exactly one vote, a unit phasor eiS/ℏe^{iS/\hbar}. Democracy among histories, with a catch: votes are complex numbers. Where neighbouring histories have wildly different actions, their phasors point in scattered directions and cancel. Only near a history where SS is stationary — where first-order wiggles cost no action — do the neighbours' phasors align and add. The classical solution is not chosen; it is what survives the interference of everything else.

Figure
Sum over histories — every path votes with a phasor, and only the stationary neighbourhood survives
start end classical history: δS = 0 every history contributes e^(iS/ħ) — one phasor each near the classical path phases agree — they add far from the classical path phases scatter — they cancel

Planck's constant sets the width of the surviving bundle. When a system's action is enormous compared with ℏ\hbar — a thrown ball racks up 1034 ℏ10^{34}\,\hbar of phase — the constructive band around the classical history is razor-thin and the motion looks perfectly deterministic. When S∼ℏS \sim \hbar — an electron in an atom — the band is wide, many histories contribute coherently, and "the trajectory" stops existing. The classical action is literally the ruler against which quantumness is measured, and the ℏ→0\hbar \to 0 limit is not a different theory but a stationary-phase approximation to this one.

Figure
The ħ ruler — rapid phase oscillation cancels everything except the stationary band
stationary band: δS = 0, phases agree rapid oscillation — neighbouring histories cancel in pairs smaller ħ ⇒ faster oscillation ⇒ narrower band history label Re e^(iS/ħ)

The same bridge can be crossed on the Hamiltonian deck. There, the classical structure that survives is not the action but the bracket: the equal-time field Poisson brackets {ϕ(x),π(y)}=δ3(x−y)\{\phi(\mathbf{x}), \pi(\mathbf{y})\} = \delta^3(\mathbf{x}-\mathbf{y}) become commutators [ϕ^,π^]=iℏ δ3[\hat\phi, \hat\pi] = i\hbar\,\delta^3, and the Noether charges that generated symmetry flows through Poisson brackets become operators generating the same transformations through commutators. Two roads — Lagrangian sum-over-histories and Hamiltonian operator algebra — one quantum field theory. Everything this track built is cargo that makes the crossing.

From least action to sum over histories

Start from the object itself. A field theory is specified by one functional of the field configuration — all of dynamics, all of symmetry, compressed into a single spacetime integral:

S[ϕ]=∫d4x  L(ϕ, ∂μϕ)(1)S[\phi] = \int \mathrm{d}^4x\;\mathcal{L}\big(\phi,\,\partial_\mu\phi\big) \qquad (1)

where: L\mathcal{L} — the Lagrangian density; ϕ\phi — the field (a scalar here; vectors and spinors generalize); d4x=dt d3x\mathrm{d}^4x = \mathrm{d}t\,\mathrm{d}^3x. Note what kind of object SS is: a machine that eats an entire field history — a function on all of spacetime — and returns one number. Classical physics asks where that number is stationary; quantum physics, we are about to see, uses the number as a phase. The classical story in two moves: vary, integrate by parts. Under ϕ→ϕ+δϕ\phi \to \phi + \delta\phi with δϕ\delta\phi vanishing on the boundary, the chain rule gives δS=∫d4x (∂L∂ϕδϕ+∂L∂(∂μϕ)∂μδϕ)\delta S = \int\mathrm{d}^4x\,\big(\tfrac{\partial\mathcal{L}}{\partial\phi}\delta\phi + \tfrac{\partial\mathcal{L}}{\partial(\partial_\mu\phi)}\partial_\mu\delta\phi\big); moving the derivative off δϕ\delta\phi by parts and demanding stationarity for arbitrary variations yields the field equation:

δSδϕ(x)=∂L∂ϕ−∂μ∂L∂(∂μϕ)=0(2)\frac{\delta S}{\delta\phi(x)} = \frac{\partial\mathcal{L}}{\partial\phi} - \partial_\mu\frac{\partial\mathcal{L}}{\partial(\partial_\mu\phi)} = 0 \qquad (2)

where: δS/δϕ\delta S/\delta\phi — the functional derivative, the field theorist's gradient. For the Klein–Gordon Lagrangian L=12(∂ϕ)2−12m2ϕ2\mathcal{L} = \tfrac12(\partial\phi)^2 - \tfrac12 m^2\phi^2 this is (□+m2)ϕ=0(\Box + m^2)\phi = 0. Keep the left-hand form in view: "δS/δϕ=0\delta S/\delta\phi = 0" is about to be reinterpreted as the location of a saddle point.

Derivation AThe propagator becomes exp(iS/ℏ) — time slicing

We prove the bridge equation (3): the quantum propagator ⟨qf∣e−iH^T/ℏ∣qi⟩\langle q_f|e^{-i\hat H T/\hbar}|q_i\rangle equals a sum over all histories weighted by eiS/ℏe^{iS/\hbar}, with SS the classical action itself. The cleanest proof is Feynman's own time-slicing argument, run for a single degree of freedom — fields follow at the end by attaching a coordinate to every point of space, since the slicing never cares how many coordinates there are. The route: chop the evolution into NN infinitesimal steps, insert position completeness between them, do one Gaussian momentum integral per step, and watch the exponents add up into ∫L dt\int L\,\mathrm{d}t. The step to watch is the Gaussian integral of Step 3, where completing the square turns pq˙−Hp\dot q - H into the Lagrangian — the Legendre transform working in reverse inside a quantum matrix element; miss that and the appearance of SS looks like magic.

Step 1. Split the quantum evolution into NN tiny steps: e−iH^T/ℏ=(e−iH^ϵ/ℏ)Ne^{-i\hat H T/\hbar} = \big(e^{-i\hat H\epsilon/\hbar}\big)^N with ϵ=T/N\epsilon = T/N. Why: no one can evaluate a finite-time propagator directly, but an infinitesimal one can be expanded to first order in ϵ\epsilon — the whole derivation is "divide and conquer" in time.

Step 2. Insert a complete set of position states 1=∫dqk ∣qk⟩⟨qk∣1 = \int\mathrm{d}q_k\,|q_k\rangle\langle q_k| between every pair of factors. Why: this converts one quantum matrix element into an integral over all intermediate positions q1,…,qN−1q_1, \dots, q_{N-1} — that is, over all discretized paths. The "sum over histories" enters here, as pure linear algebra.

Step 3. Evaluate a single step by doing the Gaussian momentum integral:

⟨qk+1∣e−iH^ϵ/ℏ∣qk⟩∝exp⁡ ⁣(iϵℏ[m2(qk+1−qkϵ)2−V(qk)])\langle q_{k+1}|e^{-i\hat H\epsilon/\hbar}|q_k\rangle \propto \exp\!\Big(\frac{i\epsilon}{\hbar}\Big[\frac{m}{2}\Big(\frac{q_{k+1}-q_k}{\epsilon}\Big)^2 - V(q_k)\Big]\Big)

Why the bracket is the Lagrangian: the momentum integral is Gaussian because H^=p^2/2m+V\hat H = \hat p^2/2m + V, and completing the square turns pq˙−Hp\dot q - H into m2q˙2−V=L\tfrac{m}{2}\dot q^2 - V = L evaluated on the segment. The Legendre transformation of this track's Hamiltonian chapters is doing the work, in reverse, inside a quantum integral.

Step 4. Multiply the NN factors and take the continuum limit: the exponents add, ∑kϵ L(qk,q˙k)→∫0TL dt=S[q]\sum_k \epsilon\,L(q_k, \dot q_k) \to \int_0^T L\,\mathrm{d}t = S[q], giving

⟨qf∣e−iH^T/ℏ∣qi⟩=∫Dq  eiS[q]/ℏ(3)\langle q_f|e^{-i\hat H T/\hbar}|q_i\rangle = \int\mathcal{D}q\;e^{iS[q]/\hbar} \qquad (3)

where: Dq\mathcal{D}q — the limit of ∏kdqk\prod_k\mathrm{d}q_k with its normalization. Why this is the theorem: the exponent assembled itself into the classical action — not an analogue of it, the very functional whose stationarity is Newton's law.

Step 5. Generalize to fields: replace the single coordinate q(t)q(t) by one coordinate per point of space, ϕ(x,t)\phi(\mathbf{x}, t), so S[q]→S[ϕ]=∫d4x LS[q] \to S[\phi] = \int\mathrm{d}^4x\,\mathcal{L} and Dq→Dϕ\mathcal{D}q \to \mathcal{D}\phi. Why nothing new is needed: a field is a mechanical system with a continuum of degrees of freedom — exactly the continuum limit built earlier in this track — and the slicing argument never cared how many coordinates there were. *Verified: standard derivation (Feynman & Hibbs; Peskin & Schroeder Ch. 9). Verifier should confirm the Gaussian integral over momenta yields the kinetic term and that the exponent assembles to SS.*

The generating machine

What does one compute with the field path integral? The observables of QFT are time-ordered correlation functions — vacuum expectation values of products of field operators — and the path integral produces them as moments of the weight eiS/ℏe^{iS/\hbar}:

⟨Ω∣ T ϕ^(x1)⋯ϕ^(xn) ∣Ω⟩=∫Dϕ  ϕ(x1)⋯ϕ(xn)  eiS[ϕ]/ℏ∫Dϕ  eiS[ϕ]/ℏ(4)\langle\Omega|\,T\,\hat\phi(x_1)\cdots\hat\phi(x_n)\,|\Omega\rangle = \frac{\int\mathcal{D}\phi\;\phi(x_1)\cdots\phi(x_n)\;e^{iS[\phi]/\hbar}}{\int\mathcal{D}\phi\;e^{iS[\phi]/\hbar}} \qquad (4)

where: TT — time ordering; ∣Ω⟩|\Omega\rangle — the interacting vacuum. The time ordering, an awkward operator instruction on the left, is automatic on the right — the sliced integral assembles operators in time order by construction. This formula is the working definition of a QFT: give me SS, and correlators (hence masses, cross sections, decay rates via LSZ reduction) follow. Rather than computing each correlator separately, couple the field to an external source J(x)J(x) — a crank handle welded onto the action:

Z[J]=∫Dϕ  exp⁡ ⁣(iℏ[S[ϕ]+∫d4x J(x) ϕ(x)])(5)Z[J] = \int\mathcal{D}\phi\;\exp\!\Big(\frac{i}{\hbar}\Big[S[\phi] + \int\mathrm{d}^4x\,J(x)\,\phi(x)\Big]\Big) \qquad (5)

where: JJ — a classical auxiliary field, set to zero at the end. Each functional derivative δ/δJ(x)\delta/\delta J(x) pulls one factor of ϕ(x)\phi(x) down from the exponent, so the whole tower of correlators lives inside one functional:

⟨0∣ T ϕ^(x1)⋯ϕ^(xn) ∣0⟩=1Z[0] 1in δnZ[J]δJ(x1)⋯δJ(xn)∣J=0(6)\langle 0|\,T\,\hat\phi(x_1)\cdots\hat\phi(x_n)\,|0\rangle = \frac{1}{Z[0]}\,\frac{1}{i^n}\,\frac{\delta^n Z[J]}{\delta J(x_1)\cdots\delta J(x_n)}\Big|_{J=0} \qquad (6)

where: the 1/in1/i^n tidies the phases. Z[J]Z[J] is to field theory what a partition function with external fields is to statistical mechanics — differentiate to interrogate. For a free (quadratic) action the integral is Gaussian and closes in one line, Z[J]=Z[0]exp⁡(−12∫J ΔF J)Z[J] = Z[0]\exp\big(-\tfrac12\int J\,\Delta_F\,J\big): all of free-field theory is the Feynman propagator ΔF\Delta_F, and perturbation theory (Feynman diagrams) is the systematic expansion of the interacting Z[J]Z[J] about that Gaussian.

Go deeper · Derivation B · 2 minStationary phase: the classical world as a saddle point
Derivation BStationary phase: the classical world as a saddle point

Now we run the bridge backwards and derive classical field theory from the quantum integral: the goal is to show that when the action is large compared with ℏ\hbar, the only configurations surviving the path integral are those obeying the Euler–Lagrange equation (2). The tool is the stationary-phase argument — functionally Taylor-expand SS about a trial configuration and examine how neighbouring histories interfere. Watch the logic of the limit: ℏ\hbar never literally goes to zero; what matters is the dimensionless ratio S/ℏS/\hbar for the phenomenon at hand, the integral's only knob. Watch also what happens to the quadratic term — it is not discarded but becomes the fluctuation determinant, the first quantum correction wrapped around the classical skeleton.

Step 1. Write the amplitude Z=∫Dϕ eiS[ϕ]/ℏZ = \int\mathcal{D}\phi\,e^{iS[\phi]/\hbar} and ask which configurations dominate when ℏ\hbar is small compared with the actions in play. Why this is the right question: ℏ\hbar never literally "goes to zero" — what varies is the size of SS for the phenomenon at hand, and S/ℏS/\hbar is the only dimensionless knob in the weight.

Step 2. Expand the action functionally about an arbitrary configuration ϕcl\phi_{cl} with fluctuation η\eta:

S[ϕcl+η]=S[ϕcl]+∫d4x δSδϕ∣cl η+12∫η δ2Sδϕ δϕ η+⋯S[\phi_{cl} + \eta] = S[\phi_{cl}] + \int\mathrm{d}^4x\,\frac{\delta S}{\delta\phi}\Big|_{cl}\,\eta + \frac12\int \eta\,\frac{\delta^2S}{\delta\phi\,\delta\phi}\,\eta + \cdots

Why a Taylor expansion is the whole tool: the integral cannot be done, but its local phase structure is entirely controlled by the first two variations.

Step 3. Apply the stationary-phase argument. If the linear term is nonzero, then along that direction the phase S/ℏS/\hbar sweeps through many cycles of 2π2\pi for tiny η\eta, and contributions cancel in pairs. Only where the linear term vanishes — δS/δϕ∣cl=0\delta S/\delta\phi|_{cl} = 0 — do neighbouring histories interfere constructively. Why the cancellation is ruthless for large S/ℏS/\hbar: a macroscopic action of 1034 ℏ10^{34}\,\hbar means a field wiggle of one part in 103410^{34} already flips the phasor — the non-stationary integrand is the most violently oscillating object in physics.

Step 4. Keep the leading terms:

Z≈eiS[ϕcl]/ℏ×(Gaussian fluctuation determinant)Z \approx e^{iS[\phi_{cl}]/\hbar} \times \big(\text{Gaussian fluctuation determinant}\big)

Why the determinant matters too: the quadratic term gives the one-loop quantum correction — the first whisper of the quantum theory around the classical skeleton; higher orders generate the full loop expansion.

Step 5. Read the result: as ℏ/S→0\hbar/S \to 0 only ϕcl\phi_{cl} survives, and ϕcl\phi_{cl} obeys exactly the Euler–Lagrange equation (2). Why this closes the loop of the entire track: the variational principle that was postulated back in Lagrangian mechanics is now a theorem — stationarity is where quantum interference permits a world to be visible. *Verified: reproduces Peskin & Schroeder's classical limit and the standard WKB/stationary-phase argument; verifier should reconstruct the functional Taylor expansion and confirm the linear term gives the EL equation.*

The Hamiltonian road and the Euclidean detour

Canonical quantization crosses the same bridge with Hamiltonian cargo. The classical field theory built the equal-time Poisson bracket {ϕ(x),π(y)}PB=δ3(x−y)\{\phi(\mathbf{x}), \pi(\mathbf{y})\}_{PB} = \delta^3(\mathbf{x}-\mathbf{y}) from the Hamiltonian formulation of fields; Dirac's prescription { , }PB→1iℏ[ , ]\{\ ,\ \}_{PB} \to \tfrac{1}{i\hbar}[\ ,\ ] promotes it to the operator statement

[ϕ^(x,t), π^(y,t)]=iℏ δ3(x−y)(7)[\hat\phi(\mathbf{x},t),\,\hat\pi(\mathbf{y},t)] = i\hbar\,\delta^3(\mathbf{x}-\mathbf{y}) \qquad (7)

where: π^=∂L/∂ϕ˙\hat\pi = \partial\mathcal{L}/\partial\dot\phi — the conjugate momentum density, exactly the object defined classically. This single commutator, imposed on each spatial point, generates the entire operator structure of free-field theory — mode expansions, creation and annihilation operators, particles. That the path integral (3) and the commutator (7) define the same theory is one of the standard equivalences of QFT; each formulation keeps a different classical structure manifest — the action and its symmetries on one side, the bracket algebra and unitarity on the other. One technical confession converts into profound physics. The Minkowski weight eiS/ℏe^{iS/\hbar} is a pure phase — the integral is oscillatory and not absolutely convergent. Rotate time, t→−iτt \to -i\tau, and the action acquires a sign-definite Euclidean form:

Z=∫Dϕ  e−SE[ϕ]/ℏ,SE=∫d4xE [12(∂Eϕ)2+V(ϕ)](8)Z = \int\mathcal{D}\phi\;e^{-S_E[\phi]/\hbar}, \qquad S_E = \int\mathrm{d}^4x_E\,\Big[\tfrac12(\partial_E\phi)^2 + V(\phi)\Big] \qquad (8)

where: SES_E — the Euclidean action, positive for stable theories. Three payoffs at once: the integral becomes a genuine (Wiener-like) measure, mathematically defensible; the weight e−SE/ℏe^{-S_E/\hbar} is a Boltzmann factor, so QFT in dd spacetime dimensions maps onto equilibrium statistical mechanics in dd spatial dimensions — with ℏ\hbar playing temperature; and the positive weight can be sampled by Monte Carlo, which is the entire basis of lattice field theory. Rare classical solutions of the Euclidean equations — instantons — carry finite SES_E and contribute e−SE/ℏe^{-S_E/\hbar}: quantum tunneling computed from a classical action.

Go deeper · Derivation C · 2 minNoether's theorem becomes a Ward identity
Derivation CNoether's theorem becomes a Ward identity

The last piece of cargo is symmetry: we derive the Ward identity (9), the quantum completion of Noether's theorem. Classically, an invariance of SS produced a conserved current; here the same manipulation performed inside the integral produces an operator identity, with contact terms as the new quantum content. The trick is to promote the global symmetry parameter to a local probe α(x)\alpha(x) and treat the transformation as a change of integration variable: renaming a dummy variable cannot change any integral, so the first-order response must vanish identically — that triviality, localized, is the whole theorem. Watch where each piece comes from — the current divergence from the action's response, the delta-function contact terms from the transformation hitting the inserted fields — and watch the assumption smuggled into Step 2, invariance of the measure Dϕ\mathcal{D}\phi, whose failure is not a technicality but the anomaly.

Step 1. Let the action be invariant under the global transformation ϕ→ϕ+α Δϕ\phi \to \phi + \alpha\,\Delta\phi with constant α\alpha, with classical Noether current jμ=∂L∂(∂μϕ)Δϕj^\mu = \tfrac{\partial\mathcal{L}}{\partial(\partial_\mu\phi)}\Delta\phi. Why start global: this is precisely the hypothesis of the classical Noether theorem — nothing quantum yet.

Step 2. Inside the path integral, perform the local version ϕ→ϕ+α(x) Δϕ\phi \to \phi + \alpha(x)\,\Delta\phi as a change of integration variable, assuming the measure Dϕ\mathcal{D}\phi is invariant. Why localizing is the trick: a change of dummy variable cannot change the value of any integral — so whatever the transformation produces must sum to zero. The spacetime-dependent α(x)\alpha(x) is the probe that converts that triviality into a local statement.

Step 3. Compute the first-order change of the action. Since constant α\alpha is a symmetry, only the gradient survives:

δS=∫d4x (∂μα) jμ=−∫d4x α(x) ∂μjμ\delta S = \int\mathrm{d}^4x\,(\partial_\mu\alpha)\,j^\mu = -\int\mathrm{d}^4x\,\alpha(x)\,\partial_\mu j^\mu

Why the current appears: this is the same integration-by-parts bookkeeping as the classical Noether derivation — the quantum context has changed nothing about the algebra.

Step 4. Apply the same shift inside a correlator with field insertions ϕ(y1)⋯ϕ(yn)\phi(y_1)\cdots\phi(y_n): invariance of the renamed integral, order by order in α(x)\alpha(x), forces

∂μ⟨jμ(x) ϕ(y1)⋯ϕ(yn)⟩=−iℏ∑aδ4(x−ya) ⟨ϕ(y1)⋯Δϕ(ya)⋯⟩(9)\partial_\mu\big\langle j^\mu(x)\,\phi(y_1)\cdots\phi(y_n)\big\rangle = -i\hbar\sum_a \delta^4(x - y_a)\,\big\langle \phi(y_1)\cdots\Delta\phi(y_a)\cdots\big\rangle \qquad (9)

where: the right-hand contact terms come from the transformation acting on the inserted fields themselves. Why the delta functions are physics, not debris: away from insertions the current is conserved as an operator equation; at insertions the contact terms say the conserved charge generates the symmetry transformation on every field — the quantum completion of Noether's theorem.

Step 5. Note the loophole deliberately left open: Step 2 assumed the measure is invariant. When it is not — as for the chiral symmetry of fermions — the identity picks up an extra term, an anomaly: a symmetry of the classical action that quantization refuses to honour. The decay π0→γγ\pi^0 \to \gamma\gamma proceeds at its observed rate precisely because of this term. Verified: matches the standard Ward–Takahashi derivation (Peskin & Schroeder Ch. 9.6, Weinberg Vol. 1 Ch. 10). Verifier should reproduce jμj^\mu and check the contact terms arise from δϕ(ya)\delta\phi(y_a).

Figure
One functional, four destinies — the classical action as the genetic code of quantum field theory
the classical action S[φ] = ∫d⁴x 𝓛 symmetries included classical field theory δS = 0 → Euler–Lagrange (the stationary saddle) path integral Z[J] = ∫𝒟φ e^(i(S+∫Jφ)/ħ) correlators → LHC amplitudes canonical quantization {φ, π} → [φ̂, π̂] = iħδ³ operators, particles Euclidean / lattice t → −iτ: Z = ∫𝒟φ e^(−S_E/ħ) stat mech, Monte Carlo, instantons vary exponentiate Legendre + brackets Wick rotate symmetries of S: Noether currents (classical) ⇄ Ward identities (quantum) one functional, four destinies

Which formula answers which question

you wantusenotes
classical field equationsδS/δϕ=0\delta S/\delta\phi = 0, eq. (2)the stationary-phase saddle
quantum amplitudeZ=∫Dϕ eiS/ℏZ = \int\mathcal{D}\phi\,e^{iS/\hbar}, eq. (3)sum over histories
all observablesZ[J]Z[J] and its δ/δJ\delta/\delta J, eqs. (5)–(6)correlators → LSZ → cross sections
operator formulation[ϕ^,π^]=iℏδ3[\hat\phi,\hat\pi] = i\hbar\delta^3, eq. (7)Dirac's bracket promotion
a well-defined integral / numericsWick-rotated e−SE/ℏe^{-S_E/\hbar}, eq. (8)lattice Monte Carlo, instantons
symmetry consequencesWard identity, eq. (9)Noether's theorem, quantized

Limiting cases

CaseResultMeaning
classical limit S≫ℏS \gg \hbarδSδϕ∣ϕcl=0\dfrac{\delta S}{\delta\phi}\Big|_{\phi_{cl}} = 0 dominatesstationary phase: classical field theory recovered exactly, with eiS[ϕcl]/ℏ×e^{iS[\phi_{cl}]/\hbar} \times fluctuations
free field (Gaussian action)Z[J]=Z[0]exp⁡(−12∫J ΔF J)Z[J] = Z[0]\exp\big(-\tfrac12\int J\,\Delta_F\,J\big)exactly solvable; all correlators from the propagator via Wick's theorem
Euclidean / statistical limitZ=∫Dϕ e−SE/ℏZ = \int\mathcal{D}\phi\,e^{-S_E/\hbar}a partition function: QFT ↔ equilibrium statistical mechanics
single-particle limitZ=∫Dx(t) eiS[x]/ℏZ = \int\mathcal{D}x(t)\,e^{iS[x]/\hbar}drop the spatial continuum: Feynman's original quantum-mechanical path integral
tunneling / instanton sectoramplitude ∼e−SE[ϕinst]/ℏ\sim e^{-S_E[\phi_{inst}]/\hbar}Euclidean saddles govern barrier penetration and vacuum decay

Worked examples

1 · The ħ ruler: how classical is a thrown ball, or an atomic electron? (confidence builder)

Estimate S/ℏS/\hbar for (a) a ball of mass 0.10 kg0.10\ \mathrm{kg} thrown at 10 m/s10\ \mathrm{m/s} for 1.0 s1.0\ \mathrm{s}, and (b) an electron of kinetic energy 1.0 eV1.0\ \mathrm{eV} crossing a 1.0 nm1.0\ \mathrm{nm} molecular junction.

Step 1. For the ball, the action along the real path is of order the kinetic energy times the time: S∼12mv2 t=12(0.10)(10)2(1.0)=5.0 J sS \sim \tfrac12 mv^2\,t = \tfrac12(0.10)(10)^2(1.0) = 5.0\ \mathrm{J\,s}.

Step 2. Divide by the quantum: S/ℏ=5.0/(1.055×10−34)=4.7×1034S/\hbar = 5.0/(1.055\times10^{-34}) = \boxed{4.7\times10^{34}} radians of phase. A neighbouring history differing by one part in 103410^{34} already flips its phasor — the constructive band around the classical trajectory is unimaginably thin. This ball is classical beyond any conceivable measurement.

Step 3. For the electron: speed v=2E/m=2(1.602×10−19)/(9.109×10−31)=5.93×105 m/sv = \sqrt{2E/m} = \sqrt{2(1.602\times10^{-19})/(9.109\times10^{-31})} = 5.93\times10^{5}\ \mathrm{m/s}; crossing time t=d/v=(1.0×10−9)/(5.93×105)=1.69×10−15 st = d/v = (1.0\times10^{-9})/(5.93\times10^{5}) = 1.69\times10^{-15}\ \mathrm{s}.

Step 4. Action: S∼E t=(1.602×10−19)(1.69×10−15)=2.7×10−34 J sS \sim E\,t = (1.602\times10^{-19})(1.69\times10^{-15}) = 2.7\times10^{-34}\ \mathrm{J\,s}, so S/ℏ=≈2.6S/\hbar = \boxed{\approx 2.6} — a handful of radians. The "sum over histories" is now genuinely plural: many histories interfere, the electron has no trajectory, and full quantum treatment is mandatory.

Step 5. The moral in one line: classicality is not a property of objects but of the ratio S/ℏS/\hbar — the same formalism covers both regimes, with stationary phase doing more or less of the work.

2 · A Euclidean saddle you can measure: tunneling in an STM (exam level)

A scanning tunneling microscope holds its tip a distance d=0.50 nmd = 0.50\ \mathrm{nm} from a surface; the vacuum gap presents an effective barrier of height V=4.5 eVV = 4.5\ \mathrm{eV} (a typical work function) to electrons at the Fermi level. Use the Euclidean-action/WKB weight e−2SE/ℏe^{-2S_E/\hbar} with SE=∫0d2mV dxS_E = \int_0^d\sqrt{2mV}\,\mathrm{d}x to estimate the tunneling suppression, and the sensitivity of the current to a 0.10 nm0.10\ \mathrm{nm} height change.

Step 1. The under-barrier momentum scale: 2mV=2(9.109×10−31)(4.5×1.602×10−19)=1.313×10−48=1.146×10−24 kg m/s\sqrt{2mV} = \sqrt{2(9.109\times10^{-31})(4.5\times1.602\times10^{-19})} = \sqrt{1.313\times10^{-48}} = 1.146\times10^{-24}\ \mathrm{kg\,m/s}.

Step 2. Euclidean action across the gap: SE=2mV d=(1.146×10−24)(0.50×10−9)=5.73×10−34 J s=5.43 ℏS_E = \sqrt{2mV}\,d = (1.146\times10^{-24})(0.50\times10^{-9}) = 5.73\times10^{-34}\ \mathrm{J\,s} = 5.43\,\hbar. This is an honest classical action — of a solution of the Euclidean equations of motion — and it is a few ℏ\hbar: the deep-quantum regime where saddle-point tunneling is both important and accurately exponential.

Step 3. Tunneling weight: e−2SE/ℏ=e−10.9=≈1.9×10−5e^{-2S_E/\hbar} = e^{-10.9} = \boxed{\approx 1.9\times10^{-5}}. Only about two electrons in a hundred thousand that strike the barrier get through — yet at femtoamp-to-nanoamp scales that is a robust, measurable current.

Step 4. Sensitivity: increasing dd by 0.10 nm0.10\ \mathrm{nm} adds Δ(2SE/ℏ)=22mV Δd/ℏ=2(1.146×10−24)(1.0×10−10)/(1.055×10−34)=2.2\Delta(2S_E/\hbar) = 2\sqrt{2mV}\,\Delta d/\hbar = 2(1.146\times10^{-24})(1.0\times10^{-10})/(1.055\times10^{-34}) = 2.2, so the current drops by e−2.2≈a factor of 9e^{-2.2} \approx \boxed{\text{a factor of } 9}. One ångström costs nearly an order of magnitude of current — which is precisely why an STM resolves individual atoms: the exponential of a classical Euclidean action is the sharpest ruler in instrumentation.

Step 5. Connect back to the formalism: what we used is the n=0n = 0 instanton-type estimate ∣amplitude∣2∼e−2SE/ℏ|{\text{amplitude}}|^2 \sim e^{-2S_E/\hbar} from eq. (8) — a saddle point of the Wick-rotated path integral. The same logic, scaled up to field space, gives false-vacuum decay rates and the θ\theta-vacuum structure of QCD.

Origins · 1 minRichard Feynman and how this was discovered
RF
Richard Feynman
1918–1988 · Princeton, Cornell, and Caltech

Feynman turned Dirac's throwaway remark into a formulation of quantum mechanics. Dirac's 1933 paper observed that eiS/ℏe^{iS/\hbar} "corresponds to" the short-time propagator; when Herbert Jehle showed Feynman the paper in a Princeton tavern in 1941, Feynman computed at the blackboard the next day and found the correspondence was an equality — composing infinitesimal steps of eiL ϵ/ℏe^{iL\,\epsilon/\hbar} reproduced the Schrödinger equation exactly. ("Dirac meant they were proportional," he later said; "I found they were equal.") His 1942 thesis, The Principle of Least Action in Quantum Mechanics, and the 1948 Reviews of Modern Physics paper built quantum theory on the sum over histories; the diagrams that organize its perturbative expansion made him, with Schwinger and Tomonaga, a 1965 Nobel laureate for quantum electrodynamics. The path integral was slow to be accepted — Feynman could not at first convince even Bohr — but it became the working language of modern field theory, from the Standard Model to lattice QCD to string theory.

History
1933
— Dirac publishes "The Lagrangian in Quantum Mechanics", noting that the quantum propagator over a short time "corresponds to" eiS/ℏe^{iS/\hbar} — and leaves it there.
1942
— Feynman's Princeton thesis turns the correspondence into an equality: quantum mechanics as a sum over histories weighted by the classical action.
1948
— The path-integral formulation appears in Reviews of Modern Physics; within two years the diagrammatic expansion revolutionizes QED.
1950
— Ward derives the identities linking symmetries to quantum correlators; Takahashi generalizes them in 1957 — Noether's theorem completes its quantum crossing.
1974
— Wilson formulates lattice gauge theory: the Euclidean path integral becomes statistical mechanics on a grid, opening QCD to computation and clarifying confinement.

Going further

Texts. Morin, Introduction to Classical Mechanics, Ch. 6 (variational principle) · Goldstein, Poole & Safko, Classical Mechanics, 3rd ed., Ch. 13 (continuous systems and fields) · Landau & Lifshitz, Mechanics, §2 (principle of least action) · Peskin & Schroeder, An Introduction to QFT, Ch. 9 (functional methods) · Weinberg, The Quantum Theory of Fields, Vol. 1, Ch. 9 · Zee, QFT in a Nutshell, Ch. I.2 (path integral formulation).

Provenance. Equations and structure follow standard references: classical action and Euler–Lagrange from Goldstein Ch. 13 and Landau–Lifshitz Vol. 1 §2; path integral, generating functional, stationary phase, Ward identities from Peskin & Schroeder Ch. 9, Weinberg Vol. 1 Ch. 9–10, and Zee Ch. I. A verifier should recompute: (1) the EL equation from δS=0\delta S=0; (2) the time-sliced derivation of eiS/ℏe^{iS/\hbar} including the momentum Gaussian; (3) the stationary-phase saddle giving the classical solution; (4) the Ward identity from a local field redefinition. All numerical constants (ℏ\hbar, ii, delta functions) are standard textbook values; no nonstandard numbers introduced. Worked-example estimates on this page (S/ℏ=4.7×1034S/\hbar = 4.7\times10^{34} for the ball and ≈2.6\approx 2.6 for the 1 eV electron; SE=5.43 ℏS_E = 5.43\,\hbar, suppression 1.9×10−51.9\times10^{-5}, and factor ≈9\approx 9 per ångström for the STM barrier) are order-of-magnitude WKB arithmetic recomputed from the stated inputs.

Builds on: The Lagrangian Density and Field Action · Euler-Lagrange Equations for Fields · Functional Derivatives and Field Variations · Hamiltonian Formulation for Fields · Poisson Brackets for Fields. Leads to: Noether's Theorem for Continuous Symmetries · Conserved Currents and Charges · The Canonical Stress-Energy Tensor · Internal Symmetries and Conserved Charge.

Does the path integral mean the particle really takes every path?

The formalism assigns an amplitude eiS/ℏe^{iS/\hbar} to every history and adds them; only the total is physical. Whether the particle "really" traverses each path is interpretation, not prediction — the same theory emerges from canonical quantization, where [ϕ^,π^]=iℏ δ3[\hat\phi, \hat\pi] = i\hbar\,\delta^3 generates everything and no path language appears at all. What is objective: when S∼ℏS \sim \hbar, no single trajectory suffices. The 1 eV electron with S/ℏ≈2.6S/\hbar \approx 2.6 genuinely has no trajectory — many histories contribute coherently, and their interference has measurable consequences. When S/ℏ∼1034S/\hbar \sim 10^{34}, the constructive band is so thin that the classical path is all anything could resolve. Safest reading: a sum over histories is how the amplitude is computed; nature only ever shows you the sum.

Every history has the same weight magnitude, so how does the classical path win?

It doesn't win by carrying a bigger weight — ∣eiS/ℏ∣=1|e^{iS/\hbar}| = 1 for the wildest history and the classical one alike. It wins by having agreeable neighbours. Near a stationary configuration, first-order wiggles cost no action, so the phasors of a whole neighbourhood point the same way and add; far from stationarity, a tiny deformation sweeps the phase through many cycles of 2π2\pi and contributions cancel in pairs. As the physical picture puts it, the classical solution is not chosen — it is what survives the interference of everything else. So resist the phrase "most probable path": no single path has a probability at all. What survives is the coherent contribution of the stationary band, whose width is set by ℏ\hbar.

How can ℏ→0\hbar \to 0 when ℏ\hbar is a fixed constant of nature?

It's shorthand, and the stationary-phase derivation is careful about this: ℏ\hbar never literally goes to zero — what varies is the size of the action for the phenomenon at hand, and S/ℏS/\hbar is the only dimensionless knob in the weight eiS/ℏe^{iS/\hbar}. A thrown ball racks up 4.7×10344.7\times10^{34} radians of phase, so a neighbouring history differing by one part in 103410^{34} already flips its phasor — the classical limit in full force, at the real value of ℏ\hbar. The 1 eV electron crossing a nanometre accumulates only a couple of radians, so many histories interfere. Same formalism, same constant; the stationary-phase approximation simply does more or less of the work. Classicality is a property of the ratio S/ℏS/\hbar, not of objects.

Is the Wick rotation t→−iτt \to -i\tau a physical claim about imaginary time, or just a trick?

It's a controlled analytic continuation, not a claim that time is secretly imaginary. Correlation functions are computed where the integral is well-behaved — in Euclidean signature, where e−SE/ℏe^{-S_E/\hbar} is a positive, Wiener-like weight — and then continued back to Minkowski, where the physics lives. But the detour earns more than convergence, which is why it feels deeper than a trick: the Euclidean weight is literally a Boltzmann factor, making QFT in dd spacetime dimensions formally identical to statistical mechanics in dd spatial dimensions; lattice QCD computes the proton mass by Monte Carlo sampling of e−SEe^{-S_E}; and Euclidean saddles — instantons — give real, measured tunneling rates ∼e−2SE/ℏ\sim e^{-2S_E/\hbar}, like the STM current in the worked example. The rotated theory is a different, rigorous handle on the same physics.

How can renaming a dummy integration variable produce an anomaly?

The renaming itself is exactly as trivial as it looks — that triviality is the engine of the Ward identity. What can fail is the hidden assumption in Step 2 of the derivation: that only the action responds to the transformation, while the measure Dϕ\mathcal{D}\phi sits inert. But Dϕ\mathcal{D}\phi is not a product of ordinary measures; it only means something after regularization, and for some symmetries — chiral rotations of fermions are the canonical case — no regulator preserves the classical invariance. The transformation then has a nontrivial Jacobian, and the Ward identity picks up an extra term: an anomaly, a symmetry of the classical action that quantization refuses to honour. This is physics, not pathology — the decay π0→γγ\pi^0 \to \gamma\gamma proceeds at its observed rate precisely because of that term.

If large S/ℏS/\hbar makes things classical, why do electrons still tunnel in an STM?

Because the classical limit suppresses tunneling only exponentially — never to zero. The rate is governed by the Euclidean action of the barrier crossing, ∣amplitude∣2∼e−2SE/ℏ|\text{amplitude}|^2 \sim e^{-2S_E/\hbar}, and everything depends on how many ℏ\hbar's that saddle costs. In the STM example, SE=2mV d=5.43 ℏS_E = \sqrt{2mV}\,d = 5.43\,\hbar, giving suppression e−10.9≈1.9×10−5e^{-10.9} \approx 1.9\times10^{-5} — tiny, but robustly measurable at nanoamp scales, and so sharp that one extra ångström of gap costs a factor of nine in current. For your coffee leaping out of its cup, SE/ℏS_E/\hbar is astronomically large, so e−2SE/ℏe^{-2S_E/\hbar} is beyond-astronomically small. Same exponential, different address on it: "classical" objects tunnel too, just at rates the universe's lifetime cannot notice.

Path integral and canonical quantization give the same theory — so which one is fundamental?

Neither, in any operationally meaningful sense: their equivalence is one of the standard results of QFT, so they are two descriptions of one theory, each keeping a different piece of classical structure manifest. The path integral wears the action and its symmetries on its sleeve — Lorentz invariance is visible, Ward identities drop out of a change of variables, the classical limit is a saddle point, and the Euclidean form hands you lattice numerics. The canonical route keeps the bracket algebra and unitarity manifest — operators, states, particles, and conservation of probability are natural there, via [ϕ^,π^]=iℏ δ3[\hat\phi, \hat\pi] = i\hbar\,\delta^3. Practitioners switch freely depending on the question. In the article's phrase: two roads, one quantum field theory — and everything classical is cargo that makes the crossing on both.

Is ∫Dϕ\int\mathcal{D}\phi actually a well-defined integral?

Not naively, and pretending otherwise costs real physics. In Minkowski signature the weight eiS/ℏe^{iS/\hbar} is a pure phase, so the integral is oscillatory and not absolutely convergent; Dϕ\mathcal{D}\phi is not a product of Lebesgue measures. Two honest definitions exist. One is constructive: the time-slicing limit of the derivation, where the integral is a limit of NN ordinary integrals as N→∞N \to \infty. The other is the Wick rotation, after which e−SE/ℏe^{-S_E/\hbar} defines a genuine measure — but even then only after regularization (a lattice, a cutoff) and careful limits. The practical moral from the warning: formal manipulations like changes of variables can silently fail, which is exactly what an anomaly is. Yet the regulated definitions are solid enough that lattice QCD extracts hadron masses to sub-percent accuracy.

Before this

Where this leads