Skip to content

Commutators, Currents, and Ward Identities

Vacuum commutator expectations and current Ward identities describe different singularities of correlation functions: the former compare Lorentzian operator orderings, while the latter determine contact terms at transformed insertions. This lesson first uses an exact scalar conformal two-point function, then derives equations of motion inside a stable regulated Euclidean integral and obtains the current identities. The final conditional pole argument introduces spontaneous breaking of a continuous global symmetry.

The current records how the action changes when a global symmetry parameter is allowed to vary in space. In a correlation function this produces contact terms at the charged insertions. The localized parameter is an arbitrary admissible infinitesimal test function; it need not vary slowly.

The bridge between the two themes is locality. Commutators tell us where an operator can causally respond. Ward identities tell us how a local symmetry variation is distributed among operator insertions. Both are statements about singularities of correlation functions.

Required background. Lesson 17 supplies ordered vacuum boundary values, distributional light-cone limits and the distinction between a raw time-ordered correlator and causal response.

Consider a nonidentity Hermitian scalar primary in a unitary conformal vacuum, with admissible scaling dimension Δ\Delta and positive normalization COC_O. Its exact Euclidean two-point function is

⟨O(xE)O(0)⟩=CO(xE2)Δ.\langle O(x_E)O(0)\rangle={C_O\over (x_E^2)^\Delta}.

This is a statement at all noncoincident separations in the exact CFT. A short-distance asymptote in a more general theory cannot determine a finite-separation correlator. For example, in D>2D>2 the leading short-distance power of a massive free scalar is massless, but it cannot be its exact correlator: away from coincidence the massive Klein–Gordon operator acting on the massless power leaves m2W+(0)≠0m^2W_+^{(0)}\ne0.

After analytic continuation to Lorentzian time, the two Wightman orderings are distributional boundary values on opposite sides of the light-cone singularity:

W+(t,x)=CO(x2−(t−i0)2)Δ=CO(ρ+i0 sgn⁡t)Δ,W_+(t,\mathbf x) ={C_O\over \left(\mathbf x^2-(t-i0)^2\right)^\Delta} ={C_O\over \left(\rho+i0\,\operatorname{sgn}t\right)^\Delta},

and

W−(t,x)=CO(x2−(t+i0)2)Δ=CO(ρ−i0 sgn⁡t)Δ.W_-(t,\mathbf x) ={C_O\over \left(\mathbf x^2-(t+i0)^2\right)^\Delta} ={C_O\over \left(\rho-i0\,\operatorname{sgn}t\right)^\Delta}.

Therefore

⟨0∣[O(t,x),O(0)]∣0⟩=W+(t,x)−W−(t,x).\langle0|[O(t,\mathbf x),O(0)]|0\rangle =W_+(t,\mathbf x)-W_-(t,\mathbf x).

See Gillioz 2023, arXiv v3, § 3.4, pp. 27–29 (PDF). His ordering is WG(x)=⟨ϕ(0)ϕ(x)⟩W_G(x)=\langle\phi(0)\phi(x)\rangle, so the present W+(x)=WG(−x)W_+(x)=W_G(-x). With his squared interval equal to ρ\rho, this translation gives the displayed lower-time prescription and the discontinuity sign below.

For spacelike separation, ρ>0\rho>0, the two boundary values are the same. The commutator vanishes:

⟨0∣[O(t,x),O(0)]∣0⟩=0(x2>t2).\boxed{ \langle0|[O(t,\mathbf x),O(0)]|0\rangle=0 \qquad (\mathbf x^2>t^2). }

This is the two-point, vacuum-expectation version of causality. Microcausality is the stronger operator statement [O(x),O(0)]=0[O(x),O(0)]=0 at spacelike separation (or the graded version for fermionic fields).

The branch-cut and causal-support diagram in Lesson 17 shows how the two banks interchange when time changes sign, and why spacelike values agree.

For non-integer Δ\Delta, the timelike discontinuity is easy to compute. When ρ<0\rho<0, write ρ=−∣ρ∣\rho=-|\rho|. On the principal branch,

(ρ+i0)−Δ=∣ρ∣−Δe−iπΔ,(ρ−i0)−Δ=∣ρ∣−Δe+iπΔ.(\rho+i0)^{-\Delta}=|\rho|^{-\Delta}e^{-i\pi\Delta}, \qquad (\rho-i0)^{-\Delta}=|\rho|^{-\Delta}e^{+i\pi\Delta}.

Taking into account the factor sgn⁡t\operatorname{sgn}t in the boundary prescription gives

W+(t,x)−W−(t,x)=−2iCOsin⁡(πΔ) sgn⁡t θ(t2−x2)(t2−x2)Δ\boxed{ W_+(t,\mathbf x)-W_-(t,\mathbf x) =-2iC_O\sin(\pi\Delta)\,\operatorname{sgn}t\, {\theta(t^2-\mathbf x^2)\over (t^2-\mathbf x^2)^\Delta} }

for non-integer Δ\Delta, away from the light cone.

The negative sign follows from the displayed ordering and principal branch. Globally the result means the distribution inherited from the analytic Wightman boundary values, including the light cone; for Δ≥1\Delta\ge1, the power times the step function is not independently an ordinary locally integrable function there. Its support satisfies

supp⁡⟨[O(x),O(0)]⟩⊆{x:t2≥x2}.\operatorname{supp}\langle[O(x),O(0)]\rangle \subseteq \{x: t^2\ge \mathbf x^2\}.

For integer Δ\Delta, the factor sin⁡(πΔ)\sin(\pi\Delta) seems to vanish. This does not mean that the commutator vanishes. Instead the discontinuity collapses to derivatives of delta functions on the light cone. For example,

1ρ+i0−1ρ−i0=−2πi δ(ρ).{1\over \rho+i0}-{1\over \rho-i0}=-2\pi i\,\delta(\rho).

This is why the free massless scalar in four dimensions has a commutator supported on the light cone rather than throughout the timelike interior.

Equal-time algebra and the ultraviolet cutoff

Section titled “Equal-time algebra and the ultraviolet cutoff”

A continuum theory can describe the long-distance limit of lattice degrees of freedom. For lattice spacing aa, the ultraviolet cutoff is of order

ΛUV∼1a.\Lambda_{\mathrm{UV}}\sim {1\over a}.

A continuum scaling law must be restricted to its scaling regime, with separations much larger than aa and the correlation length controlled when it is finite. Canonical equal-time commutators are local. For a scalar with the standard kinetic normalization, their singular support lies at coincident spatial points:

[ϕ(t,x),ϕ(t,y)]=0,[\phi(t,\mathbf x),\phi(t,\mathbf y)]=0,

and

[ϕ(t,x),π(t,y)]=iδ(d)(x−y),π=ϕ˙.[\phi(t,\mathbf x),\pi(t,\mathbf y)]=i\delta^{(d)}(\mathbf x-\mathbf y), \qquad \pi=\dot\phi.

The delta function is the continuum limit of a sharply localized lattice kernel. It is not resolved by the long-distance CFT power law. This is the same warning in a different language: the normalization and possible local improvement of contact terms are ultraviolet data that survive as distributions.

The mechanical analogy is useful but must be handled carefully. For a single nonlinear oscillator,

x¨+ω2x+λx3=0,\ddot x+\omega^2x+\lambda x^3=0,

the initial data x(0),x˙(0)x(0),\dot x(0) determine x(t)x(t). For a field, the equation of motion propagates operator data from one time slice to another. But products of fields at the same point are singular, so the equation of motion is an operator identity only after specifying renormalized composite operators and contact terms.

A stable regulated Schwinger–Dyson identity

Section titled “A stable regulated Schwinger–Dyson identity”

An explicit finite regulator makes the integration-by-parts premise precise. Take real field variables on a finite periodic Euclidean lattice and the normalized weight e−Sa/Zae^{-S_a}/Z_a, where

Sa=aD∑n[12∑μ=1D(ϕn+μ^−ϕna)2+m022ϕn2+λ04!ϕn4],λ0>0,m02≥0.\begin{aligned} S_a&=a^D\sum_n\left[ \frac12\sum_{\mu=1}^D\left(\frac{\phi_{n+\hat\mu}-\phi_n}{a}\right)^2 +\frac{m_0^2}{2}\phi_n^2+\frac{\lambda_0}{4!}\phi_n^4\right],\\ &\hspace{1em}\lambda_0>0,\qquad m_0^2\ge0. \end{aligned}

There are finitely many ordinary real integrals. The positive quartic term dominates every polynomial insertion at large field values, so the integral of ∂ϕn(e−SaO)\partial_{\phi_n}(e^{-S_a}\mathcal O) has zero field-space boundary flux. Define

Ea,n=a−D∂Sa∂ϕn=−Δaϕn+m02ϕn+λ06ϕn3,Δaϕn=1a2∑μ(ϕn+μ^+ϕn−μ^−2ϕn).\begin{aligned} E_{a,n}&=a^{-D}\frac{\partial S_a}{\partial\phi_n} =-\Delta_a\phi_n+m_0^2\phi_n+\frac{\lambda_0}{6}\phi_n^3,\\ \Delta_a\phi_n&=\frac1{a^2}\sum_\mu (\phi_{n+\hat\mu}+\phi_{n-\hat\mu}-2\phi_n). \end{aligned}

Ordinary integration by parts now proves the exact finite-dimensional identity

⟨Ea,nO⟩a=a−D⟨∂O∂ϕn⟩a.\langle E_{a,n}\mathcal O\rangle_a =a^{-D}\left\langle\frac{\partial\mathcal O}{\partial\phi_n}\right\rangle_a.

For a product of lattice fields the contacts are a−Dδn,nka^{-D}\delta_{n,n_k}, with the corresponding insertion omitted. On a bounded field-integration domain, any nonzero boundary flux must instead be retained. A real cubic action over the full real contour does not meet this positive-integral premise: its constant-field direction is unbounded below.

The continuum Schwinger–Dyson notation summarizes the corresponding identity after a continuum/composite prescription has been supplied. For a scalar field and a product of insertions O[ϕ]\mathcal O[\phi],

0=∫Dϕ δδϕ(x)(e−S[ϕ]O[ϕ]).0=\int D\phi\,{\delta\over\delta\phi(x)} \left(e^{-S[\phi]}\mathcal O[\phi]\right).

Therefore

⟨δSδϕ(x)O⟩=⟨δOδϕ(x)⟩.\boxed{ \left\langle {\delta S\over\delta\phi(x)}\mathcal O\right\rangle =\left\langle {\delta\mathcal O\over\delta\phi(x)}\right\rangle. }

For

O=ϕ(x1)⋯ϕ(xn),\mathcal O=\phi(x_1)\cdots\phi(x_n),

this becomes

⟨δSδϕ(x)ϕ(x1)⋯ϕ(xn)⟩=∑k=1nδ(D)(x−xk)⟨ϕ(x1)⋯ϕ(xk)^⋯ϕ(xn)⟩,\left\langle {\delta S\over\delta\phi(x)}\phi(x_1)\cdots\phi(x_n)\right\rangle =\sum_{k=1}^n\delta^{(D)}(x-x_k) \left\langle\phi(x_1)\cdots\widehat{\phi(x_k)}\cdots\phi(x_n)\right\rangle,

where the hat means omission. Away from x=xkx=x_k, the classical-looking equation of motion is valid inside correlation functions. At x=xkx=x_k, the contact terms are part of the identity.

For the quartic example, the formal continuum expressions are

SE=∫dDx[12(∂ϕ)2+m022ϕ2+λ04!ϕ4],S_E=\int d^D x\left[\frac12(\partial\phi)^2+\frac{m_0^2}{2}\phi^2+\frac{\lambda_0}{4!}\phi^4\right],

then

δSEδϕ(x)=−∂2ϕ(x)+m02ϕ(x)+λ06ϕ3(x),\frac{\delta S_E}{\delta\phi(x)}=-\partial^2\phi(x)+m_0^2\phi(x)+\frac{\lambda_0}{6}\phi^3(x),

so the equation

∂2ϕ=m02ϕ+λ06ϕ3\partial^2\phi=m_0^2\phi+\frac{\lambda_0}{6}\phi^3

is the bare regulated equation, not a pointwise equality of unrenormalized continuum products. Regulator removal requires counterterms and defined composite insertions; the renormalized equation and its contacts must be obtained together. The finite integral does not by itself prove a continuum limit. Compare the field-shift derivation in Schwartz 2014, § 14.7.1, pp. 273–274, and § 14.7.2, p. 275: his Lorentzian weight gives different factors of ii, while the Euclidean sign here follows directly from differentiating e−SEe^{-S_E}.

Noether currents from localizing a symmetry

Section titled “Noether currents from localizing a symmetry”

Suppose the action is invariant under a constant internal transformation

δΦk=αRkΦk,α=constant.\delta\Phi_k=\alpha R_k\Phi_k, \qquad \alpha=\text{constant}.

For an action depending on fields and their first derivatives, promote α\alpha to a smooth compactly supported function α(x)\alpha(x). To first order in its amplitude, the variation is

δS=∫dDx Jμ(x)∂μα(x).\delta S=\int d^D x\,J_\mu(x)\partial_\mu\alpha(x).

Integrating by parts gives

δS=−∫dDx α(x)∂μJμ(x),\delta S=-\int d^D x\,\alpha(x)\partial_\mu J_\mu(x),

with no boundary term for this parameter. Since constant α\alpha is an exact symmetry, only derivatives of α\alpha appear. Their coefficient is the Noether current. Higher-derivative actions can first produce higher derivatives of the parameter; integration by parts then identifies a current, with its usual improvement freedom.

For a complex scalar with Euclidean action

S=∫dDx(∂μϕ∗∂μϕ+V(∣ϕ∣2)),S=\int d^D x\left(\partial_\mu\phi^*\partial_\mu\phi+V(|\phi|^2)\right),

the global U(1)U(1) transformation is

ϕ↦eiαϕ,ϕ∗↦e−iαϕ∗.\phi\mapsto e^{i\alpha}\phi, \qquad \phi^*\mapsto e^{-i\alpha}\phi^*.

For local α(x)\alpha(x),

δϕ=iαϕ,δϕ∗=−iαϕ∗.\delta\phi=i\alpha\phi, \qquad \delta\phi^*=-i\alpha\phi^*.

The variation of the kinetic term is

δS=∫dDx i(ϕ∂μϕ∗−ϕ∗∂μϕ)∂μα.\delta S =\int d^D x\,i\left(\phi\partial_\mu\phi^*-\phi^*\partial_\mu\phi\right)\partial_\mu\alpha.

Thus one convenient Euclidean convention is

Jμ=i(ϕ∂μϕ∗−ϕ∗∂μϕ).J_\mu=i\left(\phi\partial_\mu\phi^*-\phi^*\partial_\mu\phi\right).

The conservation equation

∂μJμ=0\partial_\mu J_\mu=0

holds as an operator statement away from insertions and, equivalently, as a Ward identity inside correlators.

Let

O=Φ1(x1)Φ2(x2)⋯Φn(xn).\mathcal O=\Phi_1(x_1)\Phi_2(x_2)\cdots\Phi_n(x_n).

Use a symmetry-preserving regulator and measure, an admissible compactly supported parameter, and no anomalous Jacobian. Perform the local change of integration variables. The regulated integral does not change:

0=∫DΦ δ(e−SO)=⟨δO⟩−⟨O δS⟩.0=\int D\Phi\,\delta\left(e^{-S}\mathcal O\right) =\left\langle\delta\mathcal O\right\rangle- \left\langle\mathcal O\,\delta S\right\rangle.

Using

δS=−∫dDx α(x)∂μJμ(x),\delta S=-\int d^D x\,\alpha(x)\partial_\mu J_\mu(x),

and

δO=∑k=1nα(xk)Φ1(x1)⋯RkΦk(xk)⋯Φn(xn),\delta\mathcal O=\sum_{k=1}^n\alpha(x_k) \Phi_1(x_1)\cdots R_k\Phi_k(x_k)\cdots\Phi_n(x_n),

we obtain

0=∫dDx α(x)[∂μ⟨Jμ(x)O⟩+∑k=1nδ(D)(x−xk)⟨Φ1⋯RkΦk⋯Φn⟩].0=\int d^D x\,\alpha(x) \left[ \partial_\mu\langle J_\mu(x)\mathcal O\rangle +\sum_{k=1}^n\delta^{(D)}(x-x_k) \langle\Phi_1\cdots R_k\Phi_k\cdots\Phi_n\rangle \right].

Because α(x)\alpha(x) is arbitrary,

∂μ⟨Jμ(x)Φ1(x1)⋯Φn(xn)⟩=−∑k=1nδ(D)(x−xk)⟨Φ1⋯RkΦk⋯Φn⟩.\boxed{ \partial_\mu\langle J_\mu(x)\Phi_1(x_1)\cdots\Phi_n(x_n)\rangle =-\sum_{k=1}^n\delta^{(D)}(x-x_k) \langle\Phi_1\cdots R_k\Phi_k\cdots\Phi_n\rangle. }

For U(1)U(1) charges eke_k, with Rk=iekR_k=i e_k, this is

∂μ⟨Jμ(x)Φ1(x1)⋯Φn(xn)⟩=−i∑k=1nekδ(D)(x−xk)⟨Φ1(x1)⋯Φn(xn)⟩.\boxed{ \partial_\mu\langle J_\mu(x)\Phi_1(x_1)\cdots\Phi_n(x_n)\rangle =-i\sum_{k=1}^n e_k\delta^{(D)}(x-x_k) \langle\Phi_1(x_1)\cdots\Phi_n(x_n)\rangle. }

The right-hand side says that the current insertion detects the charge carried by the local operators. On a lattice, the exact identity uses a link current and a lattice divergence. The continuum expression presupposes the corresponding limit and contact prescription. A noninvariant measure adds its anomaly contribution.

This derivation follows Di Francesco, Mathieu and Sénéchal 1997, § 2.4.2, pp. 39–41, and § 2.4.4, pp. 43–44. Their definitions are δΦ=−iαGΦ\delta\Phi=-i\alpha G\Phi and δS=−∫jμ∂μα\delta S=-\int j_\mu\partial_\mu\alpha; the translation is R=−iGR=-iG and Jμ=−jμJ_\mu=-j_\mu, which gives precisely the contact sign above. For the charged scalar and Dμ=∂μ−ieAμD_\mu=\partial_\mu-ieA_\mu, the current defined by localizing the transformation is −δSE/δAμ-\delta S_E/\delta A_\mu at zero background. This source dictionary must accompany a change to background-field conventions.

The figure locates the contact terms at the charged insertions. Integrating over a region collects only the transformations of the operators it contains.

Region A encloses neither insertion and has zero integrated Ward contact; region B encloses only the first insertion and selects its symmetry variation

For O=Φ1(x1)Φ2(x2)\mathcal O=\Phi_1(x_1)\Phi_2(x_2), define IR=∫∂RdΣμ⟨JμO⟩I_R=\int_{\partial R}d\Sigma_\mu\langle J_\mu\mathcal O\rangle. Region AA contains no insertion, so IA=0I_A=0; region BB contains only x1x_1, so IB=−⟨(R1Φ1)Φ2⟩I_B=-\langle(R_1\Phi_1)\Phi_2\rangle. The boundaries avoid the insertions, and the stated Ward hypotheses apply. Arrows show outward surface normals, not a current-flow pattern. The Euclidean section is schematic.

As a useful check, integrate the Ward identity over a region RR. Gauss’s theorem gives

∫RdDx ∂μ⟨Jμ(x)O⟩=∫∂RdΣμ ⟨Jμ(x)O⟩.\int_R d^D x\,\partial_\mu\langle J_\mu(x)\mathcal O\rangle =\int_{\partial R} d\Sigma_\mu\,\langle J_\mu(x)\mathcal O\rangle.

If the boundary term vanishes and RR contains all insertions, then

∑k=1n⟨Φ1⋯RkΦk⋯Φn⟩=0.\sum_{k=1}^n\langle\Phi_1\cdots R_k\Phi_k\cdots\Phi_n\rangle=0.

For a U(1)U(1) correlator of charged fields, this becomes the charge-selection rule

(∑k=1nek)⟨Φ1(x1)⋯Φn(xn)⟩=0.\left(\sum_{k=1}^n e_k\right) \langle\Phi_1(x_1)\cdots\Phi_n(x_n)\rangle=0.

Thus a correlator in an invariant vacuum with no boundary flux can be nonzero only if the total charge is zero. A symmetry-breaking state or symmetry-breaking boundary condition removes this selection-rule premise; it need not be a state of nonzero total charge.

Ward identities are especially powerful in momentum space. With Euclidean momenta, define the full Fourier distribution

Gμ(q;p1,…,pn)=∫dDx∏kdDxk e+iq⋅x+i∑kpk⋅xk⟨Jμ(x)Φ1(x1)⋯Φn(xn)⟩.G_\mu(q;p_1,\ldots,p_n) =\int d^D x\prod_k d^D x_k\, e^{+iq\cdot x+i\sum_k p_k\cdot x_k} \langle J_\mu(x)\Phi_1(x_1)\cdots\Phi_n(x_n)\rangle.

This is the site-wide forward Fourier convention, so integration by parts gives ∂μ↦−iqμ\partial_\mu\mapsto-iq_\mu. Fourier transforming the position-space Ward identity therefore gives

qμGμ(q;p1,…,pn)=+∑k=1nek G(p1,…,pk+q,…,pn)\boxed{ q_\mu G_\mu(q;p_1,\ldots,p_n) =+\sum_{k=1}^n e_k\, G(p_1,\ldots,p_k+q,\ldots,p_n) }

for U(1)U(1) fields, after the common factor of ii has been canceled. A source using the opposite phase is translated by reversing every momentum, Gsite(q,pk)=G−(−q,−pk)G_{\rm site}(q,p_k)=G_{-}( -q,-p_k); it then returns to the displayed formula. Changing the generator convention also changes the displayed charge sign. The invariant content is that contracting the current momentum is equivalent to inserting the symmetry generator on each charged leg.

The momentum-space figure follows one term of this identity: a position-space contact sets the current position equal to an insertion position, so their Fourier momenta add on that leg.

Contracting the current momentum shifts the charged leg from minus p minus q to minus p or its conjugate leg from p to p plus q, producing e times the difference of two propagators

The two-point example uses incoming labels −p−q-p-q and pp on fields of charges +e+e and −e-e. A contact shifts exactly one label by qq; the underline marks that change. After stripping the common translation delta, the contraction is qμGμ=e[G(p)−G(p+q)]q_\mu G_\mu=e[G(p)-G(p+q)]. Arrows denote Fourier momentum labels of correlators, not physical charge flow. The diagram is schematic; the amputation below fixes the proper-vertex identity.

For a charged scalar and its conjugate, choose the incoming Fourier labels p1=−p−qp_1=-p-q on ϕ\phi and p2=pp_2=p on ϕ∗\phi^*. Strip the common translation delta (2π)Dδ(D)(q+p1+p2)(2\pi)^D\delta^{(D)}(q+p_1+p_2) from the three-point function and write the resulting current kernel as Gμ(p+q,p)G_\mu(p+q,p). For the rotationally invariant two-point kernel G(p)G(p), the two opposite charges then give

qμGμ(p+q,p)=e [G(p)−G(p+q)].q_\mu G_\mu(p+q,p)=e\,[G(p)-G(p+q)].

Where the full scalar propagators are invertible, define the amputated current vertex by

Gμ(p+q,p)=G(p+q)Γμ(p+q,p)G(p).G_\mu(p+q,p)=G(p+q)\Gamma_\mu(p+q,p)G(p).

Then the Ward identity becomes

qμΓμ(p+q,p)=e[G−1(p+q)−G−1(p)].\boxed{ q_\mu\Gamma_\mu(p+q,p) =e\left[G^{-1}(p+q)-G^{-1}(p)\right]. }

This formula is one of the cleanest expressions of symmetry in perturbative QFT. It says that the longitudinal part of the current vertex is not independent dynamical data; it is fixed by the inverse propagator.

If the vertex has a direction-independent regular limit at q=0q=0 and G−1G^{-1} is differentiable, expansion for arbitrary small qq gives

Γμ(p,p)=e ∂G−1(p)∂pμ,\Gamma_\mu(p,p)=e\,{\partial G^{-1}(p)\over\partial p_\mu},

in the Euclidean convention just defined. For a free scalar, G−1(p)=pE2+m2G^{-1}(p)=p_E^2+m^2 and Γμ=e(2p+q)μ\Gamma_\mu=e(2p+q)_\mu satisfy the finite-qq identity directly. A singular or direction-dependent soft vertex does not justify the limiting formula. In QED the corresponding matched Ward identity relates the vertex and fermion wavefunction renormalization factors; the renormalized charge also involves the photon normalization.

Ward identities also know about spontaneous symmetry breaking. Suppose a continuous symmetry acts nontrivially on an operator Φ\Phi, and the vacuum expectation value is nonzero:

⟨Φ⟩=v≠0.\langle\Phi\rangle=v\ne0.

Apply the Ward identity to the correlator ⟨Jμ(x)Φ(0)⟩\langle J_\mu(x)\Phi(0)\rangle. If

δΦ=ieαΦ,\delta\Phi=i e\alpha\Phi,

then

∂μ⟨Jμ(x)Φ(0)⟩=−iev δ(D)(x).\partial_\mu\langle J_\mu(x)\Phi(0)\rangle =-i e v\,\delta^{(D)}(x).

Fourier transforming gives

qμCμ(q)=ev,Cμ(q)=∫dDx e+iq⋅x⟨Jμ(x)Φ(0)⟩.q_\mu C_\mu(q)=ev, \qquad C_\mu(q)=\int d^D x\,e^{+iq\cdot x}\langle J_\mu(x)\Phi(0)\rangle.

Assume ev≠0ev\ne0. If Cμ(q)C_\mu(q) were regular at q=0q=0, the left side would vanish as q→0q\to0, contradicting this identity. Its longitudinal projection is therefore

Cμ∥(q)=ev qμqE2,qE2=∑μqμ2>0.C_\mu^{\parallel}(q)=ev\,\frac{q_\mu}{q_E^2}, \qquad q_E^2=\sum_\mu q_\mu^2>0.

Regular transverse terms may also be present, but they cannot saturate the nonzero divergence.

Under the relativistic vacuum and positive spectral assumptions stated below, the massless longitudinal contribution has the Goldstone interpretation. The spectral step combines Lorentz covariance, locality, current conservation and a nonzero symmetry variation; see Weinberg 1995, Vol. II, § 19.2, pp. 169–172. A Euclidean projection formula alone is not a proof of all the hypotheses of a Lorentzian particle theorem.

The figure compares the regularity assumption with the nonzero divergence. Its upper plot uses the scalar longitudinal coefficient FL=ev/qE2F_L=ev/q_E^2: for a reference momentum scale μ>0\mu>0, set s=qE2/μ2s=q_E^2/\mu^2 and f=μ2FL/(ev)f=\mu^2F_L/(ev). The pole gives f=1/sf=1/s, while a regular control has f=1f=1. The lower plot displays the corresponding contractions sf=1sf=1 and sf=s→0sf=s\to0. Only the pole maintains the required nonzero limit.

A nonzero Ward divergence rules out a regular current correlator at zero momentum; its normalized Euclidean longitudinal coefficient grows as one over positive squared momentum

For ev≠0ev\ne0, the Ward identity fixes Cμ∥=evqμ/qE2C_\mu^{\parallel}=evq_\mu/q_E^2. The upper curves show exact dimensionless f=1/sf=1/s and the regular control f=1f=1 on 0.1≤s≤10.1\le s\le1; the lower curves compare their contractions as s↓0s\downarrow0. The open circle marks the pole’s limiting contraction, not a value of ff at zero. Solid and dashed styles distinguish the pole and control. These are normalized Euclidean functions, not a Lorentzian two-sided pole prescription. The next lesson supplies the further spectral interpretation.

Once the current Ward identity is established, this singularity argument does not require weak coupling or a classical potential. Its use of a noninvariant vacuum is consequential: the infinite-volume limit and state selection precede removal of an infinitesimal symmetry-breaking source. A finite symmetric integral does not itself select that state. The action-based derivation above is one way of obtaining the identity; an operator derivation can start from the same local current relation.

The Lorentzian commutator expectation of two local fields is the difference between two Wightman boundary values. For a conformal power law, this difference is a branch-cut discontinuity. It vanishes in the spacelike region because there the two boundary values coincide. On the light cone, and sometimes inside it, this discontinuity supplies the spectral function from which retarded and advanced response are built.

Current conservation has an equally local interpretation. A current is obtained by promoting a global symmetry parameter to α(x)\alpha(x) and reading off the coefficient of ∂μα\partial_\mu\alpha. Inside correlation functions, the divergence of this current vanishes away from operator insertions and produces contact terms at the insertions. These contact terms are Ward identities.

In momentum space, a Ward identity says that contracting a current insertion with its momentum is equivalent to acting with the symmetry generator on each external operator. If the chosen vacuum is not invariant, the same identity forces a massless pole: the first glimpse of Goldstone physics.

Confusing time ordering with response. The raw Feynman correlator DFD_F is time ordered; the commutator expectation is the difference of two Wightman orderings. The former can be nonzero at spacelike separation. For a Hamiltonian source Hext=−∫JOH_{\rm ext}=-\int JO, the causal response kernel is R=iθ(t)⟨[O(x),O(0)]⟩\mathcal R=i\theta(t)\langle[O(x),O(0)]\rangle, while the site’s GR=−iθ(t)⟨[O(x),O(0)]⟩=−RG_R=-i\theta(t)\langle[O(x),O(0)]\rangle=-\mathcal R. These are the definitions in Lesson 17.

Dropping light-cone distributions. The off-cone formula with sin⁡(πΔ)\sin(\pi\Delta) can vanish at integer dimension while delta functions and their derivatives remain on the light cone. For noninteger dimension its global meaning also comes from the original boundary distribution.

Omitting contact terms. The equation ∂μJμ=0\partial_\mu J_\mu=0 holds away from insertions. Charged insertions contribute the specified delta-function terms.

Ignoring the regulator or Jacobian. A Ward identity follows from an admissible change of variables, with the measure and boundary conditions accounted for. A noninvariant measure adds an anomaly; a nonzero boundary flux cannot simply be discarded.

Let

W+(t,x)=C(ρ+i0sgn⁡t)Δ,W−(t,x)=C(ρ−i0sgn⁡t)Δ,ρ=x2−t2.W_+(t,\mathbf x)={C\over(\rho+i0\operatorname{sgn}t)^\Delta}, \qquad W_-(t,\mathbf x)={C\over(\rho-i0\operatorname{sgn}t)^\Delta}, \qquad \rho=\mathbf x^2-t^2.

For non-integer Δ\Delta, compute W+−W−W_+-W_- for ρ>0\rho>0 and for ρ<0\rho<0.

Solution

For ρ>0\rho>0, the two boundary values approach the positive real axis from opposite sides. Since ρ−Δ\rho^{-\Delta} has no branch cut there,

(ρ+i0)−Δ=(ρ−i0)−Δ=ρ−Δ.(\rho+i0)^{-\Delta}=(\rho-i0)^{-\Delta}=\rho^{-\Delta}.

Therefore

W+−W−=0(ρ>0).W_+-W_-=0 \qquad (\rho>0).

For ρ<0\rho<0, write ρ=−∣ρ∣\rho=-|\rho|. If t>0t>0, then

ρ+i0=∣ρ∣eiπ,ρ−i0=∣ρ∣e−iπ,\rho+i0=|\rho|e^{i\pi}, \qquad \rho-i0=|\rho|e^{-i\pi},

so

W+−W−=C∣ρ∣−Δ(e−iπΔ−eiπΔ)=−2iCsin⁡(πΔ)∣ρ∣−Δ.W_+-W_- =C|\rho|^{-\Delta}\left(e^{-i\pi\Delta}-e^{i\pi\Delta}\right) =-2iC\sin(\pi\Delta)|\rho|^{-\Delta}.

If t<0t<0, the boundary prescriptions are reversed, so the result changes sign. Thus

W+−W−=−2iCsin⁡(πΔ)sgn⁡t θ(−ρ)(−ρ)Δ.W_+-W_- =-2iC\sin(\pi\Delta)\operatorname{sgn}t\,{\theta(-\rho)\over(-\rho)^\Delta}.

Since −ρ=t2−x2-\rho=t^2-\mathbf x^2, this is

W+−W−=−2iCsin⁡(πΔ)sgn⁡t θ(t2−x2)(t2−x2)Δ.W_+-W_- =-2iC\sin(\pi\Delta)\operatorname{sgn}t\, {\theta(t^2-\mathbf x^2)\over(t^2-\mathbf x^2)^\Delta}.

Exercise 2: Integer dimension and light-cone support

Section titled “Exercise 2: Integer dimension and light-cone support”

Use the distribution identity

1x+i0−1x−i0=−2πiδ(x){1\over x+i0}-{1\over x-i0}=-2\pi i\delta(x)

to compute the commutator discontinuity for Δ=1\Delta=1.

Solution

For Δ=1\Delta=1,

W+−W−=C[1ρ+i0sgn⁡t−1ρ−i0sgn⁡t].W_+-W_-=C\left[{1\over\rho+i0\operatorname{sgn}t}-{1\over\rho-i0\operatorname{sgn}t}\right].

If t>0t>0, this is

C(1ρ+i0−1ρ−i0)=−2πiCδ(ρ).C\left({1\over\rho+i0}-{1\over\rho-i0}\right) =-2\pi iC\delta(\rho).

If t<0t<0, the i0i0 prescriptions are reversed, giving the opposite sign. Hence

W+−W−=−2πiCsgn⁡t δ(ρ).W_+-W_-=-2\pi iC\operatorname{sgn}t\,\delta(\rho).

Since ρ=x2−t2\rho=\mathbf x^2-t^2, this distribution is supported on the light cone. This is the light-cone version of the massless free-field commutator.

Exercise 3: Deriving the current Ward identity

Section titled “Exercise 3: Deriving the current Ward identity”

Derive the Ward identity

∂μ⟨Jμ(x)Φ1(x1)⋯Φn(xn)⟩=−∑k=1nδ(D)(x−xk)⟨Φ1⋯RkΦk⋯Φn⟩\partial_\mu\langle J_\mu(x)\Phi_1(x_1)\cdots\Phi_n(x_n)\rangle =-\sum_{k=1}^n\delta^{(D)}(x-x_k) \langle\Phi_1\cdots R_k\Phi_k\cdots\Phi_n\rangle

from a local change of variables in the Euclidean path integral, assuming the invariant regulator/measure and compactly supported parameter used in the text.

Solution

Let

O=Φ1(x1)⋯Φn(xn).\mathcal O=\Phi_1(x_1)\cdots\Phi_n(x_n).

Perform the infinitesimal change of variables

δΦk=α(x)RkΦk.\delta\Phi_k=\alpha(x)R_k\Phi_k.

Assume the measure is invariant. Then

0=∫DΦ δ(e−SO)=⟨δO⟩−⟨OδS⟩.0=\int D\Phi\,\delta(e^{-S}\mathcal O) =\langle\delta\mathcal O\rangle-\langle\mathcal O\delta S\rangle.

The localized symmetry variation of the action is

δS=∫dDx Jμ∂μα=−∫dDx α∂μJμ.\delta S=\int d^D x\,J_\mu\partial_\mu\alpha =-\int d^D x\,\alpha\partial_\mu J_\mu.

The insertion varies as

δO=∑k=1nα(xk)Φ1⋯RkΦk⋯Φn.\delta\mathcal O= \sum_{k=1}^n\alpha(x_k)\Phi_1\cdots R_k\Phi_k\cdots\Phi_n.

Write this as an integral over xx:

δO=∫dDx α(x)∑k=1nδ(D)(x−xk)Φ1⋯RkΦk⋯Φn.\delta\mathcal O= \int d^D x\,\alpha(x) \sum_{k=1}^n\delta^{(D)}(x-x_k) \Phi_1\cdots R_k\Phi_k\cdots\Phi_n.

Substituting into 0=⟨δO⟩−⟨OδS⟩0=\langle\delta\mathcal O\rangle-\langle\mathcal O\delta S\rangle gives

0=∫dDx α(x)[∑k=1nδ(D)(x−xk)⟨Φ1⋯RkΦk⋯Φn⟩+∂μ⟨Jμ(x)O⟩].0=\int d^D x\,\alpha(x)\left[ \sum_{k=1}^n\delta^{(D)}(x-x_k) \langle\Phi_1\cdots R_k\Phi_k\cdots\Phi_n\rangle +\partial_\mu\langle J_\mu(x)\mathcal O\rangle \right].

Since α(x)\alpha(x) is arbitrary, the integrand must vanish. This gives the Ward identity.

Let ϕ\phi have U(1)U(1) charge +1+1 and ϕ∗\phi^* charge −1-1. Use the Ward identity to show that

⟨ϕ(x1)ϕ(x2)⟩=0\langle\phi(x_1)\phi(x_2)\rangle=0

in a symmetry-invariant vacuum, assuming the integrated current has no boundary flux.

Solution

For the two-point function

O=ϕ(x1)ϕ(x2),\mathcal O=\phi(x_1)\phi(x_2),

both insertions have charge +1+1. The integrated Ward identity for a region containing both insertions and with no boundary flux gives

(e1+e2)⟨ϕ(x1)ϕ(x2)⟩=0.(e_1+e_2)\langle\phi(x_1)\phi(x_2)\rangle=0.

Here e1=e2=1e_1=e_2=1, so

2⟨ϕ(x1)ϕ(x2)⟩=0.2\langle\phi(x_1)\phi(x_2)\rangle=0.

Thus

⟨ϕ(x1)ϕ(x2)⟩=0.\langle\phi(x_1)\phi(x_2)\rangle=0.

By contrast, ⟨ϕ(x1)ϕ∗(x2)⟩\langle\phi(x_1)\phi^*(x_2)\rangle is not forbidden, because the total charge is 1−1=01-1=0.

Exercise 5: The zero-momentum Ward–Takahashi vertex

Section titled “Exercise 5: The zero-momentum Ward–Takahashi vertex”

Use the Euclidean momentum routing and stripped translation delta defined in the text. Assume that G−1G^{-1} is differentiable, the vertex has a regular direction-independent limit as q→0q\to0, and the current three-point function is written as

Gμ(p+q,p)=G(p+q)Γμ(p+q,p)G(p).G_\mu(p+q,p)=G(p+q)\Gamma_\mu(p+q,p)G(p).

Starting from the Ward–Takahashi identity

qμΓμ(p+q,p)=e[G−1(p+q)−G−1(p)],q_\mu\Gamma_\mu(p+q,p)=e\left[G^{-1}(p+q)-G^{-1}(p)\right],

show that the zero-momentum current vertex is

Γμ(p,p)=e∂G−1(p)∂pμ.\Gamma_\mu(p,p)=e{\partial G^{-1}(p)\over\partial p_\mu}.
Solution

Expand the inverse propagator for small qq:

G−1(p+q)=G−1(p)+qν∂G−1(p)∂pν+o(∣q∣).G^{-1}(p+q)=G^{-1}(p)+q_\nu{\partial G^{-1}(p)\over\partial p_\nu}+o(|q|).

Then the Ward–Takahashi identity becomes

qμΓμ(p+q,p)=eqν∂G−1(p)∂pν+o(∣q∣).q_\mu\Gamma_\mu(p+q,p) =e q_\nu{\partial G^{-1}(p)\over\partial p_\nu}+o(|q|).

Assuming the vertex is regular as q→0q\to0, set

Γμ(p+q,p)=Γμ(p,p)+o(1).\Gamma_\mu(p+q,p)=\Gamma_\mu(p,p)+o(1).

The leading terms give

qμΓμ(p,p)=eqμ∂G−1(p)∂pμ.q_\mu\Gamma_\mu(p,p) =e q_\mu{\partial G^{-1}(p)\over\partial p_\mu}.

Since this holds for arbitrary infinitesimal qμq_\mu,

Γμ(p,p)=e∂G−1(p)∂pμ.\Gamma_\mu(p,p)=e{\partial G^{-1}(p)\over\partial p_\mu}.

Suppose a continuous symmetry current satisfies

∂μ⟨Jμ(x)Φ(0)⟩=−vδ(D)(x),v≠0.\partial_\mu\langle J_\mu(x)\Phi(0)\rangle=-v\delta^{(D)}(x), \qquad v\ne0.

Show that the Fourier transform Cμ(q)C_\mu(q) cannot be regular at q=0q=0.

Solution

Define

Cμ(q)=∫dDx e+iq⋅x⟨Jμ(x)Φ(0)⟩.C_\mu(q)=\int d^D x\,e^{+iq\cdot x}\langle J_\mu(x)\Phi(0)\rangle.

The forward transform sends the derivative to −iqμ-iq_\mu. Fourier transforming the assumed identity gives

qμCμ(q)=−iv.q_\mu C_\mu(q)=-iv.

If Cμ(q)C_\mu(q) were regular at q=0q=0, then qμCμ(q)q_\mu C_\mu(q) would vanish as q→0q\to0. This contradicts the nonzero right-hand side. Therefore Cμ(q)C_\mu(q) must be singular.

Rotational invariance implies that the singular longitudinal part has the form

Cμ∥(q)=−iv qμqE2,qE2>0.C_\mu^{\parallel}(q)=-iv\,\frac{q_\mu}{q_E^2}, \qquad q_E^2>0.

Terms transverse to qμq_\mu cannot saturate the nonzero divergence. The factor −i-i follows from this exercise’s assumed contact −vδ-v\delta; the body’s different contact −ievδ-iev\delta instead gives evev. The Goldstone interpretation of the longitudinal singularity requires the vacuum and spectral hypotheses stated in the text.

  • Coleman, Sidney. “There Are No Goldstone Bosons in Two Dimensions.” Communications in Mathematical Physics 31, no. 4 (1973): 259–264. Journal record. Open PDF.
  • Di Francesco, Philippe, Pierre Mathieu, and David Sénéchal. Conformal Field Theory. Graduate Texts in Contemporary Physics. Springer, 1997. DOI.
  • Gillioz, Marc. Conformal Field Theory for Particle Physicists: From QFT Axioms to the Modern Conformal Bootstrap. SpringerBriefs in Physics. Springer, 2023. DOI. Open PDF, arXiv:2207.09474v3, 3 May 2023. Locators above refer to the preprint’s printed page labels.
  • Schwartz, Matthew D. Quantum Field Theory and the Standard Model. Cambridge University Press, 2014. DOI.
  • Weinberg, Steven. The Quantum Theory of Fields, Volume II: Modern Applications. Cambridge University Press, 1995. DOI. Locators use the printed page labels of the consulted 2012 printing.

Original QFT.org content:CC BY 4.0, unless an item supplies different terms. Third-party material retains its own terms.