# Feynman parametrization

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.

## Formulas

Richard Feynman observed that:

${\displaystyle {\frac {1}{AB}}=\int _{0}^{1}{\frac {du}{\left[uA+(1-u)B\right]^{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:

${\displaystyle \int {\frac {dp}{A(p)B(p)}}=\int dp\int _{0}^{1}{\frac {du}{\left[uA(p)+(1-u)B(p)\right]^{2}}}=\int _{0}^{1}du\int {\frac {dp}{\left[uA(p)+(1-u)B(p)\right]^{2}}}.}$

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

More generally, using the Dirac delta function ${\displaystyle \delta }$:[1]

{\displaystyle {\begin{aligned}{\frac {1}{A_{1}\cdots A_{n}}}&=(n-1)!\int _{0}^{1}du_{1}\cdots \int _{0}^{1}du_{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}du_{1}\int _{0}^{u_{1}}du_{2}\cdots \int _{0}^{u_{n-2}}du_{n-1}{\frac {1}{\left[A_{1}+u_{1}(A_{2}-A_{1})+\dots +u_{n-1}(A_{n}-A_{n-1})\right]^{n}}}.\end{aligned}}}

This formula is valid for any complex numbers A1,...,An as long as 0 is not contained in their convex hull.

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

${\displaystyle {\frac {1}{A_{1}^{\alpha _{1}}\cdots A_{n}^{\alpha _{n}}}}={\frac {\Gamma (\alpha _{1}+\dots +\alpha _{n})}{\Gamma (\alpha _{1})\cdots \Gamma (\alpha _{n})}}\int _{0}^{1}du_{1}\cdots \int _{0}^{1}du_{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 ${\displaystyle \Gamma }$ was used.[2]

## Derivation

${\displaystyle {\frac {1}{AB}}={\frac {1}{A-B}}\left({\frac {1}{B}}-{\frac {1}{A}}\right)={\frac {1}{A-B}}\int _{B}^{A}{\frac {dz}{z^{2}}}.}$

Now just linearly transform the integral using the substitution,

${\displaystyle u=(z-B)/(A-B)}$ which leads to ${\displaystyle du=dz/(A-B)}$ so ${\displaystyle z=uA+(1-u)B}$

and we get the desired result:

${\displaystyle {\frac {1}{AB}}=\int _{0}^{1}{\frac {du}{\left[uA+(1-u)B\right]^{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 ${\displaystyle {\frac {1}{A_{1}...A_{n}}}}$, we first reexpress all the factors in the denominator in their Schwinger parametrized form:

${\displaystyle {\frac {1}{A_{i}}}=\int _{0}^{\infty }ds_{i}\,e^{-s_{i}A_{i}}\ \ {\text{for }}i=1,\ldots ,n}$

and rewrite,

${\displaystyle {\frac {1}{A_{1}\cdots A_{n}}}=\int _{0}^{\infty }ds_{1}\cdots \int _{0}^{\infty }ds_{n}\exp \left(-\left(s_{1}A_{1}+\cdots +s_{n}A_{n}\right)\right).}$

Then we perform the following change of integration variables,

${\displaystyle \alpha =s_{1}+...+s_{n},}$
${\displaystyle \alpha _{i}={\frac {s_{i}}{s_{1}+\cdots +s_{n}}};\ i=1,\ldots ,n-1,}$

to obtain,

${\displaystyle {\frac {1}{A_{1}\cdots A_{n}}}=\int _{0}^{1}d\alpha _{1}\cdots d\alpha _{n-1}\int _{0}^{\infty }d\alpha \ \alpha ^{N-1}\exp \left(-\alpha \left\{\alpha _{1}A_{1}+\cdots +\alpha _{n-1}A_{n-1}+\left(1-\alpha _{1}-\cdots -\alpha _{n-1}\right)A_{n}\right\}\right).}$

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

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

${\displaystyle \int _{0}^{\infty }d\alpha \ \alpha ^{n-1}\exp(-\alpha x)={\frac {\partial ^{n-1}}{\partial (-x)^{n-1}}}\left(\int _{0}^{\infty }d\alpha \exp(-\alpha x)\right)={\frac {\left(n-1\right)!}{x^{n}}}.}$

where we have defined ${\displaystyle x=\alpha _{1}A_{1}+\cdots +\alpha _{n-1}A_{n-1}+\left(1-\alpha _{1}-\cdots -\alpha _{n-1}\right)A_{n}.}$

Substituting this result, we get to the penultimate form,

${\displaystyle {\frac {1}{A_{1}\cdots A_{n}}}=\left(n-1\right)!\int _{0}^{1}d\alpha _{1}\cdots d\alpha _{n-1}{\frac {1}{[\alpha _{1}A_{1}+\cdots +\alpha _{n-1}A_{n-1}+\left(1-\alpha _{1}-\cdots -\alpha _{n-1}\right)A_{n}]^{n}}},}$

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

${\displaystyle {\frac {1}{A_{1}\cdots A_{n}}}=\left(n-1\right)!\int _{0}^{1}d\alpha _{1}\cdots \int _{0}^{1}d\alpha _{n}{\frac {\delta \left(1-\alpha _{1}-\cdots -\alpha _{n}\right)}{[\alpha _{1}A_{1}+\cdots +\alpha _{n}A_{n}]^{n}}}.}$

Similarly, in order to derive the Feynman parametrization form of the most general case, :${\displaystyle {\frac {1}{A_{1}^{\alpha _{1}}...A_{n}^{\alpha _{n}}}}}$ one could begin with the suitable different Schwinger parametrization form of factors in the denominator, namely,

${\displaystyle {\frac {1}{A_{1}^{\alpha _{1}}}}={\frac {1}{\left(\alpha _{1}-1\right)!}}\int _{0}^{\infty }ds_{1}\,s_{1}^{\alpha _{1}-1}e^{-s_{1}A_{1}}={\frac {1}{\Gamma (\alpha _{1})}}{\frac {\partial ^{\alpha _{1}-1}}{\partial (-A_{1})^{\alpha _{1}-1}}}\left(\int _{0}^{\infty }ds_{1}e^{-s_{1}A_{1}}\right)}$

and then proceed exactly along the lines of previous case.

## Alternative Form

An alternative form of the parametrisation that is sometimes useful is

${\displaystyle {\frac {1}{AB}}=\int _{0}^{\infty }{\frac {d\lambda }{\left[\lambda A+B\right]^{2}}}.}$

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

{\displaystyle {\begin{aligned}{\frac {1}{AB}}&=\int _{0}^{1}{\frac {du}{\left[uA+(1-u)B\right]^{2}}}\\&=\int _{0}^{1}{\frac {du}{(1-u)^{2}}}{\frac {1}{\left[{\frac {u}{1-u}}A+B\right]^{2}}}\\&=\int _{0}^{\infty }{\frac {d\lambda }{\left[\lambda A+B\right]^{2}}}\\\end{aligned}}}

More generally we have

${\displaystyle {\frac {1}{A^{m}B^{n}}}={\frac {\Gamma (m+n)}{\Gamma (m)\Gamma (n)}}\int _{0}^{\infty }{\frac {\lambda ^{m-1}d\lambda }{\left[\lambda A+B\right]^{n+m}}},}$

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

This form can be useful when combining a linear denominator ${\displaystyle A}$ with a quadratic denominator ${\displaystyle 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 ${\displaystyle [-1,1]}$, leading to:

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

## References

1. ^ Weinberg, Steven (2008). The Quantum Theory of Fields, Volume I. Cambridge: Cambridge University Press. p. 497. ISBN 978-0-521-67053-1.
2. ^ Kristjan Kannike. "Notes on Feynman Parametrization and the Dirac Delta Function" (PDF). Archived from the original (PDF) on 2007-07-29. Retrieved 2011-07-24.