Feynman parametrization

Template:Short description Template:For Feynman parametrization is a technique for evaluating loop integrals which arise from Feynman diagrams with one or more loops. However, it is sometimes useful in integration in areas of pure mathematics as well. It was introduced by Julian Schwinger and Richard Feynman in 1949 to perform calculations in quantum electrodynamics.Template:Efn

Hung Cheng and T.T. Wu proved in 1987 that the sum in the Dirac delta function can be reduced to a subset of Feynman parameters. This result is known as Cheng–Wu theorem.

Formulas

Two denominators

Richard Feynman observed that

$\frac{1}{AB}={\int }_0^1\frac{du}{{[uA+(1-u)B]}^2}$

which is valid for any complex numbers A and B as long as 0 is not contained in the line segment connecting A and B. The formula helps to evaluate integrals like:

$\begin{array}{cc}\int \frac{dp}{A(p)B(p)}& =\int dp{\int }_0^1\frac{du}{{[uA(p)+(1-u)B(p)]}^2}\\ & ={\int }_0^{1}du\int \frac{dp}{{[uA(p)+(1-u)B(p)]}^2}.\end{array}$

If A(p) and B(p) are linear functions of p, then the last integral can be evaluated using substitution.

Multiple denominators

More generally, using the Dirac delta function $\delta$:

$\begin{array}{cc}\frac{1}{A_1\cdots A_n}& =(n-1)!{\int }_0^{1}d{u}_1\cdots {\int }_0^{1}d{u}_n\frac{\delta (1-\sum _{k=1}^n{u}_k)}{{\left(\sum _{k=1}^n{u}_k{A}_k\right)}^n}\\ & =(n-1)!{\int }_0^{1}d{u}_1{\int }_0^{u_1}d{u}_2\cdots {\int }_0^{u_{n-2}}d{u}_{n-1}\frac{1}{{[A_1{u}_{n-1}+A_2(u_{n-2}-u_{n-1})+\dots +A_n(1-u_1)]}^n}.\end{array}$

This formula is valid for any complex numbers A₁,...,A_n as long as 0 is not contained in their convex hull.

Even more generally, provided that $\text{Re}({\alpha }_j)>0$ for all $1\le j\le n$:

$\frac{1}{A_1^{{\alpha }_1}\cdots A_n^{{\alpha }_n}}=\frac{\mathrm{\Gamma }({\alpha }_1+\dots +{\alpha }_n)}{\mathrm{\Gamma }({\alpha }_1)\cdots \mathrm{\Gamma }({\alpha }_n)}{\int }_0^{1}d{u}_1\cdots {\int }_0^{1}d{u}_n\frac{\delta (1-\sum _{k=1}^n{u}_k)u_1^{{\alpha }_1-1}\cdots u_n^{{\alpha }_n-1}}{{\left(\sum _{k=1}^n{u}_k{A}_k\right)}^{\sum _{k=1}^n{\alpha }_k}}$

where the Gamma function $\mathrm{\Gamma }$ was used.

Derivation

$\frac{1}{AB}=\frac{1}{A-B}(\frac{1}{B}-\frac{1}{A})=\frac{1}{A-B}{\int }_B^A\frac{dz}{z^2}.$

By using the substitution $u=(z-B)/(A-B)$, we have $du=dz/(A-B)$, and $z=uA+(1-u)B$, from which we get the desired result

$\frac{1}{AB}={\int }_0^1\frac{du}{{[uA+(1-u)B]}^2}.$

In more general cases, derivations can be done very efficiently using the Schwinger parametrization. For example, in order to derive the Feynman parametrized form of $\frac{1}{A_1\mathrm{...}{A}_n}$, we first reexpress all the factors in the denominator in their Schwinger parametrized form:

$\frac{1}{A_i}={\int }_0^{\mathrm{\infty }}d{s}_i{e}^{-s_i{A}_i}\text{for }i=1,\dots ,n$

and rewrite,

$\frac{1}{A_1\cdots A_n}={\int }_0^{\mathrm{\infty }}d{s}_1\cdots {\int }_0^{\mathrm{\infty }}d{s}_n\exp (-(s_1{A}_1+\cdots +s_n{A}_n)).$

Then we perform the following change of integration variables,

$\alpha =s_1+\mathrm{...}+s_n,$

${\alpha }_i=\frac{s_i}{s_1+\cdots +s_n};i=1,\dots ,n-1,$

to obtain,

$\frac{1}{A_1\cdots A_n}={\int }_0^{1}d{\alpha }_1\cdots d{\alpha }_{n-1}{\int }_0^{\mathrm{\infty }}d\alpha {\alpha }^{n-1}\exp (-\alpha \{{\alpha }_1{A}_1+\cdots +{\alpha }_{n-1}{A}_{n-1}+(1-{\alpha }_1-\cdots -{\alpha }_{n-1})A_n\}).$

where ${\int }_0^{1}d{\alpha }_1\cdots d{\alpha }_{n-1}$ denotes integration over the region $0\le {\alpha }_i\le 1$ with ${\sum }_{i=1}^{n-1}{\alpha }_i\le 1$.

The next step is to perform the $\alpha$ integration.

${\int }_0^{\mathrm{\infty }}d\alpha {\alpha }^{n-1}\exp (-\alpha x)=\frac{{\partial }^{n-1}}{\partial (-x{)}^{n-1}}({\int }_0^{\mathrm{\infty }}d\alpha \exp (-\alpha x))=\frac{(n-1)!}{x^n}.$

where we have defined $x={\alpha }_1{A}_1+\cdots +{\alpha }_{n-1}{A}_{n-1}+(1-{\alpha }_1-\cdots -{\alpha }_{n-1})A_n.$

Substituting this result, we get to the penultimate form,

$\frac{1}{A_1\cdots A_n}=(n-1)!{\int }_0^{1}d{\alpha }_1\cdots d{\alpha }_{n-1}\frac{1}{[{\alpha }_1{A}_1+\cdots +{\alpha }_{n-1}{A}_{n-1}+(1-{\alpha }_1-\cdots -{\alpha }_{n-1})A_n{]}^n},$

and, after introducing an extra integral, we arrive at the final form of the Feynman parametrization, namely,

$\frac{1}{A_1\cdots A_n}=(n-1)!{\int }_0^{1}d{\alpha }_1\cdots {\int }_0^{1}d{\alpha }_n\frac{\delta (1-{\alpha }_1-\cdots -{\alpha }_n)}{[{\alpha }_1{A}_1+\cdots +{\alpha }_n{A}_n{]}^n}.$

Similarly, in order to derive the Feynman parametrization form of the most general case,$\frac{1}{A_1^{{\alpha }_1}\mathrm{...}{A}_n^{{\alpha }_n}}$ one could begin with the suitable different Schwinger parametrization form of factors in the denominator, namely,

$\frac{1}{A_1^{{\alpha }_1}}=\frac{1}{({\alpha }_1-1)!}{\int }_0^{\mathrm{\infty }}d{s}_1{s}_1^{{\alpha }_1-1}{e}^{-s_1{A}_1}=\frac{1}{\mathrm{\Gamma }({\alpha }_1)}\frac{{\partial }^{{\alpha }_1-1}}{\partial (-A_1{)}^{{\alpha }_1-1}}\left({\int }_0^{\mathrm{\infty }}d{s}_1{e}^{-s_1{A}_1}\right)$

and then proceed exactly along the lines of previous case.

Alternative form

An alternative form of the parametrization that is sometimes useful is

$\frac{1}{AB}={\int }_0^{\mathrm{\infty }}\frac{d\lambda }{{[\lambda A+B]}^2}.$

This form can be derived using the change of variables $\lambda =u/(1-u)$. We can use the product rule to show that $d\lambda =du/(1-u{)}^2$, then

$\begin{array}{cc}\frac{1}{AB}& ={\int }_0^1\frac{du}{{[uA+(1-u)B]}^2}\\ & ={\int }_0^1\frac{du}{(1-u{)}^2}\frac{1}{{[\frac{u}{1-u}A+B]}^2}\\ & ={\int }_0^{\mathrm{\infty }}\frac{d\lambda }{{[\lambda A+B]}^2}\end{array}$

More generally we have

$\frac{1}{A^m{B}^n}=\frac{\mathrm{\Gamma }(m+n)}{\mathrm{\Gamma }(m)\mathrm{\Gamma }(n)}{\int }_0^{\mathrm{\infty }}\frac{{\lambda }^{m-1}d\lambda }{{[\lambda A+B]}^{n+m}},$

where $\mathrm{\Gamma }$ is the gamma function.

This form can be useful when combining a linear denominator $A$ with a quadratic denominator $B$, such as in heavy quark effective theory (HQET).

Symmetric form

A symmetric form of the parametrization is occasionally used, where the integral is instead performed on the interval $[-1,1]$, leading to:

$\frac{1}{AB}=2{\int }_{-1}^1\frac{du}{{[(1+u)A+(1-u)B]}^2}.$

Further reading

• Michael E. Peskin and Daniel V. Schroeder, An Introduction To Quantum Field Theory, Addison-Wesley, Reading, 1995.
• Silvan S. Schweber, Feynman and the visualization of space-time processes, Rev. Mod. Phys, 58, p. 449 ,1986 doi:10.1103/RevModPhys.58.449
• Vladimir A. Smirnov: Evaluating Feynman Integrals, Springer, ISBN 978-3-54023933-8 (Dec., 2004).
• Vladimir A. Smirnov: Feynman Integral Calculus, Springer, ISBN 978-3-54030610-8 (Aug., 2006).
• Vladimir A. Smirnov: Analytic Tools for Feynman Integrals, Springer, ISBN 978-3-64234885-3 (Jan.,2013).
• Johannes Blümlein and Carsten Schneider (Eds.): Anti-Differentiation and the Calculation of Feynman Amplitudes, Springer, ISBN 978-3-030-80218-9 (2021).
• Stefan Weinzierl: Feynman Integrals: A Comprehensive Treatment for Students and Researchers, Springer, ISBN 978-3-030-99560-7 (Jun., 2023).