FIESTA 2: Parallelizeable multiloop numerical calculations

Abstract The program FIESTA has been completely rewritten. Now it can be used not only as a tool to evaluate Feynman integrals numerically, but also to expand Feynman integrals automatically in limits of momenta and masses with the use of sector decompositions and Mellin–Barnes representations. Other important improvements to the code are complete parallelization (even to multiple computers), high-precision arithmetics (allowing to calculate integrals which were undoable before), new integrators, Speer sectors as a strategy, the possibility to evaluate more general parametric integrals. Program summary Program title: FIESTA 2 Catalogue identifier: AECP_v2_0 Program summary URL: http://cpc.cs.qub.ac.uk/summaries/AECP_v2_0.html Program obtainable from: CPC Program Library, Queen's University, Belfast, N. Ireland Licensing provisions: GNU GPL version 2 No. of lines in distributed program, including test data, etc.: 39 783 No. of bytes in distributed program, including test data, etc.: 6 154 515 Distribution format: tar.gz Programming language: Wolfram Mathematica 6.0 (or higher) and C Computer: From a desktop PC to a supercomputer Operating system: Unix, Linux, Windows, Mac OS X Has the code been vectorised or parallelized?: Yes, the code has been parallelized for use on multi-kernel computers as well as clusters via Mathlink over the TCP/IP protocol. The program can work successfully with a single processor, however, it is ready to work in a parallel environment and the use of multi-kernel processor and multi-processor computers significantly speeds up the calculation; on clusters the calculation speed can be improved even further. RAM: Depends on the complexity of the problem Classification: 4.4, 4.12, 5, 6.5 Catalogue identifier of previous version: AECP_v1_0 Journal reference of previous version: Comput. Phys. Comm. 180 (2009) 735 External routines: QLink [1], Cuba library [2], MPFR [3] Does the new version supersede the previous version?: Yes Nature of problem: The sector decomposition approach to evaluating Feynman integrals falls apart into the sector decomposition itself, where one has to minimize the number of sectors; the pole resolution and epsilon expansion; and the numerical integration of the resulting expression. Solution method: The sector decomposition is based on a new strategy as well as on classical strategies such as Speer sectors. The sector decomposition, pole resolution and epsilon-expansion are performed in Wolfram Mathematica 6.0 or, preferably, 7.0 (enabling parallelization) [4]. The data is stored on hard disk via a special program, QLink [1]. The expression for integration is passed to the C-part of the code, that parses the string and performs the integration by one of the algorithms in the Cuba library package [2]. This part of the evaluation is perfectly parallelized on multi-kernel computers. Reasons for new version: 1. The first version of FIESTA had problems related to numerical instability, so for some classes of integrals it could not produce a result. 2. The sector decomposition method can be applied not only for integral calculation. Summary of revisions: 1. New integrator library is used. 2. New methods to deal with numerical instability (MPFR library). 3. Parallelization in Mathematica . 4. Parallelization on multiple computers via TCP-IP. 5. New sector decomposition strategy (Speer sectors). 6. Possibility of using FIESTA to for integral expansion. 7. Possibility of using FIESTA to discover poles in d. 8. New negative terms resolution strategies. Restrictions: The complexity of the problem is mostly restricted by CPU time required to perform the evaluation of the integral Running time: Depends on the complexity of the problem References: [1] http://qlink08.sourceforge.net , open source. [2] http://www.feynarts.de/cuba/ , open source. [3] http://www.mpfr.org/ , open source. [4] http://www.wolfram.com/products/mathematica/index.html .

[1]  M. Steinhauser,et al.  Three-loop static potential. , 2009, Physical review letters.

[2]  Next-to-leading mass singularities in two-loop electroweak singlet form factors , 2004, hep-ph/0401087.

[3]  M. Bergere,et al.  A theorem on asymptotic expansion of Feynman amplitudes , 1978 .

[4]  T. Binoth,et al.  Numerical evaluation of multi-loop integrals by sector decomposition , 2004 .

[5]  T. Huber,et al.  Nine-propagator master integrals for massless three-loop form factors , 2009, 0902.3512.

[6]  Two-loop electroweak next-to-leading logarithmic corrections to massless fermionic processes , 2006, hep-ph/0608326.

[7]  A. Smirnov,et al.  On the resolution of singularities of multiple Mellin–Barnes integrals , 2009, 0901.0386.

[8]  Li Lin Yang,et al.  Two-loop divergences of massive scattering amplitudes in non-abelian gauge theories , 2009, 0908.3676.

[9]  Richard J. Enbody,et al.  Using AMDAHL’s Law as a Metric to Drive Code Parallelization: Two Case Studies , 2001, Int. J. High Perform. Comput. Appl..

[10]  G. Hooft,et al.  Regularization and Renormalization of Gauge Fields , 1972 .

[11]  A. Denner,et al.  High-Energy Approximation of One-Loop Feynman Integrals , 1996 .

[12]  Christian Bogner,et al.  Operating system: Unix , 1983 .

[13]  V. Velizhanin.  Leading transcendental contribution to the four-loop universal anomalous dimension in N = 4 SYM , 2009 .

[14]  V. A. Smirnov,et al.  Hepp and Speer sectors within modern strategies of sector decomposition , 2008, 0812.4700.

[15]  M. Bergere,et al.  Bogolubov–Parasiuk theorem in the α‐parametric representation , 1976 .

[16]  A. V. Smirnov,et al.  Feynman Integral Evaluation by a Sector decomposiTion Approach (FIESTA) , 2008, Comput. Phys. Commun..

[17]  T. Binoth,et al.  Numerical evaluation of phase space integrals by sector decomposition , 2004 .

[18]  M. Tentyukov,et al.  Four-loop massless propagators: A numerical evaluation of all master integrals , 2010, 1004.1149.

[19]  G. Bell.  NNLO corrections to inclusive semileptonic B decays in the shape-function region , 2008, 0810.5695.

[20]  P. Breitenlohner,et al.  Dimensional renormalization and the action principle , 1977 .

[21]  The International Journal of High Performance Computing Applications— , 1998 .

[22]  M. Steinhauser,et al.  Fermionic contributions to the three-loop static potential , 2008, 0809.1927.

[23]  Vladimir A. Smirnov,et al.  Feynman Integral Calculus , 2006 .

[24]  V. A. Smirnov.  Analytical result for dimensionally regularized massless on shell double box , 1999 .

[26]  K. Chetyrkin,et al.  Four Loop Massless Propagators: an Algebraic Evaluation of All Master Integrals , 2010, 1004.1153.

[27]  Takahiro Ueda,et al.  A geometric method of sector decomposition , 2009, Comput. Phys. Commun..

[28]  O. I. Zavialov.  Renormalized Quantum Field Theory , 1990 .

[29]  A. Denner,et al.  Two-loop electroweak next-to-leading logarithms for processes involving heavy quarks , 2008, 0809.0800.

[30]  V. Smirnov,et al.  Three-loop on-shell Feynman integrals with two masses , 2009, 0903.4760.

[31]  V. Smirnov.  Asymptotic expansions in limits of large momenta and masses , 1990 .

[32]  G. Sterman,et al.  Threshold Resummation for Top-Pair Hadroproduction to Next-to-Next-to-Leading Log , 2009, 0907.1790.

[33]  J. Zuber,et al.  Renormalization of Feynman amplitudes and parametric integral representation , 1974 .

[34]  R. Bonciani,et al.  Two-loop QCD corrections to the heavy-to-light quark decay , 2008, 0809.4687.

[35]  Christian Bogner,et al.  Blowing up Feynman integrals , 2008, 0806.4307.

[36]  V. Smirnov,et al.  The regional strategy in the asymptotic expansion of two-loop vertex feynman diagrams , 1998, hep-ph/9812529.

[37]  V. Smirnov.  Evaluating Feynman Integrals , 2005 .

[38]  Volker Pilipp,et al.  Semi-numerical power expansion of Feynman integrals , 2008, 0808.2555.

[39]  Michal Czakon,et al.  Automatized analytic continuation of Mellin-Barnes integrals , 2005, Computer Physics Communications.

[40]  C. Duhr,et al.  An analytic result for the two-loop hexagon Wilson loop in $ \mathcal{N} = 4 $ SYM , 2009, 0911.5332.

[41]  T. Binoth,et al.  An automatized algorithm to compute infrared divergent multi-loop integrals , 2000 .

[42]  M. Steinhauser,et al.  Quark and gluon form factors to three loops. , 2009, Physical review letters.

[43]  Thomas Hahn,et al.  Cuba - a library for multidimensional numerical integration , 2004, Comput. Phys. Commun..

[44]  Master integrals for massless three-loop form factors: One-loop and two-loop insertions , 2006, hep-ph/0607185.

[45]  Gudrun Heinrich,et al.  Sector Decomposition , 2008, 0803.4177.

[46]  C. Anastasiou,et al.  Real radiation at next-to-next-to-leading order in QCD: e+e- --> 2 jets through O(alpha2s). , 2004, Physical review letters.

[47]  Ericka Stricklin-Parker,et al.  Ann , 2005 .

[48]  K. Pohlmeyer.  Large momentum behavior of the Feynman amplitudes in the φ44 ‐theory , 1982 .

[49]  R. Bonciani,et al.  Two-Loop Planar Corrections to Heavy-Quark Pair Production in the Quark-Antiquark Channel , 2009, 0906.3671.

[50]  Kirill Melnikov,et al.  Higgs-Boson production at hadron colliders: differential cross sections through next-to-next-to-leading order. , 2004, Physical review letters.

[51]  M. Beneke,et al.  Asymptotic expansion of Feynman integrals near threshold , 1997, hep-ph/9711391.

[52]  D. Maitre,et al.  Master integrals for fermionic contributions to massless three-loop form-factors , 2007, 0711.3590.

[53]  G. Amdhal,et al.  Validity of the single processor approach to achieving large scale computing capabilities , 1967, AFIPS '67 (Spring).

[54]  An Algorithm for the high-energy expansion of multi-loop diagrams to next-to-leading logarithmic accuracy , 2004, hep-ph/0408068.

[55]  G. G. Stokes.  "J." , 1890, The New Yale Book of Quotations.

[56]  Li Lin Yang,et al.  Two-loop divergences of QCD scattering amplitudes with massive partons. , 2009, Physical review letters.

[57]  Kirill Melnikov,et al.  A new method for real radiation at next-to-next-to-leading order , 2004 .

[58]  V. Smirnov.  Applied Asymptotic Expansions in Momenta and Masses , 2001 .

[59]  Kirill Melnikov,et al.  The electron energy spectrum in muon decay through (α2) , 2005 .

[60]  Y. Kiyo,et al.  ${\cal O}(\alpha\alpha_s)$ corrections to the $\gamma t\bar{t}$ vertex at the top quark threshold , 2008 .

[61]  Klaus Hepp,et al.  Proof of the Bogoliubov-Parasiuk theorem on renormalization , 1966 .

[62]  Kirill Melnikov,et al.  Fully differential Higgs boson production and the di-photon signal through next-to-next-to-leading order , 2005 .

[63]  J. B. Tausk.  Non-planar massless two-loop Feynman diagrams with four on-shell legs , 1999 .

[64]  V. Smirnov.  Problems of the strategy of regions , 1999, hep-ph/9907471.

[65]  Radiative-nonrecoil corrections of order alpha^2 (Z alpha)^5 to the Lamb shift , 2009, 0911.4078.

[66]  M. Steinhauser,et al.  Matching QCD and HQET heavy–light currents at three loops , 2009, 0911.3356.

[67]  Two-loop electroweak angular-dependent logarithms at high energies , 2003, hep-ph/0301241.

[68]  Junpei Fujimoto,et al.  New implementation of the sector decomposition on FORM , 2009, 0902.2656.