# Transverse Single-Spin Asymmetries in Proton-Proton Collisions at the AFTER@LHC Experiment in a TMD Factorisation Scheme

Authors: M. Anselmino; U. D’Alesio; S. Melis
Journal: Advances in High Energy Physics
 (2015)

---

## Abstract

The inclusive large-pT production of a single pion, jet or direct photon, and Drell-Yan processes, are considered for proton-proton collisions in the kinematical range expected for the fixed-target experiment AFTER, proposed at LHC. For all these processes, predictions are given for the transverse single-spin asymmetry, AN, computed according to a Generalised Parton Model previously discussed in the literature and based on TMD factorisation. Comparisons with the results of a collinear twist-3 approach, recently presented, are made and discussed.


---

## Body

## 1. Introduction and Formalism

Transverse Single-Spin Asymmetries (TSSAs) have been abundantly observed in several inclusive proton-proton experiments for a long time; when reaching large enough energies andpT values, their understanding from basic quark-gluon QCD interactions is a difficult and fascinating task, which has always been one of the major challenges for QCD.In fact, large TSSAs cannot be generated by the hard elementary processes, because of helicity conservation (in the massless limit) typical of QED and QCD interactions; indeed, such asymmetries were expected to vanish at high energies. Their persisting must be related to nonperturbative properties of the nucleon structure, such as parton intrinsic and orbital motion. A true understanding of the origin of TSSAs would allow a deeper understanding of the nucleon structure.Since the 1990s two different, despite being somewhat related, approaches have attempted to tackle the problem. One is based on the collinear QCD factorisation scheme and involves as basic quantities, which can generate single-spin dependences, higher-twist quark-gluon-quark correlations in the nucleon as well as higher-twist fragmentation correlators. The second approach is based on a physical, although unproven, generalisation of the parton model, with the inclusion, in the factorisation scheme, of transverse momentum dependent partonic distribution and fragmentation functions (TMDs), which also can generate single-spin dependences. The twist-3 correlations are related to moments of some TMDs. We refer to [1–9], and references therein, for more detailed account of the two approaches and possible variations, with all relevant citations. Following [10], we denote by CT-3 the first approach while the second one is, as usual, denoted by GPM.In this paper we consider TSSAs at the proposed AFTER@LHC experiment, in which high-energy protons extracted from the LHC beam would collide on a (polarised) fixed target of protons, with high luminosity. For a description of the physics potentiality of this experiment see [11] and for the latest technical details and importance for TMD studies see, for example, [12]. Due to its features the AFTER@LHC is an ideal experiment to study and understand the origin of SSAs and, in general, the role of QCD interactions in high-energy hadronic collisions; AFTER@LHC would be a polarised fixed-target experiment with unprecedented high luminosity.We recall our formalism by considering the Transverse Single-Spin AsymmetryAN, measured in pp↑→hX inclusive reactions and defined as(1)AN=dσ↑-dσ↓dσ↑+dσ↓withdσ↑,↓≡Ehdσpp↑,↓→hXd3ph,where ↑, ↓ are opposite spin orientations perpendicular to the x-z scattering plane, in the pp↑ c.m. frame. We define the ↑ direction as the +y^-axis and the unpolarised proton is moving along the +z^-direction. In such a process the only large scale is the transverse momentum pT=(ph)x of the final hadron.In the GPMAN originates mainly from two spin and transverse momentum effects, one introduced by Sivers in the partonic distributions [13, 14] and one by Collins in the parton fragmentation process [15], being all the other effects strongly suppressed by azimuthal phase integrations [16]. According to the Sivers effect the number density of unpolarised quarks q (or gluons) with intrinsic transverse momentum k⊥ inside a transversely polarised proton p↑, with three-momentum P and spin polarisation vector S, can be written as(2)f^q/p↑x,k⊥=fq/px,k⊥+12ΔNfq/p↑x,k⊥S·P^×k^⊥,where x is the proton light-cone momentum fraction carried by the quark, fq/p(x,k⊥) is the unpolarised TMD (k⊥=k⊥), and ΔNfq/p↑(x,k⊥) is the Sivers function. P^=P/|P| and k^⊥=k⊥/k⊥ are unit vectors. Notice that the Sivers function is most often denoted as f1T⊥q(x,k⊥) [17]; this notation is related to ours by [18](3)ΔNfq/p↑x,k⊥=-2k⊥mpf1T⊥qx,k⊥,where mp is the proton mass.Similarly, according to the Collins effect the number density of unpolarised hadronsh with transverse momentum p⊥ resulting in the fragmentation of a transversely polarised quark q↑, with three-momentum q and spin polarisation vector Sq, can be written as(4)D^q↑/hz,p⊥=Dh/qz,p⊥+12ΔNDq↑/hz,p⊥Sq·q^×p^⊥,where z is the parton light-cone momentum fraction carried by the hadron, Dh/q(z,p⊥) is the unpolarised TMD (p⊥=|p⊥|), and ΔNDq↑/h(z,p⊥) is the Collins function. q^=q/|q| and p^⊥=p⊥/p⊥ are unit vectors. Notice that the Collins function is most often denoted as H1⊥q(z,p⊥) [17]; this notation is related to ours by [18](5)ΔNDh/q↑z,p⊥=2p⊥zMhH1⊥qz,p⊥,where Mh is the hadron mass.According to the GPM formalism [1, 2, 16], AN can then be written as(6)AN=dσ↑-dσ↓Sivers+dσ↑-dσ↓Collinsdσ↑+dσ↓·The Collins and Sivers contributions were recently studied, respectively, in [1] and [2], and are given by(7)
d
σ
↑
-
d
σ
↓
Sivers
=
∑
a
,
b
,
c
,
d
∫
d
x
a
d
x
b
d
z
16
π
2
x
a
x
b
z
2
s
d
2
k
⊥
a
d
2
k
⊥
b
d
3
p
⊥
δ
p
⊥
·
p
^
c
·
J
p
⊥
δ
s
^
+
t
^
+
u
^
Δ
N
f
a
/
p
↑
x
a
,
k
⊥
a
cos
⁡
ϕ
a
·
f
b
/
p
x
b
,
k
⊥
b
1
2
M
^
1
0
2
+
M
^
2
0
2
+
M
^
3
0
2
a
b
→
c
d
·
D
h
/
c
z
,
p
⊥
,
(8)
d
σ
↑
-
d
σ
↓
Collins
=
∑
q
a
,
b
,
q
c
,
d
∫
d
x
a
d
x
b
d
z
16
π
2
x
a
x
b
z
2
s
d
2
k
⊥
a
d
2
k
⊥
b
d
3
p
⊥
δ
p
⊥
·
p
^
c
J
p
⊥
δ
s
^
+
t
^
+
u
^
Δ
T
q
a
x
a
,
k
⊥
a
cos
⁡
ϕ
a
+
φ
1
-
φ
2
+
ϕ
π
H
f
b
/
p
x
b
,
k
⊥
b
M
^
1
0
M
^
2
0
q
a
b
→
q
c
d
·
Δ
N
D
h
/
q
c
↑
z
,
p
⊥
.For details and full explanation of the notations in the above equations we refer to [16] (where p⊥ is denoted as k⊥C). It suffices to notice here that J(p⊥) is a kinematical factor, which at O(p⊥/Eh) equals 1. The phase factor cos⁡(ϕa) in (7) originates directly from the k⊥ dependence of the Sivers distribution [S·(P^×k^⊥), (2)]. The (suppressing) phase factor cos⁡(ϕa+φ1-φ2+ϕπH) in (8) originates from the k⊥ dependence of the unintegrated transversity distribution ΔTq, the polarized elementary interaction, and the spin-p⊥ correlation in the Collins function. The explicit expressions of φ1, φ2, and ϕπH in terms of the integration variables can be found via (60)–(63) in [16] and (35)–(42) in [19].TheM^i0’s are the three independent hard scattering helicity amplitudes describing the lowest order QCD interactions. The sum of their moduli squared is related to the elementary unpolarised cross section dσ^ab→cd; that is,(9)dσ^ab→cddt^=116πs^212∑i=13M^i02.The explicit expressions of the combinations of M^i0’s which give the QCD dynamics in (7) and (8), can be found, for all possible elementary interactions, in [16] (see also [1] for a correction to one of the product of amplitudes). The QCD scale is chosen as Q=pT.The denominator of (1) or (6) is twice the unpolarised cross section and is given in our TMD factorisation by the same expression as in (7), where one simply replaces the factor ΔNfa/p↑cos⁡(ϕa) with 2fa/p.


## 2.AN for Single Pion, Jet, and Direct Photon Production

We present here our results forAN, (1), based on our GPM scheme, (6), (7), and (8). The TMDs which enter in these equations are those extracted from the analysis of Semi-Inclusive Deep Inelastic (SIDIS) and e+e- data [20–23], adopting simple factorised forms, which we recall here. For the unpolarised TMD partonic distributions and fragmentation functions we have, respectively,(10)fq/px,k⊥=fq/pxe-k⊥2/k⊥2πk⊥2k⊥2=0.25GeV2,(11)Dh/qz,p⊥=Dh/qze-p⊥2/p⊥2πp⊥2p⊥2=0.20GeV2.The Sivers function is parameterised as(12)ΔNfq/p↑x,k⊥=2NqSxfq/pxhk⊥e-k⊥2/k⊥2πk⊥2,where(13)NqSx=NqSxαq1-xβqαq+βqαq+βqαqαqβqβq,with |NqS|≤1, and(14)hk⊥=2ek⊥Me-k⊥2/M2.Similarly, the quark transversity distribution, ΔTq(x,k⊥), and the Collins fragmentation function, ΔNDh/q↑(z,p⊥), have been parametrized as follows:(15)ΔTqx,k⊥=12NqTxfq/px+Δqxe-k⊥2/k⊥2πk⊥2,ΔNDh/q↑z,p⊥=2NqCzDh/qzhp⊥e-p⊥2/p⊥2πp⊥2,where Δq(x) is the usual collinear quark helicity distribution,(16)NqTx=NqTxaq1-xbqaq+bqaq+bqaqaqaqbq,NqCz=NqCzγq1-zδqγq+δqγq+δqγqγqδqδq,with |NqT(C)|≤1, and (17)hp⊥=2ep⊥Mce-p⊥2/Mc2.All details concerning the motivations for such a choice, the values of the parameters, and their derivation can be found in [20–23]. We do not repeat them here, but in the caption of each figure we will give the corresponding references which allow fixing all necessary values.We present our results onAN for the process pp↑→πX at the expected AFTER@LHC energy (s=115 GeV) in Figures 1–3. Following [1, 2], our results are given for two possible choices of the SIDIS TMDs and are shown as function of pT at two fixed xF values (Figure 1), as function of xF at two fixed rapidity y values (Figure 2) and as function of rapidity at one fixed pT value (Figure 3). xF is the usual Feynman variable defined as xF=2pL/s where pL=(ph)z is the z-component of the final hadron momentum. Notice that, in our chosen reference frame, a forward production, with respect to the polarised proton, means negative values of xF. The uncertainty bands reflect the uncertainty in the determinations of the TMDs and are computed according to the procedure explained in the appendix of [21]. More information can be found in the figure captions.Figure 1
Our theoretical estimates forAN versus pT at s=115 GeV, xF=-0.2 (upper plots), and xF=-0.4 (lower plots) for inclusive π± and π0 production in pp↑→πX processes, computed according to (6)–(8) of the text. The contributions from the Sivers and the Collins effects are added together. The computation is performed adopting the Sivers and Collins functions of [20, 22] (SIDIS 1-KRE, left panels) and of [21, 23] (SIDIS 2-DSS, right panels). The overall statistical uncertainty band, also shown, is the envelope of the two independent statistical uncertainty bands obtained following the procedure described in Appendix  A of [21].Figure 2
Our theoretical estimates forAN versus xF at s=115 GeV, y=-1.5 (upper plots), and y=-3.0 (lower plots) for inclusive π± and π0 production in pp↑→πX processes, computed according to (6)–(8) of the text. The contributions from the Sivers and the Collins effects are added together. The computation is performed adopting the Sivers and Collins functions of [20, 22] (SIDIS 1-KRE, left panels) and of [21, 23] (SIDIS 2-DSS, right panels). The overall statistical uncertainty band, also shown, is the envelope of the two independent statistical uncertainty bands obtained following the procedure described in Appendix  A of [21].Figure 3
Our theoretical estimates forAN versus y at s=115 GeV and pT=3 GeV, for inclusive π± and π0 production in pp↑→πX processes, computed according to (6)–(8) of the text. The contributions from the Sivers and the Collins effects are added together. The computation is performed adopting the Sivers and Collins functions of [20, 22] (SIDIS 1-KRE, left panel) and of [21, 23] (SIDIS 2-DSS, right panel). The overall statistical uncertainty band, also shown, is the envelope of the two independent statistical uncertainty bands obtained following the procedure described in Appendix  A of [21].Notice that, for both our choices of the Sivers functions, the gluon Sivers distributions are taken to be vanishing, as suggested by data [21, 24]. Gluon channels contribute instead to the unpolarised cross sections, in the denominator of (1) or (6). For the unpolarised partonic distributions we adopt the GRV98LO PDF set [25] and for the fragmentation functions the DSS set from [26] and the Kretzer (KRE) set from [27].The analogous results for the single direct photon are shown in Figures4–6 (where xF=2(pjet)z/s) and those for the single jet production in Figures 7–9  (xF=2(pγ)z/s). In these cases, obviously, there is no fragmentation process and only the Sivers effect contributes to AN, with Dh/c(z,p⊥) simply replaced by δ(z-1)δ2(p⊥) in (7) (see [2] for further details). In our leading order treatment the jet coincides with a single final parton. Notice that for a jet production we have all the same QCD subprocesses which contribute to hadron production, while for a direct photon production the basic partonic subprocesses are the Compton scattering gq(q¯)→γq(q¯) and the annihilation process qq¯→γg [28].Figure 4
Our theoretical estimates forAN versus pT at s=115 GeV, xF=-0.2 (upper plots), and xF=-0.4 (lower plots) for inclusive photon production in pp↑→γX processes, computed according to (6) and (7) of the text. Only the Sivers effect contributes. The computation is performed adopting the Sivers functions of [20] (SIDIS 1, left panels) and of [21] (SIDIS 2, right panels). The overall statistical uncertainty band, also shown, is obtained following the procedure described in Appendix  A of [21].Figure 5
Our theoretical estimates forAN versus xF at s=115 GeV, y=-1.5 (upper plots), and y=-3.0 (lower plots) for inclusive photon production in pp↑→γX processes, computed according to (6) and (7) of the text. Only the Sivers effect contributes. The computation is performed adopting the Sivers functions of [20] (SIDIS 1, left panels) and of [21] (SIDIS 2, right panels). The overall statistical uncertainty band, also shown, is obtained following the procedure described in Appendix  A of [21].Figure 6
Our theoretical estimates forAN versus y at s=115 GeV and pT=3 GeV, for inclusive photon production in pp↑→γX processes, computed according to (6)–(8) of the text. Only the Sivers effects contributes. The computation is performed adopting the Sivers functions of [20] (SIDIS 1, left panel) and of [21] (SIDIS 2, right panel). The overall statistical uncertainty band, also shown, is the envelope of the two independent statistical uncertainty bands obtained following the procedure described in Appendix  A of [21].Figure 7
Our theoretical estimates forAN versus xF at s=115 GeV, xF=-0.2 (upper plots), and xF=-0.4 (lower plots) for inclusive single jet production in pp↑→jetX processes, computed according to (6) and (7) of the text. Only the Sivers effect contributes. The computation is performed adopting the Sivers functions of [20] (SIDIS 1, left panels) and of [21] (SIDIS 2, right panels). The overall statistical uncertainty band, also shown, is obtained following the procedure described in Appendix  A of [21].Figure 8
Our theoretical estimates forAN versus xF at s=115 GeV, y=-1.5 (upper plots), and y=-3.0 (lower plots) for inclusive single jet production in pp↑→jetX processes, computed according to (6) and (7) of the text. Only the Sivers effect contributes. The computation is performed adopting the Sivers functions of [20] (SIDIS 1, left panels) and of [21] (SIDIS 2, right panels). The overall statistical uncertainty band, also shown, is obtained following the procedure described in Appendix  A of [21].Figure 9
Our theoretical estimates forAN versus y at s=115 GeV and pT=3 GeV, for inclusive single jet production in pp↑→jetX processes, computed according to (6)–(8) of the text. Only the Sivers effect contributes. The computation is performed adopting the Sivers functions of [20] (SIDIS 1, left panel) and of [21] (SIDIS 2, right panel). The overall statistical uncertainty band, also shown, is the envelope of the two independent statistical uncertainty bands obtained following the procedure described in Appendix  A of [21].


## 3.AN for Drell-Yan Processes

Drell-Yan (D-Y) processes are expected to play a crucial role in our understanding of the origin, at the partonic level, of TSSAs. For such processes, like for SIDIS processes and contrary to single hadron production, the TMD factorisation has been proven to hold, so that there is a general consensus that the Sivers effect should be visible via TSSAs in D-Y [29–32]. The widely accepted interpretation of the QCD origin of TSSAs as final or initial state interactions of the scattering partons [33] leads to the conclusion that the Sivers function has opposite signs in SIDIS and D-Y processes [34], which remains to be seen.Predictions for SiversAN in D-Y and at different possible experiments were given in [35], which we follow here.In [35] predictions were given for the p↑p→l+l-X  D-Y process in the p↑-p c.m. frame, in which one observes the four-momentum q of the final l+l- pair. Notice that q2=M2 is the large scale in the process, while qT=|qT| is the small one. In order to collect data at all azimuthal angles, one defines the weighted spin asymmetry:(18)ANsin⁡ϕγ-ϕS≡∫02πdϕγdσ↑-dσ↓sin⁡ϕγ-ϕS1/2∫02πdϕγdσ↑+dσ↓,(19)=∫dϕγ∑qeq2∫d2k⊥1d2k⊥2δ2k⊥1+k⊥2-qTΔNfq/p↑x1,k⊥1fq¯/px2,k⊥2sin⁡ϕγ-ϕS∫dϕγ∑qeq2∫d2k⊥1d2k⊥2δ2k⊥1+k⊥2-qTfq/px1,k⊥1fq¯/px2,k⊥2,where ϕγ and ϕS are, respectively, the azimuthal angle of the l+l- pair and of the proton transverse spin and we have defined (see (2))(20)ΔNfq/p↑x,k⊥≡ΔNfq/p↑x,k⊥S·P^×k^⊥=f^q/p↑x,k⊥-f^q/p↓x,k⊥.Adopting for the unpolarised TMD and the Sivers function the same expressions as in (10) and (12)–(14) allows, at O(k⊥/M), an analytical integration of the numerator and denominator of (19), resulting in a simple expression for the asymmetry ANsin⁡(ϕγ-ϕS) [35].Notice that we consider here thepp↑→l+l-X D-Y process in the p-p↑ c.m. frame. For such a process the TSSA is given by [35](21)ANsin⁡ϕγ-ϕSpp↑⟶γ∗X;xF,M,qT=-ANsin⁡ϕγ-ϕSp↑p⟶γ∗X;-xF,M,qT.Our results for the Sivers asymmetryANsin⁡(ϕγ-ϕS) at AFTER@LHC, obtained following [35], (21) and using the SIDIS extracted Sivers functionreversed in sign, are shown in Figure 10. Further details can be found in the captions of these figures.Figure 10
Our theoretical estimates forAN(ϕγ-ϕS) in D-Y processes as expected at AFTER@LHC. Our results are presented as function of M (upper plots), xF (middle plots), and x of the quark inside the polarised proton, x↑ (lower plots). The other kinematical variables are either fixed or integrated, as indicated in each figure. They are computed according to [35] and (21), adopting the Sivers functions of [20] (SIDIS 1, left panels) and of [21] (SIDIS 2, right panels),reversed in sign. The overall statistical uncertainty band, also shown, is obtained following the procedure described in Appendix  A of [21].


## 4. Comments and Conclusions

Some final comments and further details might help in understanding the importance of the measurements of the TSSAs at AFTER@LHC:(i)
Most predictions given show clear asymmetries, sufficiently large as to be easily measurable, given the expected performance of AFTER@LHC [11]. The uncertainty bands reflect the uncertainty in the extraction of the Sivers and transversity functions from SIDIS data, which are focused on small and intermediate x values (x≲0.3); in fact the bands grow larger at larger values of |xF|.(ii)
The values ofAN found for pion production can be as large as 10% for π±, while they are smaller for π0. They result from the sum of the Sivers and the Collins effects. The relative importance of the two contributions varies according to the kinematical regions and the set of distributions and fragmentation functions adopted. As a tendency, the contribution from the Sivers effect is larger than the Collins contribution with the SIDIS 1-KRE set, while the opposite is true for the SIDIS 2-DSS set.The values found here are in agreement, both in sign and qualitative magnitude, with the values found in [10] within the collinear twist-3 (CT-3) approach.(iii)
The results for single photon production are interesting; they isolate the Sivers effect and our predictions show that they can reach values of about 5%, with a reduced uncertainty band. We find positive values ofAN as the relative weight of the quark charges leads to a dominance of the u quark and the Sivers functions ΔNfu/p↑ which is positive [20, 21].Our results, obtained within the GPM, have a similar magnitude to those obtained in [3, 10], within the CT-3 approach, but havean opposite sign. Thus, a measurement of AN for a single photon production, despite being difficult, would clearly discriminate between the two approaches.(iv)
The values ofAN for single jet production, which might be interesting as they also have no contribution from the Collins effect, turn out to be very small and compatible with zero, due to a strong cancellation between the u and d quark contributions. The same result is found in [10].(v)
A measurement ofANsin⁡(ϕγ-ϕS) in D-Y processes at AFTER@LHC is the most interesting one. In such a case the TMD factorisation has been shown to be valid and the Sivers asymmetry should show the expected sign change with respect to SIDIS processes [33, 34]. Our computations, Figure 10, predict a clear asymmetry, which can be as sizeable as 10%, with a definite sign, even within the uncertainty band.Both the results of [10] and the results of this paper obtain solid non negligible values for the TSSA AN measurable at the AFTER@LHC experiment. The two sets of results are based on different approaches, respectively, the CT-3 and the GPM factorisation schemes. While the magnitude of AN is very similar in the two cases, the signs can be different, in particular, the TSSA for a direct photon production, pp↑→γX, has opposite signs in the two schemes.In this paper we have also considered azimuthal asymmetries in polarised D-Y processes, related to the Sivers effect. As explained above, in this case, due to the presence of a large and a small scale, like in SIDIS, the TMD factorisation is valid, with the expectation of an opposite sign of the Sivers function in SIDIS and D-Y processes. Also this prediction can be checked at AFTER@LHC.



---

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