# Longitudinally Polarized Photoproduction of Heavy Flavors at Next-to-Leading Order of QCD

## Abstract

We present a phenomenological study of charm quark photoproduction in longitudinally polarized lepton-hadron collisions at next-to-leading order accuracy of QCD. Our results are based on a recently developed, flexible parton-level Monte Carlo program for spin-dependent heavy flavor hadroproduction, which we extend to deal also with both direct and resolved photon contributions. The subsequent hadronization into charmed mesons is modeled in our calculations, which allows us to compare with data on double-spin asymmetries for meson production taken by the COMPASS collaboration. In general, next-to-leading order QCD corrections are found to be very significant and do not cancel in spin asymmetries. We elucidate the role of the individual hard scattering subprocesses and determine the range of parton momentum fractions predominantly probed for charm production at COMPASS. Theoretical uncertainties are estimated by varying renormalization and factorization scales and parameters controlling the hadronization of the charm quarks.

###### pacs:

13.88.+e and 12.38.Bx and 13.85.Ni## 1 Motivation and Introduction

The quest to understand the partonic structure of the nucleon spin remains to be one of the key research areas in Hadronic Physics even after more than 25 years of strenuous experimental and theoretical efforts. In particular, the polarized gluon density, defined as

(1) |

is still one of the most elusive quantities associated with the non-perturbative partonic structure of hadrons. Here, () denotes the probability of finding a gluon at a scale with light-cone momentum fraction and helicity () in a proton with helicity . The total, integrated gluon polarization,

(2) |

and a similar contribution from the sum of all quarks and antiquarks, enters the helicity sum rule of the nucleon along with the orbital angular momenta of quarks and gluons (1). The challenge is to precisely map in a wide range of in order to minimize extrapolation uncertainties in the first moment (2).

Currently, the best constraints on are derived from global QCD analyses (2); (3) which treat all available experimental probes simultaneously and consistently at a given order in the strong coupling in perturbative QCD. The availability of next-to-leading order (NLO) QCD corrections is essential for any meaningful, quantitative analysis of parton densities. In particular, single-inclusive pion (4) and jet production (5), measured in spin-dependent proton-proton collisions at BNL-RHIC (6), have started to put significant limits on the amount of gluon polarization in the nucleon (2). New, preliminary single and di-jet data from the STAR collaboration (7) show for the first time tantalizing hints for a non-zero (8); (9). Due to the given kinematics, the current probes mainly constrain in the medium-to-large region, , which is not sufficient to reliably determine its integral (2). A very significant contribution of up to one unit of , i.e., twice the proton spin, can still come from the unexplored small region (2). Narrowing down the uncertainties on and, at the same time, extending the range in continues to be the main objective of experimental efforts in the years to come. Eventually, only a future high-energy polarized electron-proton collider, such as the EIC project (10), will finally be able to quantitatively address all the remaining open questions related to the helicity structure in the small region (11).

World-data on polarized inclusive and semi-inclusive deep-inelastic scattering are pivotal in constraining the helicity quark and antiquark densities (2); (3); (8), but due to lack of sufficient kinematic coverage of present fixed-target experiments, information on the gluon density is notoriously difficult to obtain from QCD scaling violations. Viable probes of at fixed-target energies comprise of one- and two-hadron and open charm production and have been exploited by several experimental collaborations (12); (13); (14); (15). The proper theoretical description of these processes depends on the virtuality of the probing photon and is in general more involved than corresponding calculations for hadron-hadron scattering. In case of photoproduction, where a quasi-real photon is exchanged, one has to include also “resolved” contributions, where the photon fluctuates into a vector meson of the same quantum numbers before the hard scattering with partons in the proton takes place. If the virtuality of the photon is of or higher, resolved processes are sufficiently suppressed but the additional momentum scale greatly complicates the calculations of phase-space and loop integrals.

Not surprisingly, only very few calculations at NLO of QCD are available in case of polarized beams and targets (16); (17); (18), none of which for large virtualities . As a consequence, the available data sets on hadron and open charm production (12); (13); (14); (15) have not been included in global QCD analyses of helicity parton densities so far. Experiments have analyzed their data only in terms of the gluon polarization, , under certain simplifying assumptions and based on leading order (LO) matrix elements. Nonetheless, the results of these exercises, illustrated in Fig. 1, are in fairly good agreement with a NLO extraction of including the RHIC data (2). While determinations of at fixed-target energies are more involved and provide less of a constraint than collider data, they are crucial for further testing and establishing the assumed universality of helicity-dependent parton densities and hence for our understanding of the spin structure of the nucleon and QCD in general. On the one hand, lepton-nucleon scattering experiments are sensitive to different partonic hard scattering processes than jet or hadron production at RHIC and, on the other hand, the relevant momentum fractions fall within the range already probed by data.

In this paper, we present a comprehensive phenomenological analysis of longitudinally polarized photoproduction of heavy flavors at NLO of QCD. To this end, we extend our recently developed flexible parton-level Monte Carlo program for spin-dependent heavy flavor hadroproduction (19) by including all relevant subprocesses for direct (16); (17) and resolved (19); (20) photon contributions. To facilitate the comparison to data, the hadronization of the produced charm quarks into charmed mesons and, as an additional option, their semi-leptonic decays can be modeled in our calculations based on the phenomenological functions used in Ref. (21). The presented results will allow one to consistently include available and upcoming data on spin-dependent open charm photoproduction into future global QCD analyses of helicity parton densities at NLO accuracy.

As a first application, we examine in detail all aspects of open charm production in the kinematic regime accessible with the COMPASS experiment (15). We will demonstrate that NLO corrections are very significant and do not cancel in double-spin asymmetries used to extract information on the polarized gluon density . We elucidate the relevance of the individual direct and resolved partonic subprocesses and estimate theoretical uncertainties by varying renormalization and factorization scales and parameters controlling the hadronization of the produced charm quarks into experimentally observed mesons. Finally, we shall illustrate the sensitivity of the existing data on the double-spin asymmetry for charm quark photoproduction (15); (22); (23); (24) to and compare our estimates of the relevant range of momentum fractions with those obtained by the COMPASS collaboration.

We note that our flexible Monte Carlo code is capable of computing any infrared safe heavy flavor photoproduction cross section at in longitudinally polarized lepton-nucleon collisions, including correlations of the produced heavy quark pair. Our results complement and significantly extend previously existing spin-dependent NLO calculations of single-inclusive heavy quark yields based on largely analytical methods (16); (17), where the resolved contribution was neglected, any information on the partonic recoil system was lost, and most experimental cuts could not be implemented. The results presented here will be also useful for future studies of spin-dependent charm photoproduction at an EIC (10).

The outline of the paper is as follows: in Sec. II we briefly review some of the technical aspects of setting up a parton-level Monte Carlo program for heavy flavor photoproduction in polarized lepton-nucleon collisions at NLO accuracy. In Sec. III we present a detailed phenomenological study of open charm quark production at COMPASS, including the relevance of NLO corrections, the role of the different hard scattering subprocesses, an estimate of the relevant momentum fractions , and an assessment of theoretical uncertainties. We summarize our results in Sec. IV.

## 2 Technical Framework

We consider heavy quark photoproduction in longitudinally polarized lepton-nucleon collisions at NLO accuracy of QCD by consistently including for the first time both direct and resolved photon contributions. In order to compute arbitrary infrared-safe observables within flexible experimental acceptance cuts and to account for the hadronization of the produced heavy quark pair, all phase-space integrations are performed numerically with appropriate Monte Carlo techniques. To this end, we follow closely the subtraction method devised and used in Refs. (25); (26). In the following, we will only briefly outline the technical details relevant for the case of polarized photoproduction.

Assuming, as usual, factorization, the inclusive cross section for producing a heavy quark in spin-dependent lepton-proton collisions at a center-of-mass system (c.m.s.) energy can be written as a convolution,

(3) | |||||

where the subscripts in (3) label the helicity states of the lepton and proton. In analogy to Eq. (1), the denote the spin-dependent parton distribution functions (PDFs) of flavor in a lepton or proton . The sum in (3) is over all contributing partonic processes , including those with a direct photon, i.e., , with the associated polarized hard scattering cross sections. The required spin-dependent matrix elements squared at NLO accuracy in dimensional regularization for the direct and resolved photon processes can be taken from (16) and (20), respectively. Corresponding unpolarized results can be found in Refs. (25); (26); (27). At , parton can be either a gluon or a light (anti-)quark producing an additional jet or hadron which is usually not observed in experiment. denote the four-momenta of the heavy quark and antiquark with mass , i.e., . and are the renormalization and factorization scale, respectively, which are typically chosen as a combination of the hard scales characterizing the process.

The “measurement function” in Eq. (3) defines the observable, for instance, through a set of step functions implementing the experimental cuts imposed on the final-state particles and selecting a certain bin in a kinematical distribution of interest. In case of the COMPASS experiment, single-inclusive spectra of mesons, differential in transverse momentum , are measured (15); (22); (23); (24). Thus, the cross section (3) at the heavy quark-level is not yet sufficient for comparing theory with experimental results and needs to be convoluted with an additional phenomenological function modeling the hadronization of the charm quark:

(4) |

We will specify our choice for , along with all other non-perturbative inputs, in Sec. 3. Even though not relevant for the phenomenological applications considered in this paper, our flexible Monte Carlo program can also account for the semi-leptonic decays of mesons if needed.

For the integration of the fully exclusive partonic cross sections in Eq. (3), we generalize the framework used for unpolarized heavy quark production in (25); (26) to deal with singular regions of phase-space. The gist of the method is to add and subtract appropriate spin-dependent counter terms to which can be integrated analytically with respect to momenta of unresolved partons. For all resolved photon processes, , we can adopt the expressions obtained in our recent calculation of polarized heavy quark hadroproduction at (19), since these processes share the same hard scattering matrix elements. Hence, we only need to consider processes with a direct photon, , in the following.

In general, for a numerically efficient implementation of the subtraction method it is convenient (25); (26); (19) to express the three-body phase-space in terms of variables where soft and collinear singularities are particularly transparent. This is achieved by choosing , the invariant mass of the pair scaled by the available partonic c.m.s. energy squared, i.e., . In addition, one uses , the cosine of the angle between the -direction, aligned with the spatial direction of parton , and , the momentum of parton , in the c.m.s. of the incoming partons. Soft and collinear regions of phase-space are then simply associated with and , respectively, and the hard scattering matrix element for the process can be written as

(5) |

where is regular for and . The angles are used to parametrize the spatial orientation of with respect to the plane span by the other three momenta in the c.m.s. of the pair; for further details, see (25); (26); (19). We note that the genuine NLO subprocess can have only collinear singularities at .

The partonic subprocesses contributing to the direct photon cross section at , and , can be decomposed as (26)

(6) | |||||

Here, and denote the Born contribution and the one-loop corrections to the process, respectively. Analytic expressions for the virtual contributions in dimensions, with ultraviolet divergences being subtracted at the renormalization scale , can be found in Ref. (16).

In Eq. (6), denotes the soft gluon emission part of the scattering cross section, which can be obtained from the full -dimensional matrix elements squared in the limit where the phase-space integrations can be performed analytically (26). We obtain the same result for the squared three-body amplitude in the soft limit as given in Eq. (A.14) of Ref. (26) but with the Born contribution for being replaced by its spin-dependent, color-averaged counterpart, which reads in dimensions:

(7) |

where

(8) |

Here, and are the usual tree-level Mandelstam variables. Contrary to the unpolarized case, receives no contributions (16). is the electromagnetic charge of the heavy quark in units of the coupling , i.e., , and .

All processes exhibit singularities related to collinear splittings off the incoming photon () and parton () which need to be factorized into the bare photon and proton PDFs at a scale , respectively. As for the soft contributions discussed above, the kinematics collapses to the much simpler case of scattering, such that compact analytical expressions, summarized by in Eq. (6), can be obtained:

(9) | |||||

where

(10) | |||||

(11) | |||||

(12) | |||||

(13) | |||||

Here, and denote the momenta of the photon and parton , respectively. Both, the standard two-body phase-space in (9) and the Born matrix elements squared in Eqs. (10)-(13) are to be evaluated with appropriate collinear kinematics as indicated by the shift of momenta , etc; further details can be found in Refs. (25); (26); (19).

The resulting divergence in (9) assumes the form dictated by the factorization theorem, i.e., a convolution of -dimensional helicity-dependent LO splitting functions and Born matrix elements . By adding appropriate counter cross sections to (9), which to schematically read

(14) | |||||

where

(15) |

all collinear singularities can be consistently factorized into the scale evolution of the bare photon and proton PDFs, depending on whether they originate from collinear configurations involving the initial-state photon or parton . The required Born cross sections and LO are listed in the Appendix of Ref. (19), except for

(16) |

in (16) and, below, are the usual QCD color factors

The factorization scheme is fully specified by the choice of in (15) for which we take , to guarantee helicity conservation within the HVBM prescription for in dimensions (28), and otherwise. The Euler constant and , both, like the scale , artifacts of dimensional regularization, are subtracted along with the singularity. This defines the scheme in the polarized case, see, e.g., Ref. (29).

Finally, the last term in Eq. (6), , contains all the remaining, finite contributions, and the phase-space integration can be performed numerically in four dimensions. Monte Carlo integrations of (3) for different measurement functions can be done in parallel by randomly generating a sufficiently large sample of final-state configurations characterized by and to account for the possibly large cancellations among the various terms. Remnants of the regularization of soft and collinear regions of phase are the mathematical distributions appearing in the different contributions to Eq. (6). Their proper definitions through a test function, properties, and numerical treatment are discussed at length in Refs. (25); (26); (19) and need not be repeated here.

## 3 Phenomenological Studies

### 3.1 Preliminaries

Based on our Monte Carlo code described above, we present a comprehensive study of charm quark photoproduction in longitudinally polarized muon-deuterium collisions at a c.m.s. energy of relevant for the COMPASS experiment at CERN. Unless specified otherwise, we will show results for single-inclusive cross sections, differential in the transverse momentum of the observed meson.

The hadronization of the heavy quarks produced in the hard scattering is modeled using Eq. (4) with (30)

(17) |

The fragmentation function is normalized such that its first moment is unity, i.e., . For the remaining free parameter in (17) we take as the default value in all our calculations; see Table 4 in Ref. (31). To estimate the uncertainties associated with the choice of , we will vary it in the range (31). Like in the experimental analyses (15); (22); (23); (24), we adopt a cut throughout our phenomenological studies, where is the fraction of the energy of the photon taken by the meson in the laboratory frame.

In the computation of the LO and NLO unpolarized cross sections we use the LO and NLO CTEQ6 parton densities (32) and values for the strong coupling , respectively. In the polarized case, we adopt the best fit from the comprehensive global analysis performed by the DSSV group (2), which is the only set of helicity PDFs including constraints on from data. The resulting PDFs are characterized by a small in the -range predominantly probed by RHIC experiments, , with a node at (2). To study the sensitivity of charm photoproduction data to different , we also use the GRSV “standard” set (33), which has a positive that is larger in size than the one of DSSV. We note that the most recent RHIC data (7); (9) tend to prefer a polarized gluon density somewhere in between the ones obtained in the DSSV and GRSV fits (8); (9).

Equation (3) consistently includes both direct and resolved photon contributions to the spin-dependent photoproduction cross section by defining as the convolution of the polarized lepton-to-photon splitting function and the structure functions of a circularly polarized photon, i.e.,

(18) |

Nothing is known experimentally about the , which account for the hadronic structure of a photon in the resolved case where in (3). However, as will be demonstrated below, uncertainties associated with turn out to be negligible for charmed meson photoproduction at COMPASS kinematics. Contrary to inclusive hadron production (18), the resolved contribution turns out to be numerically small even when estimated with the “maximum” model for introduced in Ref. (34) and based on saturating the positivity limit at some low scale using the unpolarized set of GRV photon distributions (35) as reference. The direct part of the cross section (3) where is obtained by setting

(19) |

in Eq. (18).

The collinear emission of a quasi-real photon with low virtuality and momentum fraction off a muon with mass is given by the Weizsäcker-Williams equivalent photon spectrum which reads in the polarized case (36)

(20) | |||||

The upper limit is determined by experimental conditions. For COMPASS we take and, in addition, restrict the fraction of the muon’s momentum taken by the quasi-real photon to the range .

We use as the value of the charm quark mass for all our results. For the factorization and renormalization scales in Eq. (3) we take with as the central value and where denotes the transverse momentum of the charm quark. As is commonly done, we vary them simultaneously in the range to estimate the residual scale dependence at NLO which represents the dominant source of theoretical uncertainty and can be taken as a rough measure of yet unknown higher order corrections.

### 3.2 Numerical Results

We begin our detailed numerical studies with a discussion of the scale dependence of the unpolarized and polarized photoproduction cross sections for mesons in muon-deuterium collisions at COMPASS. Figure 2 shows the scale ambiguity from varying simultaneously in the range for the dominant direct photon contribution () only. The solid lines denote our default choice of scales, . Since the published COMPASS data (15) are divided into three bins in the energy of the produced meson, , , and , we impose the cut in Fig. 2. Results for the other two bins in are very similar and hence not shown.

As can be seen, the theoretical uncertainties due to the choice of scale are quite sizable, both in the unpolarized and in the polarized case and can be even further inflated by varying and independently as was done, for instance, in our study of heavy quark hadroproduction (19). Also, possible variations of , which we do not pursue here, would add to the theoretical error. By comparing the solid and dashed lines in Fig. 2 one can infer the relevance of the resolved photon cross section (). As was argued in Sec. 3.1, its contribution is rather small compared to the direct photon cross section and potential uncertainties due to the unknown do not matter.

Next, we study the role of the individual partonic subprocesses and their contribution to the photoproduction cross section in Eq. (3). Figure 3 displays the decomposition of the unpolarized direct photon cross section for , as shown in Fig. 2, in terms of the and processes. Again, we choose , and results for the other two bins in are very similar. As expected, the photon-gluon fusion mechanism, which is already present at LO, gives the main contribution to the cross section. The genuine NLO photon-quark channel yields a negative but small correction. The NLO result for the meson production cross section is roughly a factor of two larger than the corresponding estimate at LO as can be gathered from comparing the dashed and solid lines in Fig. 3.

A similar exercise in the polarized case is shown in Fig. 4 for two different choices of helicity PDFs. First, one notices that the results obtained with the DSSV and GRSV sets differ in sign and magnitude of the cross section, which is readily explained by the very different gluon densities in both sets. The positive definite of GRSV leads to a similar decomposition into and subprocesses as was observed in the unpolarized case in Fig. 3. Again, the cross section is strongly dominated by photon-gluon fusion, and the channel only yields a small correction. On the contrary, the oscillating of the DSSV set of helicity PDFs leads to a negative . Since the DSSV gluon is much smaller in size than the one of GRSV, the genuine NLO photon-quark contribution, which is numerically very similar in both PDF sets, is more important and yields more than a quarter of the cross section at small . Another important observation concerns the relevance of NLO corrections which appears to be very different for the DSSV and GRSV helicity PDFs. This implies that higher order QCD effects do not cancel in the experimentally relevant double-spin asymmetry,

(21) |

which we shall discuss in more detail below. Clearly, for a reliable quantitative analysis of charm photoproduction in terms of polarized PDFs, preferably as part of a global QCD fit, it is indispensable to properly include NLO corrections.

Figure 6 illustrates the uncertainties due to our choice of in the fragmentation of the charm quark into the observed meson. As in Fig. 2, we show the direct photon contribution for , but here we do not impose any cuts on the meson to allow for a better comparability with the cross sections on the charm quark level which are presented as well. Compared to the factorization scale ambiguities, the dependence of our results on is fairly moderate for both the unpolarized and the polarized cross sections. A similar observation was made in the case of hadroproduction (19). Since charm quarks lose only very little of their momentum during hadronization, i.e., is peaked at fairly large values of , the convolution (4) introduces only a rather small shift in the transverse momentum spectrum of the charm quarks. This can be inferred from the curves labeled in Fig. 6.

Another interesting question concerns the range of momentum fractions predominantly probed in the PDFs by the COMPASS data. Due to the dominance of the photon-gluon fusion process, charm photoproduction will mainly lead to a constraint on the gluon helicity distribution, which is the prime motivation for such measurements. The distribution in LO and NLO, for both the unpolarized and polarized direct photon contribution to the cross section (3) is explored in Fig. 6. No cuts are imposed in this calculation and the meson spectrum is integrated over the entire phase space. As has to be expected from the fact that the PDFs enter the cross section (3) through a convolution, a broad range of values is sampled. It turns out, however, that the mean value of , where the distribution is peaked, is fairly independent not only of the order in perturbation theory, LO or NLO, but to a large extent also of the chosen set of polarized or unpolarized PDFs. We roughly estimate the average momentum fraction to be with an error of about .

Our results differ from preliminary estimates of by the COMPASS collaboration (22); (23); (24), where NLO results have been obtained based on some parton shower Monte Carlo to approximate the phase space for the NLO matrix elements of Ref. (16). Significant differences between estimated in LO and NLO are found in this way. While their LO result for agrees with our estimate of about 0.08, their preliminary NLO result is . Since the details of the method are not yet published, it is not yet clear how these results can be compared to our full NLO calculation. We also note that once data on photoproduction processes are implemented in global QCD analyses of helicity PDFs, information on , though useful, is no longer required or relevant as the fits automatically impose the constraints from data for any given functional form assumed for the .

Next, we turn to the experimentally measured double-spin asymmetry, defined in Eq. (21), which was analyzed at LO accuracy and under certain simplifying assumptions in terms of the mean gluon polarization by the COMPASS collaboration (15); see the discussion below. Preliminary NLO estimates are also available at NLO with the hybrid method outlined above (22); (23); (24). In Fig. 7 we show for the two sets of helicity PDFs used throughout this paper, by computing the ratio of the cross sections shown in Figs. 2 and 3 for . To resolve the differences in obtained with the DSSV and GRSV PDFs, which mainly stem from , an experimental precision of at least needs to be achieved.

We also compare NLO and LO estimates of for both sets of helicity PDFs in the lower panel of Fig. 7 and find rather different patterns depending on . At small , the NLO is about a factor of two larger than the LO estimate for the DSSV set whereas a reduction by roughly the same amount is found with the GRSV PDFs. This illustrates that any approximations for the spin asymmetry, either to assume a cancellation of NLO corrections or a constant pattern independent of the choice of PDFs, are not justified and should not be used for analyzing data. Again, only a global analysis will lead to consistent results. Qualitatively very similar results have been obtained for other cuts on the energy of the observed meson, and .

Finally, we compare our calculations at NLO accuracy with the available results from the COMPASS experiment. Figure 8 shows the data in three bins of the energy of the detected meson as a function of its transverse momentum , imposing the cut . Note that instead of using the published data (15), we show new, preliminary results presented recently in Ref. (22); (23); (24). A weighted average is performed to combine the results for the three decay channels , , and listed in (22); (23); (24). We adopt the photon polarization dilution factors tabulated in Ref. (22); (23); (24) to convert the data for given in the photon-nucleon system to the double-spin asymmetries for muon-nucleon scattering computed with our Monte Carlo code.

Clearly, present experimental uncertainties are too large to discriminate between different spin-dependent gluon densities. Apart from the DSSV and GRSV “standard” sets, we also use an extreme GRSV set, characterized by a very large and negative gluon density based on setting at some low bound-state like input scale (33) for the evolution. Despite leading to a distinctively different than our two default sets, all results are compatible with data within the experimental precision.

To overcome the statistical limitations, fewer bins have to be used or all data need to be combined. COMPASS has performed such an analysis (15) yielding the result for the gluon polarization quoted in Fig. 1, which is compatible with similar extractions of from hadron production data. These kind of analyses are LO estimates, assuming, in addition, that the convolutions of the PDFs with the partonic hard cross sections in Eq. (3) can be approximated as

(22) |

Here, the “analyzing power” is evaluated at some average kinematics. As mentioned before, in a proper global analysis of helicity PDFs one does not need to resort to any of these approximations whose validity is difficult to estimate or justify. With our new Monte Carlo program it is now possible to perform such an analysis at NLO accuracy in the future. As we have demonstrated in some detail, NLO corrections are indispensable for a quantitative analysis. Most importantly, they do not cancel in the spin asymmetry as one might naively expect.

Apart from the phenomenological applications presented here, our code can be also used to estimate spin asymmetries for charm and bottom photoproduction and their impact on our knowledge of the spin structure of the nucleon at higher c.m.s. energies intended for a first polarized lepton-ion collider. Details on the EIC project can be found in (10). Such a machine will be indispensable to finally unravel the quark and gluon contributions to the nucleon spin (11).

## 4 Summary and Outlook

We have presented a flexible parton-level Monte Carlo program to compute heavy flavor distributions at NLO accuracy in longitudinally polarized lepton-nucleon collisions in the photoproduction regime. For the first time, we consistently include both direct and resolved photon contributions and found the latter to be negligibly small for charm production at COMPASS kinematics. Experimental acceptance cuts, the hadronization of the produced heavy quark pair, and, if needed, their subsequent semi-leptonic decays can be included in phenomenological applications.

Heavy flavor photoproduction receives its importance for the field of spin physics from its expected strong sensitivity to the polarized gluon density which we confirm. In general, higher order corrections are found to be sizable and strongly dependent on the chosen set of helicity PDFs. This is also true for the experimentally relevant spin asymmetry despite naive expectations that QCD corrections cancel in the ratio. Theoretical uncertainties due to the choice of the factorization scale are sizable even at next-to-leading order accuracy while ambiguities from the exact form of the charm quark fragmentation function are less important.

The results obtained in this paper allow one to include data on charm photoproduction consistently into future global analyses of helicity PDFs. We have shown that data from COMPASS can lead in principle to a constraint on the polarized gluon density at a momentum fraction of about 0.1. Currently available data will have, however, very little impact on existing fits due to the size of the experimental uncertainties which are too large to discriminate between different gluon densities.

In addition to the phenomenological studies performed in this paper, our code will be useful in assessing the physics impact of heavy quark photoproduction at a possible future polarized lepton-ion collider like the EIC project.

## Note Added

## Acknowledgements

J.R. was supported by a grant of the “Cusanuswerk”, Bonn, Germany. M.S. acknowledges support by the U.S. Department of Energy (DOE) under contract No. DE-AC02-98CH10886 and, in part, by a BNL “Laboratory Directed Research and Development Program” (LDRD 12-034). This work was supported in part by the “Bundesministerium für Bildung und Forschung”, Germany.

### References

- For a recent overview and discussions of the spin sum rule and orbital angular momentum, see the talks given at the INT workshop on “Orbital Angular Momentum in QCD”, INT, Seattle, 2012, http://www.int.washington.edu/PROGRAMS/12-49w.
- D. de Florian, R. Sassot, M. Stratmann, and W. Vogelsang, Phys. Rev. Lett. 101, 072001 (2008); Phys. Rev. D 80, 034030 (2009).
- J. Blumlein and H. Bottcher, Nucl. Phys. B 841, 205 (2010); E. Leader, A. V. Sidorov, and D. B. Stamenov, Phys. Rev. D 82, 114018 (2010).
- B. Jager, A. Schafer, M. Stratmann, and W. Vogelsang, Phys. Rev. D 67, 054005 (2003).
- B. Jager, M. Stratmann, and W. Vogelsang, Phys. Rev. D 70, 034010 (2004).
- A. Adare et al. [PHENIX Collaboration], Phys. Rev. D 76, 051106 (2007); Phys. Rev. Lett. 103, 012003 (2009); Phys. Rev. D 79, 012003 (2009); B. I. Abelev et al. [STAR Collaboration], Phys. Rev. Lett. 97, 252001 (2006); Phys. Rev. Lett. 100, 232003 (2008); L. Adamczyk et al. [STAR Collaboration], Phys. Rev. D 86, 032006 (2012).
- P. Djawotho [for the STAR Collaboration], arXiv:1106.5769; J. Phys. Conf. Ser. 295 (2011) 012061.
- D. de Florian, R. Sassot, M. Stratmann, and W. Vogelsang, talk presented at “DIS2011”, April 2011, Newport News, VA, arXiv:1108.3955; Prog. Part. Nucl. Phys. 67, 251 (2012).
- E.C. Aschenauer et al., report on “The RHIC Spin Program: Achievements and Future Opportunities”, BNL, 2012, http://www.bnl.gov/npp/docs/ RHIC-Spin-WriteUp-121105.pdf
- D. Boer et al., “INT report on EIC Science”, arXiv:1108.1713.
- E. C. Aschenauer, R. Sassot, and M. Stratmann, Phys. Rev. D 86 (2012) 054020.
- A. Airapetian et al. [HERMES Collaboration], Phys. Rev. Lett. 84, 2584 (2000); JHEP 1008, 130 (2010).
- B. Adeva et al. [Spin Muon Collaboration (SMC)], Phys. Rev. D 70, 012002 (2004).
- E. S. Ageev et al. [COMPASS Collaboration], Phys. Lett. B 633, 25 (2006); C. Adolph et al. [COMPASS Collaboration], arXiv:1202.4064.
- M. Alekseev et al. [COMPASS Collaboration], arXiv:0802.3023; Phys. Lett. B 676, 31 (2009).
- I. Bojak and M. Stratmann, Phys. Lett. B 433, 411 (1998); Nucl. Phys. B 540, 345 (1999) [Erratum-ibid. B 569, 694(E) (2000)].
- Z. Merebashvili, A. P. Contogouris, and G. Grispos, Phys. Rev. D 62, 114509 (2000) [Erratum-ibid. D 69, 019901 (2004)]; Phys. Lett. B 482, 93 (2000).
- B. Jager, M. Stratmann, and W. Vogelsang, Phys. Rev. D 68, 114018 (2003); Eur. Phys. J. C 44, 533 (2005); C. Hendlmeier, A. Schafer, and M. Stratmann, Eur. Phys. J. C 55, 597 (2008).
- J. Riedl, A. Schafer, and M. Stratmann, Phys. Rev. D 80, 114020 (2009).
- I. Bojak and M. Stratmann, Phys. Rev. D 67, 034010 (2003).
- M. Cacciari, P. Nason and R. Vogt, Phys. Rev. Lett. 95, 122001 (2005).
- K. Kurek, Habilitation Thesis, National Centre of Nuclear Research, Swierk, November 2011.
- C. Franco, Ph.D. Thesis, Universidade Tecnica de Lisboa, December 2011.
- C. Franco [for the COMPASS Collaboration], arXiv:1208.6567.
- M. L. Mangano, P. Nason, and G. Ridolfi, Nucl. Phys. B 373, 295 (1992).
- S. Frixione, M. L. Mangano, P. Nason, and G. Ridolfi, Nucl. Phys. B 412, 225 (1994).
- R. K. Ellis and P. Nason, Nucl. Phys. B 312, 551 (1989); J. Smith and W. L. van Neerven, Nucl. Phys. B 374, 36 (1992); W. Beenakker, H. Kuijf, W. L. van Neerven, and J. Smith, Phys. Rev. D 40, 54 (1989); W. Beenakker, W. L. van Neerven, R. Meng, G. A. Schuler, and J. Smith, Nucl. Phys. B 351, 507 (1991); P. Nason, S. Dawson, and R. K. Ellis, Nucl. Phys. B 327, 49 (1989) [Erratum-ibid. B335, 260 (1990)].
- G. ’t Hooft and M. J. G. Veltman, Nucl. Phys. B 44, 189 (1972); P. Breitenlohner and D. Maison, Commun. Math. Phys. 52, 11 (1977).
- R. Mertig and W. L. van Neerven, Z. Phys. C 70, 637 (1996); W. Vogelsang, Phys. Rev. D 54, 2023 (1996); Nucl. Phys. B 475, 47 (1996).
- V. G. Kartvelishvili, A. K. Likhoded, and V. A. Petrov, Phys. Lett. B 78, 615 (1978).
- For a review, see, e.g., J. Baines et al., summary report of the “Heavy Quarks Working Group” for the “HERA-LHC Workshop” proceedings, arXiv:hep-ph/0601164 and references therein.
- J. Pumplin, D. R. Stump, J. Huston, H. L. Lai, P. M. Nadolsky, and W. K. Tung, JHEP 0207, 012 (2002).
- M. Glück, E. Reya, M. Stratmann, and W. Vogelsang, Phys. Rev. D 63, 094005 (2001).
- M. Glück and W. Vogelsang, Z. Phys. C 55, 353 (1992); C 57, 309 (1993); M. Glück, M. Stratmann, and W. Vogelsang, Phys. Lett. B 337, 373 (1994); M. Stratmann and W. Vogelsang, Phys. Lett. B 386, 370 (1996).
- M. Glück, E. Reya, and A. Vogt, Phys. Rev. D 46, 1973 (1992).
- D. de Florian and S. Frixione, Phys. Lett. B 457, 236 (1999).
- C. Adolph et al. [COMPASS Collaboration], arXiv:1211.6849.