# The Gluon Sivers Distribution: Status and Future Prospects

Authors: Daniël Boer; Cédric Lorcé; Cristian Pisano; Jian Zhou
Journal: Advances in High Energy Physics
 (2015)

---

## Abstract

We review what is currently known about the gluon Sivers distribution and what are the opportunities to learn more about it. Because single transverse spin asymmetries inp↑p→πX provide only indirect information about the gluon Sivers function through the relation with the quark-gluon and tri-gluon Qiu-Sterman functions, current data from hadronic collisions at RHIC have not yet been translated into a solid constraint on the gluon Sivers function. SIDIS data, including the COMPASS deuteron data, allow for a gluon Sivers contribution of natural size expected from large Nc arguments, which is O(1/Nc) times the nonsinglet quark Sivers contribution. Several very promising processes to measure the gluon Sivers effect directly have been suggested, which besides RHIC investigations, would strongly favor experiments at AFTER@LHC and a possible future Electron-Ion Collider. Due to the inherent process dependence of TMDs, the gluon Sivers TMD probed in the various processes are different linear combinations of two universal gluon Sivers functions that have different behavior under charge conjugation and that therefore satisfy different theoretical constraints. For this reason both hadronic and DIS type of collisions are essential in the study of the role of gluons in transversely polarized protons.


---

## Body

## 1. The Sivers Function and Its Definition

The distribution of quarks and gluons in a proton (or any other spin-1/2 hadron) that is polarized transversely to its momentum need not be left-right symmetric with respect to the plane spanned by the momentum and spin directions. This asymmetry is called the Sivers effect [1]. It results in angular asymmetries of produced particles in high energy scattering processes involving a transversely polarized hadron. Experimental data in support for such a left-right asymmetry in the quark distribution was first obtained from semi-inclusive DIS process by the HERMES collaboration [2]. This review is about what is currently known about the gluonic Sivers effect distribution. The Sivers effect is of great interest theoretically as it is very sensitive to the color flow in the scattering process and to the multitude of color exchanges among initial and final states. It is the first quantity for which this has been recognized and for which color flow sensitivity can be tested unambiguously [3]. Verification of its unusual properties will provide a strong test of the formalism of Transverse-Momentum-Dependent parton distributions. It is an important quantity of nonperturbative QCD to consider both qualitatively and quantitatively. This review discusses these aspects for the gluon Sivers distribution specifically. We first start with its proper definition.The number density in momentum space of a generic parton (quark, antiquark, and gluon) inside a hadron with massM, transverse polarization ST, and momentum P can be written as(1)f^x,k⊥;ST=f1x,k⊥2-P^×k⊥·S^TMf1T⊥x,k⊥2,where f1(x,k⊥2) is the unpolarized Transverse-Momentum-Dependent (TMD) parton distribution, S^T≡ST/|ST|, and P^≡P/|P|. The function f1T⊥(x,k⊥2) describing the distortion in the distribution of unpolarized partons with light-front momentum fraction x and transverse momentum k⊥ due to the transverse polarization of the hadron is called the Sivers function. The notation used here comes from [4], but also the notation ΔNfg/h↑ from [5, 6] is sometimes used, where ΔNfg/h↑=-2(k⊥/M)f1T⊥g (analogous to the quark case [7]).The Sivers function satisfies the following positivity bound [8]:(2)k⊥Mf1T⊥x,k⊥2≤f1x,k⊥2.In [8] the operator definition of the gluon Sivers function, which was called GT=-f1T⊥g [9], was first given without gauge links. The definition including gauge links then appeared in [10, 11] (3)-P^×k⊥·S^TMf1T⊥gx,k⊥2=12f^x,k⊥;ST-f^x,k⊥;-STwith(4)f^x,k⊥;ST=δTjlxP+∫dz-d2z⊥2π2eik·zP,ST2·Tr⁡F+j0U0,zF+lzUz,0P,STz+=0,where U[a,b] is a Wilson line connecting the points a and b along a contour determined by the physical process and k+=xP+ is the fraction of parton light-front momentum. For a proper definition that is free from rapidity divergences associated with gauge links with paths (partly) along the light front, a redefinition involving the so-called soft factor is necessary [12–15]. This will however not play a significant role here and can simply be considered as implicit.


## 2. Sivers Effect,AN, and Qiu-Sterman Effect

The Sivers effect (for both quarks and gluons) was first suggested in [1] as an explanation for the large left-right single transverse spin asymmetries (AN) observed in p↑p→πX [16–23] (and similar asymmetries in K [22], in η [24], and tentatively in J/ψ [25] production). The Sivers effect was first studied phenomenologically in [5]. However, extraction of the Sivers TMD presumes all-order TMD factorization. A factorized description of the process pp→πX only applies for large transverse momentum pT of the produced pions (say for pT≳1 GeV), where in fact collinear factorization is appropriate [26], rather than TMD factorization. Although TMDs do appear in the phenomenological description of AN in [5] and subsequent studies (see [27]), that description is thus not based on a TMD factorization theorem. Rather it should be considered as an effective model description, now commonly referred to as the Generalized Parton Model (GPM). As a phenomenological approach it has proven useful in the quest to disentangle the possible underlying mechanisms of the spin asymmetries; see [27, 28] for more discussion, but the extracted “effective” TMDs may differ from the TMDs extracted from TMD-factorizing processes. What is known about the effective gluon Sivers TMD will be discussed below.In collinear factorization the single spin asymmetry (SSA) will arise at the twist-3 level [29]. In this description AN probes the (quark-gluon) Qiu-Sterman function (sometimes one factor of the coupling constant g is included in the definition of Qiu-Sterman functions, because one always encounters them multiplied by g) [30, 31](5)Tq,Fx,x=MP+∫dz-dη2πeik·zP^×S^Tj2MP,Sψ¯0·γ+F+jηzψzP,Sz+=z⊥=0,and its trigluon correlation analogues TG(f)(x,x)=TG(+)(x,x) and TG(d)(x,x)=TG(-)(x,x) [32–35](6)TG±x,x=-2MδTlmxP+2∫dz-dη2πeik·zP^×S^Tj2MP,S·C±abcFa+l0Fb+jηzFc+mzP,Sz+=z⊥=0,where C+abc=ifabc and C-abc=dabc and where the light-front gauge A+=0 has been considered for convenience.In [36] the quark-gluon Qiu-Sterman function has been related to the first transverse moment of the quark Sivers function, that is, f1T⊥(1)q(x)∝Tq,F(x,x)/M, where(7)f1T⊥1qx≡∫d2k⊥k⊥22M2f1T⊥qx,k⊥2.However, that relation was only established at tree level (beyond tree level the relation will be affected by the considered regularized definition of the Sivers TMD including its dependence on the soft factor). A similar tree level relation can be established in the gluon sector as well: f1T⊥(1)g(x)∝TG(x,x)/M (where TG that appears in this relation depends on the gauge links; see the discussion in Section 5).Another relation has been established in [37]:(8)f1T⊥qx,k⊥2~k⊥2≫M2αsMk⊥4K⊗Tq,Fx,which means that the quark-gluon Qiu-Sterman function determines the large transverse momentum tail of the quark Sivers function. Here it should be emphasized that the function K⊗Tq,F(x) consists not only of a convolution of Tq,F(x,x) but also of its derivative x∂Tq,F(x,x)/∂x and of the more general Tq,F(x,y) with y≠x. It corresponds to the fact that the evolution of Tq,F is nonautonomous and inhomogeneous; see [38]. Note that here we have discussed the nonsinglet contributions only that apply to combinations like u minus d quarks; otherwise also gluonic contributions need to be taken into account.Similarly, the tail of the gluon Sivers function is determined by several Qiu-Sterman functions [39]. It receives contributions from the quark-gluon Qiu-Sterman functions Tq,F(x,x) and Tq,F(x,y) with y≠x and from the trigluon functions TG(f/d)(x,x) (whose contribution to the gluon Sivers function depends on the gauge links; see Section 5). At small x the situation simplifies: the contributions from Tq,F(x,x) and Tq,F(x,y) with y≠x to the tail of the gluon Sivers function cancel each other [39], leaving only the trigluon correlators. Moreover, TG(d)(x,x) evolves with the same 1/x behavior at small x as the unpolarized gluon distribution and is therefore not necessarily suppressed at high energies and small values of x, whereas TG(f)(x,x) lacks this 1/x enhancement [35]. As the large pT  AN data from RHIC are generally not in the small-x region of the polarized proton, except for negative xF, any simplifications at small x should of course first be tested for validity.Information fromAN measurements at sufficiently large pT (in order to consider a collinear factorization description in the first place) can thusin principle reflect some information on Sivers functions (i.e., on the tails and perhaps also on first transverse moments), butin practice other twist-3 contributions beside the mentioned Qiu-Sterman functions, namely, chiral-odd and fragmentation function analogues, contribute to AN [28, 40–42]. From the smallness of AN in the midrapidity and backward (negative xF) regions, one would generally conclude that gluonic and sea quark contributions to the transverse single spin asymmetry are not large, but a detailed analysis is required to determine precisely the size of the various contributions.A
N for π0 production at midrapidity has been measured by the PHENIX experiment in polarized pp collisions at RHIC and was found to be consistent with zero, for pT values below 5 GeV at the permille level and for higher pT values (up to 11 GeV) at the few percent level [21, 43]. These data taken at s=200 GeV probe x values only down to x~0.006, where still a combination of Qiu-Sterman functions is expected to contribute. In [44] these π0 data were discussed, using two models for the trigluon Qiu-Sterman functions that were constrained from AN in D-meson production [45, 46]. The midrapidity π0 data are shown to mostly constrain TG(f)(x,x) (their N(x)) at low pT. The authors conclude that “both models give tiny asymmetry due to the small partonic cross sections, so the form of the three-gluon correlation functions is not much constrained by the data in this region.” From this limited model study of both π0 and D production one would conclude that TG(f)(x,x) and TG(d)(x,x) are in any case small, a permille fraction of x times the unpolarized gluon distribution. This should be investigated further with more general model forms that adhere to the correct small-x behavior and with more precise data. The experimental precision of AN can be improved much further both at RHIC and especially at the AFTER@LHC experiment [47], which would have a luminosity factor of 10–100 times higher, if not more. Such improvement and the measurements of asymmetries for many different types of produced particles are required to separately constrain or determine the various Qiu-Sterman contributions.In the GPM the smallness ofAN at midrapidity puts strong constraints on the effective gluon Sivers function. As explained, this gluon Sivers function captures the combined effect from several Qiu-Sterman contributions and may thus differ from the actual gluon Sivers function obtained from TMD-factorizing processes. In a recent GPM analysis [48], which is an updated analysis of [49], thebest fits to the PHENIX midrapidity π0AN data indeed correspond to a small effective gluon Sivers function with respect to its theoretical bound determined by the unpolarized gluon. For example, for x<0.1 it is at most only a few percent of the bound. However, the maximally allowed effective gluon Sivers function is still sizable though. Its first transverse moment is still found to be around 30% of the up quark Sivers function in the region 0.06<x<0.3, which is consistent with findings from semi-inclusive DIS for the actual gluon Sivers function and also with theoretical expectations, as discussed in the next section. In addition, it should be mentioned that this GPM analysis assumes a Gaussian k⊥ dependence, which corresponds neither to the correct power-law tail of the Sivers function, (8), nor to the unpolarized gluon distribution. Given all the caveats that come with these results, one should be careful to draw a definite conclusion about the size of the actual gluon Sivers TMD from AN data.One should also specify clearly what one calls a small gluon Sivers function. It will depend on what one compares to, that is, whether that is to the unpolarized gluon that grows very rapidly at smallx or to the up or down quark Sivers function for not too small x. At small x it becomes very important whether one discusses the f or d type contribution (see Section 5), which is an issue not addressed in the GPM studies of AN.


## 3. Sivers Asymmetry in SIDIS

The Sivers effect leads to asin⁡(ϕh-ϕS) asymmetry in semi-inclusive DIS (SIDIS) [4], which has been observed in experiments using a proton target by HERMES [2, 50] and COMPASS [51] and using 3He target by Jefferson Lab Hall A [52, 53]. The data follow to quite a good extent the expectations of a valence quark picture in the target and of favored fragmentation. In the proton case, π+ thus shows the largest asymmetry (for large z values, the asymmetry is around 4-5% or even somewhat larger when a lower cut of Q2>4 GeV2 is implemented instead of Q2>1 GeV2 [50]). The π- asymmetries are smaller and still compatible with zero. The K± asymmetries are similar to the π± asymmetries, but with larger errors. Sivers asymmetries on a deuteron target [54] are all consistent with zero. Fits to all these HERMES and COMPASS data, including the deuteron data using isospin symmetry, indicate that the Sivers function (f1T⊥) for u quarks in a proton is negative and for a d quark in a proton is positive and approximately equal in absolute value [55]. This fits the expectations from the limit of a large number of colors Nc [56, 57]:(9)f1T⊥ux,k⊥2=-f1T⊥dx,k⊥2+O1Nc.The flavor singlet combination of u and d is of the same order as the gluon contribution in Nc counting [58]. The latter is thus 1/Nc suppressed with respect to the flavor nonsinglet quark Sivers effect at not too small x (x~1/Nc) [59].Within the current accuracy, the SIDIS data do not require any sea quark or gluon contributions, which among other considerations (see Section6) led Brodsky and Gardner to conclude that the gluon Sivers function is small or even zero (“absence of gluon orbital angular momentum”) [60]. The SIDIS data from HERMES, COMPASS, and Jefferson Lab Hall A are of course at rather modest Q2 and not too small-x values, that is, in the valence region. One cannot yet draw any conclusions about the gluon Sivers function at higher Q2 and smaller values of x. Moreover, the data certainly still allow for gluon Sivers contributions of the order of 1/Nc times the valence quark Sivers functions. This is evident from the fits by Anselmino et al. [55], where the first transverse moment of the u and d Sivers functions has error bands that are at least around 30% of the central values.Note that the SSA in the “inclusive” processep→hX, where the back-scattered lepton is not observed [61, 62], does not allow for an interpretation in terms of TMDs, as the data are dominated by Q2≈0. Even for large pT the appropriate factorization would be collinear factorization and the Sivers type of asymmetry would probe the Qiu-Sterman functions instead [63], which as discussed above have some relation to the Sivers TMDs, but only via the tail or possibly via the first transverse moment. The asymmetries for pT>1 GeV are found to be at the level of 5–10% for positive hadrons. Fits will need to make clear how much room there is for a gluon Qiu-Sterman effect. Given the fact that the gluon Qiu-Sterman function does not enter at leading order in αs in this process, this room may be considerable.


## 4. Sivers Asymmetry in Other Processes

Several otherpp scattering processes to access the gluon Sivers function have been suggested over the past years: p↑p→jetjetX [64], p↑p→DX [34, 44, 65], p↑p→γX [66], p↑p→γjetX [66, 67], p↑p→γ∗X→μ+μ-X [66], p↑p→jetX (single transverse spin asymmetries in jet production measured at RHIC [68, 69] at forward rapidities (the valence region) show very small asymmetries, which is probably due to a cancellation among u and d quark contributions [70]), p↑p→πjetX [71], and p↑p→ηc/bX [39]. Several of these processes are like AN in (high-pT) pion production, which means that they deal with twist-3 collinear factorization and only provide indirect or limited information about the gluon Sivers TMD. Several other processes run into the problem of TMD factorization breaking contributions [72] and hence are not safe. In principle they do probe TMDs but as a result of TMD factorization breaking contributions, conclusions about the gluon Sivers function from their measurements cannot be drawn safely. This applies, for instance, to the process p↑p→jetjetX (measured at RHIC to be small at the few percent level [73]), which moreover suffers from cancellations between u and d contributions and between the effects of initial and final state interactions [74–76]. TMD factorization breaking would also apply to open heavy quark production: p↑p→QQ¯X, such as p↑p→D0D¯0X; compare, for example, [77]. Whether the problem also applies to double heavy quarkonium production remains to be seen, because in practice the color singlet contributions may give the dominant contribution in that case. Among the hadronic collisions the processes having one or two color singlets in the final state would in any case be the safest. One very promising example is p↑p→γjetX [67], where it depends on the rapidity of the photon and the jet, that is, on the x fraction of the parton in polarized proton, whether the gluon Sivers function dominates over the quark one or vice versa. Another very promising example is p↑p→J/ψγX, which is predominantly initiated by gluon-gluon scattering (which is an order in αs higher than the gluon contribution in p↑p→γjetX) and for which the color singlet contribution dominates over the color octet one to a large extent [78, 79]. The same applies to p↑p→J/ψJ/ψX (see the contribution by Lansberg and Shao in this special issue). AFTER@LHC would be very well-suited for studying these processes.SSA experiments could be done at AFTER@LHC where the beam of protons or lead ions of the LHC would collide with a fixed target that is transversely polarized. Suchpp↑ and Pbp↑ collisions would have a center-of-mass energy sNN of 115 and 72 GeV, respectively, and have high luminosity and good coverage in the rapidity region of the transversely polarized target (mid and large xp↑) [47]. Polarized Drell-Yan and prompt photon production studies could be done to measure the quark Sivers function very precisely, perhaps to the level that the gluon Sivers function becomes relevant, despite the large values of x in the polarized target. As mentioned γjet and J/ψγ production could be used to study the gluon Sivers effect directly, where the former would need specific selection of the rapidities. In addition, the comparison of Pbp↑→γjetX and pp↑→γjetX would give a further handle on determining the relative sizes of quark and gluon Sivers functions. Other processes, such as D-meson or J/ψ production, would allow a similar study of Qiu-Sterman functions, including the trigluon ones, which are of course interesting in their own right. See [28] for a more detailed and quantitative study of twist-3 transverse single spin asymmetries in proton-proton collisions at the AFTER@LHC experiment. All these possibilities offer a very interesting complementary opportunity or even a competitive alternative to the other existing high-energy particle physics spin projects aiming at studying the role of gluons in transversely polarized protons.In electron-proton scattering one of the most promising processes to directly probe the gluon Sivers function is open charm production,ep↑→e′cc¯X, which could ideally be studied with an Electron-Ion Collider (EIC). By selecting the charm (or bottom) quark, one effectively eliminates the subprocesses γ∗q→qg and γ∗q¯→q¯g and becomes essentially (this assumes that intrinsic charm contributions are suppressed by selecting sufficiently small-x values) sensitive to γ∗g→cc¯ and thus to the gluon Sivers function (a similar argument is used in the study of high-pT hadron pairs in muon-deuteron and muon-proton scattering [80, 81], where photon-gluon fusion is expected to dominate. The relevant asymmetry AUTsin⁡(ϕ2h-ϕS) is found to be -0.14±0.15(stat.)±0.06(syst.) at 〈xG〉=0.13 for the deuteron and -0.26±0.09(stat.)±0.08(syst.) at 〈xG〉=0.15 for the proton. For the interpretation of the data in terms of the gluon Sivers effect, Q2 and pT of each hadron need to be sufficiently large to trust factorization). Here the transverse momenta of the heavy quarks are considered to be almost back-to-back. There is no problem with TMD factorization breaking contributions of the type discussed in [72], but that does not mean that the process is as straightforward as SIDIS. Even in the case where one considers charm jets, one has to include a description of the transverse momentum distribution inside such a jet. It may be easier to consider D0D¯0 measurements (for a study of the twist-3 SSA in large pTD meson production in SIDIS, that is, ep↑→e′DX, see [33, 82]). In either case one deals with 3 TMDs. Such processes involve a different soft factor (in this case a vacuum correlator of 6 Wilson lines) compared to processes involving 2 TMDs as in SIDIS, affecting the predictability. This has been discussed at the one-loop level in [83]. The SSA in ep↑→e′D0D¯0X has been studied for some models of the gluon Sivers function in [84]; compare Section  2.3.1 of [84]. This may be the “smoking gun” process for the gluonic Sivers effect at an EIC. It should be mentioned though that it actually probes a different gluon Sivers TMD than the hadronic processes discussed above. This is discussed in the next section. It shows that hadronic processes are complementary to DIS processes.For completeness we mention that when comparing extractions of the gluon Sivers TMD from different processes, one has to take care not only of the process dependence but also of the different energy scales. Under TMD evolution from one scale to another, the transverse momentum distribution changes. For details we refer to [85–90].


## 5. Process Dependence of the Gluon Sivers Function

Once a set of processes that in principle allow probing the gluon Sivers TMD has been obtained, one still has to take into account the fact that such TMD is process-dependent. For quarks the famous overall sign change between the Sivers TMD probed in SIDIS and the one probed in Drell-Yan is expected [3, 91–93] and is currently under experimental investigation. For gluons the situation is more complicated as each gluon TMD depends on two gauge links (in the fundamental representation), so there are more possibilities [11, 94, 95]. The gauge link structure of the gluon distributions in ep→e′D0D¯0X differs from the one in, for instance, pp→γjetX (cf. [96] for the comparison at small x). Clearly, this will complicate the analysis of gluon Sivers effect which will involve more than one gluon Sivers function. In [11] it was demonstrated that any gluon Sivers function can be expressed in terms of two “universal” gluon Sivers functions:(10)f1T⊥gUx,k⊥2=∑c=12CG,cUf1T⊥gAcx,k⊥2,where the coefficients CG,c[U] are calculable for each partonic subprocess. The first transverse moments of the two distinct gluon Sivers functions are related (at least at tree level) to the two distinct trigluon Qiu-Sterman functions TG(f/d). Therefore, we will refer to the universal gluon Sivers functions as f1T⊥g(f) and f1T⊥g(d). The two functions have different behavior under charge conjugation, just like TG(f) is a matrix element of a C-even operator and TG(d) of a C-odd operator.The processep↑→e′D0D¯0X is dominated by just one partonic subprocess γg→qq¯ and thus probes the gluon Sivers function with two future-pointing (+) links [97], which is f1T⊥g(f) [11]. The process p↑p→γjetX probes the subprocesses qg→γq and qq¯→γg. If one selects kinematics such that one probes small-x values in the polarized proton, such that qg→γq dominates, then this process accesses the gluon Sivers with a future- and past-pointing link, which corresponds to f1T⊥g(d). The theoretical expectations are different for these two cases.


## 6. Theoretical Constraints on Sivers Function

Constraints on the unintegrated gluon Sivers TMDf1T⊥g(x,k⊥2) from fits have to take into account that it is theoretically possible that both the quark and the gluon Sivers TMD can have nodes in x and/or k⊥ [98, 99]. The possibility of a node in x is supported by the observation [35] that the splitting function for TG(f) is negative at small x, in analogy to the Δg case. Fits to SIDIS data (studied with a rather restrictive parameterization and in a restricted kinematic range) do not appear to require a node [100], but that does not exclude this possibility. Especially when comparing data from different kinematic regions and different processes, this option should be kept in mind. Nodes can of course have a large effect on integrals of Sivers functions, such as the first transverse moment (7) and its first Mellin moment (for parton a)(11)k⊥a=-MS^T×P^∫dxf1T⊥1ax,which is the average transverse momentum inside a transversely polarized target. The notation 〈k⊥a〉 comes from [101]. This quantity is related to the Sivers shift [102], the average transverse momentum shift orthogonal to the transverse spin direction, which is normalized to the zeroth transverse moment of the unpolarized TMD f1(0)(x)≡∫d2k⊥f1(x,k⊥2):(12)k⊥yxUTx=Mf1T⊥1x;μ,ζf10x;μ,ζ.Here only the y-component perpendicular to the transverse spin direction x is nonzero and therefore considered. Note that the Sivers shift depends in principle on the renormalization scale μ and a rapidity variable ζ but also depends on the path of the gauge link (see below). Furthermore, there is the theoretical issue whether the transverse moments converge. For that reason a regularized version using Bessel moments has been suggested in [102], which for quarks has been evaluated on the lattice in [103]. The lattice calculation of the Sivers shift for u-d is consistent with negative u and positive d Sivers functions in SIDIS, which are expectations that follow from general arguments on final state interactions [104], from a model-dependent relation to GPDs [105, 106] (see below) and from fits [107, 108].As said above,f1T⊥(1)q(x) has a (process-dependent!) relation to the Qiu-Sterman function Tq,F(x,x), which so far only has been established at tree level and for quarks. In addition, Burkardt has suggested a model-dependent relation between the integrated quantity f1T⊥(1)(x) and an integral over the GPD E(x,ξ,Δ2) at zero skewness ξ=0 (and only for quarks) [105]:(13)f1T⊥1x∝∫d2b⊥I¯b⊥∂∂byEx,b⊥2,for a nucleon polarized in the transverse x direction. Here E(x,b⊥2)≡∫(d2Δ⊥/(2π)2)e-ib⊥·Δ⊥E(x,0,-Δ⊥2) and I¯(b⊥) is called the lensing function. This relation has been obtained in models [9, 105, 106, 109]. It allows relating 〈k⊥q〉 to the anomalous magnetic moment κq associated with the quark q:(14)∫dx∫d2b⊥Eqx,b⊥2=κq,albeitin a model-dependent and (due to the different integrals involved) only qualitative way. This relation does confirm the expectations for the relative signs between the u and d Sivers functions and has been used in [110] to fit SIDIS data for the Sivers effect with the integral constrained by the anomalous magnetic moments. Interestingly, this led to a new estimation of the quark total angular momentum which turned out to be in agreement with most common GPD extractions [111–115]. The relation between 〈k⊥q〉 and κq is also at the heart of the argument by Brodsky and Gardner of why a gluon Sivers function is expected to be small. Using that κup=2κp+κn=1.673 and κdp=2κn+κp=-2.033, one sees the opposite signs reflected, but since (κp+κn)/2=0.06≪κp/n≈1.8-1.9 this suggests that there is little room for gluon contributions [60]. If (κup+κdp)/2 is taken as a measure for κgp, the latter is about 10% smaller than κqp. This would suggest that κgp is of order 1/Nc2 rather than 1/Nc, which in turn would suggest a similar additional 1/Nc suppression for the gluon Sivers function. Clearly there are various (strong!) assumptions going into this type of argument, such that the conclusion can certainly not be taken at face value. Apart from the assumptions on the relation to the gluon Sivers function, it is not clear that one can use very low-energy quantities to deduce something about the size of the gluon contributions at energies around or above 1 GeV to begin with.Burkardt derived a further constraint on the fully integrated quantity〈k⊥a〉, nowadays referred to as the Burkardt sum rule (BSR) [101], stating that the total transverse momentum of all partons in a transversely polarized proton must vanish [116]:(15)k⊥=∑a=q,q¯,gk⊥a=0.In terms of the Sivers function, the BSR takes the form [59](16)∑a=q,q¯,g∫dxf1T⊥1ax=0.Its validity has been checked explicitly in a diquark spectator model in [117]. The fits to SIDIS data from [108] at the scale Q2=2.4 GeV2 almost saturate the BSR already with the u and d quark contributions alone: (17)k⊥u=96-28+60MeV,k⊥d=-113-51+45MeV.The contributions of the sea quarks are all small and together allow the following range for the gluon contribution:(18)-10≤k⊥g≤48MeV.This means there is certainly still room for a 30% contribution from gluons with respect to the valence quarks. Of course, it should be emphasized that these values were obtained under assumptions on the k⊥ dependence, the absence of nodes, and extrapolations outside the kinematic region accessed by the SIDIS experiments.The derivation of the BSR by Burkardt [101] considers gauge links as they appear in SIDIS and involves a gluon correlator containing the antisymmetric fabc structure constant of SU(3). As shown in [94, 95], there is also a gluon correlator with the symmetric dabc structure constant. As a consequence, inclusion of gauge links in the operator definition of TMD distributions gives rise to two distinct gluon Sivers functions, f1T⊥g(f) and f1T⊥g(d) (corresponding to the (Ac) label used in [11]). However, the BSR essentially expresses transverse momentum conservation. Since the momentum operator in QCD is C-even, only the gluon Sivers function f1T⊥g(f) which is associated with a C-even operator is constrained by the BSR. The gluon Sivers function f1T⊥g(d), which is associated with a C-odd operator, is not expected to satisfy a BSR where quark and gluon contributions cancel each other. Judging from the small-x behavior of TG(d)(x,x) expected from its evolution equation, the integral of f1T⊥g(d) over x may even not converge.


## 7. Conclusions

In summary, no hard constraints on the size of the gluon Sivers function exist apart from the positivity bound, although the theoretical expectation from largeNc considerations (expected to hold approximately for not too small x) favors a 30% gluon to quark Sivers ratio, which is still completely allowed by all SIDIS data, including the COMPASS deuteron data. It may turn out that the ratio is smaller, but much smaller than 10% may in turn be considered unnaturally small. Strictly speaking, no direct conclusion about the size of the gluon Sivers function can be drawn from AN data. Like p↑p→πX, many other processes suggested in the literature to probe the gluon Sivers function actually deal with collinear factorization and as such they are sensitive to complicated linear combinations of quark-gluon and trigluon Qiu-Sterman functions (and chiral-odd and fragmentation function versions of them) rather than to Sivers functions directly. Inferring constraints on the gluon Sivers function, even on its large transverse momentum tail, must therefore be done with much care. In the Generalized Parton Model description of AN at midrapidity, theeffective gluon Sivers function is currently still allowed to be 30% of the up quark Sivers function, despite the smallness of the asymmetry. Other suggested processes that in principle probe TMDs may suffer from TMD factorization breaking contributions and any results on the gluon Sivers function from measurements of such processes cannot be trusted. This applies, for instance, to the process p↑p→jetjetX.The most promising processes that directly give access to the gluon Sivers effect arep↑p→γjetX, p↑p→J/ψγX, and ep↑→e′cc¯X. The first process can be studied at RHIC and at a polarized fixed-target experiment at LHC (AFTER@LHC), the second process also at AFTER@LHC, and the third process at a possible future Electron-Ion Collider. Due to the inherent process dependence of TMDs, the gluon Sivers TMD probed is in principle different in these processes. They can be expressed in terms of two universal gluon Sivers functions that appear in different linear combinations in different processes. Extracting and comparing these universal functions is very interesting from a theoretical point of view. The fact that a difference can exist is a consequence of the non-Abelian nature of QCD. Both functions satisfy different theoretical constraints. Although TMD factorization is expected to hold for these processes, that has not been demonstrated yet to all orders. Apart from the process dependence, there is also the issue of modified soft factors to contend with still. Nevertheless, as far as experimentally demonstrating and measuring a gluon Sivers effect in transversely polarized protons, several complementary future possibilities exist, in which AFTER@LHC can play a very important role.



---

*Source: 371396-2015-10-04.xml*