Skip to content

Integration by Parts and Equation-of-Motion Redundancy

Integration by parts (IBP) and equation-of-motion (EOM) reduction are statements about a declared action and observable, not literal identities between local densities. IBP transfers derivatives and may expose a surface term. EOM operators are removed by a local perturbative field redefinition, which preserves on-shell amplitudes through the working order only after coefficients, sources, Jacobians, and induced higher-order interactions are handled consistently.

Required background. From Operator Lists to Independent Bases defines the invariant candidate space and its quotient by redundancy relations. Local versus Integrated Operator Redundancies distinguishes equality of integrated interactions from equality of local insertions. Helpful background. Boundaries, Variations, and Well-Posed Actions explains when a variational surface term vanishes and when it is physical.

Integration by parts is an action-level relation

Section titled “Integration by parts is an action-level relation”

For ordinary scalar functions on a spacetime region MM,

Md4xAB=Md4x(μA)(μB)+MdΣμAμB.\int_M d^4x\,A\Box B =-\int_M d^4x\,(\partial_\mu A)(\partial^\mu B) +\int_{\partial M}d\Sigma_\mu\,A\partial^\mu B.

Thus ABA\Box B and (A)(B)-(\partial A)\cdot(\partial B) define the same bulk interaction only if the final surface integral is absent, fixed, or otherwise irrelevant to the target observable. They remain different local insertions. On a physical boundary, at a defect, or in an asymptotic sector with nonvanishing flux, the last term belongs to the problem and cannot be discarded.

IBP relations at a fixed EFT order can be generated systematically. Enumerate every local vector JμJ^\mu with one unit lower canonical dimension and the desired internal quantum numbers, take μJμ\partial_\mu J^\mu, expand it in the ordered candidate list, and row-reduce the resulting coefficient matrix. For gauge-covariant building blocks, use the covariant adjoint relation appropriate to the representation. Reordering covariant derivatives is a separate algebraic step because

[Dμ,Dν]Φ[D_\mu,D_\nu]\Phi

produces a field-strength insertion rather than zero.

The boundary condition is part of the relation record. Saying only “drop total derivatives” silently changes the theory whenever the surface term carries charge, edge dynamics, or a source response.

Equation-of-motion operators are field-coordinate directions

Section titled “Equation-of-motion operators are field-coordinate directions”

Let the EFT expansion be

S[Φ]=S0[Φ]+ϵS1[Φ]+ϵ2S2[Φ]+,S[\Phi]=S_0[\Phi]+\epsilon S_1[\Phi] +\epsilon^2S_2[\Phi]+\cdots,

where ϵ\epsilon denotes a power of 1/Λ1/\Lambda or another controlled small parameter. Under the local perturbative redefinition

Φi=Φi+ϵFi[Φ],\Phi_i=\Phi_i'+\epsilon F_i[\Phi'],

functional Taylor expansion gives

S[Φ]=S0[Φ]+ϵ[S1[Φ]+FiδS0δΦi]+ϵ2[S2[Φ]+FiδS1δΦi+12FiFjδ2S0δΦiδΦj]+O(ϵ3),\begin{aligned} S[\Phi] ={}&S_0[\Phi']\\ &+\epsilon\left[ S_1[\Phi'] +F_i\frac{\delta S_0}{\delta\Phi_i'} \right]\\ &+\epsilon^2\left[ S_2[\Phi'] +F_i\frac{\delta S_1}{\delta\Phi_i'} +\frac12F_iF_j \frac{\delta^2S_0}{\delta\Phi_i'\delta\Phi_j'} \right] +O(\epsilon^3), \end{aligned}

with repeated field labels including spacetime integration. Consequently, an order-ϵ\epsilon interaction proportional to FiδS0/δΦiF_i\delta S_0/\delta\Phi_i can be cancelled by a redefinition with the opposite sign. This is why the leading-order EOM generates the redundancy relation at the next order.

The second line also exposes the limitation. At order ϵ2\epsilon^2, the terms involving δS1\delta S_1 and the second variation of S0S_0 are real induced interactions. Substituting a higher-order EOM into the Lagrangian does not generally reproduce them. A field redefinition is an order-by-order change of coordinates on theory space; “set the EOM to zero” is only its first-order shorthand. This distinction and its consequences for matching are derived in Criado and Pérez-Victoria 2019, §§ 2–3, pp. 5–12, and § 5.1, pp. 16–18, Open PDF.

For a path integral with sources, the same transformation changes

JiΦi=JiΦi+ϵJiFi[Φ].J_i\Phi_i =J_i\Phi_i'+\epsilon J_iF_i[\Phi'].

It also produces the functional Jacobian det(δΦ/δΦ)\det(\delta\Phi/\delta\Phi'). For local perturbative transformations, the Jacobian can be represented by local ghost interactions and is especially simple in dimensional regularization, where the corresponding closed ghost loops vanish under the usual assumptions. That convenience is regulator- and transformation-dependent; it is not permission to omit the Jacobian without checking. Counterterms and gauge fixing must be transformed consistently as well.

The quotient map summarizes what survives after the action-level relations are imposed. The construction of operator classes modulo IBP and EOM, including the distinction between a count and explicit representatives, is developed in Henning et al. 2016, §§ 1–2, preprint pp. 1–7, Open PDF. Read the figure’s right panel first when the target is not an ordinary on-shell scattering amplitude.

Candidate invariant operators are quotiented by declared integration-by-parts, equation-of-motion, field-redefinition, and algebraic relations before normalized representatives are chosen; off-shell sources and physical boundaries require retaining extra terms.

An operator basis is a normalized section of a quotient, not the unreduced candidate list or its count. Panel (a) forms classes in Vd,qinv/R\mathcal V_{d,\mathbf q}^{\mathrm{inv}}/\mathcal R, chooses representatives, and tests spanning and independence separately. Panel (b) limits the reduction: on-shell observables with boundary conditions that remove total derivatives use the quotient, whereas off-shell Green functions, explicit sources and contact terms, or physical boundaries can require the extra operators. The diagram is schematic and not to scale.

The local reason for the off-shell exception is the Schwinger–Dyson identity. For any insertion XX and local functional FiF_i,

FiδSδΦiX=iδ(FiX)δΦi,\left\langle F_i\frac{\delta S}{\delta\Phi_i}X \right\rangle =i\left\langle \frac{\delta(F_iX)}{\delta\Phi_i} \right\rangle,

up to the stated path-integral convention. The right side consists of coincident-point contact terms. It vanishes in separated on-shell matrix elements after LSZ reduction under the equivalence-theorem hypotheses, but not in a general sourced Green function. Arzt proves the quantum EOM reduction while keeping precisely this on-shell/off-shell distinction in Arzt 1995, § 2, preprint pp. 4–7, Open PDF. This is the perturbative EFT realization of the field-coordinate equivalence established in Kamefuchi, O’Raifeartaigh, and Salam 1961, pp. 529–549.

The resulting decision is concrete:

  • For an on-shell SS-matrix element with ordinary asymptotic fields and a vanishing surface term, reduce by IBP and leading EOM and transform the Wilson coefficients.
  • For off-shell correlators, transform the sources and retain the induced contact terms, or work in an enlarged Green-function basis.
  • For a physical boundary, retain the surface operator and specify its boundary condition or boundary coupling.
  • At the next EFT order, include every interaction generated by the earlier field redefinition before reducing again.

First application: derivative scalar operators

Section titled “First application: derivative scalar operators”

Return to the four-dimensional Z2\mathbb Z_2-even scalar theory, now writing the dimension-six action explicitly:

L=12(ϕ)2λ4!ϕ4+1Λ2[c1ϕ6+c2ϕ2(ϕ)2+c3ϕ3ϕ]+O(Λ4).\begin{aligned} \mathcal L ={}&\frac12(\partial\phi)^2 -\frac{\lambda}{4!}\phi^4\\ &+\frac1{\Lambda^2}\left[ c_1\phi^6 +c_2\phi^2(\partial\phi)^2 +c_3\phi^3\Box\phi \right] +O(\Lambda^{-4}). \end{aligned}

IBP gives

ϕ2(ϕ)213ϕ3ϕ,\phi^2(\partial\phi)^2 \simeq-\frac13\phi^3\Box\phi,

so define cˉ3=c3c2/3\bar c_3=c_3-c_2/3. After this action-level step, the dimension-six interaction is c1ϕ6+cˉ3ϕ3ϕc_1\phi^6+\bar c_3\phi^3\Box\phi.

Now set

ϕ=ϕ+aΛ2(ϕ)3.\phi=\phi'+\frac{a}{\Lambda^2}(\phi')^3.

The lower-order action changes by

δL0aΛ2[(ϕ)3ϕ+λ6(ϕ)6].\delta\mathcal L_0 \simeq-\frac{a}{\Lambda^2}\left[ (\phi')^3\Box\phi' +\frac{\lambda}{6}(\phi')^6 \right].

Choosing a=cˉ3a=\bar c_3 removes the derivative operator and leaves

Lred(6)=c1,redΛ2(ϕ)6,c1,red=c1+λ18c2λ6c3.\mathcal L^{(6)}_{\mathrm{red}} =\frac{c_{1,\mathrm{red}}}{\Lambda^2}(\phi')^6, \qquad c_{1,\mathrm{red}} =c_1+\frac{\lambda}{18}c_2 -\frac{\lambda}{6}c_3.

This reproduces the relation-matrix result on the preceding page. The transformation of the dimension-six terms themselves starts at Λ4\Lambda^{-4}; it may be omitted here because the action has been declared only through Λ2\Lambda^{-2}. A calculation through Λ4\Lambda^{-4} must retain it and add the required second-order field redefinition.

The same result can be checked without using the interacting example. Take a free massive scalar with EOM E(ϕ)=(+m2)ϕE(\phi)=(\Box+m^2)\phi and the redundant interaction

OE=ϕ3(+m2)ϕ.O_E=\phi^3(\Box+m^2)\phi.

After symmetrizing over four identical external legs, its tree-level contact vertex is proportional to

i=14(pi2m2).-\sum_{i=1}^4(p_i^2-m^2).

Every external leg satisfies pi2=m2p_i^2=m^2, so the on-shell amplitude vanishes. Off shell, the same polynomial is nonzero and supplies exactly the inverse-propagator contact structure predicted by the Schwinger–Dyson identity. This is a direct observable-domain check, not merely a formal substitution.

Every reduction should carry enough information for another calculation to reconstruct the same quotient and translate coefficients into or out of it. The following semantic table is the chapter-wide minimum; later pages reuse it while specializing the relevant rows.

RecordDeclare before reductionVerification retained with the result
Field content and orderSpacetime dimension, dynamical fields, exact symmetries, charges, EFT grading, and truncationEvery candidate and relation has the declared labels and order
Flavor, Hermiticity, and CPFlavor-index ranges, conjugation rule, coefficient reality conditions, and CP conventionConjugate completion and independent real parameter count agree
Operator definitionOrdered names, explicit index contractions, derivative placement, signs, and normalization factorsEach symbolic or numerical column maps to one unambiguous operator
Renormalization dataRegulator, subtraction scheme, gauge convention when relevant, renormalization scale μ\mu, and coupling definitionsCoefficients and matrix elements use the same scheme and scale
Dimensional identitiesDimension used for Lorentz and spinor algebra, γ5\gamma_5 prescription when present, and evanescent-operator definitionsThe renormalized basis closes before any four-dimensional projection
Redundancy generatorsIBP currents and boundary conditions, lower-order EOM, field maps, and algebraic identitiesEvery relation row is reproducible from a displayed generator
Basis mapCandidate and reduced dimensions, matrix orientation, exact rank, pivots, and representative orderingNullities and ranks satisfy the quotient dimension and no pivot is tolerance-dependent
Coefficient mapDual transformation, transpose convention, finite shifts, and perturbative orderCTOC^TO is unchanged through the retained order
Implementation identitySource or notebook version, dependency versions, input hash, and output checksumA clean rerun reproduces the ordered map and checksum
Round trip and physicsForward and inverse maps on the common subspace plus one amplitude, correlator, or counting benchmarkThe round trip is the identity and the benchmark is basis independent to the stated tolerance

For the scalar example, the candidate dimension is 33, the relation rank is 22, and the reduced dimension is 11. The exact map is

(c1,c2,c3)c1+λ18c2λ6c3.(c_1,c_2,c_3) \longmapsto c_1+\frac{\lambda}{18}c_2-\frac{\lambda}{6}c_3.

Its on-shell benchmark is the vanishing four-point matrix element of OEO_E above. Its failure records are equally important: nonzero boundary flux, untransformed sources, or a requested O(Λ4)O(\Lambda^{-4}) prediction blocks the naive elimination.

Writing ϕ=\Box\phi=\cdots inside every operator. EOM replacement is a controlled field redefinition at a specific EFT order. Apply it to the higher-order sector using the lower-order action, then record the induced coefficient map.

Forgetting that IBP and derivative commutation are different. Moving a covariant derivative and swapping two covariant derivatives are separate operations; the latter can generate a field strength.

Checking only the transformed action at first order. The same transformation changes sources, counterterms, the measure, and higher-order interactions. Which pieces matter is fixed by the observable, regulator, and target accuracy.

Using an off-shell zero as an on-shell theorem. An EOM insertion equals contact terms in a Green function. It disappears from the SS matrix only after the complete equivalence-theorem argument.

Derive the scalar coefficient shift without first integrating by parts.

Solution

Substituting ϕ=ϕ+a(ϕ)3/Λ2\phi=\phi'+a(\phi')^3/\Lambda^2 directly into L0\mathcal L_0 gives

δL0=aΛ2[3(ϕ)2(ϕ)2λ6(ϕ)6].\delta\mathcal L_0 =\frac{a}{\Lambda^2}\left[ 3(\phi')^2(\partial\phi')^2 -\frac{\lambda}{6}(\phi')^6 \right].

Choose a=c2/3a=-c_2/3 to remove O2O_2. This shifts c1c_1 to c1+λc2/18c_1+\lambda c_2/18 and leaves c3c_3. A second description using the IBP relation then removes O3O_3 and shifts c1c_1 by λc3/6-\lambda c_3/6, giving the same c1,redc_{1,\mathrm{red}}.

Why does a single insertion of FδS/δϕF\,\delta S/\delta\phi need not vanish in a two-point Green function even though it is redundant for the SS matrix?

Solution

The Schwinger–Dyson identity differentiates both FF and the inserted fields. Functional derivatives of the inserted fields produce delta-function contact terms, so the off-shell correlator changes. LSZ multiplication by inverse propagators and the consistent transformation of interpolating fields remove the redundant direction from on-shell amplitudes under the theorem’s hypotheses.

  • Arzt, Christopher. “Reduced Effective Lagrangians.” Physics Letters B 342, no. 1–4 (1995): 189–195. DOI; Open PDF
  • Criado, Juan Carlos, and Manuel Pérez-Victoria. “Field Redefinitions in Effective Theories at Higher Orders.” Journal of High Energy Physics 2019, no. 3 (2019): 038. DOI; Open PDF
  • Henning, Brian, Xiaochuan Lu, Tom Melia, and Hitoshi Murayama. “Hilbert Series and Operator Bases with Derivatives in Effective Field Theories.” Communications in Mathematical Physics 347, no. 2 (2016): 363–388. DOI; Open PDF
  • Kamefuchi, S., L. O’Raifeartaigh, and A. Salam. “Change of Variables and Equivalence Theorems in Quantum Field Theories.” Nuclear Physics 28 (1961): 529–549. DOI