Planar master integrals for four-loop form factors

# Planar master integrals for four-loop form factors

Andreas von Manteuffel and Robert M. Schabinger
###### Abstract

We present the complete set of planar master integrals relevant to the calculation of three-point functions in four-loop massless Quantum Chromodynamics. Employing direct parametric integrations for a basis of finite integrals, we give analytic results for the Laurent expansion of conventional integrals in the parameter of dimensional regularization through to terms of weight eight.

\preprint

MSUHEP-19-008 institutetext: Department of Physics and Astronomy
Michigan State University, East Lansing, Michigan 48824, USA

## 1 Introduction

Due to their relevance to Drell-Yan lepton pair production Drell:1970wh () and Higgs boson production via gluon fusion Georgi:1977gs (); Wilczek:1977zn (); Shifman:1978zn (); Ellis:1979jy (); Inami:1982xt (), the basic quark and gluon form factors of massless Quantum Chromodynamics (QCD) have played a very important role in the development of the subject. For example, it has long been understood that the massless quark and gluon form factors provide a clean theoretical laboratory for the study of the dimensionally-regulated infrared singularities of perturbative scattering amplitudes in non-Abelian gauge theories Magnea:1990zb (). In particular, the cusp and collinear anomalous dimensions can be conveniently extracted from the and poles of the bare form factors Moch:2005tm ().

The four-loop cusp anomalous dimensions are especially relevant to cutting-edge analyses of Drell-Yan lepton production and gluon-fusion Higgs boson production because they are the last remaining ingredients required for a resummation of the next-to-next-to-next-to-leading Sudakov logarithms which are known only approximately Moch:2017uml (); Moch:2018wjh (). It is therefore unsurprising that the four-loop form factors and cusp anomalous dimensions of QCD have received significant recent attention vonManteuffel:2015gxa (); Henn:2016men (); Ruijl:2016pkm (); vonManteuffel:2016xki (); Lee:2016ixa (); Lee:2017mip (); Grozin:2017css (); Grozin:2018vdn (); Lee:2019zop (); Henn:2019rmi (); Bruser:2019auj (); vonManteuffel:2019wbj (). Similar studies in super Yang-Mills theory Boels:2012ew (); Boels:2015yna (); Boels:2017skl (); Boels:2017ftb () should ultimately provide a very useful cross-check on the QCD results via the principle of maximal transcendentality Kotikov:2002ab (). As a central building block for the complete calculation of the four-loop form factors of massless QCD, we extend previous results and present complete analytic expressions for all planar four-loop master integrals in this paper for the first time. The covering topologies are shown in Fig. 1.

This article is organized as follows. In Section 2, we discuss our conventions, notation, and setup. We summarize our computational method in Section 3 and our results for the most interesting integrals are provided in Section 4. Finally, we give an outlook in Section 5. We assemble the results for all 99 planar master integrals in the ancillary file ff4l-ints-pl.m on arXiv.org.

## 2 Preliminaries

In this section, we establish some notation and describe our enumeration of the planar four-loop form factor master integrals. We use Minkowskian propagators and choose an absolute normalization of

 Γ4(d/2−1)π2d (1)

for our four-loop Feynman integrals in spacetime dimensions. Our integrals depend on the virtuality, , in a trivial manner dictated by dimensional analysis. Therefore, without loss of generality, we set throughout this work. To understand our conventions, it suffices to study the four-loop generalized sunrise integral:

 =Γ4(1−ϵ)π8−4ϵ(4∏i=1∫d4−2ϵki)1k22k23k24(p1−k1+k2−k3+k4)2(p2+k1)2 ∣∣ ∣∣q2=−1 =−Γ(−3+4ϵ)Γ9(1−ϵ)Γ(5−5ϵ) =1ϵ(1576)+712304+ϵ(2681582944)+ϵ2(872675331776−548ζ3)+O(ϵ3). (2)

As a first non-trivial step, we construct a single Reduze 2 vonManteuffel:2012np (); Studerus:2009ye (); Bauer:2000cp () integral family (see Table 1) which covers all planar sectors (or topologies). To achieve this, we make highly-symmetric choices for the auxiliary propagators of the four-loop planar ladder form factor integral topology. At this stage, Reduze 2 allows for the construction of a compact sector selection encoding the minimal number of sectors for which integration by parts reductions are required. After carrying out integral reductions for all Feynman integrals with our in-house reduction code, we find just ninety-nine master integrals in ninety-seven sectors. In other words, only two of our master integral topologies are of the multi-component type. For these topologies, we prefer to work with squared propagators, marked with dotted edges at the level of graphs (see Figure 2). Finally, we remark that we use the physicists’ convention for ,

 ζ5,3=∞∑m=11m5m−1∑n=11n3≈0.0377076729848… (3)

## 3 Computational Details

In this section, we describe our computation of the master integrals. For all master integrals which diverge in four dimensions, we first derive dimensional recurrence relations Tarasov:1996br (); Lee:2012cn () and map to alternative finite integrals along the lines discussed in vonManteuffel:2014qoa (); vonManteuffel:2015gxa (). The Reduze 2 job find_finite_integrals efficiently generates a large number of finite integral candidates, which then facilitate the construction of the alternative finite basis. The change of basis to our alternative finite master integrals requires non-trivial integration by parts reductions for which we use finite field sampling and rational reconstruction vonManteuffel:2014ixa (); Hart2010 (), see also Peraro:2016wsq (). A key advantage of Finred, the private reduction program developed by one of us which realizes these ideas, is that it may be run in a highly-distributed manner on a computer cluster.

With these auxiliary integral reductions in hand, the problem reduces to one of finite Feynman integral evaluation. In all cases, the integrals are linearly reducible and accessible to HyperInt Panzer:2014caa (), a program for Feynman parametric integration, out of the box. The HyperInt program is capable of detecting and integrating out massless one-loop bubbles, a feature which pays off tremendously in a large number of cases. In fact, we find only twenty-eight integrals free of massless one-loop bubble insertions. Of these, five are two-point functions which were calculated already some time ago Baikov:2010hf (); Lee:2011jt (). This leaves us with just twenty-three non-trivial master integrals, for which we give explicit expressions through to weight eight in the next section.

Before proceeding, let us first stress one subtle point. It is essential for our workflow to avoid evaluating complicated finite four-loop Feynman integrals to excessively high orders in . In order to achieve this, it is necessary to select a basis of finite integrals for which complete weight eight information at the level of the finite basis integrals implies complete weight eight information at the level of the corresponding conventional basis integrals Schabinger:2018dyi (); vonManteuffel:2019wbj (). What is perhaps surprising is how restrictive this requirement turns out to be in some cases. For example, the Feynman integral

 \vbox\includegraphics[scale=0.175]{A_12_121295_dot}\put(−2.82,0.7)\makebox[4.3pt][l]\scalebox.8$(6−2ϵ)$=352ζ7−154ζ5ζ2+185ζ3ζ22−174ζ23−131140ζ32−754ζ5+12ζ3ζ2+110ζ22 +3ζ3+ϵ(−38320ζ5,3−292ζ5ζ3−92ζ23ζ2+835815250ζ42+267716ζ7+43ζ5ζ2+26120ζ3ζ22−35ζ23 −1943105ζ32−322ζ5−6ζ3ζ2+185ζ22+60ζ3)+O(ϵ2) (4)

appears to be the only one in this sector which is both finite in and allows for a faithful mapping of the weights in the sense just described.111For this sector, the Reduze 2 job find_finite_integrals finds nine finite integral candidates in dimensions which are not related by permutation symmetries.

## 4 Results

In this section, we present explicit expressions through to weight eight for the subset of planar four-loop three-point master integrals which do not have any massless one-loop bubble insertions. Results for the first two twelve-line topologies of Figure 1 were given, respectively, in references Henn:2016men () and Lee:2016ixa () and, very recently, results were given in reference Lee:2019zop () for the integrals of Eq. (4), Eq. (4), Eqs. (4)-(20), and Eqs. (4)-(4) below. To the best of our knowledge, the rest of the results which we present in this section are new. Curiously, it turns out that some of the eleven-line master integrals are every bit as challenging to calculate as the twelve-line master integrals in our approach. Eq. (4), for example, is actually more convenient to derive indirectly by evaluating a reducible twelve-line finite integral. In the following, we define our conventional master integrals in dimensions as described in Section 2 above and add a label to identify the topology in the conventions of Reduze 2.

 −12971152)+1ϵ3(45ζ22−143192ζ3+18364ζ2+129132304)+1ϵ2(8732ζ5+4316ζ3ζ2−74811440ζ22 +1393192ζ3−6371576ζ2−1214414608)+1ϵ(−27596ζ23+9955110080ζ32−8549192ζ5−7127288ζ3ζ2+160796ζ22 −23291576ζ3+348911152ζ2+3810473072)−20803192ζ7+593548ζ5ζ2+107340ζ3ζ22−83135576ζ23 −4475516720ζ32+38467192ζ5−153196ζ3ζ2−7639144ζ22+1900991152ζ3−461032304ζ2−37101316144 +ϵ(676960ζ5,3−237748ζ5ζ3+452916ζ23ζ2+406705331008000ζ42−160371128ζ7−1083596ζ5ζ2 −46127360ζ3ζ22+1725964ζ23+276817720160ζ32−150607192ζ5+1472932ζ3ζ2+5200432880ζ22−11001592304ζ3 −24991014608ζ2+125939414096)+O(ϵ2) (5) +253288ζ2+14596)+1ϵ3(−297160ζ22+14936ζ3−1181288ζ2−1669144)+1ϵ2(−616ζ5−13918ζ3ζ2 +362892880ζ22−7367288ζ3+674ζ2+11243144)+1ϵ(−4531144ζ23−31302110080ζ32+223116ζ5+3587144ζ3ζ2 −814811440ζ22+634748ζ3−435172ζ2−2275748)−4409764ζ7−5776ζ5ζ2−199316ζ3ζ22+910172ζ23 +2406371344ζ32−2496532ζ5−61144ζ3ζ2+99523480ζ22−88745144ζ3+1331972ζ2+382375144 +ϵ(493130ζ5,3−6038524ζ5ζ3+132136ζ23ζ2−7021967811008000ζ42+216878ζ7+998924ζ5ζ2 +3398211440ζ3ζ22−40519288ζ23−36819775040ζ32+15508748ζ5−14573ζ3ζ2−179053288ζ22+381811144ζ3 −1077724ζ2−2005247144)+O(ϵ2) (6) \vbox\includegraphics[scale=0.175]{A_11_30943}\put(−2.82,0.75)\makebox[4.3pt][l]\scalebox.6A_11_30943=1ϵ3(532ζ3)+1ϵ2(−19192ζ22)+1ϵ(47596ζ5−2512ζ3ζ2+3724ζ22−2516ζ3) −89596ζ23+243234032ζ32+1156ζ5−54ζ3ζ2−41396ζ22+8516ζ3+ϵ(12548ζ7+308548ζ5ζ2−123596ζ3ζ22 +41524ζ23+4079252ζ32−226516ζ5+1754ζ3ζ2−101596ζ22+516ζ3)+ϵ2(−17738ζ5,3−7581548ζ5ζ3 +405516ζ23ζ2+67210487403200ζ42+9104596ζ7−3103ζ5ζ2−2172ζ3ζ22+146548ζ23−67765672ζ32 (7) \vbox\includegraphics[scale=0.175]{A_11_27359}\put(−2.82,0.75)\makebox[4.3pt][l]\scalebox.6A_11_27359=1ϵ4(116ζ3)+1ϵ3(−1996ζ22)+1ϵ2(2916ζ5+23ζ3ζ2+2330ζ22+14ζ3) +1ϵ(26548ζ23+16213360ζ32−34ζ5−ζ3ζ2+340ζ22−234ζ3)−1038164ζ7+25312ζ5ζ2+252380ζ3ζ22 −20912ζ23+12599630ζ32+113ζ5+53ζ3ζ2−2639120ζ22+2134ζ3+ϵ(13920ζ5,3−1187524ζ5ζ3 +1714ζ23ζ2+426457144000ζ42+50456ζ7+171512ζ5ζ2−178320ζ3ζ22+194ζ23−28617280ζ32−2923ζ5 −313ζ3ζ2+24253120ζ22−14714ζ3)+O(ϵ2) (8) \vbox\includegraphics[scale=0.175]{A_11_14831}\put(−2.82,0.75)\makebox[4.3pt][l]\scalebox.6A_11_14831=1ϵ7(−1144)+1ϵ6(196)+1ϵ5(−19144ζ2+116)+1ϵ4(−25144ζ3+53144ζ2 −3748)+1ϵ3(−149360ζ22+125144ζ3−4972ζ2+26348)+1ϵ2(2548ζ5+3724ζ3ζ2+18536ζ22+4336ζ3 −6572ζ2−51916)+1ϵ(−7116ζ23+664995040ζ32+79148ζ5+41572ζ3ζ2−3539144ζ22−142336ζ3+140372ζ2 +838348)+38779192ζ7+635924ζ5ζ2−24941360ζ3ζ22−9035144ζ23+327831008ζ32−8638ζ5−30718ζ3ζ2 +60949720ζ22+1243336ζ3−1048172ζ2−4255748+ϵ(−378ζ5,3−78138ζ5ζ3+308972ζ23ζ2 +46907263100800ζ42−5155564ζ7+798724ζ5ζ2−67651360ζ3ζ22+43229ζ23−4906571680ζ32+1205924ζ5 −210518ζ3ζ2−140239720ζ22−8190736ζ3+6229172ζ2+6918116)+O(ϵ2) (9) +1ϵ2(2516ζ23−62695040ζ32)+1ϵ(−3881192ζ7−1074ζ5ζ2−169240ζ3ζ22)+25740ζ5,3−7298ζ5ζ3 +83916ζ23ζ2−17086001252000ζ42+O(ϵ) (10) \vbox\includegraphics[scale=0.175]{A_10_30155}\put(−2.82,0.75)\makebox[4.3pt][l]\scalebox.6A_10_30155=1ϵ5(−124ζ3)+1ϵ4(19240ζ22)+1ϵ3(−1724ζ5+1124ζ2ζ3)+1ϵ2(13ζ23 −3111260ζ32)+1ϵ(−5353384ζ7−23716ζ5ζ2+1069240ζ3ζ22)+255160ζ5,3+141712ζ5ζ3−1324ζ23ζ2 −332015984000ζ42+O(ϵ) (11) \vbox\includegraphics[scale=0.175]{A_10_30123}\put(−2.82,0.75)\makebox[4.3pt][l]\scalebox.6A_10_30123=1ϵ8(1144)+1ϵ6(18ζ2)+1ϵ5(1124ζ3)+1ϵ4(2316ζ22)+1ϵ3(414ζ5 +218ζ2ζ3)+1ϵ2(2716ζ23+13259840ζ32)+1ϵ(807164ζ7+1932ζ5ζ2−16ζ3ζ22)+32710ζ5,3 +11178ζ5ζ3−15778ζ23ζ2+554824142000ζ42+O(ϵ) (12) \vbox\includegraphics[scale=0.175]{A_10_30095}\put(−2.82,0.75)\makebox[4.3pt][l]\scalebox.6A_10_30095=1ϵ3(−516ζ3)+1ϵ2(−732ζ22+3516ζ3)+1ϵ(−12516ζ5+56ζ3ζ2+4932ζ22 −18516ζ3)−60548ζ23−7289672ζ32+87516ζ5−356ζ3ζ2−25932ζ22+87516ζ3+ϵ(−1356532ζ7 +216524ζ5ζ2−125348ζ3ζ22+423548ζ23+728996ζ32−462516ζ5+1856ζ3ζ2+122532ζ22−390516ζ3) +ϵ2(14794ζ5,3+2379524ζ5ζ3+492524ζ23ζ2−66262951201600ζ42+9495532ζ7−1515524ζ5ζ2 +877148ζ3ζ22−2238548ζ23−269693672ζ32+2187516ζ5−8756ζ3ζ2−546732ζ22+1683516ζ3) +O(ϵ3) (13) \vbox\includegraphics[scale=0.175]{A_10_27343}\put(−2.82,0.75)\makebox[4.3pt][l]\scalebox.6A_10_27343=1ϵ2(54ζ5)+1ϵ(78ζ23+1321ζ32−5ζ5)−377564ζ7−1558ζ5ζ2−740ζ3ζ22 −72ζ23−5221ζ32+20ζ5+ϵ(−192ζ5,3−23734ζ5ζ3−1834ζ23ζ2−33959600ζ42+377516ζ7+1552ζ5ζ2 +710ζ3ζ22+14ζ23+20821ζ32−80ζ5)+O(ϵ2) (14) \vbox\includegraphics[scale=0.175]{A_10_14767}\put(−2.82,0.75)\makebox[4.3pt][l]\scalebox.6A_10_14767=1ϵ7(−5576)+1ϵ6(5144)+1ϵ5(−596ζ2−536)+1ϵ4(−6796ζ3+524ζ2 +59)+1ϵ3(−269480ζ22+6724ζ3−56ζ2−209)+1ϵ2(22396ζ5−31348ζ3ζ2+269120ζ22−676ζ3 +103ζ2+809)+1ϵ(−53596ζ23−391815040ζ32−22324ζ5+31312ζ3ζ2−26930ζ22+1343ζ3−403ζ2 −3209)+307ζ7−986348ζ5ζ2−19348ζ3ζ22+53524ζ23+391811260ζ32+2236ζ5−3133ζ3ζ2+53815ζ22 −5363ζ3+1603ζ2+12809+ϵ(977960ζ5,3+1752116ζ5ζ3+531148ζ23ζ2+8076847252000ζ42−1228ζ7 +986312ζ5ζ2+19312ζ3ζ22−5356ζ23−39181315ζ32−4463ζ5+12523ζ3ζ2−215215ζ22+21443ζ3 −6403ζ2−51209)+O(ϵ2) (15) \vbox\includegraphics[scale=0.175]{A_9_30877}\put(−2.82,0.75)\makebox[4.3pt][l]\scalebox.6A_9_30877=1ϵ2(−32ζ23−433840ζ32)+1ϵ(−1574ζ7+514ζ5ζ2−15740ζ3ζ22)+46ζ5,3 +214ζ5ζ3+898ζ23ζ2−1621314200ζ42+O(ϵ) (16) \vbox\includegraphics[scale=0.175]{A_9_29899}\put(−2.82,0.75)\makebox[4.3pt][l]\scalebox.6A_9_29899=1ϵ(44164ζ7+52ζ5ζ2−740ζ3ζ22)+132ζ5,3+26ζ5ζ3+398ζ23ζ2+210395600ζ42 +O(ϵ) (17) \vbox\includegraphics[scale=0.175]{A_9_27339}\put(−2.82,0.75)\makebox[4.3pt][l]\scalebox.6A_9_27339=−4418ζ7−15ζ5ζ2−95ζ3ζ22+ϵ(365ζ5,3+51ζ5ζ3−749542625ζ42−441ζ7 −120ζ5ζ2−725ζ3ζ22)+O(ϵ2) (18) \vbox\includegraphics[scale=0.175]{A_9_14765}\put(−2.82,0.75)\makebox[4.3pt][l]\scalebox.6A_9_14765=1ϵ2(−20ζ5)+1ϵ(−2ζ23−27221ζ32−160ζ5)−16074ζ7+170ζ5ζ2 −1825ζ3ζ22−16ζ23−217621ζ32−1280ζ5+ϵ(4645ζ5,3+2432ζ5ζ3−120ζ23ζ2−1036125ζ42 −3214ζ7+1360ζ5ζ2−14565ζ3ζ22−128ζ23−1740821ζ32−10240ζ5)+O