Two-loop massive operator matrix elements for unpolarized heavy flavor production to O(ϵ)

We calculate the O(αs2) massive operator matrix elements for the twist-2 operators, which contribute to the heavy flavor Wilson coefficients in unpolarized deeply inelastic scattering in the region Q2≫m2, up to the O(ε) contributions. These terms contribute through the renormalization of the O(αs3) heavy flavor Wilson coefficients of the structure function F2(x,Q2). The calculation has been performed using light-cone expansion techniques without using the integration-by-parts method. We represent the individual Feynman diagrams by generalized hypergeometric structures, the ε-expansion of which leads to infinite sums depending on the Mellin variable N. These sums are finally expressed in terms of nested harmonic sums using the general summation techniques implemented in the Sigma package.