Elsevier Science Home
Computer Physics Communications Program Library
Full text online from Science Direct
Programs in Physics & Physical Chemistry
CPC Home

[Licence| Download | New Version Template] aeuv_v1_0.tar.gz(2483 Kbytes)
Manuscript Title: Algorithms for the symbolic integration of hyperlogarithms with applications to Feynman integrals
Authors: Erik Panzer
Program title: HyperInt
Catalogue identifier: AEUV_v1_0
Distribution format: tar.gz
Journal reference: Comput. Phys. Commun. 188(2015)148
Programming language: Maple [1], version 16 or higher.
Computer: Any that supports Maple.
Operating system: Any that supports Maple.
RAM: Highly problem dependent; from a few MiB to many GiB
Keywords: Hyperlogarithms, Polylogarithms, Symbolic integration, Computer algebra, Feynman integrals, ε-expansions.
PACS: 02.30.Gp, 02.70.Wz, 12.38.Bx.
Classification: 4.4, 5.

Nature of problem:
Feynman integrals and their ε-expansions in dimensional regularization can be expressed in the Schwinger parametrization as multi-dimensional integrals of rational functions and logarithms. Symbolic integration of such functions therefore serves a tool for the exact and direct evaluation of Feynman graphs.

Solution method:
Symbolic integration of rational linear combinations of polylogarithms of rational arguments is obtained using a representation in terms of hyperlogarithms. The algorithms exploit their iterated integral structure.

To compute multi-dimensional integrals with this method, the integrand must be linearly reducible, a criterion we state in section 4. As a consequence, only a small subset of all Feynman integrals can be addressed.

Unusual features:
The complete program works strictly symbolically and the obtained results are exact. Whenever a Feynman graph is linearly reducible, its ε-expansion can be computed to arbitrary order (subject only to time and memory restrictions) in ε, near any even dimension of space-time and for arbitrarily ε-dependent powers of propagators with integer values at ε = 0. The method is not restricted to scalar integrals and applies even to (regulated) divergent integrals.

Additional comments:
Apart from Feynman integrals, other suitable parametric integrals may be computed (or expanded in ε) as well, like for example hypergeometric functions.
An example worksheet, Manual.mw, is included. This contains an explanation of most features provided and includes plenty of examples of Feynman integral computations.

Running time:
Highly dependent on the particular problem through the number of integrations to be performed (edges of a graph), the number of remaining variables (kinematic invariants), the order in ε and the complexity of the geometry (topology of the graph). Simplest examples finish in seconds, but the time needed increases beyond any bound for sufficiently high orders in ε or graphs with many edges.

[1] Maple 161. Maplesoft, a division of Waterloo Maple Inc., Waterloo, Ontario.