• We present the matching relations of the variable flavor number scheme at next-to-leading order, which are of importance to define heavy quark partonic distributions for the use at high energy colliders such as Tevatron and the LHC. The consideration of the two-mass effects due to both charm and bottom quarks, having rather similar masses, are important. These effects have not been considered in previous investigations. Numerical results are presented for a wide range of scales. We also present the corresponding contributions to the structure function $F_2(x,Q^2)$.
  • 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.
  • We compute the two-loop QCD corrections to the heavy quark form factors in case of the vector, axial-vector, scalar and pseudo-scalar currents up to second order in the dimensional parameter $\epsilon = (4-D)/2$. These terms are required in the renormalization of the higher order corrections to these form factors.
  • 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.
  • We present the results for heavy quark form factors at two-loop order in perturbative QCD for different currents, namely vector, axial-vector, scalar and pseudo-scalar currents, up to second order in the dimensional regularization parameter. We outline the necessary computational details, ultraviolet renormalization and corresponding universal infrared structure.
  • 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.
  • 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.
  • 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.
  • 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.
  • We calculate the massive Wilson coefficients for the heavy flavor contributions to the non-singlet charged current deep-inelastic scattering structure functions $F_L^{W^+}(x,Q^2)-F_L^{W^-}(x,Q^2)$ and $F_2^{W^+}(x,Q^2)-F_2^{W^-}(x,Q^2)$ in the asymptotic region $Q^2 \gg m^2$ to 3-loop order in Quantum Chromodynamics (QCD) at general values of the Mellin variable $N$ and the momentum fraction $x$. Besides the heavy quark pair production, also the single heavy flavor excitation $s \rightarrow c$ contributes. Numerical results are presented for the charm quark contributions and consequences on the unpolarized Bjorken sum rule and Adler sum rule are discussed.
  • 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.
  • A survey is given on the status of 3-loop heavy flavor corrections to deep-inelastic structure functions at large enough virtualities $Q^2$.
  • 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.
  • We calculate the massive Wilson coefficients for the heavy flavor contributions to the non-singlet charged current deep-inelastic scattering structure function $xF_3^{W^+}(x,Q^2)+xF_3^{W^-}(x,Q^2)$ in the asymptotic region $Q^2 \gg m^2$ to 3-loop order in Quantum Chromodynamics (QCD) at general values of the Mellin variable $N$ and the momentum fraction $x$. Besides the heavy quark pair production also the single heavy flavor excitation $s \rightarrow c$ contributes. Numerical results are presented for the charm quark contributions and consequences on the Gross-Llewellyn Smith sum rule are discussed.
  • We calculate the massive flavor non-singlet Wilson coefficient for the heavy flavor contributions to the polarized structure function $g_1(x,Q^2)$ in the asymptotic region $Q^2 \gg m^2$ to 3-loop order in Quantum Chromodynamics at general values of the Mellin variable $N$ and the momentum fraction $x$, and derive heavy flavor corrections to the Bjorken sum-rule. Numerical results are presented for the charm quark contribution. Results on the structure function $g_2(x,Q^2)$ in the twist-2 approximation are also given.
  • The 3-loop heavy flavor corrections to deep-inelastic scattering are essential for consistent next-to-next-to-leading order QCD analyses. We report on the present status of the calculation of these corrections at large virtualities $Q^2$. We also describe a series of mathematical, computer-algebraic and combinatorial methods and special function spaces, needed to perform these calculations. Finally, we briefly discuss the status of measuring $\alpha_s(M_Z)$, the charm quark mass $m_c$, and the parton distribution functions at next-to-next-to-leading order from the world precision data on deep-inelastic scattering.
  • We present our most recent results on the calculation of the heavy flavor contributions to deep-inelastic scattering at 3-loop order in the large $Q^2$ limit, where the heavy flavor Wilson coefficients are known to factorize into light flavor Wilson coefficients and massive operator matrix elements. We describe the different techniques employed for the calculation and show the results in the case of the heavy flavor non-singlet and pure singlet contributions to the structure function $F_2(x,Q^2)$.
  • Contributions to heavy flavour transition matrix elements in the variable flavour number scheme are considered at 3-loop order. In particular a calculation of the diagrams with two equal masses that contribute to the massive operator matrix element $A_{gg,Q}^{(3)}$ is performed. In the Mellin space result one finds finite nested binomial sums. In $x$-space these sums correspond to iterated integrals over an alphabet containing also square-root valued letters.
  • 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.
  • We report on our latest results in the calculation of the three-loop heavy flavor contributions to the Wilson coefficients in deep-inelastic scattering in the asymptotic region $Q^2 \gg m^2$. We discuss the different methods used to compute the required operator matrix elements and the corresponding Feynman integrals. These methods very recently allowed us to obtain a series of new operator matrix elements and Wilson coefficients like the flavor non-singlet and pure singlet Wilson coefficients.
  • We consider gluonic contributions to the heavy flavor Wilson coefficients at 3-loop order in QCD with two heavy quark lines in the asymptotic region $Q^2 \gg m_{1(2)}^2$. Here we report on the complete result in the case of two equal masses $m_1 = m_2$ for the massive operator matrix element $A_{gg,Q}^{(3)}$, which contributes to the corresponding heavy flavor transition matrix element in the variable flavor number scheme. Nested finite binomial sums and iterated integrals over square-root valued alphabets emerge in the result for this quantity in $N$ and $x$-space, respectively. We also present results for the case of two unequal masses for the flavor non-singlet OMEs and on the scalar integrals ic case of $A_{gg,Q}^{(3)}$, which were calculated without a further approximation. The graphs can be expressed by finite nested binomial sums over generalized harmonic sums, the alphabet of which contains rational letters in the ratio $\eta = m_1^2/m_2^2$.
  • We calculate the massive flavor non-singlet Wilson coefficient for the heavy flavor contributions to the structure function $F_2(x,Q^2)$ in the asymptotic region $Q^2 \gg m^2$ and the associated operator matrix element $A_{qq,Q}^{(3), \rm NS}(N)$ to 3-loop order in Quantum Chromodynamics at general values of the Mellin variable $N$. This matrix element is associated to the vector current and axial vector current for the even and the odd moments $N$, respectively. We also calculate the corresponding operator matrix elements for transversity, compute the contributions to the 3-loop anomalous dimensions to $O(N_F)$ and compare to results in the literature. The 3-loop matching of the flavor non-singlet distribution in the variable flavor number scheme is derived. All results can be expressed in terms of nested harmonic sums in $N$ space and harmonic polylogarithms in $x$-space. Numerical results are presented for the non-singlet charm quark contribution to $F_2(x,Q^2)$.
  • 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.
  • We calculate the logarithmic contributions to the 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 massive operator matrix elements needed in the variable flavor number scheme. Explicit expressions are given both in Mellin-$N$ space and $z$-space.
  • We calculate the massive operator matrix element $A_{gq}^{(3)}(N)$ to 3-loop order in Quantum Chromodynamics at general values of the Mellin variable $N$. This is the first complete transition function needed in the variable flavor number scheme obtained at $O(\alpha_s^3)$. A first independent recalculation is performed for the contributions $\propto N_F$ of the 3-loop anomalous dimension $\gamma_{gq}^{(2)}(N)$.