Source author record

Jakob Ablinger

Jakob Ablinger appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

15works
8topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

15 published item(s)

preprint2021arXiv

Extensions of the AZ-algorithm and the Package MultiIntegrate

We extend the (continuous) multivariate Almkvist-Zeilberger algorithm in order to apply it for instance to special Feynman integrals emerging in renormalizable Quantum field Theories. We will consider multidimensional integrals over hyperexponential integrands and try to find closed form representations in terms of nested sums and products or iterated integrals. In addition, if we fail to compute a closed form solution in full generality, we may succeed in computing the first coefficients of the Laurent series expansions of such integrals in terms of indefinite nested sums and products or iterated integrals. In this article we present the corresponding methods and algorithms. Our Mathematica package MultiIntegrate, can be considered as an enhanced implementation of the (continuous) multivariate Almkvist Zeilberger algorithm to compute recurrences or differential equations for hyperexponential integrands and integrals. Together with the summation package Sigma and the package HarmonicSums our package provides methods to compute closed form representations (or coefficients of the Laurent series expansions) of multidimensional integrals over hyperexponential integrands in terms of nested sums or iterated integrals.

preprint2021arXiv

Solving linear difference equations with coefficients in rings with idempotent representations

We introduce a general reduction strategy that enables one to search for solutions of parameterized linear difference equations in difference rings. Here we assume that the ring itself can be decomposed by a direct sum of integral domains (using idempotent elements) that enjoys certain technical features and that the coefficients of the difference equation are not degenerated. Using this mechanism we can reduce the problem to find solutions in a ring (with zero-divisors) to search solutions in several copies of integral domains. Utilizing existing solvers in this integral domain setting, we obtain a general solver where the components of the linear difference equations and the solutions can be taken from difference rings that are built e.g., by $RΠΣ$-extensions over $ΠΣ$-fields. This class of difference rings contains, e.g., nested sums and products, products over roots of unity and nested sums defined over such objects.

preprint2016arXiv

A toolbox to solve coupled systems of differential and difference equations

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.

preprint2016arXiv

Algorithms to solve coupled systems of differential equations in terms of power series

Using integration by parts relations, Feynman integrals can be represented in terms of coupled systems of differential equations. In the following we suppose that the unknown Feynman integrals can be given in power series representations, and that sufficiently many initial values of the integrals are given. Then there exist algorithms that decide constructively if the coefficients of their power series representations can be given within the class of nested sums over hypergeometric products. In this article we will work out the calculation steps that solve this problem. First, we will present a successful tactic that has been applied recently to challenging problems coming from massive 3-loop Feynman integrals. Here our main tool is to solve scalar linear recurrences within the class of nested sums over hypergeometric products. Second, we will present a new variation of this tactic which relies on more involved summation technologies but succeeds in reducing the problem to solve scalar recurrences with lower recurrence orders. The article will work out the different challenges of this new tactic and demonstrates how they can be treated efficiently with our existing summation technologies.

preprint2015arXiv

Discovering and Proving Infinite Binomial Sums Identities

We consider binomial and inverse binomial sums at infinity and rewrite them in terms of a small set of constants, such as powers of $π$ or $\log(2)$. In order to perform these simplifications, we view the series as specializations of generating series. For these generating series, we derive integral representations in terms of root-valued iterated integrals. Using substitutions, we express the interated integrals as cyclotomic harmonic polylogarithms. Finally, by applying known relations among the cyclotomic harmonic polylogarithms, we derive expressions in terms of several constants.

preprint2014arXiv

Calculating Massive 3-loop Graphs for Operator Matrix Elements by the Method of Hyperlogarithms

We calculate convergent 3-loop Feynman diagrams containing a single massive loop equipped with twist $τ=2$ local operator insertions corresponding to spin $N$. They contribute to the massive operator matrix elements in QCD describing the massive Wilson coefficients for deep-inelastic scattering at large virtualities. Diagrams of this kind can be computed using an extended version to the method of hyperlogarithms, originally being designed for massless Feynman diagrams without operators. The method is applied to Benz- and $V$-type graphs, belonging to the genuine 3-loop topologies. In case of the $V$-type graphs with five massive propagators new types of nested sums and iterated integrals emerge. The sums are given in terms of finite binomially and inverse binomially weighted generalized cyclotomic sums, while the 1-dimensionally iterated integrals are based on a set of $\sim 30$ square-root valued letters. We also derive the asymptotic representations of the nested sums and present the solution for $N \in \mathbb{C}$. Integrals with a power-like divergence in $N$--space $\propto a^N, a \in \mathbb{R}, a > 1,$ for large values of $N$ emerge. They still possess a representation in $x$--space, which is given in terms of root-valued iterated integrals in the present case. The method of hyperlogarithms is also used to calculate higher moments for crossed box graphs with different operator insertions.

preprint2014arXiv

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.

preprint2014arXiv

The package HarmonicSums: Computer Algebra and Analytic aspects of Nested Sums

This paper summarizes the essential functionality of the computer algebra package HarmonicSums. On the one hand HarmonicSums can work with nested sums such as harmonic sums and their generalizations and on the other hand it can treat iterated integrals of the Poincare and Chen-type, such as harmonic polylogarithms and their generalizations. The interplay of these representations and the analytic aspects are illustrated by concrete examples.

preprint2013arXiv

Analytic and Algorithmic Aspects of Generalized Harmonic Sums and Polylogarithms

In recent three--loop calculations of massive Feynman integrals within Quantum Chromodynamics (QCD) and, e.g., in recent combinatorial problems the so-called generalized harmonic sums (in short $S$-sums) arise. They are characterized by rational (or real) numerator weights also different from $\pm 1$. In this article we explore the algorithmic and analytic properties of these sums systematically. We work out the Mellin and inverse Mellin transform which connects the sums under consideration with the associated Poincaré iterated integrals, also called generalized harmonic polylogarithms. In this regard, we obtain explicit analytic continuations by means of asymptotic expansions of the $S$-sums which started to occur frequently in current QCD calculations. In addition, we derive algebraic and structural relations, like differentiation w.r.t. the external summation index and different multi-argument relations, for the compactification of $S$-sum expressions. Finally, we calculate algebraic relations for infinite $S$-sums, or equivalently for generalized harmonic polylogarithms evaluated at special values. The corresponding algorithms and relations are encoded in the computer algebra package {\tt HarmonicSums}.

preprint2013arXiv

Computer Algebra Algorithms for Special Functions in Particle Physics

This work deals with special nested objects arising in massive higher order perturbative calculations in renormalizable quantum field theories. On the one hand we work with nested sums such as harmonic sums and their generalizations (S-sums, cyclotomic harmonic sums, cyclotomic S-sums) and on the other hand we treat iterated integrals of the Poincaré and Chen-type, such as harmonic polylogarithms and their generalizations (multiple polylogarithms, cyclotomic harmonic polylogarithms). The iterated integrals are connected to the nested sums via (generalizations of) the Mellin-transformation and we show how this transformation can be computed. We derive algebraic and structural relations between the nested sums as well as relations between the values of the sums at infinity and connected to it the values of the iterated integrals evaluated at special constants. In addition we state algorithms to compute asymptotic expansions of these nested objects and we state an algorithm which rewrites certain types of nested sums into expressions in terms of cyclotomic S-sums. Moreover we summarize the main functionality of the computer algebra package HarmonicSums in which all these algorithms and transformations are implemented. Furthermore, we present application of and enhancements of the multivariate Almkvist-Zeilberger algorithm to certain types of Feynman integrals and the corresponding computer algebra package MultiIntegrate.

preprint2013arXiv

Harmonic Sums, Polylogarithms, Special Numbers, and their Generalizations

In these introductory lectures we discuss classes of presently known nested sums, associated iterated integrals, and special constants which hierarchically appear in the evaluation of massless and massive Feynman diagrams at higher loops. These quantities are elements of stuffle and shuffle algebras implying algebraic relations being widely independent of the special quantities considered. They are supplemented by structural relations. The generalizations are given in terms of generalized harmonic sums, (generalized) cyclotomic sums, and sums containing in addition binomial and inverse-binomial weights. To all these quantities iterated integrals and special numbers are associated. We also discuss the analytic continuation of nested sums of different kind to complex values of the external summation bound N.

preprint2013arXiv

Massive 3-loop Ladder Diagrams for Quarkonic Local Operator Matrix Elements

3-loop diagrams of the ladder-type, which emerge for local quarkonic twist-2 operator matrix elements, are computed directly for general values of the Mellin variable $N$ using Appell-function representations and applying modern summation technologies provided by the package {\sf Sigma} and the method of hyperlogarithms. In some of the diagrams generalized harmonic sums with $ξ\in \{1,1/2,2\}$ emerge beyond the usual nested harmonic sums. As the asymptotic representation of the corresponding integrals shows, the generalized sums conspire giving well behaved expressions for large values of $N$. These diagrams contribute to the 3-loop heavy flavor Wilson coefficients of the structure functions in deep-inelastic scattering in the region $Q^2 \gg m^2$.

preprint2011arXiv

Harmonic Sums and Polylogarithms Generated by Cyclotomic Polynomials

The computation of Feynman integrals in massive higher order perturbative calculations in renormalizable Quantum Field Theories requires extensions of multiply nested harmonic sums, which can be generated as real representations by Mellin transforms of Poincaré--iterated integrals including denominators of higher cyclotomic polynomials. We derive the cyclotomic harmonic polylogarithms and harmonic sums and study their algebraic and structural relations. The analytic continuation of cyclotomic harmonic sums to complex values of $N$ is performed using analytic representations. We also consider special values of the cyclotomic harmonic polylogarithms at argument $x=1$, resp., for the cyclotomic harmonic sums at $N \rightarrow \infty$, which are related to colored multiple zeta values, deriving various of their relations, based on the stuffle and shuffle algebras and three multiple argument relations. We also consider infinite generalized nested harmonic sums at roots of unity which are related to the infinite cyclotomic harmonic sums. Basis representations are derived for weight {\sf w = 1,2} sums up to cyclotomy {\sf l = 20}.

preprint2010arXiv

A Computer Algebra Toolbox for Harmonic Sums Related to Particle Physics

In this work we present the computer algebra package HarmonicSums and its theoretical background for the manipulation of harmonic sums and some related quantities as for example Euler-Zagier sums and harmonic polylogarithms. Harmonic sums and generalized harmonic sums emerge as special cases of so-called d'Alembertian solutions of recurrence relations. We show that harmonic sums form a quasi-shuffle algebra and describe a method how we can find algebraically independent harmonic sums. In addition, we define a differentiation on harmonic sums via an extended version of the Mellin transform. Along with that, new relations between harmonic sums will arise. Furthermore, we present an algorithm which rewrites certain types of nested sums into expressions in terms of harmonic sums. We illustrate by nontrivial examples how these algorithms in cooperation with the summation package Sigma support the evaluation of Feynman integrals.