arxiv: v1 [math-ph] 18 Oct 2016

Transcription

1 Elliptic Polylogarithms and Basic Hypergeometric Functions Giampiero Passarino a, a Dipartimento di Fisica Teorica, Università di Torino, Italy INFN, Sezione di Torino, Italy Abstract arxiv: v [math-ph] 8 Oct 206 Multiple elliptic polylogarithms can be written as multiple) integrals of products of basic hypergeometric functions. The latter are computable, to arbitrary precision, using a -difference euation and -contiguous relations. Keywords: Feynman integrals; elliptic polylogaritms; basic hypergeometric functions. Introduction There is a wide class of Feynman integrals, mainly related to massless theories, which can be expressed in terms of multiple polylogarithms. More challenging are Feynman integrals, which cannot be expressed in terms of multiple polylogarithms. Evaluating this integrals one encounters elliptic generalizations of multiple) polylogarithms EP); examples can be found in Refs. [ 3] and in Refs. [4, 5]. From a more abstract point of view on multiple) elliptic polylogarithms, in particular on their analytic structure, see Refs. [6, 7] and Refs. [8 ]. An interesting problem is to find a suitable integral representation and the analytic continuation of EPs and an efficient algorithm for their numerical evaluation. 2. Elliptic polylogarithms In Refs. [ 3] the following functions, of depth two, are defined: ELi n;m x, y; ) = j= k= x j y k j n k m jk. ) For the sake of simplicity we shall assume that arg) = 0. The results derived can be extended by using analytic continuation. Following Ref. [7] we prefer to start from ELi n;m x, y; ) = x j j= j n Li m y j ), 2) where Li m z) is a generalized polylogarithm [2]. Furthermore, C x with < and y C with / R y. Thus, for x < the series converges absolutely. It is immediately seen that ELi 0;0 x, y; ) has a pole at x =. address: giampiero@to.infn.it Giampiero Passarino) Work supported by the Research Executive Agency REA) of the European Union under the Grant Agreement PITN-GA HiggsTools). October 2, 206

2 Note that Ref. [7] defines functions whose simplest example is given by reuiring < x <. The elliptic dilogarithm is defined in Ref. [5] as Ê 2 x) = n=0 and we have the following relations: where the last relation holds for z <. Ex, y; ) = x j Li y j ), 3) j Z [ ] Li 2 n x) Li 2 n x) n= Li 2 n z) = ELi 0;2, z; ), n= Λ 2 iz) = [ ] Li 2 z) Li 2 z) 2i = [Li 2 n n=0 x ) ) ] Li 2 n, 4) x ) n 2n + ) 2 iz)2n+, 5) Next, generalizations of higher depth of E.2) are also defined in Ref. [3]. They are defined by ELi [n]l ;[m] l ;[2σ] l [x] l ; [y] l ; ) =... j = j l = k =... j k j l k l k l = which are elliptic polylogarithms of depth 2l. We have introduced the abbreviation etc. The following relation holds: Furthermore, for l >, we can use x j j n... x j l j n l l y k k m... yk l k m l l l i= j i k i j l k l ) σ i, 6) [n] l = n,..., n l 7) dz dz ELi n;m x, y; ) = z ELi n ;m zx, y; ) = z ELi n;m x, zy; ). 8) 0 dz z ELi [n] l ;[m] l ;[2σ] l [x] l ; [y] l ; z) = ELi [n]l ;[m] l ;2σ +)...2σ l [x] l ; [y] l ; ), 9) ELi [n]l ;[m] l ;0...2σ l [x] l ; [y] l ; ) = ELi n ;m x, y ; ) ELi n2...n l ;m 2...m l ;2σ 2...2σ l x 2... x l ; y 2... y l ; ). 0) In Es.8) 0) we find recurrence relations where the starting point is always ELi 0;0 x, y; ), e.g. dz ELi n,n 2 ;m,m 2 ;2 x, x 2 ; y, y 2 ; ) = 0 z ELi n ;m x, y ; z) ELi n2 ;m 2 x 2, y 2 ; z), ) followed by repeated applications of E.8). In the next Section we will define basic hypergeometric series and establish the connection with elliptic polylogarithms. 2

3 3. Basic hypergeometric functions It is easy to show that ELi 0;0 x, y; ) can be written in terms of a basic hypergeometric function. Indeed, ELi 0;0 x, y; ) = = x j y j j= y j = xy y 2 φ, y; y 2 ;, x ) y j x j j= x j = xy x 2 φ, x; x 2 ;, y ), 2) where the basic hypergeometric series is defined as [ ] α,... α r ; z α ; ) n... α r ; ) [ n rφ s = β... β s ) n nn )/2] +s r z n, β ; ) n... β s+ ; ) n n=0 whit β s+ = and the -shifted factorial defined by α ;,) n = n k=0 α k). 3) The basic hypergeometric series was first introduced by Heine and was later generalized by Ramanujan. For the case of interest we will use the shorthand notation [ ] a, b; z φa, b; c;, z) = 2 φ. c Furthermore, we introduce Φx, y; ) = φ, y; y 2 ;, x ). 4) If < the φ series converges absolutely for z <. The series also converges absolutely if > and z < c / ab. With a =,b = y,c = y 2 and z = x we have that the basic hypergeometric series converges absolutely if < and x <, > and x <, which is a special case of z < c / ab. In the following we will always assume that <. Indeed, using p = / and a; ) n = a ; p) n a) n p nn )/2, 5) one obtains the following relation, φa, b; c;, z) = φ a, b ; ) ab ; p, c c pz. 6) Using E.6) we obtain and the inversion formula ELi 0;0 x, y ; ) = ELi 0;0 x, y; ) = xy ) p, p y φ py p2 ; y ; p, x, p =, 7) ) xy p ) Φ x p, y; p = ) y ELi 0;0 x p, y; p, p =. 8) The next problem is the continuation of the φ series into the complex z plane and the extension to complex inside the unit disc. We have two alternatives, analytic continuation and recursion relations. First we present an auxiliary relation that will be useful in the continuation of φ 3

4 3.. -contiguous relations There are several -contiguous relations for φ and one will be used extensively in the rest of this work, see E. A.0) of Ref. [3] see also.0 of Ref. [4] 2 ); using a shorthand notation, i.e. one derives φ = 2 φ a, b; c;, z) φc ±n ) = 2 φ a, b; c ±n ;, z ), 9) c) c)c abz)φ = c) { [ ] } 2 c c) + cb + a) + )ab z φc + ) + c a)c b)zφc +2 ). 20) 3.2. Analytic continuation A detailed discussion of the analytic continuation of basic hypergeometric series can be found in Ref. [4] Sects ) and in Ref. [6] Theorem 4.). From Ref. [4] we can use the following analytic continuation 2φ a, b; c;, z) = b, c/a; ) az, /az; ) c, b/a; ) z, /z; ) 2φ a, a/c; a/b;, c/abz) + a, c/b; ) bz, /bz; ) c, a/b; ) z, /z; ) 2φ b, b/c; b/a;, c/abz), 2) where arg z) < π, c and a/b are not integer powers of and a,b,z 0. Furthermore, a, a 2 ; ) = a ; ) a 2 ; ), a; ) = k=0 a k). 22) It is worth noting that the coefficients in E.2) are reducible to φ 0 functions, e.g. c/a; ) c/a = c; ) k ) c k = φ 0 a ;, c, 23) etc. A more convenient way to compute the -shifted factorial is given by z : ) = k=0 n=0 nn /2 ; ) n z) n, z <. 24) The condition arg z) < π must be understood as follows: The Mellin transform of φ for z R,z > 0) is a Cauchy Principal Value. With < the two series in E.2) can be used for z > z 0 with z 0 = c/ab). If z < z 0 we can use E.20); repeated applications transform c into n c until a value of n is reached for which z n = n+ c/ab) is such that z > z n. In Sect. 4.8) of Ref. [4] it is shown that extension to complex inside the unit disk is possible, provided that arg z) ω 2 ω ln z < π, 25) where ln = ω iω2. As seen in the complex mrz plane the condition is represented by a spiral of euation r = exp{ ω ω 2 θ}, z = r e θ. 26) It remains to study the convergence for =. In Ref. [7] a condition is presented so that the radius of convergence is positive; furthermore, the numbers with positive radius are densely distributed on the unit circle. Obviously, in our case a =, b = y and c = y 2 ), = reduces any ELi function to a product of polylogarithms. 2 This relation was introduced by E. Heine [5] in

5 3.3. Basic hypergeometric euation In the context of functional euations the basic hypergeometric series provides a solution to a second order - difference euation, called the basic hypergeometric euation, see Refs. [6, 8]. We will now show how φ can be computed to arbitrary precision, using a theorem proved in Ref. [9]. Theorem 3. Chen, Hou, Mu). Let f z) be a continuous function defined for z < r and d 2 an integer. Suppose that f z) = d n= with a n 0) = w n. Suppose that exists a real number M > 0 such that and a n z) f z n ), 27) a n z) w n M z, n d, 28) w d + w d + w d +... w w d <, w w d =. 29) The f z) is uniuely determined by f 0) and the functions a n z). With we have Fz) = Fz) = 2 φ a, b; c;, z) 30) a + b)z c Fz) + c abz z ) z ) Fz2 ), 3) thus, by the theorem and for < and c/ < from E.29)) we can determine uniuely Fz) by F0) and the -difference euation E.3)), i.e. we define with A i) 0 = a iz) and obtain { A i) High accuracy can be obtained by computing n+ = a iz n+ )A ) n A 2) n+ = a 2z n+ )A ) n Fz) = F0) F N z) = F0) + A i+) n, i < 2 2 lim n Ai) i= 2 i= for N high enough. Application to 2 φ, y; y 2 ;, x ) reuires y <. If > we can use the -difference euation downword. In this case Fz) = and { A i) n. 32) A i) N, 33) c + ) a + b)z F z) + z 2 c abz c abz F 2 z), 34) n+ = a iz n )A ) n A 2) n+ = a 2z n )A ) n + A i+) n, i < 2 which reuires /c <. In both cases the -difference euation determines F without limitations in the z complex plane. 5

6 To summarize, we can compute φ by using the sum of the series inside the circle of convergence and their analytic continuation) or by using the -difference euation. The advantage in using the latter is no limitation on z but < and c/ < or > and /c <. 35) To be precise, the function defined by Es.3) 35) is a meromorphic continuation of φ with simple poles located at z = n, n Z. Seen as a function of α = lnz)/2iπ) the continuation shows poles at 2iπ α = nω where ω = ln. What to do when we are outside the two regions of applicability? Instead of using analytic continuation we can do the following: assume that < but c/ >. We can use E.20); repeated applications transform c/ into n c until a value of n is reached for which n c <. Similar situation when > and /c > where repeated applications transform /c into n c until a value of n is reached for which n c <. It is worth noting that the Φ function satisfies [ 2 x ] Φx, y; ) = x + y Φx, y; ) y 2 x x Φ x 2, y; ) = x y 2 Φx, y; ) +, 36) y which are the -shift along x and y Poles Using E.2) and E.3) we conclude that ELi 0;0 x, y; ) has simple poles located at x = n, y = n, n Z. 37) Poles in the complex y-plane can also be analyzed by using E.2) and E.20). Indeed, from E.20) we obtain Φx, y; ) = [ + y x + x x 2 y x + x + ]φ, y; y 3 ;, x ) 2 3 y φ, y; y 4 ;, x ). 38) ) Using Φ x, ; = we have a simple pole of ELi 0;0 x, y; ) at y = / with residue x/ the pole at x = / has residue y/). Repeated applications of the -contiguous relation exhibit the poles at y = n,n Z. Residues can be computed according to the following chain Φx, y; ) = R x + Reg x, y; ) = R 2 x + Reg 2 x, y; ) = etc., 39) 2 where the regular part admits a Taylor expansion around x = /, x = / 2 etc. R = ) Φ, y; ) + y ) Φ, y; ), R 2 = y ) Φ, y; ) + y 2 ) Φ, y; ), = etc. 40) Using the series we obtain ) Φ k, y; = y) n=0 nk+), 4) yn+ R = y ), R = y ) y, etc. 42) 6

7 showing the following residues for ELi 0;0 x, y; ): Res x=/ n ELi 0;0 x, y; ) = y, y2 2, y3,... 3 for n =, 2, 3,... 43) The isolation of simple poles in ELi 0;0 x, y; ) is crucial in order to compute Elliptic polylogarithms of higher depth. 4. Elliptic polylogarithms of higher depth In this Section we show how to compute arbitrary elliptic polylogarithms, in particular how to identify their branch points their multi-valued component). The procedure is facilitated by the fact that both the basic hypergeometric euation and the -contiguous relation allows to isolate the simple) poles of ELi 0;0 x, y; ) with a remainder given by a + distribution. Introducing the usual i ε prescription we obtain a general recipe for computing elliptic polylogarithms on the cuts x and or y real and greater than /). In the following we discuss few explicit examples. ELi ;0 x, y; ) From E.8) we obtain ELi ;0 x, y; ) = xy dzφzx, y; ). 44) y 0 For x R the integral is defined when x < /, otherwise it is understood that x x ± iε where ε 0 +. Consider the case < x < 2. 45) From E.3) we derive t i z) Φzx, y; ) = i=,2 zx Φ zx i, y; ), 46) From E.45) it follows that no pole of the two φ functions in E.46) appears for z [0,]; therefore the two φ in the r.h.s of E.46) can be evaluated according to the strategy outlined in the previous Sections. Using ) ) t =, t 2 = y ), 47) x x we can write where the subtraction term is and x = x ± iε. Note that 0 dzφzx, y; ) = [ ] x ln x) Φ, y; ) yφ, y; ) + S, 48) S = 0 dz zx [ t i z)φ zx i, y; ) t i i=,2 x )Φ i, y; )], 49) Φ, y; ) yφ, y; ) = y. 50) If 2 < x < 3 we can iterate once more obtaining an additional ln x 2 ) etc. The function S defined in E.49) is a + distribution which has simple poles in the y-plane. The explicit result is as follows: ELi ;0 x, y; ) = ELic ;0 x, y; ) + ELir ;0 x, y; ), 7 < x < 2, 5)

8 where the cut part is while the restr is The second iteration gives ELir ;0 x, y; ) = ELic ;0 x, y; ) = y ln x), 52) xy [ F u) + y ] du + yφxu, y; ), y 0 xu F u) = x 2 u ) [ Φxu, y; ) yφ x 2 u, y; )]. 53) ELi ;0 x, y; ) = ELic 2 ;0 x, y; ) + ELir 2 ;0 x, y; ), 2 < x < 3, 54) F 2 u) = ELir 2 ;0 x, y; ) = [ + + y) ELic 2 ;0 x, y; ) = y ln x) y 2 ln x 2 ), 55) ] xy [F 2 u) F 2 x du ) y 0 xu ) y 2 + Φ x 4, y; ), F 22 u) = y ) [ Φ x 2, y; ) yφ x 3, y; )], R 2 u) = y) 2 Φ x 2, y; ) [ y + + y 2 + We continue with other examples. ELi 0; x, y; ) + F22u) F 22 ) ] x 2 xu + R 2 u), 56) 2 Φ x 2, y; ) {[ ] y + + ) + y) )} Φ x 3, y; ) xu y) ]Φ x 3, y; ) 3 ) xu + 4 Φ x 4, y; ). 57) 3 Similar results follow from E.2), i.e. by using the x, y symmetry of ELi 0;0 x, y; ). ELi 2;0 x, y; ) We can use the same derivation as before obtaining a result similar to the one in E.48) where we replace ln x) Li 2 x), 58) where Li 2 z) is the dilogarithm; for ELi n;0 x, y; ) with n > 2 we obtain polylogarithms. The explicit result is as follows: ELi 2;0 x, y; ) = ELic 2;0 x, y; ) + ELir 2;0 x, y; ), < x < 2, 59) where the cut part is ELic 2;0 x, y; ) = yli 2 x), 60) 8

9 while the rest part is ELir 2;0 x, y; ) = xy dv v [ F u) + y ] du + yφxu, y; ), y 0 v 0 xu F u) = x 2 u ) [ Φxu, y; ) yφ x 2 u, y; )], 6) When both n and m are different from zero the derivation reuires isolating x-poles and y-poles. ELi ; x, y; ) This case reuires additional work. Consider the integral dz dz 2 z 2 y Φz x, z 2 y; ). 62) 0 When y, R the integral is defined for y < /, otherwise it is understood that y y ± oε with ε 0 +. By using E.20) we obtain [ ] { [ φc) = c) c)c abu) c) +b + a)cu + )abu + c)c 2] φc) } + c a)c b)uφc 2 ), 63) where φc) = φa, b; c;, u). Next we replace a =, b = z 2 y, c = z 2 y 2 and u = z x. If / < y < / 2 we replace z 2 = /y) in the φ function of E.62) and obtain a logarithmic part as well as a subtraction part. Note that b = gives a terminating series, i.e. ELi 0,0;0,0;2 This function is defined through ln y), 64) y φ, ; 3 ;, z x ) =. 65) dz ELi 0,0;0,0;2 x, x 2 ; y, y 2 ; ) = 0 z ELi 0;0 x, y ; z) ELi 0;0 x 2, y 2 ; z). 66) The integrand has poles in the complex z-plane located at x i n+ z = y i m+ i =,2 n,m Z Poles of the first series are isolated by using the -difference euation while those in the second series are isolated by using the -contiguous relation. All poles are simple as long as none of the ratios x i /x j,y i /y j and x i /y j is eual to an integer power of. 4.. Mixed hypergeometric series When < and x < / we can write ELi n;0 x, y; ) or y < / and ELi 0;n x, y; )) as mixed hypergeometric series, e.g. ELi ;0 x, y; ) = 2,2 φ,, y;, ; y 2, 2; x ) y; ) n ) n ) n x) n = y 2, 67) ; ) n 2) n n! where a) n = Γa + n)/γa) is the Pochhammer symbol. For a previous definition of mixed hypergeometric series see Ref. [20]. 9 n=0

10 4.2. Barnes contour integrals For 0 < < we can write ELi ;0 x, y; ) = xy 2π i Γ π ds sin πs x) s s +, 68) ys+ where x < and arg x) < π. The contour of integration, denoted by Γ, runs from i to +i so that the poles at s Z lie to the right of the contour and the other poles, at s Z and s = + lny + 2mπ i)/omega with ω = ln and m Z, lie to the left and the latter are at least some ε ε 0 + ) distance away from the contour. The r.h.s. of E.68) defines an analytic function of z = x in arg z) < π. Note that E.68) can be generalized to define the analytic continuation of ELi n;0 x, y; ) and can be extended to complex inside the unit disc see E.26)) Eisenstein-Kronecker series The construction of elliptic multiple polylogarithms in Ref. [7] is largely based on the Eisenstein-Kronecker series Fξ, α, τ) defined in their Sect. 3.4; with z = eξ ), = eτ), α = eu), 69) where eξ ) = exp{2π iξ } we obtain the following relation with Φ, the basic hypergeometric series of E.4): Fξ, α, τ) = z 2π i z uz z Φu, z; ) + z uz ) Φ u, z; ). 70) The function F satisfies i.e. uasi-periodicity. Fξ +, α, τ) = Fξ, α, τ), Fξ + τ, α, τ) = Fξ, α, τ), 7) α 5. Conclusions We have established a connection between elliptic polylogarithms and basic hypergeometric functions, providing a framework for high-precision numerical evaluation of the Eli functions. Outside the region of convergence of the series the numerical evaluation uses analytic continuation via Watson s contour integral representation, basic hypergeometric euation and -contiguous relations. As an example we show in Fig. the behavior of Φx, y; ), E.4), around the poles at x = / blue curve) and x = / 2 red curve) for y = 0. and = 0.9. In Fig. 2 we show the real and the imaginary parts of Φx, y; ) for y = 0. and = i0.04. In Fig.3 we show eli ;0, defined by ELir ;0 x, y; ) = for different values of y and and / r < x < / 2 r, r = Re. xy y eli ;0 x, y; ), 72) 0

11 x 2 ) Φx, y ; ) y = 0. = 0.9 x ) Φx, y ; ) 0.37 x 0.94 Figure : Behavior of Φx, y; ), defined in E.4), around the poles at x = / blue curve) and x = / 2 red curve).

12 +25 Φx, y ; ) y = 0. = i 0.04 Im Re 0 0 x 2 Figure 2: Behavior of Φx, y; ), defined in E.4), as a function of x for y = 0. and = i

13 4 y = 0. = i 0.04 y = 0. = 0.9 y = 0.5 = 0.9 eli ; 0x, y ; ) 2 y = 0.9 = Im Re x Figure 3: The function eli ;0 x, y; ), E.72), for different values of y and and / r < x < / 2 r, r = Re. 3

14 References [] L. Adams, C. Bogner, S. Weinzierl, The iterated structure of the all-order result for the two-loop sunrise integral, J. Math. Phys. 57 3) 206) arxiv: , doi:0.063/ [2] L. Adams, C. Bogner, S. Weinzierl, The sunrise integral and elliptic polylogarithms, in: 3th DESY Workshop on Elementary Particle Physics: Loops and Legs in Quantum Field Theory LL206) Leipzig, Germany, April 24-29, 206, 206. arxiv: URL [3] L. Adams, C. Bogner, A. Schweitzer, S. Weinzierl, The kite integral to all orders in terms of elliptic polylogarithmsarxiv: [4] S. Bloch, M. Kerr, P. Vanhove, A Feynman integral via higher normal functions, Compos. Math ) arxiv: , doi:0.2/s000437x [5] S. Bloch, M. Kerr, P. Vanhove, Local mirror symmetry and the sunset Feynman integralarxiv: [6] A. Levin, Elliptic Polylogarithms: an analytic theory, Compositio Math ) ). [7] F. C. S. Brown, A. Levin, Multiple Elliptic Polylogarithms, eprint arxiv: ). [8] A. Beilison, A. Levin, Elliptic Polylogarithms, Proceedings of Symposia in Pure Mathematics, Vol ), Part 2, pp ). [9] S. Bloch, Higher regulators, algebraic K -theory and zeta-functions of elliptic curves, lecture notes, U.C.Irvine 977). [0] D. Zagier, The Bloch-Wigner-Ramakrishnan polylogarithm function, Math.Ann 286, ) 990). [] S. Bloch, P. Vanhove, The elliptic dilogarithm for the sunset graph, J. Number Theor ) arxiv: , doi: 0.06/j.jnt [2] K. S. Kolbig, Nielsen S generalized polylogarithms, SIAM J. Math. Anal ) doi:0.37/ [3] M. Ismail, C. Libis, Contiguous relations, basic hypergeometric functions, and orthogonal polynomials, Journal of Mathematical Analysis and Applications, Volume 4, Issue 2, Pages ). [4] G. Gasper, M. Rahman, Basic Hypergeometric Series, Cambridge University Press 990). [5] E. Heine, Untersuchungen über die Reihe, j. reine angew. math. 34, ). [6] S. Ruijsenaars, A Generalized Hypergeometric Function Satisfying Four Analytic Difference Euations, Comm Math Phys 999) 206: ). [7] T. Oshima, On convergence of basic hypergeometric series 205). arxiv: [8] J. Roues, Galois groups of the basic hypergeometric euations 2007). arxiv: [9] Q.-H. H. William Y. C. Chen, Y.-P. Mu, Non-Terminating Basic Hypergeometric Series and the -Zeilberger Algorithm, Proceedings of the Edinburgh Mathematical Society Series 2), Volume 5, Issue 03, pp ). [20] M. A. Khan, A. H. Khan, A note on mixed hypergeometric series, Acta Mathematica Vietnamica, Vol ). 4