Computer algebra tools for Feynman integrals and related multi-sums

In perturbative calculations, e.g., in the setting of Quantum Chromodynamics (QCD) one aims at the evaluation of Feynman integrals. Here one is often faced with the problem to simplify multiple nested integrals or sums to expressions in terms of indefinite nested integrals or sums. Furthermore, one seeks for solutions of coupled systems of linear differential equations, that can be represented in terms of indefinite nested sums (or integrals). In this article we elaborate the main tools and the corresponding packages, that we have developed and intensively used within the last 10 years in the course of our QCD-calculations

Non-planar Feynman integrals, Mellin-Barnes representations, multiple sums

The construction of Mellin-Barnes (MB) representations for non-planar Feynman diagrams and the summation of multiple series derived from general MB representations are discussed. A basic version of a new package AMBREv.3.0 is supplemented. The ultimate goal of this project is the automatic evaluation of MB representations for multiloop scalar and tensor Feynman integrals through infinite sums, preferably with analytic solutions. We shortly describe a strategy of further algebraic summation.Comment: Contribution to the proceedings of the Loops and Legs 2014 conferenc

Recent Symbolic Summation Methods to Solve Coupled Systems of Differential and Difference Equations

We outline a new algorithm to solve coupled systems of differential equations in one continuous variable $x$ (resp. coupled difference equations in one discrete variable $N$) depending on a small parameter $\epsilon$: given such a system and given sufficiently many initial values, we can determine the first coefficients of the Laurent-series solutions in $\epsilon$ if they are expressible in terms of indefinite nested sums and products. This systematic approach is based on symbolic summation algorithms in the context of difference rings/fields and uncoupling algorithms. The proposed method gives rise to new interesting applications in connection with integration by parts (IBP) methods. As an illustrative example, we will demonstrate how one can calculate the $\epsilon$-expansion of a ladder graph with 6 massive fermion lines

New Results on Massive 3-Loop Wilson Coefficients in Deep-Inelastic Scattering

We present recent results on newly calculated 2- and 3-loop contributions to the heavy quark parts of the structure functions in deep-inelastic scattering due to charm and bottom.Comment: Contribution to the Proc. of Loops and Legs 2016, PoS, in prin

Nested (inverse) binomial sums and new iterated integrals for massive Feynman diagrams

Nested sums containing binomial coefficients occur in the computation of massive operator matrix elements. Their associated iterated integrals lead to alphabets including radicals, for which we determined a suitable basis. We discuss algorithms for converting between sum and integral representations, mainly relying on the Mellin transform. To aid the conversion we worked out dedicated rewrite rules, based on which also some general patterns emerging in the process can be obtained.Comment: 13 pages LATEX, one style file, Proceedings of Loops and Legs in Quantum Field Theory -- LL2014,27 April 2014 -- 02 May 2014 Weimar, German

Calculation of Massless Feynman Integrals using Harmonic Sums

A method for the evaluation of the epsilon expansion of multi-loop massless Feynman integrals is introduced. This method is based on the Gegenbauer polynomial technique and the expansion of the Gamma function in terms of harmonic sums. Algorithms for the evaluation of nested and harmonic sums are used to reduce the expressions to get analytical or numerical results for the expansion coefficients. Methods to increase the precision of numerical results are discussed.Comment: 30 pages, 6 figures; Minor typos corrected, references added. Published in Computer Physics Communication

Recent progress on the calculation of three-loop heavy flavor Wilson coefficients in deep-inelastic scattering

We report on our latest results in the calculation of the three-loop heavy flavor contributions to the Wilson coefficients in deep-inelastic scattering in the asymptotic region $Q^2 \gg m^2$. We discuss the different methods used to compute the required operator matrix elements and the corresponding Feynman integrals. These methods very recently allowed us to obtain a series of new operator matrix elements and Wilson coefficients like the flavor non-singlet and pure singlet Wilson coefficients.Comment: 11 pages Latex, 2 Figures, Proc. of Loops and Legs in Quantum Field Theory, April 2014, Weimar, German

A Fast Approach to Creative Telescoping

In this note we reinvestigate the task of computing creative telescoping relations in differential-difference operator algebras. Our approach is based on an ansatz that explicitly includes the denominators of the delta parts. We contribute several ideas of how to make an implementation of this approach reasonably fast and provide such an implementation. A selection of examples shows that it can be superior to existing methods by a large factor.Comment: 9 pages, 1 table, final version as it appeared in the journa