epsilon : A tool to find a canonical basis of master integrals

Abstract In 2013, Henn proposed a special basis for a certain class of master integrals, which are expressible in terms of iterated integrals. In this basis, the master integrals obey a differential equation, where the right hand side is proportional to ϵ in d = 4 − 2 ϵ space–time dimensions. An algorithmic approach to find such a basis was found by Lee. We present the tool epsilon , an efficient implementation of Lee’s algorithm based on the Fermat computer algebra system as computational back end. Program summary Program Title: epsilon Program Files doi: http://dx.doi.org/10.17632/j59sy5n729.1 Licensing provisions: GPLv3 Programming language: C++ Nature of problem: For a certain class of master integrals, a canonical basis can be found in which they fulfill a differential equation with the right hand side proportional to ϵ . In such a basis the solution of the master integrals in an ϵ -expansion becomes trivial. Unfortunately, the problem of finding a canonical basis is challenging. Solution method: Algorithm by Lee [1] Restrictions: The normalization step of Lee’s algorithm will fail if the eigenvalues of the matrix residues are not of the form a + b ϵ with a , b ∈ Z . Multi-scale problems are not supported. [1] R.N. Lee, JHEP 1504 (2015) 108 [ arXiv:1411.0911 [hep-ph]].

[1]  Daniele Binosi,et al.  JaxoDraw: A graphical user interface for drawing Feynman diagrams , 2004, Comput. Phys. Commun..

[2]  Richard Kreckel,et al.  Introduction to the GiNaC Framework for Symbolic Computation within the C++ Programming Language , 2000, J. Symb. Comput..

[3]  R. Lee Reducing differential equations for multiloop master integrals , 2014, 1411.0911.

[4]  A. Kotikov Differential equation method: The Calculation of N point Feynman diagrams , 1991 .

[5]  John C. Collins,et al.  JaxoDraw: A graphical user interface for drawing Feynman diagrams. Version 2.0 release notes , 2009, Comput. Phys. Commun..

[6]  A. Kotikov NEW METHOD OF MASSIVE FEYNMAN DIAGRAMS CALCULATION , 1991 .

[7]  J. Vermaseren,et al.  Harmonic Polylogarithms , 1999, hep-ph/9905237.

[8]  A. Goncharov,et al.  Classical polylogarithms for amplitudes and Wilson loops. , 2010, Physical review letters.

[9]  J.A.M. Vermaseren,et al.  New features of FORM , 2000 .

[10]  Oleksandr Gituliar,et al.  Fuchsia : A tool for reducing differential equations for Feynman master integrals to epsilon form , 2017, Comput. Phys. Commun..

[11]  J. Henn Multiloop integrals in dimensional regularization made simple. , 2013, Physical review letters.

[12]  V. Smirnov Analytic Tools for Feynman Integrals , 2013 .

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

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

[15]  O. Gituliar Master integrals for splitting functions from differential equations in QCD , 2015, 1607.00759.

[16]  D. F. Lodato,et al.  Λ c + production in pp collisions at √s=7 TeV and in p-Pb collisions at √sNN=5.02 TeV , 2017, 1712.09581.

[17]  A. Kotikov,et al.  Differential equations method: the calculation of vertex-type Feynman diagrams , 1991 .

[18]  J. Moser,et al.  The order of a singularity in Fuchs' theory , 1959 .

[19]  Daniel Maître HPL, a Mathematica implementation of the harmonic polylogarithms , 2006, Comput. Phys. Commun..