Calculation of massive 2-loop operator matrix elements with outer gluon lines

Massive on-shell operator matrix elements and self-energy diagrams with outer gluon lines are calculated analytically at O(αs2), using Mellin-Barnes integrals and representations through generalized hypergeometric functions. This method allows for a direct evaluation without decomposing the integrals using the integration-by-parts method.