Baikov-Lee Representations Of Cut Feynman Integrals

TL;DR

This paper introduces a new framework using Baikov-Lee representation for evaluating cut Feynman integrals, enabling better understanding of their discontinuities and solutions beyond multiple polylogarithms.

Contribution

It develops a general method for cut Feynman integrals using Baikov-Lee representation, applying residue calculus, and solving differential equations beyond polylogarithms.

Findings

Reproduces the relation between cuts and discontinuities.

Provides complete solutions to differential equations for Feynman integrals.

Enables direct evaluation of complex cut Feynman integrals.

Abstract

We develop a general framework for the evaluation of $d$-dimensional cut Feynman integrals based on the Baikov-Lee representation of purely-virtual Feynman integrals. We implement the generalized Cutkosky cutting rule using Cauchy's residue theorem and identify a set of constraints which determine the integration domain. The method applies equally well to Feynman integrals with a unitarity cut in a single kinematic channel and to maximally-cut Feynman integrals. Our cut Baikov-Lee representation reproduces the expected relation between cuts and discontinuities in a given kinematic channel and furthermore makes the dependence on the kinematic variables manifest from the beginning. By combining the Baikov-Lee representation of maximally-cut Feynman integrals and the properties of periods of algebraic curves, we are able to obtain complete solution sets for the homogeneous differential equations satisfied by Feynman integrals which go beyond multiple polylogarithms. We apply our formalism to the direct evaluation of a number of interesting cut Feynman integrals.

Equations

$$I_{\mathrm{E}}=\int {\mathrm{d}}^d{k}_1\cdots \int {\mathrm{d}}^d{k}_L\prod _{\ell =1}^N{\left(Q_{\ell }^2(k_i,p_j)+m_{\ell }^2\right)}^{-{\nu }_{\ell }},$$

$$\displaystyle I_{\rm E}=\frac{\pi^{\frac{Ld}{2}}\Gamma\big{(}\nu-\frac{Ld}{2}\big{)}}{\prod_{i=1}^{N}\Gamma(\nu_{i})}\Bigg{[}\prod_{j=1}^{N}\int_{0}^{\infty}{\rm d}x_{j}\Bigg{]}\delta(1-x_{N})\,\mathcal{U}_{\rm E}^{\nu-(L+1)d/2}\mathcal{F}_{\rm E}^{Ld/2-\nu}\prod_{k=1}^{N}x_{k}^{\nu_{k}-1},$$

$$I_{\mathrm{M}}=\int {\mathrm{d}}^d{k}_1\cdots \int {\mathrm{d}}^d{k}_L\prod _{\ell =1}^N{\left(Q_{\ell }^2(k_i,p_j)-m_{\ell }^2+i0\right)}^{-{\nu }_{\ell }}$$

$$\displaystyle I_{\rm M}=\frac{i^{L}\pi^{\frac{Ld}{2}}e^{-i\pi\nu}\Gamma\big{(}\nu-\frac{Ld}{2}\big{)}}{\prod_{i=1}^{N}\Gamma(\nu_{i})}\Bigg{[}\prod_{j=1}^{N}\int_{0}^{\infty}{\rm d}x_{j}\Bigg{]}\delta(1-x_{N})\,\mathcal{U}_{\rm M}^{\nu-(L+1)d/2}\mathcal{F}_{\rm M}^{Ld/2-\nu}\prod_{k=1}^{N}x_{k}^{\nu_{k}-1},$$

$$\displaystyle I_{\rm M}=i^{L}e^{-i\pi\nu}I_{\rm E}\bigg{|}_{\omega_{i}\to-\omega_{i},\,\{p_{j}^{*}\}\to-\{p_{j}^{*}\}},$$

$$I=\int \frac{{\mathrm{d}}^d{k}_1}{i{\pi }^{d/2}}\cdots \int \frac{{\mathrm{d}}^d{k}_L}{i{\pi }^{d/2}}\prod _{\ell =1}^N{\left(Q_{\ell }^2(k_i,p_j)-m_{\ell }^2+i0\right)}^{-{\nu }_{\ell }}$$

$$G(q_1,\dots ,q_K)=\left|\begin{array}{ccc}{q}_1^2& \cdots & q_1\cdot q_K\\ & & \\ ⋮& \ddots & ⋮\\ & & \\ q_1\cdot q_K& \dots & q_K^2\end{array}\right|.$$

$$I_{\mathrm{E}}$$

$$\times {\left[G(q_1,\dots ,q_{n+L-1})\right]}^{(d-n-L)/2}\prod _{\ell =1}^N{\left(Q_{\ell }^2(q_i\cdot q_j)+m_{\ell }^2\right)}^{-{\nu }_{\ell }},$$

$$Q_1=k_1{Q}_2=k_1-p.$$

$$\displaystyle\vbox{\hbox{\includegraphics[scale={.3}]{virtoneloopbubble}}}=\underset{\mathcal{D}}{\int}\int{\rm d}(q_{1}^{2}){\rm d}(q_{1}\cdot q_{2})\frac{\pi^{\frac{3}{2}-\epsilon}\left(p^{2}q_{1}^{2}-(q_{1}\cdot q_{2})^{2}\right)^{\frac{1}{2}-\epsilon}}{\Gamma\left(\frac{3}{2}-\epsilon\right)\left(p^{2}\right)^{1-\epsilon}q_{1}^{2}\left(q_{1}^{2}-2q_{1}\cdot q_{2}+p^{2}\right)},$$

$$\displaystyle\vbox{\hbox{\includegraphics[scale={.3}]{virtoneloopbubble}}}=\int_{0}^{\infty}{\rm d}(q_{1}^{2})\int_{-\sqrt{p^{2}q_{1}^{2}}}^{\sqrt{p^{2}q_{1}^{2}}}{\rm d}(q_{1}\cdot q_{2})\frac{\pi^{\frac{3}{2}-\epsilon}\left(p^{2}q_{1}^{2}-(q_{1}\cdot q_{2})^{2}\right)^{\frac{1}{2}-\epsilon}}{\Gamma\left(\frac{3}{2}-\epsilon\right)\left(p^{2}\right)^{1-\epsilon}q_{1}^{2}\left(q_{1}^{2}-2q_{1}\cdot q_{2}+p^{2}\right)}.$$

$$I_{\mathrm{M}}$$

$$\displaystyle\qquad\quad\times\left[G(q_{1},\ldots,q_{n+L-1})\right]^{(d-n-L)/2}\prod_{\ell=1}^{N}\left(Q_{\ell}^{2}(q_{i}\cdot q_{j})+m_{\ell}^{2}\right)^{-\nu_{\ell}}\bigg{|}_{\omega_{i}\to-\omega_{i},\,\{P_{j}^{*}\}\to-\{P_{j}^{*}\}}$$

$$=\frac{i^L{\pi }^{L(3+2d-2n-L)/4}{e}^{-i\pi \nu }}{{\prod }_{r=0}^{L-1}\Gamma \left(\frac{d-n-r+1}{2}\right){\left[\bar{G}(P_1,\dots ,P_{n-1})\right]}^{(d-n)/2}}\underset{\bar{\mathcal{D}}}{\int \cdots \int }\left(\prod _{i=1}^L\prod _{j=i}^{n+L-1}\mathrm{d}(q_i\cdot q_j)\right)\times$$

$$\times {\left[\bar{G}(q_1,\dots ,q_{n+L-1})\right]}^{(d-n-L)/2}\prod _{\ell =1}^N{\left({\bar{Q}}_{\ell }^2(q_i\cdot q_j)+m_{\ell }^2\right)}^{-{\nu }_{\ell }},$$

$$\frac{1}{k^2+i0}-\frac{1}{k^2-i0}=-2\pi i\theta (k^0)\delta \left(k^2\right),$$

$$I_{{\omega }_i-\mathrm{cut}}^{(j)}=\frac{(-2\pi i{)}^{n_j}{\pi }^{L(3-2n-L)/4}{e}^{-i\pi \nu }}{{\prod }_{r=0}^{L-1}\Gamma \left(\frac{d-n-r+1}{2}\right){\left[\bar{G}(P_1,\dots ,P_{n-1})\right]}^{(d-n)/2}}\underset{{\bar{\mathcal{D}}}_{{\omega }_i-\mathrm{cut}}^{(j)}}{\int \cdots \int }\left(\prod _{k=1}^{L(L-1)/2+nL-n_j}\mathrm{d}{s}_k\right)\times$$

$$\times \mathbf{sgn}\left(\left|\frac{\partial \{{\bar{Q}}_{\ell }^2{\}}_{{\omega }_i-\mathrm{cut}}^{(j)}}{\partial \{{\bar{s}}_i\}}\right|\right)\underset{\{{\bar{s}}_i\}}{\mathrm{Res}}\left\{{\left[\bar{G}(q_1,\dots ,q_{n+L-1})\right]}^{(d-n-L)/2}\prod _{\ell =1}^N{\left({\bar{Q}}_{\ell }^2({\bar{s}}_i,s_k)+m_{\ell }^2\right)}^{-{\nu }_{\ell }}\right\},$$

$$\underset{\{a\}}{\mathrm{Res}}\left\{\frac{f(z)}{a-z}\right\}=-f(a),$$

$${\int }_{-\infty }^{\infty }\mathrm{d}z\delta (a-z)f(z)={\int }_{-\infty }^{\infty }\mathrm{d}z\delta (z-a)f(z)=f(a)$$

$${\mathrm{Disc}}_{{\omega }_i}\left(I\right)=-\sum _{j=1}^M{I}_{{\omega }_i-\mathrm{cut}}^{(j)}.$$

$${\hat{I}}_{{\omega }_i-\mathrm{cut}}^{(j)}=\frac{i^L(2\pi {)}^{n_j}{\pi }^{L(3-2n-L)/4}{e}^{-i\pi \nu }}{{\prod }_{r=0}^{L-1}\Gamma \left(\frac{d-n-r+1}{2}\right){\left[\bar{G}(P_1,\dots ,P_{n-1})\right]}^{(d-n)/2}}\underset{{\bar{\mathcal{D}}}_{{\omega }_i-\mathrm{cut}}^{(j)}}{\int \cdots \int }\left(\prod _{k=1}^{L(L-1)/2+nL-n_j}\mathrm{d}{s}_k\right)\times$$

$$\times \mathbf{sgn}\left(\left|\frac{\partial \{{\bar{Q}}_{\ell }^2{\}}_{{\omega }_i-\mathrm{cut}}^{(j)}}{\partial \{{\bar{s}}_i\}}\right|\right)\underset{\{{\bar{s}}_i\}}{\mathrm{Res}}\left\{{\left[\bar{G}(q_1,\dots ,q_{n+L-1})\right]}^{(d-n-L)/2}\prod _{\ell =1}^N{\left({\bar{Q}}_{\ell }^2({\bar{s}}_i,s_k)+m_{\ell }^2\right)}^{-{\nu }_{\ell }}\right\}$$

$${\mathrm{Disc}}_{{\omega }_i}\left(\mathcal{P}\left\{\hat{I}\right\}\hat{I}\right)=\sum _{j=1}^M\mathcal{P}\left\{{\hat{I}}_{{\omega }_i-\mathrm{cut}}^{(j)}\right\}{\hat{I}}_{{\omega }_i-\mathrm{cut}}^{(j)},$$

$$\displaystyle\vbox{\hbox{\includegraphics[scale={.3}]{virtonelooptriangle}}}=\int\frac{d^{d}k_{1}}{i\pi^{d/2}}~{}\frac{1}{[(p_{2}-k_{1})^{2}]^{\nu_{1}}[(p_{1}+k_{1})^{2}]^{\nu_{2}}[k_{1}^{2}]^{\nu_{3}}}.$$

$$=\frac{(-2\pi i{)}^2{\pi }^{-1}(-1{)}^{\nu }}{\Gamma (1-ϵ){\left(-s^2/4\right)}^{1/2-ϵ}}\mathbf{sgn}\left(\left|\frac{\partial \{s_1-2{\bar{s}}_1,s_1+2{\bar{s}}_2\}}{\partial \{{\bar{s}}_1,{\bar{s}}_2\}}\right|\right)\times$$

$$\times \underset{{\bar{\mathcal{D}}}_{s-\mathrm{cut}}}{\int }\mathrm{d}{s}_1\underset{\{{\bar{s}}_1,{\bar{s}}_2\}}{\mathrm{Res}}\left\{\frac{{\left(-{\bar{s}}_1{\bar{s}}_{2}s-s_1{s}^2/4\right)}^{-ϵ}}{{\left(s_1-2{\bar{s}}_1\right)}^{{\nu }_1}{\left(s_1+2{\bar{s}}_2\right)}^{{\nu }_2}{s}_1^{{\nu }_3}}\right\}$$

$$=\frac{2^{3-{\nu }_1-{\nu }_2}\pi (-1{)}^{{\nu }_2+{\nu }_3}i}{s^{1-ϵ}\Gamma (1-ϵ)}\times$$

$$\times \underset{{\bar{\mathcal{D}}}_{s-\mathrm{cut}}}{\int }\frac{\mathrm{d}{s}_1}{s_1^{{\nu }_3}}\underset{\{{\bar{s}}_1,{\bar{s}}_2\}}{\mathrm{Res}}\left\{\frac{{\left(4{\bar{s}}_1{\bar{s}}_2+s_{1}s\right)}^{-ϵ}}{{\left({\bar{s}}_1-s_1/2\right)}^{{\nu }_1}{\left({\bar{s}}_2+s_1/2\right)}^{{\nu }_2}}\right\}$$

$$=-\frac{2^{1+{\nu }_1-{\nu }_2}\pi \Gamma ({\nu }_1-1+ϵ)(-1{)}^{\nu }i}{s^{1-ϵ}\Gamma (ϵ)\Gamma (1-ϵ)\Gamma ({\nu }_1)}\times$$

$$\times \underset{{\bar{\mathcal{D}}}_{s-\mathrm{cut}}}{\int }\frac{\mathrm{d}{s}_1}{s_1^{{\nu }_1+{\nu }_3-1+ϵ}}\underset{\{{\bar{s}}_2\}}{\mathrm{Res}}\left\{\frac{{\bar{s}}_2^{{\nu }_1-1}{\left(s+2{\bar{s}}_2\right)}^{1-{\nu }_1-ϵ}}{{\left({\bar{s}}_2+s_1/2\right)}^{{\nu }_2}}\right\}$$

$$=-\frac{2\pi \Gamma ({\nu }_1-1+ϵ)(-1{)}^{{\nu }_3}i}{s^{2-{\nu }_2-ϵ}\Gamma (ϵ)\Gamma (1-ϵ)\Gamma ({\nu }_2)\Gamma (1+{\nu }_1-{\nu }_2)}\underset{{\bar{\mathcal{D}}}_{s-\mathrm{cut}}}{\int }\frac{\mathrm{d}{s}_1}{s_1^{{\nu }_2+{\nu }_3-1+ϵ}}\times$$

$$\times (s-s_1{)}^{2-{\nu }_1-{\nu }_2-ϵ}{}_2{F}_1\left(1-{\nu }_2,2-{\nu }_2-ϵ;1+{\nu }_1-{\nu }_2;\frac{s_1}{s}\right).$$

$$=-\frac{2\pi \Gamma ({\nu }_1-1+ϵ)(-1{)}^{{\nu }_3}i}{s^{2-{\nu }_2-ϵ}\Gamma (ϵ)\Gamma (1-ϵ)\Gamma ({\nu }_2)\Gamma (1+{\nu }_1-{\nu }_2)}{\int }_0^s\frac{\mathrm{d}{s}_1}{s_1^{{\nu }_2+{\nu }_3-1+ϵ}}\times$$

$$\times (s-s_1{)}^{2-{\nu }_1-{\nu }_2-ϵ}{}_2{F}_1\left(1-{\nu }_2,2-{\nu }_2-ϵ;1+{\nu }_1-{\nu }_2;\frac{s_1}{s}\right).$$

Full text

TCDMATH-17-12

1 Introduction

As has long been known, the evaluation of multi-loop Feynman integrals is an important component of many high-precision collider physics calculations. Even at the two-loop level, both purely-virtual and cut Feynman integrals often provide a remarkable computational challenge. This is particularly true if one proceeds analytically, and, as a consequence, a number of specialized techniques have been developed to aid in the evaluation of such integrals. Most direct analytic integration techniques have traditionally employed some variant of the well-known Feynman (or Schwinger) parametric representation Feynman:1949zx as a starting point in the purely-virtual case, due to the fact that it makes the dependence of the integral on Lorentz-invariant quantities manifest. The aim of this paper is to supply an analogous framework for cuts which is suitable not only for the direct calculation of real radiative master integrals, but also for the maximally-cut integrals relevant to the study of purely-virtual Feynman integrals in the method of differential equations Kotikov:1990kg ; Kotikov:1991hm ; Kotikov:1991pm ; Bern:1992em ; Bern:1993kr ; Remiddi:1997ny ; Gehrmann:1999as . As we shall see, our work is a natural out-growth of earlier work on the subject of integral reduction which we apply in a novel way to evaluate various types of cut Feynman integrals.

It is unfortunately the case that phenomenologically-relevant multi-loop integration by parts reductions Tkachov:1981wb ; Chetyrkin:1981qh often require dedicated effort and substantial resources to compute. Although it is probably fair to say that most recent higher-order perturbative calculations rely on some variant of Laporta’s algorithm Laporta:2001dd ; vonManteuffel:2012np ; vonManteuffel:2014ixa ; Smirnov:2014hma , several other interesting and useful algorithms have been worked out and implemented over the years (see e.g. Gorishnii:1989gt ; Baikov:1996iu ; Gluza:2010ws ; Lee:2012cn ; Ruijl:2015aca ; Ita:2015tya ; Larsen:2015ped ; Ueda:2016sxw ; Abreu:2017idw ; Abreu:2017xsl ). However, even very different approaches, such as the one advocated by Baikov Baikov:1996iu , sometimes turn out to have ramifications for Laporta’s method as well. Some time ago, the derivation of the Baikov formula for purely-virtual Feynman integrals was clarified by Lee and then used to great effect as a generating function for integration by parts relations Lee:2010wea . The product of his analysis, what we shall hereafter refer to as the Baikov-Lee representation, will be of great interest to us in this work. Larsen and Zhang argued in Larsen:2015ped that evaluating the integration by parts relations generated by the Baikov representation on the support of various cuts dramatically improves the approach to integral reduction originally advocated by Gluza, Kajda, and Kosower Gluza:2010ws . In fact, from their work and a similar study by Ita Ita:2015tya , one can readily guess that the Baikov-Lee representation ought to offer a useful starting point for the evaluation of cut Feynman integrals.

Even though cut Feynman integrals play every bit as important a role as purely-virtual Feynman integrals in the perturbative computation of collider observables, methods for their direct evaluation have not been as thoroughly developed.111Certainly, we do not wish to suggest that such techniques do not exist. For instance, reference Anastasiou:2013srw describes some pertinent traditional and state-of-the-art direct integration techniques in great detail. In fact, a setup which bears at least some rudimentary resemblance to the one discussed in this paper was developed for tree-level cross section calculations in $d=4$ long ago Kumar:1970cr . It is entirely possible that, at least in part, this state of affairs has persisted due to the fact that it is not completely trivial to write down a Baikov representation for Feynman integrals directly in Minkowski space. The issue is that the standard derivation of the Baikov-Lee representation (see e.g. Grozin:2011mt for a detailed exposition) relies heavily upon Euclidean geometric intuition which does not immediately generalize. In this paper, we write down a simple recipe for the analytical continuation of the Euclidean Baikov-Lee formula by drawing an analogy to the more familiar situation that one encounters in the derivation of the Feynman parametric representation. To pass from uncut to cut propagators, we use sequential applications of Cauchy’s residue theorem to implement a natural generalization of Cutkosky’s cutting rule Cutkosky:1960sp ; Lee:2012te .222Due to the fact that the generalized Cutkosky rule of reference Lee:2012te is actually written in terms of delta distributions and their derivatives, certain subtleties apply. To be consistent, one must first integrate out a complete set of scalar products localized by the distributions before attempting any of the more non-trivial integrations over scalar products. We would like to thank Ruth Britto for emphasizing this point to us. Finally, the integration region is determined by analyzing the analytical structure of the integrand and applying the available constraints from physics. At each step of the calculation, one integrates between branch points of the integrand with respect to the current variable of integration. This procedure allows one to write down Baikov-Lee representations for a wide class of cut Feynman integrals. Crucially, this approach makes the dependence of the cut integral on the kinematic invariants of the problem manifest and eliminates the need to set up a convenient reference frame to carry out integrations over the components of the cut loop momenta.333To appreciate this point, we invite the reader to compare and contrast the maximally-cut Feynman integral calculation which appears in both references Primo:2016ebd and Frellesvig:2017aai .

As mentioned above, the Baikov-Lee representation for cut Feynman integrals also has important applications to the more indirect method of differential equations. Writing differential equations with respect to the available kinematic parameters allows for the complete determination of a large class of Feynman integrals in terms of an appropriate set of iterated integrals, order-by-order in the parameter of dimensional regularization, $\epsilon$. For simple families of Feynman integrals with relatively few ratios of scales, the function space associated to Feynman integrals is spanned by iterated integrals with rational integrating factors. This set of iterated integrals, comprised of the well-known multiple polylogarithms LDanilevsky , has been studied and popularized by many authors (Goncharov:1998kja ; Remiddi:1999ew ; Gehrmann:2000zt to name a few). In the polylogarithmic case, it is always possible, by suitably choosing the basis of Feynman integrals Kotikov:2010gf ; Henn:2013pwa ; Henn:2013tua , to reduce the problem to one which admits an elementary formal solution in terms of Chen iterated integrals Chen:1977oja . Using the symbol-coproduct calculus Goncharov:2005sla ; Brown:2011ik ; Duhr:2011zq ; Duhr:2012fh , performance-optimized solutions which can be readily interfaced with a Monte Carlo integration program may be constructed for the purposes of phenomenology (see e.g. Gehrmann:2015ora ; vonManteuffel:2015msa ; Bonciani:2015eua ; vonManteuffel:2017myy ). Even through to weight four, this is not always straightforward to do in practice, despite the fact that the weight four function space has been studied extensively and is in principle well-understood Frellesvig:2016ske .

For general Feynman integrals, far more complicated analytic structures may appear; even simple-looking two-loop integrals which depend on sufficiently many kinematic variables may already involve elliptic polylogarithms and related functions BrownLevin ; Bloch:2013tra ; Bloch:2014qca ; Adams:2014vja ; Adams:2015gva ; Adams:2016xah . In this context, maximally-cut Feynman integrals play an important role because they satisfy the homogeneous differential equations for the associated uncut Feynman integrals Lee:2012te ; Primo:2016ebd .444Although, we will leave a detailed discussion for future work, there is nothing stopping us from applying this technique also to the differential equations satisfied by multi-scale cut Feynman integrals in the reverse-unitarity method Anastasiou:2002yz ; Anastasiou:2003yy ; Anastasiou:2003ds . This property is particularly useful when Feynman integrals cannot be expressed in terms of multiple polylogarithms, and iterated integrals over special functions need to be considered Caffo:1998du ; Laporta:2004rb ; Adams:2015ydq ; Bloch:2016izu ; Remiddi:2016gno ; Bonciani:2016qxi ; vonManteuffel:2017hms . In fact, higher-order differential equations appear in non-polylogarithmic cases, and, at the present time, no general solution algorithm is known. Nevertheless, it was observed by Primo and Tancredi in Primo:2016ebd that, upon setting $\epsilon$ to zero, maximally-cut Feynman integrals can often be computed in closed form, allowing one to find at least a single homogeneous solution to the higher-order differential equations under consideration by direct integration. Provided that a complete, linearly independent set of homogeneous solutions can be found, the full solution can finally be determined using the variation of parameters technique. In our opinion, it is of great importance to supply an algorithm which comes up with not just one, but rather a complete set of homogeneous solutions. We provide a general prescription555A few days prior to the appearance of this paper, we became aware of a recent preprint, Primo:2017ipr , which discusses many of the same technical issues for non-polylogarithmic Feynman integrals. In fact, the maximally-cut case was also discussed recently in the Baikov approach by yet another group Bosma:2017ens , with, however, a different set of physical problems in mind. and show that it allows for the straightforward construction of complete sets of homogeneous solutions for the non-polylogarithmic examples of Section 4.

The plan of this paper is as follows. In Section 2, we explain in detail how to write down Baikov-Lee representations for Feynman integrals cut in a single kinematic channel and how to work out an explicit description of the relevant integration domain. Our focus will be on the physical case relevant to cross section calculations because the procedure for maximally-cut Feynman integrals is closely analogous and has already been discussed in reference Frellesvig:2017aai . Although we do not have a computer program which finds the integration limits for arbitrarily complicated cut Feynman integrals, we have a solid conceptual understanding which could in the future lead us to an explicit algorithm. In Section 3, we go through a number of well-studied one- and two-loop examples of Feynman integrals cut in a single kinematic channel in order to give the reader a feeling for how explicit computations typically proceed when one adopts a Baikov-Lee representation as the starting point. To the best of our knowledge, our treatment of classical cut Feynman integrals in the Baikov-Lee representation is new and effectively extends the work of Frellesvig and Papadopoulos Frellesvig:2017aai beyond the maximally-cut case. Although we work with generic values of the spacetime dimension for pedagogical purposes, our calculations strongly suggest that, in practice, an expansion in $\epsilon$ under the integral sign must be the way to go for all but the very simplest of cut Feynman integrals.

In Section 4, we move on to maximally-cut Feynman integrals which evaluate to complete elliptic integrals. We discuss a general solution strategy applicable to many problems of practical interest and then demonstrate the general procedure by focusing on examples which are suitable for exposition. We emphasize in particular the utility of integrating out one loop at a time, as this generically leads to simpler Baikov-Lee representations. Finally, we conclude in Section 5 and outline our plans for future research. We also include a number of appendices for pedagogical purposes and cross-checks. To streamline the exposition in the body of the paper, we summarize a number of purely mathematical results from the theory of hypergeometric-like functions in Appendix A. In Appendix B, we reproduce physical-region results from the literature for the uncut versions of the integrals considered in Section 3. This allows the reader to easily verify our results using the classical relation between discontinuities and cuts in a given kinematic channel (see e.g. tHooft:1973wag ; Gehrmann-DeRidder:2003pne ; Abreu:2014cla ). Finally, in Appendix C, we evaluate a simple uncut Feynman integral using the Baikov-Lee setup to help the less familiar reader understand the relation between our prescriptions for cut Feynman integrals and the usual prescriptions for Baikov’s method in the purely-virtual case.

2 General formalism

In this section, we define our notation, recall some results from the literature, and explain how we generalize the Baikov-Lee representation to the case of cut Feynman integrals.

2.1 Preliminaries

Let us begin by discussing our notation for purely-virtual, $L$-loop Feynman integrals and recalling some useful facts about them. For the direct integration of purely-virtual Feynman integrals, a very common starting point is the Feynman (or Schwinger) parametric representation (see e.g. Smirnov:2004ym for a detailed exposition). In many cases, it is convenient to write down the Feynman parametric representation in Euclidean space, treating all $n$ external momenta, $\{p_{i}\}$, on an equal footing by taking them all to be outgoing. In the most general case Smirnov:2010hn , it suffices to consider Feynman integrals of the form

[TABLE]

where $Q_{\ell}(k_{i},p_{j})$ denotes the momentum of the $\ell$-th propagator and the $N$ propagators in (1) are linearly independent. It is often the case that one can profitably work with the Feynman representation for the all-plus metric and ultimately obtain results which differ from the Minkowski space results in an appropriate Euclidean kinematic region only by trivial phases.666Of course, certain assumptions must be satisfied. For a more in-depth discussion, see reference Panzer:2015ida .

This approach has been used to great effect in recent years by Brown, Panzer, and others Brown:2008um ; Brown:2009ta ; Panzer:2013cha ; Panzer:2014gra ; Panzer:2014caa ; vonManteuffel:2014qoa ; Panzer:2015ida ; vonManteuffel:2015gxa , culminating recently in an impressive calculation of the six-loop $\beta$ function in $\phi^{4}$ theory Kompaniets:2016hct . Working through the details of the straightforward derivation (see e.g. Smirnov:2004ym ), one finds the all-plus Feynman parametrization

[TABLE]

where $\mathcal{U}_{\rm E}$ and $\mathcal{F}_{\rm E}$ are respectively the first and second Symanzik polynomials Bogner:2010kv in Euclidean space and $\nu=\sum_{i=1}^{N}\nu_{i}$. However, for the evaluation of Feynman integrals relevant to the computation of collider observables, it is arguably more natural to work in Minkowski space from the very beginning, considering Feynman integrals of the form

[TABLE]

with a momentum flow suitable for the description of a scattering experiment.

As is well-known, the derivation of Eq. (2) goes through with minor modifications if one works directly in Minkowski space. The mostly-minus Feynman parameter representation has the form

• [TABLE]

where $\mathcal{U}_{\rm M}$ and $\mathcal{F}_{\rm M}$ are the first and second Symanzik polynomials in Minkowski space. For our subsequent analysis of the Baikov-Lee representation, it is important to note that the functional dependence of $I_{\rm E}$ and $I_{\rm M}$ on the external kinematics is nearly identical.777At this juncture, it is of critical importance to clarify that, strictly speaking, this statement is not true. Obviously, a vector in Euclidean space which squares to zero is identically zero, whereas this is not the case in Minkowski space. However, one may simply write a formal expression for a Euclidean Feynman integral with the squares of certain momenta set to zero, remembering that, to be rigorous, one would have to work out the connection between Euclidean and Minkowski space representations with fake external masses and then set them to zero after the fact. Given some spanning set of external kinematic invariants, $\{\omega_{1},\ldots,\omega_{n(n-1)/2}\}$, constructed along the lines described in Bern:1992em , we can straightforwardly obtain one from the other,

[TABLE]

provided that we remember the $+i0$ prescription for the external kinematic invariants and define $\{p_{j}^{*}\}$ to be the set of external momenta which happen to be incoming in the physical kinematics of interest. To understand the above relation, recall that $\mathcal{U}_{\rm M}=\mathcal{U}_{\rm E}$ and that one can generate the Minkowski space function $\mathcal{F}_{\rm M}$ from $\mathcal{F}_{\rm E}$ by flipping the signs of all generalized Mandelstam variables and external masses which appear and, subsequently, appropriately adjusting the signs of the external momenta which must now be regarded as incoming.888Naturally, we have assumed that Euclidean and Minkowskian generalized Mandelstam invariants are defined in the usual way. For instance, one would have $t=(p_{1}+p_{3})^{2}$ in Euclidean space but $t=(p_{1}-p_{3})^{2}$ in the usual physical kinematics for $2\to 2$ scattering.

For what concerns the explicit examples discussed in the following sections, we essentially adopt the conventions of reference vonManteuffel:2014qoa . That is to say, for our actual calculations, we consider Minkowskian purely-virtual Feynman integrals of the form

[TABLE]

in physical kinematics. Note that we do not include factors in the measure to prevent the Euler-Mascheroni constant from appearing in $\epsilon$-expanded expressions because, in this paper, we either study cut Feynman integrals at $\mathcal{O}(\epsilon^{0})$ or to all orders in $\epsilon$. To simplify our discussion later on, it is also convenient to introduce complex-conjugated purely-virtual Feynman integrals, $I^{\dagger}$, where the $+i0$ prescription becomes a $-i0$ prescription and the $i\pi^{d/2}$ factors in (6) above are replaced by factors of $-i\pi^{d/2}$. The maximally-cut examples of Section 4 are far less sensitive to such details because overall phases make no difference at all if the only goal is to produce a valid solution to a given homogeneous differential equation. For the sake of definiteness, we will use the same normalization conventions in both Sections 3 and 4.

2.2 The Euclidean Baikov-Lee representation and its analytical continuation

To write the Euclidean Baikov-Lee formula succinctly, let us first recall that the Gram determinant on the $K$ linearly independent vectors $\{q_{i}\}$ is given by

[TABLE]

If we let $q_{i}$ be an element of the set $\{k_{1},\ldots,k_{L},P_{1},\ldots,P_{n-1}\}$,999Here, $\{P_{i}\}$ is nothing but a convenient permutation of the set of independent external momenta, $\{p_{i}\}$. This formulation is convenient because it is often desirable to eliminate a momentum other than the $n$-th. Our treatment of the $s$-cut of the massless one-loop box integral in Section 3.2 clearly illustrates this point. the Baikov-Lee representation of the purely-virtual Feynman integral $I_{\rm E}$ defined above in Eq. (1) is then

[TABLE]

where $\mathcal{D}$ is the domain of integration.

Even if one works in Euclidean space, finding an explicit description of the integration domain is in general a non-trivial task. To understand how this works in practice, let us consider the evaluation of the one-loop bubble with no internal masses in the Baikov-Lee approach. This example will both illustrate a general strategy for the determination of the integration region (briefly discussed in Frellesvig:2017aai ) and give the reader a sense as to why it is more convenient to integrate purely-virtual Feynman integrals using Feynman parameters as a starting point. The arguments advanced in this section are generally applicable, but it may be quite challenging to work out the details in examples with many integration variables and/or rich analytic structures. For future applications, we expect tools for the explicit solution of systems of inequalities such as the Reduce routine of Mathematica to play an important role.

For the one-loop bubble, a possible routing of the propagator momenta is

[TABLE]

In this case, (2.2) becomes

[TABLE]

where we have set $d=4-2\epsilon$. At this stage, it is important to note that the form of Eq. (2) and the definitions of $\mathcal{U}_{\rm E}$ and $\mathcal{F}_{\rm E}$ (see e.g. Bogner:2010kv ) together imply that Euclidean Feynman integrals are positive definite if all input kinematic variables are positive definite. This is a powerful analytic constraint from the Baikov-Lee point of view and effectively determines the shape of the integration region. The Baikov polynomial $p^{2}q_{1}^{2}-(q_{1}\cdot q_{2})^{2}$ inside the integrand above depends on the variables of integration and is raised to a non-integer power. This is a generic feature of Baikov-Lee calculations. The point is that, in order to keep the solution real-valued, one must consider e.g. $q_{1}\cdot q_{2}$ to lie between the branching points $\pm\sqrt{p^{2}q_{1}^{2}}$.

Apart from the obvious positivity of $q_{1}^{2}$, there are no further constraints on the variables of integration in this case, and we arrive at

[TABLE]

Although, the Baikov polynomial is no more than quadratic in the scalar products involving the loop momenta, the situation may become more complicated once the first scalar product is integrated out. In favorable cases, it is possible to find an analogous integration variable at each step of the calculation. However, it is not guaranteed that all polynomial structures remaining in the integrand after some number of integration steps have at most quadratic dependence on the remaining variables of integration. Although it is probably clear already, let us emphasize that the two-fold integral above is far, far more complicated than the trivial one-fold integral which one finds in the Feynman parametric approach to this simple problem. The evaluation of (10) involves non-trivial hypergeometric function identities, and we refer the interested reader to Appendix C for a detailed discussion.

As we shall see, the procedure described above for the limits of integration is conceptually even simpler for Feynman integrals cut in a single kinematic channel. In such cases, one can also make use of the fact that the region of integration is bounded; cut Feynman integrals of this type will be closely related to the real radiative master integrals for some physical decay or scattering process which has a finite amount of energy and momentum in the initial state Gehrmann-DeRidder:2003pne . The story is otherwise analogous to what was described above for purely-virtual Euclidean Feynman integrals because, up to phase, one again has a natural positivity condition. That is to say, in favorable cases, one can integrate each variable between branching points of the current integrand and then assign a positive orientation to the integration contour for the current variable (i.e. one must integrate $q_{1}\cdot q_{2}$ from $-\sqrt{p^{2}q_{1}^{2}}$ to $\sqrt{p^{2}q_{1}^{2}}$ in the above example, not vice versa). We have applied these ideas to explicitly evaluate a variety of Feynman integrals cut in a single kinematic channel at one, two, and three loops.

Of course, before defining the cut Baikov-Lee representation which we will study throughout the rest of this paper, we first need to analytically continue Eq. (2.2). The key idea is to recognize that (2) and (2.2) are nothing but two different integral representations of the same function, $I_{\rm E}$. Since our recipe to pass from $I_{\rm E}$ to $I_{\rm M}$, Eq. (5), does not depend at all on the details of the Feynman representation, it is natural to apply it to Eq. (2.2) as well,101010Due to the fact that the Baikov-Lee representation utilizes scalar product integration variables which have non-trivial dependence on the external momenta, it is not obvious that one can work in this way. However, we have found experimentally that this prescription does in fact make sense provided that one allows $\{P_{j}^{*}\}\to-\{P_{j}^{*}\}$ to act on the relevant scalar product integration variables as well. thereby obtaining a putative physical, Minkowski space version of the Baikov-Lee representation,

[TABLE]

where $\bar{G}(P_{1},\ldots,P_{n-1})$, $\bar{\mathcal{D}}$, $\bar{G}(q_{1},\ldots,q_{n+L-1})$, and $\bar{Q}_{\ell}(q_{i}\cdot q_{j})$ denote the various objects which appear in Eq. (2.2) after the replacements prescribed by (5) have been implemented.

2.3 The cut Baikov-Lee representation and unitarity

The key ingredient missing from the discussion so far is the generalized Cutkosky cutting rule written down by Lee and Smirnov in reference Lee:2012te . In a nutshell, they suggest that one can treat cut Feynman integrals with propagator denominators raised to powers greater than one by simply differentiating both sides of Cutkosky’s relation,

[TABLE]

an appropriate number of times with respect to $k^{2}$.111111We thank Gil Paz for pointing out that the distributional identity behind the generalized cutting rule was available in textbooks on the subject (e.g. GelfandShilov:1964 ) long before the appearance of reference Lee:2012te . To avoid digressing into a lengthy discussion of distributional calculus,121212The direct integration of delta distributions and their derivatives is straightforward Folland:1992 , but it seems somewhat less convenient from the perspective of implementation in a computer algebra system. it is convenient to actually define our Baikov-Lee representation for Feynman integrals cut in a single kinematic channel using the familiar language of residue calculus. The idea is that, up to a possible overall sign, the process of putting some number of propagators on the mass shell is completely equivalent to performing sequential residue computations which localize a subset of the scalar product integration variables. Our logic is similar to that of reference Abreu:2017ptx , except that, for our purposes, we find it more natural to repeatedly apply Cauchy’s residue theorem to Eq. (2.2). Due to the fact that we will use the main result of this section for the explicit examples discussed in the following sections, we find it natural to work with the absolute normalization of Eq. (6) in what follows.

For the sake of discussion, suppose that a particular unitarity cut of $I$ in the $\omega_{i}$ channel, say the $j$-th out of $M$, puts $n_{j}$ propagators on shell. By assumption, these propagators are linearly independent and there must therefore exist a subset of the scalar products depending on the loop momenta which one can sequentially integrate out using $n_{j}$ applications of the residue theorem. To simplify our notation, let $\{\bar{s}_{i}\}$ be the subset of the scalar product integration variables to be localized by the cut propagators, $\{s_{k}\}$ be the set of $L(L-1)/2+nL-n_{j}$ scalar product integration variables left over,131313In some cases, such as that of the one-loop double-cut bubble integral, the set $\{s_{k}\}$ is actually empty. and $\{\bar{Q}_{\ell}\}^{(j)}_{\omega_{i}{\rm-cut}}$ be the momenta of the cut propagators. Note that, at this stage, the order of both $\{\bar{s}_{i}\}$ and $\{\bar{Q}_{\ell}\}^{(j)}_{\omega_{i}{\rm-cut}}$ should be fixed to match the order in which the localization of the propagators and associated scalar products will be implemented by the residue theorem.

Finally, we define the Baikov-Lee representation of the $j$-th Feynman integral cut in the $\omega_{i}$ channel to be

[TABLE]

where the $\mathbf{sgn}$ function returns the sign of its argument and $\mathbf{Res}$ denotes the sequence of residue computations which localizes the $\{\bar{s}_{i}\}$. The presence of the $\mathbf{sgn}$ factor is necessary because sequential residue computations do differ from sequential localizations implemented with delta distributions and their derivatives in one important aspect. The issue is that

[TABLE]

but

[TABLE]

for arbitrary test functions $f(z)$ regular at $z=a$. In fact, one must include the sign of the Jacobian factor in (2.3) above, or the definition yields nonsensical results which may vary depending upon precisely what momentum routing is chosen for the cut Feynman integral under consideration.

We explicitly evaluate a number of cut Feynman integrals in Section 3 using Eq. (2.3) as a starting point, and we find that our definition is consistent in all cases. In particular, we find the expected unitarity relation between sums of cut Feynman integrals and the direct discontinuities141414We define the direct discontinuity of a Feynman integral in Appendix B. of their purely-virtual counterparts tHooft:1973wag ,

[TABLE]

In fact, we have successfully used our framework to study a large number of other examples of comparable complexity at one, two, and three loops. However, since we have no proof that our formulation is equivalent to the usual one where one considers all Feynman integrals to be embedded in an ambient generalized scalar field theory, it is important to write down a cut integral definition and associated unitarity relation along the lines of Gehrmann-DeRidder:2003pne ; Abreu:2014cla . If $\hat{I}=i^{L}I$, then

[TABLE]

and

[TABLE]

where $\mathcal{P}\left\{\hat{I}\right\}$ and $\mathcal{P}\left\{\hat{I}_{\omega_{i}{\rm-cut}}^{(j)}\right\}$ are scalar field theory phase factors defined in e.g. Abreu:2014cla .

For the purposes of this paper, we employ Eq. (2.3) with the attitude that it streamlines the exposition in Section 3 and makes it easier for the reader to check our analysis using the results from the literature collected in Appendix B (i.e. it allows us to forget about the annoying Feynman graph-dependent phase factors on both sides of (18)). Before leaving this section, let us emphasize that one can also employ the formalism discussed above to treat maximally-cut Feynman integrals in an analogous fashion. Interesting examples of maximally-cut Feynman integrals will be discussed in Section 4.

3 Discontinuities from cuts: one- and two-loop examples

To get a feeling for the ideas put forward in Section 2, we now consider a number of illustrative one- and two-loop examples. In the spirit of reference Abreu:2014cla , we compute the $s$-channel cuts of selected Feynman integrals using our cut Baikov-Lee representation, Eq. (2.3), and demonstrate that, in all cases, our results match the predictions of the optical theorem (i.e. the predictions obtained by using Appendix B to compute the direct discontinuities on the left-hand side of Eq. (16)).

3.1 The one-external-mass one-loop triangle

In this section, we consider the $s$-channel cut of the one-external-mass one-loop triangle with positive integer propagator exponents,

[TABLE]

This example clearly demonstrates the applicability of our formalism to propagators of higher multiplicity, and its elementary nature should give the reader ample opportunity to adjust to our notation. Clearly, there is only one cut Feynman integral which needs to be evaluated. We have

[TABLE]

in the cut Baikov-Lee representation of Eq. (2.3), where $s_{1}=k_{1}^{2}$, $\bar{s}_{1}=k_{1}\cdot p_{2}$, $\bar{s}_{2}=k_{1}\cdot p_{1}$, and $s=(p_{1}+p_{2})^{2}$.

Eq. (3.1) can be conveniently rewritten as

[TABLE]

to make manifest the fact that its right-hand side is purely imaginary for $s>0$ and to streamline the applications of the residue theorem which follow. Using the principle of mathematical induction, it is straightforward to evaluate the symbolic derivatives which enter into the residue calculations. Reading the list of barred integration variables from left to right, we find

[TABLE]

At this stage of the calculation, we have to analyze the integrand to determine the $s_{1}$ integration domain. By studying the form of Eq. (3.1), we see immediately that the integrand has precisely two branching points, at $s_{1}=0$ and at $s_{1}=s$. From the general discussion in Section 2.2, it therefore follows that

[TABLE]

For $\epsilon$ such that $\mathfrak{Re}(\nu_{1}+\nu_{2}+\epsilon)<3$ and $\mathfrak{Re}(\nu_{2}+\nu_{3}+\epsilon)<2$, the above integral converges and we find that

[TABLE]

Under an additional assumption, Eq. (3.1) can be simplified using the Saalschütz summation formula. For $\nu_{2}>1$, Eq. (86) implies that

[TABLE]

Actually, the principle of analytical continuation allows us to conclude that Eq. (25) is not only valid for arbitrary positive integer propagator exponents as desired, but that it even furnishes a definition of the cut Feynman integral for generic complex values of the propagator exponents. In fact, this analytical continuation of the above result will prove useful later on in Section 3.3. Finally, one can readily check using (25) and the $s>0$ evaluation of (19) given in Eq. (91) that

[TABLE]

3.2 The massless one-loop box

For our next example, we consider the $s$-channel cut of the massless one-loop box integral,

[TABLE]

Although it is again the case that just one cut Feynman integral needs to be evaluated, it is useful to study this example because it illustrates the applicability of our framework to multi-scale problems. In fact, it is not obvious to us that this cut can be computed to all orders in $\epsilon$ using traditional cut parameterizations. We shall see that it is also not entirely straightforward in the Baikov-Lee approach; the cut Feynman integral of interest here contains structures which bear a remarkable resemblance to those which appeared during the evaluation of certain three-loop, single-scale cut Feynman integrals Li:2014bfa .

From Eq. (2.3), we obtain

[TABLE]

for our cut Baikov-Lee representation after performing some trivial manipulations to massage the expressions into a more convenient form. In Eq. (28), we have made the definitions $s_{1}=k_{1}\cdot p_{1}$, $s_{2}=k_{1}\cdot p_{4}$, $\bar{s}_{1}=k_{1}^{2}$, $\bar{s}_{2}=k_{1}\cdot p_{2}$, $s=(p_{1}+p_{2})^{2}$, and $u=(p_{1}-p_{4})^{2}$. Carrying out the residue computations, we find

[TABLE]

Next, we integrate out the variable $s_{2}$. As explained in Section 2.2, the integration runs between the real zeros of the polynomial from the last line of Eq. (3.2),

[TABLE]

In fact, by observing that the $s_{2}$ integration is completely analogous to the integration with respect to $q_{1}\cdot q_{2}$ carried out in Appendix C, we can already anticipate from the branch cut structure of the $s_{2}^{\pm}$ (with respect to $s_{1}$) that the $s_{1}$ integration will run between [math] and $s/2$. We therefore readily obtain

[TABLE]

where we have made the change of variables $s_{1}=s/2\,x$ in the last line.

At this point, the above result may be rewritten to exhibit a ${}_{2}F_{1}$ of argument $\frac{4z}{(1+z)^{2}}$, where $z$ can be chosen to be either $\sqrt{\frac{1-x}{\frac{s+u}{-u}x}}$ or its reciprocal. From this observation, we see that a simple strategy to eliminate the square roots appearing in Eq. (3.2) is to split the integral at the point where $z=1$ and then apply a quadratic hypergeometric function transformation which is valid for $|z|<1$ (i.e. Eq. (84)) to both terms. Carrying out these steps, we find

[TABLE]

Now that the square root structures have been eliminated, we can deal with the two terms in (3.2) above by mapping them to linear combinations of known generalized Euler integrals. The first step is to make the change of variables $x=\frac{y}{\frac{s+u}{-u}+y}$ in the first integral and the change of variables $x=\frac{1}{1+\frac{s+u}{-u}y}$ in the second integral to bring them into the generalized Euler form. In fact, after making these transformations, one can immediately evaluate the second integral by applying integration formula (89):

[TABLE]

In our opinion, it is most convenient to deal with the remaining integral by replacing ${}_{2}F_{1}\left(1,1+\epsilon;1-\epsilon;y\right)$ with its integral representation, Eq. (80), and then applying Eq. (81) to integrate out the variable $y$ with an Appell series. This yields

• [TABLE]

In this case, the sum of the second and third parameters of the $F_{1}$ is equal to the fourth parameter and reduction formula (88) therefore immediately leads to

[TABLE]

Finally, we can map the remaining integral onto a linear combination of standard Euler integrals via connection formula (83),

[TABLE]

where we have obtained the second equality by evaluating the integral on the first line of Eq. (3.2) with the help of (80). Applying Eq. (82) to the final integral remaining on the right-hand side of (3.2), we see that the cut Feynman integral evaluates to

[TABLE]

Note that the final integral on the right-hand side of (3.2) exactly cancels the ${}_{3}F_{2}$ term, thereby removing all dependence on the generalized hypergeometric series from the result. One can check using (37) and the physical region ($s>0$ and $-s<u<0$) evaluation of (27) given in Eq. (92) that

[TABLE]

3.3 The one-external-mass six-line two-loop double triangle

Our final example will be the $s$-channel cut of the one-external-mass six-line two-loop double triangle,

[TABLE]

The $s$-channel cut of this integral was chosen because three cut Feynman integrals contribute to it and one of these integrals has an integration domain which is non-trivial to determine. This example will show the reader what Baikov-Lee computations look like beyond one loop, where irreducible scalar products and real-virtual contributions come into play for the first time. As before, we begin with Eq. (2.3). This time, however, we must enumerate the distinguishable cut Feynman integrals which contribute. Let us consider the triple cut first and the double cut second. The conjugate of the double-cut contribution can be obtained from the double cut without any additional calculation.151515By virtue of the $+i0$ in Eq. (3), the double cut contribution and its conjugate are distinguishable.

After performing some trivial manipulations, we obtain a cut Baikov-Lee representation of the form

[TABLE]

for the triple-cut contribution. In Eq. (40), we have made the definitions $s_{1}=k_{1}\cdot p_{1}$, $s_{2}=k_{1}^{2}$, $s_{3}=k_{2}\cdot p_{1}$, $s_{4}=k_{2}\cdot p_{2}$, $\bar{s}_{1}=k_{2}^{2}$, $\bar{s}_{2}=k_{1}\cdot k_{2}$, $\bar{s}_{3}=k_{1}\cdot p_{2}$, and $s=(p_{1}+p_{2})^{2}$. Carrying out the residue computations, we find

[TABLE]

We begin our treatment of the unbarred integration variables with $s_{4}$. Following the discussion of Section 2.2, we obtain our integration domain by studying the zeros of polynomial structures inside the integrand. Due to the fact that the polynomial

[TABLE]

is raised to the power $-1/2-\epsilon$, the integrand obtained above in (41) has the branch points

[TABLE]

in $s_{4}$. As in the previous example, we can actually deduce the limits of integration on the remaining variables by simply looking at the branch cut structure of $s_{4}^{\pm}$.161616It is worth pointing out that our logic here is similar to that used to formulate the compatibility graphs method for purely-virtual Feynman integrals Brown:2009ta . Namely, we do not expect new singularity structures to arise beyond those which are already encoded in the analytical structure of the initial integrand. First, the polynomial structures $s_{2}(2s_{1}+s_{2})$ and $s_{3}(s_{3}-s_{1})$ under the radical tell us that the limits of integration for $s_{2}$ are either $[0,-2s_{1}]$ or $[-2s_{1},0]$ and that the limits of integration for $s_{3}$ are either $[0,s_{1}]$ or $[s_{1},0]$. The only question that remains is whether $s_{1}$, what we shall consider to be the final variable of integration, is positive or negative.171717All possible limits of integration for $s_{2}$ and $s_{3}$ force one of the $s_{1}$ integration limits to be [math]. In the physical region, the polynomial $2s_{1}+s$ under the radical implies that the variable $s_{1}$ is a negative number which runs between $-s/2$ and [math]. Finally, we conclude that the $s_{2}$ integration runs between [math] and $-2s_{1}$ and that the $s_{3}$ integration runs between $s_{1}$ and [math]. The upshot is that

[TABLE]

In this case, all of the remaining integrations are elementary. The $s_{4}$ integration is closely analogous to the integration with respect to $q_{1}\cdot q_{2}$ carried out in Appendix C and the other integrations may be carried out with the help of a computer algebra system such as Mathematica. We find

[TABLE]

for the triple-cut contribution.

Our next task is to calculate the double-cut contribution. In fact, we can write down the answer immediately by recycling calculations that we have already carried out. We require only the result obtained in Section 3.1 for the $s$-channel cut of the one-external-mass one-loop triangle integral together with the well-known result for its purely-virtual counterpart. The key observation is that, on the support of the double cut, the virtual part of the cut one-external-mass six-line two-loop double triangle is precisely a purely-virtual one-loop triangle with external mass $k_{1}^{2}$. From Eq. (91), we see that181818Here, the $+i0$ prescription for the propagator plays a crucial role.

[TABLE]

It therefore follows that we can treat the integration over $k_{1}$ using Eq. (25) with propagator exponents $\nu_{1}=1$, $\nu_{2}=1$, and $\nu_{3}=2+\epsilon$. The desired result is

[TABLE]

The final cut Feynman integral of interest is the conjugate double-cut contribution. The result may be simply obtained by taking the complex conjugate of the virtual part of the double-cut contribution. We have

[TABLE]

Finally, the sum of the three contributions to the $s$-channel cut, Eqs. (3.3), (3.3), and (3.3), may be rewritten as

[TABLE]

Using (48) and the $s>0$ evaluation of (39) given in Eq. (93), the reader can check that

[TABLE]

4 Maximally-cut Feynman integrals and differential equations

As mentioned in the introduction, it was observed in Primo:2016ebd that maximally-cut Feynman integrals Bern:2004ky ; Britto:2004nc in the $\epsilon\to 0$ limit solve the homogeneous parts of the differential equations satisfied by their uncut counterparts. In non-polylogarithmic cases, it was furthermore shown explicitly that maximally-cut Feynman integrals in the $\epsilon\to 0$ limit may often be computed by direct integration. However, complete solution sets of their homogeneous differential equations were obtained indirectly by using a convenient mathematical property of the complete elliptic integral of the first kind. Subsequently, it was shown in Frellesvig:2017aai that Baikov’s method may be applied to maximally-cut Feynman integrals and offers conceptual advantages relative to the traditional approach to maximally-cut Feynman integrals. However only one solution of their higher-order homogeneous differential equations was provided.

We explain in this section how one can directly obtain the necessary, complete solution sets for a class of interesting non-polylogarithmic Feynman integrals. As we shall see, the key idea is to again employ our cut Baikov-Lee representation, Eq. (2.3), but to allow for more general integration domains than have so far been considered. For all of the classical unitarity cuts studied in Section 3, a unique solution was obtained via a generalized phase-space volume computation. However, the criteria that we employed to determine the integration domain must now be generalized further to allow for multiple solutions, due to the fact that the differential equations satisfied by non-polylogarithmic Feynman integrals are of order greater than one. Although the primary focus of this section will be on maximally-cut Feynman integrals which evaluate to generic complete elliptic integrals, the ideas discussed here apply to maximally-cut Feynman integrals which evaluate to complete hyperelliptic integrals as well Carlson:2003 .

Let us consider the following family of elliptic curves,

[TABLE]

where $t$ is a parameter and $\{a_{i}(t)\}$ is the set of branch points of $y(x,t)$, with $a_{i}(t)\neq a_{j}(t)$ for $i\neq j$. Now, for some Feynman integral, $I$, let us suppose that we are able to obtain a one-fold integral representation for its maximal cut, $\bar{I}$, of the form

[TABLE]

where $R(x,y(x,t),t)$ is a rational function of its arguments. If $R(x,y(x,t),t)$ has no poles in $x$, it is then the case that a complete set of solutions to the homogeneous part of the differential equations satisfied by $I$ may be obtained by considering

[TABLE]

for $i\neq j$. If poles are present, one must also consider closed contour integrals around each pole,

[TABLE]

to find all possible solutions. In (53), $\gamma_{k}$ denotes a closed contour which encircles the $k$-th pole of $R(x,y(x,t),t)$, but no other pole or branch point of the integrand. Note that, quite generically, the set of solutions obtained in this manner will actually be overcomplete.

By using properties of elliptic curves we can argue that (52) and (53) represent a complete set of solutions. Integrals $\bar{I}(a_{i}(t),a_{j}(t))$ for $i\neq j$ and $\bar{I}(\gamma_{k})$ are periods of the elliptic curve $y(x,t)$ Kontsevich:2001 and the given homogeneous differential equation for $I$ with respect to $t$ is nothing but the associated Picard-Fuchs equation.191919It has been clear for a long time that Picard-Fuchs equations play a very important role in the theory of Feynman integrals (see e.g. MullerStach:2011ru ; MullerStach:2012mp ). By construction, these differential equations are the same for every period, and a complete set of periods provides a complete set of solutions to the Picard-Fuchs equation (for a comprehensive review of the subject of periods see Kontsevich:2001 and the references therein).

The prescription described above represents a substantial generalization of that described in Section 2.2. When considering a cut Feynman integral associated with the unitarity cut in some physical kinematic channel, one naturally expects to obtain a result which is real-valued up to an overall phase; in this context, however, no such constraint applies and the integration domain is no longer uniquely determined. Indeed, Eq. (52) implies that, in the elliptic case, we can integrate the maximal cut of $I$ over six distinct domains and this obviously leads to some solutions which possess both real and imaginary parts. If the integrand has poles one must also include solutions of the form (53). This would be relevant, for example, when considering a complete elliptic integral of the third kind.

In the absence of poles, complete elliptic integrals admit a simple description as closed contour integrals which wrap the torus Carlson:2003 . As claimed above, it is clear from this point of view that our solutions cannot form a linearly independent set. To see this, recall that the fundamental group of the torus, $\mathbb{Z}\times\mathbb{Z}$, is isomorphic to the first homology group (Hurewicz’s theorem Hatcher:2002 ). This means that, in the absence of poles, the torus admits just two independent cycles for us to integrate along. We can therefore conclude that four of the six functions generated by applying the prescription given in (52) above are actually spurious and may be disposed of. In general, one must also check whether contour integrals of the form (53) around different poles of the integrand yield linearly dependent results.

The above discussion generalizes and systematizes the analysis of, e.g., Adams:2013kgc to generic complete elliptic integrals of the form (51). Moreover, it is possible to generalize it to curves of higher genus, i.e. when the square of (50) is a polynomial of degree greater than four, by considering the set of periods over the relevant higher-genus Riemann surface. The application of these techniques to curves of genus greater than one goes beyond the scope of the present paper but will likely play a role in future calculations. In the following, we consider illustrative examples taken from the virtual corrections to Higgs + jet with exact top mass dependence.

4.1 A Higgs + jet non-polylogarithmic three-point function

As a first example, we consider the maximal cut of the two-loop crossed form factor,

[TABLE]

where we have made the definitions,

[TABLE]

The evaluation of this Feynman integral is relevant to the calculation of the non-planar part of the two-loop virtual corrections to Higgs + jet with exact top mass dependence. Its maximal cut was considered in Primo:2016ebd , where it was evaluated with a traditional cut parametrization in an effort to obtain a solution to the homogeneous part of the associated system of differential equations. Using Eq. (2.3), we arrive at the following one-fold integral representation of the maximally-cut Feynman integral,

[TABLE]

where we have made the definitions $s_{1}=k_{1}\cdot p_{1}$ and $s=(p_{1}+p_{2})^{2}$, and we have chosen to work in the physical region above threshold where $s>p_{2}^{2}>0$ and $m^{2}>0$.

We now take the first step towards finding a complete set of solutions to the homogeneous differential equation for the uncut integral by computing the maximal cut in $d=4$. We have,

[TABLE]

where

[TABLE]

Following our general prescription, concrete results are obtained by integrating between branch points of the integrand. In other words, we obtain six possible solutions by picking distinct pairs of elements from the set of branch points of the integrand,

[TABLE]

where we have introduced the convenient shorthand

[TABLE]

in (59). Note that, in the kinematic region that we are working in, the elements of (59) are real-valued and ordered from smallest to largest.

As discussed above, the solutions to our homogeneous differential equations will be complete elliptic integrals and these may be thought of as periods of the torus. The torus admits two linearly independent cycles and we therefore expect to find just two linearly independent periods. As we shall see, it is convenient to take

[TABLE]

to be our independent basis elements.

That $f_{1}$ and $f_{2}$ are actually independent periods may be seen by writing them in a standard form. For $f_{1}$, this is easily achieved by making the appropriate analytical continuation of Eq. (58) and then changing variables according to Primo:2016ebd

[TABLE]

The result is

[TABLE]

in the physical kinematic region of interest.202020In Eqs. (63) and (64), $K(z)$ is the complete elliptic integral of the first kind,

$\displaystyle K(z)=\int_{0}^{1}\frac{{\rm d}t}{\sqrt{(1-t^{2})(1-z\,t^{2})}}\,.$

Similar considerations lead to

[TABLE]

again in the kinematic region of interest. The other possible solutions may be written as linear combinations of $f_{1}$ and $f_{2}$ and it is clear that some of them will have both real and imaginary parts. For example, we have from (59) and the explicit expressions for $f_{1}$ and $f_{2}$ that

[TABLE]

where $f_{1}$ is purely imaginary and $f_{2}$ is purely real in the region where $s>p_{2}^{2}>0$ and $m^{2}>0$.

Finally, let us point out that we have explicitly checked that the maximal cut calculated in this section solves the associated homogeneous differential equation. A more non-trivial Higgs + jet four-point function is discussed in the next section as a further application of our ideas. The following calculation demonstrates the utility of integrating out cut loops one at a time. As pointed out in Frellesvig:2017aai , working in this way often allows one to write down a more compact integrand than in the straightforward all-at-once approach to the construction of a cut Baikov-Lee representation adopted so far in this paper.

4.2 A Higgs + jet non-polylogarithmic four-point function

We consider the following two-loop Higgs + jet integral, denoted in Bonciani:2016qxi as $f^{A}_{66}$,

[TABLE]

with propagators,

[TABLE]

As usual, the kinematics is

[TABLE]

The maximal cut of this Feynman integral was also considered more recently in both references Primo:2016ebd and Frellesvig:2017aai . For the purposes of our analysis in this section, it is convenient to work in the kinematic region where $s>p_{4}^{2}>0$, $s>4m^{2}>0$, and $p_{4}^{2}-s>t$.

As we are considering a four-point function at two loops, there are nine scalar product integration variables, and it would therefore seem that we must consider a three-fold integral representation of the maximal cut. Fortunately, one can obtain a one-fold integral representation of the maximal cut by proceeding recursively loop-by-loop. First, we integrate out the one-loop triangle subintegral defined by the propagators $D_{4}$, $D_{5}$, and $D_{6}$,

[TABLE]

by localizing $k_{2}^{2}$, $k_{2}\cdot p_{3}$, and $k_{2}\cdot p_{4}$. To do so, we evaluate the maximal cut of (69) using Eq. (2.3), our cut Baikov-Lee representation. Of course, the maximal cut of a one-loop triangle involves no non-trivial integrations and one immediately finds

[TABLE]

after carrying out the residue computations.

We now integrate out the remaining loop by localizing $k_{1}^{2}$, $k_{1}\cdot p_{1}$, and $k_{1}\cdot p_{4}$ with the remaining cut conditions. A moment’s thought reveals that Eq. (2.3) may still be straightforwardly applied to loop-by-loop Baikov-Lee calculations; the only difference is that the results of the previous loop integrations appear in the current integrand. In other words, if we make the definitions $s_{1}=k_{1}\cdot p_{3}$, $\bar{s}_{1}=k_{1}^{2}$, $\bar{s}_{2}=k_{1}\cdot p_{1}$, and $\bar{s}_{3}=k_{1}\cdot p_{4}$, we have

[TABLE]

and can immediately write

[TABLE]

where

[TABLE]

in the region of interest.

We now turn to the problem of finding a complete set of homogeneous solutions to the differential equations for the uncut integral, proceeding as described at the beginning of Section 4. Possible solutions are obtained by integrating between the branch points of the integrand,

[TABLE]

where we have set

[TABLE]

In the kinematic region that we have chosen to work in, the elements of (74) are real-valued and ordered from smallest to largest.

As guaranteed by the form of Eqs. (72) and (73), we again find just two linearly independent solutions,

[TABLE]

and

[TABLE]

That these two solutions are actually linearly independent follows from the explicit formulas written in terms of complete elliptic integrals of the first kind,

[TABLE]

and

[TABLE]

which may be derived by making appropriate changes of variables in Eqs. (76) and (77).212121To derive (78), one must analytically continue the second square root structure in Eq. (73) above. We have explicitly checked that $g_{1}$ and $g_{2}$ satisfy the appropriate homogeneous second-order differential equations.

5 Conclusions

In this paper, we formulated and studied cut Baikov-Lee representations, both for Feynman integrals cut in a single kinematic channel and for maximally-cut Feynman integrals. For a wide class of interesting problems, our framework provides a convenient setup for the explicit computation of cut Feynman integrals. It makes the dependence on the Lorentz-invariant kinematic variables manifest and may be used directly or in conjunction with other methods such as sector decomposition Binoth:2000ps ; Bogner:2007cr . Although some elements of our analysis in Section 2 relied upon physically-motivated plausibility arguments and experimentation, we subsequently presented a substantial amount of evidence in Sections 3 and 4 that our master formula, Eq. (2.3), is correct.222222Although we have employed Eq. (2.3) throughout this paper, Eq. (2.3) is actually on more solid ground from the theoretical point of view. The difference is that we have included all of the phases in Eq. (2.3) which one would find by considering the relevant Feynman integrals to be Feynman graphs inside of an appropriate generalized scalar field theory. It is not clear to us why Eq. (2.3) works as well as it does.

It would be very interesting in future work to consider still more non-trivial examples232323When a reference evaluation is not available, it is important to check analytical results numerically. We therefore note that the recently-released program pySecDec Borowka:2017idc should allow for the evaluation of a wide class of cut Feynman integrals numerically (up to some fixed order in $\epsilon$) with moderate user input. such as the $s$-channel cut of the massless two-loop non-planar double box; as discussed in vonManteuffel:2014qoa , examples for which one cannot avoid imaginary parts on the mass shell require special care and may be instructive. In an effort to remove as many superfluous assumptions as possible from the formulation given in Section 2, it would of course also be desirable to put the theoretical foundations of the cut Baikov-Lee representation on a firmer footing. Although we have employed the familiar language of classical complex analysis throughout this work, it might be interesting to reformulate our findings in more modern language along the lines of Abreu:2017ptx ; Abreu:2017enx . It is unclear to us, however, that such a reformulation will immediately lead to clarifications.

In fact, there exist several interesting classes of cut Feynman integrals which were not discussed in this work at all. First of all, it would be interesting to study representative sequential cuts of the type discussed in reference Abreu:2014cla . One should also check whether a cut-discontinuity relation of the type discussed in Abreu:2014cla also exists for “crossed” sequential cuts such as the $s$-channel $+$ $t$-channel cut of the massless one-loop box of Section 3.2. Although success is less certain, it might be interesting to use the Baikov-Lee formalism developed in this work to study iterated cuts in a single channel. As a start, one could consider the Feynman integral analog of the double two-particle cuts at two loops discussed in Bern:2000dn . For such iterated cuts, it is not obvious that a cut-discontinuity relation exists at all, and it would therefore be interesting to take a fresh look at the problem using our Baikov-Lee machinery. Finally, it goes almost without saying that we would very much like to apply our techniques to the evaluation of the master integrals relevant to the current generation of phenomenologically-important unsolved problems in perturbative quantum field theory.

Acknowledgments

The authors would especially like to thank Ruth Britto for many interesting discussions, support, and comments on the manuscript. We also gratefully acknowledge an illuminating discussion with Stefan Müller-Stach and thank him for reading excerpts from our manuscript. RMS would like to thank Sven-Olaf Moch for an interesting discussion which, at least in part, inspired the author to pursue this line of research. This project has received funding from the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation programme under grant agreement No 647356 (CutLoops). Our figures were generated using Jaxodraw Binosi:2003yf , based on AxoDraw Vermaseren:1994je .

Appendix A Mathematical relations for hypergeometric-like functions

In this appendix, we review some well-known facts used in this paper about hypergeometric functions and their generalizations. First of all, let us recall the usual integral representation of the hypergeometric function which provides the analytical continuation of the hypergeometric series in non-exceptional cases. For $|{\rm arg}(1-z)|<\pi$ and $\mathfrak{Re}(c)>\mathfrak{Re}(a)>0$, we have Lebedev:1965

[TABLE]

In fact, for $|{\rm arg}(1-z_{1})|<\pi$, $|{\rm arg}(1-z_{2})|<\pi$, and $\mathfrak{Re}(c)>\mathfrak{Re}(a)>0$, a completely analogous formula holds for the Appell $F_{1}$ function Schlosser:2013 :

[TABLE]

We also encounter the generalized hypergeometric function ${}_{3}F_{2}$, which, for the purposes of this paper, may be defined via the integral representation

[TABLE]

which is valid for $|{\rm arg}(1-z)|<\pi$ and $\mathfrak{Re}(b_{2})>\mathfrak{Re}(a_{3})>0$ Lebedev:1965 .

In Lebedev Lebedev:1965 , one also finds a very clear discussion of both linear and quadratic transformations of the hypergeometric function ${}_{2}F_{1}$. Of particular interest to us are,

[TABLE]

[TABLE]

and

[TABLE]

Eq. (83) is valid for non-integral $a-b$, $|{\rm arg}(1-z)|<\pi$, and $|{\rm arg}(-z)|<\pi$, whereas Eq. (84) is valid for $2b\neq-1,-3,-5,\dots$ and, crucially, $\left|z\right|<1$. Eq. (85) is valid so long as the condition $|{\rm arg}(1-z)|<\pi$ is satisfied.

We also require some reduction identities, two for the generalized hypergeometric function ${}_{3}F_{2}$ and one for the Appell $F_{1}$ function. The Saalschütz summation formula Slater:1966 ,

[TABLE]

applies if $b_{1}+b_{2}-a_{1}-a_{2}-a_{3}=1$ and one element of $\{a_{1},a_{2},a_{3}\}$ is a negative integer. A more non-trivial summation formula involving two ${}_{3}F_{2}$ functions on the left-hand side is WolframFunctions

[TABLE]

Eq. (87) is valid for $\mathfrak{Re}(1+b_{1}-a_{1}-a_{2})>0$ and may be verified using the techniques described in Koornwinder:1998 . Finally, the Appell $F_{1}$ function collapses to a ${}_{2}F_{1}$ if $c=b_{1}+b_{2}$ Schlosser:2013 ,

[TABLE]

The standard integral representation given above for the generalized hypergeometric function ${}_{3}F_{2}$ is only one of a number of Euler integrals involving ${}_{2}F_{1}$ which may be evaluated using the ${}_{3}F_{2}$ series. Many evaluations of such generalized Euler integrals are given in reference Letessier:1988 . The main result of interest to us is

[TABLE]

which is valid for $|{\rm arg}(1-z)|<\pi$, $\mathfrak{Re}(\gamma)>0$, $\mathfrak{Re}(\rho)>0$, and $\mathfrak{Re}(\gamma+\rho-\alpha-\beta)>0$ Gradshteyn:2007 .

Appendix B All-order-in-$\epsilon$ results for selected purely-virtual Feynman integrals

In this appendix, we collect some useful results from the Feynman integral literature. All results which follow are presented in the normalization of Eq. (6) and are valid to all orders in the parameter of dimensional regularization, $\epsilon$, for generic phase-space points in physical kinematics. In order to compare with the cut calculations performed in Section 3, we must explain how to parse the ${\rm Disc}$ operation introduced in Section 2. The direct discontinuity of Feynman integral $I$ in the $s$-channel is simply

[TABLE]

where $\{v_{j}\}$ denotes the set of variables (the parameter of dimensional regularization, generalized Mandelstam variables, and, in general, internal masses) that $I$ is a function of.

B.1 The one-external-mass one-loop triangle

The one-external-mass one-loop triangle with generic propagator exponents is given by Smirnov:2004ym

[TABLE]

to all orders in $\epsilon$, where $s=(p_{1}+p_{2})^{2}>0$ and $\nu=\sum_{i=1}^{3}\nu_{i}$.

B.2 The massless one-loop box

The massless one-loop box is given by vanNeerven:1985xr ; Smirnov:2004ym

[TABLE]

to all orders in $\epsilon$, where $s=(p_{1}+p_{2})^{2}$ and $u=(p_{1}-p_{4})^{2}$. In the physical region, $s>0$ and $-s<u<0$.

B.3 The one-external-mass six-line two-loop double triangle

The purely-virtual counterpart of the one-external-mass six-line two-loop double triangle integral studied in Section 3 is given by

[TABLE]

to all orders in $\epsilon$, where $s=(p_{1}+p_{2})^{2}>0$. To our knowledge, Eq. (93) was first derived by van Neerven vanNeerven:1985xr . His idea was to first calculate all $s$-channel cuts and then deduce the associated purely-virtual result using unitarity. Consequently, to obtain (93) without referring to cuts, it was necessary for us to evaluate the integral ourselves using Feynman parameters. This exercise is elementary and may be carried out using the loop-by-loop integration strategy suggested in Gonsalves:1983nq .

Appendix C Baikov-Lee for the purely-virtual one-external-mass one-loop bubble

In this appendix, we complete the calculation of the Euclidean purely-virtual one-loop bubble integral with no internal masses which was initiated in Section 2. Our point of departure will be Eq. (10),

[TABLE]

The first step is to map the domain of the first integration variable, $q_{1}\cdot q_{2}$, onto the unit interval. This can be achieved straightforwardly by making the change of variables $q_{1}\cdot q_{2}=2\sqrt{p^{2}q_{1}^{2}}z-\sqrt{p^{2}q_{1}^{2}}$.242424A variable change of this form often allows one to recognize the definite integrals which arise from all-orders-in-$\epsilon$ Feynman integral calculations as Euler integrals of hypergeometric type (see e.g. Gehrmann-DeRidder:2003pne ). We arrive at

[TABLE]

By comparing Eq. (95) to the form of Eq. (80), it is now obvious that the $z$ integral may be evaluated in terms of the ${}_{2}F_{1}$ series.

The one-fold integral that remains,

[TABLE]

is most naturally evaluated by splitting the integral at the point $q_{1}^{2}=p^{2}$ and then mapping both the integral from [math] to $p^{2}$ and the integral from $p^{2}$ to $\infty$ onto the unit interval. This will allow for the simultaneous application of quadratic transformation (85) to both integrals. Making the change of variables $q_{1}^{2}=p^{2}x$ in the first integral and the change of variables $q_{1}^{2}=p^{2}/x$ in the second integral, we find

[TABLE]

At this stage, we can straightforwardly evaluate both integrals using Eq. (82). The result obtained in this manner,

[TABLE]

is correct but far more complicated than it needs to be. In this case, we can simplify the result by applying summation formula (87) to eliminate the second ${}_{3}F_{2}$ series in Eq. (98) above:

[TABLE]

Needless to say, the above result agrees with what one obtains (far more easily) using the Feynman representation.252525One can see that this is the case by consulting a standard text such as Smirnov Smirnov:2004ym .