Contact Terms and Renormalized Operator Products

Renormalizing each local insertion separately determines a product only away from coincident points. When $x\to y$, new divergent subgraphs can contain both marked vertices. Their counterterms are distributions supported on the coincidence diagonal, such as $\delta^{(d)}(x-y)[O_c](y)$ and its derivatives. These contact terms do not alter separated-point correlators, but they are indispensable in integrated observables, repeated source derivatives, Ward identities, and operator-product expansions.

This page classifies that diagonal freedom and works the simplest example: the product of two $\phi^2/2$ insertions is divergent even in a free four-dimensional scalar theory. It also fixes the chapter-wide matrix convention record used on the following mixing and evolution pages.

Required background. Renormalized Composite-Operator Insertions defines one insertion through source differentiation. Coincident Products and Contact Terms supplies the distributional distinction between a separated product and its extension to a diagonal.

Helpful background. Contact Terms, Equal-Time Commutators, and Schwinger Terms explains their symmetry role. Free-Field OPE Preview separates the ordinary short-distance expansion from its interacting renormalization.

Extensions on the coincidence diagonal

Suppose the time-ordered product

$$t_{ab}(x-y) = T\,[O_a](x)[O_b](y)$$

has already been renormalized for $x\ne y$. A renormalized product is an extension of this distribution to $x=y$. If two extensions $t_{ab}^{(1)}$ and $t_{ab}^{(2)}$ agree away from the diagonal, their difference has support only at $x=y$. Locality therefore gives

$$t_{ab}^{(1)}(x-y)-t_{ab}^{(2)}(x-y) = \sum_{c,\alpha} c_{ab}^{c,\alpha}\, \partial^\alpha\delta^{(d)}(x-y)\,[O_c](y).$$

The sum is finite at any declared power-counting order. Lorentz symmetry, internal quantum numbers, exchange symmetry, ghost number, and other exact identities filter the allowed $O_c$ and multi-index $\alpha$. If the operators have engineering dimensions $\Delta_a$ and $\Delta_b$, dimensional analysis requires

$$d+|\alpha|+\Delta_c \le \Delta_a+\Delta_b,$$

with any difference supplied by masses or other dimensionful parameters. In a mass-independent homogeneous sector, only terms with equal total dimension occur. A cutoff can expose power-divergent mixing into lower-dimensional operators; dimensional regularization can hide those powers without changing the symmetry classification.

This extension problem is not optional notation. A distribution such as $1/(x^2)^2$ is well defined for $x\ne0$ in four dimensions but is not locally integrable at the origin. Its Fourier transform has a logarithmic ultraviolet divergence. Choosing a subtraction is precisely choosing an extension, and changing its finite part adds a multiple of $\delta^{(4)}(x)$.

Nonlinear source counterterms

The source method organizes all diagonals at once. For sources $s_a(x)$ coupled to $[O_a](x)$, locality permits

$$\begin{aligned} S_{\rm ct}[s] =S_{\rm ct}[0] &+\int d^dx\,s_a A_{ab}O_b\\ &+\frac12\int d^dx \sum_{a,b,c,\alpha,\beta} k_{ab}^{c,\alpha\beta} (\partial^\alpha s_a)(\partial^\beta s_b)O_c +\mathcal O(s^3). \end{aligned}$$

One source derivative sees only the first line. Two derivatives of the source-quadratic terms produce derivatives of $\delta^{(d)}(x-y)$ multiplying local operators. Thus

$$\begin{aligned} \frac{\delta^2W_{\rm R}}{\delta s_a(x)\delta s_b(y)} ={}& \left\langle T\,[O_a](x)[O_b](y) \right\rangle_{\rm R,c}\\ &+\sum_{c,\alpha} K_{ab}^{c,\alpha}\, \partial^\alpha\delta^{(d)}(x-y) \left\langle[O_c](y)\right\rangle_{\rm R}. \end{aligned}$$

Integration by parts moves derivatives between the sources, the delta distribution, and $O_c$; a convention must be fixed before comparing coefficients. Commutativity of functional derivatives constrains the exchange $a,x\leftrightarrow b,y$. Symmetry identities impose further relations. Arbitrarily choosing each contact coefficient independently can therefore make the source functional nonintegrable or violate a Ward identity.

Collins exhibits the same structure graphically: separate renormalizations of two insertions remove subgraphs around either marked vertex, while an overall subgraph containing both vertices generates a new local counterterm Collins 1984/2023, §§ 6.2.2–6.3, pp. 145–149.

Free scalar example: a divergent product of finite insertions

Let

$$O(x)=\frac12\phi^2(x)$$

in a free massive Euclidean scalar theory in four dimensions. At $x\ne0$, Wick’s theorem gives the connected product

$$\left\langle O(x)O(0)\right\rangle_c = \frac12\Delta(x)^2.$$

In the massless short-distance limit,

$$\Delta(x)=\frac1{4\pi^2x^2}, \qquad \left\langle O(x)O(0)\right\rangle_c = \frac1{32\pi^4x^4}.$$

The radial integral near the origin behaves as

$$\int_0^\eta dr\,r^3\frac1{r^4} = \int_0^\eta\frac{dr}{r},$$

so the singularity is logarithmic. The one-insertion operator $O$ is perfectly finite in the free theory; the product is not.

In $d=4-2\epsilon$, the Fourier transform is half the scalar bubble from the preceding page:

$$\begin{aligned} \Pi_0(p^2) &= \frac12\mu^{2\epsilon} \int\frac{d^{4-2\epsilon}k}{(2\pi)^{4-2\epsilon}}\, \frac1{(k^2+m^2)[(k+p)^2+m^2]}\\ &= \frac1{32\pi^2} \left[ \frac1{\bar\epsilon} -\int_0^1dx\, \ln\frac{m^2+x(1-x)p^2}{\mu^2} \right] +\mathcal O(\epsilon). \end{aligned}$$

Minimal subtraction defines

$$\Pi_{\rm R}(p^2) = \Pi_0(p^2)-\frac1{32\pi^2\bar\epsilon} = -\frac1{32\pi^2} \int_0^1dx\, \ln\frac{m^2+x(1-x)p^2}{\mu^2}.$$

The subtracted constant in momentum space is a $\delta^{(4)}(x)$ counterterm in position space. A position-space extension that makes the same logarithmic scale dependence explicit is

$$\left[\frac1{x^4}\right]_{\rm R} = -\frac14\Box \left[ \frac{\ln(\mu^2x^2)}{x^2} \right].$$

For $x\ne0$, applying $\Box$ recovers $1/x^4$. As a distribution, its scale derivative is local:

$$\mu\frac{d}{d\mu} \left\langle O(x)O(0)\right\rangle_{\rm R,c} = \frac1{16\pi^2}\delta^{(4)}(x).$$

The same result follows by differentiating $\Pi_{\rm R}$ with respect to $\mu$. This agreement checks the normalization and sign of the contact term.

A second valid scheme can add a finite constant $c$:

$$\left\langle O(x)O(0)\right\rangle_{\rm R,c}' = \left\langle O(x)O(0)\right\rangle_{\rm R,c} +c\,\delta^{(4)}(x).$$

All correlators with $x\ne0$ agree. The integrated susceptibility, zero-momentum source derivative, and local Ward identities can change unless the same finite contact convention is translated everywhere. “It vanishes at separated points” is therefore not a reason to discard $c$.

Contact terms in Ward identities

For a current $j^\mu$ generating an infinitesimal transformation $\delta$, a time-ordered Ward identity has the schematic distributional form

$$\partial_\mu \left\langle T\,j^\mu(x)O_1(x_1)\cdots O_n(x_n) \right\rangle = \sum_{r=1}^n \delta^{(d)}(x-x_r) \left\langle T\,O_1\cdots(\delta O_r)\cdots O_n \right\rangle +\mathcal A(x;\{x_r\}),$$

where overall factors of $i$ depend on the real-time or Euclidean convention. The delta terms are required: differentiating the time ordering and localizing the transformation produces them. The term $\mathcal A$ is reserved for a possible anomalous breaking after all admissible local counterterms have been considered.

Three statements must be kept distinct:

1. a symmetry-required contact term implements the transformation of another insertion;
2. a finite contact redefinition changes the local representative while preserving separated-point data;
3. a nontrivial anomaly is a consistent local breaking that cannot be removed by an allowed redefinition.

Calling every delta-supported term an anomaly confuses these categories. Conversely, dropping contact terms can manufacture an apparent violation of the Ward identity. The complete source-dependent identity, including nonlinear source counterterms, is the appropriate check.

Contact terms and the operator-product expansion

At separated points, an operator-product expansion has the form

$$[O_a](x)[O_b](0) \sim \sum_c C_{ab}^{\;\;c}(x,\mu)[O_c](0), \qquad x\ne0.$$

Its coefficient functions encode nonlocal short-distance dependence. When the product is integrated through $x=0$, the $C_{ab}^{\;\;c}$ themselves require distributional extensions. Two extensions can differ by derivatives of delta functions, exactly as above. Contact terms therefore belong to the renormalized product and to integrated OPE applications, but they are not ordinary values of the coefficient functions at nonzero separation.

Zimmermann proved perturbative short-distance expansions using renormalized normal products and made their finite redefinition freedom systematic Zimmermann 1973, pp. 570–601. This page uses the source-functional language because it makes multiple insertions and Ward identities visible in the same object.

Mixing convention record

This semantic table is the chapter’s required convention record. Later pages specialize it; no matrix result is accepted unless the applicable rows are declared and checked.

Entry	Chapter declaration	Where it enters	Required check
Operator orientation	$O$ is a column	Every matrix equation	Restore explicit indices and verify left action.
Renormalization direction	$O_0=Z_OO$	Pole subtraction and finite basis maps	Re-expand insertion vertices and cancel every pole.
Anomalous dimension	$\gamma=Z_O^{-1}\mu\,dZ_O/d\mu$	Operator RG equation	Differentiate $O_0$ and recover $\mu\,dO/d\mu=-\gamma O$.
Coefficient pairing	$\mathcal L_{\rm eff}\supset C^{\mathsf T}O$	Matching and RG evolution	Verify $\mu\,d(C^{\mathsf T}O)/d\mu=0$.
Coefficient flow	$\mu\,dC/d\mu=\gamma^{\mathsf T}C$	Wilson-coefficient running	A nonsymmetric test matrix must expose a missing transpose.
Physical sector	Representatives and quotient are named	Observable matrix elements	Show independence of allowed redundant shifts.
Equation-of-motion sector	Included when off-shell closure requires it	Off-shell Green functions and field redefinitions	Verify disappearance only in the stated on-shell or quotient projection.
Total derivatives	Retained for nonforward kinematics	Momentum-transfer matrix elements	Check the insertion momentum before dropping them.
BRST-exact sector	Included in gauge-fixed closure when applicable	Gauge-variant off-shell renormalization	Project to cohomology only after the full matrix closes.
Evanescent sector	Declared in dimensionally continued tensor or spinor bases	Pole subtraction beyond four dimensions	Include its pole-times-vanishing finite feedback.
Identity and lower dimension	Included when quantum numbers and dimensions allow	Vacuum terms and power-divergent mixing	Test zero-leg functions and regulator scaling.
Contact sector	Recorded separately from one-insertion $Z_O$	Coincident products and source derivatives	Compare separated points with integrated identities.
Protection	Exact identity, anomaly assumption, normalization, and improvement freedom stated	Currents and stress tensors	Check the complete Ward identity, not one zero matrix entry.
Finite basis change	$O'=BO$, $C'=B^{-\mathsf T}C$	Scheme or basis translation	Perform an exact coefficient–operator round trip.

Common pitfalls

Multiplying two one-insertion factors. The product $Z_{O_a}^{-1}Z_{O_b}^{-1}O_{0a}(x)O_{0b}(y)$ removes subgraphs localized at either insertion. It does not remove a subgraph containing both vertices, whose counterterm lives on $x=y$.

Inferring contacts from separated data. No measurement restricted to $x\ne y$ fixes a finite delta term. A normalization of an integrated source derivative, a Ward identity, or another local convention is required.

Dropping derivatives of delta functions by integration by parts. Moving a derivative changes which source or operator it acts on. The operation is legitimate only after boundary conditions and the source convention are fixed.

Calling a contact term unphysical. Individual finite contact coefficients are scheme dependent, but consistent contact terms affect integrated identities and are necessary for scheme-independent complete predictions.

Exercises

1. Classify the contact terms in the product of two scalar operators of dimension two in four dimensions, assuming a mass-independent scheme and no other quantum numbers.

Solution

The product has dimension four. A term $\partial^\alpha\delta^{(4)}(x)O_c$ has dimension $4+|\alpha|+\Delta_c$. Equality requires $|\alpha|=0$ and $\Delta_c=0$. Thus only

$$c\,\delta^{(4)}(x)\mathbf1$$

is allowed. Masses or a power-divergent regulator can enlarge the lower-dimensional bookkeeping, but they must be declared explicitly.

2. Verify the scale derivative of the differential extension of $1/x^4$.

Solution

Differentiate:

$$\mu\frac{d}{d\mu} \left[-\frac14\Box\frac{\ln(\mu^2x^2)}{x^2}\right] = -\frac12\Box\frac1{x^2}.$$

In four-dimensional Euclidean distribution theory,

$$\Box\frac1{x^2}=-4\pi^2\delta^{(4)}(x).$$

Hence the derivative is $2\pi^2\delta^{(4)}(x)$. Multiplying by the Wick-contraction factor $1/(32\pi^4)$ gives $\delta^{(4)}(x)/(16\pi^2)$.

3. Explain how omitting the transformation contact in a current Ward identity can mimic an anomaly.

Solution

The divergence of the time-ordered current product contains delta functions at every transformed insertion. If one compares the divergence only with a separated-point conservation equation, those required local terms appear as an unexplained breaking. An anomaly can be claimed only after the transformation contacts and all admissible local finite counterterms have been included and the remaining breaking satisfies the consistency condition but is not removable.

Continue to Operator Mixing and Renormalization Matrices to construct a closed one-insertion sector without confusing it with the contact sector. For theorem-first distribution extensions, use Scaling Degree and Distribution Extension and Epstein–Glaser Induction.

References

• Collins, John C. 1984; open-access reissue 2023. Renormalization: An Introduction to Renormalization, the Renormalization Group and the Operator-Product Expansion. Cambridge Monographs on Mathematical Physics. Cambridge University Press. DOI and Open PDF.

• Zimmermann, Wolfhart. 1973. “Normal Products and the Short Distance Expansion in the Perturbation Theory of Renormalizable Interactions.” Annals of Physics 77 (1–2): 570–601. DOI.