---
title: Feynman Residues in Quantum Field Theory
url: https://www.emergentmind.com/topics/feynman-residues
type: topic
---

# Feynman Residues in Quantum Field Theory

“Feynman residues” denotes several related constructions attached to singularities of Feynman integrals, propagators, and amplitude representations. In perturbative quantum field theory, the term can refer to residues of loop-energy poles that define cuts, residues of higher-order poles produced by raised propagators, residues of propagator poles that carry spectral weight or LSZ normalization data, and residue-type contributions appearing in Mellin–Barnes, Grassmannian, twisted-cohomology, or configuration-space formulations [1702.03163] [1903.02286] [2312.03157] [1710.07165]. The common principle is localization on singular loci—on-shell hypersurfaces, complex poles, Landau varieties, or exceptional divisors—while the operational meaning depends on the framework.

## 1. Cut residues and multivariate on-shell localization

A precise residue-theoretic definition of a cut Feynman integral is given at one loop by interpreting a cut as a multivariate residue of the loop integrand on the locus where selected propagators are put on shell [1702.03163]. For a one-loop scalar integral
\[
I_n^D=\int \omega_n^D,\qquad \omega_n^D= \frac{e^{\gamma_E\epsilon}}{i\pi^{D/2}} \frac{d^Dk}{D_1\cdots D_n},
\]
with propagators \(D_j=(k-q_j)^2-m_j^2+i0\), and a cut set \(C\subseteq\{1,\dots,n\}\), the composed residue is
\[
\mathrm{Res}_C[\omega_n^D]=\mathrm{Res}_{S_1\cdots S_c}[\omega_n^D],
\]
where \(S_i=\{D_i=0\}\). The associated cut integral is expressed as
\[
\mathcal{C}_C I_n = (2\pi i)^{-\lceil c/2\rceil}\int_{\Gamma_C}\omega_n^D, \qquad \Gamma_C=\delta_C S_\bot.
\]
Here \(\delta_C\) is the iterated Leray coboundary, and the cut contour is the tubular cycle around the on-shell variety.

This residue definition makes the geometry of the cut explicit. The cut variety \(S_C=\bigcap_{j\in C}S_j\) is a sphere \(S_C\cong S^{D-c}\), and the cut contour is its Leray coboundary [1702.03163]. First-type Landau singularities arise at finite loop momentum and are characterized by the vanishing of the modified Cayley determinant \(Y_C\). Second-type singularities are associated with pinches at infinity and are controlled by Gram determinants in the compactified description. A notable consequence is that cuts associated with second-type singularities can be written as specific combinations of ordinary cuts, leading to linear relations among distinct cuts of the same integral.

Multidimensional residues also provide a direct method for extracting master-integral coefficients when propagators occur with generic powers [1112.4136]. For
\[
\omega = \frac{h(x)\,dx}{g_1(x)\cdots g_n(x)}, \qquad dx = dx_1\wedge\cdots\wedge dx_n,
\]
the local residue at an isolated common zero \(p\) is
\[
\operatorname{Res}_p \frac{h(x)\,dx}{g_1(x)\cdots g_n(x)} = \frac{1}{(2\pi i)^n}\int_{T_g(\varepsilon)} \frac{h(x)\,dx}{g_1(x)\cdots g_n(x)}.
\]
When \(g_i(x)=(x_i-p_i)^{a_i+1}\), the residue becomes a derivative:
\[
\operatorname{Res}_p \frac{h(x)\,dx}{\prod_{i=1}^n (x_i-p_i)^{a_i+1}}
= \frac{1}{a_1!\cdots a_n!}
\left.
\frac{\partial^{a_1+\cdots+a_n}h}
{\partial x_1^{a_1}\cdots \partial x_n^{a_n}}
\right|_{x=p}.
\]
At one loop this allows direct extraction of the scalar box coefficient from a quadruple cut, bypassing iterative IBP reduction [1112.4136].

## 2. Higher-order poles from raised propagators

Starting from two loops, self-energy insertions on internal lines generate propagators raised to higher powers, so the integrand develops higher-order poles [1903.02286]. If an internal line is represented by
\[
D_j = k_j^2 - m^2 + i\delta,
\]
a self-energy insertion produces factors such as \(1/D_j^2\), and at higher loops even higher powers can appear. Analytically, such integrals can be reduced by IBP identities, but in numerical approaches the residue of the higher-order pole becomes a direct computational obstacle.

The standard residue formula for a pole of order \(\nu\) is
\[
\mathrm{res}\!\left(f,z_0\right) = \frac{1}{(\nu-1)!}
\left.
\left(\frac{d}{dz}\right)^{\nu-1}
\Big[(z-z_0)^\nu f(z)\Big]
\right|_{z=z_0}.
\]
For \(\nu>1\), derivatives of the full integrand are required, and in the multivariate case the problem becomes still more cumbersome [1903.02286]. This is precisely the situation encountered in loop-tree duality and numerical unitarity when a raised propagator goes on shell.

For renormalised amplitudes in the on-shell scheme, the problematic residue can be made to vanish already at the integrand level by combining the self-energy graph with a suitable counterterm graph [1903.02286]. The key device is the on-shell projection
\[
p^\flat = \big(\mathrm{sign}(E)\sqrt{\vec p^{\,2}+m^2},\,\vec p\big), \qquad (p^\flat)^2=m^2,
\]
together with “flattened” denominators such as
\[
D_1^\flat=\left(k+\frac12 p^\flat\right)^2-m^2, \qquad
D_2^\flat=\left(k-\frac12 p^\flat\right)^2-m^2.
\]
The loop integrand and its counterterm are engineered so that their sum vanishes quadratically as the raised propagator goes on shell:
\[
\lim_{k^2\to m^2}(\text{Self-energy integrand}) = \mathcal O\!\left((E-E^\flat)^2\right).
\]
In that situation the double pole and the accompanying subleading pole cancel, and the residue of the raised propagator vanishes [1903.02286].

The scope of the construction is scheme dependent. The quadratic on-shell cancellation is tied to the on-shell renormalisation scheme and does not hold in \(\overline{\mathrm{MS}}\), because finite scheme-changing terms remain nonzero in the on-shell limit [1903.02286]. The same logic extends to quark and gluon self-energies in QCD, although for massive quarks and gluons some terms vanish only linearly or are proportional to \(p^\mu p^\nu\) or \(p^\flat{}^\mu p^\flat{}^\nu\); these terms are harmless when contracted with physical gauge-invariant quantities.

## 3. Nested residues, causal denominators, and Landau geometry

In the Loop-Tree Duality framework, multiloop contour integration is organized into iterated residues in the loop-energy variables [2010.12971]. For a generic rational integrand with quadratic denominators,
\[
I=\int\prod_{i=1}^L \frac{dx_i}{2\pi i}\, f(x_1,\dots,x_L),
\]
the primitive variables are complexified successively, not simultaneously. This produces residue chains with an important dichotomy: displaced poles cancel, while the surviving terms are the nested residues.

The paper proves that contributions from displaced poles cancel pairwise [2010.12971]. Those poles are generated when the position of a pole in one variable depends on the residue already taken in another variable. After explicit Laurent expansion, the relevant iterated residues have equal magnitude and opposite sign. The surviving nested residues encode the physically relevant information and are naturally mapped onto nondisjoint on-shell states. A central structural consequence is that unphysical singularities vanish, and the final expressions can be written using only causal denominators.

This cancellation pattern is closely related to Landau geometry. At one loop, first-type singularities are tied to finite-momentum pinch configurations and second-type singularities to pinch configurations at infinity [1702.03163]. In a more algebraic formulation of integral reduction, the maximal-cut Landau locus in the Baikov representation,
\[
J^\Gamma_{\text{Landau}}=\langle B\rangle+\sum_{i\in\text{ISPs}(\Gamma)}\langle \partial_i B\rangle + J_{\text{cut}^\Gamma},
\]
controls the critical syzygies used in reduction [2512.05869]. The decomposition
\[
\text{CSyz}(\Gamma)=\sum_{i=0}^{n_{\text{IR}}}\text{FittSyz}(U_i^\Gamma)
\]
expresses critical syzygies as a sum over irreducible components of the Landau locus. The paper describes the associated saturation and ideal-quotient operations as residue-like, in the sense that they remove unwanted singular contributions and retain the components relevant for reduction [2512.05869].

A plausible implication is that several apparently distinct residue constructions—cut residues, nested LTD residues, and Landau-controlled syzygies—are different localizations of the same singular-support data. The supplied literature does not collapse these viewpoints into a single formalism, but it repeatedly identifies the singular locus as the organizing object.

## 4. Propagator-pole residues, spectral weights, and LSZ factors

In many-body Green’s function theory, the residue at a propagator pole is the spectral weight carried by that pole [2312.03157]. The exact one-particle Green’s function satisfies Dyson’s equation
\[
\mathbf G(\omega) = \mathbf G^{(0)}(\omega) + \mathbf G^{(0)}(\omega)\,\mathbf\Sigma(\omega)\,\mathbf G(\omega),
\]
or equivalently
\[
\mathbf G(\omega) = \bigl[\omega\mathbf 1-\mathbf\epsilon-\mathbf\Sigma(\omega)\bigr]^{-1}.
\]
If \(\omega_q\) is a pole, the residue is
\[
F(\omega_q) \equiv \operatorname{Res}_{\omega_q} G_{qq}(\omega)
= \left[
1-\mathbf U_q^\dagger
\left(\frac{\partial \mathbf\Sigma(\omega)}{\partial\omega}\right)_{\omega_q}
\mathbf U_q
\right]^{-1},
\]
and in the diagonal approximation
\[
F(\omega_q) =
\left[
1-\left.\frac{\partial \Sigma_{qq}(\omega)}{\partial\omega}\right|_{\omega_q}
\right]^{-1}.
\]
For the exact Green’s function, these residues satisfy the sum rules quoted in the paper [2312.03157].

Finite-order Feynman–Dyson perturbation theory does not preserve this structure uniformly [2312.03157]. The second-order self-energy is described as mostly physical: principal roots have residues close to \(1\), while satellite roots have near-zero residues. At odd perturbation orders, however, the self-energy can acquire the wrong concave or convex shape within a frequency bracket, yielding complex roots or real roots with residues outside the physical range \([0,1]\). At higher even orders, numerous phantom poles can appear with essentially zero residues. The paper attributes the nonconvergence to the nonanalyticity of the rational-function exact Green’s function at many frequencies, and argues that Padé approximants can largely restore the correct pole and residue structure [2312.03157].

A distinct but related notion appears in LSZ reduction for mixed propagators. For systems with non-diagonal fermionic or scalar propagators, the pole part factorizes at each complex pole [1710.07165]. For scalars,
\[
\widetilde{G}(p,-p) = i\,\zeta\, [p^2-m^2]^{-1}\,\zeta^{\,T} +\text{[non-pole part]},
\]
while for Majorana fermions
\[
\hat{G}(p) = \hat\zeta\,[p^2-m^2]^{-1}[\,\slashed p + m\,]\,\hat\zeta^{\,T} +\text{[non-pole part]}.
\]
The matrices \(\zeta\) and \(\hat\zeta\) are the “square-rooted” residues. The paper gives all-orders prescriptions for these factors in arbitrary renormalization schemes for Majorana fermions, Dirac fermions, generic mixing fermions, and scalars [1710.07165]. In the stable-particle case they match the usual LSZ wave-function normalization logic; for unstable particles they remain useful for defining effective couplings and resonant amplitudes.

## 5. Hypergeometric, Grassmannian, and Yangian residue structures

Residues also organize exact function spaces of Feynman integrals. In conformal integral bootstrap, the \(D\)-dimensional box integral with generic propagator powers is fixed by Yangian symmetry to a specific linear combination of Appell \(F_4\) functions [1912.05561]. In Mellin–Barnes form,
\[
\phi_4=\frac{N_4}{(2\pi i)^2}\int_{\kappa_4+i\mathbb R^2} dz_1\wedge dz_2\,\omega_4,
\]
with
\[
\omega_4=u^{z_1}v^{z_2}\Gamma_{-z_1}\Gamma_{-z_2}\Gamma_{1-\gamma-z_1}\Gamma_{1-\gamma'-z_2}\Gamma_{\alpha+z_1+z_2}\Gamma_{\beta+z_1+z_2}.
\]
Closing the contour in compatible cones yields
\[
\phi_4 = N_4\sum_{\vec z^*\in R_i}\operatorname*{res}_{\vec z=\vec z^*}\omega_4.
\]
In cone \(R_{\mathrm I}\), the fundamental residue reproduces the Appell \(F_4\) coefficient structure exactly [1912.05561]. The same paper argues that Yangian invariance and Mellin–Barnes residue calculus are two descriptions of the same analytic structure: the Yangian PDEs determine the coefficient system, while the Mellin–Barnes contour expresses those coefficients as sums over residues.

In four-dimensional undeformed kinematics, the conformal box collapses to the Bloch–Wigner function [1912.05561]. In this sense, the residue sum is not merely a computational tool but an analytic mechanism selecting the physically relevant branch from a larger solution space.

A different residue geometry appears in the Grassmannian formulation of scattering amplitudes. Leading singularities are obtained as residues of the Grassmannian integral over \(G(n,k+2)\), after enough minors are set to zero to localize the integral [1012.4136]. The localization dimension is
\[
d = k(n-k-4).
\]
Composite residues arise when vanishing of one minor forces further factorization. The coordinates \((w_a,t_a)\) introduced in the paper place composite and non-composite residues on equal footing and make residue theorems more uniform [1012.4136].

For generalized biadjoint amplitudes \(m_n^{(3)}\), certain \((n-5)\)-dimensional residues are even larger objects: they reproduce the entire ordinary biadjoint partial amplitude \(m_n^{(2)}\) rather than a single diagram [2204.01743]. The main theorem states
\[
{\rm Res}[m^{(3)}_n]_{(\,t_{45\ldots n}=0,\;t_{56\ldots n}=0,\;\ldots,\;s_{n-2,n-1,n}=0\,)} = m_n^{(2)}.
\]
In the \(n=6\) example, \(m_6^{(3)}\) has \(48\) generalized Feynman-diagram terms, but only \(14\) survive at the chosen residue, and \(14\) is exactly the number of planar cubic diagrams in \(m_6^{(2)}\) [2204.01743]. Here “Feynman residue” denotes a residue whose value is a full tree-level amplitude for fixed planar ordering.

## 6. Renormalization residues, higher residue pairings, and period invariants

In configuration-space renormalization, residues appear as Poincaré residues of the pulled-back Feynman form on a wonderful compactification of the graph configuration space [1012.5485]. For a local form
\[
\omega=\frac{g(z)\,dz_1\wedge\cdots\wedge dz_n}{f(z)},
\]
the Poincaré residue along \(f=0\) is
\[
\operatorname{Res}[\omega] =
(-1)^{i-1}\left.
\frac{g(z)\,dz_1\wedge\cdots\wedge \widehat{dz_i}\wedge\cdots\wedge dz_n}{\partial f/\partial z_i}
\right|_{f=0}.
\]
If \(E_y\) is the exceptional divisor associated with a divergent subgraph \(y\), regularization is achieved by replacing the singular part of the real integration cycle with a Leray coboundary, and the ambiguity is measured by
\[
2\pi i \int_{\sigma_y} \operatorname{Res}\big[T^*(\omega_\Gamma)\big].
\]
For primitive logarithmically divergent graphs, the pulled-back form has a simple pole along the deepest exceptional divisor; with logarithmic subdivergences, iterated residues appear on intersections of exceptional divisors indexed by nests of divergent subgraphs [1012.5485].

A formally different but conceptually adjacent construction is provided by higher residue pairings in the twisted-cohomology description of Feynman integrals [1910.11852]. After Schwinger parametrization, a generic \(L\)-loop integral takes the form
\[
I_{\nu_1,\dots,\nu_P} =\int_{\mathbb R_+^P} e^{\varepsilon W}\,\varphi_{\nu_1,\dots,\nu_P},
\]
with twisted differential
\[
\nabla_{dW}=d+\varepsilon\, dW\wedge.
\]
The relevant intersection number has an \(\varepsilon^{-1}\) expansion
\[
\langle \varphi_-|\varphi_+\rangle_{dW}
=\sum_{k=0}^\infty \varepsilon^{-k}\,(\varphi_-|\varphi_+)_{dW,k},
\]
whose coefficients are Saito’s higher residue pairings. The leading term is the ordinary Grothendieck residue
\[
(\varphi_-|\varphi_+)_{dW,0} =
\operatorname{Res}_{dW=0}
\left(
\frac{\widehat\varphi_-\,\widehat\varphi_+\,d^m z}
{\partial_1W\,\partial_2W\cdots \partial_mW}
\right).
\]
Because the connection matrix is polynomial in \(\varepsilon\), finitely many higher residue pairings suffice to reconstruct the exact differential equations near four dimensions [1910.11852].

Residue language also appears in the arithmetic side of perturbation theory. Feynman periods are the renormalization-group independent parts of logarithmically divergent, subdivergence-free graphs and are equivalently the \(\frac1\varepsilon\) residues of those graphs in dimensional regularization [2206.10460]. The period is defined by
\[
P_G=\frac{I_G(Q)}{(Q^2)^{-\omega_G}}, \qquad
\omega_G=\sum_{e\in E_G}\nu_e-\frac{D}{2}h_1(G).
\]
Using graphical functions and conformal four-point integrals, all subdivergence-free Feynman periods in \(\phi^3\) theory up to six loops and \(561\) of \(607\) Feynman periods at seven loops were computed [2206.10460]. In \(\phi^4\) graph theory, the extended graph permanent was introduced as an infinite sequence of residues from prime order finite fields and shown to be preserved by completion/decompletion, planar duality, and the Schnetz twist [1704.06350]. This suggests a broad residue hierarchy ranging from analytic pole coefficients to arithmetic graph invariants.

A recurring misconception is that “residue” always means a simple coefficient extracted from a one-variable pole. In the literature summarized here, residues can instead be multivariate Leray residues, higher-order derivatives at repeated poles, spectral weights of propagator roots, factorized LSZ wave-function data, Mellin–Barnes lattice sums, Poincaré residues on exceptional divisors, or \(\frac1\varepsilon\) coefficients of renormalized amplitudes. The unifying content is not a single formula but the localization of Feynman-theoretic information on singular structures.

Source: https://www.emergentmind.com/topics/feynman-residues