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 → ∞, 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 w = 1,2 sums up to cyclotomy l = 20.Dedicated to Martinus Veltman on the occasion of his 80th birthday.

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 ±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 HarmonicSums.

We consider finite iterated generalized harmonic sums weighted by the binomial 2k k in numerators and denominators. A large class of these functions emerges in the calculation of massive Feynman diagrams with local operator insertions starting at 3-loop order in the coupling constant and extends the classes of the nested harmonic, generalized harmonic, and cyclotomic sums. The binomially weighted sums are associated by the Mellin transform to iterated integrals over square-root valued alphabets. The values of the sums for N → ∞ and the iterated integrals at x = 1 lead to new constants, extending the set of special numbers given by the multiple zeta values, the cyclotomic zeta values and special constants which emerge in the limit N → ∞ of generalized harmonic sums. We develop algorithms to obtain the Mellin representations of these sums in a systematic way. They are of importance for the derivation of the asymptotic expansion of these sums and their analytic continuation to N ∈ C. The associated convolution relations are derived for real parameters and can therefore be used in a wider context, as, e.g., for multi-scale processes. We also derive algorithms to transform iterated integrals over root-valued alphabets into binomial sums. Using generating functions we study a few aspects of infinite (inverse) binomial sums. C 2014 AIP Publishing LLC. [http://dx

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 Sigma and the method of hyperlogarithms. In some of the diagrams generalized harmonic sums with ξ ∈ {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 ≫ m 2 .

The pure singlet asymptotic heavy flavor corrections to 3-loop order for the deep-inelastic scattering structure function F 2 (x, Q 2 ) and the corresponding transition matrix element A (3),PS Qq in the variable flavor number scheme are computed. In Mellin-N space these inclusive quantities depend on generalized harmonic sums. We also recalculate the complete 3-loop pure singlet anomalous dimension for the first time. Numerical results for the Wilson coefficients, the operator matrix element and the contribution to the structure function F 2 (x, Q 2 ) are presented.