#### Abstract

We present results for transverse single-spin asymmetries in proton-proton collisions at kinematics relevant for AFTER, a proposed fixed-target experiment at the Large Hadron Collider. These include predictions for pion, jet, and direct photon production from analytical formulas already available in the literature. We also discuss specific measurements that will benefit from the higher luminosity of AFTER, which could help resolve an almost 40-year puzzle of what causes transverse single-spin asymmetries in proton-proton collisions.

#### 1. Introduction

Transverse single-spin asymmetries (TSSAs), denoted by , have been a fundamental observable since the mid-1970s to test perturbative quantum chromodynamics (pQCD). Such measurements were first conducted at FermiLab, where large effects were found in [1]. These results contradicted the naïve collinear parton model, which said that should be extremely small [2], and doubts were raised as to whether pQCD can actually describe these reactions [2]. However, in the 1980s it was shown that if one went beyond the parton model and included collinear twist-3 (CT3) quark-gluon-quark correlations in the nucleon, substantial TSSAs could be generated [3, 4]. In the 1990s this CT3 approach was worked out in more detail for proton-proton collisions, first for direct photon production [5–7] and then for pion production [8]. Over the last decade, several other analyses furthered the development of this formalism; see [9–19] and references therein. During the same time, another mechanism was also put forth to explain TSSAs in proton-proton collisions. This approach involves the Sivers [20, 21], Collins [22], and Boer-Mulders [23, 24] transverse momentum dependent (TMD) functions and became known as the Generalized Parton Model (GPM); see [25–29] and references therein. (We mention that since most likely a rigorous factorization formula involving TMD functions does not hold for single-inclusive processes (which have only one scale), the GPM can only be considered a phenomenological model.) In addition to all of this theoretical work, many experimental measurements of have been performed at proton-(anti)proton accelerators [30–42]. Most of the experimental data in the more negative region has come in the form of light-hadron asymmetries , for example, , , , with the exception of the jet asymmetry measured a few years ago at the Relativistic Heavy Ion Collider (RHIC) by the DY Collaboration [40]. (Throughout the paper we will use the convention , where is the momentum of the outgoing particle, and the transversely polarized proton moves along the -axis. That is, means large momentum fractions of the parton probed inside the transversely polarized proton. This setup causes to be opposite in sign to the one used in collider experiments (like those at RHIC).) Plans are also in place to measure the direct photon asymmetry at RHIC by both the PHENIX Collaboration and the STAR Collaboration [43–45].

Although much progress has been made in understanding TSSAs, there is not a definitive answer on what their origin might be. In the CT3 approach it was assumed for many years that a soft-gluon pole (SGP) chiral-even quark-gluon-quark () matrix element called the Qiu-Sterman (QS) function was the main cause of [8, 11]. However, this led to a so-called “sign-mismatch” between the QS function and the TMD Sivers function extracted from semi-inclusive deep-inelastic scattering (SIDIS) [46]. This issue could not be resolved through more flexible parameterizations of the Sivers function [47]. Moreover, the authors of [48] argued, by looking at data on the target TSSA in inclusive DIS [49, 50], that cannot be the main source of . This observation led us last year in [51] to analyze by including not only the QS function but also the fragmentation mechanism, whose analytical formula was first fully derived in [18] (the so-called “derivative term” was first computed in [52]) (see also [17, 53, 54] for fragmentation terms in other processes). We found in this situation for the first time in pQCD that one can fit all RHIC high transverse momentum pion data very well without any sign-mismatch issue. Furthermore, we showed that a simultaneous description of TSSAs in , SIDIS, and is possible. Nevertheless, more work must be done to confirm/refute this explanation and its predictions. We mention that, in the GPM, one cannot draw a definitive conclusion as to whether the Sivers or Collins mechanism is the main cause of [28, 29]. (In principle the Boer-Mulders function and gluon Sivers function can also contribute to the GPM formalism, but these pieces have not been analyzed in the literature.) This is due to the theoretical error bands being too large, since the associated TMD functions are mostly unconstrained in the large- regime covered by the data [28, 29]. For a detailed discussion of the GPM formalism and its predictions for the AFTER experiment, see [55].

In addition, in order to have a complete knowledge of TSSAs, it is important to have a “clean” extraction of the QS function from observables like and that do not have any fragmentation contributions. (We will ignore photons coming from fragmentation [56], which can be largely suppressed by using isolation cuts.) (For recent analyses of in collisions, see [57, 58].) This is necessary in order to help resolve the sign-mismatch issue and better understand the role of rescattering effects in the nucleon. The jet asymmetry has been studied in [11, 46, 59, 60] and the direct photon asymmetry has been investigated in [6, 11, 12, 56, 59–64]. It is important to point out that other contributions to and exist besides the one from the (SGP chiral-even) QS function. These other pieces include (i) soft-fermion pole (SFP) chiral-even functions, (ii) SGP and SFP chiral-odd functions, and (iii) SGP trigluon functions. For the numerical analyses in [59, 64] show that (i) is negligible for while the study in [64] draws a similar conclusion for (ii) as does the work in [62] for (iii). That is, for , the QS function dominates the asymmetry. We mention that, at present in the GPM, is predicted to have the* opposite sign* to that from the CT3 approach [29]. Therefore, as was emphasized in [64], this observable could allow us for the first time to clearly distinguish between the two frameworks as well as learn about the process dependence of the Sivers function [65], which is a feature of this nonperturbative object that is crucial to our current understanding of TMD functions.

For the conclusions as to which piece dominates are not as clear. The study in [66, 67] provides evidence that (ii) should be small in the whole -region. The work in [59] shows that the same is most likely true for (i) but that analysis suffers from the sign-mismatch issue. Also, in [19] there is an indication that (iii) could be significant. Therefore, it will be necessary to reassess the impact of (i) and (iii) on . Nevertheless, one can gain insight into these other terms by looking at the contribution from the QS function and comparing it with data.

Given the open issues that still remain, it is an opportune time for the Large Hadron Collider (LHC) to produce data on TSSAs in proton-proton collisions via the AFTER experiment. These measurements will not only add to the data from FermiLab, AGS, and RHIC but also, through the high luminosity of the experiment [68, 69], probe certain features that remain ambiguous. For example, the behavior of at large pion transverse momentum appears to fall off very slowly (or is even flat), a feature which the theory says should persist to high [29, 51, 55, 59] (see also [70] in the context of production). However, the data from RHIC [71] has too large error bars (or not enough statistics) in this high- region to ascertain whether or not this is true. Also, measured by DY [40] has large error bars as becomes more negative, which makes it difficult to determine whether or not the QS function alone can describe that data. Moreover, as previously mentioned, has never been measured before, yet it could be a tremendous opportunity to learn about the process dependence of the Sivers function and distinguish between the CT3 and GPM frameworks. Already, PHENIX and STAR plan to carry out such experiments [43–45].

Therefore, in this paper we give predictions within the CT3 formalism for , , and at AFTER@LHC kinematics. (For related work on charmonium and bottomonium production we refer to [72, 73].) Since the relevant analytical formulas already exist within the literature, in Section 2 we focus on the phenomenology and refer the reader to the appropriate papers on the underlying theory. These numerical results are summarized in Section 3, and there we highlight again how AFTER can offer unique insight into TSSAs in proton-proton collisions, which is a truly fundamental observable to test pQCD at higher twist.

#### 2. Pion, Photon, and Jet TSSAs at AFTER

We start first with , where we follow our numerical work in [51]. (We also refer the reader to [11, 13, 18] for more formal discussions of the relevant analytical formulas.) There we took into account the contribution from the QS function and the fragmentation term. The former has a model-independent relation to the Sivers function [74], while the latter involves three nonperturbative CT3 fragmentation functions (FFs): , , and . Of these, has a model-independent connection to the Collins function [18, 52] and can be written in terms of the other two through a QCD equation-of-motion relation [18]. In Figures 1 and 2 we provide predictions for neutral and charged pion production at AFTER based on our fit in [51]. One sees in Figure 1 from versus that the magnitude of the asymmetry can be anywhere from ~ and from versus that it increases with more negative (center-of-mass) rapidity . (Recall the relation between , , and : , so versus ( versus ) at fixed () implies a running in ().) One also notices that turns over at more negative values, which was also observed in some of the STAR data [33, 34, 39]. In Figure 2, where we show versus , one sees that the asymmetry is flat or falls off very slowly as increases, a feature that had also been measured by STAR [71]. It will be important to establish with more precision if this flatness persists at higher- values, say 12–15 GeV, and AFTER, with its much higher luminosity, will be in a position to make such a measurement.

**(a)**

**(b)**

**(c)**

**(a)**

**(b)**

**(c)**

We next look at and , which do not receive contributions from FFs. As we discussed in Section 1, the former may receive nonnegligible contributions from terms other than the QS function, while for the latter we recently showed in [64] that the QS function is the dominant piece to that asymmetry. (All of the analytical expressions for can be found in [64] while those for are determined simply by setting (the unpolarized FF) to in the equations for given in [11, 14, 19, 66, 67].) (We note the analytical formulas for the piece involving chiral-odd functions are new from [64], while those involving chiral-even functions were derived before in the literature, and the relevant references are cited therein.) However, given that the other pieces for are not reliably known, for that asymmetry we will only look at the contribution from the QS function using its relation to the Sivers function, while for we adopt our work in [64]. In Figures 3 and 4 we show results for jet and photon production at AFTER. We see that is very small, although we caution the reader that the Sivers function (which we use as input for the QS function) is mostly unconstrained in the large- region, and when this uncertainty is taken into account, one could obtain a measurable asymmetry [60]. Also, as we mentioned, there is the potential for (chiral-even) SFP and/or trigluon functions to make an impact. Therefore, in order to determine if the Sivers function alone can describe , along with the current data from DY, we need more precise data in the far backward region, which should be possible at AFTER. (We note that STAR has preliminary data on electromagnetic “jets” that could also be helpful [75].)

**(a)**

**(b)**

**(c)**

**(a)**

**(b)**

**(c)**

Unlike the jet asymmetry, could be on the order of at less negative and more negative (see Figures 3(a) and 3(b)) or smaller and less negative (see Figure 4). Both of these observations are consistent with the behavior of as a function of rapidity (see Figure 3(c)), where the asymmetry peaks at (with GeV), which corresponds to . Since the QS function is the dominant source of the asymmetry, we can have “clean” access to it. We state again that the GPM framework at present predicts to be* positive* [29]. Therefore, a clear nonzero signal for this observable would help to distinguish between the CT3 and GPM formalisms. However, we emphasize that should data contradict the predictions of the GPM, this does not invalidate the results obtained for TMD observables that are based on rigorous TMD factorization proofs. Also, since we use the Sivers function from SIDIS as our input for the QS function, we can learn about the predicted process dependence of the Sivers function.

#### 3. Summary and Outlook

In this paper we have discussed TSSAs in single-inclusive pion, jet, and photon production from proton-proton collisions, that is, , at kinematics relevant for the proposed AFTER@LHC experiment. These asymmetries have been fundamental observables to test pQCD at higher twist for close to 40 years, and much work has been performed on both the theoretical and experimental sides. Nevertheless, issues still remain as to the origin of these TSSAs, which makes a measurement of at the LHC via the AFTER experiment timely. For we have found that AFTER should expect (absolute) asymmetries on the order of as a function of and increasing as the rapidity becomes more negative. Also, the dependence of still falls off slowly and flattens out at high . For we predict a very small asymmetry, but we must remember that uncertainties in the Sivers function could allow for a measurable observable [60] and also that other contributions (like chiral-even SFP and trigluon) could make an impact. Lastly, for we expect asymmetries on the order of and decreasing with more negative and increasing . These are opposite in sign to the ones predicted from the GPM [29].

Even though these observables have been (or are planned to be) measured at RHIC, AFTER has the ability, through its much higher luminosity, to not only supplement the RHIC data but also provide important information on still unknown issues. For example, it will be key to determine if stays flat at higher-, say to 12–15 GeV, like theory predicts [29, 51, 55, 59, 70] and STAR has evidence for [71]. Also, higher statistics should allow for more precise measurements of at more negative , which will be necessary to determine if the QS function is the sole source of that asymmetry. Moreover, has never been measured before and provides the opportunity to clearly distinguish between the CT3 and GPM frameworks and learn about the process dependence of the Sivers function. Given the questions that remain as to the origin of TSSAs, which has been unresolved for almost 40 years, AFTER could provide valuable data on these observables.

#### Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

#### Acknowledgments

This work has been supported by the Grant-in-Aid for Scientific Research from the Japanese Society of Promotion of Science under Contract no. 26287040 (Y. Koike), the National Science Foundation under Contract no. PHY-1205942 (K. Kanazawa and A. Metz), and the RIKEN BNL Research Center (D. Pitonyak).