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

We calculate the O(α2s) 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(α3s) 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. (orig.)