Internal Reduction method for computing Feynman Integrals

A new approach to compute Feynman Integrals is presented. It relies on an integral representation of a given Feynman Integral in terms of simpler ones. Using this approach, we present, for the first time, results for a certain family of non-planar five-point two-loop Master Integrals with one external off-shell particle, relevant for instance for H + 2 jets production at the LHC, in both Euclidean and physical kinematical regions. ∗ Electronic address: costas.papadopoulos@cern.ch † Electronic address: christopher.wever@tum.de ar X iv :1 91 0. 06 27 5v 1 [ he pph ] 1 4 O ct 2 01 9


I. INTRODUCTION
Precise theoretical predictions are nowadays a vital element of discoveries in Physics, from gravitational wave astronomy [1] to high-energy physics [2,3]. The coming LHC Run 3 and the High Luminosity LHC Run scheduled after it, require the most precise theoretical predictions in order to fully exploit the machine's potential [4]. In the future, the FCC (Future Circular Collider) project will also further boost the demands in the direction of precision [5].
To be precise, next-to-next-to-leading order (NNLO) accuracy is the next level needed for physics analysis at the LHC, for the vast majority of QCD dominated scattering processes (see [6] and references therein). Mixed QCD-electroweak effects are also important for processes and observables involving gauge bosons productions [7][8][9], whereas when the NNLO corrections turn out to be large, even higher loop calculations (N 3 LO) are required [10].
Over the last years, NNLO QCD corrections for most of the 2 → 2 processes, including two-jet, top-pair and gauge bosons production, have been completed and already used in phenomenological and experimental studies [11]. This is the result of an intense theoretical work, which can be schematically classified in two frontiers: the frontier of NNLO radiative corrections and the frontier of two-loop amplitude computation. Building on the progress in these two frontiers, first results on NNLO QCD corrections for 2 → 3 processes start to emerge [12][13][14][15].
Two-loop amplitude computations require the reduction of the scattering matrix element in terms of basis integrals, usually referred to as Master Integrals (MI). Traditional reduction techniques based on integration-by-part identities [16][17][18] (IBP) are now more and more replaced by integrand reduction methods, following the one-loop paradigm [19]. Results for five-point two-loop amplitudes, relevant for three-jet/photon, V, H + 2 jets production have been recently presented [15]. A remarkable contradistinction with the NLO case, is that the basis of Master Integrals at two loops is still incomplete.
Multi-loop Master Integrals have been investigated for many years now. The most successful method to obtain analytic expressions and accurate numerical estimates of multi-scale multi-loop Feynman Integrals is, for the time being, the differential equations approach [20][21][22][23][24]. With the introduction of the canonical form of the differential equations [25], a major step towards the understanding of the mathematical structure of Feynman Integrals and sub-sequently of the scattering amplitudes has been achieved. The complexity of two-loop Feynman Integrals is determined by the number of internal massive propagators and the number of external particles, i.e. the total number of independent "kinematical" scales involved.
Feynman Integrals with a relatively small number of scales satisfy canonical differential equations and can be expressed in terms of multiple (or Goncharov) polylogarithms [26][27][28], a class of functions that have been well understood by now. Moreover, in the last couple of years, new mathematical structures [29][30][31][32][33] (elliptic polylogarithms) have been studied in order to obtain analytic insight of more complicated Feynman Integrals. With a complete basis of two-loop Master Integrals, it is hoped that an automation of NNLO calculations for arbitrary scattering processes can be achieved in the near future.
Five-point two-loop Master Integrals determine the current frontier of this endeavour.
The computation of all planar and non-planar five-point two-loop Master Integrals with massless internal propagators and light-like external momenta, has been recently completed [34,35]. The next step on this path of computing the five-point two-loop Master Integrals would be those with one of the external legs being off-shell. The planar and nonplanar topologies corresponding to these Master Integrals are shown in Fig. 1 and Fig. 2 respectively. Based on a simplified differential equations approach, we have also computed and expressed in terms of Goncharov polylogarithms, all Master Integrals for the first nontrivial planar family of five-point two-loop Master Integrals with massless internal propagators and one external particle carrying a space-or time-like momentum, P 1 in Fig. 1, as well as the full set of planar five-point two-loop massless Master Integrals with light-like external momenta [36]. In this paper we introduce a new idea, that we call "Internal Reduction", and present for the first time results for the non-planar family of five-point two-loop Master Integrals with massless internal propagators and one external particle carrying a space-or time-like momentum, N 1 in Fig. 2.
In section II, we briefly introduce the idea and present some key examples of its applicability. Section III gives a detailed presentation of the computation of the non-planar

II. INTERNAL REDUCTION
As two-loop Feynman Integrals with many scales (many external legs and/or internal masses) become more and more difficult to be expressed in analytic form, it is very welcome to investigate the possibility to compute them in terms of simpler ones, for which analytic results are available. In this respect, a Feynman parameter is introduced, in order to appropriately combine internal propagators of a multi-loop Feynman Integral. The simplest case is to combine two neighbouring propagators with the same loop momentum By appropriately choosing the propagators to be combined, the resulting Feynman Integral that appears under the integral of the introduced Feynman parameter x, i.e. the integrand, which contains one internal line less corresponds then either to a simpler topology with the same number of external legs or to a Feynman Integral with less external lines. The task then is straightforward: using known results for the integrand with one fewer internal line, the original Feynman Integral can be computed as an one-dimensional integration, done preferably analytically or in the worst case numerically. In the latter case, one expects that numerical integration shall be much more efficient, in comparison with more traditional approaches based on sector decomposition or Mellin-Barnes representation. We will call this process Internal Reduction, since the integrand contains the squared propagator on the r.h.s. in Eq. (1) and will typically itself be reducible to Master Integrals of the topology of which the integrand belongs to. Since this topology would contain one propagator less, its IBP reduction would be typically much easier than the reduction of integrals in the sector of the the original Feynman integral, that contains instead the two propagators on the l.h.s. of It is usually preferable to combine neighbouring propagators separated by an external leg carrying a light-like momentum, since in that case the new internal propagator is still massless. In several case, the reduced integral can be easily expressed in an analytic form using standard techniques, as for instance IBP identities, differential equations, etc. To illustrate the way Eq. (1) works in that case, a rather trivial example, shown graphically below k k + p 1 has been considered in detail in [37].
If the two neighbouring propagators are separated by an external leg carrying a time-like or a space-like momentum, the reduced graph will exhibit a non-zero internal mass. To illustrate how the method works in that case, we consider the following example, In Eq. (3), p 2 1 = m 1 , p 2 2 = 0, p 2 12 = m 3 and the internal massive propagator is given by Using IBP identities we can straightforwardly express the integrand appearing in the right-hand side of the Eq. (3) in terms of known Master Integrals [38], Now using the analytic results of [38] for the above Master Integrals, the integrand in Eq. (3) is, up to an overall factor of (M 2 ) −2 , with d = 4 − 2 , a function of the dimensionless Using the known properties of multiple polylogarithms G, the integrand can be expressed in terms of multiple polylogarithms with argument x and letters that depend on m 1 , m 3 and x. The integrand in the right-hand side of Eq. (3) is then easily integrable in the interval x ∈ [0, 1] and a straightforward integration results in the l.h.s. of Eq. (3), that, after replacingx → x above, agrees with the known results of ref. [39]. In fact, although results are shown above up to O( ) for the sake of brevity, we have checked the result up to order O( 2 ), namely up to transcendental weight four. Notice that, in this case, despite the fact that the Internal Reduction produces an integral with an internal mass, solutions of the sunrise graph, obtained from the differential equations approach, do allow for a straightforward analytic evaluation of the two-loop triangle, using Eq. (3).

III. APPLICATION TO PENTABOX INTEGRALS
In this section we will apply the method discussed in the previous section to a family of planar and nonplanar pentaboxes with one off-shell leg. In particular we will compute the The planar family of pentabox integrals of which Fig. 3 is the top sector is defined as while the nonplanar family of pentabox integrals in Fig. 4 is defined as The pentabox integrals in the two families above depend on six variables that we define as, We note that the MIs inside the planar family P have been computed previously in [36], but we decided to recompute the three in the top sector using the steps explained in this paper since it serves as a nontrivial check of the method. The top sector, shown in Fig. 3, contains the following three integrals, The nonplanar family N P contains 4 sectors, shown in Fig. 4 with ten previously unknown integrals that we choose as follows, All other MI in the family N P lie inside the doublebox families D, N D in Eqs. (14), (15) below and have been computed before [40][41][42][43].
We may now use the Internal Reduction method discussed in the previous section to relate the scalar planar pentabox MI to a planar double-box integral as shown below.
Similarly, the scalar nonplanar pentabox MI can be related to the non-planar double-box family as shown below.
The planar and nonplanar double box integral that appear inside the integral on the r.h.s.
in Eqs. (12) and (13) can be IBP reduced to a set of MIs that are known analytically and have been computed in [40,42,43] and [41][42][43] respectively. Their family of integrals is defined as The four external momenta p 1 , p 2 , p 3 , p 4 of the planar doublebox integrals may be expressed in terms of the five-point external momenta of the pentabox and similarly goes for the invariants, For the nonplanar doublebox, we flipped p 3 and p 4 in Eqs. (16), (17) (cf. Eq. (13)).
For the evaluation of double-box integrals we have used a common basis of MI taken from refs. [40,41] and two different parametrizations, namely the one given in [40,41] and the second one from [42]. In both cases the basis elements satisfy the following canonical equation with M (k) being constant matrices consisting of rational numbers. In the former case the so-called letters, α k , are given by α = {x, y, z, 1 + x, 1 − y, 1 − z, 1 + xy, z − y, 1 + y(1 + x) − z, xy + z, for the planar family, and α = {x, 1 + x, 1 − y, y, 1 + xy, 1 + x(1 + y − z), 1 − z, z − y, 1 + y − z 1 + y + xy − z, z, xy + z, 1 + x + xy − xz, 1 + xz}, for the non-planar family, with x, y, z defined by whereas in the second case the letters are given by α k ≡x − β k , with for the planar family, and for the non-planar family, withx, S 12 , S 23 , q defined by In both cases the solution of the equation (18) Eqs. (12), (13)) at x f = 1, where G i,(a) 111121100,res does not depend on x f but only on the five-point invariants in Eq (9) through Eq. (17). We note that the pentabox MI in Eqs. (10), (11) have been chosen such that their corresponding integrands have at most a (1−x f ) n−a singularity shown in Eq. (25) with n ≥ −1. The resummed expression for the integrand can be computed by making use of the differential equations of the corresponding doublebox MIs, Eq. (18). We may expand the differential equations around x f = 1, The The remaining integral in Eq. (27) is now integrable in the interval x f ∈ [0, 1] and may be safely expanded in . Its evaluation may be performed numerically with any integration routine that links to GiNaC to evaluate the GPLs of the four-point integrand G i 111121100 . For the second parametrization of the double-box integrals, we opt to replace the integration over x f with that over thex variable by using the following relation Eqs. (12) and (13).
In the ancillary files we provide the integrand for the ten nonplanar pentabox MI in Eq. (11) in terms of a canonical basis. For the canonical basis we have plugged in the two different solutions, described above. The first one can be used more efficiently in the physical region, where the five-point momenta q 1,2 are incoming, while the second parametrization in the Euclidean region. We have checked numerically that both solutions produce the same results in all regions. Furthermore, we also provide the resummed expression G i res around x f = 1 (x = 1) (as in Eq. (25)) for the canonical MI of the doublebox families i = D, N D, plus their expanded versions in , that is needed in Eq. (27).

IV. RESULTS AND CHECKS
In this section we discuss the integration over the Feynman parameter x f (x). The simplest way is to evaluate the integral in Eq. (27) numerically in for example Mathematica [45].
We provide a Mathematica notebook in the ancillary files where the reader can load the ingre-

dients of Eq. (27) and evaluate the integral over x f . A very naive integration in Mathematica
with NIntegrate takes about the order of ten minutes to compute all ten integrals on a home laptop to at least 5 digits of precision in a specific point.
We compare in Table I our results of the numerical integration for the ten new nonplanar pentaboxes with Secdec. The Secdec error is given for the last significant digit in brackets.
The error of our analytic result is smaller than the given digits and is left out in the table.
For the planar MI we compared our results in this paper with the results computed in [36] and we found excellent agreement, both in the Euclidean and physical region. We have chosen the following Euclidean and physical points, We managed to obtain numerical results in the Euclidean region with Secdec and we find very good agreement as seen from the table. We were however not able to obtain numerical results in the physical region with Secdec.
The numerical integration takes the longest time around x f = 1 (x = 1). Therefore one may consider speeding the integration up drastically by splitting the integration interval in The former evaluates in Mathematica in a manner of seconds, while the latter may be performed analytically by expanding the integrand around x f = 1 (x = 1) as explained in the previous section. In this way the numerical integration could in principle be performed in a few seconds on a home laptop.
One other way of speeding up the evaluation of the integrands is by performing the integrations analytically. For the first parametrization, this may be done by making a variable substitution that rationalizes the square-roots that appear in x f . For example, in the parametrization of ref. [40] the square-root that appears in the integration variable x f is and can be rationalized with a new variable v as follows The Goncharov polylogs that one finds have all rational letters in the new variable v. For the second parametrization [42] this is already achieved by the transformation given by Eq. (28).
By expressing the Goncharov polylogs in the form G(w 1 , . . . , w n ; v) (G(w 1 , . . . , w n ;x)) (e.g. using HyperInt [46]), with weights w i that do not depend on v (x), the integral over v (x) can be performed analytically and the result expressed again in terms of Goncharov polylogs.
On the other hand, we expect that analytic results for the canonical basis of the pentabox integral families with one off-shell leg, which is not known yet, will be much more manageable. It is therefore an interesting question how the method of Internal Reduction presented in this paper, can be used to relate directly elements of canonical bases of different families, i.e. those of pentabox and double-box integral families. We leave these considerations for a future publication.

V. CONCLUSIONS
In this paper we have presented a new method called Internal Reduction. By judiciously combining two propagators through a Feynman parameter representation one may reduce multi-scale (e.g. five-point) Feynman graphs to simpler Feynman graphs with less scales (e.g. four-point). Since there is already an expanding database of MI with up to five external legs that are known analytically, we may use this method to recycle these known MIs by relating new unknown multi-leg Feynman integrals to one-fold integrals of these known lower scale MIs.
As a proof of concept, we applied this method to a set of planar and nonplanar twoloop five-point MI with one off-shell leg. The planar family has been computed before and served as a non-trivial check of our implementation. The nonplanar family on the other hand contains ten new MI that are not known. In this paper we derived one-fold integral representations for these ten new MI by using the Internal Reduction method. We compared our results with Secdec in the Euclidean region and found excellent agreement. We have also presented results for these integrals in one of the six physical regions. The extension to the other five physical regions is left for future work. Together with this paper we include the results for the ten new nonplanar integrals in ancillary files. They are expressed in terms of one-fold integrals of Goncharov polylogs that can be evaluated efficiently with Ginac.
For the future, it would be interesting to apply this method to the other topologies of the five-point MI shown in Figs 1 and 2, providing thus the missing ingredients for NNLO QCD computations concerning for instance V, H + 2 jets production at the LHC. Besides this, one of the possible extensions of the Internal Reduction method is to reduce more internal lines introducing more Feynman parameters. In conclusion, the idea put forward in this paper, namely to relate through integration over Feynman (or possibly externally introduced) parameters more complicated integrals to simpler ones, for arbitrary multi-loop MI families, deserves further investigation.