Massless planar double box
The content on this page was written by AI under human supervision.
The planar two-loop double box with seven massless propagators and four massless on-shell legs, expanded in the dimensional regulator from the leading pole through the finite part. Smirnov gave its closed form in 1999 as harmonic polylogarithms of a single ratio of invariants. It is recomputed here without using that answer.
The integral
The double box is a scalar Feynman integral of the 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 $e$ of the graph $G$, the $l_\ell$ are the $L$ loop momenta, $q_e$ is the momentum flowing through line $e$ and $m_e$ its mass; here $L=2$ and every $m_e$ vanishes. With no internal mass to set a scale, the integral depends on the two Mandelstam invariants, apart from an overall power of $s$, only through the dimensionless ratio
$$x=\frac{t}{s},\qquad s=(p_1+p_2)^2,\quad t=(p_2+p_3)^2.$$
The quantity computed is the Laurent expansion in $\varepsilon$ of the normalized integral $I$, obtained from $I_G$ by dividing out the loop-measure factor $(i\pi^{d/2})^2$, the powers $(-s)^{-2\varepsilon}$ and $e^{-2\gamma_E\varepsilon}$ (with $\gamma_E$ Euler's constant) and fixing the overall normalization so that the coefficient of $\varepsilon^{-4}$ is the leading singularitythe rational prefactor of the integral, computed by putting every propagator on shell; dividing by it leaves a function of uniform transcendental weight $1/(-2s^2t)$. The expansion runs from order $\varepsilon^{-4}$ through the finite part, as a function of $x$ in the Euclidean region $s,t\lt0$, where $x\gt0$ and $I$ is real.
At a glance
- Process or family: planar two-loop double box, seven massless propagators, four massless on-shell legs
- Loops and legs: two loops; four-point function of $x=t/s$
- Master integrals: 8
- Function class: polylog; harmonic polylogarithms of $x$, uniform weight at every order
- Singular points or alphabet: letters $\{x,\,1+x\}$; branch points $x=0$ and $x=-1$
- Status: known result, rederived
- Paper: Polylogarithmic, elliptic, and Calabi–Yau Feynman integrals from a hybrid bootstrap, section 4.2
Smirnov evaluated this integral analytically in 1999 as an expansion in $\varepsilon$ whose coefficients are harmonic polylogarithmsiterated integrals of the kernels $dx/x$, $dx/(1+x)$ and $dx/(1-x)$; only the first two occur here of $x$, the first analytic result for a massless on-shell double box. The non-planar double box of Tausk followed the same year. The integral has since become a standard example for differential equations in canonical form. The family has eight master integralsa finite set of integrals of the family in terms of which every other member is a rational linear combination, by integration-by-parts identities; in the canonical basis of Henn (2013) the differential equation in $x$ has logarithmic singularities only at $x=0$ and $x=-1$, and the coefficient of $\varepsilon^j$ is a pure functionan iterated integral of uniform transcendental weight with constant rational coefficients, no rational or algebraic functions of the kinematics multiplying it of uniform transcendental weight $4+j$.
Both features, the two singular points and the uniform weight, follow from the graph alone, before any coefficient is computed. The maximal cutthe integral with every propagator put on shell, which exposes the algebraic prefactor and the geometry is rational, so the leading singularity is the rational prefactor divided out above and no elliptic curve or higher period appears. Dividing it out leaves functions built from the two letters $x$ and $1+x$, which vanish at $x=0$ (where $t=0$) and at $x=-1$ (where the third Mandelstam invariant $u=-s-t$ vanishes). Three conditions then fix which iterated integrals (words in the two letters) may appear, stated in terms of the symbolthe tensor of letters that encodes the iterated-integral structure of a polylogarithmic function; integrability is the condition that such a tensor comes from an actual function as reviewed by Duhr: the first letter of every word must be a letter whose vanishing is a Landau singularity of the graph (the first-entry condition, which both letters satisfy here); the symbol must be integrable; and discontinuities may occur only in the $s$ and $t$ channels, with none in $u$. Together with their products with zeta values of complementary weight, the words that satisfy these conditions span an ansatz of dimension $(1,2,5,11,23)$ at weights $0$ through $4$, so twenty-three unknown rational coefficients at weight four. Fitting that ansatz to high-precision numerical values of the integral at Euclidean points with an integer-relation algorithm (PSLQ) determines every coefficient as an exact rational number, weight by weight.
Closed form
Every order of the expansion is the leading singularity times a rational combination of harmonic polylogarithms of $x$ and their products with zeta values. Smirnov's function $K(x,\varepsilon)=\sum_{j\ge-4}\varepsilon^{j}K^{(j)}(x)$ is related to $I$ by $I=\frac{1}{-2s^2t}\cdot\frac{K(x,\varepsilon)}{-4}$, and through the finite part
$$ \begin{aligned} K^{(-4)}&=-4,\qquad K^{(-3)}=5\,G(0;x),\qquad K^{(-2)}=-4\,G(0,0;x)+15\,\zeta_2,\\ K^{(-1)}&=-4\,G(0,0,0;x)+4\,G(-1,0,0;x)-33\,\zeta_2\,G(0;x)+12\,\zeta_2\,G(-1;x) +\tfrac{65}{3}\,\zeta_3,\\ K^{(0)}&=32\,G(0,0,0,0;x)-16\,G(-1,0,0,0;x)-20\,G(0,-1,0,0;x)+4\,G(-1,-1,0,0;x)\\ &\quad+72\,\zeta_2\,G(0,0;x)-20\,\zeta_2\,G(-1,0;x)-60\,\zeta_2\,G(0,-1;x)+12\,\zeta_2\,G(-1,-1;x) \\ &\quad-\tfrac{88}{3}\,\zeta_3\,G(0;x)+4\,\zeta_3\,G(-1;x)+87\,\zeta_4 \end{aligned} $$
where $G(a_1,\dots,a_n;x)$ is the Goncharov iterated integral with indices $a_i=0$ and $a_i=-1$ standing for the integration kernels $dx/x$ and $dx/(1+x)$; with these two indices only, the $G$ functions are the harmonic polylogarithms of Smirnov's answer. In the weight count a word of length $n$ times $\zeta_k$ carries weight $n+k$.
Checks
Checks
| point | closed form | independent value | digits |
|---|---|---|---|
| $x=2/5$, order $\varepsilon^{0}$ | $\;\;144.4685602\ldots$ | $\;\;144.4685602\ldots$ | 90 |
| $x=2/5$, order $\varepsilon^{-1}$ | $\;\;144.4970281\ldots$ | $\;\;144.4970281\ldots$ | 97 |
| $x=2/5$, order $\varepsilon^{-2}$ | $\;\;64.04595960\ldots$ | $\;\;64.04595960\ldots$ | 104 |
| $x=3/4$, order $\varepsilon^{0}$ | $-7.244798182\ldots$ | $-7.244798182\ldots$ | 89 |
| $x=3/4$, order $\varepsilon^{-1}$ | $\;\;37.46122097\ldots$ | $\;\;37.46122097\ldots$ | 97 |
| $x=3/4$, order $\varepsilon^{-2}$ | $\;\;31.33814886\ldots$ | $\;\;31.33814886\ldots$ | 104 |
The entries are the Laurent coefficients in $\varepsilon$ of $e^{-2\gamma_E\varepsilon}K(x,\varepsilon)/x$; the independent values are auxiliary-mass-flow evaluations (AMFlow) of the double box at the two Euclidean points $s=-1$, $t=-x$, which were not used in the fit.
Evaluator
Evaluator
- massless-dbox-evaluate.py: the $\varepsilon$ expansion of the massless double box from $\varepsilon^{-4}$ through $\varepsilon^0$ at any Euclidean $x=t/s\gt0$, compared with the Checks values.
python3 massless-dbox-evaluate.py --dps 110 - The polylogarithmic-suite bundle (zip, 473 KB) — the data file the evaluator reads (words, exact rational coefficients, normalization and auxiliary-mass-flow values), the companion evaluators of the same suite, vendored helper modules and checksums; the zip also contains the script, so run it from the unzipped folder.
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
| tool | role |
|---|---|
| Landau Alphabet | fixes the two-letter alphabet $\{x,\,1+x\}$ from the Landau singularities of the graph |
| Ansatzer | imposes first entry, integrability and the absence of a $u$-channel cut on the word space |
| PSLQ | integer-relation fits that determine the exact rational coefficients weight by weight |
| AMFlow | high-precision values at Euclidean points: the input of the fits and the independent evaluations in the Checks table |
| GPLEval | arbitrary-precision evaluation of the harmonic polylogarithms |
Same family
- The outer-mass double box has the same graph with a common mass on the six perimeter lines and a massless rung, and is finite in four dimensions: a single pure function of weight four whose alphabet contains three square roots.
- The central-mass double box puts the mass on the shared rung instead and keeps the perimeter massless; its alphabet is seven rational letters in $s$ and $t$, with no square root.
- The three-loop ladder has the same two-letter alphabet $\{x,\,1+x\}$ one loop higher: a three-point ladder with nine massless lines and two off-shell legs, whose finite part has weight six.
The paper
- Polylogarithmic, elliptic, and Calabi–Yau Feynman integrals from a hybrid bootstrap (PDF) — Matthew D. Schwartz.
References
| V. A. Smirnov, Analytical result for dimensionally regularized massless on-shell double box, Phys. Lett. B 460 (1999) 397–404 [arXiv:hep-ph/9905323] | the closed form in harmonic polylogarithms that is rederived here |
| J. B. Tausk, Non-planar massless two-loop Feynman diagrams with four on-shell legs, Phys. Lett. B 469 (1999) 225 [arXiv:hep-ph/9909506] | the non-planar companion of this integral, solved the same year |
| J. M. Henn, Multiloop integrals in dimensional regularization made simple, Phys. Rev. Lett. 110 (2013) 251601 [arXiv:1304.1806] | the canonical form of the differential equation, with this family as the worked example: eight master integrals, letters $x$ and $1+x$ |
| C. Duhr, Mathematical aspects of scattering amplitudes, TASI 2014 lectures, in: Journeys Through the Precision Frontier (WSP, 2015) 419–476 [arXiv:1411.7538] | symbols, the first-entry condition and integrability, the conditions imposed on the words |
| A. B. Goncharov, Multiple polylogarithms and mixed Tate motives (2001) [arXiv:math/0103059] | the iterated integrals $G(a_1,\dots,a_n;x)$ in which the result is written |
| 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 values in the fit and in the Checks table |