arxiv: v2 [hep-ph] 12 Jul 2017

Transcription

1 MaPhy-AvH/7-5 MITP/7-38 arxiv:75.895v [hep-ph] Jul 7 Analytic continuation and numerical evaluation of the kite integral and the equal mass sunrise integral Christian Bogner a, Armin Schweitzer a and Stefan Weinzierl b a Institut für Physik, Humboldt-Universität zu Berlin, D - 99 Berlin, Germany b PRISMA Cluster of Excellence, Institut für Physik, Johannes Gutenberg-Universität Mainz, D Mainz, Germany Abstract We study the analytic continuation of Feynman integrals from the kite family, expressed in terms of elliptic generalisations of multiple) polylogarithms. Expressed in this way, the Feynman integrals are functions of two periods of an elliptic curve. We show that all what is required is just the analytic continuation of these two periods. We present an explicit formula for the two periods for all values of t R. Furthermore, the nome q of the elliptic curve satisfies over the complete range in t the inequality q, where q = is attained only at the singular points t {m,9m, }. This ensures the convergence of the q-series expansion of the ELi-functions and provides a fast and efficient evaluation of these Feynman integrals.

2 Introduction Precision calculations in high-energy particle physics require the evaluation of Feynman loop integrals. Within an analytical approach towards Feynman loop integrals one computes first the Feynman integrals in terms of a specific class of transcendental function in a particular kinematic region usually the Euclidean region). In a second step one studies the analytic continuation into the full kinematic region and the numerical evaluation of the transcendental functions. This procedure has been successful for a wide class of Feynman integrals, which evaluate to multiple polylogarithms [ 4]. Numerical algorithms are available to evaluate multiple polylogarithms for arbitrary complex-valued arguments [5 ]. It is well-known that not all Feynman integrals can be expressed in terms of multiple polylogarithms. The simplest counter-example is given by the two-loop equal mass sunrise integral. By now, this integral has been studied extensively in the literature [ 3]. This is of course justified by the fact that the sunrise integral is the first Feynman integral which cannot be expressed in terms of multiple polylogarithms. Further examples of Feynman integrals not expressible in terms of multiple polylogarithms are discussed in [3 37]. In this paper we are interested in a larger family of Feynman integrals, the family of kite integrals [3 33]. By the term family of Feynman integrals we mean Feynman integrals which differ only in the powers of their propagators, including the case, where one or more exponents are zero. In this case the corresponding propagator is absent and the Feynman integral reduces to a sub-topology. The kite integral contains as a sub-topology the equal-mass sunrise integral. Like the sunrise integral, the kite integral cannot be expressed in terms of multiple polylogarithms. Both the kite integral and the sunrise integral are two-point functions, depending on a variable t = p equal to the incoming or outgoing) momentum squared. In two recent publications we showed that the equal-mass sunrise integral [7] and the kite integral [33] can be expressed in a neighbourhood of t = to all orders in the dimensional regularisation parameter ε in terms of elliptic generalisations of multiple) polylogarithms. These generalisations are denoted as ELi-functions and are functions of several variables. Of particular interest is the dependence on one particular variable denoted by q. Feynman integrals are related to periods of algebraic varieties [38, 39]. In the case of the family of kite integrals the non-trivial algebraic variety which prohibits an evaluation in terms of multiple polylogarithms is an elliptic curve. The variable q is the nome of the elliptic curve. By the modularity theorem, every elliptic curve over Q has a modular parametrisation and in [9] one of the authors showed that one may express the family of kite integrals to all orders in the dimensional regularisation parameter ε as iterated integrals of modular forms. The ELi-functions give the q-series expansion of these iterated integrals of modular forms. At t = we have q = and therefore q < in a neighbourhood of t =. In this neighbourhood the ELi-functions provide a simple and convenient evaluation of the Feynman integrals of the kite family. In this paper we consider the analytic continuation and the numerical evaluation of the integrals of the kite family in the complete kinematic region t R. We note that in a recent paper [4] the analytic continuation of the ELi-functions has been discussed. It is one of the results of the present paper, that for the family of kite integrals the analytic continuation of the ELi-functions is not needed. The nome q is defined in terms of the two periods ψ and ψ of the elliptic curve by q=expiπψ /ψ ). All what is needed is the analytic continuation of the two periods ψ and

3 ψ. We present an explicit formula for the two periods for all values of t R. It turns out that once the two periods ψ and ψ are properly defined, we have q for all t R. and q = is attained only at the singular points t {m,9m, }. As a second main result of this paper we show that the same expressions in terms of the ELi-functions which we found for the Feynman integrals of the kite family in the region around t = hold for all t R. The ELi-functions provide therefore a fast and efficient way to evaluate these integrals over the complete kinematic range. Of course, Feynman integrals may also be evaluated by purely numerical methods [4 5]. A flexible numerical method to check Feynman integrals at individual kinematic points is given by sector decomposition [5 54]. We compare numerical results from the evaluation of the Feynman integrals from the kite family expressed in terms of ELi-functions with numerical results from the SeCDec program [55 58]. We find perfect agreement in all kinematic regions, including the regions close to the thresholds. As already mentioned, we only need to continue analytically the two periods ψ and ψ. Taking Feynman s i-prescription into account, the two periods are continuous functions of t. We express them in terms of complete elliptic integrals of the first kind. The complete elliptic integral of the first kind Kz) has a branch cut in the variable z at [, [. If we cross this branch cut, we have to compensate the discontinuity of Kz) by taking into account the monodromy around z =. This paper is organised as follows: In section we briefly recall the conventions as arguments of functions approach a branch cut for mathematical software on the one hand and Feynman s i- prescription in physics on the other hand. In section 3 we study a family of elliptic curves E t and their periods. We express the periods in terms of complete elliptic integrals of the first kind, such that the periods are continuous functions as t varies continuously. Equipped with the appropriate definition of the periods we discuss the analytic continuation of the Feynman integrals of the kite family in section 4. The correct definition of the periods requires to take a monodromy matrix at t = m into account. Section 5 is devoted to the derivation of this monodromy matrix. In section 6 we show numerical results for three Feynman integrals from the kite family. Finally, our conclusions are given in section 7. Appendix A gives information on our conventions regarding the roots of the cubic polynomial of the Weierstrass normal form. In appendix B we summarise the definition of the ELi-functions and the definition of the E-functions, the latter being linear combinations of the former. Appendix C reviews the algorithm for the numerical computation of the complete elliptic integrals of the first kind based on the arithmetic-geometric mean. Conventions We will encounter mathematical functions of a complex variable z, like z or the complete elliptic integral Kz). These functions have branch cuts. Let us clarify, which values to assign to these functions on the branch cuts. The standard convention for mathematical software is as follows: Implementations shall map a cut so the function is continuous as the cut is approached coming around the finite endpoint of the cut in a counter clockwise direction [59]. With this convention, cuts on the positive real axis are continuous to the lower complex half-plane while 3

4 cuts on the negative real axis are continuous to the upper complex half-plane. However, this is not always what we want in physics. In physics the analytic continuation is dictated by Feynman s i-prescription, where we substitute a real) variable t by t t+ i. The symbol +i denotes an infinitesimal positive imaginary part. The small imaginary part overrides the mathematical convention above. To give an example we have t = { i t, t >, t, t, t i = { i t, t >, t, t. ) In the sequel we will always assume that a small positive imaginary part is added to the variable t. Furthermore we will deal with equations of the form y = x e )x e )x e 3 )x e 4 ), ) and we would like to express y as the square root of the right-hand side. The square root of the right-hand side may be viewed as a multi-valued function of x, taking two possible value which differ by a sign. We are interested in a single-valued and continuous) function of x. Of course this cannot be done on the entire complex plane C, but only on the complex plane minus some cuts. Let us first assume that the roots e, e, e 3 and e 4 are real and ordered as e < e < e 3 < e 4. 3) Let us choose the cuts to be given by the line segment from e to e and by the line segment from e 3 to e 4. Then we may choose y as a single-valued and continuous function onc\[e,e ] [e 3,e 4 ]). It is not too difficult to show that one possible choice is given by the other choice is given by These functions have branch cuts at y = x e x e x e3 x e4, 4) y = x e x e x e3 x e4. 5) e i + λ, λ R. 6) On the interval [e,e 3 ] we have a superposition of the branch cuts starting at e 3 and e 4, making the function continuous there. The same happens on the interval ],e ], where we have a superposition of all four branch cuts. It is worth noting that the naive guess y = x e )x e )x e 3 )x e 4 ) 7) does in general not give a single-valued continuous function of x on C\[e,e ] [e 3,e 4 ]). The cuts of eq. 7) are in general algebraic functions of x and determined by x e )x e )x e 3 )x e 4 ) λ =, λ R. 8) 4

5 Finally, let us remove the assumption e i R and let us consider the general case e i C. We continue to denote by [e i,e j ] the line segment from e i to e j, now in the complex plane. Then we may express y as a single-valued and continuous function onc\[e,e ] [e 3,e 4 ]) through x e x e x e3 x e4 y = ±e e )e 4 e 3 ). 9) e e e e e 4 e 3 e 4 e 3 3 The elliptic curve We consider a family of elliptic curves E t P C) given in the chart z = by the Weierstrass normal form with g t) = g 3 t) = y = 4x 3 g t)x g 3 t), ) 3m µ 8 t ) 3m 6 3m 4 t+ 9m t t 3), 3m 4 6µ + 6m t t ) 9m 8 36m 6 t+ 3m 4 t m t 3 +t 4). ) As discussed previously in the literature, this elliptic curve is associated with the equal-mass sunrise integral [ 9]. Eq. ) is obtained by transforming the equation where F =, ) F = x x x 3 t+ m x + x + x 3 )x x + x x 3 + x 3 x ). 3) denotes the second graph polynomial of the equal-mass sunrise integral to the Weierstrass normal form. Graph polynomials are reviewed in [6]. We may factor the cubic polynomial on the righthand side of eq. ). We then obtain y = 4x e )x e )x e 3 ), with e + e + e 3 =, 4) where the roots are given by e = t 4µ 4 + 6m t+ 3m m t ) 3 9m t ) ), e = t 4µ 4 + 6m t+ 3m 4 3 m t ) 3 9m t ) ), e 3 = t 4µ 4 m t 6m 4). 5) In appendix A we discuss the motivation for defining the roots as in eq. 5). The roots e, e and e 3 vary with t. This is shown in fig.. At the singular points t {,m,9m, } two of the three roots coincide. In detail we have 5

6 roots, real part 4 Reei) 4 e e e t/m roots, imaginary part 4 Imei) 4 e 3 e e t/m Figure : The variation of the roots e, e and e 3 of the elliptic curve y = 4x e )x e )x e 3 ) with t. The roots e and e acquire an imaginary part for m < t < 9m. At t = the roots e and e 3 coincide. At t = m and at t = 9m the roots e and e coincide. At t = the roots e and e 3 coincide. In the interval m t 9m the real parts of the roots e and e coincide. 6

7 t m 9m roots e = e 3 e = e e = e e = e 3 The modulus k and the complementary modulus k of the elliptic curve are defined by k = e 3 e e e, k = k = e e 3 e e. 6) The modulus k and the complementary modulus k appear as arguments of the complete elliptic integral of the first kind Kk), defined by Kk) = dt t ) k t ). 7) It is clear from the definition that Kk) is only a function of k and some authors prefer therefore the notation Kk ). In this paper we follow the standard conventions [6] and define Kk) as in eq. 7). The complete elliptic integral of the first kind Kk) has branch cuts on the real axis at ], ] and[, [ in the complex k-plane. The function Kk )=Kk) has then a branch cut at [, [ in the complex k -plane. Let us denote by δ and δ two cycles on the elliptic curve E t which generate the homology group H E t,z). The differential forms η = dx y, η = xdx y 8) are generators of the cohomology group H dr E t). The periods of the elliptic curve are P i j = δ i η j, i, j {,}. 9) We are in particular interested in the periods ψ = P and ψ = P, involving η. These periods satisfy the differential equation [ d dt + t + t m + ) d t 9m dt + m )] 3t + 4t m ) + t 9m ψ i =. ) ) The periods P and P satisfy a slightly different differential equation.) Eq. ) is called the Picard-Fuchs equation. We define the two cycles δ and δ such that the periods ψ and ψ are given for t < by ψ = e 3 e dx y, ψ = 7 e 3 e dx y. )

8 Imt) I C II C III C 9 IV m 9m Ret) Figure : The path for the analytic continuation in the variable t. Feynman s i-prescription avoids the singular points at, m and 9m as shown in the figure. We take y = x e x e x e3. ) The path of integration in x-space is such that we have an infinitesimal small negative imaginary part for x. In other words, the integration path is below all cuts. Working out the integrals one finds in the region t < 4µ 4iµ ψ = Kk), ψ m t) 3 4 9m t) = K k ). 3) 4 m t) 3 4 9m t) 4 The periods ψ and ψ vary continuously with t. We study the variation of these periods along the path shown in fig. ). We divide the path into seven pieces: Four line segments and three small semi-circles. The line segments are characterised by Region I : < Ret) <, Region II : < Ret) < m, Region III : m < Ret) < 9m, Region IV : 9m < Ret) <. 4) The three semi-circles C, C and C 9 encircle the points, m and 9m, respectively. Everywhere on the path we have Imt)>. This implements Feynman s i-prescription. Of course, the path is equivalent to a straight line parallel to the real axis with a small imaginary part. However, it is advantageous to discuss the individual pieces separately, in particular the small semi-circles. We may extend the path to a closed) path on the Riemann sphere by adding a semi-circle at infinity. As already mentioned, the periods ψ and ψ are continuous functions of t. In eq. 3) we expressed the two periods for t in region I in terms of complete elliptic integrals of the first kind. The complete elliptic integral Kk) has branch cuts on the real axis. In terms of the variable k the branch cut is on the positive real axis given by the interval [, [. In fig. 3) we sketch the path in the complex k -space, as t varies along the path shown in fig. ). For t = +i we start in k -space at k = i, continuing below the real axis until we reach for t = m + i the point k = i. The small semi-circle C around t = is harmless and can be deformed 8

9 Imk ) III II C I IV Rek ) C 9 C Figure 3: The path in the complex k -space, as t varies along the path of fig. ). In the regions IV, I and II the path in k -space is at an infinitesimal distance below the real axis. The path crosses the branch cut [, ] in the region C. In region III we have Rek )=/. 9

10 away. However, the semi-circle C around t = m is essential. It corresponds to a three-quarter circle in k -space and crosses the branch cut [, [. The path continues in region III, where we have Rek ) = /. The semi-circle around t = 9m brings the path close to the real axis with Rek )>. In region IV the path continues at an infinitesimal distance below the real axis back to k = i. The three-quarter circle in fig. 3) corresponding to C is equivalent to the missing quarter-circle in the clockwise direction and a full circle in the anti-clockwise direction. In this way the monodromy of Kk) enters the expression for ψ. For ψ we study the path in k -space, as t varies along the path of fig. ). Due to the relation k + k =, 5) the path in k -space is simply obtained from the path in k -space by reflection on the point /. The path in k -space is sketched in fig. 4). Let us stress that the path in k -space does not cross the branch cut[, ]. In expressing the periods ψ and ψ for all values of t along the path of fig. ) in terms of complete elliptic integrals of the first kind Kk) and Kk ) we have to take the monodromy at t = m into account. The monodromy relation is derived in section 5. The result is as follows: For all values of t R the periods ψ and ψ are expressed as ψ t+ i) ψ t+ i) ) = 4µ m t i) 3 4 9m t i) 4 γ t ikk t+ i)) Kkt+ i)) ), 6) where the matrix γ t is given by γ t = ), < t < m, ), m < t <, 7) and encodes the monodromy. Note that in region II the elliptic integral Kk ) is evaluated above the branch cut, while in region IV the elliptic integral Kk) is evaluated below the branch cut. In both cases we may use for k ], [ the formula [6] [ Kk± i) = ) )] K ± ik k. 8) k k We note that this relation can also be used for an analytical continuation by fixing the constants in a small neighbourhood around the singular points [3, 35]. However, we will pursue a different and simpler) path here. Having defined the periods ψ and ψ in eq. 6), we introduce the standard definitions of the ratio of the two periods τ and the nome q of the elliptic curve: τ = ψ ψ, q = e iπτ. 9)

11 Imk ) C C 9 IV I C II Rek ) III Figure 4: The path in the complex k -space, as t varies along the path of fig. ). In the regions IV, I and II the path in k -space is at an infinitesimal distance above the real axis. In region III we have Rek )=/.

12 path in τ.4. Imτ) Reτ).5 path in q t/m ],] t/m [ : ] t/m [ : 9] t/m [9, [.5 Imq).5.5 Req).5 Figure 5: The path in τ-space and in q-space, as t varies along the path of fig. ). The value t = corresponds to τ = and q =, the value t = corresponds to τ = i and q =, the value t = m corresponds to τ= and q=, the value t = 9m corresponds to τ= 3 and q= i 3.

13 In fig. 5) we show the path in τ-space and in q-space as t varies along the path of fig. ). We have τt =± )=, τt = )=i, τt = m )=, τt = 9m )= 3. 3) For all values t along the path of fig. ) we have We have Imτ)>, q <. 3) Imτ)=, q =. 3) for t {m,9m,± }, i.e. without Feynman s i-prescription. There is a bijection between all points t on the path of fig. ) and the points τ on the path in τ-space of fig. 5). In the direction t τ the mapping is given by eq. 3) τ = ψ t) ψ t). 33) In the reverse direction we have t = 9m ητ)4 η ) 3τ 4 η6τ) 4 η ) τ 4 ητ) 4 η3τ) 4, 34) where ητ) denotes Dedekind s eta function, defined by ητ) = e iπτ n= e πinτ ). 35) We have checked numerically that eq. 34) is the inverse mapping to eq. 33) for all t on the path of fig. ). 4 Analytic continuation of Feynman integrals In this section we discuss the analytic continuation of the Feynman integrals associated to the kite family, i.e. the kite integral and all sub-topologies, including in particular the equal mass sunrise integral. We will see that this is trivial, once the periods ψ and ψ are defined by eq. 6). The Feynman integrals of the kite family depend on the variable t = p. 36) Feynman s i-prescription instructs us to add a small imaginary part t t + i where necessary. We are interested in the values of these integrals as t ranges over the real numbers. Using integration-by-parts identities [6, 63] we may express all integrals from this family as 3

14 linear combinations of a few master integrals. For the kite family there are eight master integrals [3, 33], which we denote as I =I,I,...,I 8 ). We will follow the notation of [33], where the definition of I -I 8 is given. The eight master integrals I for the kite system satisfy a system of differential equations in t of Fuchsian type [64 7], which has singularities at {,m,9m, }. 37) The system of differential equations reads µ d [ dt µ I = t A + µ t m A + µ t 9m A 9 ] I, 38) where A, A and A 9 are 8 8-matrices with entries of the form a+bε, where a,b Q and ε denotes the dimensional regularisation parameter. The explicit expressions are given in [33]. The system of differential equations in eq. 38) holds for all t C. In [33] we solved the system of differential equations in the Euclidean region i.e. region I) by performing a change of variables from t to the nome q. For this change of variables the relation between t and τ is given by t = fτ), fτ)= 9m ητ)4 η ) 3τ 4 η6τ) 4 η ) τ 4 ητ) 4 η3τ) 4. 39) The relation between τ and q is as usual q=expiπτ). In terms of the variable q the system of differential equations 4) becomes q d dq I = g, A + g, A + g,9 A 9 ) I, 4) where g,, g, and g,9 are modular forms of modular weight for the congruence subgroup Γ ), defined in [9]. We have seen in eq. 34) that the validity of eq. 39) is not restricted to region I, but holds for all values t on the path of fig. ). Thus eq. 4) holds for all values of q on the path of fig. 5). We already know a solution for this system of differential equations: In [33] we solved these differential equations in terms of ELi-functions in the region I. The solution extends to all values t on the path of fig. ). Furthermore from eq. 3) it follows that q < for t R\{m,9m, }. 4) This ensures the convergence of the q-series, which we use to express our results. Let us give an example: The first term of the Laurent expansion in the dimensional regularisation parameter ε of the equal-mass sunrise integral in two space-time dimensions reads S ),t) = 3ψ iπ { Li r 3 ) Li r 3 )+ELi ; r 3 ; ; q) ELi ; r 3 ; ; q ) }, 4) 4

15 Imk ) Imk ) Imk ) Rek ) Rek ) Rek ) C Figure 6: Deformation of the path C in k -space into a quarter-circle in the clockwise direction and a small full circle around in the counter-clockwise direction. where r 3 = expπi/3) and the functions ELi n;m x,y,q) are defined by ELi n;m x;y;q) = j= k= x j y k j n k m q jk. 43) Eq. 4) holds for all values of t, the periods ψ and ψ and in turn the nome q) are computed by eq. 6). Let us summarise: With the definition of the periods ψ and ψ given in eq. 6) the differential equation 4) holds for all values of t R + i. Thus any solution It), originally defined in one particular region say region I), extends to all values t R+i. Furthermore we have for t R+i q <, 44) ensuring the convergence of the ELi-series. At the singular points t {m,9m, } we have q =.) In other words: The analytic continuation of It) could not be simpler! 5 Monodromy and the Picard-Lefschetz theorem In this section we derive the monodromy relation of eq. 6) and eq. 7). This is a standard application of the Picard-Lefschetz theorem and we follow the textbooks [73, 74]. We may deform the path C in fig. 3) in k -space into a quarter-circle in the clockwise direction and a small full circle around in the anti-clockwise direction, as shown in fig. 6. For the monodromy we have to study the contribution from the small full circle around. Let us set λ=k. We recall that we may write the complete elliptic integral of the first kind equally well as Kk) = λ dt t λ t t. 45) 5

16 Thus, Kk) gives a quarter-period of the elliptic curve in Legendre form The roots of the polynomial px) are E λ : y = px), px)=xx λ)x ). 46), λ,,. 47) Let us now study the behaviour of Kk) as λ moves in a small circle around. It is slightly simpler to consider instead of the family of elliptic curves E λ the family of elliptic curves E ϕ given by where the roots are given by E ϕ : y = px,ϕ), px,ϕ) = xx e ϕ))x e ϕ)), 48) e ϕ)= re iϕ, e ϕ)=+re iϕ. 49) We study this family for small positive r and ϕ [,π]. We consider the periods P ϕ) = δ dx y, P ϕ) = δ dx y, y = x x e ϕ) x e ϕ), 5) where the two cycles δ and δ form a basis of H E ϕ,z). The orientation of the two cycles is such that for ϕ=we have P ) = e ) dx y = e ) dx y, P ) = e ) e ) dx y, 5) where the path of integration in x-space is such that we have an infinitesimal small negative imaginary part for x. We note that and px,) = px,π) = px,π), 5) e ) = e π) = e π), e ) = e π) = e π). 53) Thus, under a half-turn the polynomial px,ϕ) transforms into itself, however the roots e ϕ) and e ϕ) exchange their roles. Let us now discuss the effect of a half-turn on the periods P ϕ) and P ϕ). As px,ϕ) transforms into itself, we may express the periods P π) and P π) as a linear combination of the periods P ) and P ). As ϕ ranges over the interval[,π] the cycles of integration for the two periods transform as shown in fig. 7). We recall that we may think of an elliptic curve as two copies of an Riemann sphere, each sphere with two cuts. The elliptic 6

17 e ) e ) e ) e ) δ ) δ ) δ π) δ π) Figure 7: The left picture shows the cycles δ and δ for ϕ=, the right picture shows the cycles δ and δ for ϕ=π. Solid lines corresponds to a path on one Riemann sheet, dashed lines to a path on the other Riemann sheet. a) b) e ) e ) e ) e ) δ π) c) d) e ) e ) e ) e ) e) e ) e ) Figure 8: Deformation of the cycle δ π). In the step from c) to d) we bring the closed path around e ) from one Riemann sheet to the other Riemann sheet. This changes the orientation. The last figure shows that δ π)=δ ) δ ). 7

18 3 p 4 5 Figure 9: The kite graph. Solid lines correspond to massive propagators of mass m, dashed lines correspond to massless propagators. curve is obtained by gluing together the two spheres at the cuts with opposite orientations. From fig. 7) it is clear that P π) = P ). 54) We may deform the cycle of integration δ π) for P π) as shown in fig. 8). Thus P π) = P ) P ). 55) Eq. 55) is an application of the Picard-Lefschetz theorem. In the example above δ is a vanishing cycle. The Picard-Lefschetz theorem states that δ π) = δ ) δ ) δ )) δ ), 56) whereδ a δ b ) denotes the intersection number of the cycles δ a and δ b [73]. Combining two half-turns we obtain the monodromy relation of eq. 6) and eq. 7) from P π)=p ) P ), P π)=p ). 57) 6 Numerical results In this section we consider three examples of Feynman integrals from the family of kite integrals. The family of kite integrals is given in D-dimensional Minkowski space by I ν ν ν 3 ν 4 ν 5 D, p,m,µ ) = ) ν 345 µ ) ν 345 D d D k d D k with the propagators iπ D iπ D D ν Dν Dν 3 3 Dν 4 4 Dν 5 5, 58) D = k m, D = k, D 3 =k k ) m, D 4 =k p), D 5 =k p) m 59) and ν 345 = ν + ν + ν 3 + ν 4 + ν 5. The kite graph is shown in fig. 9). The internal momenta are denoted by k and k, the internal mass by m and the external momentum by p. We denote 8

19 Figure : The Feynman graphs of the examples considered in this section. A dot on a propagator indicates, that this propagator is raised to the power two. t = p. In the following we will suppress the dependence of the integrals on the mass m and the scale µ and we write I ν ν ν 3 ν 4 ν 5 D,t) = I ν ν ν 3 ν 4 ν 5 D,t,m,µ ). 6) The three examples which we consider in this section are the kite integral I 4 ε,t), the equal-mass sunrise integral in ε space-time dimensions I ε,t) and the integral bubble squared ) I 4 ε,t). The Feynman graphs for these integrals are shown in fig. ). All three Feynman integrals start in the ε-expansion at Oε ). In [33] we presented an algorithm to express all integrals from the kite family to all orders in the dimensional regularisation parameter in terms of ELi-functions, which are generalisations of multiple) polylogarithms towards the elliptic case. For the kite integral and the sunrise integral the occurrence of ELi-functions is generic, the bubble integral squared may be expressed in terms of harmonic polylogarithms [75, 76]. All polylogarithms occurring in the integrals from the kite family may be expressed in terms of ELi-functions. We include the trivial example of the bubble integral squared to show explicitly that our method includes the case where the ELi-functions degenerate to multiple polylogarithms. Expressed in terms of ELi-functions we have I 4 ε,t) = [ µ t ln t )] m + Oε) 6) = 9 µ4 [ E; t ;; q) E ; r 6 ;; q) ] + Oε), I ε,t) = 3ψ [ ) ] π Cl + E ; r 3 ; ; q) + Oε), π 3 I 4 ε,t) = [G,,;y) µ G,,;y)+ π t 6 G;y) ) ] π +7Cl E ; r 3 ;; q)+7e,;,; r 3,r 3 ;, ; q) + Oε). 3 The notation is as follows: We introduced the dimensionless variable y=t/m. The symbol r n 9

20 denotes the n-th root of unity r n = exp ) πi. 6) n The E-functions are linear combinations of the ELi-functions. Both are defined in appendix B. Cl denotes the Clausen function, defined by Cl ϕ) = i [ Li e iϕ ) Li e iϕ )]. 63) The harmonic polylogarithms G;y), G,,;y) and G,,;y) may be expressed in terms of ELi-functions. We have G;y)=3 [ E ; ;; q) E ; r 6 ;; q) ], 64) G,,;y)=9 [ E,;,;4, ;,; q) E,;,;4,r 6 ;,; q) E,;,;4 r 6, ;,; q)+ E,;,;4 r 6,r 6 ;,; q) ] 36 [ E,,;,,;, r 3,, ;,,; q) E,,;,,;, r 3,,r 6 ;,,; q) E,,;,,;, r 3,r 6, ;,,; q)+ E,,;,,;, r 3,r 6,r 6 ;,,; q) ], G,,;y)=9 [ E,;,;, ;,; q) E,;,;,r 6 ;,; q) E,;,; r 6, ;,; q)+ E,;,; r 6,r 6 ;,; q) ] 36 [ E,,;,,;,,r 3, ;,,; q) E,,;,,;,,r 3,r 6 ;,,; q) E,,;,,;, r 6,r 3, ;,,; q)+ E,,;,,;, r 6,r 3,r 6 ;,,; q) ]. In the following we set µ = m, express all integrals entirely in terms of ELi-functions and plot the real and the imaginary part of the ε -term of I 4 ε,t), I ε,t), I 4 ε,t) 65) as a function of y = t/m. Note that the expressions in eq. 6) are valid over the complete kinematic range t R. The periods ψ and ψ and in turn the nome q) are computed according to eq. 6). We compare our results to numerical results obtained from the SecDec-program [55 58]. The comparison is shown in fig. ) for y [ 5,5]. This range includes in particular the thresholds t = m and t = 9m. We find perfect agreement over the full range y R. The results for the bubble integral squared I i.e. the first two plots in fig. )) show that the numerical evaluation of this integral in terms of ELi-functions correctly reproduces the result of this integral, which alternatively may be expressed in terms of multiple polylogarithms. We note that truncating the q-series at Oq ) gives already a very good approximation, including the thresholds. In order to quantify this, we show in fig. ) for the kite integral the relative difference of the Oq )-approximation with the Oq )-approximation, defined by I I = I 4,t) q I 4,t) q I 4,t) q. 66)

21 8 bubble squared, real part our result SecDec 8 bubble squared, imaginary part our result SecDec ReI) ImI) t/m t/m 5 8 sunrise, real part our result SecDec 8 sunrise, imaginary part our result SecDec 6 6 ReI) 4 ImI) t/m t/m kite, real part our result SecDec 6 5 kite, imaginary part our result SecDec 4 4 ReI) 3 ImI) t/m t/m 5 Figure : The real and the imaginary parts of the ε -term of the Feynman integrals I 4 ε,t), I ε,t) and I 4 ε,t) as y = t/m ranges over the interval [ 5,5]. The dashed vertical lines indicate the thresholds at t = m and t = 9m.

22 .8 kite, approximation error, linear scale kite, approximation error, log scale I/I.4 I/I t/m t/m 5 Figure : The relative difference of the Oq )-approximation with the Oq )-approximation for the kite integral I. The left plot shows the relative difference on a linear scale, the right plot shows the relative difference on a logarithmic scale. We evaluate the integral in the interval t/m [ 5,5] at steps δt/m =., excluding the singular points t = m and t = 9m. Thus the closest points to the singular points are t =.98m, t =.m, t = 8.98m and 9.m. We see that the relative difference is rather small over a wide range of t, taking a maximum I 5% around t = 9m, where the kite integral is continuous. To see the behaviour at the threshold t = m, we show in the right plot of fig. ) the relative difference on a log-scale. The relative difference around the threshold t = m is I 6. The integrals I and I show a similar behaviour. Again, the maximum of the relative difference for t/m [ 5,5] is reached around t = 9m, being of the order of %. The ELi-functions provide therefore a fast and efficient way to evaluate these integrals over the complete kinematic range with the exception of small neighbourhoods around the singular points t {m,9m, }, where q =. 7 Conclusions In this paper we studied the analytic continuation and the numerical evaluation of Feynman integrals from the kite family. In previous publications we showed that these integrals can be expressed to all orders in the dimensional regularisation parameter ε in a neighbourhood of t = in terms of elliptic generalisations of multiple) polylogarithms, denoted as ELi-functions. In this paper we showed that these expressions hold in the full kinematic range t R after analytic continuation of the two periods. The ELi-functions are power series in the variable q. Furthermore we showed that q holds in the full kinematic range t R and q = is attained only at the singular points t {m,9m, }. Therefore the q-series expansion of the ELi-functions provide a fast and efficient way to evaluate these integrals over a wide kinematic range. We compared results from this method to numerical results from the SecDec program and found perfect agreement.

23 Acknowledgements C.B. and A.S. thank Konrad Schultka for useful discussions. S.W. would like to thank Hubert Spiesberger for useful discussions. C.B. is supported by Deutsche Forschungsgemeinschaft under the project BO45/-. A Convention regarding the roots of the cubic equation In section 3 we defined the roots e, e, e 3 of the cubic polynomial 4x 3 g x g 3 in eq. 5). The roots e and e were given by e = t 4µ 4 + 6m t+ 3m m t ) 3 9m t ) ), e = 4µ 4 t + 6m t+ 3m 4 3 m t ) 3 9m t ) ). 67) In this appendix we discuss the consequences of an alternative convention for the roots of the cubic equation, where e and e are replaced by ) ẽ = 4µ 4 t + 6m t+ 3m m t) 3 9m t), ) ẽ = 4µ 4 t + 6m t+ 3m 4 3 m t) 3 9m t), 68) respectively. Let us first note that for t < m and 9m < t one finds that e agrees with ẽ and e agrees with ẽ. Therefore the differences between the convention in eq. 67) and the one of eq. 68) are restricted to the region m < t < 9m. Let us first consider e and e. We recall that we assume that a infinitesimal small positive imaginary part is added to the variable t. We have m t i ) 3 = i t m ) 3 for t > m. 69) Thus Ime ) >, Ime ) <, for m < t < 9m. 7) Let us now consider ẽ and ẽ. Working to first order in the infinitesimal small imaginary part one obtains m t i) 3 9m t i) = m t) 3 9m t) 4m t) 7m t)i. Thus Imẽ ) { <, m < t < 7m, >, 7m < t < 9m, Imẽ ) { >, m < t < 7m, <, 7m < t < 9m, 7) 3

24 In the interval m < t < 7m one has e = e and ẽ = e. This just exchanges the roles of e and e. However, at t = 7m the roots ẽ and ẽ are not continuous functions of t. For this reason it is better to use the definitions of the roots as given in eq. 67). The discontinuity in t at 7m has the following origin: For 7m < t < 9m the imaginary part coming from 9m t i) is larger and has opposite sign as the imaginary part coming from m t i). Similar considerations apply to the algebraic prefactor in eq. 3), where we express the periods ψ and ψ in terms of complete elliptic integrals of the first kind. We defined the algebraic prefactor as and not as 4µ, 7) m t) 3 4 9m t) 4 4µ [ ] m t) 3. 73) 9m 4 t) Both definitions agree for t < m, but the latter leads to spurious discontinuities in the analytic continuation. B Elliptic generalisations of polylogarithms In this appendix we define the ELi-functions and the E-functions. The latter are just linear combinations of the ELi-functions. Let us start with the ELi-functions. These are functions of l+ ) variables x,..., x l, y,..., y l, q and3l ) indices n,..., n l, m,..., m l, o,..., o l. For l = we set ELi n;m x;y;q) = j= k= x j j n y k k m q jk. 74) For l > we define ELi n,...,n l ;m,...,m l ;o,...,o l x,...,x l ;y,...,y l ;q)= =... j = j l = k =... k l = x j j n... x j l l j n l l y k k m... yk l l k m l l q j k j l k l l j i k i j l k l ) o i i=. 75) We have the relations ELi n ;m x ;y ;q)eli n,...,n l ;m,...,m l ;o,...,o l x,...,x l ;y,...,y l ;q)= = ELi n,n,...,n l ;m,m,...,m l ;,o,...,o l x,x,...,x l ;y,y,...,y l ;q) 76) 4

25 and q dq q ELi n,...,n l ;m,...,m l ;o,o,...,o l x,...,x l ;y,...,y l ;q ) = = ELi n,...,n l ;m,...,m l ;o +),o,...,o l x,...,x l ;y,...,y l ;q). 77) It will be convenient to introduce abbreviations for certain linear combinations, which occur quite often. We define a prefactor c n and a sign s n, both depending on an index n by c n = [+i)+ i) )n ]= For l = we define the linear combinations {, n even, i, n odd, s n = ) n = {, n even,, n odd. 78) More explicitly, we have E n;m x;y;q) = c n+m i E n;m x;y;q) = [ ELin;m x;y;q) s n+m ELi n;m x ;y ;q )]. 79) i [ ELin;m x;y;q) ELi n;m x ;y ;q )], n+m even, ELi n;m x;y;q)+eli n;m x ;y ;q ), For l > we proceed as follows: For o = we set E n,...,n l ;m,...,m l ;,o,...,o l x,...,x l ;y,...,y l ;q)= For o > we set recursively n+m odd. 8) = E n ;m x ;y ;q)e n,...,n l ;m,...,m l ;o,...,o l x,...,x l ;y,...,y l ;q). 8) E n,...,n l ;m,...,m l ;o,o,...,o l x,...,x l ;y,...,y l ;q)= = q dq q E n,...,n l ;m,...,m l ;o ),o,...,o l x,...,x l ;y,...,y l ;q ). 8) The E-functions are linear combinations of the ELi-functions with the same indices. More concretely, an E-function of depth l can be expressed as a linear combination of l ELi-functions. We have E n,...,n l ;m,...,m l ;o,...,o l x,...,x l ;y,...,y l ;q)= 83) [ l c n j +m j ) t =... j sn j= i j +m j ]ELi n,...,n l ;m,...,m l ;o,...,o l x s t,...,xs t l l ;y s ) t,...,ys t l l ;q. t = t l = 5

26 C The arithmetic-geometric mean In this appendix we review the numerical evaluation of the complete elliptic integral of the first kind with the help of the arithmetic-geometric mean. Let a and b be two complex numbers. For n N one sets The sign of the square root is chosen such that [77] a n+ = a n+ b n ), b n+ = ± a n b n. 84) a n+ b n+ a n+ + b n+, 85) and in case of equality one demands in addition ) bn+ Im a n+ >. 86) The sequences a n ) and b n ) converge to a common limit lim a n = lim b n = agma,b ), 87) n n known as the arithmetic-geometric mean. The complete elliptic integral of the first kind is given by References Kk) = π agmk,), k = [] A. B. Goncharov, Math. Res. Lett. 5, ). [] A. B. Goncharov, ), math.ag/359. k. 88) [3] J. M. Borwein, D. M. Bradley, D. J. Broadhurst, and P. Lisonek, Trans. Amer. Math. Soc. 353:3, 97 ), math.ca/9945. [4] S. Moch, P. Uwer, and S. Weinzierl, J. Math. Phys. 43, 3363 ), hep-ph/83. [5] J. Vollinga and S. Weinzierl, Comput. Phys. Commun. 67, 77 5), hep-ph/459. [6] T. Gehrmann and E. Remiddi, Comput. Phys. Commun. 4, 96 ), hep-ph/773. [7] T. Gehrmann and E. Remiddi, Comput. Phys. Commun. 44, ), hep-ph/55. [8] D. Maitre, Comput. Phys. Commun. 74, 6), hep-ph/575. [9] D. Maitre, Comput. Phys. Commun. 83, 846 ), arxiv:hep-ph/735. 6

27 [] H. Frellesvig, D. Tommasini, and C. Wever, JHEP 3, 89 6), arxiv: [] D. J. Broadhurst, J. Fleischer, and O. Tarasov, Z.Phys. C6, ), arxiv:hepph/ [] F. A. Berends, M. Buza, M. Böhm, and R. Scharf, Z.Phys. C63, 7 994). [3] S. Bauberger, M. Böhm, G. Weiglein, F. A. Berends, and M. Buza, Nucl.Phys.Proc.Suppl. 37B, ), arxiv:hep-ph/ [4] S. Bauberger, F. A. Berends, M. Böhm, and M. Buza, Nucl.Phys. B434, ), arxiv:hep-ph/ [5] S. Bauberger and M. Böhm, Nucl.Phys. B445, 5 995), arxiv:hep-ph/95. [6] M. Caffo, H. Czyz, S. Laporta, and E. Remiddi, Nuovo Cim. A, ), arxiv:hepth/9858. [7] S. Laporta and E. Remiddi, Nucl. Phys. B74, 349 5), hep-ph/466. [8] B. A. Kniehl, A. V. Kotikov, A. Onishchenko, and O. Veretin, Nucl. Phys. B738, 36 6), arxiv:hep-ph/535. [9] S. Groote, J. G. Körner, and A. A. Pivovarov, Annals Phys. 3, 374 7), arxiv:hepph/5686. [] S. Groote, J. Körner, and A. Pivovarov, Eur.Phys.J. C7, 85 ), arxiv: [] D. H. Bailey, J. M. Borwein, D. Broadhurst, and M. L. Glasser, J. Phys. A4, 53 8), arxiv:8.89. [] S. Müller-Stach, S. Weinzierl, and R. Zayadeh, Commun. Num. Theor. Phys. 6, 3 ), arxiv:.436. [3] L. Adams, C. Bogner, and S. Weinzierl, J. Math. Phys. 54, 533 3), arxiv:3.74. [4] S. Bloch and P. Vanhove, J. Numb. Theor. 48, 38 5), arxiv: [5] L. Adams, C. Bogner, and S. Weinzierl, J. Math. Phys. 55, 3 4), arxiv: [6] L. Adams, C. Bogner, and S. Weinzierl, J. Math. Phys. 56, 733 5), arxiv: [7] L. Adams, C. Bogner, and S. Weinzierl, J. Math. Phys. 57, 334 6), arxiv: [8] S. Bloch, M. Kerr, and P. Vanhove, 6), arxiv:

28 [9] L. Adams and S. Weinzierl, 7), arxiv: [3] E. Remiddi and L. Tancredi, Nucl.Phys. B88, 343 4), arxiv: [3] A. Sabry, Nucl. Phys. 33, 4 96). [3] E. Remiddi and L. Tancredi, Nucl. Phys. B97, 4 6), arxiv:6.48. [33] L. Adams, C. Bogner, A. Schweitzer, and S. Weinzierl, J. Math. Phys. 57, 3 6), arxiv: [34] S. Bloch, M. Kerr, and P. Vanhove, Compos. Math. 5, 39 5), arxiv: [35] A. Primo and L. Tancredi, 7), arxiv: [36] R. Bonciani et al., JHEP, 96 6), arxiv: [37] A. von Manteuffel and L. Tancredi, 7), arxiv: [38] S. Bloch, H. Esnault, and D. Kreimer, Comm. Math. Phys. 67, 8 6), math.ag/5. [39] C. Bogner and S. Weinzierl, J. Math. Phys. 5, 43 9), arxiv: [4] G. Passarino, European Physical Journal C 77, 77 7), arxiv:6.67. [4] Z. Nagy and D. E. Soper, JHEP 9, 55 3), hep-ph/387. [4] C. Anastasiou, S. Beerli, and A. Daleo, JHEP 5, 7 7), hep-ph/738. [43] W. Gong, Z. Nagy, and D. E. Soper, Phys. Rev. D79, 335 9), arxiv: [44] S. Becker, C. Reuschle, and S. Weinzierl, JHEP, 3 ), arxiv:.487. [45] S. Becker, C. Reuschle, and S. Weinzierl, JHEP 7, 9 ), arxiv:5.96. [46] S. Becker and S. Weinzierl, Phys.Rev. D86, 749 ), arxiv: [47] S. Becker and S. Weinzierl, Eur.Phys.J. C73, 3 3), arxiv:.59. [48] S. Catani, T. Gleisberg, F. Krauss, G. Rodrigo, and J.-C. Winter, JHEP 89, 65 8), arxiv: [49] R. J. Hernandez-Pinto, G. F. R. Sborlini, and G. Rodrigo, JHEP, 44 6), arxiv: [5] S. Buchta, G. Chachamis, P. Draggiotis, and G. Rodrigo, Eur. Phys. J. C77, 74 7), arxiv:

29 [5] G. F. R. Sborlini, F. Driencourt-Mangin, R. Hernandez-Pinto, and G. Rodrigo, JHEP 8, 6 6), arxiv: [5] T. Binoth and G. Heinrich, Nucl. Phys. B585, 74 ), hep-ph/43. [53] C. Bogner and S. Weinzierl, Comput. Phys. Commun. 78, 596 8), arxiv: [54] A. V. Smirnov and V. A. Smirnov, JHEP 5, 4 9), arxiv:8.47. [55] S. Borowka et al., Comput. Phys. Commun. 96, 47 5), arxiv: [56] S. Borowka, J. Carter, and G. Heinrich, Comput. Phys. Commun. 84, 396 3), arxiv:4.45. [57] S. Borowka and G. Heinrich, Comput. Phys. Commun. 84, 55 3), arxiv: [58] S. C. Borowka, Evaluation of multi-loop multi-scale integrals and phenomenological twoloop applications, PhD thesis, Munich, Tech. U., 4, arxiv: [59] The British Standards Institution, The C Standard Wiley, 3). [6] C. Bogner and S. Weinzierl, Int. J. Mod. Phys. A5, 585 ), arxiv: [6] F. Olver, D. Lozier, R. Boisvert, and C. Clark, NIST Handbook of Mathematical Functions Cambridge University Press, ). [6] F. V. Tkachov, Phys. Lett. B, 65 98). [63] K. G. Chetyrkin and F. V. Tkachov, Nucl. Phys. B9, 59 98). [64] A. V. Kotikov, Phys. Lett. B54, 58 99). [65] A. V. Kotikov, Phys. Lett. B67, 3 99). [66] E. Remiddi, Nuovo Cim. A, ), hep-th/9788. [67] T. Gehrmann and E. Remiddi, Nucl. Phys. B58, 485 ), hep-ph/9939. [68] M. Argeri and P. Mastrolia, Int. J. Mod. Phys. A, ), arxiv: [69] S. Müller-Stach, S. Weinzierl, and R. Zayadeh, Commun.Math.Phys. 36, 37 4), arxiv: [7] J. M. Henn, Phys. Rev. Lett., 56 3), arxiv: [7] J. M. Henn, J. Phys. A48, 53 5), arxiv:4.96. [7] L. Adams, E. Chaubey, and S. Weinzierl, Phys. Rev. Lett. 8, 46 7), arxiv:

30 [73] J. Carlson, S. Müller-Stach, and C. Peters, Period Mappings and Period Domains Cambridge University Press, 3). [74] R. C. Hwa and V. L. Teplitz, Homology and Feynman Integrals W. A. Benjamin, 966). [75] J. A. M. Vermaseren, Int. J. Mod. Phys. A4, ), hep-ph/9868. [76] E. Remiddi and J. A. M. Vermaseren, Int. J. Mod. Phys. A5, 75 ), hep-ph/ [77] D. A. Cox, Enseign. Math. 3, ). 3