This paper was converted on www.awesomepapers.org from LaTeX by an anonymous user.
Want to know more? Visit the Converter page.

gbsn

aainstitutetext: CAS Key Laboratory of Theoretical Physics, Institute of Theoretical Physics, Chinese Academy of Sciences, Beijing 100190, Chinabbinstitutetext: School of Fundamental Physics and Mathematical Sciences, Hangzhou Institute for Advanced Study and ICTP-AP, UCAS, Hangzhou 310024, Chinaccinstitutetext: Peng Huanwu Center for Fundamental Theory, Hefei 230026, Chinaddinstitutetext: School of Physical Sciences, University of Chinese Academy of Sciences, No.19A Yuquan Road, Beijing 100049, China

Notes on conformal integrals: Coulomb branch amplitudes, magic identities and bootstrap

Song He [email protected] a    Xuhang Jiang [email protected] a,d    Jiahao Liu [email protected] a,d    Yao-Qi Zhang [email protected]
Abstract

We study multi-loop conformal integrals for four-point correlators of planar 𝒩=4{\cal N}=4 super-Yang-Mills theory, and in particular those contributing to Coulomb branch amplitudes in the ten-dimensional lightlike limit, where linear combinations of such integrals are determined by the large R-charge octagons exactly known from integrability. Exploiting known results for integrands, we review those combinations of dual conformal invariant (DCI) integrals that must evaluate to determinants of ladders, generalizing the simplest cases of Basso-Dixon fishnet integrals; in this way, we summarize all-loop predictions for the integrands (which are extracted from ff-graphs) contributing to components of Coulomb branch amplitudes, such as next-to-fishnet integrals. Moreover, this exercise produces new “magic identities”, i.e. certain combinations of DCI integrals equal zero, and we enumerate and simplify such identities up to six loops explicitly.

On the other hand, most of these individual integrals have not been computed beyond three loops, and as a first step we consider a bootstrap program for DCI integrals based on their leading singularities and the space of pure functions. We bootstrap the 33 non-trivial DCI integrals for four-loop Coulomb branch amplitudes (providing an independent verification of the four-loop magic identity), which all take remarkably simple form as weight-88 single-valued harmonic polylogarithms. We also compute all leading singularities and a large portion of the pure functions for the 3434 DCI integrals contributing to five-loop amplitudes, where not only some integrals evaluate to functions beyond harmonic polylogarithms but they also contain lower-weight pieces individually.

1 Introduction

Recent years have seen enormous progress in the study of scattering amplitudes and other physical quantities in quantum field theory, especially in the most symmetric of all four-dimensional gauge theories, the 𝒩=4{\cal N}=4 supersymmetric Yang-Mills (SYM) theory in the planar limit. The two complementary sides of the story are the development of powerful new tools for computations, and consequently the discovery of hidden simplicity and mathematical structures (which usually leads to new computational tools etc.). In this story, correlation functions of half-BPS operators (the simplest local operators in the theory) have played a key role and attracted lots of interests because they stand at the crossroads of perturbative calculations, integrability and conformal bootstrap.

There are two seemingly unrelated relations between half-BPS correlators and scattering amplitudes in two different theories (see the nice review Heslop:2022xgp and references therein). At strong coupling they are dual to type-IIB supergravity amplitudes in string theory on AdS×5S5{}_{5}\times S^{5} background via AdS/CFT; at weak coupling, we have the famous duality relating lightlike limits of correlators to null polygonal Wilson loops Alday:2010zy ; Adamo:2011dq and equivalently (the square of) planar scattering amplitudes in 𝒩=4{\cal N}=4 SYM Alday:2007hr ; Drummond:2007aua ; Brandhuber:2007yx ; Bern:2008ap ; Drummond:2008aq ; Eden:2010zz ; Arkani-Hamed:2010zjl ; Mason:2010yk ; Eden:2010ce ; Caron-Huot:2010ryg ; Eden:2011yp ; Eden:2011ku . Furthermore, this “triality” holds at the integrand level, and has provided important insights and rich data for both correlators and amplitudes: the \ell-loop planar integrand of the simplest four-point correlator of the stress-tensor supermultiplet already contains the \ell-loop integrand of four-point amplitudes and (1)(\ell-1)-loop five-point amplitudes etc. Eden:2010zz ; Ambrosio:2013pba ; Heslop:2018zut , in various lightlike limits. In fact the most accurate perturbative data for amplitude integrands are obtained from that of the four-point correlators, and, in the other direction, the relation to squared amplitudes also imposes powerful constraints on planar correlators. Following earlier works Gonzalez-Rey:1998wyj ; Eden:1998hh ; Eden:1999kh ; Eden:2000mv ; Bianchi:2000hn , the planar integrand of four-point correlator was determined for three to seven loops in Eden:2011we ; Eden:2012tu ; Ambrosio:2013pba (the integrated results are known up to three loops Drummond:2013nda , and for results of higher Kaluza-Klein modes and higher points, see Chicherin:2015edu ; Chicherin:2018avq ; Caron-Huot:2021usw ; Bargheer:2022sfd ; Caron-Huot:2023wdh ). The key breakthrough was the discovery of a hidden permutation symmetry for the integrand of the correlator Eden:2011we , which when combined with correlator/amplitude duality has allowed a powerful graphical bootstrap program which determine four-point integrands (for both correlators and amplitudes) up to ten loops Bourjaily:2011hi ; Bourjaily:2015bpz ; Bourjaily:2016evz . Very recently, a new, universal properties has been found for the correlators which greatly improves the graphical bootstrap: this single rule has very efficiently determined correlator (and amplitude) integrands up to eleven He:2024cej and even twelve loops toappear:2 !

Although the planar loop integrand for four-point correlators and amplitudes has been determined to impressively high orders, very little is known about integrated results in general. The one- and two-loop correlators involve conformal integrals which are dual to the well-known box and double-box integrals respectively, but already the three-loop case needs rather non-trivial integrals which were bootstrapped in Drummond:2013nda . These integrals contain new leading singularities and new transcendental functions beyond the simplest (single-valued) harmonic polylogarithms Remiddi:1999ew ; Brown:2004ugm ; Dixon:2012yy ; Schnetz:2013hqa , which, at the symbol level Goncharov:2010jf ; Duhr:2011zq ; Duhr:2012fh means new symbol letters beyond the simplest alphabet. Our original motivation was to similarly bootstrap all (3232) conformal integrals at four loops: while we have succeeded for the majority of these integrals, the remaining few turn out to be very difficult to determine (including two integrals that contain elliptic pieces toappear:1 ). Therefore, except for some families of special cases such as ladders Usyukina:1993ch and fishnet integrals Basso:2017jwq ; Basso:2021omx ; Aprile:2023gnh , it remains an important open question how to systematically compute higher-loop, four-point conformal integrals analytically.

On the other hand, in a remarkable paper Caron-Huot:2021usw , the authors have discovered a hidden ten-dimensional conformal symmetry which not only allows one to package all higher-charge four-point correlators into the 10d correlator (whose integrands can be expressed in terms of ff-graphs), but also provide exact, integrated results for 10d lightlike limit of the correlators known as large R-charge correlator or the octagon Coronado:2018ypq . As explained in Caron-Huot:2021usw , the latter has a nice physical interpretation as the four-point Coulomb branch amplitudes (with non-vanishing v.e.v. for scalar fields), thus via the duality to the octagon these amplitudes (with different R-charge components) have been solved to all orders, which are given by determinants of ladder integrals. As we will review shortly, the interplay between (10d) correlator integrands using ff-graphs, 10d lightlike limit as Coulomb branch amplitudes and the exact results provided by the octagon proves extremely powerful for constraining both the integrand and integrated results for dual conformal invariant (DCI) integrals contributing to Coulomb branch amplitudes.

In the first part of this paper we will summarize some “data” obtained from these constraints. In sec. 2 we will discuss how the 10d limit can be used to determine coefficients of some ff-graphs for the correlator integrand: already the consistency with possible 10d limit implies vanishing of a large class of such coefficients (extending previous zeros from square rules Bourjaily:2016evz ), and we also obtain intriguing connections between some ff-graphs and special DCI integrals such as Basso-Dixon fishnet and generalizations such as next-to fishnet Caron-Huot:2021usw . We have checked all these constraints using both explicit results for ff-graphs up to eleven loops He:2024cej and even higher-loop checks from square rules etc.. Moreover, in sec. 3 we will continue the study of the so-called “magic identities” initiated in Caron-Huot:2021usw ; it is well known that many DCI integrals are related to each other: already the original two-loop “magic identity” for double-box integral Drummond:2006rz implies numerous higher-loop relations, e.g. three-loop “tennis-court” integral equals the ladder integral; by exploiting the fact that certain combinations of DCI integrals equal to various components of Coulomb branch amplitudes which are given by determinants of ladders, the authors of Caron-Huot:2021usw have derived more non-trivial “magic identities” at four and five loops. Note that these are highly non-trivial identities among DCI integrals which cannot be explained from other method such as integration-by-parts identities etc.. We simplify such identities by subtracting known “zeros” which can be understood from lower-loop identities, e.g. they can be these ladder-type identities (originated from double-box case) or those generated by “inverse boxing” of lower-loop ones. In this way, we enumerate such non-trivial, simplified magic identities through six loops, which shows how these seemingly different DCI integrals can be closely related to each other. We will also comment on periods of conformal integrals Broadhurst:1995km ; Schnetz:2008mp ; Brown:2009ta (or the underlying ff-graphs) which can serve as consistency checks of such magic identities.

All these studies have shown that, although it seems too difficult to compute all these conformal integrals for four loops and above, those DCI integrals contributing to Coulomb branch amplitudes may be relatively tamed. Indeed in sec. 4 we will attempt to bootstrap such individual DCI integrals at four and five loops. Recall that for four loops, out of the 3232 conformal integrals, only 88 contribute to Coulomb branch amplitudes where 55 of them all equal to ladder integral (by the original, ladder-type magic identity) Drummond:2006rz . As shown in Caron-Huot:2021usw , the remaining 33 non-trivial DCI integrals (dubbed d2d2, f2f2 and ff) form a new magic identity, and our main target is to bootstrap them individually, which then provide an independent check for the identity. Note that the d2d_{2} integral and f2f_{2} integral have been obtained before: the former was bootstrapped in Drummond:2013nda  111In fact d2d_{2} and f2f_{2} integral can be directly computed by the package HyperlogProcedures Schnetz:2013hqa ; Borinsky:2022lds (but not for ff integral), which proves to be extremely useful for these computations. and the latter is the simplest 2×22\times 2 fishnet integral Basso:2017jwq . We obtain all these integrals analytically in terms of single-valued harmonic polylogarithms (SVHPL) Brown:2004ugm ; Drummond:2012bg ; Schnetz:2013hqa ; Brown:2015ztw : we first analyze their leading singularities and the accompanying function space consists of (parity-odd or even) SVHPL at weight-88; by imposing boundary conditions such as asymptotic expansion of these integrals Chicherin:2018avq , it becomes very straightforward to localize these functions, which all take remarkably simple form in terms of SVHPL functions.

We then move to five loops: there are 3434 DCI integrals but lots of them are either known (e.g. equal to ladder integral) or can be computed using HyperlogProcedures, but we find that some of these integrals are more complicated than expected. They can evaluate to functions that are no longer SVHPL but (single-valued) multiple polylogarithms Chavez:2012kn ; Drummond:2013nda ; Schnetz:2021ebf (in analogy with “hard” integral at three loops or general conformal integrals at four loops). More importantly, some of these integrals contain lower-weight pieces, e.g. one can easily see that some of them actually give weight-99 functions. We compute all leading singularities of these DCI integrals, which then in principle allows us to bootstrap their maximal-transcendental part, and we also present some of the accompanying pure functions (both for maximal, weight-1010 cases and for some weight-99 parts). We end with conclusions and some open equations in sec. 5.

2 From ff-graphs to Coulomb-branch amplitudes/octagons

We consider 𝒩=4\mathcal{N}=4 super Yang-Mills with gauge group SU(Nc)\mathrm{SU}(N_{c}). There are 66 real scalars, ϕI(x)\phi^{I}(x) where I=1,2,,6I=1,2,\ldots,6 in the fundamental representation of the internal symmetry group. By contracting with the co-ordinates yIy_{I} satisfying yIyI=0y_{I}y^{I}=0,

ϕI(x)ϕ(x,y)=yIϕI(x),\phi^{I}(x)\rightarrow\phi(x,y)=y_{I}\phi^{I}(x), (1)

the half BPS operators fall into symmetric traceless representations of SO(6)\mathrm{SO}(6) which are then represented as the products of ϕ(x,y)\phi(x,y). In this paper, we focus on the single trace operator with general charge kik_{i},

𝒪ki(xi)=1ki(24π2Nc)ki2Tr[ϕ(x,y)ki],\mathcal{O}_{k_{i}}(x_{i})=\frac{1}{k_{i}}\left(\frac{2}{4\pi^{2}N_{c}}\right)^{\frac{k_{i}}{2}}\mathrm{Tr}\left[\phi(x,y)^{k_{i}}\right], (2)

where ki2k_{i}\geq 2 since we are considering the SU(Nc)\mathrm{SU}(N_{c}) group. When ki=2k_{i}=2, we come back to the lowest charge half-BPS operator which is also called the stress-tensor multiplet since it contains the stress-tensor and the on-shell Lagrangian.

We are interested in the connected correlator of four such operators in the planar limit, perturbatively we have

Nc2𝒪k1𝒪k2𝒪k3𝒪k4c=Gk1k2k3k4free+=1Gk1k2k3k4()+O(1/Nc2),N_{c}^{2}\langle\mathcal{O}_{k_{1}}\mathcal{O}_{k_{2}}\mathcal{O}_{k_{3}}\mathcal{O}_{k_{4}}\rangle_{\rm c}=G^{\rm free}_{k_{1}k_{2}k_{3}k_{4}}+\sum_{\ell=1}^{\infty}G^{(\ell)}_{k_{1}k_{2}k_{3}k_{4}}+O(1/N_{c}^{2}), (3)

where the loop corrections Gk1k2k3k4()G^{(\ell)}_{k_{1}k_{2}k_{3}k_{4}} can be computed using the Lagrangian insertion method, i.e. the (4+)(4+\ell)-point correlator with four external operators and \ell chiral Lagrangians at Born level Intriligator:1998ig ; Eden:2011we :

Gk1k2k3k4()\displaystyle G^{(\ell)}_{k_{1}k_{2}k_{3}k_{4}} =(g2)!d4x5π2d4x4+nπ2𝒢k1k2k3k4()\displaystyle=\frac{(-g^{2})^{\ell}}{\ell!}\int\frac{d^{4}x_{5}}{\pi^{2}}\cdots\frac{d^{4}x_{4+n}}{\pi^{2}}\mathcal{G}^{(\ell)}_{k_{1}k_{2}k_{3}k_{4}} (4)
=(g2)!d4x5π2d4x4+nπ2𝒪k1𝒪k2𝒪k3𝒪k4(x5)(x4+)(0).\displaystyle=\frac{(-g^{2})^{\ell}}{\ell!}\int\frac{d^{4}x_{5}}{\pi^{2}}\cdots\frac{d^{4}x_{4+n}}{\pi^{2}}\langle\mathcal{O}_{k_{1}}\mathcal{O}_{k_{2}}\mathcal{O}_{k_{3}}\mathcal{O}_{k_{4}}\mathcal{L}(x_{5})\cdots\mathcal{L}(x_{4+\ell})\rangle^{(0)}.

For the four-point correlator, the non-renormalization theorem Eden:2000bk states that the loop integrand 𝒢k1k2k3k4()\mathcal{G}^{(\ell)}_{k_{1}k_{2}k_{3}k_{4}} can factor out a universal prefactor

𝒢k1k2k3k4=R1234(2x122x132x142x232x242x342)k1k2k3k4(),\mathcal{G}_{k_{1}k_{2}k_{3}k_{4}}=R_{1234}\left(2x_{12}^{2}x_{13}^{2}x_{14}^{2}x_{23}^{2}x_{24}^{2}x_{34}^{2}\right)\mathcal{H}^{(\ell)}_{k_{1}k_{2}k_{3}k_{4}}, (5)

where the k1k2k3k4()\mathcal{H}^{(\ell)}_{k_{1}k_{2}k_{3}k_{4}} is the integrand for the reduced correlator and

R1234=d132d242x132x242+d12d23d34d14(x132x242x122x342x142x232)+(12)+(14),R_{1234}=d_{13}^{2}d_{24}^{2}x_{13}^{2}x_{24}^{2}+d_{12}d_{23}d_{34}d_{14}\left(x_{13}^{2}x_{24}^{2}-x_{12}^{2}x_{34}^{2}-x_{14}^{2}x_{23}^{2}\right)+(1\leftrightarrow 2)+(1\leftrightarrow 4), (6)

where dij=yij2xij2d_{ij}=-\frac{y_{ij}^{2}}{x_{ij}^{2}}. Notice that the operator 𝒪ki(xi,yi)\mathcal{O}_{k_{i}}(x_{i},y_{i}) carries kik_{i} charges and mass dimension kik_{i} while the prefactor is degree 22 of both yiy_{i} and xix_{i}, therefore, the reduced integrand carries ki2k_{i}-2 charges and dimension ki+2k_{i}+2.

2.1 Review: four-point correlators and Coulomb branch amplitudes

ff-graphs for correlator integrands

Especially, for the lowest charge (ki=2k_{i}=2) half-BPS operator, the reduced integrand 2222()\mathcal{H}^{(\ell)}_{2222} is just a rational function of xij2x_{ij}^{2} which has at most simple poles as a consequence of OPE Green:2020eyj and has conformal weight 4 with respect to each xix_{i}. More importantly, it is invariant under arbitrary permutation of x1,,x4+x_{1},\cdots,x_{4+\ell} i,e, the full permutation S4+S_{4+\ell} instead of only S4×SS_{4}\times S_{\ell} Eden:2011we . The full permutation symmetry indicates that it is useful to package the \ell-loop integrands into (4+)(4+\ell)-point ff-graphs where the solid line between vertex ii and jj represents a pole 1/xij21/x_{ij}^{2} and the dashed line denotes a numerator xij2x_{ij}^{2}. Each ff-graph (after performing Expand[] in Mathematica) represents a permutation invariant sum of terms, divided by the order of the automorphism group of the ff-graph |aut(f)|\left|\mathrm{aut}(f)\right|, so that the coefficient of any term is equal to 1. For example when 3\ell\leq 3 Eden:2012tu ,

2222(1)=[Uncaptioned image]=11i<j5xij2,\displaystyle\mathcal{H}^{(1)}_{2222}=\includegraphics[align=c,scale={0.3}]{graph/f1.pdf}=\frac{1}{\prod_{1\leq i<j\leq 5}x_{ij}^{2}}, (7)
2222(2)=[Uncaptioned image]=148x122x342x562+S6permutations1i<j6xij2,\displaystyle\mathcal{H}^{(2)}_{2222}=\includegraphics[align=c,scale={0.3}]{graph/f2.pdf}=\frac{1}{48}\frac{x_{12}^{2}x_{34}^{2}x_{56}^{2}+S_{6}\,\text{permutations}}{\prod_{1\leq i<j\leq 6}x_{ij}^{2}},
2222(3)=[Uncaptioned image]=120(x122)2(x342x452x562x672x372)+S7permutations1i<j7xij2.\displaystyle\mathcal{H}^{(3)}_{2222}=\includegraphics[align=c,scale={0.3}]{graph/f3.pdf}=\frac{1}{20}\frac{(x_{12}^{2})^{2}\left(x_{34}^{2}x_{45}^{2}x_{56}^{2}x_{67}^{2}x_{37}^{2}\right)+S_{7}\,\text{permutations}}{\prod_{1\leq i<j\leq 7}x_{ij}^{2}}.

Generally, 2222()\mathcal{H}^{(\ell)}_{2222} can be written as linear combination of all ff-graphs with (4+)(4+\ell)-vertecies 2222()=i=1𝒩ci()fi()\mathcal{H}^{(\ell)}_{2222}=\sum_{i=1}^{\mathcal{N}_{\ell}}c_{i}^{(\ell)}f_{i}^{(\ell)} where 𝒩\mathcal{N}_{\ell} denotes the number of (4+)(4+\ell)-point ff-graphs and ci()c_{i}^{(\ell)} are some rational coefficients. Considering the cusp limit where two consecutive separations x122=x232=0x_{12}^{2}=x_{23}^{2}=0, these coefficients are fixed up to =11\ell=11 He:2024cej and even =12\ell=12 toappear:2 .

According to (4)-(7), we obtain the explicit expression of the correlator, written as the so-called conformal integrals,

G2222(1)\displaystyle G^{(1)}_{2222} =2g2R1234×g1234,\displaystyle=-2g^{2}R_{1234}\times g_{1234},
G2222(2)\displaystyle G^{(2)}_{2222} =2g4R1234(h12;34+h34;12+h14;23+h23;14+h13;24+h24;13\displaystyle=2g^{4}R_{1234}\bigg{(}h_{12;34}+h_{34;12}+h_{14;23}+h_{23;14}+h_{13;24}+h_{24;13}
+12(x122x342+x132x242+x142x232)[g1234]2),\displaystyle\qquad\qquad\left.\qquad+\frac{1}{2}\left(x_{12}^{2}x_{34}^{2}+x_{13}^{2}x_{24}^{2}+x_{14}^{2}x_{23}^{2}\right)[g_{1234}]^{2}\right),

where the one-loop box and two-loop ladder are defined as

g1234\displaystyle g_{1234} =d4x5x152x252x352x452,\displaystyle=\int\frac{d^{4}x_{5}}{x_{15}^{2}x_{25}^{2}x_{35}^{2}x_{45}^{2}},
h13;24\displaystyle h_{13;24} =x242d4x5d4x6(x152x252x452)x562(x262x362x462).\displaystyle=x_{24}^{2}\int\frac{d^{4}x_{5}d^{4}x_{6}}{(x_{15}^{2}x_{25}^{2}x_{45}^{2})x_{56}^{2}(x_{26}^{2}x_{36}^{2}x_{46}^{2})}.

We can also identify the integral to a Feynman diagram in either momentum or position space. For example, the two-loop ladder integral can be drawn as

h13;24[Uncaptioned image],h_{13;24}\Leftrightarrow\includegraphics[align=c]{graph/mompos.pdf}, (8)

where the position space diagram (gray line) is dual to the momentum space diagram (black line). And the conformal integral h13;24h_{13;24} is equivalent to the double-box integral with 4 external legs.

At =3\ell=3, there are 55 conformal integrals from (7).

Refer to caption
Figure 1: 5 conformal integrals at =3\ell=3.
L13;24\displaystyle L_{13;24} =a=57d4xa(x242)2(π2)3d4x5d4x6d4x7(x152x252x452)x562(x262x462)x672(x272x372x472)\displaystyle=\int\prod_{a=5}^{7}d^{4}x_{a}\frac{\left(x_{24}^{2}\right)^{2}}{\left(\pi^{2}\right)^{3}}\int\frac{d^{4}x_{5}d^{4}x_{6}d^{4}x_{7}}{\left(x_{15}^{2}x_{25}^{2}x_{45}^{2}\right)x_{56}^{2}\left(x_{26}^{2}x_{46}^{2}\right)x_{67}^{2}\left(x_{27}^{2}x_{37}^{2}x_{47}^{2}\right)} (9)
T13;24\displaystyle T_{13;24} =a=57d4xax242(π2)3d4x5d4x6d4x7x172(x152x252)(x162x462)(x272x372x472)x562x572x672\displaystyle=\int\prod_{a=5}^{7}d^{4}x_{a}\frac{x_{24}^{2}}{\left(\pi^{2}\right)^{3}}\int\frac{d^{4}x_{5}d^{4}x_{6}d^{4}x_{7}x_{17}^{2}}{\left(x_{15}^{2}x_{25}^{2}\right)\left(x_{16}^{2}x_{46}^{2}\right)\left(x_{27}^{2}x_{37}^{2}x_{47}^{2}\right)x_{56}^{2}x_{57}^{2}x_{67}^{2}}
E12;34\displaystyle E_{12;34} =a=57d4xax232x242(π2)3d4x5d4x6d4x7x162(x152x252x352)x562(x262x362x462)x672(x172x272x472)\displaystyle=\int\prod_{a=5}^{7}d^{4}x_{a}\frac{x_{23}^{2}x_{24}^{2}}{\left(\pi^{2}\right)^{3}}\int\frac{d^{4}x_{5}d^{4}x_{6}d^{4}x_{7}x_{16}^{2}}{\left(x_{15}^{2}x_{25}^{2}x_{35}^{2}\right)x_{56}^{2}\left(x_{26}^{2}x_{36}^{2}x_{46}^{2}\right)x_{67}^{2}\left(x_{17}^{2}x_{27}^{2}x_{47}^{2}\right)}
H12;34\displaystyle H_{12;34} =a=57d4xax142x232x342(π2)3d4x54d6d4x7x572(x152x252x352x452)x562(x362x462)x672(x172x272x372x472)\displaystyle=\int\prod_{a=5}^{7}d^{4}x_{a}\frac{x_{14}^{2}x_{23}^{2}x_{34}^{2}}{\left(\pi^{2}\right)^{3}}\int\frac{d^{4}x_{5}^{4}d_{6}d^{4}x_{7}x_{57}^{2}}{\left(x_{15}^{2}x_{25}^{2}x_{35}^{2}x_{45}^{2}\right)x_{56}^{2}\left(x_{36}^{2}x_{46}^{2}\right)x_{67}^{2}\left(x_{17}^{2}x_{27}^{2}x_{37}^{2}x_{47}^{2}\right)}
(gh)13;24\displaystyle(gh)_{13;24} =a=57d4xax132(x242)2(π2)3d4x5x152x252x352x452d4x6d4x7(x162x262x462)x672(x272x372x472)\displaystyle=\int\prod_{a=5}^{7}d^{4}x_{a}\frac{x_{13}^{2}\left(x_{24}^{2}\right)^{2}}{\left(\pi^{2}\right)^{3}}\int\frac{d^{4}x_{5}}{x_{15}^{2}x_{25}^{2}x_{35}^{2}x_{45}^{2}}\int\frac{d^{4}x_{6}d^{4}x_{7}}{\left(x_{16}^{2}x_{26}^{2}x_{46}^{2}\right)x_{67}^{2}\left(x_{27}^{2}x_{37}^{2}x_{47}^{2}\right)}

We remark that at =1,2,,5\ell=1,2,\cdots,5, there are 1,2,5,32,2551,2,5,32,255 conformal integrals, but the number of DCI integrals contributing to Coloumb branch amplitudes are smaller: already at two and three loops only 11 (double-box) and 22 (ladder and tennis-court) contribute, while at four and five loops, only 88 and 3434 DCI integrals contribute respectively. These integrals are defined in Appendix A. There are 229,1873,66928343229,1873,\cdots 66928343 such DCI integrals (with non-zero coefficients) contribute to Coulomb branch amplitudes at =6,7,,11\ell=6,7,\cdots,11 Bourjaily:2016evz ; He:2024cej .

For higher charge correlators, the reduced integrand can be written as a sum of different RR-charge structures

k1k2k3k4()=bijki=2+jbij{bij}()×1i<j4(dij)bij.\mathcal{H}^{(\ell)}_{k_{1}k_{2}k_{3}k_{4}}=\sum_{\begin{subarray}{c}b_{ij}\\ k_{i}=2+\sum_{j}b_{ij}\end{subarray}}\mathcal{F}^{(\ell)}_{\{b_{ij}\}}\times\prod_{1\leq i<j\leq 4}\left(d_{ij}\right)^{b_{ij}}. (10)

For example, 2222()=F{0,0,0,0,0,0}()\mathcal{H}_{2222}^{(\ell)}=F_{\{0,0,0,0,0,0\}}^{(\ell)}. The higher-charge integrand was bootstraped up to =5\ell=5 in Chicherin:2015edu ; Chicherin:2018avq . Surprisingly, the integrals that appear for higher-charge correlators are exactly those of the stress-tensor multiplet only different by their coefficients. In fact, Caron-Huot:2021usw indicates that all correlators with different R-charges come from a single master four-point correlator 𝖦(xij2,yij2)\mathsf{G}(x_{ij}^{2},y_{ij}^{2}) of the master operator

𝖦(xij2,yij2)=𝖮(x1,y1)𝖮(x2,y2)𝖮(x3,y3)𝖮(x4,y4),with𝖮(x,y)k=2Ok(x,y).\mathsf{G}(x_{ij}^{2},y_{ij}^{2})=\langle\mathsf{O}(x_{1},y_{1})\mathsf{O}(x_{2},y_{2})\mathsf{O}(x_{3},y_{3})\mathsf{O}(x_{4},y_{4})\rangle,\quad\text{with}\;\mathsf{O}(x,y)\equiv\sum_{k=2}^{\infty}\mathrm{O}_{k}(x,y). (11)

Similar to (5), the loop integrand of the master correlator can pull out a universal factor

𝖦()=(g2)!R1234(2x122x132x142x232x242x342)dx54π2dx4+n4π2𝖧()(xij2,yij2),\mathsf{G}^{(\ell)}=\frac{(-g^{2})^{\ell}}{\ell!}R_{1234}\left(2x_{12}^{2}x_{13}^{2}x_{14}^{2}x_{23}^{2}x_{24}^{2}x_{34}^{2}\right)\int\frac{dx_{5}^{4}}{\pi^{2}}\cdots\frac{dx_{4+n}^{4}}{\pi^{2}}\mathsf{H}^{(\ell)}(x_{ij}^{2},y_{ij}^{2}), (12)

where yi=0y_{i}=0 for i5i\geq 5. As indicated in Caron-Huot:2021usw , 𝖧()(xij2,yij2)\mathsf{H}^{(\ell)}(x_{ij}^{2},y_{ij}^{2}) is a ten-dimensional conformal invariant object, i.e. a function of the distance in D=10D=10 Xij2xij2+yij2=(1dij)xij2X_{ij}^{2}\equiv x_{ij}^{2}+y_{ij}^{2}=(1-d_{ij})x_{ij}^{2}, and also possesses the full permutation symmetry S4+S_{4+\ell}. A given correlator can be obtained by extracting the coefficient of the correct power of yiy_{i},

k1k2k3k4()(xij2,yij2)=coefficient of (i=14βiki2) in 𝖧()(Xij2)|yij2βiβjyij2,\mathcal{H}^{(\ell)}_{k_{1}k_{2}k_{3}k_{4}}(x_{ij}^{2},y_{ij}^{2})=\mbox{coefficient of $\left(\prod_{i=1}^{4}\beta_{i}^{k_{i}-2}\right)$ in }\mathsf{H}^{(\ell)}(X_{ij}^{2})\big{|}_{y_{ij}^{2}\to\beta_{i}\beta_{j}y^{2}_{ij}}, (13)

For example, at =1\ell=1

𝖧(1)=11i<j5xij2=11i<j5xij2(1dij)=1ξ4g1234{bij}1i<j4(dij)bij,\mathsf{H}^{(1)}=\frac{1}{\prod_{1\leq i<j\leq 5}x_{ij}^{2}}=\frac{1}{\prod_{1\leq i<j\leq 5}x_{ij}^{2}(1-d_{ij})}=\frac{1}{\xi_{4}}g_{1234}\sum_{\{b_{ij}\}}\prod_{1\leq i<j\leq 4}\left(d_{ij}\right)^{b_{ij}}, (14)

where ξ4=x122x232x342x142x134x244\xi_{4}=x_{12}^{2}x_{23}^{2}x_{34}^{2}x_{14}^{2}x_{13}^{4}x_{24}^{4}. Therefore, all component correlators are equal and given by the one-loop box function (with unit coefficient).

From (13), the 10d correlator can be obtained by uplifting the lowest charge correlator in D=4D=4. The only possible issue is that certain linear combination of conformal integrals might vanish at D=4D=4 due to gram identity, which could otherwise contribute to the 10d correlator. However, as verified in Caron-Huot:2021usw up to 7 loops, no nontrivial gram identities exist.

Coulomb-branch amplitudes from 10d lightlike limit

One of the most important properties of the simplest half-BPS correlator is the amplitude/correlator duality Eden:2010ce . To be more explicit, the amplitude/correlator duality is realized in two steps. First, the light-like limit of the correlator is proportional to a light-like polygon Wilson loop in the adjoint representation Alday:2010zy ; Adamo:2011dq where the later one equals the square of the Wilson loop in the fundamental representation in the planar limit. Secondly, in the planar limit, the scattering amplitude is dual to the polygonal Wilson loop in the fundamental representation Alday:2007hr ; Drummond:2007aua ; Brandhuber:2007yx ; Bern:2008ap ; Drummond:2008aq ; Eden:2010zz ; Arkani-Hamed:2010zjl ; Mason:2010yk ; Caron-Huot:2010ryg ; Eden:2011yp ; Eden:2011ku . At 44-point

limxi,i+120𝒢2222(xi)𝒢2222free=A(xi)2,\lim_{x_{i,i+1}^{2}\to 0}\frac{\mathcal{G}_{2222}(x_{i})}{\mathcal{G}_{2222}^{\rm free}}=A(x_{i})^{2}, (15)

where A(xi)A(x_{i}) is the planar four-point amplitude, divided by the tree-level amplitude. Equivalently, for the \ell-loop correction

limxi,i+1202ξ42222()=(A(xi)2)(),\lim_{x_{i,i+1}^{2}\to 0}2\xi_{4}\mathcal{H}_{2222}^{(\ell)}=(A(x_{i})^{2})^{(\ell)}, (16)

where ξ4\xi_{4} is defined as earlier and (A(xi)2)()==0()A(xi)()A(xi)()(A(x_{i})^{2})^{(\ell)}=\sum_{\ell^{\prime}=0}^{\ell}\binom{\ell}{\ell^{\prime}}A(x_{i})^{(\ell^{\prime})}A(x_{i})^{(\ell-\ell^{\prime})}. Graphically, the \ell-loop amplitudes A(xi)()=A(xi)()×A(xi)(0)A(x_{i})^{(\ell)}=A(x_{i})^{(\ell)}\times A(x_{i})^{(0)} term can be obtained from the ff-graph by taking inequivalent four-face, i.e. four-cycles which encompass no internal points. The loop integrands are the so-called dual conformal invariant (DCI) integrands.

A similar question can be asked about the 10d light-like limit of the master correlator (11), i.e. Xi,i+120di,i+11yij2xij2X^{2}_{i,i+1}\to 0\Leftrightarrow d_{i,i+1}\to 1\Leftrightarrow y_{ij}^{2}\to-x_{ij}^{2}. In fact, as shown in Caron-Huot:2021usw , the light-like limit in D=10D=10 corresponds to the large R charge limit of the master correlator. Under the large charge limit, the master correlator becomes the square of a certain octagon form factor 𝕆\mathbb{O} Coronado:2018ypq obtained by gluing two hexagons Basso:2015zoa together

Resdi,i+1=1G|integrand=𝕆×𝕆|integrand,\mathrm{Res}_{d_{i,i+1}=1}G\big{|}_{\rm integrand}=\mathbb{O}\times\mathbb{O}\big{|}_{\rm integrand}, (17)

which is similar to the correlator/Wilson loop duality. From the amplitude side, a natural extension to D=10D=10 is the scattering amplitude on the Coulomb branch since the nonzero yi,i4y_{i},i\leq 4 coordinates can be interpreted as vacuum expectation values for the scalar fields of the 𝒩=4\mathcal{N}=4 SYM theory

𝕆(x,y)𝕆free(x,y)=M(x,y),\frac{\mathbb{O}(x,y)}{\mathbb{O}^{\text{free}}(x,y)}=M(x,y), (18)

where MM is a massive four-particle amplitude in the Coulomb branch222this is a special case of Coulomb branch amplitude where yi2=0y_{i}^{2}=0. normalized so that Mfree=1M^{\text{free}}=1, with external momenta and masses are piμxi,i+1μ,mi2yi,i+12p_{i}^{\mu}\equiv x^{\mu}_{i,i+1},m_{i}^{2}\equiv y_{i,i+1}^{2}. Thus, the light-like condition Xi,i+12=0X_{i,i+1}^{2}=0 becomes the on-shell condition pi2+mi2=0p_{i}^{2}+m_{i}^{2}=0. In addition, the generic external mass mi20m_{i}^{2}\neq 0 also keeps the integrals off-shell, meaning there are no divergences of the integrals. Therefore, (17) and (18) hold both at the integrand and integrated level.

Note that the octagon is known for arbitrary values of the coupling from integrability Belitsky:2020qrm ; Belitsky:2020qir , and the Coulomb branch amplitude can be computed from the master correlator perturbatively, (18) produces some all-loop predictions for the correlators at the integrand level as well as some new results for integrals. Before that, we first discuss the DCI integrals from light-like limit on the right hand side of (18) at weak coupling.

DCI integrals from light-like limit and ladder-type magic identity

The Coulomb branch amplitude can be obtained by taking the D=10D=10 light-like limit of the master correlator. For example, at 3\ell\leq 3.

M\displaystyle M =1g2X132X242g1234\displaystyle=1-g^{2}{X_{13}^{2}X_{24}^{2}}g_{1234}
+g4[X132(X242)2x242h13;24+(X132)2X242x132h24;13]\displaystyle\qquad+g^{4}\left[\frac{{X_{13}^{2}(X_{24}^{2})^{2}}}{x_{24}^{2}}h_{13;24}+\frac{{(X_{13}^{2})^{2}X_{24}^{2}}}{x_{13}^{2}}h_{24;13}\right]
g6[X132(X242)3(x242)2L13;24+(X132)3X242(x132)2L24;13+X132(X242)2x242(T13;24+T13;42)\displaystyle\qquad-g^{6}\left[\frac{{X_{13}^{2}(X_{24}^{2})^{3}}}{(x_{24}^{2})^{2}}L_{13;24}+\frac{{(X_{13}^{2})^{3}X_{24}^{2}}}{(x_{13}^{2})^{2}}L_{24;13}+\frac{{X_{13}^{2}(X_{24}^{2})^{2}}}{x_{24}^{2}}(T_{13;24}+T_{13;42})\right.
+(X132)2X242x132(T24;13+T24;31)]+𝒪(g8)\displaystyle\qquad\qquad\quad\left.+\frac{{(X_{13}^{2})^{2}X_{24}^{2}}}{x_{13}^{2}}(T_{24;13}+T_{24;31})\right]+\mathcal{O}(g^{8})

We can see that Ma,bM_{a,b} can be naturally grouped by the power of X132(1d13)X_{13}^{2}\sim(1-d_{13}) and X242(1d24)X_{24}^{2}\sim(1-d_{24})

M=1+a,b1(1d13)a(1d24)bMa,b,M=1+\sum_{a,b\geq 1}(1-d_{13})^{a}(1-d_{24})^{b}M_{a,b}, (19)

where Ma,b=Mb,aM_{a,b}=M_{b,a} at the integrated level. From the ff-graph perspective, the \ell-loop integrand Ma,b()M_{a,b}^{(\ell)} can be obtained by taking the four-faces with specific (a,b)(a,b)-structure in Fig. 2.

Refer to caption
Figure 2: Any ff-graph with (a,b)(a,b)-structure (on the 44-face labeled by 1,2,3,41,2,3,4) contributes to Ma,bM_{a,b}.

When a=1a=1 or b=1b=1, the dashed line becomes the solid line. By taking light-like limit of certain ff-graph, we focus on the rational functions corresponding to Fig. 2 and denote the remaining part as f~\tilde{f}

fx122x232x342x142x134x244(x132)a2(x242)b2x122x232x342x142x134x244f~(1d13)a(1d24)bMa,b().f\rightarrow x_{12}^{2}x_{23}^{2}x_{34}^{2}x_{14}^{2}x_{13}^{4}x_{24}^{4}\frac{(x_{13}^{2})^{a-2}(x_{24}^{2})^{b-2}}{x_{12}^{2}x_{23}^{2}x_{34}^{2}x_{14}^{2}x_{13}^{4}x_{24}^{4}}\tilde{f}\rightarrow(1-d_{13})^{a}(1-d_{24})^{b}M_{a,b}^{(\ell)}. (20)

Therefore, an ff-graph contributes to Ma,bM_{a,b} if and only if it has the (a,b)(a,b)-structure.

Notice that when a=b=1a=b=1, the (1,1)(1,1)-structure is just the k=4k=4 divergence structure in Bourjaily:2015bpz , which means that no ff-graphs can contribute to M1,1M_{1,1} except for the non-planar ff-graph at =1\ell=1, which gives the box integral. When a=1,b=2a=1,b=2, this is the double triangle structure discussed in He:2024cej , which means that almost all ff -graphs contribute to M1,2M_{1,2}. In fact, one can check that except for the 4 ff-graphs without double-triangle structure that occur at =8,11,12\ell=8,11,12, all ff-graphs contribute to M1,2M_{1,2}. The other extremal case is when a,ba,b is large. We will discuss about this case in the next subsection.

At the integrated level, some of these DCI integrals are related by the ladder-type magic identity Drummond:2006rz . This kind of magic identity stems from the fact that the ss-channel double box at =2\ell=2 equals the tt-channel one.

[Uncaptioned image]=[Uncaptioned image]\includegraphics[scale={0.5},align=c]{graph/ladders}=\includegraphics[scale={0.5},align=c]{graph/laddert} (21)

This equality tells us has all loop generalization by adding box to the momentum space Feynman diagram (or the so-called sling-shot rule in position space). For example at =3\ell=3, it tells us the ladder integral equals to the tennis-court integral.

[Uncaptioned image]=[Uncaptioned image]\includegraphics[scale={0.5},align=c]{graph/ladder3}=\includegraphics[scale={0.5},align=c]{graph/tsc} (22)

We can also translate the “adding box” procedure to position space, which turn out to be the “slingshot” rule Drummond:2006rz 333Slingshot rule can only generate part of the ladder-type identities. We need to consider relation like ~()=~()~()\tilde{\mathcal{I}}^{(\ell)}=\tilde{\mathcal{I}}^{(\ell^{\prime})}\tilde{\mathcal{I}}^{(\ell-\ell^{\prime})}, which has slingshot rule as its special case when =1\ell^{\prime}=1.. To be explicit, we start from the magic identity 1()=2()==m(){\mathcal{I}}^{(\ell)}_{1}={\mathcal{I}}^{(\ell)}_{2}=\ldots={\mathcal{I}}_{m_{\ell}}^{(\ell)} at \ell-loop whose external points are 1,2,3,41,2,3,4 in order. We pick up one of the external points, for example, the point 1, then multiply i(){\mathcal{I}}^{(\ell)}_{i} by the slingshot factor x242x122x142x+5,12\frac{x_{24}^{2}}{x_{12}^{2}x_{14}^{2}x_{\ell+5,1}^{2}}

j(+1)=[Uncaptioned image]×i(),i=1,2,,m,{\mathcal{I}}^{(\ell+1)}_{j}=\includegraphics[scale={1.2},align=c]{graph/shot.pdf}\times{\mathcal{I}}^{(\ell)}_{i},\quad i=1,2,\ldots,m_{\ell}, (23)

where the vertex 11 becomes an internal point of j(+1){\mathcal{I}}_{j}^{(\ell+1)} and the external points are now +5,2,3,4\ell+5,2,3,4. All the j(+1){\mathcal{I}}_{j}^{(\ell+1)} are eqaul at the integrated level. At the ff-graph level, we can see that the “slingshot” rule can be translated to the inverse square rule we discuss later.

2.2 Fishnet and next-to-fishnet from ff-graphs

We first discuss some all-loop structures from (18) at the integrand level. From integrability, (18) tells us that the leading order of Ma,bM_{a,b} starts from (g2)ab(g^{2})^{ab}. However, it is not natural from the ff-graph perspective. A simple combinatorial argument only shows that an upper bound for (a,b)(a,b)-structure in a (4+)(4+\ell)-point ff-graph is a+b+min(a,b)+2a+b+\mathrm{min}(a,b)\leq\ell+2. For example, the 66-loop ff-graph in Fig. 3 has a (2,4)(2,4)-structure, and its coefficient vanishes from square rule Bourjaily:2016evz .

Refer to caption
Figure 3: A 66-loop ff-graph with (2,4)(2,4) structure; its coefficient vanishes due to the k=4k=4 divergence structure.

Here we conjecture that the coefficient of a (4+)(4+\ell)-point ff-graph with (a,b),ab>(a,b),ab>\ell structure must vanish.444In fact, Ma,b(>ab)M_{a,b}^{(\ell>ab)} vanishes exactly at the integrated level, and we found that they vanish at the integrand level up to =11\ell=11. We list the number of ff-graphs that will naively contribute to Ma,b,ab>M_{a,b},ab>\ell up to =11\ell=11 in Table 1. Up to =10\ell=10, all zeros can be understood by the square rule, however, at =11\ell=11, there are 88 vanishing ci(11)c_{i}^{(11)}’s that the square rule can not explain.

\ell (a,b)(a,b) (2,4) (2,5) (2,6) (2,7) (2,8) (2,9) (3,3) (3,4) (3,5) (3,6) (3,7)
6 1
7 4 1
8 9 1 4 2
9 120 9 1 1 2
10 174 12 1 738 59 2
11 4049 212 12 1 21764 1845 58 2
Table 1: Numbers of \ell-loop ff-graphs potentially /contributing to Ma,bM_{a,b} where ab>ab>\ell.

At the leading order Ma,b(=ab)M_{a,b}^{(\ell=ab)}, the integrability tells us Ma,b(ab)M_{a,b}^{(ab)} only has the contribution from the fishnet diagram Basso:2017jwq .

Mab(ab)=(1)a(ba)×[Uncaptioned image]M_{ab}^{(ab)}=(-1)^{a(b-a)}\times\includegraphics[align=c]{graph/fishnet.pdf} (24)

However, from the ff-graph perspective, there are more than one ff-graphs contributing to Ma,b(=ab)M_{a,b}^{(\ell=ab)}, while only one of them does not vanish. For example, at =6\ell=6, there are 3 ff-graphs with (2,3)(2,3)-structure, but only the middle one in Fig. 4 has a nonzero coefficient.

Refer to caption
Refer to caption
Refer to caption
Figure 4: At =6\ell=6, there are 3 ff-graphs with (2,3)(2,3) structures: only the middle one has nonzero coefficient, while the other two vanish due to the square rule.

For M2,4(8),M3,3(9),M2,5(10)M_{2,4}^{(8)},M_{3,3}^{(9)},M_{2,5}^{(10)}, there are 60,135,242060,135,2420 ff-graphs contributing to it respectively, but only the one producing the fishnet diagram is non-vanishing. Although up to =10\ell=10, the square rule accounts for all the zeros from this fishnet constraint, there might be some zeros that the square rule can not explain at higher loop. And the coefficient of the unique ff-graph contributing to the fishnet is (1)a(ba)(-1)^{a(b-a)}.

At the subleading order Ma,b(=ab+1)M_{a,b}^{(\ell=ab+1)}, as indicated in Caron-Huot:2021usw , the DCI integrand contributing to Ma,b(=ab+1)M_{a,b}^{(\ell=ab+1)} can be obtained by deforming one of the 2×22\times 2 sublattice of the regular fishnet. From the ff-graph perspective, all (5+ab)(5+ab)-point ff-graphs with (a,b)(a,b)-structures can contribute, but only those with nonzero coefficients will really contribute. The non-vanishing ff-graphs contributing to the subleading order can be obtained by two rules, i.e. the inverse square rule and the squeeze rule. The coefficient of other (ab+1)(ab+1)-loop ff-graphs with (a,b)(a,b) structures vanish.

The inverse square rule is related to the square rule; it maps a \ell-loop ff-graph to a (+1)(\ell+1)-loop one. To be explicit, starting from the =ab\ell=ab-loop ff-graph contributing to the fishnet diagram, we find a 44-cycle labeled as 1,2,3,41,2,3,4 in it which does not correspond to the (a,b)(a,b)-structure, then multiply the one-loop box factor x132x242x1,+52x2,+52x3,+52x4,+52\frac{x_{13}^{2}x_{24}^{2}}{x_{1,\ell+5}^{2}x_{2,\ell+5}^{2}x_{3,\ell+5}^{2}x_{4,\ell+5}^{2}} to obtain the new (5+)(5+\ell)-point ff-graph.

[Uncaptioned image]×[Uncaptioned image][Uncaptioned image].\includegraphics[align=c,scale={2}]{graph/box.pdf}\times\includegraphics[align=c,scale={2}]{graph/isf.pdf}\longrightarrow\includegraphics[align=c,scale={2}]{graph/boxresult.pdf}. (25)

The coefficient of the new graphs equal to the original one. For example, from the 44-loop ff-graph contributing to the 2×22\times 2 fishnet, the inverse square rule gives the following two ff-graphs contributing to M2,2(5)M_{2,2}^{(5)}.

[Uncaptioned image][Uncaptioned image]\displaystyle\includegraphics[scale={0.4},align=c]{graph/M224h1.pdf}\Rightarrow\includegraphics[scale={0.4},align=c]{graph/M2251.pdf} (26)
[Uncaptioned image][Uncaptioned image]\displaystyle\includegraphics[scale={0.4},align=c]{graph/M224h2.pdf}\Rightarrow\includegraphics[scale={0.4},align=c]{graph/M2252.pdf}

Meanwhile, the squeeze rule also starts from a 44-cycle which does not correspond to the (a,b)(a,b)-structure of the ff-graph, but multiply a different factor x122x342x1,+52x2,+52x3,+52x4,+52\frac{x_{12}^{2}x_{34}^{2}}{x_{1,\ell+5}^{2}x_{2,\ell+5}^{2}x_{3,\ell+5}^{2}x_{4,\ell+5}^{2}}.

[Uncaptioned image]×[Uncaptioned image][Uncaptioned image].\includegraphics[align=c,scale={2}]{graph/box.pdf}\times\includegraphics[align=c,scale={2}]{graph/sqf.pdf}\longrightarrow\includegraphics[align=c,scale={2}]{graph/sqresult.pdf}. (27)

However, the squeeze rule can result in non-planar ff-graphs, and we only pick up the planar ones. From the 44-loop ff-graph contributing to the 2×22\times 2 fishnet, the squeeze rule gives the following ff-graph which also contribute to M2,2(5)M_{2,2}^{(5)}.

[Uncaptioned image][Uncaptioned image],\includegraphics[scale={0.4},align=c]{graph/M224h2.pdf}\Rightarrow\includegraphics[scale={0.4},align=c]{graph/M2253.pdf}, (28)

One can check that when applied to the ff-graph contributing to the fishnet integral, the coefficient of the ff-graph obtained by squeezing equals the minus of the original one, which can be checked up to =11\ell=11.

3 Magic identities for DCI integrals

In this section, we study magic identities for DCI integrals by equating their combinations (Ma,bM_{a,b}) with octagon results as determinants of ladders. We will see that once we obtain a magic identity at a given loop, it implies an infinite family of higher-loop identities, which can in turn be used to simplify new magic identities. We enumerate such simplified magic identities through six loops.

3.1 Review: from determinants of ladders to magic identities

Now we move onto the Octagons at week coupling, which appears on the left hand side of (18). The hexagonalization Fleury:2016ykk ; Eden:2016xvg procedure allows us to uplift the free correlator to the correlator at the finite coupling. In this procedure, the octagon can be decomposed into two hexagons by inserting a complete set of states along a mirror cut. This cut can be considered as stretching between the hexagons connected by the bridge of length ll.

𝕆=𝕆0+l=1(d13)l𝕆l+(d24)l𝕆l,\mathbb{O}=\mathbb{O}_{0}+\sum_{l=1}^{\infty}(d_{13})^{l}\mathbb{O}_{l}+(d_{24})^{l}\mathbb{O}_{l}, (29)

where the 𝕆l\mathbb{O}_{l} can be expressed as a sum over the number of particles nn Coronado:2018ypq

𝕆l(z,z¯,d13,d24)=1+n=1(1d13d24)nIn,l(z,z¯).\mathbb{O}_{l}(z,\bar{z},d_{13},d_{24})=1+\sum_{n=1}^{\infty}(1-d_{13}d_{24})^{n}I_{n,l}(z,\bar{z}). (30)

where the z,z¯z,\bar{z} is defined via two cross ratio u,vu,v

ux122x342x132x242zz¯,vx142x232x132x242(1z)(1z¯).u\equiv\frac{x_{12}^{2}x_{34}^{2}}{x_{13}^{2}x_{24}^{2}}\equiv z\bar{z},\,v\equiv\frac{x_{14}^{2}x_{23}^{2}}{x_{13}^{2}x_{24}^{2}}\equiv(1-z)(1-\bar{z}). (31)

At weak coupling, the nn-particle contribution In,lI_{n,l} is related to the determinant of ladder integral as a generalization of the Basso-Dixon fishnet Basso:2017jwq . Recall the integrated result of pp-loop ladder integral is Usyukina:1993ch

Fp(z,z¯)=(1)j=p2pj![log(zz1z¯z¯1)]2pjp!(jp)!(2pj)![Lij(zz1)Lij(z¯z¯1)zz¯].F_{p}(z,\bar{z})=(-1)\sum_{j=p}^{2p}\frac{j!\left[-\log(\frac{z}{z-1}\frac{\bar{z}}{\bar{z}-1})\right]^{2p-j}}{p!(j-p)!(2p-j)!}\left[\frac{\text{Li}_{j}(\frac{z}{z-1})-\text{Li}_{j}(\frac{\bar{z}}{\bar{z}-1})}{z-\bar{z}}\right]. (32)

The determinant of ladder integrals is defined as

Fi1,i2,,in=m=1n1im!(im1)!|fi1fi21finn+1fi1+1fi2finn+2fi1+n1fi2+n2fin|with fp=p!(p1)!Fp.F_{i_{1},i_{2},\cdots,i_{n}}=\prod_{m=1}^{n}\frac{1}{i_{m}!(i_{m}-1)!}\begin{vmatrix}f_{i_{1}}&f_{i_{2}-1}&\cdots&f_{i_{n}-n+1}\\ f_{i_{1}+1}&f_{i_{2}}&\cdots&f_{i_{n}-n+2}\\ \vdots&\vdots&\ddots&\vdots\\ f_{i_{1}+n-1}&f_{i_{2}+n-2}&\dots&f_{i_{n}}\\ \end{vmatrix}\qquad\text{with }f_{p}=p!(p-1)!F_{p}. (33)

In,lI_{n,l} in (30) can be written as linear combination of fishnet integrals perturbatively

In,l==l(n+l)(g2)In,l()==l(n+l)(g2)(1)nl×i1++in=with (ip+1ip)2 and i1>lcl;{i}nFi1,i2,,inI_{n,l}=\sum_{\ell=l(n+l)}^{\infty}(-g^{2})^{\ell}I_{n,l}^{(\ell)}=\sum_{\ell=l(n+l)}^{\infty}(-g^{2})^{\ell}(-1)^{nl}\times\sum_{\begin{subarray}{c}i_{1}+\cdots+i_{n}=\ell\\ \text{with }(i_{p+1}-i_{p})\geq 2\text{ and }i_{1}>l\end{subarray}}c_{l;\{i\}_{n}}F_{i_{1},i_{2},\cdots,i_{n}} (34)

where cl;{i}c_{l;\{i\}} are some coefficients. Although these coefficients do not have a closed form generally, the leading and subleading contribution of In,lI_{n,l} is conjectured in Caron-Huot:2021usw

In,l(n(n+l))=(1)nlF1+l,3+l,,2n1+l,\displaystyle I^{(n(n+{l}))}_{n,{l}}=(-1)^{nl}F_{1+{l},3+{l},\cdots,2n-1+{l}}, (35)
In,l(n(n+l)+1)=(1)nl2(2n1+l)F1+l,3+l,,2n3+l,2n1+l+1,\displaystyle I^{(n(n+{l})+{1})}_{n,{l}}=(-1)^{nl}2(2n-1+l)F_{1+{l},3+{l},\cdots,2n-3+l,2n-1+l+1},

where the F1+l,3+l,,2n1+lF_{1+{l},3+{l},\cdots,2n-1+{l}} is just the (n+l)×n(n+l)\times n fishnet integral. And for later use, I2,0(6)=28F1,5+185F2,4I_{2,0}^{(6)}=28F_{1,5}+\frac{18}{5}F_{2,4}.

Notice that 𝕆free=1+l=1(d13)l+l=1(d24)l\mathbb{O}^{\text{free}}=1+\sum_{l=1}^{\infty}(d_{13})^{l}+\sum_{l=1}^{\infty}(d_{24})^{l}, we can collect the dependence of d13,d24d_{13},d_{24} of 𝕆𝕆free\frac{\mathbb{O}}{\mathbb{O}^{\text{free}}} as

𝕆𝕆free=1+a,b1(1d13)a(1d24)bn=min(a,b)a+b1l0l+na,bCa,bn,lIn,l,\frac{\mathbb{O}}{\mathbb{O}^{\text{free}}}=1+\sum_{a,b\geq 1}(1-d_{13})^{a}(1-d_{24})^{b}\sum_{n=\min(a,b)}^{a+b-1}\sum_{\begin{subarray}{c}l\geq 0\\ l+n\geq a,b\end{subarray}}^{\infty}C^{n,l}_{a,b}I_{n,l}, (36)

where

Ca,bn,l=(1)a+bn11+δl,0[(n1a1)(a+l1a+bn1)+(ab)].C^{n,l}_{a,b}=\frac{(-1)^{a+b-n-1}}{1+\delta_{l,0}}\left[\binom{n-1}{a-1}\binom{a+l-1}{a+b-n-1}+(a\leftrightarrow b)\right]. (37)

We list the result of independent Ma,bM_{a,b}’s up to =6\ell=6 in Table 2.

[Uncaptioned image]
Table 2: Result of Ma,b()M_{a,b}^{(\ell)} up to =6\ell=6. The colored highlights and arrows indicate a diagonal structure, which is explained in sec. 3.2. The yellow and blue highlight means the results can be obtained by ladder-type magic identities. The arrows mean the “inverse-boxing” of lower loop magic identity can be used to simplify the higher loop ones.

Formally, the Ma,b()M_{a,b}^{(\ell)} can be written as

Ma,b()=i=1𝒩ci()a,bfi(),M_{a,b}^{(\ell)}=\sum_{i=1}^{\mathcal{N}_{\ell}}c_{i}^{(\ell)}\mathbb{P}_{a,b}f_{i}^{(\ell)}, (38)

where the operator a,b\mathbb{P}_{a,b} is defined as finding all (a,b)(a,b)-structures in a given ff-graph. By combining (38) and (36), we obtain the following identity

i=1𝒩ci()a,bfi()=n=min(a,b)a+b1l0l+na,bCa,bn,lIn,l().\sum_{i=1}^{\mathcal{N}_{\ell}}c_{i}^{(\ell)}\mathbb{P}_{a,b}f_{i}^{(\ell)}=\sum_{n=\min(a,b)}^{a+b-1}\sum_{\begin{subarray}{c}l\geq 0\\ l+n\geq a,b\end{subarray}}^{\infty}C^{n,l}_{a,b}I_{n,l}^{(\ell)}. (39)

The above equation not only gives us some all-loop predictions about coefficients of ff-graphs, but also teaches us new magic identities among multi-loop conformal integrals (listed through six loops).

3.2 Magic identities through six loops

In this subsection, we derive new relations between DCI integrals from (39) and give some results to all loops.

For M1,()M_{1,\ell}^{(\ell)}, it is not hard to see only one ff-graph contributes, since it’s the only graph containing a (1,)(1,\ell)-structure. Thus, M1,()=FM_{1,\ell}^{(\ell)}=F_{\ell} for any \ell. (The yellow highlight in Table. 2.)

M1,()=[Uncaptioned image]=F.M_{1,\ell}^{(\ell)}=\includegraphics[scale={0.4},align=c]{graph/M1L.pdf}=F_{\ell}. (40)

Then for M1,1()M_{1,\ell-1}^{(\ell)}, the contributing integrals are obtained by adding a “rung” to an (l1)(l-1)-loop ladder, which all equal to the ll-loop ladder integral through ladder-type magic identity. This means M1,1()=(1)FM_{1,\ell-1}^{(\ell)}=(\ell-1)F_{\ell}. (The blue highlight in Table. 2.) Take =3\ell=3 and =4\ell=4 for example:

M1,2(3)=[Uncaptioned image]+Dih.=2F3,M_{1,2}^{(3)}=\includegraphics[scale={0.45},align=c]{graph/M123.pdf}+\text{Dih.}=2F_{3}, (41)
M1,3(4)=[Uncaptioned image]+[Uncaptioned image]+Dih.=3(4)+𝔭2a4(4)=3F4,M_{1,3}^{(4)}=\includegraphics[scale={0.45},align=c]{graph/M1341.pdf}+\includegraphics[scale={0.45},align=c]{graph/M1342.pdf}+\text{Dih.}=\mathcal{I}_{3}^{(4)}+\mathfrak{p}_{2a}\,\mathcal{I}_{4}^{(4)}=3F_{4}, (42)

where “Dih.” means all the in-equivalent dihedral permutation images of 4 external legs. The definition of the integral is present in Appendix A and the permutation operator is defined as 𝔭2aXX+(13)\mathfrak{p}_{2a}X\equiv X+(1{\leftrightarrow}3). For later convenience, we also define 𝔭2bXX+(24)\mathfrak{p}_{2b}X\equiv X+(2{\leftrightarrow}4) and 𝔭4𝔭2a𝔭2b\mathfrak{p}_{4}\equiv\mathfrak{p}_{2a}\mathfrak{p}_{2b}.

Starting from M1,2()M_{1,\ell-2}^{(\ell)} (blue arrows in Table. 2), we have contributions beyond ladder-equivalent integrals, which results in the non-trivial magic identities we are interested in. However, we can still “count” and “enumerate” the ladder-equivalent contributions for M1,2()M_{1,\ell-2}^{(\ell)}. In this case, we should add two rungs to an (l2)(l-2)-loop ladder while maintaining the (1,2)(1,\ell-2)-structure. It is easy to see that the ladder contribution to M1,2()M_{1,\ell-2}^{(\ell)} is C2FC_{\ell}^{2}F_{\ell}. Explicitly, for =4\ell=4 and =5\ell=5, we have:

M1,2(4)[Uncaptioned image]+[Uncaptioned image]+Dih.=𝔭41(4)+𝔭2a7(4)=6F4,M_{1,2}^{(4)}\supset\includegraphics[scale={0.45},align=c]{graph/M1241.pdf}+\includegraphics[scale={0.45},align=c]{graph/M1242.pdf}+\text{Dih.}=\mathfrak{p}_{4}\,\mathcal{I}_{1}^{(4)}+\mathfrak{p}_{2a}\,\mathcal{I}_{7}^{(4)}=6F_{4}, (43)
M1,35\displaystyle M_{1,3}^{5} [Uncaptioned image]+[Uncaptioned image]+[Uncaptioned image]+[Uncaptioned image]+[Uncaptioned image]+Dih.\displaystyle\supset\includegraphics[scale={0.45},align=c]{graph/M1352.pdf}+\includegraphics[scale={0.45},align=c]{graph/M1351.pdf}+\includegraphics[scale={0.45},align=c]{graph/M1355.pdf}+\includegraphics[scale={0.45},align=c]{graph/M1354.pdf}+\includegraphics[scale={0.45},align=c]{graph/M1353.pdf}+\text{Dih.}
=𝔭411(5)+𝔭2b13(5)+15(5)+17(5)+𝔭2a19(5)=10F5.\displaystyle=\mathfrak{p}_{4}\,\mathcal{I}_{11}^{(5)}+\mathfrak{p}_{2b}\,\mathcal{I}_{13}^{(5)}+\mathcal{I}_{15}^{(5)}+\mathcal{I}_{17}^{(5)}+\mathfrak{p}_{2a}\,\mathcal{I}_{19}^{(5)}=10F_{5}. (44)

The first non-trivial magic identity starts from =4\ell=4 and has been given in Caron-Huot:2021usw . Note that M1,2(2)+M2,2(2)=6F4M_{1,2}^{(2)}+M_{2,2}^{(2)}=6F_{4}, and the ladder-equivalent contribution in (43) exactly corresponds to 6F46F_{4}. Thus we are left with

0=[Uncaptioned image][Uncaptioned image][Uncaptioned image]+Dih.=𝔭2a2(4)𝔭2b5(4)6(4).0=\includegraphics[scale={0.45},align=c]{graph/DCI4l2.pdf}-\includegraphics[scale={0.45},align=c]{graph/DCI4l5.pdf}-\includegraphics[scale={0.45},align=c]{graph/DCI4l6.pdf}+\text{Dih.}=\mathfrak{p}_{2a}\,\mathcal{I}^{(4)}_{2}-\mathfrak{p}_{2b}\,\mathcal{I}^{(4)}_{5}-\mathcal{I}^{(4)}_{6}. (45)

By attaching box (we call this operation “inverse-boxing”) to this magic identity, we get higher loop magic identities for M1,2()M_{1,\ell-2}^{(\ell)}. And they are exactly magic identities stemming from M1,2(5)M_{1,\ell-2}^{(\ell\geq 5)} after subtracting all ladder-equivalent integrals. For =5\ell=5, the new magic identity reads

0=[Uncaptioned image][Uncaptioned image][Uncaptioned image]+Dih.=𝔭4(2(5)30(5))𝔭2a8(5).0=\includegraphics[scale={0.45},align=c]{graph/magicM1351.pdf}-\includegraphics[scale={0.45},align=c]{graph/magicM1352.pdf}-\includegraphics[scale={0.45},align=c]{graph/magicM1353.pdf}+\text{Dih.}=\mathfrak{p}_{4}\left(\mathcal{I}_{2}^{(5)}-\mathcal{I}_{30}^{(5)}\right)-\mathfrak{p}_{2a}\,\mathcal{I}_{8}^{(5)}. (46)

For =6\ell=6, we can group the integrals by places where the boxes are attached.

0=\displaystyle 0= ([Uncaptioned image][Uncaptioned image][Uncaptioned image]+Dih.)\displaystyle\left(\includegraphics[scale={0.45},align=c]{graph/magic611.pdf}-\includegraphics[scale={0.45},align=c]{graph/magic612.pdf}-\includegraphics[scale={0.45},align=c]{graph/magic613.pdf}+\text{Dih.}\right) (47)
+\displaystyle+ ([Uncaptioned image][Uncaptioned image][Uncaptioned image]+Dih.)\displaystyle\left(\includegraphics[scale={0.45},align=c]{graph/magic621.pdf}-\includegraphics[scale={0.45},align=c]{graph/magic622.pdf}-\includegraphics[scale={0.45},align=c]{graph/magic623.pdf}+\text{Dih.}\right)

And it is easy to generalize this kind of magic identity to higher loops by inverse-boxing.

Now we arrive at M1,3()M_{1,\ell-3}^{(\ell)} (orange arrows in Table. 2). Let us firstly look at =5\ell=5. In Table. 2, notice that M1,2(5)+M2,2(5)=20F5M_{1,2}^{(5)}+M_{2,2}^{(5)}=20F_{5}, where all the ladder-equivalent integrals contribute 20F520F_{5}. Subtracting them gives the 55-loop magic identity already mentioned in Caron-Huot:2021usw . We rewrite it in our notation as

0\displaystyle 0 =[Uncaptioned image][Uncaptioned image]+[Uncaptioned image]+[Uncaptioned image][Uncaptioned image]\displaystyle=\includegraphics[scale={0.45},align=c]{graph/DCI5l33.pdf}-\includegraphics[scale={0.45},align=c]{graph/DCI5l7.pdf}+\includegraphics[scale={0.45},align=c]{graph/DCI5l14.pdf}+\includegraphics[scale={0.45},align=c]{graph/DCI5l1.pdf}-\includegraphics[scale={0.45},align=c]{graph/DCI5l25.pdf}
+[Uncaptioned image]+[Uncaptioned image][Uncaptioned image][Uncaptioned image]+[Uncaptioned image]\displaystyle+\includegraphics[scale={0.45},align=c]{graph/DCI5l3.pdf}+\includegraphics[scale={0.45},align=c]{graph/DCI5l12.pdf}-\includegraphics[scale={0.45},align=c]{graph/DCI5l24.pdf}-\includegraphics[scale={0.45},align=c]{graph/DCI5l6.pdf}+\includegraphics[scale={0.45},align=c]{graph/DCI5l14.pdf}
+[Uncaptioned image][Uncaptioned image][Uncaptioned image][Uncaptioned image]\displaystyle+\includegraphics[scale={0.45},align=c]{graph/DCI5l18.pdf}-\includegraphics[scale={0.45},align=c]{graph/DCI5l26.pdf}-\includegraphics[scale={0.45},align=c]{graph/DCI5l28.pdf}-\includegraphics[scale={0.45},align=c]{graph/DCI5l29.pdf}
+[Uncaptioned image][Uncaptioned image][Uncaptioned image]+Dih.\displaystyle+\includegraphics[scale={0.45},align=c]{graph/DCI5l34.pdf}-\includegraphics[scale={0.45},align=c]{graph/DCI5l27.pdf}-\includegraphics[scale={0.45},align=c]{graph/DCI5l5.pdf}+\text{Dih.}
=33(5)7(5)+14(5)+1(5)25(5)+13(5)+12(5)24(5)6(5)+4(5)\displaystyle=\,\mathcal{I}^{(5)}_{33}-\mathcal{I}^{(5)}_{7}+\mathcal{I}^{(5)}_{14}+\mathcal{I}^{(5)}_{1}-\mathcal{I}^{(5)}_{25}+\mathcal{I}^{(5)}_{13}+\mathcal{I}^{(5)}_{12}-\mathcal{I}^{(5)}_{24}-\mathcal{I}^{(5)}_{6}+\mathcal{I}^{(5)}_{4}
+18(5)26(5)28(5)29(5)+34(5)27(5)5(5)+Dih.,\displaystyle\hskip 10.00002pt+\mathcal{I}^{(5)}_{18}-\mathcal{I}^{(5)}_{26}-\mathcal{I}^{(5)}_{28}-\mathcal{I}^{(5)}_{29}+\mathcal{I}^{(5)}_{34}-\mathcal{I}^{(5)}_{27}-\mathcal{I}^{(5)}_{5}+\text{Dih.}, (48)

where the diagrams only show the denominators, and the corresponding numerators are defined explicitly in Appendix A. We will use the inverse-boxing of the above equation to simplify a magic identity at =6\ell=6. Explicitly, after subtracting the ladder-equivalent integrals from M1,3(6)+M2,3(6)=35F6M_{1,3}^{(6)}+M_{2,3}^{(6)}=35F_{6} 555The ladder-equivalent integrals contribute 36F636F_{6}, that’s why the LHS. of (49) is nonzero., we get an identity involving 125 integrals (37 dihedral seeds). Applying inverse-boxing from different direction to (3.2), i.e. attaching box as in (46), we will have a 6-loop identity. After subtracting it from the 125 integral identity, we have the following nice identity:

F6=\displaystyle-F_{6}= [Uncaptioned image][Uncaptioned image]+[Uncaptioned image]+[Uncaptioned image]\displaystyle-\includegraphics[scale={0.45},align=c]{graph/magic61.pdf}-\includegraphics[scale={0.45},align=c]{graph/magic62.pdf}+\includegraphics[scale={0.45},align=c]{graph/magic63.pdf}+\includegraphics[scale={0.45},align=c]{graph/magic64.pdf} (49)
+[Uncaptioned image]+[Uncaptioned image][Uncaptioned image]+[Uncaptioned image]+Dih.\displaystyle+\includegraphics[scale={0.45},align=c]{graph/magic65.pdf}+\includegraphics[scale={0.45},align=c]{graph/magic66.pdf}-\includegraphics[scale={0.45},align=c]{graph/magic67.pdf}+\includegraphics[scale={0.45},align=c]{graph/magic68.pdf}+\text{Dih.}

which only involves 8 integral dihedral seeds, where the dashed lines correspond to numerators. From 70F52F2,4=M1,2(6)+M2,2(6)70F_{5}-2F_{2,4}=M_{1,2}^{(6)}+M_{2,2}^{(6)}, we will have another magic identity which is too long to be explicitly shown. The last comment is that, from the above analysis, we can detect a diagonal structure shown in Table. 2, which indicates the inverse-boxing of magic identity stemming from lower loop can be used to simplify higher loop ones when they lie on the same diagonal.

3.3 Comments on periods and integrated correlators

Before proceeding, we briefly comment on the periods of conformal integrals, which have been studied quite extensively in the literature Broadhurst:1995km ; Schnetz:2008mp ; Brown:2009ta ; Brown:2013gia ; Schnetz:2013hqa ; Panzer:2016snt ; Panzer:2014caa ; Borinsky:2020rqs ; Borinsky:2022lds partly due to the connection to integrated correlator Binder:2019jwn ; Chester:2019pvm ; Chester:2020vyz ; Dorigoni:2021guq ; Dorigoni:2022zcr ; Paul:2022piq ; Wen:2022oky ; Brown:2023zbr . Recall that for an integrand I()I^{(\ell)} with (4+)(4+\ell) vertex666Here I()I^{(\ell)} denotes the integrand of DCI integrals, it should be distinguished from In,l(l)I_{n,l}^{(l)} in the octagons, which always has two indices in the subscript., its period is defined as integrating over all spacetime coordinates with specific measure Wen:2022oky ; Brown:2023zbr

𝒫(I())=1π2(+1)d4x1d4x+4vol(SO(2,4))I()\mathcal{P}(I^{(\ell)})=\frac{1}{\pi^{2(\ell+1)}}\int\frac{\mathrm{d}^{4}x_{1}\cdots\mathrm{d}^{4}x_{\ell+4}}{\mathrm{vol}(SO(2,4))}I^{(\ell)} (50)

where the vol(SO(2,4))\mathrm{vol}(SO(2,4)) allow us to fix three xix_{i}. Since we are integrating over all spacetime coordinates, periods of all conformal integrals coming from one ff-graph are equal. So we only need to consider periods of ff-graphs, which is related to the periods of integrals by

𝒫(fi())=(4+)!|aut(fi())|𝒫(Ii()),\mathcal{P}(f_{i^{\prime}}^{(\ell)})=\frac{(4+\ell)!}{\left|\mathrm{aut}(f_{i^{\prime}}^{(\ell)})\right|}\mathcal{P}(I_{i}^{(\ell)}), (51)

where Ii()I_{i}^{(\ell)} is one of the conformal integrals from the ff-graph fi()f_{i^{\prime}}^{(\ell)}.

As perhaps the simplest constraints, the periods of integrals related by the original magic identity are all equal. Now we consider the “slingshot” (23) rule at the ff-graph level: recall that integrals Ii(){I}_{i}^{(\ell)} can be obtained from certain ff-graph fi()f_{i^{\prime}}^{(\ell)} by multiplying x122x232x342x412x132x242x_{12}^{2}x_{23}^{2}x_{34}^{2}x_{41}^{2}x_{13}^{2}x_{24}^{2}, (23) can be translated to

fj(+1)\displaystyle f_{j^{\prime}}^{(\ell+1)} =1x+5,22x232x342x4,+52x+5,32x242×x242x122x412x+5,12×(x122x232x342x412x132x242)fi()\displaystyle=\frac{1}{x_{\ell+5,2}^{2}x_{23}^{2}x_{34}^{2}x_{4,\ell+5}^{2}x_{\ell+5,3}^{2}x_{24}^{2}}\times\frac{x_{24}^{2}}{x_{12}^{2}x_{41}^{2}x_{\ell+5,1}^{2}}\times\left(x_{12}^{2}x_{23}^{2}x_{34}^{2}x_{41}^{2}x_{13}^{2}x_{24}^{2}\right)f_{i^{\prime}}^{(\ell)} (52)
=x132x242x+5,12x+5,22x+5,32x+5,42fi(),\displaystyle=\frac{x_{13}^{2}x_{24}^{2}}{x_{\ell+5,1}^{2}x_{\ell+5,2}^{2}x_{\ell+5,3}^{2}x_{\ell+5,4}^{2}}f_{i^{\prime}}^{(\ell)},

which is just the inverse square rule (25). It should be noted that the four-cycle 12341-2-3-4 in fif_{i^{\prime}} correspond to DCI integrals involved in the \ell-loop magic identity. The periods of all fj(+1)f^{(\ell+1)}_{j^{\prime}} are equal up to the symmetry factor.

Moreover, by integrating over all xix_{i} in (39), we obtain some constraints between the periods of ff-graphs.

i=1𝒩ci()|aut(fi())||a,b|𝒫(fi())=n=min(a,b)a+b1l0l+na,bCa,bn,l𝒫(In,l()),\sum_{i=1}^{\mathcal{N}_{\ell}}c_{i}^{(\ell)}\left|\mathrm{aut}(f_{i}^{(\ell)})\right|\left|\mathbb{P}_{a,b}\right|\mathcal{P}(f_{i}^{(\ell)})=\sum_{n=\min(a,b)}^{a+b-1}\sum_{\begin{subarray}{c}l\geq 0\\ l+n\geq a,b\end{subarray}}^{\infty}C^{n,l}_{a,b}\mathcal{P}(I_{n,l}^{(\ell)}), (53)

where |a,b|\left|\mathbb{P}_{a,b}\right| is the number of (a,b)(a,b)-structures in the ff-graph. For example, at four loops, 2(4)\mathcal{I}^{(4)}_{2}, 5(4)\mathcal{I}^{(4)}_{5} and 6(4)\mathcal{I}^{(4)}_{6} in (45) comes from f1(4),f2(4),f2(4)f^{(4)}_{1},f^{(4)}_{2},f^{(4)}_{2}, thus we have

2×8×𝒫(f1(4))2×16×𝒫(f2(4))=16×𝒫(f2(4)),2\times 8\times\mathcal{P}(f_{1}^{(4)})-2\times 16\times\mathcal{P}(f_{2}^{(4)})=16\times\mathcal{P}(f_{2}^{(4)}), (54)

which can be easily verified since 𝒫(f1(4))=8!8×252ζ9,𝒫(f2(4))=8!16×168ζ9\mathcal{P}({f_{1}^{(4)}})=\frac{8!}{8}\times 252\zeta_{9},\mathcal{P}({f_{2}^{(4)}})=\frac{8!}{16}\times 168\zeta_{9}. Similarly, the constraints between periods of five-loop ff-graphs reduce all seven periods to three independent ones Wen:2022oky . For example, the identity M2,2(5)=6F1,4M_{2,2}^{(5)}=-6F_{1,4} in Table. 2 gives a linear relation (see also (73)) among 𝒫(f1(5))\mathcal{P}(f_{1}^{(5)}), 𝒫(f5(5))\mathcal{P}(f_{5}^{(5)}) and 𝒫(f7(5))\mathcal{P}(f_{7}^{(5)}) (the order of five-loop ff-graphs is defined as in Eden:2012tu ):

89!𝒫(f5(5))+2×129!𝒫(f7(5))8×29!𝒫(f1(5))=5×924ζ11.-\frac{8}{9!}\mathcal{P}(f_{5}^{(5)})+2\times\frac{12}{9!}\mathcal{P}(f_{7}^{(5)})-8\times\frac{2}{9!}\mathcal{P}(f_{1}^{(5)})=-5\times 924\zeta_{11}. (55)

This relation supplements the constraint from localization result of integrated correlator presented in Wen:2022oky . For periods of 26 six-loop ff-graphs with non-zero coefficients, the relations between Coulomb branch amplitudes and octagons also play a key role in determining some of them. Here we take M1,3(6)+M2,3(6)=35F6M_{1,3}^{(6)}+M_{2,3}^{(6)}=35F_{6} as an example, which is simplified to (49). Considering the corresponding ff-graphs on the RHS and the integrated result of FpF_{p} is given by C2p+2p+1ζ2p+1C_{2p+2}^{p+1}\zeta_{2p+1} , thus we obtain the following constraint on the integrated result:

4×110!𝒫(f25(6))4×210!𝒫(f17(6))+2×210!𝒫(f8(6))+2×210!𝒫(f2(6))+2×410!𝒫(f1(6))\displaystyle-4\times\frac{1}{10!}\mathcal{P}(f_{25}^{(6)})-4\times\frac{2}{10!}\mathcal{P}(f_{17}^{(6)})+2\times\frac{2}{10!}\mathcal{P}(f_{8}^{(6)})+2\times\frac{2}{10!}\mathcal{P}(f_{2}^{(6)})+2\times\frac{4}{10!}\mathcal{P}(f_{1}^{(6)}) (56)
+2×210!𝒫(f23(6))2×210!𝒫(f6(6))+410!𝒫(f14(6))=3432ζ13,\displaystyle+2\times\frac{2}{10!}\mathcal{P}(f_{23}^{(6)})-2\times\frac{2}{10!}\mathcal{P}(f_{6}^{(6)})+\frac{4}{10!}\mathcal{P}(f_{14}^{(6)})=-3432\zeta_{13},

which provides a non-trivial linear relation between 𝒫(f1(6))\mathcal{P}(f_{1}^{(6)}) and other periods. It supplements the constraints derived from simplification of correlators by substitution dijγiγjd_{ij}\to\gamma_{i}\gamma_{j}, as described in Brown:2023zbr . It also serves as a consistency check of our results since we already know the periods of these ff-graphs from Brown:2023zbr .

4 Bootstrapping DCI integrals at four loops and beyond

In this section, we turn our attention to individual DCI integrals (as opposed to those special combinations studied above), which have only been systematically studied up to three loops Drummond:2013nda . Starting at four loops, although some of these DCI integrals are rather trivial, e.g. equivalent to the ladder integral (or generalized ladders Usyukina:1993ch ; Broadhurst:2010ds ; Drummond:2012bg ), the remaining ones have not been computed with a few exceptions such as fishnet integrals Basso:2017jwq ; Basso:2021omx ; Aprile:2023gnh , e.g. the four-loop combination M2,2M_{2,2} consists of a single integral denoted as f2{\cal I}^{f2} in Caron-Huot:2021usw . As we have mentioned, some four-loop conformal integrals involve elliptic cuts which prevent us from computing all of them, but it turns out to be much easier for DCI integrals at four and even higher loops.

Based on our analysis of leading singularities and the very restrictive function space, we will see that it is straightforward to bootstrap the 33 four-loop DCI integrals that are not equal to ladder integral, which also provides an independent check for the four-loop magic identity. Moreover, we will see that the same bootstrap method can be applied to five-loop DCI integrals. We will again enumerate all possible leading singularities and also obtain some results for the pure functions, though the latter become more complicated: there are also lower-weight pieces for some of these integrals, and in general they involve functions beyond harmonic polylogarithms.

4.1 All four-loop DCI integrals

Among the magic identities we have studied in sec. 3, there is a simple enough but nontrivial one

d2(u,v)+d2(v,u)+f2(u,v)=f(u,v)+f(v,u).\mathcal{I}^{d2}(u,v)+\mathcal{I}^{d2}(v,u)+\mathcal{I}^{f2}(u,v)=\mathcal{I}^{f}(u,v)+\mathcal{I}^{f}(v,u). (57)

Here, to maintain consistency with the notation in Caron-Huot:2021usw , we rename these function as follows.

d25(4),f26(4),f2(4).\mathcal{I}^{d2}\equiv\mathcal{I}^{(4)}_{5},\,\mathcal{I}^{f2}\equiv\mathcal{I}^{(4)}_{6},\,\mathcal{I}^{f}\equiv\mathcal{I}^{(4)}_{2}. (58)

The arguments of f(u,v)\mathcal{I}^{f}(u,v) specify the order of external legs. For example, f(v,u)\mathcal{I}^{f}(v,u) indicates external legs 11 and 33 are exchanged in the original definition of f(u,v)\mathcal{I}^{f}(u,v). While the analytic expressions for d2\mathcal{I}^{d2} and f2\mathcal{I}^{f2} (the fishnet) have already been calculated and involve only SVHPL functions Drummond:2013nda ; Eden:2016dir , f(u,v)\mathcal{I}^{f}(u,v) remains unknown. (57) gives the combination of f(u,v)\mathcal{I}^{f}(u,v) and f(v,u)\mathcal{I}^{f}(v,u), offering constraints on bootstrapping the analytic form of f\mathcal{I}^{f}. The above magic identity suggests that f\mathcal{I}^{f} should also be expressible in terms of SVHPLs with same prefactors in d2\mathcal{I}^{d2} and f2\mathcal{I}^{f2}. Additionally, the asymptotic expansion of f\mathcal{I}^{f} in various kinematic regions, as provided in Chicherin:2018avq , can be used in the bootstrap process.

Before proceeding with the bootstrap, we review some properties of SVHPL and analyze the leading singularities of f\mathcal{I}^{f}. The symbol alphabet of SVHPLs consists of {z,z¯,1z,1z¯}\{z,\bar{z},1-z,1-\bar{z}\}, and SVHPLs satisfy the property that they have no branch cuts in the whole complex plane. However, this does not imply that they are purely real or imaginary. Nevertheless, we can construct such basis which have definite parity under exchanging zz and z¯\bar{z}:

parity odd: odd(z,z¯)=odd(z¯,z),\displaystyle\mathcal{L}^{odd}(z,\bar{z})=-\mathcal{L}^{odd}(\bar{z},z), (59)
parity even: even(z,z¯)=even(z¯,z)\displaystyle\mathcal{L}^{even}(z,\bar{z})=\mathcal{L}^{even}(\bar{z},z)

This basis has been constructed up to weight 8 using HyperlogProcedure. The number of independent SVHPL basis functions for each weight nn is summarized in Table. 3.

weight 1 2 3 4 5 6 7 8
# of parity odd 0 1 2 6 12 28 56 120
# of parity even 2 3 6 10 20 36 72 136
# of total 2 4 8 16 32 64 128 256
Table 3: The counting of SVHPL up to weight 8.

It has been proven that the total number is 2n2^{n} at weight nn, while the number of parity even is 2L1+2L122^{L-1}+2^{\lfloor\frac{L-1}{2}\rfloor} Brown:2004ugm ; Schnetz:2013hqa ; Alday:2024ksp . This basis can be directly constructed using the method described in Schnetz:2013hqa and is implemented in the package HyperlogProcedures up to very high weights. Although the expressions for d2\mathcal{I}^{d2} and f2\mathcal{I}^{f2} are available in the literature, as mentioned earlier, they can also be bootstrapped using the asymptotic expansions provided in Chicherin:2018avq . The ansatz for d2\mathcal{I}^{d2} or f2\mathcal{I}^{f2} can be set up as

ansd2=iodd,iw=8Δ,ansf2=ieven,iw=8Δ2,\mathcal{I}^{d2}_{ans}=\sum_{i}\frac{\mathcal{L}^{w=8}_{odd,i}}{\Delta},\quad\mathcal{I}^{f2}_{ans}=\sum_{i}\frac{\mathcal{L}^{w=8}_{even,i}}{\Delta^{2}}, (60)

where Δ\Delta is defined as

Δzz¯=12u2v+u22uv+v2.\Delta\equiv z-\bar{z}=\sqrt{1-2u-2v+u^{2}-2uv+v^{2}}. (61)

Note that w=8\mathcal{L}^{w=8}_{\bullet} indicates all SVHPL functions with total weight 8 which contain functions like ζ3w=5\zeta_{3}\mathcal{L}^{w=5} where w=5\mathcal{L}^{w=5} consists of weight-5 SVHPL777This counting of w=n\mathcal{L}^{w=n} in Tab. 3 at weight nn does not include terms such as ζmw=nm,m>0\zeta_{m}\mathcal{L}^{w=n-m},m>0.. Using the asymptotic expansions around (u,v)(0,1)(u,v)\sim(0,1), (1/u,v/u)(0,1)(1/u,v/u)\sim(0,1) and (1/v,u/v)(0,1)(1/v,u/v)\sim(0,1), the results are determined to be888In the definition of d2\mathcal{I}^{d2} and f2\mathcal{I}^{f2} of Chicherin:2018avq , the external legs are ordered as 1,3,2,41,3,2,4. Therefore, in the expressions for d2\mathcal{I}^{d2} and f\mathcal{I}^{f}, zz and z¯\bar{z} should be replaced by 1/z1/z and 1/z¯1/\bar{z}, and an overall factor zz¯z\bar{z} should be factored out to get the expressions with order 1,2,3,41,2,3,4.

d2(vu,1u)=1Δ(Iz,1,0,1,0,1,1,0,1,0Iz,1,0,1,1,0,1,0,1,015ζ5Iz,1,0,1,04418ζ7Iz,1,0),\displaystyle\mathcal{I}^{d2}(\frac{v}{u},\frac{1}{u})=\frac{1}{\Delta}\left(\mathrm{I}_{z,1,0,1,0,1,1,0,1,0}-\mathrm{I}_{z,1,0,1,1,0,1,0,1,0}-15\zeta_{5}\mathrm{I}_{z,1,0,1,0}-\frac{441}{8}\zeta_{7}\mathrm{I}_{z,1,0}\right), (62)
f2(1u,vu)=4zz¯Δ2(Iz,0,1,1,0,1,1,1,1,0Iz,0,1,1,1,1,0,1,1,0Iz,1,0,1,1,0,1,1,1,0+Iz,1,0,1,1,1,1,0,1,0\displaystyle\mathcal{I}^{f2}(\frac{1}{u},\frac{v}{u})=\frac{4z\bar{z}}{\Delta^{2}}\Big{(}\mathrm{I}_{z,0,1,1,0,1,1,1,1,0}-\mathrm{I}_{z,0,1,1,1,1,0,1,1,0}-\mathrm{I}_{z,1,0,1,1,0,1,1,1,0}+\mathrm{I}_{z,1,0,1,1,1,1,0,1,0}
+Iz,1,1,0,1,1,0,1,1,0Iz,1,1,0,1,1,1,1,0,0Iz,1,1,1,0,1,1,0,1,0+Iz,1,1,1,1,0,1,1,0,0+2ζ3Iz,0,1,1,1,1,0\displaystyle+\mathrm{I}_{z,1,1,0,1,1,0,1,1,0}-\mathrm{I}_{z,1,1,0,1,1,1,1,0,0}-\mathrm{I}_{z,1,1,1,0,1,1,0,1,0}+\mathrm{I}_{z,1,1,1,1,0,1,1,0,0}+2\zeta_{3}\mathrm{I}_{z,0,1,1,1,1,0}
+4ζ3Iz,1,0,1,1,1,04ζ3Iz,1,1,1,0,1,02ζ3Iz,1,1,1,1,0,0+10ζ5Iz,0,1,1,010ζ5Iz,1,1,0,06ζ32Iz,1,1,0),\displaystyle+4\zeta_{3}\mathrm{I}_{z,1,0,1,1,1,0}\!-\!4\zeta_{3}\mathrm{I}_{z,1,1,1,0,1,0}\!-\!2\zeta_{3}\mathrm{I}_{z,1,1,1,1,0,0}\!+\!10\zeta_{5}\mathrm{I}_{z,0,1,1,0}\!-\!10\zeta_{5}\mathrm{I}_{z,1,1,0,0}\!-\!6\zeta_{3}^{2}\mathrm{I}_{z,1,1,0}\Big{)},

We apply the function Iz,,0\mathrm{I}_{z,...,0} defined within HyperlogProcedures as a basis for SVHPLs. Iz,,0\mathrm{I}_{z,...,0} is defined recursively as follows Schnetz:2013hqa :

Iz,an,an1,,a1,0=svdzzanIz,an1,,a1,0,an=0,1;\displaystyle\mathrm{I}_{z,a_{n},a_{n-1},\ldots,a_{1},0}=\int_{sv}\frac{\mathrm{d}z}{z-a_{n}}\mathrm{I}_{z,a_{n-1},\ldots,a_{1},0},\,a_{n}=0,1; (63)
Iz,0=1;Iz,0,0=logzz¯,Iz,1,0=log(1z)(1z¯).\displaystyle\mathrm{I}_{z,0}=1;\mathrm{I}_{z,0,0}=\log z\bar{z},\,\mathrm{I}_{z,1,0}=\log(1-z)(1-\bar{z}).

sv\int_{sv} is the single-valued integration defined by Eq. (2.53) in Schnetz:2013hqa , which not only is an integration of zz, but also accounts for the antiholomorphic part in z¯\bar{z} and the boundary condition which requires the integration vanishes when zz approaches 0 (reguralized around z=0z=0). It keeps the result single-valued. This definition can be taken as the single-valued version for HPL functions (here, we adopt the definition within HyperlogProcedures for consistency, which is slightly different from the standard definition in Remiddi:1999ew )

iz,an,an1,,a1,0=0zdzzaniz,an1,,a1,0,an=0,±1;\displaystyle\mathrm{i}_{z,a_{n},a_{n-1},\ldots,a_{1},0}=\int_{0}^{z}\frac{\mathrm{d}z^{\prime}}{z^{\prime}-a_{n}}\mathrm{i}_{z,a_{n-1},\ldots,a_{1},0},\,a_{n}=0,\pm 1; (64)
iz,0=1.\displaystyle\mathrm{i}_{z,0}=1.

where iz,,0\mathrm{i}_{z,...,0} is undertood as already regularized around z=0z=0. SVHPLs can be generally expressed as products of HPLs:

Iz,a,0=i,jci,jiz,ai,0iz¯,aj,0.\mathrm{I}_{z,\vec{a},0}=\sum_{i,j}c_{i,j}\mathrm{i}_{z,\vec{a}_{i},0}\mathrm{i}_{\bar{z},\vec{a}_{j},0}. (65)

For instance, Iz,1,0,0\mathrm{I}_{z,1,0,0} is a weight-2 function given by iz,1,0iz¯,0,0+iz,1,0,0+iz¯,0,1,0\mathrm{i}_{z,1,0}\mathrm{i}_{\bar{z},0,0}+\mathrm{i}_{z,1,0,0}+\mathrm{i}_{\bar{z},0,1,0}. iz,1,0,0\mathrm{i}_{z,1,0,0} is a HPL function with the symbol z(1z)z\otimes(1-z). (62) agrees with direct computations using HyperlogProcedures. For example, d2\mathcal{I}^{d2} can be computed in HyperlogProcedures using the following syntax:

    int:={[{1,5},1],[{2,6},1],[{3,5},1],[{3,7},1],[{5,7},1],[{5,8},1],
    [{7,8},1],[{6,7},1],[{6,8},1]}:
    Id2:=GraphicalFunction(int,[1,3,2],4,0,[z,zz]);
where the first line defines the integrand using edges of corresponding graph. The basic element is {edge,weight}. For example, [{1,5},1] represents x152x_{15}^{2} in the denominator. A weight of -1 would indicate a numerator term. The second line computes the integral using the graph function method described in Schnetz:2013hqa . We have set x4x_{4} to infinity and fixed x1=0,x3=1,x2=zx_{1}=0,x_{3}=1,x_{2}=z999For more details on the function’s usage and the meaning of every argument, refer to the documentation of HyperlogProcedures.. zz is a shorthand for z¯\bar{z}. f2\mathcal{I}^{f2} is a fishnet integral whose expression can be obtained by the determinant of ladders FnF_{n} Basso:2017jwq which has a rather compact form expressed using Iz,,0\mathrm{I}_{z,...,0}:
Fn(uv,1v)=(1)nΔ×(Iz,0,,0n,1,0,,0n1,0Iz,0,,0n1,1,0,,0n,0).F_{n}(\frac{u}{v},\frac{1}{v})=\frac{(-1)^{n}}{\Delta}\times\big{(}\mathrm{I}_{z,{\scriptscriptstyle\underbrace{0,...,0}_{n}},1,{\scriptscriptstyle\underbrace{0,...,0}_{n-1}},0}-\mathrm{I}_{z,{\scriptscriptstyle\underbrace{0,...,0}_{n-1}},1,{\scriptscriptstyle\underbrace{0,...,0}_{n}},0}\big{)}. (66)

The expression of all fishnets have been implemented into HyperlogProcedures.

Although d2\mathcal{I}^{d2} and f2\mathcal{I}^{f2} can be computed directly in this manner, f\mathcal{I}^{f} cannot be evaluated using HyperlogProcedures101010f2\mathcal{I}^{f2} is a fishnet, and its result is pre-implemented in HyperlogProcedures Borinsky:2022lds . Thus, it is not explicitly calculated. In this sense, our result for f\mathcal{I}^{f} can serve as a seed for HyperlogProcedures. . This necessitates the usage of the bootstrap method. Now we bootstrap f\mathcal{I}^{f} in the same way as above. To construct an ansatz for f\mathcal{I}^{f}, we first analyze its leading singularities (or equivalently, construct dlogd\log forms for the integrand).

f=d4x8d4x7d4x6d4x5x182x372x242x522x532x572x582x612x642x672x682x712x722x832x842x782\mathcal{I}^{f}=\int\mathrm{d}^{4}x_{8}\mathrm{d}^{4}x_{7}\mathrm{d}^{4}x_{6}\mathrm{d}^{4}x_{5}\frac{x_{18}^{2}x_{37}^{2}x_{24}^{2}}{x_{52}^{2}x_{53}^{2}x_{57}^{2}x_{58}^{2}x_{61}^{2}x_{64}^{2}x_{67}^{2}x_{68}^{2}x_{71}^{2}x_{72}^{2}x_{83}^{2}x_{84}^{2}x_{78}^{2}} (67)

It (2(4)\mathcal{I}_{2}^{(4)}) is depicted in Appendix. A. After first integrating out x5x_{5} and x6x_{6}, the remaining integrand takes the form

I78f=x182x372x242x172x272x382x482x782λ2378λ1478.I^{f}_{78}=\frac{x_{18}^{2}x_{37}^{2}x_{24}^{2}}{x_{17}^{2}x_{27}^{2}x_{38}^{2}x_{48}^{2}x_{78}^{2}\lambda_{2378}\lambda_{1478}}. (68)

where λijkl=det(xmn2)m,n=i,j,k,l\lambda_{ijkl}=\sqrt{\det(x_{mn}^{2})_{m,n=i,j,k,l}} represents the leading singularity of a four-mass box. Now the question is how many leading singularities can be derived from I78fI^{f}_{78}. Since the integrand is symmetric of x7x_{7} and x8x_{8}, we can first cut x7x_{7} and then x8x_{8}. The cut for x7x_{7} and the corresponding remaining terms are summarized in Table. 4.

cut for x7x_{7} remaining term cut for x8x_{8} leading singularity
x172,x272,x782,λ1478x_{17}^{2},x_{27}^{2},x_{78}^{2},\lambda_{1478} x242x282x382x482λ1248\displaystyle\frac{x_{24}^{2}}{x_{28}^{2}x_{38}^{2}x_{48}^{2}\lambda_{1248}} x282,x382,x482,λ1248x_{28}^{2},x_{38}^{2},x_{48}^{2},\lambda_{1248} 1Δ\displaystyle\frac{1}{\Delta}
x172,x272,x782,λ2378x_{17}^{2},x_{27}^{2},x_{78}^{2},\lambda_{2378} 0111111The cut vanishes due to numerator x372x_{37}^{2}. Same reason for other 0’s in the table. 0 0
x172,x782,λ2378,λ1478x_{17}^{2},x_{78}^{2},\lambda_{2378},\lambda_{1478} x242x282x482λ148{23}\displaystyle\frac{x_{24}^{2}}{x_{28}^{2}x_{48}^{2}\lambda_{148\{23\}}} x182,x282,x382,x482x_{18}^{2},x_{28}^{2},x_{38}^{2},x_{48}^{2} 1Δ\displaystyle\frac{1}{\Delta}
x272,x782,λ2378,λ1478x_{27}^{2},x_{78}^{2},\lambda_{2378},\lambda_{1478} 0 0 0
x172,x272,λ2378,λ1478x_{17}^{2},x_{27}^{2},\lambda_{2378},\lambda_{1478} x182x232x242x282x382x482λ12{348}\displaystyle\frac{x_{18}^{2}x_{23}^{2}x_{24}^{2}}{x_{28}^{2}x_{38}^{2}x_{48}^{2}\lambda_{12\{348\}}} x182,x282,x382,x482\displaystyle x_{18}^{2},x_{28}^{2},x_{38}^{2},x_{48}^{2} 1Δ\displaystyle\frac{1}{\Delta} uv1Δ2\displaystyle\frac{u-v-1}{\Delta^{2}}
Table 4: Cuts for I78fI_{78}^{f}. λ148{23}2=[x182(x382x242x282x342)x482(x382x122x282x132)]24x142x232x182x282x382x482\lambda^{2}_{148\{23\}}=[x_{18}^{2}(x_{38}^{2}x_{24}^{2}\!-\!x_{28}^{2}x_{34}^{2})\!-\!x_{48}^{2}(x_{38}^{2}x_{12}^{2}\!-\!x_{28}^{2}x_{13}^{2})]^{2}-4x_{14}^{2}x_{23}^{2}x_{18}^{2}x_{28}^{2}x_{38}^{2}x_{48}^{2} and λ12{348}2=[x122(x182x282x342x142x342x282x182x482x232)(x142x282x182x242)(x232x182x282x132)]24x124x142x232x182x282x382x482\lambda_{12\{348\}}^{2}=[x_{12}^{2}(x_{18}^{2}x_{28}^{2}x_{34}^{2}-x_{14}^{2}x_{34}^{2}x_{28}^{2}-x_{18}^{2}x_{48}^{2}x_{23}^{2})-(x_{14}^{2}x_{28}^{2}-x_{18}^{2}x_{24}^{2})(x_{23}^{2}x_{18}^{2}-x_{28}^{2}x_{13}^{2})]^{2}-4x_{12}^{4}x_{14}^{2}x_{23}^{2}x_{18}^{2}x_{28}^{2}x_{38}^{2}x_{48}^{2}.

Two λ\lambdas come from the Gram determinants:

λ148{23}2G(x71μ,x74μ,x78μ,x382x72μx282x73μ)|cut,\displaystyle\lambda^{2}_{148\{23\}}\equiv G(x_{71}^{\mu},x_{74}^{\mu},x_{78}^{\mu},x_{38}^{2}x_{72}^{\mu}-x_{28}^{2}x_{73}^{\mu})|_{cut},\, (69)
cut: x712=x742=x782=x382x722x282x732=0.\displaystyle\text{cut: }x_{71}^{2}=x_{74}^{2}=x_{78}^{2}=x_{38}^{2}x_{72}^{2}-x_{28}^{2}x_{73}^{2}=0.
λ12{348}2G(x71μ,x72μ,x142x78μx182x74μ,x232x78μx282x73μ)|cut,\displaystyle\lambda_{12\{348\}}^{2}\equiv G(x_{71}^{\mu},x_{72}^{\mu},x_{14}^{2}x_{78}^{\mu}-x_{18}^{2}x_{74}^{\mu},x_{23}^{2}x_{78}^{\mu}-x_{28}^{2}x_{73}^{\mu})|_{cut},\,
cut: x712=x722=x142x782x182x742=x232x782x282x732=0.\displaystyle\text{cut: }x_{71}^{2}=x_{72}^{2}=x_{14}^{2}x_{78}^{2}-x_{18}^{2}x_{74}^{2}=x_{23}^{2}x_{78}^{2}-x_{28}^{2}x_{73}^{2}=0.

Using the magic identity (57) and the leading singularities analysis in Table. 4, we propose the following ansatz for f\mathcal{I}^{f}:

ansf=iodd,iw=8Δ+(uv1)jeven,jw=8Δ2.\mathcal{I}_{ans}^{f}=\sum_{i}\frac{\mathcal{L}^{w=8}_{odd,i}}{\Delta}+(u-v-1)\sum_{j}\frac{\mathcal{L}^{w=8}_{even,j}}{\Delta^{2}}. (70)

In this case, we find that the asymptotic expansion alone suffices to determine all parameters. However, we can also directly apply the functional constraints from (57), which must hold at the function level without asymptotic expansion:

1Δ(iodd,iw=8(z,z¯)iodd,iw=8(1z,1z¯))=d2(z,z¯)+d2(1z,1z¯),\displaystyle\frac{1}{\Delta}\left(\sum_{i}\mathcal{L}^{w=8}_{odd,i}(z,\bar{z})-\sum_{i}\mathcal{L}^{w=8}_{odd,i}(1-z,1-\bar{z})\right)=\mathcal{I}^{d2}(z,\bar{z})+\mathcal{I}^{d2}(1-z,1-\bar{z}), (71)
1Δ2((uv1)jeven,jw=8(z,z¯)+(vu1)jeven,jw=8(1z,1z¯))=f2(z,z¯).\displaystyle\frac{1}{\Delta^{2}}\left((u\!-\!v\!-\!1)\sum_{j}\mathcal{L}^{w=8}_{even,j}(z,\bar{z})+(v\!-\!u\!-\!1)\sum_{j}\mathcal{L}^{w=8}_{even,j}(1-z,1-\bar{z})\right)=\mathcal{I}^{f2}(z,\bar{z}).

Note that to apply the above identity, the external points must be labeled consistently to ensure that zz and z¯\bar{z} are defined in the same way. While the magic identity alone cannot fix all the parameters, it reduces the number of free parameters to 63. This significantly simplifies the bootstrap process when combined with asymptotic expansion results. For example, applying the magic identity constraint first reduces the required number of asymptotic expansions from four to two. Although more complicated than d2\mathcal{I}^{d2} and f2\mathcal{I}^{f2}, the explicit result for f\mathcal{I}^{f} can still fits into a few lines (with the usual relabeling of external points as x2,x1,x3,x4x_{2},x_{1},x_{3},x_{4}, which simplifies the expression):

f(uv,1v)=uv1Δ2(2Iz,1,0,0,0,0,1,0,0,02Iz,1,0,0,1,0,0,0,0,02Iz,0,0,0,0,1,0,0,1,0\displaystyle\mathcal{I}^{f}(\frac{u}{v},\frac{1}{v})=\frac{u-v-1}{\Delta^{2}}\Bigg{(}2\mathrm{I}_{z,1,0,0,0,0,1,0,0,0}-2\mathrm{I}_{z,1,0,0,1,0,0,0,0,0}-2\mathrm{I}_{z,0,0,0,0,1,0,0,1,0} (72)
+2Iz,0,0,1,0,0,0,0,1,02Iz,0,1,0,0,0,0,1,0,0+2Iz,0,1,0,0,1,0,0,0,02Iz,0,0,1,0,0,1,0,0,0+2Iz,0,0,0,1,0,0,1,0,0)\displaystyle+2\mathrm{I}_{z,0,0,1,0,0,0,0,1,0}\!-\!2\mathrm{I}_{z,0,1,0,0,0,0,1,0,0}\!+\!2\mathrm{I}_{z,0,1,0,0,1,0,0,0,0}\!-\!2\mathrm{I}_{z,0,0,1,0,0,1,0,0,0}\!+\!2\mathrm{I}_{z,0,0,0,1,0,0,1,0,0}\Bigg{)}
1Δ(40ζ5Iz,0,0,1,020ζ5Iz,0,1,0,08ζ3Iz,0,0,0,0,1,0+8ζ3Iz,0,0,0,1,0,0+4ζ3Iz,0,1,0,0,1,0\displaystyle-\frac{1}{\Delta}\Bigg{(}\!40\zeta_{5}\mathrm{I}_{z,0,0,1,0}\!-\!20\zeta_{5}\mathrm{I}_{z,0,1,0,0}\!-\!8\zeta_{3}\mathrm{I}_{z,0,0,0,0,1,0}+8\zeta_{3}\mathrm{I}_{z,0,0,0,1,0,0}+4\zeta_{3}\mathrm{I}_{z,0,1,0,0,1,0}
4ζ3Iz,0,1,0,1,0,02Iz,0,0,0,0,1,0,0,1,0+2Iz,0,0,0,0,1,0,1,0,0+2Iz,0,0,1,0,0,0,0,1,02Iz,0,0,1,0,0,0,1,0,0\displaystyle-4\zeta_{3}\mathrm{I}_{z,0,1,0,1,0,0}-2\mathrm{I}_{z,0,0,0,0,1,0,0,1,0}+2\mathrm{I}_{z,0,0,0,0,1,0,1,0,0}+2\mathrm{I}_{z,0,0,1,0,0,0,0,1,0}-2\mathrm{I}_{z,0,0,1,0,0,0,1,0,0}
+2Iz,0,1,0,0,0,1,0,0,0Iz,0,1,0,0,1,0,1,0,02Iz,0,1,0,1,0,0,0,0,0+Iz,0,1,0,1,0,0,1,0,02Iz,1,0,0,0,0,1,0,0,0\displaystyle+2\mathrm{I}_{z,0,1,0,0,0,1,0,0,0}-\mathrm{I}_{z,0,1,0,0,1,0,1,0,0}-2\mathrm{I}_{z,0,1,0,1,0,0,0,0,0}+\mathrm{I}_{z,0,1,0,1,0,0,1,0,0}-2\mathrm{I}_{z,1,0,0,0,0,1,0,0,0}
+2Iz,1,0,0,1,0,0,0,0,0).\displaystyle+2\mathrm{I}_{z,1,0,0,1,0,0,0,0,0}\Bigg{)}.

This has also been verified numerically using pySecDec Borowka:2017idc with a precision of 𝒪(106)\mathcal{O}(10^{-6}).

4.2 Leading singularities and some results for five-loop DCI integrals

Among the 34 five-loop DCI integrals defined in Appendix A, 23 can be directly computed using HyperlogProcedures. There are several identities among these integrals, and they will be discussed first. Next, we study integrals which are unknown but can be calculated through boxing, i.e. a second-order differential equation Drummond:2010cz . The second-order differential operator of z,z¯z,\bar{z} acting on the integrals will relate them to lower-loop DCI integrals. In contrast, the inverse boxing is the technique for directly solving the second-order differential equation through single-valued integration of lower-loop integrals Schnetz:2013hqa ; Borinsky:2022lds . Two of the unknown integrals can be related to the four-loop integrals f\mathcal{I}^{f} calculated in the previous section. They can no longer be expressed using only SVHPLs and they are not of uniform weight as well, due to the appearance of a leading singularity, 1/Δ21/\Delta^{2}, in the four-loop DCI integrals121212In the inverse boxing of f\mathcal{I}^{f}, 1/Δ21/\Delta^{2} will result in the single-valued integration of the following type: dzdz¯zz¯(1z)(1z¯)(zz¯)×pure functions\int\frac{\mathrm{d}z\mathrm{d}\bar{z}}{z\bar{z}(1-z)(1-\bar{z})(z-\bar{z})}\times\text{pure functions}, the measure of which is not of dlog\mathrm{d}\log type for z,z¯z,\bar{z}. Thus, the integration will result in non-uniform-weight part.. Although the rest unknown integrals can be studied by bootstrap in principle, we have less control over both the corresponding ansatz and required constraints, due to the possible loss of uniform-weight property and the loss of higher-rank asymptotic expansions. Here, we only take a first step in studying all leading singularities of the 34 DCI integrals and present them in Table. 5131313For those integrals that are not of uniform weight, the leading singularity analysis only accounts for the coefficients of the weight-10 part.. Finally, there is a five-loop relation (M2,2(5)=6F1,4M_{2,2}^{(5)}=-6F_{1,4} in Table. 2 ) as presented in Caron-Huot:2021usw which is similar to the four-loop one:

5(5)(u,v)+234(5)(u,v)4[27(5)(u,v)+27(5)(v,u)]=6F1,4=6F1F4+F2F3,-\mathcal{I}_{5}^{(5)}(u,v)+2\mathcal{I}_{34}^{(5)}(u,v)-4\left[\mathcal{I}_{27}^{(5)}(u,v)+\mathcal{I}_{27}^{(5)}(v,u)\right]=-6F_{1,4}=-6F_{1}F_{4}+F_{2}F_{3}, (73)

where 34(5)\mathcal{I}_{34}^{(5)} which only involves SVHPLs is known from HyperlogProcedures. This magic identity will give additional constraints when bootstrapping 27(5)\mathcal{I}_{27}^{(5)} and 5(5)\mathcal{I}_{5}^{(5)}. And it indicates that 27(5)\mathcal{I}_{27}^{(5)} and 5(5)\mathcal{I}_{5}^{(5)} may be expressed using only SVHPLs.

In the 23 integrals that can be computed using HyperlogProcedures, 14 equal to the five-loop ladder F5F_{5}. They satisfy the original magic identity discussed in Drummond:2006rz and are listed as follows.

9(5)=10(5)=11(5)=13(5)=15(5)=16(5)=17(5)=19(5)=20(5)=21(5)=22(5)=23(5)=31(5)=32(5).\displaystyle\mathcal{I}^{(5)}_{9}\!=\!\mathcal{I}^{(5)}_{10}\!=\!\mathcal{I}^{(5)}_{11}\!=\!\mathcal{I}^{(5)}_{13}\!=\!\mathcal{I}^{(5)}_{15}\!=\!\mathcal{I}^{(5)}_{16}\!=\!\mathcal{I}^{(5)}_{17}\!=\!\mathcal{I}^{(5)}_{19}\!=\!\mathcal{I}^{(5)}_{20}\!=\!\mathcal{I}^{(5)}_{21}\!=\!\mathcal{I}^{(5)}_{22}\!=\!\mathcal{I}^{(5)}_{23}\!=\!\mathcal{I}^{(5)}_{31}\!=\!\mathcal{I}^{(5)}_{32}. (74)
9(5)(uv,1v)=Iz,0,0,0,0,1,0,0,0,0,0,0Iz,0,0,0,0,0,1,0,0,0,0,0Δ.\displaystyle\mathcal{I}_{9}^{(5)}(\frac{u}{v},\frac{1}{v})=\frac{\mathrm{I}_{z,0,0,0,0,1,0,0,0,0,0,0}-\mathrm{I}_{z,0,0,0,0,0,1,0,0,0,0,0}}{\Delta}.

The other two identities are

24(5)(1v,uv)=25(5)(1v,uv)=26(5)(uv,1v)\displaystyle\mathcal{I}_{24}^{(5)}(\frac{1}{v},\frac{u}{v})=\mathcal{I}_{25}^{(5)}(\frac{1}{v},\frac{u}{v})=\mathcal{I}_{26}^{(5)}(\frac{u}{v},\frac{1}{v}) (75)
=70ζ7Iz,0,1,0,08ζ5Iz,0,1,0,0,1,0+20ζ5Iz,0,1,0,1,0,0Iz,0,1,0,0,0,1,0,0,1,0,0+Iz,0,1,0,0,1,0,0,0,1,0,0Δ,\displaystyle=\frac{70\zeta_{7}\mathrm{I}_{z,0,1,0,0}\!-\!8\zeta_{5}\mathrm{I}_{z,0,1,0,0,1,0}+20\zeta_{5}\mathrm{I}_{z,0,1,0,1,0,0}-\mathrm{I}_{z,0,1,0,0,0,1,0,0,1,0,0}+\mathrm{I}_{z,0,1,0,0,1,0,0,0,1,0,0}}{\Delta},

and

28(5)(1v,uv)=29(5)(1v,uv)=30(5)(1v,uv)\displaystyle\mathcal{I}_{28}^{(5)}(\frac{1}{v},\frac{u}{v})=\mathcal{I}_{29}^{(5)}(\frac{1}{v},\frac{u}{v})=\mathcal{I}_{30}^{(5)}(\frac{1}{v},\frac{u}{v}) (76)
=1Δ(42ζ7Iz,0,0,1,020ζ5Iz,0,0,1,0,0,0+12ζ5Iz,0,0,1,0,1,0+4ζ3Iz,0,0,1,0,0,1,0,0\displaystyle=\frac{1}{\Delta}\Big{(}-42\zeta_{7}\mathrm{I}_{z,0,0,1,0}-20\zeta_{5}\mathrm{I}_{z,0,0,1,0,0,0}+12\zeta_{5}\mathrm{I}_{z,0,0,1,0,1,0}+4\zeta_{3}\mathrm{I}_{z,0,0,1,0,0,1,0,0}
4ζ3Iz,0,0,1,0,1,0,0,0Iz,0,0,1,0,0,1,0,1,0,0,0+Iz,0,0,1,0,1,0,0,1,0,0,0).\displaystyle-4\zeta_{3}\mathrm{I}_{z,0,0,1,0,1,0,0,0}-\mathrm{I}_{z,0,0,1,0,0,1,0,1,0,0,0}+\mathrm{I}_{z,0,0,1,0,1,0,0,1,0,0,0}\Big{)}.

They take the form of generalized ladders Drummond:2012bg . The first identity can be proved by directly applying the magic identity to the two-loop subdiagrams. The second identity can be proved by performing twice the boxing for x3x_{3} and x4x_{4}. In this process, we have used the the same fact that any DCI integral is invariant under the exchange of legs, 14,231\leftrightarrow 4,2\leftrightarrow 3. Once the seeds are provided, they can “grow” in various ways by the inverse boxing. Therefore, these three identities have the same origin. In the remaining three integrals that can be computed by HyperlogProcedures, 8(5)\mathcal{I}_{8}^{(5)} can be related to f2\mathcal{I}^{f2} in the previous section through boxing and 33(5)\mathcal{I}_{33}^{(5)} can be related to some lower-loop DCI integral as well, but it is not any one that we discussed before. Interestingly, it is a purely lower-weight integral which takes the form:

33(5)(u,v)=20ζ5Δ(\displaystyle\mathcal{I}_{33}^{(5)}(u,v)=\frac{20\zeta_{5}}{\Delta}\Big{(} Iz,0,0,1,0,0Iz,0,0,1,1,0Iz,0,1,0,0,0+Iz,0,1,0,1,0Iz,1,0,1,0,0+Iz,1,0,1,1,0\displaystyle\mathrm{I}_{z,0,0,1,0,0}-\mathrm{I}_{z,0,0,1,1,0}-\mathrm{I}_{z,0,1,0,0,0}+\mathrm{I}_{z,0,1,0,1,0}-\mathrm{I}_{z,1,0,1,0,0}+\mathrm{I}_{z,1,0,1,1,0} (77)
+Iz,1,1,0,0,0Iz,1,1,0,1,0).\displaystyle+\mathrm{I}_{z,1,1,0,0,0}-\mathrm{I}_{z,1,1,0,1,0}\Big{)}.

It is a weigh-9 function. This agrees with the leading singularity analysis for the integrand, since the multivariate residue of the integrand is 0. At last, 34(5)\mathcal{I}_{34}^{(5)} is computed with a unique leading singularity 1/Δ/(uv)1/\Delta/(u-v).

Next, we study three unknown integrals, 2(5),4(5),6(5)\mathcal{I}_{2}^{(5)},\mathcal{I}_{4}^{(5)},\mathcal{I}_{6}^{(5)}, which can also be related to some lower-loop integrals through boxing. 2(5)\mathcal{I}_{2}^{(5)} and 4(5)\mathcal{I}_{4}^{(5)} are the same with each other. Because through the boxing they are related to two identical DCI integrals with legs 14,231\leftrightarrow 4,2\leftrightarrow 3 exchanged, which is f\mathcal{I}^{f} calculated in the previous section. Then 4(5)\mathcal{I}_{4}^{(5)}(2(5)\mathcal{I}_{2}^{(5)}) is obtained by the single-value integration and its boundary values provided in Chicherin:2018avq . 6(5)\mathcal{I}_{6}^{(5)} is related to a four-loop integral that we have not encountered before. It is depicted in Fig. 5.

Refer to caption
Figure 5: The four-loop DCI integral after acting boxing on 6(5)\mathcal{I}_{6}^{(5)}. It is not one of the four-loop integrals that originate from taking four cycles of ff-graphs.
integrals leading singularities (of weight-10)
1(5)\mathcal{I}_{1}^{(5)} 1u+vΔ2\frac{1-u+v}{\Delta^{2}}, 1Δ\frac{1}{\Delta},1(uv)Δ-\frac{1}{(u-v)\Delta}
3(5)\mathcal{I}_{3}^{(5)} 1u+vΔ2\frac{1-u+v}{\Delta^{2}},1Δ\frac{1}{\Delta}
5(5)\mathcal{I}_{5}^{(5)} 1Δ2\frac{1}{\Delta^{2}},1(uv)Δ-\frac{1}{(u-v)\Delta}
7(5)\mathcal{I}_{7}^{(5)} 1Δ\frac{1}{\Delta},1(uv)Δ-\frac{1}{(u-v)\Delta}
12(5)=18(5)\mathcal{I}_{12}^{(5)}=\mathcal{I}_{18}^{(5)} 1u+vΔ2\frac{1-u+v}{\Delta^{2}},1Δ\frac{1}{\Delta}
14(5)\mathcal{I}_{14}^{(5)} 1u+vΔ2\frac{1-u+v}{\Delta^{2}},1+uvΔ2\frac{1+u-v}{\Delta^{2}},1Δ\frac{1}{\Delta}
27(5)\mathcal{I}_{27}^{(5)} 1Δ2\frac{1}{\Delta^{2}}
33(5)\mathcal{I}_{33}^{(5)} -
34(5)\mathcal{I}_{34}^{(5)} 1(uv)Δ-\frac{1}{(u-v)\Delta}
the rest 1Δ\frac{1}{\Delta}
Table 5: The leading analysis for five-loop DCI integrals. For those integrals that are not of uniform weight, the leading singularities correspond to the coefficients of weight-10 part of the full results. It should also be noted that 33(5)\mathcal{I}_{33}^{(5)} is a purely lower-weight integral whose form is given in (77).

The result of 2(5)\mathcal{I}_{2}^{(5)} and 4(5)\mathcal{I}_{4}^{(5)} are provided as follows:

2(5)(1v,uv)=4(5)(1v,uv)=1Δ(4Iz,0,z¯,0,0,0,0,1,0,0,1,04Iz,0,1,0,0,1,0,0,0,0,1,0\displaystyle\mathcal{I}_{2}^{(5)}(\frac{1}{v},\frac{u}{v})=\mathcal{I}_{4}^{(5)}(\frac{1}{v},\frac{u}{v})=\frac{1}{\Delta}\Big{(}-4\mathrm{I}_{z,0,\bar{z},0,0,0,0,1,0,0,1,0}-4\mathrm{I}_{z,0,1,0,0,1,0,0,0,0,1,0} (78)
+4Iz,0,z¯,1,0,0,0,0,1,0,0,04Iz,0,z¯,1,0,0,1,0,0,0,0,04Iz,0,z¯,0,1,0,0,0,0,1,0,0+4Iz,0,z¯,0,0,0,1,0,0,1,0,0\displaystyle+4\mathrm{I}_{z,0,\bar{z},1,0,0,0,0,1,0,0,0}-4\mathrm{I}_{z,0,\bar{z},1,0,0,1,0,0,0,0,0}-4\mathrm{I}_{z,0,\bar{z},0,1,0,0,0,0,1,0,0}+4\mathrm{I}_{z,0,\bar{z},0,0,0,1,0,0,1,0,0}
+4Iz,0,z¯,0,1,0,0,1,0,0,0,0+Iz,0,0,0,0,1,0,0,1,0,0,0+3Iz,0,0,1,0,0,1,0,0,0,0,0+4Iz,0,z¯,0,0,1,0,0,0,0,1,0\displaystyle+4\mathrm{I}_{z,0,\bar{z},0,1,0,0,1,0,0,0,0}+\mathrm{I}_{z,0,0,0,0,1,0,0,1,0,0,0}+3\mathrm{I}_{z,0,0,1,0,0,1,0,0,0,0,0}+4\mathrm{I}_{z,0,\bar{z},0,0,1,0,0,0,0,1,0}
+4Iz,0,1,0,0,0,0,1,0,0,1,04Iz,0,z¯,0,0,1,0,0,1,0,0,02Iz,0,0,1,0,0,0,0,1,0,0,02Iz,0,0,0,1,0,0,1,0,0,0,0+\displaystyle+4\mathrm{I}_{z,0,1,0,0,0,0,1,0,0,1,0}-4\mathrm{I}_{z,0,\bar{z},0,0,1,0,0,1,0,0,0}-2\mathrm{I}_{z,0,0,1,0,0,0,0,1,0,0,0}-2\mathrm{I}_{z,0,0,0,1,0,0,1,0,0,0,0}+
56ζ7Iz,0,0,1,0Iz,0,0,0,0,0,1,0,1,0,0,02Iz,0,0,1,0,0,0,1,0,0,0,0+4ζ3(Iz,0,0,0,0,0,1,0,0Iz,0,0,0,0,1,0,0,0\displaystyle 56\zeta_{7}\mathrm{I}_{z,0,0,1,0}-\mathrm{I}_{z,0,0,0,0,0,1,0,1,0,0,0}-2\mathrm{I}_{z,0,0,1,0,0,0,1,0,0,0,0}+4\zeta_{3}\left(\mathrm{I}_{z,0,0,0,0,0,1,0,0}-\mathrm{I}_{z,0,0,0,0,1,0,0,0}\right.
+Iz,0,0,1,0,0,1,0,0Iz,0,0,1,0,1,0,0,0)+4ζ5(3Iz,0,0,0,0,1,010Iz,0,0,0,1,0,0+3Iz,0,0,1,0,1,0\displaystyle\left.+\!\mathrm{I}_{z,0,0,1,0,0,1,0,0}\!-\!\mathrm{I}_{z,0,0,1,0,1,0,0,0}\right)\!+\!4\zeta_{5}\!\left(3\mathrm{I}_{z,0,0,0,0,1,0}\!-\!10\mathrm{I}_{z,0,0,0,1,0,0}\!+\!3\mathrm{I}_{z,0,0,1,0,1,0}\!\right.
12Iz,0,1,0,0,1,0)Iz,0,0,1,0,0,1,0,1,0,0,0+2Iz,0,0,0,1,0,0,0,1,0,0,0+Iz,0,0,0,0,1,0,0,0,0,0,0\displaystyle\left.-12\mathrm{I}_{z,0,1,0,0,1,0}\right)-\mathrm{I}_{z,0,0,1,0,0,1,0,1,0,0,0}+2\mathrm{I}_{z,0,0,0,1,0,0,0,1,0,0,0}+\mathrm{I}_{z,0,0,0,0,1,0,0,0,0,0,0}
Iz,0,0,0,0,0,1,0,0,0,0,0+Iz,0,0,1,0,1,0,0,0,0,0,0+Iz,0,0,1,0,1,0,0,1,0,0,0)\displaystyle-\mathrm{I}_{z,0,0,0,0,0,1,0,0,0,0,0}+\mathrm{I}_{z,0,0,1,0,1,0,0,0,0,0,0}+\mathrm{I}_{z,0,0,1,0,1,0,0,1,0,0,0}\Big{)}
+2Δ(2Iz,1,0,0,0,0,1,0,0,1,02Iz,1,0,0,1,0,0,0,0,1,024ζ5Iz,1,0,0,1,0+(z+z¯)(Iz,0,0,0,0,1,0,0,1,0,0\displaystyle+\frac{2}{\Delta}\Bigg{(}2\mathrm{I}_{z,1,0,0,0,0,1,0,0,1,0}-2\mathrm{I}_{z,1,0,0,1,0,0,0,0,1,0}-24\zeta_{5}\mathrm{I}_{z,1,0,0,1,0}+(z+\bar{z})\Big{(}\mathrm{I}_{z,0,0,0,0,1,0,0,1,0,0}
Iz,0,0,0,1,0,0,1,0,0,0Iz,0,0,1,0,0,0,0,1,0,0+Iz,0,0,1,0,0,1,0,0,0,0+Iz,0,1,0,0,0,0,1,0,0,0\displaystyle-\mathrm{I}_{z,0,0,0,1,0,0,1,0,0,0}-\mathrm{I}_{z,0,0,1,0,0,0,0,1,0,0}+\mathrm{I}_{z,0,0,1,0,0,1,0,0,0,0}+\mathrm{I}_{z,0,1,0,0,0,0,1,0,0,0}
Iz,0,1,0,0,1,0,0,0,0,0Iz,1,0,0,0,0,1,0,0,1,0+Iz,1,0,0,1,0,0,0,0,1,0+12ζ5(Iz,1,0,0,1,0Iz,0,0,0,1,0)))\displaystyle-\mathrm{I}_{z,0,1,0,0,1,0,0,0,0,0}-\mathrm{I}_{z,1,0,0,0,0,1,0,0,1,0}+\mathrm{I}_{z,1,0,0,1,0,0,0,0,1,0}+12\zeta_{5}(\mathrm{I}_{z,1,0,0,1,0}-\mathrm{I}_{z,0,0,0,1,0})\Big{)}\Bigg{)}
+2(Iz,0,0,0,0,1,0,0,1,0,0Iz,0,0,0,1,0,0,1,0,0,0Iz,0,0,1,0,0,0,0,1,0,0+Iz,0,0,1,0,0,1,0,0,0,0\displaystyle+2\left(\mathrm{I}_{z,0,0,0,0,1,0,0,1,0,0}-\mathrm{I}_{z,0,0,0,1,0,0,1,0,0,0}-\mathrm{I}_{z,0,0,1,0,0,0,0,1,0,0}+\mathrm{I}_{z,0,0,1,0,0,1,0,0,0,0}\right.
+Iz,0,1,0,0,0,0,1,0,0,0Iz,0,1,0,0,1,0,0,0,0,0Iz,1,0,0,0,0,1,0,0,1,0+Iz,1,0,0,1,0,0,0,0,1,0\displaystyle+\mathrm{I}_{z,0,1,0,0,0,0,1,0,0,0}-\mathrm{I}_{z,0,1,0,0,1,0,0,0,0,0}-\mathrm{I}_{z,1,0,0,0,0,1,0,0,1,0}+\mathrm{I}_{z,1,0,0,1,0,0,0,0,1,0}
+2Iz,z¯,0,0,0,0,1,0,0,1,02Iz,z¯,0,0,0,1,0,0,1,0,02Iz,z¯,0,0,1,0,0,0,0,1,0+2Iz,z¯,0,0,1,0,0,1,0,0,0\displaystyle+2\mathrm{I}_{z,\bar{z},0,0,0,0,1,0,0,1,0}-2\mathrm{I}_{z,\bar{z},0,0,0,1,0,0,1,0,0}-2\mathrm{I}_{z,\bar{z},0,0,1,0,0,0,0,1,0}+2\mathrm{I}_{z,\bar{z},0,0,1,0,0,1,0,0,0}
+2Iz,z¯,0,1,0,0,0,0,1,0,02Iz,z¯,0,1,0,0,1,0,0,0,02Iz,z¯,1,0,0,0,0,1,0,0,0+2Iz,z¯,1,0,0,1,0,0,0,0,0)\displaystyle\left.+2\mathrm{I}_{z,\bar{z},0,1,0,0,0,0,1,0,0}-2\mathrm{I}_{z,\bar{z},0,1,0,0,1,0,0,0,0}-2\mathrm{I}_{z,\bar{z},1,0,0,0,0,1,0,0,0}+2\mathrm{I}_{z,\bar{z},1,0,0,1,0,0,0,0,0}\right)
+24ζ5(Iz,1,0,0,1,0Iz,0,0,0,1,0)+168ζ9.\displaystyle+24\zeta_{5}(\mathrm{I}_{z,1,0,0,1,0}-\mathrm{I}_{z,0,0,0,1,0})+168\zeta_{9}.

The first part is the weight-10 part of the result and the remaining part is lower-weight (weight-9) terms.

Then we perform the leading singularity analysis for all 34 five-loop DCI integrals and list the results in Table. 5. It should be noted that, for those integrals that are not of uniform weight, the leading singularities given only correspond to the weight-10 part of the full results. Among all the unknown integrals, 12(5)=18(5)\mathcal{I}_{12}^{(5)}=\mathcal{I}_{18}^{(5)} due to the magic identity existing in the two-loop sub-diagrams.

For completeness, we record analytic results for all integrals we computed/bootstrapped (including most interesting ones like 2(5)\mathcal{I}_{2}^{(5)}, 4(5)\mathcal{I}_{4}^{(5)} and 2(4)=f{\cal I}_{2}^{(4)}={\cal I}^{f}) in an ancillary file intresults.nb along with this paper.

5 Discussions and Outlook

In this paper, we have continued explorations for multi-loop integrands and integrated results for four-point correlators and amplitudes, and especially those DCI integrals contributing to Coulomb-branch amplitudes in the ten-dimensional lightlike limit Caron-Huot:2021usw . At integrand level, we identified graphical structures in ff-graphs which contribute to component Ma,bM_{a,b}, which in turn provides all-loop predictions for coefficients of some (4+)(4+\ell)-point ff-graphs containing (a,b)(a,b)-structure with ab1ab\geq\ell-1. By comparing the Coulomb-branch amplitudes with octagons at integrated level, we obtain three infinite families of magic identities, which originates from those prototypes for M1,2(4),M1,2(5)M_{1,2}^{(4)},M_{1,2}^{(5)} and M1,3(6)M_{1,3}^{(6)}. For individual DCI integrals, we have bootstrapped the 33 non-trivial ones at four loops, in terms of nice weight-88 SVHPL functions; at five loops, 2525 of 3434 DCI integrals can be computed either by using HyperlogProcedures or by solving “boxing” differential equations, and we have analyzed all leading singularities and obtained functions beyond SVHPL at both weight 1010 and 99. Our preliminary studies have opened up several new avenues for future investigations.

We have obtained some constraints on the coefficients of ff-graphs from integrability, and it would be interesting to understand them by using the pinching limit He:2024cej (and other graphical constraints Bourjaily:2016evz ), and even to improve the efficiency of graphical bootstrap from these inputs. In general, it would be highly desirable to see whether the octagon could tell us more about higher-loop ff-graphs. Perhaps we could first understand the rule for coefficients of next-to-fishnet, which seems to be tip of an iceberg of such constructive rules for certain ff-graphs. Besides, all the discussions here are restricted to four-point Coulomb branch amplitudes, but it would be very interesting if one could produce integrands of higher-point amplitudes with more general 10d lightlike limits.

For integrated results, it would be highly desirable if one could systematically classify these (infinite families of) “irreducible” magic identities provided by octagons at higher loops: is there some all-loop structure to be discovered here, and can one understand them from the point of view of conformal/Feynman integrals e.g. via some representations such as Mellin-Barnes? Any progress in this direction may shed new lights into mathematical structures and relations for multi-loop DCI integrals. Note that these Coulomb branch amplitudes from 10d lightlike limit differ from those considered in e.g. Arkani-Hamed:2023epq ; Flieger:2025ekn , which have been partly motivated by geometric structures, and it is tempting to ask if they are related. Also it would be interesting to explore all-loop infrared divergences when we take the massless limit as discussed in Caron-Huot:2021usw . Another important direction concerns periods of ff-graphs, which are strongly constrained by magic identities as well as exact results from integrated correlators etc.; in particular, an interesting open question is if there exists closed-formula for periods associated with general determinants of ladders from octagons.

We have illustrated the power of bootstrap method for both four and five-loop DCI integrals, based on their leading singularities and function space of SVHPL (or its extensions). It would be interesting to reorganize the integrands such that individual integrals all have uniform, maximal transcendentality (weight 22\ell), which would greatly simplify such bootstrap program for higher loop integrals. Knowing the analytic results for DCI integrals at four and higher loops would be very useful for computing other amplitudes or correlators that can be expressed in terms of these integrals, and we leave these potential applications of our results to future investigations.

Recall that our original motivation was to study all four-loop conformal integrals, extending the classic results in Drummond:2013nda at =3\ell=3. Indeed by looking at the 3030 conformal integrals at =4\ell=4 (excluding the two containing elliptic pieces), we have bootstrapped most of them in terms of SVHPL or SVMPL functions, but there are still a few integrals that we cannot fix due to lack of boundary data. It would be highly desirable to determine these 3030 integrals, which can be applied to the computation of other physical quantities. Last but not least, although we have focused on the general goal of bootstrapping individual conformal integrals, a similar bootstrap program can also be set up for the four-point correlator without solving all individual integrals, where e.g. through four loops we do know all possible leading singularities and accompanying function spaces, but more physical constraints/boundary data are needed (perhaps from integrability or other inputs). It would be extremely interesting to study this problem in the future.

gbsn

Acknowledgements.
It is our pleasure to thank Jacob Bourjaily, Yu-tin Huang, Chia-Kai Kuo, Canxin Shi, Yichao Tang, Qinglin Yang for collaborations on related topics and/or helpful discussions. The work of SH is supported by the National Natural Science Foundation of China under Grant No. 12225510, 11935013, 12047503, 12247103, and by the New Cornerstone Science Foundation through the XPLORER PRIZE.

Appendix A The definition of four and five-loop DCI integrals

Here we present the definition of DCI integrals at =4,5\ell=4,5 used in the main text . They have all been normalized by x132x242x_{13}^{2}x_{24}^{2}.

1(4)\displaystyle\mathcal{I}^{(4)}_{1} =[Uncaptioned image]\displaystyle=\includegraphics[scale={0.4},align=c]{graph/DCI4l1.pdf} =[Uncaptioned image]\displaystyle=\includegraphics[scale={0.4},align=c]{graph/posi4l1.pdf} =a=58d4xax132x182x244x472x162x172x272x282x382x452x462x482x562x572x582x672x782\displaystyle=\int\prod_{a=5}^{8}d^{4}x_{a}\frac{x_{13}^{2}x_{18}^{2}x_{24}^{4}x_{47}^{2}}{x_{16}^{2}x_{17}^{2}x_{27}^{2}x_{28}^{2}x_{38}^{2}x_{45}^{2}x_{46}^{2}x_{48}^{2}x_{56}^{2}x_{57}^{2}x_{58}^{2}x_{67}^{2}x_{78}^{2}}
2(4)\displaystyle\mathcal{I}^{(4)}_{2} =[Uncaptioned image]\displaystyle=\includegraphics[scale={0.4},align=c]{graph/DCI4l2.pdf} =[Uncaptioned image]\displaystyle=\includegraphics[scale={0.4},align=c]{graph/posi4l2.pdf} =a=58d4xax132x182x244x372x162x172x262x282x352x382x452x472x572x582x672x682x782\displaystyle=\int\prod_{a=5}^{8}d^{4}x_{a}\frac{x_{13}^{2}x_{18}^{2}x_{24}^{4}x_{37}^{2}}{x_{16}^{2}x_{17}^{2}x_{26}^{2}x_{28}^{2}x_{35}^{2}x_{38}^{2}x_{45}^{2}x_{47}^{2}x_{57}^{2}x_{58}^{2}x_{67}^{2}x_{68}^{2}x_{78}^{2}}
3(4)\displaystyle\mathcal{I}^{(4)}_{3} =[Uncaptioned image]\displaystyle=\includegraphics[scale={0.4},align=c]{graph/DCI4l3.pdf} =[Uncaptioned image]\displaystyle=\includegraphics[scale={0.4},align=c]{graph/posi4l3.pdf} =a=58d4xax132x246x782x182x262x272x282x372x452x472x482x562x572x582x672x682\displaystyle=\int\prod_{a=5}^{8}d^{4}x_{a}\frac{x_{13}^{2}x_{24}^{6}x_{78}^{2}}{x_{18}^{2}x_{26}^{2}x_{27}^{2}x_{28}^{2}x_{37}^{2}x_{45}^{2}x_{47}^{2}x_{48}^{2}x_{56}^{2}x_{57}^{2}x_{58}^{2}x_{67}^{2}x_{68}^{2}}
4(4)\displaystyle\mathcal{I}^{(4)}_{4} =[Uncaptioned image]\displaystyle=\includegraphics[scale={0.4},align=c]{graph/DCI4l4.pdf} =[Uncaptioned image]\displaystyle=\includegraphics[scale={0.4},align=c]{graph/posi4l4.pdf} =a=58d4xax132x246x382x172x262x272x282x352x362x452x472x482x562x582x682x782\displaystyle=\int\prod_{a=5}^{8}d^{4}x_{a}\frac{x_{13}^{2}x_{24}^{6}x_{38}^{2}}{x_{17}^{2}x_{26}^{2}x_{27}^{2}x_{28}^{2}x_{35}^{2}x_{36}^{2}x_{45}^{2}x_{47}^{2}x_{48}^{2}x_{56}^{2}x_{58}^{2}x_{68}^{2}x_{78}^{2}}
5(4)\displaystyle\mathcal{I}^{(4)}_{5} =[Uncaptioned image]\displaystyle=\includegraphics[scale={0.4},align=c]{graph/DCI4l5.pdf} =[Uncaptioned image]\displaystyle=\includegraphics[scale={0.4},align=c]{graph/posi4l5.pdf} =a=58d4xax132x244x152x272x282x372x452x462x562x582x672x682x782\displaystyle=\int\prod_{a=5}^{8}d^{4}x_{a}\frac{x_{13}^{2}x_{24}^{4}}{x_{15}^{2}x_{27}^{2}x_{28}^{2}x_{37}^{2}x_{45}^{2}x_{46}^{2}x_{56}^{2}x_{58}^{2}x_{67}^{2}x_{68}^{2}x_{78}^{2}}
6(4)\displaystyle\mathcal{I}^{(4)}_{6} =[Uncaptioned image]\displaystyle=\includegraphics[scale={0.4},align=c]{graph/DCI4l6.pdf} =[Uncaptioned image]\displaystyle=\includegraphics[scale={0.4},align=c]{graph/posi4l6.pdf} =a=58d4xax134x244x172x182x262x272x352x362x452x482x562x582x672x782\displaystyle=\int\prod_{a=5}^{8}d^{4}x_{a}\frac{x_{13}^{4}x_{24}^{4}}{x_{17}^{2}x_{18}^{2}x_{26}^{2}x_{27}^{2}x_{35}^{2}x_{36}^{2}x_{45}^{2}x_{48}^{2}x_{56}^{2}x_{58}^{2}x_{67}^{2}x_{78}^{2}}
7(4)\displaystyle\mathcal{I}^{(4)}_{7} =[Uncaptioned image]\displaystyle=\includegraphics[scale={0.4},align=c]{graph/DCI4l7.pdf} =[Uncaptioned image]\displaystyle=\includegraphics[scale={0.4},align=c]{graph/posi4l7.pdf} =a=58d4xax132x184x244x152x162x172x272x282x382x452x482x562x582x672x682x782\displaystyle=\int\prod_{a=5}^{8}d^{4}x_{a}\frac{x_{13}^{2}x_{18}^{4}x_{24}^{4}}{x_{15}^{2}x_{16}^{2}x_{17}^{2}x_{27}^{2}x_{28}^{2}x_{38}^{2}x_{45}^{2}x_{48}^{2}x_{56}^{2}x_{58}^{2}x_{67}^{2}x_{68}^{2}x_{78}^{2}}
8(4)\displaystyle\mathcal{I}^{(4)}_{8} =[Uncaptioned image]\displaystyle=\includegraphics[scale={0.4},align=c]{graph/DCI4l8.pdf} =[Uncaptioned image]\displaystyle=\includegraphics[scale={0.4},align=c]{graph/posi4l8.pdf} =a=58d4xax132x248x182x252x262x272x282x372x452x462x472x482x562x572x682\displaystyle=\int\prod_{a=5}^{8}d^{4}x_{a}\frac{x_{13}^{2}x_{24}^{8}}{x_{18}^{2}x_{25}^{2}x_{26}^{2}x_{27}^{2}x_{28}^{2}x_{37}^{2}x_{45}^{2}x_{46}^{2}x_{47}^{2}x_{48}^{2}x_{56}^{2}x_{57}^{2}x_{68}^{2}}

These 88 DCI integrals at =4\ell=4 were denoted as (e),(d2),(d),(c),(f),(f2),(b),(a)(e),(d2),(d),(c),(f),(f2),(b),(a) in Figure 14 of Caron-Huot:2021usw . Note that the original magic identity implies 1(4)=3(4)=4(4)=7(4)=8(4)\mathcal{I}^{(4)}_{1}=\mathcal{I}^{(4)}_{3}=\mathcal{I}^{(4)}_{4}=\mathcal{I}^{(4)}_{7}=\mathcal{I}^{(4)}_{8}, while the new magic identity involve 2,5,6{\cal I}_{2},{\cal I}_{5},{\cal I}_{6} as we discussed in eq (57).

1(5)\displaystyle\mathcal{I}^{(5)}_{1} =[Uncaptioned image]\displaystyle=\includegraphics[scale={0.4},align=c]{graph/DCI5l1.pdf} =[Uncaptioned image]\displaystyle=\includegraphics[scale={0.4},align=c]{graph/posi5l1.pdf} =a=59d4xax132x182x244x392x672x162x192x272x292x372x382x462x482x562x572x582x592x682x692x782x792\displaystyle=\int\prod_{a=5}^{9}d^{4}x_{a}\frac{x_{13}^{2}x_{18}^{2}x_{24}^{4}x_{39}^{2}x_{67}^{2}}{x_{16}^{2}x_{19}^{2}x_{27}^{2}x_{29}^{2}x_{37}^{2}x_{38}^{2}x_{46}^{2}x_{48}^{2}x_{56}^{2}x_{57}^{2}x_{58}^{2}x_{59}^{2}x_{68}^{2}x_{69}^{2}x_{78}^{2}x_{79}^{2}}
2(5)\displaystyle\mathcal{I}^{(5)}_{2} =[Uncaptioned image]\displaystyle=\includegraphics[scale={0.4},align=c]{graph/DCI5l2.pdf} =[Uncaptioned image]\displaystyle=\includegraphics[scale={0.4},align=c]{graph/posi5l2.pdf} =a=59d4xax132x246x372x892x192x262x282x292x352x382x452x472x492x572x582x672x682x692x782x792\displaystyle=\int\prod_{a=5}^{9}d^{4}x_{a}\frac{x_{13}^{2}x_{24}^{6}x_{37}^{2}x_{89}^{2}}{x_{19}^{2}x_{26}^{2}x_{28}^{2}x_{29}^{2}x_{35}^{2}x_{38}^{2}x_{45}^{2}x_{47}^{2}x_{49}^{2}x_{57}^{2}x_{58}^{2}x_{67}^{2}x_{68}^{2}x_{69}^{2}x_{78}^{2}x_{79}^{2}}
3(5)\displaystyle\mathcal{I}^{(5)}_{3} =[Uncaptioned image]\displaystyle=\includegraphics[scale={0.4},align=c]{graph/DCI5l3.pdf} =[Uncaptioned image]\displaystyle=\includegraphics[scale={0.4},align=c]{graph/posi5l3.pdf} =a=59d4xax132x172x244x382x492x152x182x282x292x362x392x452x462x472x572x582x672x692x782x792x892\displaystyle=\int\prod_{a=5}^{9}d^{4}x_{a}\frac{x_{13}^{2}x_{17}^{2}x_{24}^{4}x_{38}^{2}x_{49}^{2}}{x_{15}^{2}x_{18}^{2}x_{28}^{2}x_{29}^{2}x_{36}^{2}x_{39}^{2}x_{45}^{2}x_{46}^{2}x_{47}^{2}x_{57}^{2}x_{58}^{2}x_{67}^{2}x_{69}^{2}x_{78}^{2}x_{79}^{2}x_{89}^{2}}
4(5)\displaystyle\mathcal{I}^{(5)}_{4} =[Uncaptioned image]\displaystyle=\includegraphics[scale={0.4},align=c]{graph/DCI5l4.pdf} =[Uncaptioned image]\displaystyle=\includegraphics[scale={0.4},align=c]{graph/posi5l4.pdf} =a=59d4xax132x244x272x392x482x192x252x282x292x362x382x462x472x492x572x582x592x672x682x782x792\displaystyle=\int\prod_{a=5}^{9}d^{4}x_{a}\frac{x_{13}^{2}x_{24}^{4}x_{27}^{2}x_{39}^{2}x_{48}^{2}}{x_{19}^{2}x_{25}^{2}x_{28}^{2}x_{29}^{2}x_{36}^{2}x_{38}^{2}x_{46}^{2}x_{47}^{2}x_{49}^{2}x_{57}^{2}x_{58}^{2}x_{59}^{2}x_{67}^{2}x_{68}^{2}x_{78}^{2}x_{79}^{2}}
5(5)\displaystyle\mathcal{I}^{(5)}_{5} =[Uncaptioned image]\displaystyle=\includegraphics[scale={0.4},align=c]{graph/DCI5l5.pdf} =[Uncaptioned image]\displaystyle=\includegraphics[scale={0.4},align=c]{graph/posi5l5.pdf} =a=59d4xax134x244x682x792x182x192x272x282x362x372x462x492x562x572x582x592x672x692x782x892\displaystyle=\int\prod_{a=5}^{9}d^{4}x_{a}\frac{x_{13}^{4}x_{24}^{4}x_{68}^{2}x_{79}^{2}}{x_{18}^{2}x_{19}^{2}x_{27}^{2}x_{28}^{2}x_{36}^{2}x_{37}^{2}x_{46}^{2}x_{49}^{2}x_{56}^{2}x_{57}^{2}x_{58}^{2}x_{59}^{2}x_{67}^{2}x_{69}^{2}x_{78}^{2}x_{89}^{2}}
6(5)\displaystyle\mathcal{I}^{(5)}_{6} =[Uncaptioned image]\displaystyle=\includegraphics[scale={0.4},align=c]{graph/DCI5l6.pdf} =[Uncaptioned image]\displaystyle=\includegraphics[scale={0.4},align=c]{graph/posi5l6.pdf} =a=59d4xax132x182x244x492x172x192x262x292x352x452x472x482x562x582x682x692x782x792x892\displaystyle=\int\prod_{a=5}^{9}d^{4}x_{a}\frac{x_{13}^{2}x_{18}^{2}x_{24}^{4}x_{49}^{2}}{x_{17}^{2}x_{19}^{2}x_{26}^{2}x_{29}^{2}x_{35}^{2}x_{45}^{2}x_{47}^{2}x_{48}^{2}x_{56}^{2}x_{58}^{2}x_{68}^{2}x_{69}^{2}x_{78}^{2}x_{79}^{2}x_{89}^{2}}
7(5)\displaystyle\mathcal{I}^{(5)}_{7} =[Uncaptioned image]\displaystyle=\includegraphics[scale={0.4},align=c]{graph/DCI5l7.pdf} =[Uncaptioned image]\displaystyle=\includegraphics[scale={0.4},align=c]{graph/posi5l7.pdf} =a=59d4xax132x192x244x382x172x182x262x272x362x392x482x492x562x572x582x592x692x782x892\displaystyle=\int\prod_{a=5}^{9}d^{4}x_{a}\frac{x_{13}^{2}x_{19}^{2}x_{24}^{4}x_{38}^{2}}{x_{17}^{2}x_{18}^{2}x_{26}^{2}x_{27}^{2}x_{36}^{2}x_{39}^{2}x_{48}^{2}x_{49}^{2}x_{56}^{2}x_{57}^{2}x_{58}^{2}x_{59}^{2}x_{69}^{2}x_{78}^{2}x_{89}^{2}}
8(5)\displaystyle\mathcal{I}^{(5)}_{8} =[Uncaptioned image]\displaystyle=\includegraphics[scale={0.4},align=c]{graph/DCI5l8.pdf} =[Uncaptioned image]\displaystyle=\includegraphics[scale={0.4},align=c]{graph/posi5l8.pdf} =a=59d4xax132x246x392x192x262x282x292x362x372x452x472x492x572x582x592x672x682x892\displaystyle=\int\prod_{a=5}^{9}d^{4}x_{a}\frac{x_{13}^{2}x_{24}^{6}x_{39}^{2}}{x_{19}^{2}x_{26}^{2}x_{28}^{2}x_{29}^{2}x_{36}^{2}x_{37}^{2}x_{45}^{2}x_{47}^{2}x_{49}^{2}x_{57}^{2}x_{58}^{2}x_{59}^{2}x_{67}^{2}x_{68}^{2}x_{89}^{2}}
9(5)\displaystyle\mathcal{I}^{(5)}_{9} =[Uncaptioned image]\displaystyle=\includegraphics[scale={0.4},align=c]{graph/DCI5l9.pdf} =[Uncaptioned image]\displaystyle=\includegraphics[scale={0.4},align=c]{graph/posi5l9.pdf} =a=59d4xax132x244x392x482x792x192x282x292x372x382x452x472x492x562x572x592x672x682x692x782x892\displaystyle=\int\prod_{a=5}^{9}d^{4}x_{a}\frac{x_{13}^{2}x_{24}^{4}x_{39}^{2}x_{48}^{2}x_{79}^{2}}{x_{19}^{2}x_{28}^{2}x_{29}^{2}x_{37}^{2}x_{38}^{2}x_{45}^{2}x_{47}^{2}x_{49}^{2}x_{56}^{2}x_{57}^{2}x_{59}^{2}x_{67}^{2}x_{68}^{2}x_{69}^{2}x_{78}^{2}x_{89}^{2}}
10(5)\displaystyle\mathcal{I}^{(5)}_{10} =[Uncaptioned image]\displaystyle=\includegraphics[scale={0.4},align=c]{graph/DCI5l10.pdf} =[Uncaptioned image]\displaystyle=\includegraphics[scale={0.4},align=c]{graph/posi5l10.pdf} =a=59d4xax132x172x182x244x492x152x162x192x282x292x382x462x472x482x562x572x592x672x782x792x892\displaystyle=\int\prod_{a=5}^{9}d^{4}x_{a}\frac{x_{13}^{2}x_{17}^{2}x_{18}^{2}x_{24}^{4}x_{49}^{2}}{x_{15}^{2}x_{16}^{2}x_{19}^{2}x_{28}^{2}x_{29}^{2}x_{38}^{2}x_{46}^{2}x_{47}^{2}x_{48}^{2}x_{56}^{2}x_{57}^{2}x_{59}^{2}x_{67}^{2}x_{78}^{2}x_{79}^{2}x_{89}^{2}}
11(5)\displaystyle\mathcal{I}^{(5)}_{11} =[Uncaptioned image]\displaystyle=\includegraphics[scale={0.4},align=c]{graph/DCI5l11.pdf} =[Uncaptioned image]\displaystyle=\includegraphics[scale={0.4},align=c]{graph/posi5l11.pdf} =a=59d4xax132x246x472x892x192x272x282x292x382x452x462x482x492x562x572x592x672x682x782x792\displaystyle=\int\prod_{a=5}^{9}d^{4}x_{a}\frac{x_{13}^{2}x_{24}^{6}x_{47}^{2}x_{89}^{2}}{x_{19}^{2}x_{27}^{2}x_{28}^{2}x_{29}^{2}x_{38}^{2}x_{45}^{2}x_{46}^{2}x_{48}^{2}x_{49}^{2}x_{56}^{2}x_{57}^{2}x_{59}^{2}x_{67}^{2}x_{68}^{2}x_{78}^{2}x_{79}^{2}}
12(5)\displaystyle\mathcal{I}^{(5)}_{12} =[Uncaptioned image]\displaystyle=\includegraphics[scale={0.4},align=c]{graph/DCI5l12.pdf} =[Uncaptioned image]\displaystyle=\includegraphics[scale={0.4},align=c]{graph/posi5l12.pdf} =a=59d4xax132x182x244x392x492x152x192x272x292x372x382x452x462x482x562x592x682x692x782x792x892\displaystyle=\int\prod_{a=5}^{9}d^{4}x_{a}\frac{x_{13}^{2}x_{18}^{2}x_{24}^{4}x_{39}^{2}x_{49}^{2}}{x_{15}^{2}x_{19}^{2}x_{27}^{2}x_{29}^{2}x_{37}^{2}x_{38}^{2}x_{45}^{2}x_{46}^{2}x_{48}^{2}x_{56}^{2}x_{59}^{2}x_{68}^{2}x_{69}^{2}x_{78}^{2}x_{79}^{2}x_{89}^{2}}
13(5)\displaystyle\mathcal{I}^{(5)}_{13} =[Uncaptioned image]\displaystyle=\includegraphics[scale={0.4},align=c]{graph/DCI5l13.pdf} =[Uncaptioned image]\displaystyle=\includegraphics[scale={0.4},align=c]{graph/posi5l13.pdf} =a=59d4xax132x246x392x482x172x272x282x292x352x382x452x462x472x492x562x582x682x692x792x892\displaystyle=\int\prod_{a=5}^{9}d^{4}x_{a}\frac{x_{13}^{2}x_{24}^{6}x_{39}^{2}x_{48}^{2}}{x_{17}^{2}x_{27}^{2}x_{28}^{2}x_{29}^{2}x_{35}^{2}x_{38}^{2}x_{45}^{2}x_{46}^{2}x_{47}^{2}x_{49}^{2}x_{56}^{2}x_{58}^{2}x_{68}^{2}x_{69}^{2}x_{79}^{2}x_{89}^{2}}
14(5)\displaystyle\mathcal{I}^{(5)}_{14} =[Uncaptioned image]\displaystyle=\includegraphics[scale={0.4},align=c]{graph/DCI5l14.pdf} =[Uncaptioned image]\displaystyle=\includegraphics[scale={0.4},align=c]{graph/posi5l14.pdf} =a=59d4xax132x172x182x244x392x152x162x192x262x282x372x382x452x472x572x592x682x692x782x792x892\displaystyle=\int\prod_{a=5}^{9}d^{4}x_{a}\frac{x_{13}^{2}x_{17}^{2}x_{18}^{2}x_{24}^{4}x_{39}^{2}}{x_{15}^{2}x_{16}^{2}x_{19}^{2}x_{26}^{2}x_{28}^{2}x_{37}^{2}x_{38}^{2}x_{45}^{2}x_{47}^{2}x_{57}^{2}x_{59}^{2}x_{68}^{2}x_{69}^{2}x_{78}^{2}x_{79}^{2}x_{89}^{2}}
15(5)\displaystyle\mathcal{I}^{(5)}_{15} =[Uncaptioned image]\displaystyle=\includegraphics[scale={0.4},align=c]{graph/DCI5l15.pdf} =[Uncaptioned image]\displaystyle=\includegraphics[scale={0.4},align=c]{graph/posi5l15.pdf} =a=59d4xax132x192x246x392x162x172x272x282x292x352x382x452x462x492x582x592x672x692x792x892\displaystyle=\int\prod_{a=5}^{9}d^{4}x_{a}\frac{x_{13}^{2}x_{19}^{2}x_{24}^{6}x_{39}^{2}}{x_{16}^{2}x_{17}^{2}x_{27}^{2}x_{28}^{2}x_{29}^{2}x_{35}^{2}x_{38}^{2}x_{45}^{2}x_{46}^{2}x_{49}^{2}x_{58}^{2}x_{59}^{2}x_{67}^{2}x_{69}^{2}x_{79}^{2}x_{89}^{2}}
16(5)\displaystyle\mathcal{I}^{(5)}_{16} =[Uncaptioned image]\displaystyle=\includegraphics[scale={0.4},align=c]{graph/DCI5l16.pdf} =[Uncaptioned image]\displaystyle=\includegraphics[scale={0.4},align=c]{graph/posi5l16.pdf} =a=59d4xax132x244x382x494x182x282x292x372x392x452x462x472x482x562x572x592x682x692x792x892\displaystyle=\int\prod_{a=5}^{9}d^{4}x_{a}\frac{x_{13}^{2}x_{24}^{4}x_{38}^{2}x_{49}^{4}}{x_{18}^{2}x_{28}^{2}x_{29}^{2}x_{37}^{2}x_{39}^{2}x_{45}^{2}x_{46}^{2}x_{47}^{2}x_{48}^{2}x_{56}^{2}x_{57}^{2}x_{59}^{2}x_{68}^{2}x_{69}^{2}x_{79}^{2}x_{89}^{2}}
17(5)\displaystyle\mathcal{I}^{(5)}_{17} =[Uncaptioned image]\displaystyle=\includegraphics[scale={0.4},align=c]{graph/DCI5l17.pdf} =[Uncaptioned image]\displaystyle=\includegraphics[scale={0.4},align=c]{graph/posi5l17.pdf} =a=59d4xax132x246x894x192x272x282x292x382x462x482x492x562x572x582x592x682x692x782x792\displaystyle=\int\prod_{a=5}^{9}d^{4}x_{a}\frac{x_{13}^{2}x_{24}^{6}x_{89}^{4}}{x_{19}^{2}x_{27}^{2}x_{28}^{2}x_{29}^{2}x_{38}^{2}x_{46}^{2}x_{48}^{2}x_{49}^{2}x_{56}^{2}x_{57}^{2}x_{58}^{2}x_{59}^{2}x_{68}^{2}x_{69}^{2}x_{78}^{2}x_{79}^{2}}
18(5)\displaystyle\mathcal{I}^{(5)}_{18} =[Uncaptioned image]\displaystyle=\includegraphics[scale={0.4},align=c]{graph/DCI5l18.pdf} =[Uncaptioned image]\displaystyle=\includegraphics[scale={0.4},align=c]{graph/posi5l18.pdf} =a=59d4xax132x194x244x382x162x172x182x252x282x352x392x462x492x582x592x672x692x782x792x892\displaystyle=\int\prod_{a=5}^{9}d^{4}x_{a}\frac{x_{13}^{2}x_{19}^{4}x_{24}^{4}x_{38}^{2}}{x_{16}^{2}x_{17}^{2}x_{18}^{2}x_{25}^{2}x_{28}^{2}x_{35}^{2}x_{39}^{2}x_{46}^{2}x_{49}^{2}x_{58}^{2}x_{59}^{2}x_{67}^{2}x_{69}^{2}x_{78}^{2}x_{79}^{2}x_{89}^{2}}
19(5)\displaystyle\mathcal{I}^{(5)}_{19} =[Uncaptioned image]\displaystyle=\includegraphics[scale={0.4},align=c]{graph/DCI5l19.pdf} =[Uncaptioned image]\displaystyle=\includegraphics[scale={0.4},align=c]{graph/posi5l19.pdf} =a=59d4xax132x246x394x152x252x282x292x362x372x382x452x472x492x592x672x682x692x792x892\displaystyle=\int\prod_{a=5}^{9}d^{4}x_{a}\frac{x_{13}^{2}x_{24}^{6}x_{39}^{4}}{x_{15}^{2}x_{25}^{2}x_{28}^{2}x_{29}^{2}x_{36}^{2}x_{37}^{2}x_{38}^{2}x_{45}^{2}x_{47}^{2}x_{49}^{2}x_{59}^{2}x_{67}^{2}x_{68}^{2}x_{69}^{2}x_{79}^{2}x_{89}^{2}}
20(5)\displaystyle\mathcal{I}^{(5)}_{20} =[Uncaptioned image]\displaystyle=\includegraphics[scale={0.4},align=c]{graph/DCI5l20.pdf} =[Uncaptioned image]\displaystyle=\includegraphics[scale={0.4},align=c]{graph/posi5l20.pdf} =a=59d4xax132x194x244x482x162x172x182x272x292x392x452x462x492x562x582x592x682x782x792x892\displaystyle=\int\prod_{a=5}^{9}d^{4}x_{a}\frac{x_{13}^{2}x_{19}^{4}x_{24}^{4}x_{48}^{2}}{x_{16}^{2}x_{17}^{2}x_{18}^{2}x_{27}^{2}x_{29}^{2}x_{39}^{2}x_{45}^{2}x_{46}^{2}x_{49}^{2}x_{56}^{2}x_{58}^{2}x_{59}^{2}x_{68}^{2}x_{78}^{2}x_{79}^{2}x_{89}^{2}}
21(5)\displaystyle\mathcal{I}^{(5)}_{21} =[Uncaptioned image]\displaystyle=\includegraphics[scale={0.4},align=c]{graph/DCI5l21.pdf} =[Uncaptioned image]\displaystyle=\includegraphics[scale={0.4},align=c]{graph/posi5l21.pdf} =a=59d4xax132x244x394x782x192x282x292x362x372x382x472x492x562x572x582x592x672x682x792x892\displaystyle=\int\prod_{a=5}^{9}d^{4}x_{a}\frac{x_{13}^{2}x_{24}^{4}x_{39}^{4}x_{78}^{2}}{x_{19}^{2}x_{28}^{2}x_{29}^{2}x_{36}^{2}x_{37}^{2}x_{38}^{2}x_{47}^{2}x_{49}^{2}x_{56}^{2}x_{57}^{2}x_{58}^{2}x_{59}^{2}x_{67}^{2}x_{68}^{2}x_{79}^{2}x_{89}^{2}}
22(5)\displaystyle\mathcal{I}^{(5)}_{22} =[Uncaptioned image]\displaystyle=\includegraphics[scale={0.4},align=c]{graph/DCI5l22.pdf} =[Uncaptioned image]\displaystyle=\includegraphics[scale={0.4},align=c]{graph/posi5l22.pdf} =a=59d4xax132x248x892x172x262x272x282x292x392x452x472x482x492x562x582x592x682x692x782\displaystyle=\int\prod_{a=5}^{9}d^{4}x_{a}\frac{x_{13}^{2}x_{24}^{8}x_{89}^{2}}{x_{17}^{2}x_{26}^{2}x_{27}^{2}x_{28}^{2}x_{29}^{2}x_{39}^{2}x_{45}^{2}x_{47}^{2}x_{48}^{2}x_{49}^{2}x_{56}^{2}x_{58}^{2}x_{59}^{2}x_{68}^{2}x_{69}^{2}x_{78}^{2}}
23(5)\displaystyle\mathcal{I}^{(5)}_{23} =[Uncaptioned image]\displaystyle=\includegraphics[scale={0.4},align=c]{graph/DCI5l23.pdf} =[Uncaptioned image]\displaystyle=\includegraphics[scale={0.4},align=c]{graph/posi5l23.pdf} =a=59d4xax132x248x392x172x262x272x282x292x352x362x452x472x482x492x562x592x692x782x892\displaystyle=\int\prod_{a=5}^{9}d^{4}x_{a}\frac{x_{13}^{2}x_{24}^{8}x_{39}^{2}}{x_{17}^{2}x_{26}^{2}x_{27}^{2}x_{28}^{2}x_{29}^{2}x_{35}^{2}x_{36}^{2}x_{45}^{2}x_{47}^{2}x_{48}^{2}x_{49}^{2}x_{56}^{2}x_{59}^{2}x_{69}^{2}x_{78}^{2}x_{89}^{2}}
24(5)\displaystyle\mathcal{I}^{(5)}_{24} =[Uncaptioned image]\displaystyle=\includegraphics[scale={0.4},align=c]{graph/DCI5l24.pdf} =[Uncaptioned image]\displaystyle=\includegraphics[scale={0.4},align=c]{graph/posi5l24.pdf} =a=59d4xax132x244x492x152x252x292x372x462x472x482x562x592x682x692x782x792x892\displaystyle=\int\prod_{a=5}^{9}d^{4}x_{a}\frac{x_{13}^{2}x_{24}^{4}x_{49}^{2}}{x_{15}^{2}x_{25}^{2}x_{29}^{2}x_{37}^{2}x_{46}^{2}x_{47}^{2}x_{48}^{2}x_{56}^{2}x_{59}^{2}x_{68}^{2}x_{69}^{2}x_{78}^{2}x_{79}^{2}x_{89}^{2}}
25(5)\displaystyle\mathcal{I}^{(5)}_{25} =[Uncaptioned image]\displaystyle=\includegraphics[scale={0.4},align=c]{graph/DCI5l25.pdf} =[Uncaptioned image]\displaystyle=\includegraphics[scale={0.4},align=c]{graph/posi5l25.pdf} =a=59d4xax132x244x892x192x262x292x382x452x482x572x582x592x672x682x692x782x792\displaystyle=\int\prod_{a=5}^{9}d^{4}x_{a}\frac{x_{13}^{2}x_{24}^{4}x_{89}^{2}}{x_{19}^{2}x_{26}^{2}x_{29}^{2}x_{38}^{2}x_{45}^{2}x_{48}^{2}x_{57}^{2}x_{58}^{2}x_{59}^{2}x_{67}^{2}x_{68}^{2}x_{69}^{2}x_{78}^{2}x_{79}^{2}}
26(5)\displaystyle\mathcal{I}^{(5)}_{26} =[Uncaptioned image]\displaystyle=\includegraphics[scale={0.4},align=c]{graph/DCI5l26.pdf} =[Uncaptioned image]\displaystyle=\includegraphics[scale={0.4},align=c]{graph/posi5l26.pdf} =a=59d4xax132x192x244x172x182x252x262x352x482x492x562x592x672x692x782x792x892\displaystyle=\int\prod_{a=5}^{9}d^{4}x_{a}\frac{x_{13}^{2}x_{19}^{2}x_{24}^{4}}{x_{17}^{2}x_{18}^{2}x_{25}^{2}x_{26}^{2}x_{35}^{2}x_{48}^{2}x_{49}^{2}x_{56}^{2}x_{59}^{2}x_{67}^{2}x_{69}^{2}x_{78}^{2}x_{79}^{2}x_{89}^{2}}
27(5)\displaystyle\mathcal{I}^{(5)}_{27} =[Uncaptioned image]\displaystyle=\includegraphics[scale={0.4},align=c]{graph/DCI5l27.pdf} =[Uncaptioned image]\displaystyle=\includegraphics[scale={0.4},align=c]{graph/posi5l27.pdf} =a=59d4xax134x244x492x152x162x252x292x382x392x462x472x482x562x592x672x782x792x892\displaystyle=\int\prod_{a=5}^{9}d^{4}x_{a}\frac{x_{13}^{4}x_{24}^{4}x_{49}^{2}}{x_{15}^{2}x_{16}^{2}x_{25}^{2}x_{29}^{2}x_{38}^{2}x_{39}^{2}x_{46}^{2}x_{47}^{2}x_{48}^{2}x_{56}^{2}x_{59}^{2}x_{67}^{2}x_{78}^{2}x_{79}^{2}x_{89}^{2}}
28(5)\displaystyle\mathcal{I}^{(5)}_{28} =[Uncaptioned image]\displaystyle=\includegraphics[scale={0.4},align=c]{graph/DCI5l28.pdf} =[Uncaptioned image]\displaystyle=\includegraphics[scale={0.4},align=c]{graph/posi5l28.pdf} =a=59d4xax132x192x244x152x182x252x262x392x482x492x562x572x672x692x782x792x892\displaystyle=\int\prod_{a=5}^{9}d^{4}x_{a}\frac{x_{13}^{2}x_{19}^{2}x_{24}^{4}}{x_{15}^{2}x_{18}^{2}x_{25}^{2}x_{26}^{2}x_{39}^{2}x_{48}^{2}x_{49}^{2}x_{56}^{2}x_{57}^{2}x_{67}^{2}x_{69}^{2}x_{78}^{2}x_{79}^{2}x_{89}^{2}}
29(5)\displaystyle\mathcal{I}^{(5)}_{29} =[Uncaptioned image]\displaystyle=\includegraphics[scale={0.4},align=c]{graph/DCI5l29.pdf} =[Uncaptioned image]\displaystyle=\includegraphics[scale={0.4},align=c]{graph/posi5l29.pdf} =a=59d4xax132x192x244x162x182x282x292x392x472x492x562x572x582x592x672x682x792\displaystyle=\int\prod_{a=5}^{9}d^{4}x_{a}\frac{x_{13}^{2}x_{19}^{2}x_{24}^{4}}{x_{16}^{2}x_{18}^{2}x_{28}^{2}x_{29}^{2}x_{39}^{2}x_{47}^{2}x_{49}^{2}x_{56}^{2}x_{57}^{2}x_{58}^{2}x_{59}^{2}x_{67}^{2}x_{68}^{2}x_{79}^{2}}
30(5)\displaystyle\mathcal{I}^{(5)}_{30} =[Uncaptioned image]\displaystyle=\includegraphics[scale={0.4},align=c]{graph/DCI5l30.pdf} =[Uncaptioned image]\displaystyle=\includegraphics[scale={0.4},align=c]{graph/posi5l30.pdf} =a=59d4xax132x246x182x262x282x292x392x452x472x492x562x572x582x672x682x792\displaystyle=\int\prod_{a=5}^{9}d^{4}x_{a}\frac{x_{13}^{2}x_{24}^{6}}{x_{18}^{2}x_{26}^{2}x_{28}^{2}x_{29}^{2}x_{39}^{2}x_{45}^{2}x_{47}^{2}x_{49}^{2}x_{56}^{2}x_{57}^{2}x_{58}^{2}x_{67}^{2}x_{68}^{2}x_{79}^{2}}
31(5)\displaystyle\mathcal{I}^{(5)}_{31} =[Uncaptioned image]\displaystyle=\includegraphics[scale={0.4},align=c]{graph/DCI5l31.pdf} =[Uncaptioned image]\displaystyle=\includegraphics[scale={0.4},align=c]{graph/posi5l31.pdf} =a=59d4xax132x196x244x152x162x172x182x282x292x392x452x492x562x592x672x692x782x792x892\displaystyle=\int\prod_{a=5}^{9}d^{4}x_{a}\frac{x_{13}^{2}x_{19}^{6}x_{24}^{4}}{x_{15}^{2}x_{16}^{2}x_{17}^{2}x_{18}^{2}x_{28}^{2}x_{29}^{2}x_{39}^{2}x_{45}^{2}x_{49}^{2}x_{56}^{2}x_{59}^{2}x_{67}^{2}x_{69}^{2}x_{78}^{2}x_{79}^{2}x_{89}^{2}}
32(5)\displaystyle\mathcal{I}^{(5)}_{32} =[Uncaptioned image]\displaystyle=\includegraphics[scale={0.4},align=c]{graph/DCI5l32.pdf} =[Uncaptioned image]\displaystyle=\includegraphics[scale={0.4},align=c]{graph/posi5l32.pdf} =a=59d4xax132x2410x192x252x262x272x282x292x382x452x462x472x482x492x562x582x672x792\displaystyle=\int\prod_{a=5}^{9}d^{4}x_{a}\frac{x_{13}^{2}x_{24}^{10}}{x_{19}^{2}x_{25}^{2}x_{26}^{2}x_{27}^{2}x_{28}^{2}x_{29}^{2}x_{38}^{2}x_{45}^{2}x_{46}^{2}x_{47}^{2}x_{48}^{2}x_{49}^{2}x_{56}^{2}x_{58}^{2}x_{67}^{2}x_{79}^{2}}
33(5)\displaystyle\mathcal{I}^{(5)}_{33} =[Uncaptioned image]\displaystyle=\includegraphics[scale={0.4},align=c]{graph/DCI5l33.pdf} =[Uncaptioned image]\displaystyle=\includegraphics[scale={0.4},align=c]{graph/posi5l33.pdf} =a=59d4xax132x244x162x282x292x372x462x472x562x572x582x592x692x782x892\displaystyle=\int\prod_{a=5}^{9}d^{4}x_{a}\frac{x_{13}^{2}x_{24}^{4}}{x_{16}^{2}x_{28}^{2}x_{29}^{2}x_{37}^{2}x_{46}^{2}x_{47}^{2}x_{56}^{2}x_{57}^{2}x_{58}^{2}x_{59}^{2}x_{69}^{2}x_{78}^{2}x_{89}^{2}}
34(5)\displaystyle\mathcal{I}^{(5)}_{34} =[Uncaptioned image]\displaystyle=\includegraphics[scale={0.4},align=c]{graph/DCI5l34.pdf} =[Uncaptioned image]\displaystyle=\includegraphics[scale={0.4},align=c]{graph/posi5l34.pdf} =a=59d4xax134x244x182x192x262x282x352x362x452x492x572x592x672x682x782x792\displaystyle=\int\prod_{a=5}^{9}d^{4}x_{a}\frac{x_{13}^{4}x_{24}^{4}}{x_{18}^{2}x_{19}^{2}x_{26}^{2}x_{28}^{2}x_{35}^{2}x_{36}^{2}x_{45}^{2}x_{49}^{2}x_{57}^{2}x_{59}^{2}x_{67}^{2}x_{68}^{2}x_{78}^{2}x_{79}^{2}}

At =5\ell=5, the original magic identity implies 9(5)=10(5)=11(5)=13(5)=15(5)=16(5)=17(5)=19(5)=20(5)=21(5)=22(5)=23(5)=31(5)=32(5)\mathcal{I}^{(5)}_{9}=\mathcal{I}^{(5)}_{10}=\mathcal{I}^{(5)}_{11}=\mathcal{I}^{(5)}_{13}=\mathcal{I}^{(5)}_{15}=\mathcal{I}^{(5)}_{16}=\mathcal{I}^{(5)}_{17}=\mathcal{I}^{(5)}_{19}=\mathcal{I}^{(5)}_{20}=\mathcal{I}^{(5)}_{21}=\mathcal{I}^{(5)}_{22}=\mathcal{I}^{(5)}_{23}=\mathcal{I}^{(5)}_{31}=\mathcal{I}^{(5)}_{32}.

References

  • (1) P. Heslop, The SAGEX Review on Scattering Amplitudes, Chapter 8: Half BPS correlators, J. Phys. A 55 (2022), no. 44 443009, [arXiv:2203.13019].
  • (2) L. F. Alday, B. Eden, G. P. Korchemsky, J. Maldacena, and E. Sokatchev, From correlation functions to Wilson loops, JHEP 09 (2011) 123, [arXiv:1007.3243].
  • (3) T. Adamo, M. Bullimore, L. Mason, and D. Skinner, A Proof of the Supersymmetric Correlation Function / Wilson Loop Correspondence, JHEP 08 (2011) 076, [arXiv:1103.4119].
  • (4) L. F. Alday and J. M. Maldacena, Gluon scattering amplitudes at strong coupling, JHEP 06 (2007) 064, [arXiv:0705.0303].
  • (5) J. M. Drummond, G. P. Korchemsky, and E. Sokatchev, Conformal properties of four-gluon planar amplitudes and Wilson loops, Nucl. Phys. B 795 (2008) 385–408, [arXiv:0707.0243].
  • (6) A. Brandhuber, P. Heslop, and G. Travaglini, MHV amplitudes in N=4 super Yang-Mills and Wilson loops, Nucl. Phys. B 794 (2008) 231–243, [arXiv:0707.1153].
  • (7) Z. Bern, L. J. Dixon, D. A. Kosower, R. Roiban, M. Spradlin, C. Vergu, and A. Volovich, The Two-Loop Six-Gluon MHV Amplitude in Maximally Supersymmetric Yang-Mills Theory, Phys. Rev. D 78 (2008) 045007, [arXiv:0803.1465].
  • (8) J. M. Drummond, J. Henn, G. P. Korchemsky, and E. Sokatchev, Hexagon Wilson loop = six-gluon MHV amplitude, Nucl. Phys. B 815 (2009) 142–173, [arXiv:0803.1466].
  • (9) B. Eden, G. P. Korchemsky, and E. Sokatchev, From correlation functions to scattering amplitudes, JHEP 12 (2011) 002, [arXiv:1007.3246].
  • (10) N. Arkani-Hamed, J. L. Bourjaily, F. Cachazo, S. Caron-Huot, and J. Trnka, The All-Loop Integrand For Scattering Amplitudes in Planar N=4 SYM, JHEP 01 (2011) 041, [arXiv:1008.2958].
  • (11) L. J. Mason and D. Skinner, The Complete Planar S-matrix of N=4 SYM as a Wilson Loop in Twistor Space, JHEP 12 (2010) 018, [arXiv:1009.2225].
  • (12) B. Eden, G. P. Korchemsky, and E. Sokatchev, More on the duality correlators/amplitudes, Phys. Lett. B 709 (2012) 247–253, [arXiv:1009.2488].
  • (13) S. Caron-Huot, Notes on the scattering amplitude / Wilson loop duality, JHEP 07 (2011) 058, [arXiv:1010.1167].
  • (14) B. Eden, P. Heslop, G. P. Korchemsky, and E. Sokatchev, The super-correlator/super-amplitude duality: Part I, Nucl. Phys. B 869 (2013) 329–377, [arXiv:1103.3714].
  • (15) B. Eden, P. Heslop, G. P. Korchemsky, and E. Sokatchev, The super-correlator/super-amplitude duality: Part II, Nucl. Phys. B 869 (2013) 378–416, [arXiv:1103.4353].
  • (16) R. G. Ambrosio, B. Eden, T. Goddard, P. Heslop, and C. Taylor, Local integrands for the five-point amplitude in planar N=4 SYM up to five loops, JHEP 01 (2015) 116, [arXiv:1312.1163].
  • (17) P. Heslop and V.-V. Tran, Multi-particle amplitudes from the four-point correlator in planar 𝒩\mathcal{N} = 4 SYM, JHEP 07 (2018) 068, [arXiv:1803.11491].
  • (18) F. Gonzalez-Rey, I. Y. Park, and K. Schalm, A Note on four point functions of conformal operators in N=4 superYang-Mills, Phys. Lett. B 448 (1999) 37–40, [hep-th/9811155].
  • (19) B. Eden, P. S. Howe, C. Schubert, E. Sokatchev, and P. C. West, Four point functions in N=4 supersymmetric Yang-Mills theory at two loops, Nucl. Phys. B 557 (1999) 355–379, [hep-th/9811172].
  • (20) B. Eden, P. S. Howe, C. Schubert, E. Sokatchev, and P. C. West, Simplifications of four point functions in N=4 supersymmetric Yang-Mills theory at two loops, Phys. Lett. B 466 (1999) 20–26, [hep-th/9906051].
  • (21) B. Eden, C. Schubert, and E. Sokatchev, Three loop four point correlator in N=4 SYM, Phys. Lett. B 482 (2000) 309–314, [hep-th/0003096].
  • (22) M. Bianchi, S. Kovacs, G. Rossi, and Y. S. Stanev, Anomalous dimensions in N=4 SYM theory at order g**4, Nucl. Phys. B 584 (2000) 216–232, [hep-th/0003203].
  • (23) B. Eden, P. Heslop, G. P. Korchemsky, and E. Sokatchev, Hidden symmetry of four-point correlation functions and amplitudes in N=4 SYM, Nucl. Phys. B 862 (2012) 193–231, [arXiv:1108.3557].
  • (24) B. Eden, P. Heslop, G. P. Korchemsky, and E. Sokatchev, Constructing the correlation function of four stress-tensor multiplets and the four-particle amplitude in N=4 SYM, Nucl. Phys. B 862 (2012) 450–503, [arXiv:1201.5329].
  • (25) J. Drummond, C. Duhr, B. Eden, P. Heslop, J. Pennington, and V. A. Smirnov, Leading singularities and off-shell conformal integrals, JHEP 08 (2013) 133, [arXiv:1303.6909].
  • (26) D. Chicherin, J. Drummond, P. Heslop, and E. Sokatchev, All three-loop four-point correlators of half-BPS operators in planar 𝒩\mathcal{N} = 4 SYM, JHEP 08 (2016) 053, [arXiv:1512.02926].
  • (27) D. Chicherin, A. Georgoudis, V. Gonçalves, and R. Pereira, All five-loop planar four-point functions of half-BPS operators in 𝒩=4\mathcal{N}=4 SYM, JHEP 11 (2018) 069, [arXiv:1809.00551].
  • (28) S. Caron-Huot and F. Coronado, Ten dimensional symmetry of 𝒩\mathcal{N} = 4 SYM correlators, JHEP 03 (2022) 151, [arXiv:2106.03892].
  • (29) T. Bargheer, T. Fleury, and V. Gonçalves, Higher-point integrands in 𝒩=4\mathcal{N}=4 super Yang-Mills theory, SciPost Phys. 15 (2023), no. 2 059, [arXiv:2212.03773].
  • (30) S. Caron-Huot, F. Coronado, and B. Mühlmann, Determinants in self-dual 𝒩\mathcal{N} = 4 SYM and twistor space, JHEP 08 (2023) 008, [arXiv:2304.12341].
  • (31) J. L. Bourjaily, A. DiRe, A. Shaikh, M. Spradlin, and A. Volovich, The Soft-Collinear Bootstrap: N=4 Yang-Mills Amplitudes at Six and Seven Loops, JHEP 03 (2012) 032, [arXiv:1112.6432].
  • (32) J. L. Bourjaily, P. Heslop, and V.-V. Tran, Perturbation Theory at Eight Loops: Novel Structures and the Breakdown of Manifest Conformality in N=4 Supersymmetric Yang-Mills Theory, Phys. Rev. Lett. 116 (2016), no. 19 191602, [arXiv:1512.07912].
  • (33) J. L. Bourjaily, P. Heslop, and V.-V. Tran, Amplitudes and Correlators to Ten Loops Using Simple, Graphical Bootstraps, JHEP 11 (2016) 125, [arXiv:1609.00007].
  • (34) S. He, C. Shi, Y. Tang, and Y.-Q. Zhang, The Cusp Limit of Correlators and A New Graphical Bootstrap for Correlators/Amplitudes to Eleven Loops, arXiv:2410.09859.
  • (35) J. L. Bourjaily, S. He, C. Shi, and Y. Tang, to appear, 2025. to appear.
  • (36) E. Remiddi and J. A. M. Vermaseren, Harmonic polylogarithms, Int. J. Mod. Phys. A 15 (2000) 725–754, [hep-ph/9905237].
  • (37) F. C. S. Brown, Polylogarithmes multiples uniformes en une variable, Compt. Rend. Math. 338 (2004), no. 7 527–532.
  • (38) L. J. Dixon, C. Duhr, and J. Pennington, Single-valued harmonic polylogarithms and the multi-Regge limit, JHEP 10 (2012) 074, [arXiv:1207.0186].
  • (39) O. Schnetz, Graphical functions and single-valued multiple polylogarithms, Commun. Num. Theor. Phys. 08 (2014) 589–675, [arXiv:1302.6445].
  • (40) A. B. Goncharov, M. Spradlin, C. Vergu, and A. Volovich, Classical Polylogarithms for Amplitudes and Wilson Loops, Phys. Rev. Lett. 105 (2010) 151605, [arXiv:1006.5703].
  • (41) C. Duhr, H. Gangl, and J. R. Rhodes, From polygons and symbols to polylogarithmic functions, JHEP 10 (2012) 075, [arXiv:1110.0458].
  • (42) C. Duhr, Hopf algebras, coproducts and symbols: an application to Higgs boson amplitudes, JHEP 08 (2012) 043, [arXiv:1203.0454].
  • (43) S. He, Y.-t. Huang, and C.-K. Kuo, to appear, 2025. to appear.
  • (44) N. I. Usyukina and A. I. Davydychev, Exact results for three and four point ladder diagrams with an arbitrary number of rungs, Phys. Lett. B 305 (1993) 136–143.
  • (45) B. Basso and L. J. Dixon, Gluing Ladder Feynman Diagrams into Fishnets, Phys. Rev. Lett. 119 (2017), no. 7 071601, [arXiv:1705.03545].
  • (46) B. Basso, L. J. Dixon, D. A. Kosower, A. Krajenbrink, and D.-l. Zhong, Fishnet four-point integrals: integrable representations and thermodynamic limits, JHEP 07 (2021) 168, [arXiv:2105.10514].
  • (47) F. Aprile and E. Olivucci, Multipoint fishnet Feynman diagrams: Sequential splitting, Phys. Rev. D 108 (2023), no. 12 L121902, [arXiv:2307.12984].
  • (48) F. Coronado, Perturbative four-point functions in planar 𝒩=4\mathcal{N}=4 SYM from hexagonalization, JHEP 01 (2019) 056, [arXiv:1811.00467].
  • (49) J. M. Drummond, J. Henn, V. A. Smirnov, and E. Sokatchev, Magic identities for conformal four-point integrals, JHEP 01 (2007) 064, [hep-th/0607160].
  • (50) D. J. Broadhurst and D. Kreimer, Knots and numbers in Phi**4 theory to 7 loops and beyond, Int. J. Mod. Phys. C 6 (1995) 519–524, [hep-ph/9504352].
  • (51) O. Schnetz, Quantum periods: A Census of phi**4-transcendentals, Commun. Num. Theor. Phys. 4 (2010) 1–48, [arXiv:0801.2856].
  • (52) F. C. S. Brown, On the periods of some Feynman integrals, arXiv:0910.0114.
  • (53) M. Borinsky and O. Schnetz, Recursive computation of Feynman periods, JHEP 22 (2022) 291, [arXiv:2206.10460].
  • (54) J. M. Drummond, Generalised ladders and single-valued polylogarithms, JHEP 02 (2013) 092, [arXiv:1207.3824].
  • (55) F. Brown and O. Schnetz, Single-valued multiple polylogarithms and a proof of the zig–zag conjecture, J. Number Theor. 148 (2015) 478–506.
  • (56) F. Chavez and C. Duhr, Three-mass triangle integrals and single-valued polylogarithms, JHEP 11 (2012) 114, [arXiv:1209.2722].
  • (57) O. Schnetz, Generalized single-valued hyperlogarithms, arXiv:2111.11246.
  • (58) K. A. Intriligator, Bonus symmetries of N=4 superYang-Mills correlation functions via AdS duality, Nucl. Phys. B 551 (1999) 575–600, [hep-th/9811047].
  • (59) B. Eden, A. C. Petkou, C. Schubert, and E. Sokatchev, Partial nonrenormalization of the stress tensor four point function in N=4 SYM and AdS / CFT, Nucl. Phys. B 607 (2001) 191–212, [hep-th/0009106].
  • (60) M. B. Green and C. Wen, Maximal U(1)Y-violating n-point correlators in 𝒩\mathcal{N} = 4 super-Yang-Mills theory, JHEP 02 (2021) 042, [arXiv:2009.01211].
  • (61) B. Basso, S. Komatsu, and P. Vieira, Structure Constants and Integrable Bootstrap in Planar N=4 SYM Theory, arXiv:1505.06745.
  • (62) A. V. Belitsky and G. P. Korchemsky, Octagon at finite coupling, JHEP 07 (2020) 219, [arXiv:2003.01121].
  • (63) A. V. Belitsky and G. P. Korchemsky, Crossing bridges with strong Szegő limit theorem, JHEP 04 (2021) 257, [arXiv:2006.01831].
  • (64) T. Fleury and S. Komatsu, Hexagonalization of Correlation Functions, JHEP 01 (2017) 130, [arXiv:1611.05577].
  • (65) B. Eden and A. Sfondrini, Tessellating cushions: four-point functions in 𝒩\mathcal{N} = 4 SYM, JHEP 10 (2017) 098, [arXiv:1611.05436].
  • (66) F. Brown, Single-valued Motivic Periods and Multiple Zeta Values, SIGMA 2 (2014) e25, [arXiv:1309.5309].
  • (67) E. Panzer and O. Schnetz, The Galois coaction on ϕ4\phi^{4} periods, Commun. Num. Theor. Phys. 11 (2017) 657–705, [arXiv:1603.04289].
  • (68) E. Panzer, Algorithms for the symbolic integration of hyperlogarithms with applications to Feynman integrals, Comput. Phys. Commun. 188 (2015) 148–166, [arXiv:1403.3385].
  • (69) M. Borinsky, Tropical Monte Carlo quadrature for Feynman integrals, Ann. Inst. H. Poincare D Comb. Phys. Interact. 10 (2023), no. 4 635–685, [arXiv:2008.12310].
  • (70) D. J. Binder, S. M. Chester, S. S. Pufu, and Y. Wang, 𝒩\mathcal{N} = 4 Super-Yang-Mills correlators at strong coupling from string theory and localization, JHEP 12 (2019) 119, [arXiv:1902.06263].
  • (71) S. M. Chester, Genus-2 holographic correlator on AdS5× S5 from localization, JHEP 04 (2020) 193, [arXiv:1908.05247].
  • (72) S. M. Chester, M. B. Green, S. S. Pufu, Y. Wang, and C. Wen, New modular invariants in 𝒩\mathcal{N} = 4 Super-Yang-Mills theory, JHEP 04 (2021) 212, [arXiv:2008.02713].
  • (73) D. Dorigoni, M. B. Green, and C. Wen, Exact properties of an integrated correlator in 𝒩\mathcal{N} = 4 SU(N) SYM, JHEP 05 (2021) 089, [arXiv:2102.09537].
  • (74) D. Dorigoni, M. B. Green, and C. Wen, Exact results for duality-covariant integrated correlators in 𝒩=4\mathcal{N}=4 SYM with general classical gauge groups, SciPost Phys. 13 (2022), no. 4 092, [arXiv:2202.05784].
  • (75) H. Paul, E. Perlmutter, and H. Raj, Integrated correlators in 𝒩\mathcal{N} = 4 SYM via SL(2, \mathbb{Z}) spectral theory, JHEP 01 (2023) 149, [arXiv:2209.06639].
  • (76) C. Wen and S.-Q. Zhang, Integrated correlators in 𝒩\mathcal{N} = 4 super Yang-Mills and periods, JHEP 05 (2022) 126, [arXiv:2203.01890].
  • (77) A. Brown, P. Heslop, C. Wen, and H. Xie, Integrated correlators in 𝒩=4\mathcal{N}=4 SYM beyond localisation, Phys. Rev. Lett. 132 (2024), no. 10 101602, [arXiv:2308.07219].
  • (78) D. J. Broadhurst and A. I. Davydychev, Exponential suppression with four legs and an infinity of loops, Nucl. Phys. B Proc. Suppl. 205-206 (2010) 326–330, [arXiv:1007.0237].
  • (79) B. Eden and V. A. Smirnov, Evaluating four-loop conformal Feynman integrals by D-dimensional differential equations, JHEP 10 (2016) 115, [arXiv:1607.06427].
  • (80) L. F. Alday and T. Hansen, Single-valuedness of the AdS Veneziano amplitude, JHEP 08 (2024) 108, [arXiv:2404.16084].
  • (81) S. Borowka, G. Heinrich, S. Jahn, S. P. Jones, M. Kerner, J. Schlenk, and T. Zirke, pySecDec: a toolbox for the numerical evaluation of multi-scale integrals, Comput. Phys. Commun. 222 (2018) 313–326, [arXiv:1703.09692].
  • (82) J. M. Drummond, J. M. Henn, and J. Trnka, New differential equations for on-shell loop integrals, JHEP 04 (2011) 083, [arXiv:1010.3679].
  • (83) N. Arkani-Hamed, W. Flieger, J. M. Henn, A. Schreiber, and J. Trnka, Coulomb Branch Amplitudes from a Deformed Amplituhedron Geometry, Phys. Rev. Lett. 132 (2024), no. 21 211601, [arXiv:2311.10814].
  • (84) W. Flieger, J. Henn, A. Schreiber, and J. Trnka, Two-loop four-point amplitudes on the Coulomb branch of 𝒩=4{\mathcal{N}}=4 super Yang-Mills, arXiv:2501.09454.