difference equations

We describe a new method of calculation of generic multiloop master integrals based on the numerical solution of systems of difference equations in one variable. We show algorithms for the construction of the systems using integration-by-parts identities and methods of solutions by means of expansions in factorial series and Laplace transformation. We also describe new algorithms for the identification of master integrals and the reduction of generic Feynman integrals to master integrals, and procedures for generating and solving systems of differential equations in masses and momenta for master integrals. We apply our method to the calculation of the master integrals of massive vacuum and self-energy diagrams up to three loops and of massive vertex and box diagrams up to two loops. Implementation in a computer program of our approach is described. Important features of the implementation are: the ability to deal with hundreds of master integrals and the ability to obtain very high precision results expanded at will in the number of dimensions.

[1]  Irene A. Stegun,et al.  Handbook of Mathematical Functions. , 1966 .

[2]  R. Stuart,et al.  On the precise determination of the Fermi coupling constant from the muon lifetime , 1999, hep-ph/9904240.

[3]  Deep-inelastic structure functions at two loops , 1999, hep-ph/9912355.

[4]  G. Jikia,et al.  Analytic evaluation of two-loop renormalization constants of enhanced electroweak strength in the Higgs sector of the Standard Model , 1996, hep-ph/9609447.

[5]  F. Tkachov,et al.  Integration by parts: The algorithm to calculate β-functions in 4 loops , 1981 .

[6]  J. B. Tausk,et al.  Two-loop self-energy diagrams with different masses and the momentum expansion , 1993 .

[7]  P. Baikov,et al.  Equivalence of recurrence relations for Feynman integrals with the same total number of external and loop momenta , 2000, hep-ph/0001192.

[8]  F. Tkachov.  A theorem on analytical calculability of 4-loop renormalization group functions , 1981 .

[9]  E. Remiddi,et al.  The analytical value of the electron (g − 2) at order α3 in QED , 1996 .

[10]  Differential Equations for Two-Loop Four-Point Functions , 1999, hep-ph/9912329.

[11]  E. N.,et al.  The Calculus of Finite Differences , 1934, Nature.

[12]  K. Knopp.  Theory and Application of Infinite Series , 1990 .

[13]  The analytical value of the electron (g-2) at order alpha^3 in QED , 1996, hep-ph/9602417.

[14]  P. Baikov.  Explicit solutions of the 3-loop vacuum integral recurrence relations , 1996 .

[15]  R. Roskies,et al.  Hyperspherical approach to quantum electrodynamics: sixth-order magnetic moment , 1974 .

[16]  J. Fleischer.  Erratum to: Two-loop self-energy master integrals on shell , 1999 .

[17]  D. J. Broadhurst.  Massive 3-loop Feynman diagrams reducible to SC$^*$ primitives of algebras of the sixth root of unity , 1999 .

[18]  M. Steinhauser,et al.  Automatic Computation of Feynman Diagrams , 1999 .

[19]  A. Kotikov.  Differential equations method. New technique for massive Feynman diagram calculation , 1991 .