Cuts of Feynman Integrals in Baikov representation

Hjalte Frellesvig; Costas G. Papadopoulos

arXiv:1701.07356·hep-ph·May 24, 2017

Cuts of Feynman Integrals in Baikov representation

Hjalte Frellesvig, Costas G. Papadopoulos

PDF

TL;DR

This paper introduces a systematic method using Baikov representation to compute cuts of Feynman integrals in $d$ dimensions, aiding in understanding the functions needed for their analytical expressions.

Contribution

It presents a novel systematic approach to compute Feynman integral cuts within the Baikov representation framework, extending the analytical tools available in quantum field theory.

Findings

01

Method for computing cuts in Baikov representation

02

Potential to identify function classes for full integrals

03

Framework applicable in $d$-dimensional space

Abstract

Based on the Baikov representation, we present a systematic approach to compute cuts of Feynman Integrals, appropriately defined in $d$ dimensions. The information provided by these computations may be used to determine the class of functions needed to analytically express the full integrals.

Figures1

Click any figure to enlarge with its caption.

Equations174

F_{α_{1} ... α_{N}} = \int (i = 1 \prod L \frac{d ^{d} k _{i}}{i π ^{d /2}}) \frac{1}{D _{1}^{α_{1}} ... D _{N}^{α_{N}}}

F_{α_{1} ... α_{N}} = \int (i = 1 \prod L \frac{d ^{d} k _{i}}{i π ^{d /2}}) \frac{1}{D _{1}^{α_{1}} ... D _{N}^{α_{N}}}

D_{a} = i = 1 \sum L j = i \sum M A_{_{a}}^{ij} s_{ij} + f_{a} = i = 1 \sum L j = i \sum L A_{_{a}}^{ij} k_{i} \cdot k_{j} + i = 1 \sum L j = L + 1 \sum M A_{_{a}}^{ij} k_{i} \cdot p_{j - L} + f_{a}, a = 1, \dots, N

D_{a} = i = 1 \sum L j = i \sum M A_{_{a}}^{ij} s_{ij} + f_{a} = i = 1 \sum L j = i \sum L A_{_{a}}^{ij} k_{i} \cdot k_{j} + i = 1 \sum L j = L + 1 \sum M A_{_{a}}^{ij} k_{i} \cdot p_{j - L} + f_{a}, a = 1, \dots, N

F_{α_{1} ... α_{N}} = C_{N}^{L} (G (p_{1}, ..., p_{E}))^{(- d + E + 1) /2} \int \frac{d x _{1} ... d x _{N}}{x _{1}^{α_{1}} ... x _{N}^{α_{N}}} P_{N}^{L} (x_{1} - f_{1}, ..., x_{N} - f_{N})^{(d - M - 1) /2}

F_{α_{1} ... α_{N}} = C_{N}^{L} (G (p_{1}, ..., p_{E}))^{(- d + E + 1) /2} \int \frac{d x _{1} ... d x _{N}}{x _{1}^{α_{1}} ... x _{N}^{α_{N}}} P_{N}^{L} (x_{1} - f_{1}, ..., x_{N} - f_{N})^{(d - M - 1) /2}

C_{N}^{L} = \frac{π ^{- L (L - 1) /4 - L E /2}}{\prod _{i = 1}^{L} Γ ( \frac{d - M + i}{2} )} det (A_{ij}^{a})

C_{N}^{L} = \frac{π ^{- L (L - 1) /4 - L E /2}}{\prod _{i = 1}^{L} Γ ( \frac{d - M + i}{2} )} det (A_{ij}^{a})

P_{N}^{L}\left({{x_{1}},{x_{2}},...,{x_{N}}}\right)={{{{{{G\left({{k_{1}},...,{k_{L}},{p_{1}},...,{p_{E}}}\right)}}}}}\Big{|}_{{s_{ij}}=\sum\limits_{a=1}^{N}{A_{ij}^{a}{x_{a}}\,\,\&\,\,\,{s_{ji}}={s_{ij}}}}}

P_{N}^{L}\left({{x_{1}},{x_{2}},...,{x_{N}}}\right)={{{{{{G\left({{k_{1}},...,{k_{L}},{p_{1}},...,{p_{E}}}\right)}}}}}\Big{|}_{{s_{ij}}=\sum\limits_{a=1}^{N}{A_{ij}^{a}{x_{a}}\,\,\&\,\,\,{s_{ji}}={s_{ij}}}}}

O_{ij} P_{N}^{L} = 0

O_{ij} P_{N}^{L} = 0

j \leq L (q_{j} = k_{j}) O_{ij} = d δ_{ij} + a = 1 \sum N b = 1 \sum N m = 1 \sum M A_{a}^{mi} A_{mj}^{b} (1 + δ_{mi}) (x_{b} - f_{b}) \frac{\partial}{\partial x _{a}}

j \leq L (q_{j} = k_{j}) O_{ij} = d δ_{ij} + a = 1 \sum N b = 1 \sum N m = 1 \sum M A_{a}^{mi} A_{mj}^{b} (1 + δ_{mi}) (x_{b} - f_{b}) \frac{\partial}{\partial x _{a}}

j > L (q_{j} = p_{j - L}) O_{ij} = a = 1 \sum N (m = 1 \sum L b = 1 \sum N A_{a}^{mi} A_{mj}^{b} (1 + δ_{mi}) (x_{b} - f_{b}) + m = L + 1 \sum M A_{a}^{mi} s_{mj}) \frac{\partial}{\partial x _{a}}

j > L (q_{j} = p_{j - L}) O_{ij} = a = 1 \sum N (m = 1 \sum L b = 1 \sum N A_{a}^{mi} A_{mj}^{b} (1 + δ_{mi}) (x_{b} - f_{b}) + m = L + 1 \sum M A_{a}^{mi} s_{mj}) \frac{\partial}{\partial x _{a}}

x_{1} = k^{2} - m_{1}^{2}, x_{2} = (k + p_{1})^{2} - m_{2}^{2}, \dots, x_{N} = (k + p_{1} + \dots + p_{N - 1})^{2} - m_{N}^{2}

x_{1} = k^{2} - m_{1}^{2}, x_{2} = (k + p_{1})^{2} - m_{2}^{2}, \dots, x_{N} = (k + p_{1} + \dots + p_{N - 1})^{2} - m_{N}^{2}

F_{α_{1} \dots α_{N}} = C_{N}^{1} G (p_{1}, \dots, p_{N - 1})^{(N - d) /2} \int \frac{d x _{1} ... d x _{N}}{x _{1}^{α_{1}} ... x _{N}^{α_{N}}} P_{N}^{1}^{(d - N - 1) /2}

F_{α_{1} \dots α_{N}} = C_{N}^{1} G (p_{1}, \dots, p_{N - 1})^{(N - d) /2} \int \frac{d x _{1} ... d x _{N}}{x _{1}^{α_{1}} ... x _{N}^{α_{N}}} P_{N}^{1}^{(d - N - 1) /2}

C_{N}^{1} = \frac{π ^{- (N - 1) /2}}{Γ ( \frac{d - N + 1}{2} )} (\frac{1}{2})^{N - 1}

C_{N}^{1} = \frac{π ^{- (N - 1) /2}}{Γ ( \frac{d - N + 1}{2} )} (\frac{1}{2})^{N - 1}

F_{α_{1} ... α_{N - 1} 0}

F_{α_{1} ... α_{N - 1} 0}

x_{N}^{-} \int x_{N}^{+} d x_{N} P_{N}^{1}^{(d - N - 1) /2}

x_{N}^{-} \int x_{N}^{+} d x_{N} P_{N}^{1}^{(d - N - 1) /2}

x_{1}

x_{1}

x_{4}

x_{7}

F_{0111001} = C_{3}^{2} (G (p_{1}, p_{2}))^{- (d - 3) /2} \int \frac{d x _{1} \dots d x _{7}}{x _{2} x _{3} x _{4} x _{7}} (P_{3}^{2})^{(d - 5) /2}

F_{0111001} = C_{3}^{2} (G (p_{1}, p_{2}))^{- (d - 3) /2} \int \frac{d x _{1} \dots d x _{7}}{x _{2} x _{3} x _{4} x _{7}} (P_{3}^{2})^{(d - 5) /2}

F_{0111001} = \int \frac{d ^{d} k _{1}}{i π ^{d /2}} \frac{d ^{d} k _{2}}{i π ^{d /2}} \frac{1}{x _{2} x _{3} x _{4} x _{7}} = C_{2}^{1} \int \frac{d ^{d} k _{1}}{i π ^{d /2}} \frac{1}{x _{2} x _{3}} (x_{1} + m_{1}^{2})^{- (d - 2) /2} \int \frac{d x _{4} d x _{7}}{x _{4} x _{7}} (P_{2}^{1})^{(d - 3) /2}

F_{0111001} = \int \frac{d ^{d} k _{1}}{i π ^{d /2}} \frac{d ^{d} k _{2}}{i π ^{d /2}} \frac{1}{x _{2} x _{3} x _{4} x _{7}} = C_{2}^{1} \int \frac{d ^{d} k _{1}}{i π ^{d /2}} \frac{1}{x _{2} x _{3}} (x_{1} + m_{1}^{2})^{- (d - 2) /2} \int \frac{d x _{4} d x _{7}}{x _{4} x _{7}} (P_{2}^{1})^{(d - 3) /2}

F_{0111001} = C_{2}^{1} C_{3}^{1} (G (p_{1}, p_{2}))^{- (d - 3) /2} \int \frac{d x _{1} d x _{2} d x _{3} d x _{4} d x _{7}}{x _{2} x _{3} x _{4} x _{7}} (x_{1} + m_{1}^{2})^{- (d - 2) /2} (P_{3}^{1})^{(d - 5) /2} (P_{2}^{1})^{(d - 3) /2}

F_{0111001} = C_{2}^{1} C_{3}^{1} (G (p_{1}, p_{2}))^{- (d - 3) /2} \int \frac{d x _{1} d x _{2} d x _{3} d x _{4} d x _{7}}{x _{2} x _{3} x _{4} x _{7}} (x_{1} + m_{1}^{2})^{- (d - 2) /2} (P_{3}^{1})^{(d - 5) /2} (P_{2}^{1})^{(d - 3) /2}

\frac{\partial}{\partial X} F_{α_{1} ... α_{N}}

\frac{\partial}{\partial X} F_{α_{1} ... α_{N}}

b \frac{\partial P _{N}^{L}}{\partial X} + a \sum c_{a} \frac{\partial P _{N}^{L}}{\partial x _{a}} = 0

b \frac{\partial P _{N}^{L}}{\partial X} + a \sum c_{a} \frac{\partial P _{N}^{L}}{\partial x _{a}} = 0

\frac{\partial}{\partial X} F_{α_{1} ... α_{N}}

\frac{\partial}{\partial X} F_{α_{1} ... α_{N}}

\frac{\partial}{\partial X} F_{α_{1} ... α_{N}}

\frac{\partial}{\partial X} F_{α_{1} ... α_{N}}

\frac{1}{x _{1}^{α^{'}_{1}} ... x _{N}^{α^{'}_{N}}}

\frac{1}{x _{1}^{α^{'}_{1}} ... x _{N}^{α^{'}_{N}}}

\frac{\partial}{\partial X} F_{α_{1} ... α_{N}} = i \sum R_{i} F_{α_{1}^{(i)} ... α_{N}^{(i)}}

\frac{\partial}{\partial X} F_{α_{1} ... α_{N}} = i \sum R_{i} F_{α_{1}^{(i)} ... α_{N}^{(i)}}

F_{α_{1} ... α_{N}} ∣_{n \times cut} \equiv C_{N}^{L} (G)^{(- d + E + 1) /2} (a = n + 1 \prod N \int d x_{a}) c = 1 \prod n x_{c = 0} \oint d x_{c} \frac{1}{x _{1}^{α_{1}} \dots x _{N}^{α_{N}}} P_{N}^{L}^{(d - M - 1) /2}

F_{α_{1} ... α_{N}} ∣_{n \times cut} \equiv C_{N}^{L} (G)^{(- d + E + 1) /2} (a = n + 1 \prod N \int d x_{a}) c = 1 \prod n x_{c = 0} \oint d x_{c} \frac{1}{x _{1}^{α_{1}} \dots x _{N}^{α_{N}}} P_{N}^{L}^{(d - M - 1) /2}

\frac{\partial}{\partial X _{j}} F_{i} = l = 1 \sum I M_{i l}^{(j)} F_{l}

\frac{\partial}{\partial X _{j}} F_{i} = l = 1 \sum I M_{i l}^{(j)} F_{l}

\frac{\partial}{\partial X _{j}} F_{i} ∣_{n \times cut} = l = 1 \sum I M_{i l}^{(j)} F_{l} ∣_{n \times cut}

\frac{\partial}{\partial X _{j}} F_{i} ∣_{n \times cut} = l = 1 \sum I M_{i l}^{(j)} F_{l} ∣_{n \times cut}

F_{α_{1} \dots α_{N}} = \int \frac{d ^{d} k}{i π ^{d /2}} \frac{1}{D _{1}^{α_{1}} \dots D _{N}^{α_{N}}}

F_{α_{1} \dots α_{N}} = \int \frac{d ^{d} k}{i π ^{d /2}} \frac{1}{D _{1}^{α_{1}} \dots D _{N}^{α_{N}}}

F_{α_{1} \dots α_{N}} = \frac{π ^{- E /2}}{i Γ (( d - E ) /2 )} \int \frac{λ ^{(d - E - 2) /2}}{D _{1}^{α_{1}} \dots D _{N}^{α_{N}}} d^{E} k_{∣∣} d λ

F_{α_{1} \dots α_{N}} = \frac{π ^{- E /2}}{i Γ (( d - E ) /2 )} \int \frac{λ ^{(d - E - 2) /2}}{D _{1}^{α_{1}} \dots D _{N}^{α_{N}}} d^{E} k_{∣∣} d λ

F_{α_{1} \dots α_{N}} = \frac{i π ^{- E /2}}{Γ (( d - E ) /2 )} \frac{det ( A ^{- 1} )}{- G} \int \frac{λ ^{(d - E - 2) /2}}{x _{1}^{α_{1}} \dots x _{N}^{α_{N}}} d^{E + 1} x

F_{α_{1} \dots α_{N}} = \frac{i π ^{- E /2}}{Γ (( d - E ) /2 )} \frac{det ( A ^{- 1} )}{- G} \int \frac{λ ^{(d - E - 2) /2}}{x _{1}^{α_{1}} \dots x _{N}^{α_{N}}} d^{E + 1} x

Peer Reviews

No public reviews on file for this paper yet. If you reviewed it on a platform where reviews are public (OpenReview, ICLR, NeurIPS, ICML), you can paste yours below so the community can read it here.

Videos

No videos yet. Explain this paper in a talk, walkthrough, or lecture? Add one.

Full text

††institutetext: Institute of Nuclear and Particle Physics, NCSR ‘Demokritos’, Agia Paraskevi, 15310, Greece

Cuts of Feynman Integrals in Baikov representation

Hjalte Frellesvig and Costas G. Papadopoulos

[email protected]

Abstract

Based on the Baikov representation, we present a systematic approach to compute cuts of Feynman Integrals, appropriately defined in $d$ dimensions. The information provided by these computations may be used to determine the class of functions needed to analytically express the full integrals.

Keywords:

Feynman integrals, QCD, NLO and NNLO calculations

1 Introduction

It is almost seventy years from the time Feynman Integrals (FI) were first introduced Feynman:1949zx ; Dyson:1949bp ; Dyson:1949ha and forty-five years since the dimensional regularisation tHooft:1972tcz set up the framework for an efficient use of loop integrals in computing scattering matrix elements, and still the frontier of multi-scale multi-loop integral calculations (maximal both in number of scales and number of loops) is determined by the planar five-point two-loop on-shell massless integrals Gehrmann:2015bfy ; Papadopoulos:2015jft , recently computed111Complete results, including physical region kinematics, are presented in Papadopoulos:2015jft . Notice that numerical codes, like for instance SecDec Borowka:2015mxa , can reproduce analytic results only at Euclidean region kinematics; results for physical region kinematics are not supported due to poor numerical convergence.. On the other hand, in order to keep up with the increasing experimental accuracy as more data is collected at the LHC, more precise theoretical predictions and higher loop calculations are required Andersen:2014efa .

In the last years our understanding of the reduction of one-loop amplitudes to a set of Master Integrals (MI), a minimal set of FI that form a basis, either based on unitarity methods Bern:1994cg ; Bern:1994zx ; Berger:2008sj or at the integrand level via the OPP method Ossola:2006us ; Ossola:2008xq , has drastically changed the way one-loop calculations are preformed resulting in many fully automated numerical tools (some reviews on the topic are AlcarazMaestre:2012vp ; Ellis:2011cr ; vanDeurzen:2015jmn ), making the next-to-leading order (NLO) approximation the default precision for theoretical predictions at the LHC. In the recent years, progress has been made also towards the extension of these reduction methods for two-loop amplitudes at the integral Gluza:2010ws ; Kosower:2011ty ; CaronHuot:2012ab ; Johansson:2012zv ; Johansson:2012sf ; Johansson:2013sda ; Sogaard:2013fpa ; Ita:2015tya ; Johansson:2015ava ; Mastrolia:2016dhn as well as the integrand Mastrolia:2011pr ; Badger:2012dp ; Mastrolia:2012wf ; Badger:2013gxa ; Papadopoulos:2013hra ; Badger:2015lda level. Two-loop MI are defined using the integration by parts (IBP) identities Chetyrkin:1981qh ; Tkachov:1981wb ; Laporta:2001dd , an indispensable tool beyond one loop. Contrary to the one-loop case, where MI have been known for a long time already 'tHooft:1978xw , a complete library of MI at two-loops is still missing. At the moment this is the main obstacle to obtain a fully automated NNLO calculation framework similar to the one-loop one, that will satisfy the precision requirements at the LHC Andersen:2014efa .

Many methods have been introduced in order to compute MI Smirnov:2012gma . The overall most successful one, is based on expressing the FI in terms of an integral representation over Feynman parameters, involving the two well-known Symanzik Polynomials $U$ and $F$ Bogner:2010kv . The introduction of the sector decomposition Hepp:1966eg ; Roth:1996pd ; Binoth:2000ps ; Binoth:2003ak ; Bogner:2007cr method resulted in a powerful computational framework for the numerical evaluation of FI, see for instance SecDec Borowka:2015mxa . An alternative is based on Mellin-Barnes representation Smirnov:1999gc ; Tausk:1999vh , implemented in Czakon:2005rk 222See also https://mbtools.hepforge.org. Nevertheless, the most successful method to calculate multi-scale multi-loop FI is, for the time being, the differential equations (DE) approach Kotikov:1990kg ; Kotikov:1991pm ; Bern:1992em ; Remiddi:1997ny ; Gehrmann:1999as , which has been used in the past two decades to calculate various MI at two-loops and beyond. Following the work of refs. Goncharov:1998kja ; Remiddi:1999ew ; Goncharov:2001iea , there has been a building consensus that the so-called Goncharov Polylogarithms (GPs) form a functional basis for many MI. The so-called canonical form of DE, introduced by Henn Henn:2013pwa , manifestly results in MI expressed in terms of GPs 333For an alternative method in the single scale case see also ref. Ablinger:2015tua . Nevertheless the reduction of a given DE to a canonical form is by no means fully understood. First of all, despite recent efforts Lee:2014ioa ; Meyer:2016slj ; Gituliar:2017vzm , and the existence of sufficient conditions that a given MI can be expressed in terms of GPs, no criterion, with practical applicability, that is at the same time necessary and sufficient has been introduced so far. Moreover, it is well known that when for instance enough internal masses are introduced, MI are not anymore expressible in terms of GPs, and in fact a new class of functions involving elliptic integrals is needed Adams:2015gva ; Bonciani:2016qxi .

In this paper we are studying another representation of FI introduced by Baikov Baikov:1996iu ; Baikov:1996rk ; Smirnov:2003kc ; Lee:2010wea ; Grozin:2011mt ; Lee:2013hzt . As we will see in Section 2, where we present a review of its derivation, it has several nice features, including its conceptual simplicity, a direct factorisation of kinematics and loop ‘topology’, incorporation of IBP identities in a straightforward manner. We also present, to the best of our knowledge for the first time, a consistent definition of the integration limits which will be important for the computation of cut integrals in $d$ dimensions. In Section 3 we elaborate on the loop-by-loop approach within the Baikov representation, that has a minimal number of integration variables for a given FI. In Section 4 we present a novel approach to obtain DE from the Baikov representation. In Section 5 we introduce the definition of the cut integral in Baikov representation, that satisfies the same DE and the same IBP identities as the uncut one Henn:2014qga , and we conjecture that computing the corresponding maximally cut integral Primo:2016ebd we may have a necessary and sufficient condition for the expression of the uncut integral in terms of GPs and when applied to the whole family of MI on the possibility to obtain a canonical form. Finally in Appendix A we present an alternative derivation of the Baikov representation for one- and two-loop FI and in Appendix B we collect several examples of maximally cut integrals in Baikov representation that support our findings.

2 The Baikov representation

In this section we introduce the Baikov representation following ref. Lee:2010wea ; Grozin:2011mt . An $L$ -loop Feynman Integral with $E+1$ external lines can be written in the form

[TABLE]

with $N=\frac{{L\left({L+1}\right)}}{2}+LE$ , $\alpha_{i}$ arbitrary integers, and $D_{a}$ , $a=1,...,N$ , inverse Feynman propagators, $P^{2}-M^{2}$ , where $P$ represents, collectively, a linear combination of loop and external momenta and $M$ internal masses, as dictated by the Feynman Integral in consideration.

To be more specific, let us define $q_{i}=k_{i},(i=1,...,L)$ the loop momenta and $q_{L+i}=p_{i},(i=1,...,E)$ , the independent external momenta, $M=L+E$ and $s_{ij}=q_{i}\cdot q_{j}$ . Then

[TABLE]

where $f_{a}$ depend on external kinematics and internal masses. $A_{{}_{a}}^{ij}$ can be understood as an $N\times N$ matrix, with $a$ running obviously from $1$ to $N$ and with $(ij)$ taking also $N$ values as $i=1,\ldots,L$ and $j=i,\ldots,M$ . The elements of the matrix $A_{{}_{a}}^{ij}$ are integer numbers taken from the set $\left\{-2,-1,0,+1,+2\right\}$ . This matrix is characteristic of the corresponding Feynman graph and can, in a loose sense, be associated with the ‘topology’ of the graph. Then, by projecting each of the loop momenta $q_{i}=k_{i},(i=1,...,L)$ with respect to the space spanned by the external momenta involved plus a transverse component (for details see Grozin:2011mt ), we may write

[TABLE]

with

[TABLE]

and

[TABLE]

with $G$ representing the Gram determinant, $G\left({{q_{1}},\ldots,{q_{n}}}\right)=\det\left({{q_{i}}\cdot{q_{j}}}\right)$ and $A_{ij}^{a}$ is the inverse of the topology matrix $A_{{}_{a}}^{ij}$ . An alternative derivation of the Baikov representation for one- and two-loop FI is given in Appendix A.

Integration-by-parts identities can easily be accommodated in the Baikov representation. The generators of the IBP identities can be cast into the form

[TABLE]

with the operators $O_{ij}$ given by ( $i=1,\ldots,L$ )

[TABLE]

and

[TABLE]

The derivation of the Baikov representation can easily be implemented in a computer algebra code444A Mathematica script, Baikov.m, is provided as an attachment.

We conclude this section by elaborating on the limits of the $x_{a}-$ integrations in Eq. (3). In order to simplify the discussion, let us start with a generic one-loop configuration defined by

[TABLE]

Then consider the generic integral ( $\alpha_{i}\geq 0$ ),

[TABLE]

It is easy to verify that ${P^{1}_{N}}$ is a polynomial that is quadratic in the variables $x_{a}$ Smirnov:2003kc , and that obviously when $\alpha_{N}=0$ , the external momentum $p_{N-1}$ decouples, so that

[TABLE]

where ${P^{1}_{N}}\left({x_{N}^{+}}\right)={P^{1}_{N}}\left({x_{N}^{-}}\right)=0$ and

[TABLE]

using $P_{N}^{1}=\frac{1}{4}G\left({{p_{1}},\ldots,{p_{N-2}}}\right)\left({x_{N}^{+}-{x_{N}}}\right)\left({{x_{N}}-x_{N}^{-}}\right)$ and ${\left({x_{N}^{+}-x_{N}^{-}}\right)^{2}}=16\frac{{G\left({{p_{1}},\ldots,{p_{N-1}}}\right)}}{{G{{\left({{p_{1}},\ldots,{p_{N-2}}}\right)}^{2}}}}P_{N-1}^{1}$ . This can be repeated straightforwardly for all variables except $x_{1}=k^{2}-m_{1}^{2}$ whose integration limits are simply derived from the $k-$ modulus integration limits. The generalisation to the two-loop case is straightforward, with the integration at each step performed over the $x-$ variables involving a given external momentum, and the last ones derived by the corresponding $k_{1}-$ and $k_{2}-$ modulus integration limits. We have checked both analytically and numerically that the limits, as defined above, reproduce the known results for several examples at one and two loops.

3 The loop-by-loop approach

At one loop the number of Baikov variables and the number of propagators of a generic Feynman Integral can be the same. At higher loops this is not the case anymore. Consider for instance the double-box, with $E=3$ independent external momenta and $L=2$ , giving $N=9$ and $M=5$ . In that case we have to integrate over nine variables whereas the scalar double-box graph has only seven denominators. Of course this mismatch is the origin of the well known irreducible scalar products (ISP) appearing for $L\geq 2$ . It is eventually desirable to have a representation with the minimal number of integration variables. This can be achieved by considering the projection over the space spanned by the external momenta of each loop integration momentum separately. As an example to illustrate the concept, consider the two-loop three-point Feynman Integral of Figure 1. The standard Baikov representation is based on the following set of seven inverse propagators

[TABLE]

For the integral ${F_{0111001}}$ , in the loop-by-loop approach one can consider the $k_{1}$ and $k_{2}$ integrations separately. Starting with $k_{2}$ integration it is easy to see that we have a two-point function with external momentum $k_{1}$ .

[TABLE]

Then for the $k_{1}$ integration we have a three-point integral with two external momenta $p_{1}$ and $p_{2}$ .

[TABLE]

Notice that the same result can be obtained from Eq. (12) integrating out $x_{5}$ and $x_{6}$ , as detailed at the end of the previous section. Nevertheless, the loop-by-loop approach offers an alternative way of obtaining a minimal number of integration variables.

4 Deriving differential equations

It is instructive to study the way differential equations with respect to external kinematics and masses can be obtained in Baikov representation. Differential equations are usually written in terms of external kinematical invariants, $s_{ij}=(p_{i}+p_{j})^{2}$ and internal masses, $m_{i}^{2}$ . In the standard approach, since the integral in the momentum-space representation is not an explicit function of the kinematical invariants, derivatives with respect to external momenta, $p_{j}^{\mu}\frac{\partial}{{\partial p_{i}^{\mu}}}$ , are used. In Baikov representation though, the dependence on external kinematical invariants and internal masses is explicit. Indeed in Eq. (3), it is easy to identify two terms that depends on the external kinematics and/or masses, namely the overall factor $G\left(p_{1},\ldots,p_{E}\right)^{\left({-d+E+1}\right)/2}$ and the Baikov polynomial $P_{N}^{L}$ itself. The differentiation of the first factor causes no problem since the result is expressed in terms of the original integral. For the Baikov polynomial this is not so, since the derivative introduces a different integrand that is not directly expressible in terms of FI, Eq. (1). To be more specific, let us denote by $X$ a generic kinematical variable, for instance a Mandelstam invariant $X=\left(p_{i}+p_{j}\right)^{2}$ or an internal mass $X=m_{i}^{2}$ . Then

[TABLE]

where $G$ is used for ${G\left({{p_{1}},...,{p_{E}}}\right)}$ . Based on the fact that the derivatives $\frac{{\partial P_{N}^{L}}}{{\partial X}},\frac{{\partial P_{N}^{L}}}{{\partial{x_{a}}}}$ are polynomials in $x_{a}$ , the idea is to turn the derivative with respect to $X$ into derivatives with respect to $x_{a}$ . This can be achieved by the equation, known as the syzygy equation Larsen:2015ped ; Georgoudis:2016wff ,

[TABLE]

with $b$ and $c_{a}$ being polynomials in $x_{a}$ .

Assuming that a solution of this equation has been found such that $b$ is independent of $x_{a}$ (eventually depending on external kinematics and internal masses and not identical to zero), we have

[TABLE]

Then integrating by parts the second term in the rhs of the above equation and assuming that surface terms are vanishing (a standard assumption through Baikov representation) we get

[TABLE]

The term in the curly bracket is easily seen to be a sum of terms of the form

[TABLE]

The powers $\alpha^{\prime}_{a}$ depend on the actual form of the solution of the syzygy equation, Eq. (16). The result is as expected

[TABLE]

with coefficients $R_{i}$ that are rational functions of the space-time dimension $d$ , the external kinematics and the internal masses. The rhs of the above equation contains integrals that are in general not MI. We have verified in numerous examples, that after applying a standard IBP reduction to MI for the rhs of the above equation, the resulting differential equations for the MI are the same as those obtained with the standard approach. It is still interesting to note that the initial form, Eq. (19), is generally not.

5 Cutting Feynman Integrals

Cutting FI in the Baikov representation has a very natural definition. Indeed we define an $n-$ cut as follows

[TABLE]

where the Baikov variables $\{x_{a}:a=1,...,N\}$ have been divided in two subsets, containing $n$ cut propagators and $(N-n)$ uncut ones. The cut operation defined above is operational in any space-time dimension $d$ and for any FI given by Eq. (1). Notice that the definition of the cut, Eq. (20), is not identical to the traditional unitarity cut, see for instance Section 8.4 of ref. Veltman:1994wz , due to the lack of the $\theta$ -function constraint on the energy, and therefore it is not directly related to the discontinuity of the FI Abreu:2014cla ; Abreu:2015zaa .

Let us now consider a set of MI, ${F_{i}}\equiv{F_{\alpha_{1}^{\left(i\right)}...\alpha_{N}^{\left(i\right)}}}$ , $i=1,...,I$ , satisfying a system of DE, with respect to variables $X_{j}$ ,

[TABLE]

with matrices $M^{(j)}$ depending on kinematical variables, internal masses, and the space-time dimension, $d$ . Since the derivation of DE in Section 4 is insensitive to the cut operation, as defined in Eq. (20), we may immediately write555Care should be taken in defining the DE so that no symmetries of MI are used that may be violated by the corresponding $n-$ cut.

[TABLE]

with $F|_{n\times\text{cut}}$ representing an arbitrary $n-$ cut: in other words, the cut integrals satisfy the same DE as the uncut ones666See also ref. Anastasiou:2002yz ; Lee:2012te ; Larsen:2015ped for related considerations.. Of course for a given $n-$ cut many of the MI that are not supported on the corresponding cut vanish identically. Nevertheless, Eq. (22) remains valid. Especially for the maximally cut integrals defined so that $n$ is equal to the number of propagators (with $\alpha_{i}>0$ ) of the integral, all integrals not supported on the cut vanish and the resulting DE is restricted to its homogeneous part. Evaluating the maximally cut MI provides therefore a solution to the homogeneous equation Henn:2014qga ; Primo:2016ebd . Non-maximally cut integrals, on the other hand, can resolve non-homogenous parts of the DE as well Henn:2014qga .

One important implication is that cut and uncut integrals, although very different in many respects, as for instance their structure in $\epsilon-$ expansion ( $\epsilon\equiv(4-d)/2$ ), they are expressed in terms of the same class of functions777See also related discussion in ref. CaronHuot:2012ab , section 3.4.1.. This is particularly important if we want to know a priori if a system of DE can be solved, for instance, in terms of Goncharov Polylogarithms, or if the solution contains a larger class of functions including, for instance Elliptic Integrals.

In Appendix B we have collected several results of maximally cut MI: in B.1 we study a double-box with a massive loop, and find that it is expressible in terms of Polylogarithmic functions; in B.2 we study two-loop sunset graph and show that only the fully massive one is elliptic888By elliptic we mean that it is expressible only in terms of Elliptic Integrals.; in B.3 we show how results can be obtained beyond $d=4$ and verified that certain maximally cut integrals as well as certain combinations of maximally cut integrals are expressed in terms of GPs in exactly the same way as their uncut counterparts; in B.4 we show that the cut of the elliptic box-triangle integral Bonciani:2016qxi is elliptic as well, and finally in B.5 we study the elliptic double-box Bonciani:2016qxi , and verify that the Elliptic Integrals only enter through its sub-topologies.

6 Discussion and Outloook

In this paper we have studied properties of Feynman integrals in Baikov representation. We have shown how to determine the limits of integration and how to obtain DE with respect to external kinematics and internal masses. We have introduced also a loop-by-loop approach in constructing the Baikov representation, so that for certain FI a smaller number of integration variables is obtained. Then we provided a definition of a cut integral, operational in $d$ dimensions, and show that a cut integral satisfies the same system of DE as the uncut, original integral. We have shown how to compute the simplest, i.e. maximally, cut integral in Baikov representation, and give explicit results for several cases.

Based on the fact that cut integrals satisfy the same system of DE as the full, uncut integrals we have verified that their analytic expressions are given in terms of the same class of functions, such as Goncharov Polylogarithms or Elliptic Integrals. We have therefore arrived at the conclusion that in a family of MI satisfying a given system of DE, the study of the maximally cut integrals for all its members can provide a necessary and sufficient criterion for the existence of a canonical form of the DE, and in the case when such a canonical form does not exist, it provides solutions of the homogeneous parts of the system of DE (see also ref. Henn:2014qga ; Primo:2016ebd ). An application of these ideas to non-planar pentabox integrals will be discussed elsewhere.

Baikov representation is well suited for these considerations, drastically simplifying the computation of cut integrals for arbitrary external momenta and internal masses. It is still an open question if it can also be used to actually compute the MI. To this end, an algorithm, allowing the resolution of singularities in $\epsilon$ , needs to be devised. It remains to be seen if this is possible and more importantly what kind of integral representations for the individual terms in this expansion such an algorithm produces.

Acknowledgements

We gratefully acknowledge fruitful discussions with P. A. Baikov, V. A. Smirnov, M. Czakon, R. N. Lee, E. Panzer, C. Anastasiou, K. J. Larsen, J. Henn and C. Duhr during various stages of this project. This research was supported through the Initial Training Network HiggsTools under contract PITN-GA-2012-316704.

Appendix A Alternative derivation of the Baikov representation

In this appendix we will show an alternative derivation of the Baikov representation of Feynman integrals, complementary to the one given in section 2 in the main text.

We start from Eq. (1) for the one-loop case

[TABLE]

The integrand $(D_{1}^{\alpha_{1}}\cdots D_{N}^{\alpha_{N}})^{-1}$ depends on $E$ independent external momenta, which allows us to split the integration into an $E$ -dimensional “parallel” subspace $k_{||}$ and a $(d-E)$ -dimensional “orthogonal” subspace $k_{\perp}$ . Notice that the integrand depends on the orthogonal directions only through $k_{\perp}^{2}$ , allowing us to perform the angular part of the integral over the orthogonal space

[TABLE]

where we define $\lambda\equiv k_{\perp}^{2}$ .

We change now to a different set of variables $\varsigma_{i}\equiv k\cdot p_{i}$ which is equivalent to $s_{1,E+i}$ from section 2. Defining $G$ as the Gram determinant of the external momenta, it is not hard to realize that $d^{E}k_{||}=-(-G)^{-1/2}d^{E}\varsigma$ , and we now want to change the integration to the set of Baikov variables $x$ . Using the same argument as in section 2, we get that $d^{E}k_{||}d\lambda=\det(A^{-1})d^{E+1}x$ where $A$ is the matrix defined in Eq. (2), and we see that $\det(A^{-1})=\pm 2^{-E}$ depending on the sign of the external momenta in the definition of the loop momenta. Putting this together, yields

[TABLE]

If additionally we use that $\lambda=P/G$ where $P$ – the Baikov polynomial – is given as the Gram determinant of the full set of momenta $\{q\}=\{k,p_{1},\ldots,p_{E}\}$ , we see that Eq. (25) is equivalent to Eq. (3) in the one-loop case. The above derivation can be straightforwardly generalised for any numbers of loops.

The loop-by-loop case

We may now go through the same procedure as above, loop by loop, for a two-loop process. A traditional Baikov representation of a two-loop Feynman integral with $E$ independent external momenta would involve $3+2E$ integrations, but in most cases it is possible to get the number further down. Each individual loop will in general be dependent only on a subset of the $E$ external momenta. Let us denote with $k_{2}$ the loop momentum of the loop with the smallest such subset, and the size of that subset $E_{2}-1$ . This means that the integrand of this loop depends on $E_{2}$ different momenta including $k_{1}$ . Performing the transverse angular part of the $k_{2}$ integration individually with the method leading up to Eq. (24), gives

[TABLE]

where $\lambda_{22}={k_{2}}_{\perp}^{2}$ .

The integrand now depends on all the external momenta, so the transverse angular part of the $k_{1}$ integration may be done, leaving $E+E_{2}+2$ integrations in total. So only in the cases where both loops depend on $E$ external momenta (so $E_{2}=E+1$ ) the number of integrations will be equal to that of the standard Baikov approach, otherwise it will be smaller.

Changing integration variables from $k_{||}$ to $\varsigma$ and from $\varsigma$ and $\lambda$ to the Baikov variables $x$ , our result for the two-loop case in the loop-by-loop approach is written as

[TABLE]

with $G_{1}$ the Gram determinant of the $E$ external momenta and $G_{2}$ the Gram determinant of the $E_{2}$ different momenta including $k_{1}$ .

Appendix B Examples of cuts performed in Baikov variables

In the appendix we will go through several examples of maximally cut Feynman integrals.

B.1 Double-box with massive loop

As a first example we will do a double-box with one massive loop. That integral is given as

[TABLE]

with

[TABLE]

where $x_{8}$ is needed for later. Using the loop-by-loop approach, Eq. (30), we may perform the angular parts of the integral, leaving integrals over the eight Baikov variables

[TABLE]

The various quantities are999This is the only example in which these quantities will be written out in full

[TABLE]

and $s={\left({{p_{1}}+{p_{2}}}\right)^{2}},t={\left({{p_{1}}+{p_{4}}}\right)^{2}},s>0,t<0,s+t>0,{m^{2}}>0$ .

Inserting all this, cutting the seven propagators, expanding in $\epsilon$ , and renaming $x_{8}$ to $z$ , gives

[TABLE]

where

[TABLE]

and $r_{\mp}$ are the two roots of $P(z)$ . The result is

[TABLE]

In fact we have computed also the order $\epsilon$ which is given in terms of weight 2 functions, but since it is quite extensive we give the result for another MI, defined as follows,

[TABLE]

whose expression is simpler,

[TABLE]

suggesting that the solution is expressible in terms of logarithmic/polylogarithmic functions.

B.2 The sunset

We will here go through the same considerations for the well-studied sunset-integral Adams:2015gva . The integral is given as

[TABLE]

with

[TABLE]

Using the loop-by-loop approach, we get

[TABLE]

with

[TABLE]

$\lambda(s_{1},s_{2},s_{3})\equiv s_{1}^{2}+s_{2}^{2}+s_{3}^{2}-2s_{1}s_{2}-2s_{2}s_{3}-2s_{3}s_{1}$ , is the Källén function and $s=p^{2}$ . Performing the triple cut gives Remiddi:2016gno ,

[TABLE]

Expanding in $\epsilon$

[TABLE]

it is easily seen that this is an integral over the square-root of a quartic polynomial Remiddi:2016gno , which yields elliptic integrals, and such functions are indeed present in the result for the full sunset integral as well. On the other hand, if any of the masses vanish, the square root partially factorises, and the result is expressible in terms of polylogarithmic functions.

Notice, that the extra cut over $x_{4}$ ( $m_{2}>m_{3}$ ), after expanding in $\epsilon$ ,

[TABLE]

yields a rational term at order $\mathcal{O}(\epsilon^{0})$ , which seems to contradict our previous findings, Eq. (45), as well as the known result for the full integral Adams:2015gva . Nevertheless, it is easy to verify, starting from Eq. (42), that $F_{\text{sunset}|4\times\text{cut}}=0$ and that the result of Eq. (46) is just an artefact of expanding in $\epsilon$ before cutting.

B.3 The bubble-triangle

As mentioned in Section 1, the introduction of the idea of canonical bases Henn:2013pwa caused a renaissance in the derivation of analytical expressions of Feynman integrals. We will not go through the details of what a canonical basis implies, merely say that no general algorithm for obtaining such a basis exists. In ref. Henn:2014qga it was implied that canonical integrals have maximal cuts which equal a number (rather than a function of kinematical variables). This criterion significantly limits the search space, and has been a helpful guideline in many cases.

Yet for some cases the maximal cut derived in the traditional way in four dimensions, does not provide the full information obtainable from such a criterion. For an example lets us re-examine the triangle with bubble insertion discussed in section 3. We will here parametrize it as

[TABLE]

with

[TABLE]

Using the loop-by-loop approach of eq. (30), we may turn this into a five-fold integral

[TABLE]

The set of canonical differential equations used in ref. Henn:2014lfa contains two integrals with the topology of this integral, and they are given as ( $s={\left({{p_{1}}+{p_{2}}}\right)^{2}}$ )

[TABLE]

where $R_{12}$ denotes the Källén function

[TABLE]

Performing the cut in four dimensions will only catch the leading term in $\epsilon$ , and thus not the second term in $I_{2}$ . On the contrary, the definition of maximal cut given in Eq. (20), allows us to compute the cut integral in $d$ dimensions and for arbitrary powers of propagators in the integrand. The result is given by

[TABLE]

where $z=x_{5}$ and $r_{\pm}$ the roots of the polynomial: $p_{1}^{2}p_{2}^{2}-p_{1}^{2}z-p_{2}^{2}z+sz+z^{2}$ . These integrals evaluate to hypergeometric functions, and as the second integral diverges in the $d\rightarrow 4$ limit, the integration does not commute with the $\epsilon$ expansion. Thus one has to evaluate them in $d$ dimensions and then expand the resulting functions, for instance using the HypExp package Huber:2007dx , to the desired order in $\epsilon$ . The result for the integrals defined in Eq. (50) is given by

[TABLE]

where the variables $x$ , $y$ have been introduced, defined through $s\equiv p_{1}^{2}(1+x)(1+xy)$ and $p_{2}^{2}\equiv p_{1}^{2}x^{2}y$ .

In order to explicitly show the uniform transcendentallity property of these functions, we present the results expanded up to order $\epsilon^{3}$ ,

[TABLE]

with $N_{\epsilon}=e^{2\gamma_{E}\epsilon}(p_{1}^{2})^{\epsilon}x^{\epsilon}(x+1)^{\epsilon}(xy+1)^{\epsilon}$ .

B.4 The elliptic box-triangle

Here we will go through the generalized cut of the elliptic box-triangle studied in ref. Bonciani:2016qxi where it is denoted $I^{A}_{66}$ , and in ref. Primo:2016ebd . It is given as

[TABLE]

with

[TABLE]

where $x_{7}$ has been introduced for later convenience. The kinematics is such that

[TABLE]

With the loop-by-loop approach as described in Section 3, we see that the $k_{2}$ integral has a transverse space with dimension $2$ , as it is formed by $p_{4}$ and $k_{1}$ , while the $k_{1}$ integral has a transverse space formed by $p_{1}$ , $p_{2}$ , and $p_{4}$ , which has dimension $3$ . Thus we can write $F_{\text{box-triangle}}$ as a $7$ -dimensional integral over the seven Baikov variables given above.

[TABLE]

Here $\lambda_{22}$ is a function of all the Baikov variables, and $\lambda_{11}$ and $G_{2}$ of those not containing $k_{2}$ . We may then do the hexa-cut of the integral, yielding an integral over one variable $z=x_{7}$ . That integral is well-behaved in the $d=4$ limit, and thus the integration and the $\epsilon$ expansion commute, giving

[TABLE]

where $C=i/(4\pi^{3})$ and

[TABLE]

and $r_{\mp}$ are the two roots of $F_{2}(z)$ which are given as

[TABLE]

Following Byrd:1971 , we may perform the integral with the result

[TABLE]

where $K(k^{2})$ is the complete elliptic integral of the first kind, and where

[TABLE]

which holds in the physical region with $m^{2}>0$ , $s>0$ , $p_{4}^{2}>0$ , $t<0$ , $u<0$ .

B.5 The elliptic double-box

Here we will go through the generalized cut of the elliptic double-box studied in ref. Bonciani:2016qxi where it is denoted $I^{A}_{70}$ . It is given as

[TABLE]

with

[TABLE]

The kinematics is as in the previous example.

Here each of the two integrals over the loop momenta $k_{1}$ and $k_{2}$ are over transverse spaces with dimensionality three, giving an eight-dimensional integral after the angular components have been integrated out:

[TABLE]

As before we may perform a cut of all the propagators – an hepta-cut in this case, yielding an integral over the last variable variable $z=x_{8}$

[TABLE]

where $C=1/(4\pi^{4})$ and

[TABLE]

and $r_{\mp}$ are the two roots of $f(z)$ which are given as

[TABLE]

That integral may be done and yields

[TABLE]

verifying that the elliptic character of the uncut integral stems from the non-homogeneous part of the corresponding differential equation through its elliptic sub-topologies Bonciani:2016qxi .

Bibliography80

The reference list from the paper itself. Each links out to its DOI / PubMed record.

1(1) R. P. Feynman, Space - time approach to quantum electrodynamics , Phys. Rev. 76 (1949) 769–789 . · doi ↗
2(2) F. J. Dyson, The Radiation theories of Tomonaga, Schwinger, and Feynman , Phys. Rev. 75 (1949) 486–502 . · doi ↗
3(3) F. J. Dyson, The S matrix in quantum electrodynamics , Phys. Rev. 75 (1949) 1736–1755 . · doi ↗
4(4) G. ’t Hooft and M. J. G. Veltman, Regularization and Renormalization of Gauge Fields , Nucl. Phys. B 44 (1972) 189–213 . · doi ↗
5(5) T. Gehrmann, J. M. Henn and N. A. Lo Presti, Analytic form of the two-loop planar five-gluon all-plus-helicity amplitude in QCD , Phys. Rev. Lett. 116 (2016) 062001 , [ 1511.05409 ]. · doi ↗
6(6) C. G. Papadopoulos, D. Tommasini and C. Wever, The Pentabox Master Integrals with the Simplified Differential Equations approach , JHEP 04 (2016) 078 , [ 1511.09404 ]. · doi ↗
7(7) S. Borowka, G. Heinrich, S. P. Jones, M. Kerner, J. Schlenk and T. Zirke, Sec Dec-3.0: numerical evaluation of multi-scale integrals beyond one loop , Comput. Phys. Commun. 196 (2015) 470–491 , [ 1502.06595 ]. · doi ↗
8(8) J. R. Andersen et al., Les Houches 2013: Physics at Te V Colliders: Standard Model Working Group Report , 1405.1067 .