
We present the matching relations of the variable flavor number scheme at
nexttoleading 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 twomass 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 twomass 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 twoloop QCD corrections to the heavy quark form factors in
case of the vector, axialvector, scalar and pseudoscalar currents up to
second order in the dimensional parameter $\epsilon = (4D)/2$. These terms are
required in the renormalization of the higher order corrections to these form
factors.

We present the twomass QCD contributions to the pure singlet operator matrix
element at three loop order in xspace. 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 twoloop order in
perturbative QCD for different currents, namely vector, axialvector, scalar
and pseudoscalar 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 twomass contribution to the 3loop 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 3loop master integrals for heavy quark correlators and the
3loop QCD corrections to the $\rho$parameter. They obey nonfactorizing
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, squareroot 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 LambertEisenstein
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 3loop order, the massive Wilson coefficients for deepinelastic
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 twomass 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 nonsinglet case, as well as
for transversity and the matrix element $A_{gq}^{(3)}$. We also calculate the
twomass 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 twomass calculation, the corresponding contributions to the 3loop
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 computeralgebraic form.

We calculate the unpolarized twist2 threeloop splitting functions
$P_{qg}^{(2)}(x)$ and $P_{gg}^{(2,\rm N_F)}(x)$ and the associated anomalous
dimensions using massive threeloop 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 nonsinglet charged current deepinelastic 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
3loop 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 3loop contributions to
the heavy quark parts of the structure functions in deepinelastic scattering
due to charm and bottom.

A survey is given on the status of 3loop heavy flavor corrections to
deepinelastic 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 MellinBarnes representations, on
difference ring algorithms for symbolic summation, on an optimized version of
the multivariate AlmkvistZeilberger 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
productsum expressions. This structural result is based on new results of our
difference ring theory. In the cases discussed we deal with iterative sum and
integralsolutions over general alphabets. The final results are expressed in
terms of special sums, forming quasishuffle 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 nonsinglet charged current deepinelastic 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 3loop 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 GrossLlewellyn Smith sum rule are
discussed.

We calculate the massive flavor nonsinglet 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 3loop 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 sumrule. Numerical results are
presented for the charm quark contribution. Results on the structure function
$g_2(x,Q^2)$ in the twist2 approximation are also given.

The 3loop heavy flavor corrections to deepinelastic scattering are
essential for consistent nexttonexttoleading 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,
computeralgebraic 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 nexttonexttoleading order from the world
precision data on deepinelastic scattering.

We present our most recent results on the calculation of the heavy flavor
contributions to deepinelastic scattering at 3loop 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 nonsinglet 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 3loop 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 squareroot
valued letters.

The pure singlet asymptotic heavy flavor corrections to 3loop order for the
deepinelastic 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 3loop 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 threeloop heavy
flavor contributions to the Wilson coefficients in deepinelastic 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 nonsinglet
and pure singlet Wilson coefficients.

We consider gluonic contributions to the heavy flavor Wilson coefficients at
3loop 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 squareroot 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 nonsinglet 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 nonsinglet 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 3loop 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 3loop anomalous dimensions to $O(N_F)$ and
compare to results in the literature. The 3loop matching of the flavor
nonsinglet 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
nonsinglet 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 3loop
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 rootvalued 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 deepinelastic scattering in the asymptotic region $Q^2 \gg m^2$ to 3loop
order in the fixedflavor 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 3loop
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 3loop anomalous dimension
$\gamma_{gq}^{(2)}(N)$.