Two-loop Sudakov form factor

The content on this page was written by AI under human supervision.

The planar two-loop ladder vertex of the massless Sudakov form factor, with six massless lines, two on-shell legs and one off-shell leg. Such integrals were computed by Gonsalves (1983) and Kramer and Lampe (1987). The ladder is recomputed without using the published answer and given in closed form for arbitrary dimension $d$, as a combination of Gamma functions.

The integral

Feynman diagram of the two-loop Sudakov ladder vertex. A triangle of thin black massless lines has a wavy off-shell leg at its apex, labeled q = p_1 + p_2 with q squared equal to s. A horizontal thin massless rung joins the two slanted sides of the triangle, so that the lower part of the drawing is a four-sided box. Straight thin external legs p_1 and p_2, with momentum arrows pointing inward, enter the two bottom corners, and the label p_1 squared = p_2 squared = 0 sits below the base. All six internal lines are drawn thin: every propagator is massless.
Thin black lines are massless propagators; the wavy line is the off-shell leg $q=p_1+p_2$ with $q^2=s$, and the straight external lines are the on-shell legs $p_1$ and $p_2$ with $p_1^2=p_2^2=0$. The upper triangle and the lower box share the rung, and $s$ is the only scale of the integral.

The ladder in the figure is the planar two-loop vertex integral that enters the two-loop massless Sudakov form factor, with two massless on-shell legs, $p_1^2=p_2^2=0$, and one off-shell leg carrying $s=q^2=(p_1+p_2)^2$. It is a scalar Feynman integral of the standard form

$$I_G \;=\; \int \prod_{\ell=1}^{L} d^d l_\ell \;\prod_{e\in G}\frac{1}{q_e^2-m_e^2}\,,\qquad d=4-2\varepsilon,$$

where the product runs over the internal lines of the graph $G$, $L$ is the number of loops and $l_\ell$ the loop momenta, $q_e$ is the momentum flowing through line $e$ and $m_e$ its mass; here $L=2$, the graph has six lines and every $m_e$ vanishes. With the two loop momenta written $k_1$ and $k_2$, the seven propagators of its integral family are

$$D_1=k_1^2,\quad D_2=(k_1+p_1)^2,\quad D_3=(k_1+p_1+p_2)^2,\quad D_4=k_2^2,\quad D_5=(k_2+p_1)^2,\quad D_6=(k_2+p_1+p_2)^2,\quad D_7=(k_1-k_2)^2,$$

and $I[a_1,\dots,a_7]{}(s)$ denotes the integral of $\prod_j D_j^{-a_j}$ with the measure $d^dk_1\,d^dk_2/(i\pi^{d/2})^2$, that is, the form above with a factor $1/(i\pi^{d/2})$ for each loop, and no factor $e^{\varepsilon\gamma_E}$ (with $\gamma_E$ Euler's constant) or $(4\pi)^{\varepsilon}$. Two loops and two external momenta give seven scalar products, so for any six-line graph of the family one of the seven $D_j$ is an irreducible numerator; for the ladder that numerator is $D_2$, which is absent from the graph. The triangle is formed by $D_1$, $D_3$ and the rung $D_7$, the box by $D_4$, $D_5$, $D_6$ and $D_7$. The quantity computed is the ladder $I[1,0,1,1,1,1,1]{}(s)$ exactly in $d$.

At a glance

Planar two-loop vertex integrals of this kind were obtained in closed form long ago, by Gonsalves (1983) for the on-shell quark form factor and by Kramer and Lampe (1987) for massless QCD, and the master integralsa finite basis of integrals of the family in terms of which every other member is a rational combination in $d$ and the kinematics of the family are given as Euler $\Gamma$ functions by Gehrmann, Huber and Maître (2005) in their computation of the two-loop quark and gluon form factors. With a single scale the only logarithm that can appear is $\log(-s)$, and the whole $s$ dependence is the overall power $(-s)^{-2-2\varepsilon}$. An integration-by-parts reduction leaves three master integrals, each a product of one-loop bubble and triangle integrals and therefore a ratio of $\Gamma$ functions: the sunrise $I[0,0,1,1,0,0,1]$, the bubble-inserted triangle $I[0,1,0,1,0,1,1]$ and the product of two $q$-channel bubbles $I[1,0,1,1,0,1,0]$, which are the masters $A_3$, $A_4$ and $A_2^{2}$ of Gehrmann, Huber and Maître. The family has no six-line master integral, so the ladder itself reduces to a rational combination in $d$ of these three $\Gamma$-function masters and its $\varepsilon$ expansion can be carried to any order.

Closed form

With $G(a,b)$ the one-loop massless bubble with propagator powers $a$ and $b$, and $T(a,b,c)$ the one-loop massless triangle with two on-shell legs, $b$ being the power of the line between those legs and $a$, $c$ the powers of the other two,

$$G(a,b)=\frac{\Gamma(a+b-\tfrac d2)\,\Gamma(\tfrac d2-a)\,\Gamma(\tfrac d2-b)}{\Gamma(a)\,\Gamma(b)\,\Gamma(d-a-b)},\qquad T(a,b,c)=\frac{\Gamma(a+b+c-\tfrac d2)\,\Gamma(\tfrac d2-a-b)\,\Gamma(\tfrac d2-b-c)}{\Gamma(a)\,\Gamma(c)\,\Gamma(d-a-b-c)},$$

the ladder is, exactly in $d$,

$$I[1,0,1,1,1,1,1](s)=(-s)^{-2-2\varepsilon}\Big[-\frac{3}{\varepsilon}\,G(1,2)\,T(1,1+\varepsilon,1)+\frac{(1-2\varepsilon)^2}{\varepsilon^2}\,G(1,1)^2+\frac{3(1-3\varepsilon)}{\varepsilon^2}\,G(1,2)\,G(1+\varepsilon,1)\Big].$$

The three terms are the dotted bubble-inserted triangle, the square of the one-loop bubble and the dotted sunrise, each multiplied by its reduction coefficient; a dot (the index 2 in $G(1,2)$) marks a squared propagator, and each dotted integral is a rational multiple in $d$ of the corresponding master. At $s=-1$ the ladder multiplied by $e^{2\varepsilon\gamma_E}$ has the expansion

$$e^{2\varepsilon\gamma_E}\,I(-1)=\frac{1}{4\varepsilon^4}+\frac{5\pi^2}{24\varepsilon^2}+\frac{29\zeta_3}{6\varepsilon}+\frac{3\pi^4}{32}+\varepsilon\Big(\frac{329\zeta_5}{10}-\frac{107\pi^2\zeta_3}{36}\Big)+\varepsilon^2\Big(\frac{1723\pi^6}{60480}-\frac{833\zeta_3^2}{18}\Big)+\varepsilon^3\Big(\frac{3149\zeta_7}{14}-\frac{211\pi^2\zeta_5}{12}-\frac{73\pi^4\zeta_3}{48}\Big)+O(\varepsilon^4).$$

Here $I(-1)$ is the ladder at $s=-1$ and $\zeta_n$ are the Riemann zeta values. The $\varepsilon^{-3}$ pole is absent, and the coefficient of $\varepsilon^j$ has uniform transcendental weightthe integer grade carried by a period: $\pi$ and $\log$ have weight one, $\zeta_n$ weight $n$, and weights add under multiplication $j+4$ and is a rational combination of products of even powers of $\pi$ and the odd zeta values $\zeta_3,\zeta_5,\zeta_7,\ldots$. At any other Euclidean point the expansion is that of $s=-1$ multiplied by the expansion of $(-s)^{-2-2\varepsilon}$.

Checks

Checks
pointclosed formindependent valuedigits
$s=-1$, order $\varepsilon^{0}$$\;\;3.813564111\ldots$$\;\;3.813564111\ldots$110
$s=-1$, order $\varepsilon^{1}$$-8.349416389\ldots$$-8.349416389\ldots$109
$s=-1$, order $\varepsilon^{6}$$-677.9909672\ldots$$-677.9909672\ldots$78
$s=-2$, order $\varepsilon^{0}$$\;\;0.3603200583\ldots$$\;\;0.3603200583\ldots$109
$s=-2$, order $\varepsilon^{1}$$-2.859473574\ldots$$-2.859473574\ldots$109
$s=-2$, order $\varepsilon^{6}$$-80.72591179\ldots$$-80.72591179\ldots$78

The independent values are auxiliary-mass-flow evaluations (AMFlow) of the Laurent coefficients of the ladder $I(s)$ itself, without the factor $e^{2\varepsilon\gamma_E}$ of the expansion above, at $s=-1$ and $s=-2$.

Evaluator

Evaluator

Python 3 with mpmath; the script reads its data file from the unzipped bundle folder, so run it there. A full run takes a few seconds on a laptop.

Tools
toolrole
Kiraintegration-by-parts reduction of the ladder onto its three $\Gamma$-function master integrals
PSLQinteger-relation recognition of the Laurent coefficients at $s=-1$ as rational combinations of zeta values
AMFlowthe independent numerical evaluations in the Checks table

Same family

The nearest relatives of the vertex among the massless polylogarithmic integrals computed by the same method are two other ladders. The massless planar double box is the two-loop four-point ladder with every line massless, whose expansion coefficients are harmonic polylogarithms of one ratio of invariants, and the three-loop ladder is a three-point function of two ratios of invariants given at weight six in classical polylogarithms. The one-loop massless pentagon is the other massless integral of the set, a five-point function on a 21-letter alphabet.

The paper

References

R. J. Gonsalves, Dimensionally regularized two-loop on-shell quark form-factor, Phys. Rev. D 28 (1983) 1542 [doi:10.1103/PhysRevD.28.1542]an early closed-form computation of the two-loop on-shell vertex integrals in dimensional regularization
G. Kramer and B. Lampe, Integrals for two loop calculations in massless QCD, J. Math. Phys. 28 (1987) 945 [doi:10.1063/1.527586]closed forms for the integrals of two-loop calculations in massless QCD, including the planar vertex
T. Gehrmann, T. Huber and D. Maître, Two-loop quark and gluon form-factors in dimensional regularisation, Phys. Lett. B 622 (2005) 295–302 [arXiv:hep-ph/0507061]the three master integrals of the family, $A_2^2$, $A_3$ and $A_4$, as $\Gamma$ functions exact in $d$
X. Liu and Y.-Q. Ma, AMFlow: a Mathematica package for Feynman integrals computation via auxiliary mass flow (2022) [arXiv:2201.11669]the auxiliary-mass-flow method used for the numerical evaluations in the Checks table

← back to Feynman diagrams