A numerical test of differential equations for one- and two-loop sunrise diagrams using configuration space techniques

Transcription

1 Eur. Phys. J. C 7:85 DOI.4/epjc/s5--85-z Special Article - Tools for Experiment and Theory A numerical test of differential equations for one- and two-loop sunrise diagrams using configuration space techniques S. Groote,,a, J.G. Körner, A.A. Pivovarov 3,4 Füüsika Instituut, Tartu Ülikool, Tähe 4, 5 Tartu, Estonia Institut für Physik, Johannes-Gutenberg-Universität, Staudinger Weg 7, 5599 Mainz, Germany 3 Institute for Nuclear Research, Russian Academy of Sciences, Moscow 73, Russia 4 Department Physik, Universität Siegen, Walter-Flex-Str. 3, 5768 Siegen, Germany Received: June / Revised: 5 July / Published online: 8 July The Authors. This article is published with open access at Springerlink.com Abstract We use configuration space methods to write down one-dimensional integral representations for one- and two-loop sunrise diagrams also called Bessel moments which we use to numerically check on the correctness of the second order differential equations for one- and two-loop sunrise diagrams that have recently been discussed in the literature. Introduction Sunrise-type diagrams have been under investigation since many years. Exact analytical results can be obtained only for special mass or kinematic configurations such as for the equal or zero mass cases or for the threshold region. For example, threshold expansions of the non-degenerate massive two-loop sunrise diagram have been studied in Refs. [, ]. The construction of differential equations for the corresponding correlator function provides some hope that by solving these differential equations, a general analytical solution can be obtained. Recently, mathematical methods were used to construct the coefficients of such a differential equation in a systematic way [3]. This work supplements the work of Kotikov [4] and Remiddi et al. [5 7] onthe same subject. While traditionally the correlator is calculated in momentum space, configuration space techniques allow for a surprisingly simple solution for sunrise-type diagrams: The correlator in configuration space is just a product of single propagators which in turn can be expressed by modified Bessel functions of the second kind. Transforming back to momentum space, one ends up with a one-dimensional a groote@thep.physik.uni-mainz.de integral over Bessel functions, known as Bessel moments [8, 9]. Asoutlinedin a series of papers [, 6], the corresponding one-dimensional integral can be easily integrated numerically for an arbitrary number of propagators with different masses in any space-time dimension. Therefore, configuration space techniques can be used to numerically check the differential equations for the correlator function obtained by other means. This will be detailed in this note. The paper is organized as follows: In Sect. we introduce the configuration space techniques which will be used in Sect. 3 to check the differential equations for one-loop sunrise-type diagrams. In Sect. 4 we check the differential equations for the two-loop sunrise diagrams for the equal mass case, while in Sect. 5 we will deal with nondegenerate cases. Our conclusions can be found in Sect. 6. Even though the configuration space techniques are well suited to treat general D 4 space-time dimensions, we will mainly deal with the case of D space-time dimensions in this paper. For reasons of simplicity, throughout this paper we work in the Euclidean domain. The transition to the Minkowskian domain can be obtained as usual by a Wick rotation or, equivalently, by replacing p p. Configuration space techniques In configuration space, the n-loop n-particle irreducible correlation function Πx T jxj connecting the space-time points and x is given by the product of the propagators, n+ Πx Dx,m i, i

2 Page of 6 Eur. Phys. J. C 7:85 where the free propagator of a particle with mass m in D- dimensional Euclidean space-time is given by Dx,m e ip x d D p π D p + m mxλ K λ mx π λ+ x λ 3 D λ +. K λ z is the McDonald function modified Bessel function of the second kind. Note that p and x in the integral expression of Eq. 3 ared-dimensional Lorentz vectors, i.e. p x p μ x μ, while the quantity x in the rightmost expression of Eq. 3 and, therefore, also in the argument of the propagator denotes the absolute value x x μ x μ. In the limit mx atfixedx, the propagator simplifies to Dx, π D e ip x d D p p Ɣλ 4π λ+, 4 xλ where Ɣλ is Euler s Gamma function. In this note we write the n-particle irreducible correlator function in Euclidean momentum space. The momentum space n-particle irreducible correlator function is given by the Fourier transform of the n-particle irreducible correlator function Πx in configuration space, Πp Πxe ip x d D x. 5 As a product of propagators, Πx in 5 depends only on the absolute value x x μ x μ. Therefore, one proceeds by first integrating the exponential factor over the D dimensional hypersphere. We write d D x x D d D ˆxdx where d D ˆx denotes the D dimensional integration measure over the D dimensional hypersphere. The integration of the exponential factor over the unit sphere gives λ J λ px. 6 px e ip x d D ˆx π λ+ J λ z is the Bessel function of the first kind. As before, p and x on the right hand side of Eq. 6 stand for the absolute values p p μ p μ and x x μ x μ. Therefore, the correlator in momentum space depends only on the absolute value of the momentum, px λ Πp π λ+ J λ pxπxx λ+ dx. 7 This is the central formula for our numerical verification of the correctness of the differential equations. 3 The one-loop case In Ref. [5], Remiddi explains how to obtain the differential equation for the one-loop sunrise-type diagram with arbitrary masses and dimensions. By applying the integrationby-parts technique to the correlator in momentum space, recurrence relations can be obtained. Finally, Euler s theorem for homogeneous functions connects the loose ends of the iterative steps involving partial derivative with respect to p. We have numerically checked all these steps and have found numerical consistency up to Stokes contributions due to surface terms in integer space-times dimensions. To be more precise, the integral d D k π D k μ v μ k + m p k + m for v k,p or a linear combination of both leads to a surface term which can be assumed to vanish up to Stokes contributions. The integration-by-parts technique consists in calculating the integral explicitly and one then expresses the result in terms of scalar integrals Sα,α : 8 d D k π D k + m α p k + m. 9 α The two resulting recurrence relations read DSα,α + α m Sα +,α Sα,α + α p + m + m Sα,α + Sα,α + Sα,α, α p m + m Sα +,α Sα +,α + Sα,α + α p + m m Sα,α + Sα,α + + Sα,α. Equation can be replaced by Eq. with the two lines interchanged, DSα,α + α m Sα,α + Sα,α + α p + m + m Sα +,α Sα +,α Sα,α. It is obvious that Eq. is reproduced as difference of Eq. and Eq.. Therefore, one has to check only Eq.. We will perform this numerical check for the parameter choice α α which is relevant for the differential equation, and for D space-time dimensions. As mentioned in Ref. [3], even though other dimensions are feasible, this choice avoids singular contributions and serves for the simplest integrand. The equation to be checked is m S, + p + m + m S, S, + S,. 3 Starting from d k S, π k + m p k + m J pxk m xk m xdx, 4 π the integrals Sα,α with higher integer values of α i can be obtained as partial derivatives with respect to the masses,

3 Eur. Phys. J. C 7:85 Page 3 of 6 S, p m m p + m m + m m + S,. m S, Because of J z J z, Eq. can be translated to d px [ pxj pxk m xk m x k m π k + m p k + m + J pxm xk m xk m x d k + J pxk m xm xk m x π k + m p k + m S, 5 J pxk m xk m x ] dx x. and accordingly We checked on the latter relation with an even better precision of the order of 3. The differential equation S, S,. 6 m m In addition one has d k S, π k + m 4πm, S, d k π p k + m 4πm. 7 The derivative can be expressed by K z K z. The higher order integrals in the configuration space representation are given by S, S, m m 4πm J pxm xk m xk m xx dx, 8 S, S, m m 4πm J pxk m xm xk m xx dx. Therefore, Eq. 3 in the configuration space representation reads J px [ m x m xk m xk m x + px + m x + m x K m xm xk m x m x K m xk m x ] dx x. 9 We were able to check this equation numerically for different values of p as function of m and m. The 3D-plot in MATHEMATICA shows stochastic fluctuations around the expected value of of the order of 9. Euler s theorem of homogeneous functions leads to the differential equation p + m + m p + m m S, p p + m + m S, in Ref. [5] is obtained by inserting the recurrence relations into Euler s differential equation. Using the configuration space representation, Eq. reads [ px px + m x + m x px + m x m x pxj px + px + m x + m x ] J px K m xk m x dx x. 3 This equation could be checked with a precision of the order of 8. While Euler s differential equation can be derived from general principles also for the configuration space representation, the recurrence relations can be derived only via the momentum space representation. If one does not use this technique, it remains unclear why such integral identities exist for general parameters p, m and m. In order to check whether one can derive further relations by using integral identities in configuration space, we have used the general ansatz J px [ A px, m x,m xk m xk m x + A px, m x,m xm xk m xk m x + A px, m x,m xk m xm xk m x ] dx x, 4 where A i p, m,m a i p +a i m +a im i,,. By choosing random values for p, m and m, and solving the resulting system of equations, one obtains a, a a, a + a, a a, a a, a a, 5 a a, a + a, a a,

4 Page 4 of 6 Eur. Phys. J. C 7:85 where a is an arbitrary parameter. This leads to Eq. 9 and J px [ m x m x K m xk m x + px + m x m x m x K m xk m x px m x + m x K m xm xk m x ] dx x, 6 which is the difference of Eq. 9 and the same equation with m and m interchanged. We conclude that no more recurrence relations can be found that go beyond Eq The two-loop case with equal masses The differential equation for the two-loop degenerate sunrise diagram has been given in Ref. [7]. It reads p p + m p + 9m d dp + 34 Dp Dm p + 9Dm 4 d dp + D 3 D 4p D + 4m S,, 3 D 4 π, 7 which simplifies to p p + m p + 9m d dp + 3p 4 + m p + 9m 4 d dp + p + 3m S,, 3 8π 8 in D space-time dimensions. We write S,, in a form which is easily adapted to the non-degenerate mass case to be discussed later on. One has S,, π Differentiation of Eq. 9 gives d S,, dp d S,, p dp J pxk m xk m xk m 3 xx dx. 9 π p pxj pxk m xk m x K m 3 xx dx 3 and d dp S,, 4π p 4 [ pxj px + px J px J px ] K m xk m xk m 3 xx dx. 3 Returning to the degenerate mass case the differential equation 8 can be translated to 3 [ p 4p + m p + 9m pxj px + px J px J px 3p 4 + m p + 9m 4 pxj px + 4p p + 3m ] J px K mx 3 xdx, 3 where J z J z J z J z/ is used. Because the result contains the second derivative of the Bessel function, one can use Bessel s differential equation z J λ z + zj λ z + z λ J λ z 33 for λ to compactify the result, 3 4p [ 4p p + 3m p + m p + 9m px J px p p + 5m J px ] K mx 3 xdx. 34 These results have been checked with a precision of the order of 8. 5 The two-loop case with arbitrary masses The final section of this paper is devoted to the second order differential equation, derived for the two-loop sunrise diagram with arbitrary masses in Ref. [3]. After adjusting the normalization, the differential equation can be written as [ p p d dp + p p d dp + p p ] S,, p 3 p 4π, 35 where the coefficients p i t i,,, 3 are given by

5 Eur. Phys. J. C 7:85 Page 5 of 6 p t t t m + m + m 3 t m + m + m 3 t m m + m 3 t m + m m 3 p 3 p [ p p p 4 pxj px + px 3t M t M + M, 36 J px J px p t 9t 6 3M t M + 7M t 4 p p 8M M + 44M t 3 p pxj px + 4p p ] J px K m xk m xk m 3 xx dx 3M 4 36M M 4M t [ 8M 5 + 4M 4 6M 3 96M 4p 3 p p p p 4 px J px + 44M t M 6 6M 5 + 5M ] 4 p p M M 4 M 3 6M, 37 p pxj px p t 3t 5 7M t 4 + M + 6M t 3 + 6M 3 4M t 5M 4 8M 3 + 6M 8M t + M 5 3M 4 + M 3 + 8M 3 M, 38 p 3 t 8t 4 + 4M t 3 + 4M 4M t + 8M 3 + 8M + 48M t + M 4 + 8M 3 M 8M + ct,m,m,m 3 ln m /μ K m xk m xk m 3 xx dx, 4 where p t p t + p t/t. Using different values for p and m 3, in terms of m and m we obtain a 3D-plot with MATHEMATICA which shows again stochastic fluctuations of the order of 4. In the course of our numerical checks we were able to identify two typos in the coefficients of ct,m,m,m 3 in the preprint version of Ref. [3] which we have corrected. 6 Conclusions + ct,m,m 3,m ln m /μ + ct,m 3,m,m ln m 3 /μ, 39 and where M λ λ λ 3 m σλ m σλ m σλ3 3 4 σ are monomial symmetric polynomials in m, m and m 3 and where ct,m,m,m 3 m + m + 3 m t 3 + 6m 4 3m4 3m 4 3 7m m 7m m 3 + 4m m 3 t + 6m 6 + 3m6 + 3m6 3 + m4 m + m4 m 3 8m m4 8m m4 3 3m4 m 3 3m m4 3 t + m 8 m8 m8 3 5m6 m 5m6 m 3 + m m6 + m m m6 m 3 + 4m m m 4 m4 + 3m4 m4 3 6m4 m4 3 + m4 m m 3 m m4 m 3 m m m4 3 4 for details, cf. Ref. [3]. Using Eqs. 9, 3 and 3, one obtains Using configuration space techniques, we were able to check numerically the differential equations for sunrise-type diagrams found in the literature. The precision of our numerical test is still quite moderate, but gives sufficient confidence in the validity of the differential equations derived by other means. For example, the introduction of artificial typos in the coefficients of the differential equations are easily discovered. More rigorous tests would require the use of more stable integration routines than those provided by MATHE- MATICA. For the future we hope to find independent routes to discover further relations between Bessel moments which may lead to generalizations of the present findings to cases involving three-loop or even higher order sunrise-type diagrams. Acknowledgements We want to thank Stefan Weinzierl and Anatoly Kotikov for helpful discussions and Volodya Smirnov for encouragement. S.G. acknowledges the support by the Estonian target financed Project No. 856s9, by the Estonian Science Foundation under grant No. 8769, and by the Deutsche Forschungsgemeinschaft DFG under No. KO 69/4-. A.A.P. acknowledges partial support by the RFFI grant No. --8-a. Open Access This article is distributed under the terms of the Creative Commons Attribution Noncommercial License which permits any noncommercial use, distribution, and reproduction in any medium, provided the original authors and source are credited.

6 Page 6 of 6 Eur. Phys. J. C 7:85 References. A.I. Davydychev, V.A. Smirnov, Nucl. Phys. B 554, S. Groote, A.A. Pivovarov, Nucl. Phys. B 58, S. Müller-Stach, S. Weinzierl, R. Zayadeh, Commun. Number Theory Phys. 6, 4. A.V. Kotikov, Phys. Lett. B 54, E. Remiddi, Nuovo Cimento A, M. Caffo, H. Czyż, S. Laporta, E. Remiddi, Nuovo Cimento A, S. Laporta, E. Remiddi, Nucl. Phys. B 74, D.H. Bailey, J.M. Borwein, D. Broadhurst, M.L. Glasser, arxiv:8.89 [hep-th] 9. D. Broadhurst, arxiv:8.483 [hep-th]. E. Mendels, Nuovo Cimento A 45, S. Groote, J.G. Körner, A.A. Pivovarov, Phys. Lett. B 443, S. Groote, J.G. Körner, A.A. Pivovarov, Nucl. Phys. B 54, S. Groote, J.G. Körner, A.A. Pivovarov, Eur. Phys. J. C, R. Delbourgo, M.L. Roberts, J. Phys. A 36, S. Groote, J.G. Körner, A.A. Pivovarov, Eur. Phys. J. C 36, S. Groote, J.G. Körner, A.A. Pivovarov, Ann. Phys. 3, 374 7