A toolbox to solve coupled systems of differential and difference equations

Abstract

We present algorithms to solve coupled systems of linear differential equations, arising in the calculation of massive Feynman diagrams with local operator insertions at 3-loop order, which do {\it not} request special choices of bases. Here we assume that the desired solution has a power series representation and we seek for the coefficients in closed form. In particular, if the coefficients depend on a small parameter $\ep$ (the dimensional parameter), we assume that the coefficients themselves can be expanded in formal Laurent series w.r.t.\ $\ep$ and we try to compute the first terms in closed form. More precisely, we have a decision algorithm which solves the following problem: if the terms can be represented by an indefinite nested hypergeometric sum expression (covering as special cases the harmonic sums, cyclotomic sums, generalized harmonic sums or nested binomial sums), then we can calculate them. If the algorithm fails, we obtain a proof that the terms cannot be represented by the class of indefinite nested hypergeometric sum expressions. Internally, this problem is reduced by holonomic closure properties to solving a coupled system of linear difference equations. The underlying method in this setting relies on decoupling algorithms, difference ring algorithms and recurrence solving. We demonstrate by a concrete example how this algorithm can be applied with the new Mathematica package \texttt{SolveCoupledSystem} which is based on the packages \texttt{Sigma}, \texttt{HarmonicSums} and \texttt{OreSys}. In all applications the representation in $x$-space is obtained as an iterated integral representation over general alphabets, generalizing Poincaré iterated integrals.

References in corpus (11)

• A Computer Algebra Toolbox for Harmonic Sums Related to Particle Physics
• A Refined Difference Field Theory for Symbolic Summation
• Calculating Three Loop Ladder and V-Topologies for Massive Operator Matrix Elements by Computer Algebra
• The 3-Loop Non-Singlet Heavy Flavor Contributions and Anomalous Dimensions for the Structure Function $F_2(x,Q^2)$ and Transversity
• The 3-Loop Pure Singlet Heavy Flavor Contributions to the Structure Function $F_2(x,Q^2)$ and the Anomalous Dimension
• The O(α_s^3 T_F^2) Contributions to the Gluonic Operator Matrix Element
• The O(α_s^3) Heavy Flavor Contributions to the Charged Current Structure Function xF_3(x,Q^2) at Large Momentum Transfer
• Parameterized Telescoping Proves Algebraic Independence of Sums
• A Difference Ring Theory for Symbolic Summation
• Recent Symbolic Summation Methods to Solve Coupled Systems of Differential and Difference Equations
• Complexity Estimates for Two Uncoupling Algorithms

Cited by in corpus (2)

• 3-Loop Corrections to the Heavy Flavor Wilson Coefficients in Deep-Inelastic Scattering
• Three loop QCD corrections to heavy quark form factors