The Polarized Two-Loop Massive Pure Singlet Wilson Coefficient for Deep-Inelastic Scattering

We calculate the polarized massive two--loop pure singlet Wilson coefficient contributing to the structure functions $g_1(x,Q^2)$ analytically in the whole kinematic region. The Wilson coefficient contains Kummer--elliptic integrals. We derive the representation in the asymptotic region $Q^2 \gg m^2$, retaining power corrections, and in the threshold region. The massless Wilson coefficient is recalculated. The corresponding twist--2 corrections to the structure function $g_2(x,Q^2)$ are obtained by the Wandzura--Wilczek relation. Numerical results are presented.


Introduction
Analytic expressions in perturbative Quantum Chromodynamics, as for partonic sub-system scattering cross sections and Wilson coefficients, provide an excellent basis for both numerical studies, fits to precision data and provide analytic insight into the structure of Feynman integrals. Analytic results also allow to derive important limiting cases and make it easier to incorporate resummations in specific kinematic regions. With a growing number of scales and loops for the respective processes it becomes more difficult to obtain analytic results, although there has been significant progress in analytic integration methods recently; for a survey cf. [1]. The precise knowledge of these corrections is of importance to measure the polarized parton densities in high energy collisions and to determine, related to it, the strong coupling constant α s (M Z ) and the heavy quark masses, cf. [2].
The first two-loop QCD heavy flavor corrections to the polarized structure function g 1 (x, Q 2 ) have been calculated in [3] in the asymptotic region Q 2 m 2 , where Q 2 denotes the virtuality of the exchanged photon and m the mass of the heavy quark. The asymptotic two-loop QCD corrections have been recalculated in [4,5]. In [3] the region of low values of Q 2 has been modeled by an ansatz. The leading threshold resummation for the gluonic contributions has been studied in [6]. The complete two-loop polarized heavy flavor Wilson coefficient in the non-singlet case has been calculated analytically in the tagged flavor case in [3] and for the complete contribution to the structure function g 1 (x, Q 2 ) in [7], also completing former work on the polarized Bjorken sum rule in [8]. Numerical results for the polarized two-loop heavy flavor case have been given in [9] recently. Finally, in the non-singlet case the asymptotic contributions have been calculated to three-loop order analytically in [10,11].
In the present paper we follow Ref. [12] in the unpolarized case and calculate the polarized pure singlet two-loop heavy flavor corrections for the structure function g 1 (x, Q 2 ) in the whole kinematic range analytically. We also compute the corresponding massless contributions, which have first been calculated in [13] and later in [14]. Since the calculations are carried out using dimensional regularization in D = 4 + ε dimensions one may work in the Larin-scheme [15] 1 and perform a finite renormalization to the M-scheme [17,18] afterwards. We derive both the result in the asymptotic range Q 2 m 2 , see also Refs. [3,5], and in the threshold region. Numerical results are presented. Various technical aspects of the present calculation can be found in Ref. [12] already.
The paper is organized as follows. In Section 2 we summarize basic relations for the polarized deep-inelastic scattering cross section. In Section 3 the result for the massless pure singlet Wilson coefficient C (2),PS g 1 is presented. The recalculation of the massless Wilson coefficient is necessary, since in Ref. [13] different schemes have been used in part. The corresponding massive Wilson coefficient is calculated in Section 4. The corresponding results for the twist-2 contributions to the structure function g 2 (x, Q 2 ) can be obtained by using the Wandzura-Wilczek relation [19], as has been shown for the massless quarkonic [20,21] and gluonic [22] cases, for diffractive scattering [23], non-forward scattering [24], and the target mass corrections [25,27]. Limiting cases are studied in Section 5 and numerical results are presented in Section 6. The conclusions are given in Section 7. Some Mellin convolutions appearing due to renormalization are listed in the Appendix.

The Deep-inelastic Scattering Cross Section
The scattering cross sections for deep-inelastic charged lepton scattering of polarized nucleons 1 For other γ 5 schemes see Refs. [16].
are obtained polarizing the incoming lepton longitudinally and the target nucleon either longitudinally or transversally, resulting into the spin 4-vectors S L and S T , (1) S T = M (0, cos(β), sin(β); 0) (2) in the nucleon rest frame. One has S L .p = S T .p = 0, with p the nucleon 4-momentum. The scattering cross sections are given by, cf. e.g. [25,26], for pure virtual photon exchange. Here S denotes the cms energy of the process, M is the nucleon mass, λ the degree of lepton polarization, α = e 2 /(4π) is the fine structure constant, Q 2 = −q 2 denotes the photon virtuality and x = Q 2 /(Sy), y = l.q/p.q are the Bjorken variables with l the incoming charged lepton and proton momenta, S = (p + l) 2 and φ is the azimuthal angle of the final state lepton, which can be integrated over in the case of longitudinal polarization. In the following we will present a series of relations in Mellin-N space for convenience. The respective quantities in momentum-fraction z-space are related to those in Mellin-space by the transformation The structure function g 1 (N, Q 2 ) is given in the twist-two approximation using the factorization theorems [28] by Here denotes the singlet distribution, G(N ) the gluon distribution, ∆(N ) the flavor non-singlet distribution e k labels the electric charge of the kth light quark.
The Mellin transform of the structure function g 2 is related to that of g 1 by Note that in the massive pure singlet case the support of both structure functions in limited by 0 < x < 1/(1 + 4m 2 /Q 2 ) due to the production of two heavy quarks. The different steps in the renormalization and factorization of the polarized massless Wilson coefficients have been described in [13,29] and for the massive Wilson coefficients in [5] using the Larin-scheme. In the present case the finite renormalization moving to the M-scheme only affects the massless or massive Wilson coefficient by adding the term −z

The Massless Wilson Coefficient
The Feynman diagrams contributing to the polarized massless two-loop Wilson coefficient are shown in Figure 1. forming the amplitude squared G µν , where the indices refer to the coupling of the virtual photon. Here all quark lines are massless. The massless resp. massive Wilson coefficients are obtained following Ref. [13], Eqs. (3.7-3.18). The corresponding phase space integrals have been given in Ref. [12]. We apply the Larin-scheme [15] in which the contraction of the free indices of the two appearing Levi-Civita tensors have to be performed in D dimensions.
The unrenormalized two-loop massless pure singlet Wilson coefficient reads in Mellin-N spacê g 1 ,g +c (2),PS whereâ s =ĝ 2 s /(4π) 2 denotes the unrenormalized strong coupling constant, S ε = exp[ ε 2 (γ E − ln(4π))] the spherical factor and γ E the Euler-Mascheroni constant,c (k) i and a (1) g 1 ,g are the expansion coefficients of the one-loop Wilson coefficient witĥ given by the Feynman diagrams in Figure 2, where all quark lines are massless. One obtains The contributing splitting functions [30][31][32][33][34] are where C F = (N 2 c − 1)/2N c , T F = 1/2 and N c = 3 for the gauge group SU (3) c of Quantum Chromodynamics. In the following we will use the convention The harmonic polylogarithms are given by [35] and the letters f c read We use the shorthand notation H a (z) ≡ H a in case of the argument z. The harmonic polylogarithms are dual to the harmonic sums [36,37] by the Mellin transformation (5).
In the Larin-scheme we obtain performing the phase space integrations as has been outlined in [12,29].
At O(a 2 s ) neither the renormalization of the heavy quark mass nor the coupling constant contributes in case of the massive or massless pure singlet Wilson coefficient. The poles in ε in Eq. (21) are due to collinear singularities only, which have to be factorized. One may proceed as follows. The unfactorized quarkonic Wilson coefficients for the structure function g 1 ,Ĉ NS,S 1,q in Mellin-space are given byĈ The pure singlet contribution is obtained bŷ with The factorized massless pure singlet two-loop Wilson coefficient C is given by where Here we set µ F = µ, and work with a single scale only for the factorization and renormalization scale; a s = g 2 s /(4π) is the running coupling constant and the spherical factor S ε , as usually, is set to one at the end of the calculation. Note that the splitting function P (1),PS qq is correctly obtained, cf. [33,34], despite working in the Larin-scheme, cf. also [18].
The massless Wilson coefficient in the M-scheme is obtained by the following finite renormalization with [17] z (2) cf. [5,18]. C in z-space is given by

The Massive Wilson Coefficient
The kinematic domain for the massive Wilson coefficient is given by The unrenormalized two-loop massive pure singlet Wilson coefficient reads in Mellin-N spacê The contributing Feynman diagrams are shown in Figure 1, where now the outgoing quark lines with momenta k 1 and k 2 are taken massive, see also Figure 2. Hereh given by the diagrams in Figure 2 now with massive quark lines. The expansion coefficients are given by Here β denotes the velocity of the produced heavy quarks, Since the two heavy quarks do not induce collinear divergences the mass factorization in the massive case readsĤ Here the Mellin convolution of two functions in z-space is given by We find Identifying the renormalization and factorization scale, µ = µ F , we finally obtain Note that in the pure singlet case neither the heavy quark mass nor the coupling constant is renormalized at two-loop order. The massive pure singlet Wilson coefficient is obtained as a four-fold integral over two angular and two energy variables, cf. [12] for details of the calculation. These integrals are systematically turned into iterative integrals. This process leads to a set of letters, through which these integrals are defined, see also [41]. We use the code [42] in Mathematica, which also uses the routine DSolveRational of the package HolonomicFunctions [43]; see [44,45] for the general theory underlying [42]. We also refer to [46] for the simpler case when no singularities are present at the endpoints of integration, which, however, does not apply here.

The Asymptotic and Threshold Expansions
The complete expressions calculated in Section 4 allow now to perform the asymptotic expansion for Q 2 m 2 and the threshold expansion for β 1.
In the asymptotic limit Q 2 m 2 the first expansion coefficients of the polarized massive pure singlet Wilson coefficient read setting µ 2 = Q 2 H (2),PS,M with the polynomials In this expansion the Kummer-elliptic integrals turn into harmonic polylogarithms. The leading term, which is free of power corrections of O((m 2 /Q 2 ) k ), k ∈ N, k ≥ 1, can be predicted using the representation of the massive Wilson coefficient by massive operator matrix elements (OMEs), cf. [3,5,47,48], and massless Wilson coefficients, Here the massless Wilson coefficientC cf. [3,5,47,48]; for its renormalization see Ref. [49]. The constant part of the unrenormalized polarized OME a (2),PS Qq is given by [3,5] a (2),PS Qq in z-space. The calculation of A (2),PS Qq is performed in the Larin-scheme. One has to apply the tensor decomposition method, however, to obtain the correct result. These aspects are discussed in Ref. [5] in detail.
The threshold expansion of the Wilson coefficients for β 1 is given by

Numerical Results
Let us now illustrate the analytic results numerically. In Figure 3 the two-loop heavy flavor Wilson coefficient H    For large values of Q 2 these results approach the asymptotic result for H 2,PS g 1 ,q . In the small x region this Wilson coefficient is negative.
Next we study the ratios comparing the full (57) and the asymptotic results,H, (74) for the leading term in Figure 4.
For H 2,PS g 1 ,q the asymptotic expansion agrees with the full calculation up to Q 2 /m 2 ≡ χ = 10 to about 2% for z = 10 −4 , χ = 40 for z = 10 −2 and χ = 200 for z = 1/2. However, the Wilson coefficients are very small already for the last value. Similar to the ratio of the full and asymptotic Wilson coefficient we define the ratio whereg (2),PS 1,q is the structure function obtained by using the expansion of the respective Wilson coefficient up to the desired level. The corresponding results are depicted in Figure 5. We use the parameterization of the parton distributions Ref. [50] at NLO with the corresponding values of α s (Q 2 ) at NNLO [51] to compare to previous non-singlet results in [10]. Demanding an agreement within ±2% for g     [50] . In Figure 7 we illustrate the ratios Eq. (85) as a function of x for different values of Q 2 for g PS 1 comparing the asymptotic result to the full result. For a better visibility and to avoid to depict zero transitions in the denominator we separate the small x and large x part into two plots. The corrections behave widely flat in x for larger values of Q 2 and develop some profile for Q 2 < 100 GeV 2 .
In Figure 8 we depict the ratio of the full result over the O((m 2 /Q 2 ) 2 ) improved asymptotic results for g PS 1 as a function of x for a series of Q 2 -values, again separating the small x and the large x ranges because of zero transitions for this ratio. For Q 2 > ∼ 100GeV 2 the ratios are rather flat and are close to one. The line for Q 2 = 100GeV 2 for x > 0.5 deviates from one by more than 5%. Larger deviations are found for Q 2 = 50GeV 2 , where the 5% marging is only met for x < 3 · 10 −3 . We limited the expansion to terms of ∼ O((m 2 /Q 2 ) 2 ) in the present paper, but higher order terms can be given straigtforwardly. The expanded expressions do also allow direct Mellin transforms and provide a suitable analytic basis for Mellin-space programmes.

Conclusions
We have calculated the massless and massive polarized two-loop pure singlet Wilson coefficients for deep-inelastic scattering in analytic form. The calculation has been performed in the Larinscheme, with a final finite renormalization to the M-scheme, which has been introduced in Refs. [17,18]. The massless Wilson coefficient can be expressed by the harmonic polylogarithms in z-space and harmonic sums in Mellin-N space. In the massive case the polarized two-loop pure singlet Wilson coefficient is also given by iterative integrals, however, of a more general kind, the Kummer-elliptic integrals, here based on an alphabet of 12 letters, cf. [12]. From the expansion of the massive Wilson coefficient in the region Q 2 m 2 one obtains the asymptotic result, which can be given in terms of a massive OME and the massless Wilson coefficient, cf. [5]. In the region of lower values of Q 2 and larger values of x, the power corrections to the massive two-loop Wilson coefficient are essential. From the available analytic result one can construct the series in m 2 /Q 2 analytically. Since the deep-inelastic process is usually considered only for virtualities Q 2 > ∼ 5 GeV 2 , this series gives the proper numerical representation in case of the charm-quark corrections retaining a relatively small number of terms. The latter representation has the advantage that it can be transformed into Mellin space directly, since the expansion coefficients are given in terms of harmonic polylogarithms in z-space.

A Contributing Expressions due to Renormalization
In the following we list the Mellin-convolutions, which occurred in Eq. (57). These are convolutions with leading order splitting functions, referring to the parameter κ = m 2 /Q 2 .