BZPEERPapers, people, evidence and research pathsDiscover
Research tools
Network
Opportunities
Account
Source author record
Carsten Schneider appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.
Catalog footprint
Research graph
Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.
BZPEER is loading the nearby papers, people, topics and institutions for this page.
Published work
This is an introduction to, and invitation to read, a series of review articles on scattering amplitudes in gauge theory, gravity, and superstring theory. Our aim is to provide an overview of the field, from basic aspects to a selection of current (2022) research and developments.
The analytic integration and simplification of multi-loop Feynman integrals to special functions and constants plays an important role to perform higher order perturbative calculations in the Standard Model of elementary particles. In this survey article the most recent and relevant computer algebra and special function algorithms are presented that are currently used or that may play an important role to perform such challenging precision calculations in the future. They are discussed in the context of analytic zero, single and double scale calculations in the Quantum Field Theories of the Standard Model and effective field theories, also with classical applications. These calculations play a central role in the analysis of precision measurements at present and future colliders to obtain ultimate information for fundamental physics.
We present recent computer algebra methods that support the calculations of (multivariate) series solutions for (certain coupled systems of partial) linear differential equations. The summand of the series solutions may be built by hypergeometric products and more generally by indefinite nested sums defined over such products. Special cases are hypergeometric structures such as Appell-functions or generalizations of them that arise frequently when dealing with parameter Feynman integrals.
We present a complete algorithm that computes all hypergeometric solutions of homogeneous linear difference equations and rational solutions of parameterized linear difference equations in the setting of $ΠΣ^*$-fields. More generally, we provide a flexible framework for a big class of difference fields that is built by a tower of $ΠΣ^*$-field extensions over a difference field that satisfies certain algorithmic properties. As a consequence one can compute all solutions in terms of indefinite nested sums and products that arise within the components of a parameterized linear difference equation, and one can find all hypergeometric solutions that are defined over the arising sums and products of a homogeneous linear difference equation.
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.
Using symbolic summation tools in the setting of difference rings, we prove a two-parametric identity that relates rational approximations to $ζ(4)$.
In (S.B. Ekhad and D. Zeilberger, 2020) an exciting case study has been initiated in which experimental mathematics and symbolic computation are utilized to discover new properties concerning the so-called Absent-Minded Passengers Problem. Based on these results, Doron Zeilberger raised some challenging tasks to gain further probabilistic insight. In this note we report on this enterprise. In particular, we demonstrate how the computer algebra packages of RISC can be used to carry out the underlying heavy calculations.
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.
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.
Recently, $RΠΣ^*$-extensions have been introduced which extend Karr's $ΠΣ^*$-fields substantially: one can represent expressions not only in terms of transcendental sums and products, but one can work also with products over primitive roots of unity. Since one can solve the parameterized telescoping problem in such rings, covering as special cases the summation paradigms of telescoping and creative telescoping, one obtains a rather flexible toolbox for symbolic summation. This article is the continuation of this work. Inspired by Singer's Galois theory of difference equations we will work out several alternative characterizations of $RΠΣ^*$-extensions: adjoining naively sums and products leads to an $RΠΣ^*$-extension iff the obtained difference ring is simple iff the ring can be embedded into the ring of sequences iff the ring can be given by the interlacing of $ΠΣ^*$-extensions. From the viewpoint of applications this leads to a fully automatic machinery to represent indefinite nested sums and products in such $RΠΣ^*$-rings. In addition, we work out how the parameterized telescoping paradigm can be used to prove algebraic independence of indefinite nested sums. Furthermore, one obtains an alternative reduction tactic to solve the parameterized telescoping problem in basic $RΠΣ^*$-extensions exploiting the interlacing property.
A summation framework is developed that enhances Karr's difference field approach. It covers not only indefinite nested sums and products in terms of transcendental extensions, but it can treat, e.g., nested products defined over roots of unity. The theory of the so-called $RΠΣ^*$-extensions is supplemented by algorithms that support the construction of such difference rings automatically and that assist in the task to tackle symbolic summation problems. Algorithms are presented that solve parameterized telescoping equations, and more generally parameterized first-order difference equations, in the given difference ring. As a consequence, one obtains algorithms for the summation paradigms of telescoping and Zeilberger's creative telescoping. With this difference ring theory one obtains a rigorous summation machinery that has been applied to numerous challenging problems coming, e.g., from combinatorics and particle physics.
We present an algebraic framework to represent indefinite nested sums over hypergeometric expressions in difference rings. In order to accomplish this task, parts of Karr's difference field theory have been extended to a ring theory in which also the alternating sign can be expressed. The underlying machinery relies on algorithms that compute all solutions of a given parameterized telescoping equation. As a consequence, we can solve the telescoping and creative telescoping problem in such difference rings.
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.
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.
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.
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 $ε$: given such a system and given sufficiently many initial values, we can determine the first coefficients of the Laurent-series solutions in $ε$ 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 $ε$-expansion of a ladder graph with 6 massive fermion lines.
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}.
We employ computer algebra algorithms to prove a collection of identities involving Bessel functions with half-integer orders and other special functions. These identities appear in the famous Handbook of Mathematical Functions, as well as in its successor, the DLMF, but their proofs were lost. We use generating functions and symbolic summation techniques to produce new proofs for them.
Parameterized telescoping (including telescoping and creative telescoping) and refined versions of it play a central role in the research area of symbolic summation. Karr introduced 1981 $ΠΣ$-fields, a general class of difference fields, that enables one to consider this problem for indefinite nested sums and products covering as special cases, e.g., the ($q$--)hypergeometric case and their mixed versions. This survey article presents the available algorithms in the framework of $ΠΣ$-extensions and elaborates new results concerning efficiency.
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$.
A large class of Feynman integrals, like e.g., two-point parameter integrals with at most one mass and containing local operator insertions, can be transformed to multi-sums over hypergeometric expressions. In this survey article we present a difference field approach for symbolic summation that enables one to simplify such definite nested sums to indefinite nested sums. In particular, the simplification is given -if possible- in terms of harmonic sums, generalized harmonic sums, cyclotomic harmonic sums or binomial sums. Special emphasis is put on the developed packages Sigma, EvaluateMultiSums and SumProduction that assist in the task to perform these simplifications completely automatically for huge input expressions.
In this survey article we present difference field algorithms for symbolic summation. Special emphasize is put on new aspects in how the summation problems are rephrased in terms of difference fields, how the problems are solved there, and how the derived results in the given difference field can be reinterpreted as solutions of the input problem. The algorithms are illustrated with the Mathematica package \SigmaP\ by discovering and proving new harmonic number identities extending those from (Paule and Schneider, 2003). In addition, the newly developed package \texttt{EvaluateMultiSums} is introduced that combines the presented tools. In this way, large scale summation problems for the evaluation of Feynman diagrams in QCD (Quantum ChromoDynamics) can be solved completely automatically.
Given a Feynman parameter integral, depending on a single discrete variable $N$ and a real parameter $ε$, we discuss a new algorithmic framework to compute the first coefficients of its Laurent series expansion in $ε$. In a first step, the integrals are expressed by hypergeometric multi-sums by means of symbolic transformations. Given this sum format, we develop new summation tools to extract the first coefficients of its series expansion whenever they are expressible in terms of indefinite nested product-sum expressions. In particular, we enhance the known multi-sum algorithms to derive recurrences for sums with complicated boundary conditions, and we present new algorithms to find formal Laurent series solutions of a given recurrence relation.
The $O(α_s^3 n_f T_F^2 C_{A,F})$ terms to the massive gluonic operator matrix elements are calculated for general values of the Mellin variable $N$. These twist-2 matrix elements occur as transition functions in the variable flavor number scheme at NNLO. The calculation uses sum-representations in generalized hypergeometric series turning into harmonic sums. The analytic continuation to complex values of $N$ is provided.
We continue to investigate which polynomials can possibly occur as factors in the denominators of rational solutions of a given partial linear difference equation. In an earlier article we had introduced the distinction between periodic and aperiodic factors in the denominator, and we gave an algorithm for predicting the aperiodic ones. Now we extend this technique towards the periodic case and present a refined algorithm which also finds most of the periodic factors.
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}.
We investigate which polynomials can possibly occur as factors in the denominators of rational solutions of a given partial linear difference equation (PLDE). Two kinds of polynomials are to be distinguished, we call them /periodic/ and /aperiodic/. The main result is a generalization of a well-known denominator bounding technique for univariate equations to PLDEs. This generalization is able to find all the aperiodic factors of the denominators for a given PLDE.
The scalar form factor of the nucleon and related physical quantities are investigated in the framework of the semibosonized SU(3) Nambu-Jona-Lasinio soliton model. We take into account the rotational $1/N_c$ corrections and linear $m_s$ corrections. The strangeness content of the nucleon in the scalar form factor is discussed in detail. In particular, it is found that the $m_s$ corrections play an essential role of reducing the $\langle N | \bar{s} s | N \rangle$ arising from the leading order and rotational $1/N_c$ contributions. We obtain the $σ_{πN} (0)=40.80\;\mbox{MeV}$, $Δσ= σ_{πN} (2m^{2}_π)-σ_{πN} (0) = 18.18\;\mbox{MeV}$ and $\langle r^2\rangle^{S}_{N} = 1.50\;\mbox{fm}^2$. The results are in a remarkable agreement with empirical data analyzed by Gasser, Leutwyler, and Sainio~\cite{gls}.