Search

The unpolarized and polarized massive operator matrix elements $A_{Qg}^{(3)}$ and $\Delta A_{Qg}^{(3)}$ contain first-order factorizable and non-first-order factorizable contributions in the determining difference or differential equations of their master integrals. We compute their first-order factorizable contributions in the single heavy mass case for all contributing Feynman diagrams. Moreover, we present the complete color-$\zeta$ factors for the cases in which also non-first-order factorizable contributions emerge in the master integrals, but cancel in the final result as found by using the method of arbitrary high Mellin moments. Individual contributions depend also on generalized harmonic sums and on nested finite binomial and inverse binomial sums in Mellin $N$-space, and correspondingly, on Kummer-Poincar\'e and square-root valued alphabets in Bjorken-$x$ space. We present a complete discussion of the possibilities of solving the present problem in $N$-space analytically and we also discuss the limitations in the present case to analytically continue the given $N$-space expressions to $N \in \mathbb{C}$ by strict methods. The representation through generating functions allows a well synchronized representation of the first-order factorizable results over a 17-letter alphabet. We finally obtain representations in terms of iterated integrals over the corresponding alphabet in $x$-space, also containing up to weight {\sf w = 5} special constants, which can be rationalized to Kummer-Poincar\'e iterated integrals at special arguments. The analytic $x$-space representation requires separate analyses for the intervals $x \in [0,1/4], [1/4,1/2], [1/2,1]$ and $x > 1$. We also derive the small and large $x$ limits of the first-order factorizable contributions., Comment: 58 pages, 4 Figures

We calculate the massless unpolarized Wilson coefficients for deeply inelastic scattering for the structure functions $F_2(x,Q^2), F_L(x,Q^2), x F_3(x,Q^2)$ in the $\overline{\sf MS}$ scheme and the polarized Wilson coefficients of the structure function $g_1(x,Q^2)$ in the Larin scheme up to three--loop order in QCD in a fully automated way based on the method of arbitrary high Mellin moments. We work in the Larin scheme in the case of contributing axial--vector couplings or polarized nucleons. For the unpolarized structure functions we compare to results given in the literature. The polarized three--loop Wilson coefficients are calculated for the first time. As a by--product we also obtain the quarkonic three--loop anomalous dimensions from the $O(1/\varepsilon)$ terms of the unrenormalized forward Compton amplitude. Expansions for small and large values of the Bjorken variable $x$ are provided., Comment: 159 pages Latex

We report on the calculation of the three-loop polarized and unpolarized flavor non-singlet and the polarized singlet anomalous dimensions using massless off-shell operator matrix elements in a gauge-variant framework. We also reconsider the unpolarized two-loop singlet anomalous dimensions and correct errors in the foregoing literature., Comment: 12 pages Latex, Contributeion to the Proceedings of LL2022

We calculate the unpolarized and polarized two--loop massless off--shell operator matrix elements in QCD to $O(\varepsilon)$ in the dimensional parameter in an automated way. Here we use the method of arbitrary high Mellin moments and difference ring theory, based on integration-by-parts relations. This method also constitutes one way to compute the QCD anomalous dimensions. The presented higher order contributions to these operator matrix elements occur as building blocks in the corresponding higher order calculations up to four--loop order. All contributing quantities can be expressed in terms of harmonic sums in Mellin--$N$ space or by harmonic polylogarithms in $z$--space. We also perform comparisons to the literature., Comment: 101 pages Latex

Hypergeometric structures in single and multiscale Feynman integrals emerge in a wide class of topologies. Using integration-by-parts relations, associated master or scalar integrals have to be calculated. For this purpose it appears useful to devise an automated method which recognizes the respective (partial) differential equations related to the corresponding higher transcendental functions. We solve these equations through associated recursions of the expansion coefficient of the multivalued formal Taylor series. The expansion coefficients can be determined using either the package {\tt Sigma} in the case of linear difference equations or by applying heuristic methods in the case of partial linear difference equations. In the present context a new type of sums occurs, the Hurwitz harmonic sums, and generalized versions of them. The code {\tt HypSeries} transforming classes of differential equations into analytic series expansions is described. Also partial difference equations having rational solutions and rational function solutions of Pochhammer symbols are considered, for which the code {\tt solvePartialLDE} is designed. Generalized hypergeometric functions, Appell-,~Kamp\'e de F\'eriet-, Horn-, Lauricella-Saran-, Srivasta-, and Exton--type functions are considered. We illustrate the algorithms by examples., Comment: 55 pages, several anc. files

Future high luminosity polarized deep--inelastic scattering experiments will improve both the knowledge of the spin sub--structure of the nucleons and contribute further to the precision determination of the strong coupling constant, as well as, reveal currently yet unknown higher twist contributions in the polarized sector. For all these tasks to be performed, it is necessary to know the QCD leading twist scaling violations of the measured structure functions. Here an important ingredient consists in the polarized singlet anomalous dimensions and splitting functions in QCD. We recalculate these quantities to three--loop order in the M--scheme by using the traditional method of space--like off--shell massless operator matrix elements, being a gauge--dependent framework. Here one obtains the anomalous dimensions without referring to gravitational currents, needed when calculating them using the forward Compton amplitude. We also calculate the non--singlet splitting function $\Delta P_{\rm qq}^{(2), \rm s, NS}$ and compare the final results to the literature, also including predictions for the region of small values of Bjorken $x$., Comment: 28 pages Latex, 1 figure

We calculate the unpolarized and polarized three--loop anomalous dimensions and splitting functions $P_{\rm NS}^+, P_{\rm NS}^-$ and $P_{\rm NS}^{\rm s}$ in QCD in the $\overline{\sf MS}$ scheme by using the traditional method of space--like off shell massless operator matrix elements. This is a gauge--dependent framework. For the first time we also calculate the three--loop anomalous dimensions $P_{\rm NS}^{\rm \pm tr}$ for transversity directly. We compare our results to the literature., Comment: 40 pages Latex

We compute the logarithmic contributions to the polarized massive Wilson coefficients for deep-inelastic scattering in the asymptotic region $Q^2 \gg m^2$ to 3-loop order in the fixed-flavor number scheme and present the corresponding expressions for the polarized massive operator matrix elements needed in the variable flavor number scheme. The calculation is performed in the Larin scheme. For the massive operator matrix elements $A_{qq,Q}^{(3),\rm PS}$ and $A_{qg,Q}^{(3),\rm S}$ the complete results are presented. The expressions are given in Mellin-$N$ space and in momentum fraction $z$-space., Comment: 86 pages Latex

An automated treatment of iterated integrals based on letters induced by real-valued quadratic forms and Kummer--Poincar\'e letters is presented. These quantities emerge in analytic single and multi--scale Feynman diagram calculations. To compactify representations, one wishes to apply general properties of these quantities in computer-algebraic implementations. We provide the reduction to basis representations, expansions, analytic continuation and numerical evaluation of these quantities., Comment: 14 pages LATEX, 1 anc. file

We calculate the polarized massive operator matrix element $A_{gq}^{(3)}(N)$ to 3-loop order in Quantum Chromodynamics analytically at general values of the Mellin variable $N$ both in the single- and double-mass case in the Larin scheme. It is a transition function required in the variable flavor number scheme at $O(\alpha_s^3)$. We also present the results in momentum fraction space., Comment: 26 pages Latex

We compute the two-mass contributions to the polarized massive operator matrix element $A_{gg,Q}^{(3)}$ at third order in the strong coupling constant $\alpha_s$ in Quantum Chromodynamics analytically. These corrections are important ingredients for the matching relations in the variable flavor number scheme and for the calculation of Wilson coefficients in deep--inelastic scattering in the asymptotic regime $Q^2 \gg m_c^2, m_b^2$. The analytic result is expressed in terms of nested harmonic, generalized harmonic, cyclotomic and binomial sums in $N$-space and by iterated integrals involving square-root valued arguments in $z$ space, as functions of the mass ratio. Numerical results are presented. New two--scale iterative integrals are calculated., Comment: 59 Latex, 2 figures

We calculate the massive polarized three-loop pure singlet operator matrix element $A_{Qq}^{(3), \rm PS}$ in the single mass case in the Larin scheme. This operator matrix element contributes to the massive polarized three-loop Wilson coefficient $H_{Qq}^{(3),\rm PS}$ in deep-inelastic scattering and constitutes a three-loop transition matrix element in the variable flavor number scheme. We provide analytic results in Mellin $N$ and in $x$ space and study the behaviour of this operator matrix element in the region of small and large values of the Bjorken variable $x$., Comment: 25 pages Latex, 2 Figures

We present the two-mass QCD contributions to the polarized pure singlet operator matrix element at three loop order in $x$-space. These terms are relevant for calculating the polarized structure function $g_1(x,Q^2)$ at $O(\alpha_s^3)$ as well as for the matching relations in the variable flavor number scheme and the polarized heavy quark distribution functions at the same order. The result for the operator matrix element is given in terms of generalized iterated integrals. These integrals depend on the mass ratio through the main argument, and the alphabet includes square--root valued letters., Comment: 21 pages LATEX, 3 Figures. arXiv admin note: substantial text overlap with arXiv:1711.06717

The post-Newtonian and post-Minkowskian solutions for the motion of binary mass systems in gravity can be derived in terms of momentum expansions within effective field theory approaches. In the post-Minkowskian approach the expansion is performed in the ratio $G_N/r$, retaining all velocity terms completely, while in the post-Newtonian approach only those velocity terms are accounted for which are of the same order as the potential terms due to the virial theorem. We show that it is possible to obtain the complete post-Minkowskian expressions completely algorithmically, under most general purely mathematical conditions from a finite number of velocity terms and illustrate this up to the third post-Minkowskian order given in \cite{Bern:2019crd}., Comment: 12 pages Latex

We calculate all contributions $\propto T_F$ to the polarized three-loop anomalous dimensions in the M-scheme using massive operator matrix elements and compare to results in the literature. This includes the complete anomalous dimensions $\gamma_{qq}^{(2),\rm PS}$ and $\gamma_{qg}^{(2)}$. We also obtain the complete two-loop polarized anomalous dimensions in an independent calculation. While for most of the anomalous dimensions the usual direct computation methods in Mellin $N$-space can be applied since all recurrences factorize at first order, this is not the case for $\gamma_{qg}^{(2)}$. Due to the necessity of deeper expansions of the master integrals in the dimensional parameter $\varepsilon = D-4$, we had to use the method of arbitrary high moments to eliminate elliptic contributions in intermediate steps. 4000 moments were generated to determine this anomalous dimension and 2640 moments turned out to be sufficient. As an aside, we also recalculate the contributions $\propto T_F$ to the three-loop QCD $\beta$-function., Comment: 40 pages Latex, 2 Figures

A summary of the calculation of the color-planar and complete light quark contributions to the massive three-loop form factors is presented. Here a novel calculation method for the Feynman integrals is used, solving general uni-variate first order factorizable systems of differential equations. We also present predictions for the asymptotic structure of these form factors., Comment: 8 pages; To appear in the contributions of DAE-BRNS 2018, 10th - 14th December 2018, Chennai, India

Higher order calculations in perturbative Quantum Field Theories often produce coupled linear systems of differential equations which factorize to first order. Here we present an algorithm to solve such systems in terms of iterated integrals over an alphabet the structure of which is implied by the coefficient matrix of the given system. We apply this method to calculate the master integrals in the color-planar and complete light quark contributions to the three-loop massive form factors., Comment: 9 pages; To appear in the proceedings of ACAT 2019, 10th - 15th March 2019, Saas Fee, Switzerland

We present an algorithm which allows to solve analytically linear systems of differential equations which factorize to first order. The solution is given in terms of iterated integrals over an alphabet where its structure is implied by the coefficient matrix of the differential equations. These systems appear in a large variety of higher order calculations in perturbative Quantum Field Theories. We apply this method to calculate the master integrals of the three--loop massive form factors for different currents, as an illustration, and present the results for the vector form factors in detail. Here the solution space emerging is given by the cyclotomic harmonic polylogarithms and their associated special constants. No special basis representation of the master integrals is needed. The algorithm can be applied as well to more general cases factorizing at first order, which are based on more general alphabets, iterated integrals and associated constants., Comment: 39 pages Latex, 19 figures, 2 style files, and anc-files

We present the FORTRAN-code HPOLY.f for the numerical calculation of harmonic polylogarithms up to w = 8 at an absolute accuracy of $\sim 4.9 \cdot 10^{-15}$ or better. Using algebraic and argument relations the numerical representation can be limited to the range $x \in [0, \sqrt{2}-1]$. We provide replacement files to map all harmonic polylogarithms to a basis and the usual range of arguments $x \in ]-\infty,+\infty[$ to the above interval analytically. We also briefly comment on a numerical implementation of real valued cyclotomic harmonic polylogarithms., Comment: 19 pages LATEX, 3 Figures, ancillary data

We describe the analytic calculation of the master integrals required to compute the two-mass three-loop corrections to the $\rho$ parameter. In particular, we present the calculation of the master integrals for which the corresponding differential equations do not factorize to first order. The homogeneous solutions to these differential equations are obtained in terms of hypergeometric functions at rational argument. These hypergeometric functions can further be mapped to complete elliptic integrals, and the inhomogeneous solutions are expressed in terms of a new class of integrals of combined iterative non-iterative nature., Comment: 14 pages Latex, 7 figures, to appear in the Proceedings of "Loops and Legs in Quantum Field Theory - LL 2018", 29 April - 4 May 2018, PoS

We calculate the two-mass QCD contributions to the massive operator matrix element $A_{gg,Q}$ at $\mathcal{O} (\alpha_s^3)$ in analytic form in Mellin $N$- and $z$-space, maintaining the complete dependence on the heavy quark mass ratio. These terms are important ingredients for the matching relations of the variable flavor number scheme in the presence of two heavy quark flavors, such as charm and bottom. In Mellin $N$-space the result is given in the form of nested harmonic, generalized harmonic, cyclotomic and binomial sums, with arguments depending on the mass ratio. The Mellin inversion of these quantities to $z$-space gives rise to generalized iterated integrals with square root valued letters in the alphabet, depending on the mass ratio as well. Numerical results are presented., Comment: 99 pages LATEX, 2 Figures

Various of the single scale quantities in massless and massive QCD up to 3-loop order can be expressed by iterative integrals over certain classes of alphabets, from the harmonic polylogarithms to root-valued alphabets. Examples are the anomalous dimensions to 3-loop order, the massless Wilson coefficients and also different massive operator matrix elements. Starting at 3-loop order, however, also other letters appear in the case of massive operator matrix elements, the so called iterative non-iterative integrals, which are related to solutions based on complete elliptic integrals or any other special function with an integral representation that is definite but not a Volterra-type integral. After outlining the formalism leading to iterative non-iterative integrals,we present examples for both of these cases with the 3-loop anomalous dimension $\gamma_{qg}^{(2)}$ and the structure of the principle solution in the iterative non-interative case of the 3-loop QCD corrections to the $\rho$-parameter., Comment: 13 pages LATEX, 2 Figures

We present the two-mass QCD contributions to the pure singlet operator matrix element at three loop order in x-space. These terms are relevant for calculating the structure function $F_2(x,Q^2)$ at $O(\alpha_s^3)$ as well as for the matching relations in the variable flavor number scheme and the heavy quark distribution functions at the same order. The result for the operator matrix element is given in terms of generalized iterated integrals that include square root letters in the alphabet, depending also on the mass ratio through the main argument. Numerical results are presented., Comment: 28 papges Latex, 3 figures

We calculate the two-mass contribution to the 3-loop vacuum polarization of the gluon in Quantum Chromodynamics at virtuality $p^2 = 0$ for general masses and also present the analogous result for the photon in Quantum Electrodynamics., Comment: 5 pages Latex

We calculate 3-loop master integrals for heavy quark correlators and the 3-loop QCD corrections to the $\rho$-parameter. They obey non-factorizing differential equations of second order with more than three singularities, which cannot be factorized in Mellin-$N$ space either. The solution of the homogeneous equations is possible in terms of convergent close integer power series as $_2F_1$ Gau\ss{} hypergeometric functions at rational argument. In some cases, integrals of this type can be mapped to complete elliptic integrals at rational argument. This class of functions appears to be the next one arising in the calculation of more complicated Feynman integrals following the harmonic polylogarithms, generalized polylogarithms, cyclotomic harmonic polylogarithms, square-root valued iterated integrals, and combinations thereof, which appear in simpler cases. The inhomogeneous solution of the corresponding differential equations can be given in terms of iterative integrals, where the new innermost letter itself is not an iterative integral. A new class of iterative integrals is introduced containing letters in which (multiple) definite integrals appear as factors. For the elliptic case, we also derive the solution in terms of integrals over modular functions and also modular forms, using $q$-product and series representations implied by Jacobi's $\vartheta_i$ functions and Dedekind's $\eta$-function. The corresponding representations can be traced back to polynomials out of Lambert--Eisenstein series, having representations also as elliptic polylogarithms, a $q$-factorial $1/\eta^k(\tau)$, logarithms and polylogarithms of $q$ and their $q$-integrals. Due to the specific form of the physical variable $x(q)$ for different processes, different representations do usually appear. Numerical results are also presented., Comment: 68 pages LATEX, 10 Figures

Starting at 3-loop order, the massive Wilson coefficients for deep-inelastic scattering and the massive operator matrix elements describing the variable flavor number scheme receive contributions of Feynman diagrams carrying quark lines with two different masses. In the case of the charm and bottom quarks, the usual decoupling of one heavy mass at a time no longer holds, since the ratio of the respective masses, $\eta = m_c^2/m_b^2 \sim 1/10$, is not small enough. Therefore, the usual variable flavor number scheme (VFNS) has to be generalized. The renormalization procedure in the two--mass case is different from the single mass case derived in \cite{Bierenbaum:2009mv}. We present the moments $N=2,4$ and $6$ for all contributing operator matrix elements, expanding in the ratio $\eta$. We calculate the analytic results for general values of the Mellin variable $N$ in the flavor non-singlet case, as well as for transversity and the matrix element $A_{gq}^{(3)}$. We also calculate the two-mass scalar integrals of all topologies contributing to the gluonic operator matrix element $A_{gg}$. As it turns out, the expansion in $\eta$ is usually inapplicable for general values of $N$. We therefore derive the result for general values of the mass ratio. From the single pole terms we derive, now in a two-mass calculation, the corresponding contributions to the 3-loop anomalous dimensions. We introduce a new general class of iterated integrals and study their relations and present special values. The corresponding functions are implemented in computer-algebraic form., Comment: 99 pages Latex, 12 figures, 2 style files

We calculate the unpolarized twist-2 three-loop splitting functions $P_{qg}^{(2)}(x)$ and $P_{gg}^{(2,\rm N_F)}(x)$ and the associated anomalous dimensions using massive three-loop operator matrix elements. While we calculate $P_{gg}^{(2,\rm N_F)}(x)$ directly, $P_{qg}^{(2)}(x)$ is computed from 1200 even moments, without any structural prejudice, using a hierarchy of recurrences obtained for the corresponding operator matrix element. The largest recurrence to be solved is of order 12 and degree 191. We confirm results in the foregoing literature., Comment: 38 pages Latex, 2 Figures

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 print

Three loop ladder and $V$-topology diagrams contributing to the massive operator matrix element $A_{Qg}$ are calculated. The corresponding objects can all be expressed in terms of nested sums and recurrences depending on the Mellin variable $N$ and the dimensional parameter $\varepsilon$. Given these representations, the desired Laurent series expansions in $\varepsilon$ can be obtained with the help of our computer algebra toolbox. Here we rely on generalized hypergeometric functions and Mellin-Barnes representations, on difference ring algorithms for symbolic summation, on an optimized version of the multivariate Almkvist-Zeilberger algorithm for symbolic integration, and on new methods to calculate Laurent series solutions of coupled systems of differential equations. The solutions can be computed for general coefficient matrices directly for any basis also performing the expansion in the dimensional parameter in case it is expressible in terms of indefinite nested product-sum expressions. This structural result is based on new results of our difference ring theory. In the cases discussed we deal with iterative sum- and integral-solutions over general alphabets. The final results are expressed in terms of special sums, forming quasi-shuffle algebras, such as nested harmonic sums, generalized harmonic sums, and nested binomially weighted (cyclotomic) sums. Analytic continuations to complex values of $N$ are possible through the recursion relations obeyed by these quantities and their analytic asymptotic expansions. The latter lead to a host of new constants beyond the multiple zeta values, the infinite generalized harmonic and cyclotomic sums in the case of $V$-topologies., Comment: 110 pages Latex, 4 Figures

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_{Qq}^{(3), \sf PS}$ 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., Comment: 85 pages Latex, 14 Figures, 2 style files

We consider finite iterated generalized harmonic sums weighted by the binomial $\binom{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 \rightarrow \infty$ 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 \rightarrow \infty$ 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 \in \mathbb{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., Comment: 62 pages Latex, 1 style file

The $O(\alpha_s^3 T_F^2 C_F (C_A))$ contributions to the transition matrix element $A_{gg,Q}$ relevant for the variable flavor number scheme at 3--loop order are calculated. The corresponding graphs contain two massive fermion lines of equal mass leading to terms given by inverse binomially weighted sums beyond the usual harmonic sums. In $x$-space two root-valued letters contribute in the iterated integrals in addition to those forming the harmonic polylogarithms. We outline technical details needed in the calculation of graphs of this type, which are as well of importance in the case of two different internal massive lines., Comment: 39 papges LATEX, 8 Figures

A survey is given on mathematical structures which emerge in multi-loop Feynman diagrams. These are multiply nested sums, and, associated to them by an inverse Mellin transform, specific iterated integrals. Both classes lead to sets of special numbers. Starting with harmonic sums and polylogarithms we discuss recent extensions of these quantities as cyclotomic, generalized (cyclotomic), and binomially weighted sums, associated iterated integrals and special constants and their relations., Comment: 10 pages, Proceedings ACAT 2013

Recent results on the calculation of 3-loop massive operator matrix elements in case of one and two heavy quark masses are reported. They concern the $O(n_f T_F^2 C_{F,A})$ and $O(T_F^2 C_{F,A})$ gluonic corrections, two-mass quarkonic moments, and ladder- and Benz-topologies. We also discuss technical aspects of the calculations., Comment: 12 pages, Proc. Loops and Legs in Quantum Field Theory - 11th DESY Workshop on Elementary Particle Physics, April 15-20, 2012, Wernigerode, Germany

Two-point Feynman parameter integrals, with at most one mass and containing local operator insertions in $4+\ep$-dimensional Minkowski space, can be transformed to multi-integrals or multi-sums over hyperexponential and/or hypergeometric functions depending on a discrete parameter $n$. Given such a specific representation, we utilize an enhanced version of the multivariate Almkvist--Zeilberger algorithm (for multi-integrals) and a common summation framework of the holonomic and difference field approach (for multi-sums) to calculate recurrence relations in $n$. Finally, solving the recurrence we can decide efficiently if the first coefficients of the Laurent series expansion of a given Feynman integral can be expressed in terms of indefinite nested sums and products; if yes, the all $n$ solution is returned in compact representations, i.e., no algebraic relations exist among the occurring sums and products., Comment: 14 pages, Proceedings, Loops and Legs in Quantum Field Theory 2012, Wernigerode,D; PoS(2012)

A big class of Feynman integrals, in particular, the coefficients of their Laurent series expansion w.r.t.\ the dimension parameter $\ep$ can be transformed to multi-sums over hypergeometric terms and harmonic sums. In this article, we present a general summation method based on difference fields that simplifies these multi--sums by transforming them from inside to outside to representations in terms of indefinite nested sums and products. In particular, we present techniques that assist in the task to simplify huge expressions of such multi-sums in a completely automatic fashion. The ideas are illustrated on new calculations coming from 3-loop topologies of gluonic massive operator matrix elements containing two fermion lines, which contribute to the transition matrix elements in the variable flavor scheme., Comment: 9 pages Latex, Contr. Proc. RADCOR 2011

The contributions $\propto n_f$ to the $O(\alpha_s^3)$ massive operator matrix elements describing the heavy flavor Wilson coefficients in the limit $Q^2 \gg m^2$ are computed for the structure function $F_2(x,Q^2)$ and transversity for general values of the Mellin variable $N$. Here, for two matrix elements, $A_{qq,Q}^{\sf PS}(N)$ and $A_{qg,Q}(N)$, the complete result is obtained. A first independent computation of the contributions to the 3--loop anomalous dimensions $\gamma_{qg}(N)$, $\gamma_{qq}^{\sf PS}(N$ and $\gamma_{qq}^{\sf NS,(TR)}(N)$ is given. In the computation advanced summation technologies for nested sums over products of hypergeometric terms with harmonic sums have been used. For intermediary results generalized harmonic sums occur, while the final results can be expressed by nested harmonic sums only., Comment: 29 pages, 1 style file

We report on results for the heavy flavor contributions to $F_2(x,Q^2)$ in the limit $Q^2\gg m^2$ at {\sf NNLO}. By calculating the massive $3$--loop operator matrix elements, we account for all but the power suppressed terms in $m^2/Q^2$. Recently, the calculation of fixed Mellin moments of all $3$--loop massive operator matrix elements has been finished. We present new all--$N$ results for the $O(n_f)$--terms, thereby confirming the corresponding parts of the $3$--loop anomalous dimensions. Additionally, we report on first genuine $3$--loop results of the ladder--type diagrams for general values of the Mellin variable $N$., Comment: 8 pages LATEX, 4 figures, 1 style file

By symbolic summation methods based on difference fields we present a general strategy that transforms definite multi-sums, e.g., in terms of hypergeometric terms and harmonic sums, to indefinite nested sums and products. We succeeded in this task with all our concrete calculations of 2--loop and 3--loop massive single scale Feynman diagrams with local operator insertion., Comment: 6 pages 2 figures