arxiv:hep-ph/ v1 28 Mar 2006

Transcription

1 DESY-06-0 SFB/CPP-06- Hypergeometric representation of the two-loop equal mass sunrise diagram arxiv:hep-ph/0607v 8 Mar 006 O.V. Tarasov Deutsches Elektronen - Synchrotron DESY Platanenallee 6, D-578 Zeuthen, Germany Oleg.Tarasov@desy.de Abstract A recurrence relation between equal mass two-loop sunrise diagrams differing in dimensionality by is derived and it s solution in terms of Gauss F and Appell s F hypergeometric functions is presented. For arbitrary space-time dimension d the imaginary part of the diagram on the cut is found to be the F hypergeometric function with argument proportional to the maximum of the Kibble cubic form. The analytic expression for the threshold value of the diagram in terms of the hypergeometric function F of argument / is given. On leave of absence from JINR, 4980 Dubna (Moscow Region), Russian Federation.

2 Introduction The evaluation of radiative corrections for modern high precision particle physics is becoming a more and more demanding task. Without inventing new mathematical methods and new computer algorithms the progress in calculating multi-loop, multi-leg Feynman diagrams depending on several momentum and mass scales will be not possible. An important class of radiative corrections comes from self-energy type of Feynman diagrams, which also occur in evaluating vertex, box and higher multileg diagrams. At the two-loop level different approaches for calculating self-energy diagrams are available []. A general algorithm for reduction of propagator type of diagrams to a minimal set of master integrals was proposed in []. This recurrence relations algorithm has been implemented in computer packages in [] and [4],[5]. At present the most advanced package available for calculating two-loop self-energy diagrams with arbitrary massive particles was written by Stephen Martin and David Robertson [5]. It includes procedures for numerical evaluation of master integrals with arbitrary masses and also a database of analytically known master integrals. Integral representation for master integrals with arbitrary masses in four dimensional space time was proposed in [6]. Despite intensive efforts by many authors not all two-loop self-energy integrals with a mass are known analytically. Even the imaginary part of the simplest sunrise self-energy diagram with three equal mass propagators was not known for arbitrary space-time dimension d until now. The two-loop sunrise integral with three equal masses was investigated in many publications [7]- [5]. Small and large momentum expansions of this integral for arbitrary space-time dimension d can be found in []. It s threshold expansion was given in [4]. A numerical procedure for evaluating the sunrise integral was described in []. The very latest effort of an analytic calculation of the diagram, by using the differential equation approach [6], was undertaken in Ref. [7]. It is the purpose of this paper to describe a new method and to present the analytic result for the equal mass two-loop sunrise master integral. To accomplish our goal we use the method of evaluation of master integrals by dimensional recurrences proposed in [8]. The application of this method to one-loop integrals was presented in [9] and [0]. In the present publication we extend the method to two-loop integrals. As was already discovered in the one-loop case, the solutions of dimensional recurrences are combinations of hypergeometric functions. The knowledge of the hypergeometric representation of an integral means that we possess the most complete mathematical information available. This information can be effectively used in several respects. First, through analytic continuation formulae, the hypergeometric functions valid in one kinematic domain can be re-expressed in a different kinematic region. Second, these hypergeometric functions often have integral representation themselves, in which an expansion in ε = (4 d)/ can be made, yielding expressions in logarithms, dilogarithms, elliptic integrals, etc.. Since very similar hypergeometric functions come from different kind of Feynman integrals the ε expansion derived in solving one problem can be used in other applications. Essential progress in the ε expansion of hypergeometric functions encountered in evaluating Feynman diagrams was achieved in [],[], [] and []. Third, because the hypergeometric series is convergent and well behaved in a particular region of kinematical variables, it can be numerically evaluated [], [4]. In addition a hypergeometric representation allows an asymptotic expansion of the integral in terms of ratios of different Gram determinants or ratios of momentum and mass scales which can provide fast numerical convergence of the result. Our paper is organized as follows. In Sec., we present the relevant difference equation connecting sunrise integrals with dimensionality differing by as well as the differential equation for this integral. In Sec.. the method for finding the full solution of the dimensional recurrency is elaborated. Explicit expression for the sunrise integral in terms of Appell s function F and Gauss hypergeometric function F is constructed in Sec. 4 and in Sec. 5 the differential equation approach and the method of dimensional recurrences is compared. In the Appendix some useful formulae for the hypergeometric

3 functions F, F and F are given together with their integral representations. Difference and differential equations for the sunrise integral The generic two-loop self-energy type diagram in d dimensional Minkowski space with three equal mass propagators is given by the following integral: d d (ν k d d k,ν,ν ) (iπ d/ ) (k m ) ν ((k k ) m ) ν ((k q) m ) ν. (.) For integer values of ν j the integrals (.) can be expressed in terms of only three basis integrals (d) (,,), J(d) (,,) and J(d) (0,,) = (T (m )) where T (d) d (m d k ) = ( [iπ d ] k m = Γ d ) m d. (.) The relation connecting d and d dimensional integrals (ν,ν,ν ) follows from the relationship given in Ref. []: J (d ) (ν,ν,ν ) = ν ν (ν +,ν +,ν ) + ν ν (ν +,ν,ν + ) + ν ν (ν,ν +,ν + ). (.) Relation (.) taken at ν = ν = ν = and ν =, ν = ν = gives two equations. To simplify these equations we use the recurrence relations proposed in []. From these two equations by shifting d d + two more relations follow. They are used to exclude (,,) from one of the relations, so that we obtain a difference equation for the master integral where z (d + )(d )(d + 4)(d + ) 4m 4 ( z)( 4z + 9z )z(d )d 4m 8 ( z) ( 9z) (,,) J(d) : = z[(z + )(7z + 8z )d 4z( + 9z)d 48z ]m d+ Γ J (d+4) J (d+) ( d ), (.4) z = m q. (.5) The integral satisfies also a second order differential equation []. Taking the second derivative of with respect to mass gives d dm dm J(d) (,,) = 6J(d) (,,) + 6J(d) (,,). (.6) Again using the recurrence relations from [], the integrals on the r.h.s can be reduced to the same three basis integrals. Using (,,) = d dm J(d) (,,) (.7) from (.6) we obtain: ( z)( 9z)z d dz z[9z (d 4) + 0z(d ) + 8 d] dj(d) dz + (d )[z(d + 4) + d 4] = zm (d 6) Γ ( d ). (.8) The differential equation (.8) will be used in Sec.4 to find the momentum dependence of arbitrary periodic constants in the solution of the difference equation (.4).

4 Solution of the dimensional recurrency Equation (.4) is a second order inhomogeneous equation with polynomial coefficients in d. The full solution of this equation is given by (see Ref. [5] and references therein): where p = p + w a(d) a + w b(d) b, (.9) b is a fundamental system of solutions of the associated homogeneous equation and w a (d), w b (d) are arbitrary periodic functions of d satisfying relations: w a (d + ) = w a (d), w b (d + ) = w b (d). (.0) is a particular solution of (.4), a,j(d) The order of the polynomials in d of the associated homogeneous difference equation can be reduced by making the substitution ( ) Γ d = ( ) ( ) Γ d Γ d. (.) The associated homogeneous equation for takes the simpler form Putting 6z 7m 8 ( z) ( 9z) J(d+4) d( z)( 4z + 9z )z 7m 4 ( z) ( 9z) d = k ε, we transform Eq.(.) to a standard form J (d+) (d )(d 4) 6 = 0. (.) y (k) = ρ k J (k ε), (.) Aρ y (k+) + (B + C k)ρy (k+) (α + k)(β + k)y (k) = 0, (.4) where A = 6z 7m 8 ( z) ( 9z), B = 4ε ( z)( 4z + 9z )z 7 m 4 ( z) ( 9z), C = B ε, α = ε, β = ε, (.5) and ρ is for the time being, an arbitrary constant. In order to get Eq.(.4) into a more convenient form, we will define three parameters ρ, x and γ by the equations Aρ = x( x), Bρ = γ (α + β + )x, Cρ = x. (.6) These have the solution x = Cρ = ( 9z) ( + z) = q (q 9m ) (q + m ), (.7) ρ = = 7 m 4 ( z) ( 9z) 4A + C 4 z( + z) = 7 4 m (q m ) (q 9m ) (q + m ), (.8) γ = Bρ + (α + β + )x = ε, (.9) and Eq. (.4) can accordingly be written in the form x( x)y (k+) + [( x)k + γ (α + β + )x]y (k+) (α + k)(β + k)y (k) = 0. (.0)

5 The fundamental system of solutions of this equation consists of two hypergeometric functions [5]. For example, in the case when x < (large q ) the solutions are y (k) = ( ) k Γ(α + k)γ(β + k) Γ(α + β γ + k + ) F (α + k,β + k,α + β γ + k + ; x), y (k) = Γ(α + β γ + k) ( x) k F (γ α,γ β,γ α β + k; x). (.) Once we know the solutions of the homogeneous equation, a particular solution p can be obtained by using Lagrange s method of variation of parameters. Lagrange method for finding a particular solution is well described in [5]. The application of the method is straightforward but tedious. Explicit result will be given in the next section. It is interesting to note that the argument of the Gauss hypergeometric function is related to the maximum of the Kibble cubic form [6]: Φ(s,t,u) = stu (s + t + u)m (m + q ) + m 4 (m + q ), (.) provided that the following condition is satisfied: s + t + u = q + m. (.) The maximal value Φ max = 7 q (q 9m ) occurs at s = t = u = (q + m ) and we see that the kinematical variable (.7 ) can be written as x = Φ(s,t,u) stu s=t=u= (q +m ). (.4) This observation may be useful in finding the characteristic variable in the general mass case [7]. Also one can try to apply the method described above to find the imaginary part of the sunrise integral in the general mass case in arbitrary space-time dimension. 4 Explicit analytic expression for To find the full solution of Eq. (.4) we assume that q is large. The region of large momentum squared corresponds to x and therefore as a fundamental system of solutions of the homogeneous equation we take y (k) and y (k). According to (.), (.) and (.) the solution of the associated homogeneous difference equation will be of the form ( ( ) ( ),h = w (z) Γ d ) Γ d Γ d [ d ( ) ( ) ( ) Γ d Γ d ρ d e iπ d Γ d F, d ; ] d ; x ( ) Γ d [ ρ d + w (z) ( ) ( ) Γ d Γ d F, ; ] ( x) d d ; x. (4.5) The arbitrary periodic functions w (z) and w (z) can be determined either from the d asymptotics or using the differential equation (.8). Substituting (4.5) into (.8) we obtain two simple equations z( z)( + z)( 9z) dw (z) ( + 6z 9z )w (z) = 0, dz z( + z)( 9z) dw (z) + ( z)w (z) = 0. (4.6) dz 4

6 Both equations are independent of d and their solutions w (z) = κ z ( + z) ( 9z) ( z), w (z) = κ z ( + z)( 9z), (4.7) determine the periodic functions up to integration constants κ, κ which we fix from the first two terms of the large momentum expansion of presented in []: [ ] Γ( + ε)γ = m 4ε Γ ( ε) ( + ε) zγ ( + ε)γ( ε) ( z)ε + 6Γ ( ε) Γ( ε) ( z)ε + O(z). (4.8) The application of Lagrange s method of finding a particular solution gives ( p = zmd 6 ( + z) Γ d ) ( F,, d, d,d ; ) zr,r, (4.9) where F is the Appell function [8] defined as R = 4 z ( + z), (4.0) F (α,β,β,γ,γ ;x,y) = k,l=0 (α) k+l (β) k (β ) l (γ) k (γ ) l x k y l, x + y <, (4.) k! l! and (a) n = Γ(a + n)/γ(a) denotes the Pochhammer symbol. Collecting all contributions, setting d = 4 ε, applying Euler s transformation for the first F function in (4.5) we obtain the following solution of the difference equation (.4): = 6Γ ( ε)γ ( + ε)( z) ε ( z) ε [ m 4ε F, ; Γ( ε)( + z) ε; + Γ( + ε)γ ( ε)( z) ε ( 9z) ε m 4ε Γ( ε)z( + z) + zm 4ε ( + z) Γ ( + ε) F F [, ; ε; 7( z) ] z ( + z) (,, ε, ε, ε; zr,r 7( z) ] z ( + z) ). (4.) This is our main result. The imaginary part of on the cut comes from the first two terms in (4.). The analytic continuation of the two F functions gives the relatively simple expression: Im = 4z π πm 4ε [ ( 9z) ( ) Γ ε Γ ( ε) ( + z) 08z ] ε F [, ; ε; ( 9z) ] ( + z). (4.) At d = 4 for the imaginary part this verifies the result of []. Expanding at q = 9m we reproduced the singular and finite in ε parts of several terms of the on-threshold asymptotic expansion presented in Ref. [4]. We found the following integral representation for Appell s F function in (4.) F (,, ε, ε, ε; ) zr,r = Γ( ε) Γ ( ε )( + z) 0 dt [t( t)] ε [,ε; (4zt + z + L) F ε; 5 ] 4tz + z L, (4.4) 4tz + z + L

7 where L = (4zt z) 4z. (4.5) This integral representation can be used for the ε expansion of the F function. Integral representation for the F functions convenient for ε expansion is given in the Appendix. The analytic continuation of valid near the singular points q = 0,m,9m can be directly obtained from (4.) by performing the analytic continuation of the hypergeometric functions. Explicit formulae of the analytic continuations of in terms of Olsson functions [9] as well as the ε expansion of the result we are planning to present in a separate publication [0]. Using (4.) we can find the on-threshold value of the integral. In formula (4.) at q = 9m the imaginary part of the first term cancels the imaginary part of the second term. The real part coming from the F terms cancels the term with F which comes from the Appell function F at the threshold ( F, ε,, ε, ε; 4, ) = 4 [ 6( ε), + ε, + ( ε) F ε; + ε, ε; ] ( ) + 4Γ( ε)γ( ε)γ + ε [ Γ( ε) ( ) ( ) ε Γ ε Γ 5 ε F + ε, ε; 5 ε; ]. (4.6) The cancellation of the F functions happens due to the fact that [ ε, F + ε; 5 ε; ] ( ) ( ) = 5 ε Γ 5 4 ε Γ 7 4 ε ( ) ( ). ε Γ 4 ε Γ 5 ε (4.7) Adding contributions from different hypergeometric functions gives a rather simple expression for the diagram at q = 9m = q =9m = Γ [ (ε), + ε, ( ε)( ε) F ε; Γ ( + ε) ( ε)( ε) + ε, ε; { ε + 9 4ε π + O(ε) The first several terms in the ε expansion are in agreement with the result of Ref.[4]. 5 Conclusions ] }. (4.8) Here we would like to add several remarks, which underline the most important points which follow from the results of this paper. Using a new method, for the first time, we were able to obtain an analytic expression for the two-loop sunrise, self-energy diagram with equal mass propagators. Until now all attempts to find such a result with other methods failed for this integral. This clearly demonstrates that the method of dimensional recurrences is a powerful tool for calculating Feynman integrals. In our opinion there is a deep reason why, for example, with the differential equation method [6] used in [7] an explicit formula could not be found. It turns out that even the associated homogeneous differential equation is rather complicated. Relevant for this case is the Heun equation [] with four regular singular points, located at q = 0,m,9m,. In general the reduction of the Heun equation to the hypergeometric equation is a complicated mathematical problem [] which is not completely solved until now. 6

8 At the same time the associated homogeneous difference equation for is rather simple, and admits a reduction to a hypergeometric type of equation with linear coefficients. In fact this is a rather general situation. Kinematical singularities of Feynman integrals are located on complicated manifolds. In the case when the differential equations are of the first order there are no problems to solve them. However, to solve a second or higher order differential equations in general will be a problem because of complicated structure of the kinematical singularities. The location of the singularities of Feynman integrals with respect to the space time dimension d is well known. This has been used for a rather evident rescaling of the integral by ratios of Γ functions which allowed us to significantly reduce the order of the polynomial coefficients in the difference equation as we have seen in Sec.. Finally this simplification allowed us to obtain the explicit result. We expect that a further development of the method we used in the present paper will help to find analytic results for other more complicated Feynman integrals. 6 Appendix The series representation for the Appell function F : F (α,β,β,γ,γ ;x,y) = k=0 (α) k (β) k (γ) k k! x k F (α + k,β ;γ ;y) (6.9) is convenient for the analytic continuations and also for the evaluation of F at some particular values of their arguments. The most frequently used integral representations for F is F (α,β,β,γ,γ ;w,z) = 0 0 Γ(γ)Γ(γ ) Γ(β)Γ(β )Γ(γ β)γ(γ β ) du dv u β v β ( u) γ β ( v) γ β ( uw vz) α. (6.40) For the special parameters which appear in the explicit result (4.) we have: ( F, ε, ), ε, ε;x,y = Γ( ε)γ( ε) ( ) ( ) Γ ε Γ 0 0 dudv [u( u)( v)] ε. (6.4) v ( uw vz) Furthermore, we found the following specific relation between the F function and yet another Appell s function F [ F, ; 7( z) ] ( ) z ( + z) ε; ( + z) = ( z) F, + ε, + ε, ε; 4z ( + 4z z), ( z). (6.4) To our knowledge such a relation has not been found sofar in the mathematical literature. The integral representation for the F function reads ( ) F, + ε, + ε, ε;w,z Γ( ε) du = ( ) ( ) Γ Γ ε [( u)( wu)( zu)] ε (6.4) 0 u and is convenient for the ε expansion. This F function can be considered as a generating function of a new generalization of elliptic integrals which may appear in evaluating Feynman integrals. Acknowledgments. I am very thankful to Fred Jegerlehner for useful remarks and carefully reading the manuscript. This work was supported by DFG Sonderforschungsbereich Transregio

9 References [] G. Passarino and S. Uccirati, Nucl.Phys. B 69 (00) 97; G. Weiglein, R. Mertig, R. Scharf and M. Böhm, in New computing techniques in physics research, ed. D. Perret-Gallix (World Scientific, Singapore, 99) p.67. [] O. V. Tarasov, Nucl. Phys. B 50 (997) 455. [] R. Mertig and R. Scharf, Comput. Phys. Commun. (998) 65. [4] S. P. Martin, Phys. Rev. D 68 (00) [5] S. P. Martin and D. G. Robertson, Comput. Phys. Commun. 74 (006). [6] S. Bauberger and M. Bohm, Nucl. Phys. B 445 (995) 5. [7] E. Mendels, Nuovo Cim. A45, 87 (978). [8] D. J. Broadhurst, Z. Phys. C 47, 5 (990). [9] C. Ford, I. Jack and D. R. T. Jones, Nucl. Phys. B 87, 7 (99) [Erratum-ibid. B 504, 55 (997)]. [0] A. I. Davydychev and J. B. Tausk, Nucl. Phys. B 97, (99). [] D. J. Broadhurst, J. Fleischer and O. V. Tarasov, Z. Phys. C 60 (99) 87. [] F. A. Berends, M. Buza, M. Bohm and R. Scharf. Z. Phys. C 6 (994) 7, [] J. Fleischer, F. Jegerlehner, O. V. Tarasov and O. L. Veretin, Nucl. Phys. B 59 (999) 67 [Erratum-ibid. B 57 (000) 5]. [4] A. I. Davydychev and V. A. Smirnov, Nucl. Phys. B 554 (999) 9. [5] B. A. Kniehl, A. V. Kotikov, A. Onishchenko and O. Veretin, Nucl. Phys. B 78 (006) 06. [6] A.V. Kotikov, Phys. Lett. B 54 (99) 58; A.V. Kotikov, Mod. Phys. Lett. A 6 (99) 677. [7] S. Laporta and E. Remiddi, Nucl. Phys. B 704 (005) 49. [8] O. V. Tarasov, Phys. Rev. D 54 (996) [9] O. V. Tarasov, Nucl. Phys. Proc. Suppl. 89 (000) 7. [0] J. Fleischer, F. Jegerlehner and O. V. Tarasov, Nucl. Phys. B 67 (00) 0. [] A.I. Davydychev, M.Yu. Kalmykov, Nucl. Phys. B (Proc. Suppl.) 89 (000) 8; Nucl. Phys. B605 (00) 66; A.I. Davydychev, hep-th/99080; Phys. Rev. D6 (000) 08770; A.I. Davydychev, M.Yu. Kalmykov, Proc. Workshop CPP00, Tokyo, Japan, November 00, KEK Proceedings 00-, p. 69 (hep-th/00); S. Moch, P. Uwer, S. Weinzierl, J. Math. Phys. 4 (00) 6; F. Jegerlehner, M.Yu. Kalmykov, O. Veretin, Nucl. Phys. B658 (00) 49; S. Weinzierl, J. Math. Phys. 45 (004) 656; A.I. Davydychev, M.Yu. Kalmykov, Nucl. Phys. B699 (004) ; M.Yu. Kalmykov, Nucl. Phys. B (Proc. Suppl.) 5 (004) 80. 8

10 [] M. Y. Kalmykov, arxiv:hep-th/ [] S. Bauberger, F. Berends, M. Böhm and M. Buza, Nucl. Phys. B44, 8 (995). [4] M.Yu. Kalmykov, A. Sheplyakov, Comput. Phys. Commun. 7 (005) 45. [5] L. M. Milne-Thomson, The Calculus of Finite Differences, Chelsea publishing company, New- York, 98. [6] T.W.B. Kibble, Phys. Rev. 7, 59 (960). [7] A. I. Davydychev and R. Delbourgo, J. Phys. A 7 (004) 487. [8] P. Appell and J. Kampé de Fériet Fonctions hypergeometriques et hyperspériques. Paris: Gauthier Villars (96). [9] Olsson P. O. M. A hypergeometric function of two variables of importance in perturbation theory. II, Arkiv för Fysik, 964. Band 9, No. 8. P ; Olsson P. O. M. Hypergeometric function of two variables of importance in perturbation theory. I, Arkiv för Fysik, 965. Band 0, No. 4. P [0] O. V. Tarasov, in preparation. [] A. Erdèly et. al.,higher Transcendental Functions (McGraw-Hill, New York, 95), Vol., p.57. [] R.S. Maier, On Reducing the Heun Equation to the Hypergeometric Equation, Journal of Differential Equations, v. (005) 7. 9