Skip to content

Feynman and Schwinger Parameters

Parameterization converts a product of propagators into one denominator raised to a higher power. Schwinger parameters make Gaussian loop integration and graph homogeneity transparent; Feynman parameters quotient out their common scale and leave an integral over a simplex. Neither transformation discards the causal boundary value.

Required background. Anatomy of a Loop Integral supplies the denominators, powers, routing, and i0i0 data that are being combined.

For ReA>0\operatorname{Re}A>0 and Reν>0\operatorname{Re}\nu>0, the Euclidean Schwinger identity is

1Aν=1Γ(ν)0dααν1eαA.\frac{1}{A^\nu}=\frac{1}{\Gamma(\nu)} \int_0^\infty \mathrm d\alpha\,\alpha^{\nu-1}e^{-\alpha A}.

In a Lorentzian integral it is used after the Feynman i0i0 has selected a Wick rotation, or with the corresponding oscillatory representation and its convergence factor. Dropping that prescription before rotating loses the sheet information.

For two denominator factors, introduce α1=ρx\alpha_1=\rho x and α2=ρ(1x)\alpha_2=\rho(1-x). The Jacobian is ρ\rho, and the ρ\rho integral gives

1AaBb=Γ(a+b)Γ(a)Γ(b)01dxxa1(1x)b1[xA+(1x)B]a+b.\frac{1}{A^aB^b} =\frac{\Gamma(a+b)}{\Gamma(a)\Gamma(b)} \int_0^1\mathrm d x\, \frac{x^{a-1}(1-x)^{b-1}} {[xA+(1-x)B]^{a+b}}.

For nn factors this becomes an integral over xj0x_j\ge0 with jxj=1\sum_jx_j=1. The gamma-function normalization is fixed by the Dirichlet integral; checking the special case A=BA=B immediately recovers A(a+b)A^{-(a+b)}.

The Schwinger and Feynman constructions, including the projective nature of the latter, are derived in Weinzierl 2022, §§2.5.2–2.5.3, pp. 43–55.

Apply the two-denominator identity to

D1=2m12+i0,D2=(+p)2m22+i0.D_1=\ell^2-m_1^2+i0, \qquad D_2=(\ell+p)^2-m_2^2+i0.

With weight xx on D2D_2, their combination is

(1x)D1+xD2=(+xp)2Δ(x)+i0,(1-x)D_1+xD_2 =(\ell+xp)^2-\Delta(x)+i0,

where

Δ(x)=(1x)m12+xm22x(1x)p2.\Delta(x)=(1-x)m_1^2+xm_2^2-x(1-x)p^2.

The shift k=+xpk=\ell+xp is valid in a translation-invariant regulator. Thus

B(p2)=μ2ϵ01dxddk(2π)d1[k2Δ(x)+i0]2.B(p^2)=\mu^{2\epsilon}\int_0^1\mathrm d x \int\frac{\mathrm d^d k}{(2\pi)^d} \frac{1}{[k^2-\Delta(x)+i0]^2}.

This representation makes the denominator geometry visible. In the Euclidean region p2<0p^2<0 with positive masses, Δ(x)>0\Delta(x)>0 throughout the interval. As p2p^2 is continued, a zero can enter the integration domain. At the normal threshold the quadratic has a double zero, and beyond it the logarithm generated by the loop integration samples Δi0\Delta-i0 on the lower side of its cut.

The double-zero conditions Δ(x)=Δ(x)=0\Delta(x_*)=\Delta'(x_*)=0 are solved, for positive masses and x(0,1)x_*\in(0,1), by

x=m1m1+m2,p2=(m1+m2)2.x_* = \frac{m_1}{m_1+m_2}, \qquad p^2=(m_1+m_2)^2.

The algebraic companion (m1m2)2(m_1-m_2)^2 belongs to a different parameter or sheet configuration and is not the normal positive-parameter threshold. The domain condition is what prevents a root of the discriminant from being mistaken for an accessible physical pinch.

The same bubble calculation and its massless specialization appear in Weinzierl 2022, §2.5.3, pp. 46–52.

After Gaussian integration of an LL-loop scalar graph with N=jνjN=\sum_j\nu_j, the projective form has the schematic structure

IΓ ⁣(NLd2)xj0 ⁣ΩUN(L+1)d/2(Fi0U)NLd/2jxjνj1.I\propto \Gamma\!\left(N-\frac{Ld}{2}\right) \int_{x_j\ge0}\!\Omega\, \frac{\mathcal U^{N-(L+1)d/2}} {(\mathcal F-i0\,\mathcal U)^{N-Ld/2}} \prod_j x_j^{\nu_j-1}.

Ω\Omega denotes the projective parameter measure in the simplex gauge jxj=1\sum_jx_j=1. U\mathcal U is homogeneous of degree LL and records spanning-tree information; F\mathcal F is homogeneous of degree L+1L+1 and contains masses and external invariants. In this gauge the momentum-space denominator k2F/U+i0k^2-\mathcal F/\mathcal U+i0 produces Fi0U\mathcal F-i0\,\mathcal U. In a fully homogeneous projective expression the infinitesimal is instead written i0Ujxj-i0\,\mathcal U\sum_jx_j; some references absorb it into the definition of F\mathcal F. Overall prefactors depend on the chosen measure normalization, but the homogeneities and boundary value do not. Parameter boundaries xj=0x_j=0 correspond to contracted lines and are natural places to inspect UV or IR subregions. Interior zeros of F\mathcal F can participate in threshold pinches.

A complementary derivation and concise statement of these polynomials is given in Abreu, Britto, and Duhr 2022, §1.2, pp. 5–8.

Shifting before regulating. Completing the square is a change of variables in a defined regulated integral. For a cutoff that is not translation invariant, the shifted integration region contributes a surface term.

Reading every zero as a physical singularity. A zero of the parameter polynomial is only part of the pinch analysis. The domain, derivative conditions, i0i0 prescription, and possible numerator cancellation still matter.

Removing the common scale twice. Either retain all positive Schwinger parameters or fix one projective condition such as xj=1\sum x_j=1. Mixing the two without the correct Jacobian changes the normalization.

  1. Set a=b=1a=b=1 and A=BA=B in the Feynman identity. The parameter integral is one and the result is 1/A21/A^2.
  2. For m1=m2=mm_1=m_2=m, minimize Δ(x)=m2x(1x)p2\Delta(x)=m^2-x(1-x)p^2. Its first zero occurs at x=1/2x=1/2 and p2=4m2p^2=4m^2, the two-particle threshold.
  • Abreu, Samuel, Ruth Britto, and Claude Duhr. “The SAGEX Review on Scattering Amplitudes, Chapter 3: Mathematical Structures in Feynman Integrals.” Journal of Physics A: Mathematical and Theoretical 55 (2022): 443004. doi:10.1088/1751-8121/ac87de.
  • Weinzierl, Stefan. Feynman Integrals. Cham: Springer, 2022. doi:10.1007/978-3-030-99558-4. Open PDF.