Nucleon Structure Physics by -Proton Collisions

Ming Xiong Liu

Physics Division, Los Alamos National Laboratory, Los Alamos, NM, 87545, USA

E-mail: [email protected] (Received April 23, 2019)

Today, the nucleons ( and ) make up most of the visible matter in the universe, from the nucleus to the galaxy. Understanding the fundamental constituents of the nucleon and the internal parton dynamics is of great interest of modern physics. However, describing the nucleon structure with the fundamental degree of freedom of and in the QCD framework remains as one of the most outstanding challenges in modern high energy nuclear and particle physics. In the last decades, largely due to the worldwide collaborative effort and advances in beam acceleration and target polarization technologies, tremendous progress has been made in the study of nucleon structures with polarized high energy leptons and hadron beams. The complementarity of lepton scattering and pure hadronic interaction is critical in the study of nucleon structure and QCD. In this presentation, I highlight a few recent achievements in the study of the nucleon structures with high energy (polarized) hadronic scatterings at the Relativistic Heavy Ion Collider (RHIC) at Brookhaven National Laboratory, at COMPASS at CERN as well as the SeaQuest experiments at Fermilab. Future prospects of new opportunities at above and other facilities worldwide are also discussed briefly. KEYWORDS: Proton, , spin, QCD

1. Introduction

Understanding the structure of the nucleon at the fundamental level with and ’s degree of freedom is one of the goals of modern physics. In the high-energy nucleon reactions, the nucleon can be described as a collection of partons, whose distributions inside the nucleon are generally universal and process independent, the so called parton distribution functions (PDFs). The reaction rate can be calculated in the perturbative QCD (pQCD) framework as a sum of the factorized individual parton-parton interactions. PDFs are fundamental properties of the nucleon. At the leading twist in the collinear factorization framework, there are three kinds of parton f(x, Q2) distribution functions: 1) the unpolarized parton distribution functions, AAAB8XicbVBNS8NAEJ3Ur1q/qh69LBahgpSkCnosevHYgv3ANpbNdtMu3WzC7kYsof/CiwdFvPpvvPlv3LQ5aOuDgcd7M8zM8yLOlLbtbyu3srq2vpHfLGxt7+zuFfcPWiqMJaFNEvJQdjysKGeCNjXTnHYiSXHgcdr2xjep336kUrFQ3OlJRN0ADwXzGcHaSPd++ekMNR6qp4V+sWRX7BnQMnEyUoIM9X7xqzcISRxQoQnHSnUdO9JugqVmhNNpoRcrGmEyxkPaNVTggCo3mV08RSdGGSA/lKaERjP190SCA6UmgWc6A6xHatFLxf+8bqz9KzdhIoo1FWS+yI850iFK30cDJinRfGIIJpKZWxEZYYmJNiGlITiLLy+TVrXinFeqjYtS7TqLIw9HcAxlcOASanALdWgCAQHP8ApvlrJerHfrY96as7KZQ/gD6/MHLYePSA== ; 2) the parton 2 h (x, Q2) f(x, Q ) AAAB83icbVBNS8NAEJ3Ur1q/qh69LBahgpSkFfRY9OKxBfsBbSyb7aZdutmE3Y1YQv+GFw+KePXPePPfuGlz0NYHA4/3ZpiZ50WcKW3b31ZubX1jcyu/XdjZ3ds/KB4etVUYS0JbJOSh7HpYUc4EbWmmOe1GkuLA47TjTW5Tv/NIpWKhuNfTiLoBHgnmM4K1kfrjgVN+ukDNh+p5YVAs2RV7DrRKnIyUIENjUPzqD0MSB1RowrFSPceOtJtgqRnhdFbox4pGmEzwiPYMFTigyk3mN8/QmVGGyA+lKaHRXP09keBAqWngmc4A67Fa9lLxP68Xa//aTZiIYk0FWSzyY450iNIA0JBJSjSfGoKJZOZWRMZYYqJNTGkIzvLLq6RdrTi1SrV5WarfZHHk4QROoQwOXEEd7qABLSAQwTO8wpsVWy/Wu/WxaM1Z2cwx/IH1+QNbHY/u helicity distribution functions, AAAB+nicbVBNS8NAEN34WetXqkcvi0WoICWpgh6LevDYgv2ANpbNdtIu3WzC7kYtsT/FiwdFvPpLvPlvTNoctPXBwOO9GWbmuSFnSlvWt7G0vLK6tp7byG9ube/smoW9pgoiSaFBAx7ItksUcCagoZnm0A4lEN/l0HJHV6nfugepWCBu9TgExycDwTxGiU6knlnoXgPXBHulxxNcv6sc53tm0SpbU+BFYmekiDLUeuZXtx/QyAehKSdKdWwr1E5MpGaUwyTfjRSEhI7IADoJFcQH5cTT0yf4KFH62AtkUkLjqfp7Iia+UmPfTTp9oodq3kvF/7xOpL0LJ2YijDQIOlvkRRzrAKc54D6TQDUfJ4RQyZJbMR0SSahO0kpDsOdfXiTNStk+LVfqZ8XqZRZHDh2gQ1RCNjpHVXSDaqiBKHpAz+gVvRlPxovxbnzMWpeMbGYf/YHx+QNq2pIl ; and 3) the transversity functions, 1 . Here x represents the parton light-cone momentum fraction, and Q2 denotes the momentum transfer under the reaction. In the three-dimensional partonic picture of the nucleon, one can go one step further beyond the above one-dimensional PDFs, describing the parton distribution functions in the (2+1) dimension, including the parton’s transverse momentum kT, in addition to the longitudinal momentum fraction x, thus they are Transverse-Momentum-Dependent (TMD) functions. At leading twist, there are eight TMD PDFs describing the partons’ longitudinal and transverse motion distributions and their correlations with the nucleon spin and momentum directions. All these functions can be probed in high-energy hadron-hadron or hadron-lepton scattering experiments, for recent reviews and some early work, see [1] and references therein.

With the advance of accelerator technology, the world’s first high energy polarized proton collider was built at the Brookhaven National Laboratory in New York, USA, in late 1990s, and has been in operation since 2001. Since then, tremendous progress has been made in studying the proton internal structures and QCD dynamics using longitudinally and transversely polarized proton-proton collisions at high energies at the Relativistic Heavy Ion Collider (RHIC). Fig. 1 shows the layout of the world’s first high-energy polarized proton accelerator complex at RHIC.

Figure 1. LH: The layout of the world’s first high energy polarized proton + proton collider at the Brookhaven National Laboratory in New York. The proton beam has been polarized up to about 60% either in longitudinal or transverse direction relative to the beam momentum direction, at center of mass energies from 62 to 510GeV. RH: High energy particles are produced in the polarized proton + proton collisions at RHIC. The polarized PDFs have be derived from the measured spin- dependent cross-sections in the pQCD framework. QCD is the theory of strong interactions. Particle production in the high energy polarized proton + proton collisions can be calculated as a sum of the product of the universal PDFs, f(x1) and f(x2) of the incoming partons and the partonic scattering cross sections that can be calculated in the perturbative QCD framework:

x +x h +h +X f(x ) f(x )ˆ 1 2! 1 2 AAACPHicbVA9TyMxEPXyGcJX7ihpLCIkEFK0i5CgRHfNlUEQiJQN0azjzVrY65U9Cxet8sOuuR9BR0VDAUK01DibFHw9yfLzmzfyzIsyKSz6/p03Mzs3v7BYWaour6yurdd+/Dy3OjeMt5iW2rQjsFyKlLdQoOTtzHBQkeQX0dXvcf3imhsrdHqGw4x3FQxSEQsG6KRe7TS0YqCAukvReOdvL9iloUahuC2f+7thAlhMXKPLwhnoHnU6DY0YJAjG6BualGri1D3aHvVqdb/hl6BfSTAldTJFs1e7Dfua5YqnyCRY2wn8DLsFGBRM8lE1zC3PgF3BgHccTcFN1y3K5Ud02yl9GmvjToq0VN93FKCsHarIORVgYj/XxuJ3tU6O8VG3EGmWI0/Z5KM4lxQ1HSdJ+8JwhnLoCDAj3KyUJWCAocu76kIIPq/8lZzvNwK/EZwc1I9/TeOokE2yRXZIQA7JMflDmqRFGPlH7skjefL+ew/es/cysc54054N8gHe6xtj1qu7 1 2 ⇠ ⌦ Using above theoretical tool, one can derive the universal parton distribution function f(x) with

the experimentally measured cross section, AAAB7nicbVBNSwMxEJ34WetX1aOXYBE8ld0q6LHoxWMF+wHtUrJptg1NskuSFcrSH+HFgyJe/T3e/Ddm2z1o64OBx3szzMwLE8GN9bxvtLa+sbm1Xdop7+7tHxxWjo7bJk41ZS0ai1h3Q2KY4Iq1LLeCdRPNiAwF64STu9zvPDFteKwe7TRhgSQjxSNOiXVSp2/4SJLyoFL1at4ceJX4BalCgeag8tUfxjSVTFkqiDE930tskBFtORVsVu6nhiWETsiI9RxVRDITZPNzZ/jcKUMcxdqVsniu/p7IiDRmKkPXKYkdm2UvF//zeqmNboKMqyS1TNHFoigV2MY4/x0PuWbUiqkjhGrubsV0TDSh1iWUh+Avv7xK2vWaf1mrP1xVG7dFHCU4hTO4AB+uoQH30IQWUJjAM7zCG0rQC3pHH4vWNVTMnMAfoM8f1DCPOw== , and the theoretically calculated partonic cross- section, AAAB9HicbVBNS8NAEN34WetX1aOXxSJ4KkkV9Fj04rGC/YAmlM120y7d3cTdSaGE/g4vHhTx6o/x5r9x0+agrQ8GHu/NMDMvTAQ34Lrfztr6xubWdmmnvLu3f3BYOTpumzjVlLVoLGLdDYlhgivWAg6CdRPNiAwF64Tju9zvTJg2PFaPME1YIMlQ8YhTAlYK/BGBzDd8KMms3K9U3Zo7B14lXkGqqECzX/nyBzFNJVNABTGm57kJBBnRwKlgs7KfGpYQOiZD1rNUEclMkM2PnuFzqwxwFGtbCvBc/T2REWnMVIa2UxIYmWUvF//zeilEN0HGVZICU3SxKEoFhhjnCeAB14yCmFpCqOb2VkxHRBMKNqc8BG/55VXSrte8y1r94arauC3iKKFTdIYukIeuUQPdoyZqIYqe0DN6RW/OxHlx3p2PReuaU8ycoD9wPn8ArWWSCA== ˆ.

In the following, I present a few highlights selected from recent experimental and theoretical studies of the (polarized) nucleon structures using high energy (polarized) hadronic scatterings at BNL, CERN and Fermilab.

2. The proton spin puzzle and the polarized parton distribution functions

The nucleon’s quark and gluon helicity distributions, Df(x, Q2), are universal property of the nucleon. They are derived from high-energy polarized nucleon-nucleon or polarized nucleon- lepton scattering data. Before 1987, it was believed through some naïve nucleon quark models that the proton’s spin is mostly coming from the sum of quarks’ spin. When EMC collaboration at CERN carried out the first measurement of quark’s spin contribution to the proton spin in the deep- inelastic polarized muon-proton scattering in 1987, the results socked the world – the sum of the quark’s spin can only account for about 30% of the total proton spin, known as the “spin crisis”, see [2] and references therein. Several new polarized DIS experiments were conducted following the EMC’s discovery and further confirmed the surprising results. Where does the rest of the proton spin come from?

In QCD, the proton helicity sum rule is given by [1],

AAAB8HicbVBNSwMxEM3Wr1q/qh69BIvgqexWQY9FBT1WsB/SLiWbZtvQJLsks0JZ+iu8eFDEqz/Hm//GbLsHbX0w8Hhvhpl5QSy4Adf9dgorq2vrG8XN0tb2zu5eef+gZaJEU9akkYh0JyCGCa5YEzgI1ok1IzIQrB2MrzO//cS04ZF6gEnMfEmGioecErDSY++GCSD4ttQvV9yqOwNeJl5OKihHo1/+6g0imkimgApiTNdzY/BTooFTwaalXmJYTOiYDFnXUkUkM346O3iKT6wywGGkbSnAM/X3REqkMRMZ2E5JYGQWvUz8z+smEF76KVdxAkzR+aIwERginH2PB1wzCmJiCaGa21sxHRFNKNiMshC8xZeXSatW9c6qtfvzSv0qj6OIjtAxOkUeukB1dIcaqIkokugZvaI3RzsvzrvzMW8tOPnMIfoD5/MHg1aPjw== where ⌃AAAB+HicbVBNS8NAEN34WetHox69LBbBU0mqoMeiHjxWtB/QhLLZbtqlm03YnQg19Jd48aCIV3+KN/+NmzYHbX0w8Hhvhpl5QSK4Bsf5tlZW19Y3Nktb5e2d3b2KvX/Q1nGqKGvRWMSqGxDNBJesBRwE6yaKkSgQrBOMr3O/88iU5rF8gEnC/IgMJQ85JWCkvl3xbpgAgr17PowILvftqlNzZsDLxC1IFRVo9u0vbxDTNGISqCBa91wnAT8jCjgVbFr2Us0SQsdkyHqGShIx7Wezw6f4xCgDHMbKlAQ8U39PZCTSehIFpjMiMNKLXi7+5/VSCC/9jMskBSbpfFGYCgwxzlPAA64YBTExhFDFza2YjogiFExWeQju4svLpF2vuWe1+t15tXFVxFFCR+gYnSIXXaAGukVN1EIUpegZvaI368l6sd6tj3nrilXMHKI/sD5/ADMpkiI= denotes the sum of the quark spin contribution, and G is the gluon’s spin contribution, and L + L represents the contributions from quarks’ and gluons’ orbital motion

AAAB8XicbVBNS8NAEJ3Ur1q/qh69LBZBEEpSBT0WvXjooYL9wDaEzXbTLt1s4u5GKKH/wosHRbz6b7z5b9y0OWjrg4HHezPMzPNjzpS27W+rsLK6tr5R3Cxtbe/s7pX3D9oqSiShLRLxSHZ9rChngrY005x2Y0lx6HPa8cc3md95olKxSNzrSUzdEA8FCxjB2kgPDe8RnaGGNyx55YpdtWdAy8TJSQVyNL3yV38QkSSkQhOOleo5dqzdFEvNCKfTUj9RNMZkjIe0Z6jAIVVuOrt4ik6MMkBBJE0JjWbq74kUh0pNQt90hliP1KKXif95vUQHV27KRJxoKsh8UZBwpCOUvY8GTFKi+cQQTCQztyIywhITbULKQnAWX14m7VrVOa/W7i4q9es8jiIcwTGcggOXUIdbaEILCAh4hld4s5T1Yr1bH/PWgpXPHMIfWJ8/iiWPhQ== q g inside the proton. Could gluon and sea-quark’s spin make up the rest of the missing proton spin? It was highly motivated theoretically, see [2] and references therein. Since DIS is not sensitive to the (polarized) gluon and flavor-identified sea quark distributions, a new experimental tool was highly desired. This led to the development and construction of the world’s first polarized proton collider at the Relativistic Heavy Ion Collider at BNL, to directly probe the gluon and sea quark’s spin contributions to the proton spin.

2.1 Gluon polarization inside the proton The polarized gluon distribution function (more exactly, the gluon helicity distribution function) is a fundamental quantity of the proton properties and has long been considered as a primary candidate for the missing piece required to solve the proton spin puzzle. RHIC has provided a great advantage over previous DIS by directly probing the gluon distribution at the leading order in the high energy QCD process. The double longitudinal spin asymmetry ALL of a hard process, like the inclusive high-energy jet or high pT hadron production in p+p collisions, is g(x) directly coupled to the polarized gluon distribution function AAAB8nicbVBNSwMxEM3Wr1q/qh69BItQL2W3Cnos6sFjBVsL26Vk02wbmk2WZFYsS3+GFw+KePXXePPfmLZ70NYHA4/3ZpiZFyaCG3Ddb6ewsrq2vlHcLG1t7+zulfcP2kalmrIWVULpTkgME1yyFnAQrJNoRuJQsIdwdD31Hx6ZNlzJexgnLIjJQPKIUwJW8rs3TADBg+rTaa9ccWvuDHiZeDmpoBzNXvmr21c0jZkEKogxvucmEGREA6eCTUrd1LCE0BEZMN9SSWJmgmx28gSfWKWPI6VtScAz9fdERmJjxnFoO2MCQ7PoTcX/PD+F6DLIuExSYJLOF0WpwKDw9H/c55pREGNLCNXc3orpkGhCwaZUsiF4iy8vk3a95p3V6nfnlcZVHkcRHaFjVEUeukANdIuaqIUoUugZvaI3B5wX5935mLcWnHzmEP2B8/kDJdyQgg== at the leading order in pQCD, thus providing the ultimate sensitivity to the polarized gluon distribution.

Evidences of the first definite non-zero gluon polarization came from the precision

measurements of double longitudinal spin asymmetry AAAAB7XicbVA9SwNBEJ2LXzF+RS1tFoNgFe5iQMuojUWKCOYDkiPsbfaSNXu7x+6eEI78BxsLRWz9P3b+GzfJFZr4YODx3gwz84KYM21c99vJra1vbG7ltws7u3v7B8XDo5aWiSK0SSSXqhNgTTkTtGmY4bQTK4qjgNN2ML6d+e0nqjST4sFMYupHeChYyAg2Vmpd99N6fdovltyyOwdaJV5GSpCh0S9+9QaSJBEVhnCsdddzY+OnWBlGOJ0WeommMSZjPKRdSwWOqPbT+bVTdGaVAQqlsiUMmqu/J1IcaT2JAtsZYTPSy95M/M/rJia88lMm4sRQQRaLwoQjI9HsdTRgihLDJ5Zgopi9FZERVpgYG1DBhuAtv7xKWpWyd1Gu3FdLtZssjjycwCmcgweXUIM7aEATCDzCM7zCmyOdF+fd+Vi05pxs5hj+wPn8AUBBjuo= LL from the year 2009 high statistic RHIC 0 data, in the inclusive jet production (STAR) and the inclusive ⇡AAAB7HicbVBNS8NAEJ3Ur1q/qh69LBbBU0mqoMeiF48VTFtoY9lsN+3SzSbsToRS+hu8eFDEqz/Im//GbZuDtj4YeLw3w8y8MJXCoOt+O4W19Y3NreJ2aWd3b/+gfHjUNEmmGfdZIhPdDqnhUijuo0DJ26nmNA4lb4Wj25nfeuLaiEQ94DjlQUwHSkSCUbSS303Fo9srV9yqOwdZJV5OKpCj0St/dfsJy2KukElqTMdzUwwmVKNgkk9L3czwlLIRHfCOpYrG3AST+bFTcmaVPokSbUshmau/JyY0NmYch7Yzpjg0y95M/M/rZBhdBxOh0gy5YotFUSYJJmT2OekLzRnKsSWUaWFvJWxINWVo8ynZELzll1dJs1b1Lqq1+8tK/SaPowgncArn4MEV1OEOGuADAwHP8ApvjnJenHfnY9FacPKZY/gD5/MHdOyOcw== production (PHENIX), Fig. 2.

STAR: PRL 115, 092002 (2015) PHENIX: PRD 90, 012007 (2014)

Figure 2. First precision data from RHIC show the first evidence of non-zero gluon polarization in the proton.

Several theory groups carried out new next-to-leading order (NLO) global QCD analyses to include the latest RHIC data and extracted the gluon spin contribution. Fig. 3 (LH) shows one G =0.2+0.06 example from DSSV global fit, AAACB3icbVDLSsNAFJ3UV62vqEtBBosgiCGpYt0IRQVdVrAPaGOYTKft0MkkzEyEErJz46+4caGIW3/BnX/jtM1CWw8M93DOvdy5x48Ylcq2v43c3PzC4lJ+ubCyura+YW5u1WUYC0xqOGShaPpIEkY5qSmqGGlGgqDAZ6ThDy5HfuOBCElDfqeGEXED1OO0SzFSWvLM3fYVYQrBa3gObat0nxzaln2aesmRruXUM4u6jgFniZORIshQ9cyvdifEcUC4wgxJ2XLsSLkJEopiRtJCO5YkQniAeqSlKUcBkW4yviOF+1rpwG4o9OMKjtXfEwkKpBwGvu4MkOrLaW8k/ue1YtU9cxPKo1gRjieLujGDKoSjUGCHCoIVG2qCsKD6rxD3kUBY6egKOgRn+uRZUi9ZzrFVuj0pVi6yOPJgB+yBA+CAMqiAG1AFNYDBI3gGr+DNeDJejHfjY9KaM7KZbfAHxucP5Q+WIw== 0 . 07 , in the x range from 0.05 to 0.2 that is covered by the RHIC data, and Fig. 3 (RH) shows how the future RHIC and EIC data will further constrain the polarized gluon contribution to the proton spin. For details, see [3] and references therein.

PRL 113, 012001

Figure 3. LH: DSSV NLO global fit indicates non-zero gluon spin contribution in the x range from 0.05 to 0.2. RH: Projected sensitivity to the gluon spin contribution from the future RHIC and EIC data [3].

2.2 Sea quark polarization For the first, flavor identified sea quark polarized PDFs are directly probed at leading order using the weak boson W+/- productions in longitudinally polarized p+p collisions at RHIC. Although DIS data have provided high precision determination of the quark polarized PDFs at the intermediate x-range, flavor identified sea-quark polarization distributions are poorly constrained since DIS can’t directly access the quark flavor and only derived using additional not-well determined parton fragmentation functions. With center of mass energy as high as 510GeV, RHIC has provided a new unique experimental tool to directly determine the flavor identified sea quark polarization through the measurements of spin-dependent W+/- production in the longitudinally polarized p+p collisions.

At RHIC energy (500 or 510GeV in the center of mass), weak bosons are predominantly produced through the quark and anti-quark annihilation process, shown in Fig. 4. The weak coupling of W+/- to left-handed quarks and right-handed anti-quarks directly selects the incoming quarks’ helicity, and the charge of the produced W boson naturally determines the incoming quarks’ flavor,

Figure 4. At RHIC energy, weak boson W's are predominantly produced through quark and anti-quark annihilation. q(x ) q¯ (x ) +/-

AAAB9HicbVDLSgNBEOyNrxhfUY9ehgQxIoTdeNBj0IvHCOYB2SXMTibJkNnZzcxscFnyF4IXD4p49WO85W+cPA4aLWgoqrrp7vIjzpS27amVWVvf2NzKbud2dvf2D/KHRw0VxpLQOgl5KFs+VpQzQeuaaU5bkaQ48Dlt+sPbmd8cU6lYKB50ElEvwH3BeoxgbSTP9bFMR2eT0mOnct7JF+2yPQf6S5wlKVYL7sXTtJrUOvkvtxuSOKBCE46Vajt2pL0US80Ip5OcGysaYTLEfdo2VOCAKi+dHz1Bp0bpol4oTQmN5urPiRQHSiWBbzoDrAdq1ZuJ/3ntWPeuvZSJKNZUkMWiXsyRDtEsAdRlkhLNE0MwkczcisgAS0y0ySlnQnBWX/5LGpWyc1mu3Js0bmCBLJxAAUrgwBVU4Q5qUAcCI3iGV3izxtaL9W59LFoz1nLmGH7B+vwG7/yUhg== where AAAB7XicbVDLSgNBEOyNrxhfUY9ehgQhIoTdeNBj0IvHCOYByRJmJ7PJmNmZdWZWXEL+wYMeFPHq/3jL3zh5HDRa0FBUddPdFcScaeO6Eyezsrq2vpHdzG1t7+zu5fcPGlomitA6kVyqVoA15UzQumGG01asKI4CTpvB8GrqNx+o0kyKW5PG1I9wX7CQEWys1LgvPXa9k26+6JbdGdBf4i1IsVronD5Pqmmtm//q9CRJIioM4VjrtufGxh9hZRjhdJzrJJrGmAxxn7YtFTii2h/Nrh2jY6v0UCiVLWHQTP05McKR1mkU2M4Im4Fe9qbif147MeGFP2IiTgwVZL4oTDgyEk1fRz2mKDE8tQQTxeytiAywwsTYgHI2BG/55b+kUSl7Z+XKjU3jEubIwhEUoAQenEMVrqEGdSBwB0/wCm+OdF6cd+dj3ppxFjOH8AvO5zfAWZGP 1 are 0 2 are the light quark and anti-quark that annihilate into W through the

+ u (x )+d¯ (x ) W d (x )+¯u (x ) W

AAACEnicbVDLSsNAFJ34rPUVdelmsAgthZJUQZdFNy5cVLEPaGKYTKbt0MkkzEzUEvoNbvwVNy4UcevKnX/jtM1CWw9cOJxzL/fe48eMSmVZ38bC4tLyympuLb++sbm1be7sNmWUCEwaOGKRaPtIEkY5aSiqGGnHgqDQZ6TlD87HfuuOCEkjfqOGMXFD1OO0SzFSWvLMUuJdFh88uwTL0PGRSIORd62Fagk6gvb6CgkR3cPWbdkzC1bFmgDOEzsjBZCh7plfThDhJCRcYYak7NhWrNwUCUUxI6O8k0gSIzxAPdLRlKOQSDedvDSCh1oJYDcSuriCE/X3RIpCKYehrztDpPpy1huL/3mdRHVP3ZTyOFGE4+mibsKgiuA4HxhQQbBiQ00QFlTfCnEfCYSVTjGvQ7BnX54nzWrFPqpUr44LtbMsjhzYBwegCGxwAmrgAtRBA2DwCJ7BK3gznowX4934mLYuGNnMHvgD4/MHJ9Gb6g== following processes, L 1 R 2 ! and AAACEnicbVDLSsNAFJ3UV62vqEs3g0VoEUtSBV0W3bhwUcU+oIlhMpm2QyeTMDNRS+g3uPFX3LhQxK0rd/6N0zYLrR64cDjnXu69x48Zlcqyvozc3PzC4lJ+ubCyura+YW5uNWWUCEwaOGKRaPtIEkY5aSiqGGnHgqDQZ6TlD87GfuuWCEkjfq2GMXFD1OO0SzFSWvLMcuBdlO49uwz3oeMjkSYj70oL1TJ0BO31FRIiuoOtmwPPLFoVawL4l9gZKYIMdc/8dIIIJyHhCjMkZce2YuWmSCiKGRkVnESSGOEB6pGOphyFRLrp5KUR3NNKALuR0MUVnKg/J1IUSjkMfd0ZItWXs95Y/M/rJKp74qaUx4kiHE8XdRMGVQTH+cCACoIVG2qCsKD6Voj7SCCsdIoFHYI9+/Jf0qxW7MNK9fKoWDvN4siDHbALSsAGx6AGzkEdNAAGD+AJvIBX49F4Nt6M92lrzshmtsEvGB/fKcmb7A== L 1 R 2 . At the leading order, the +/- ! longitudinal single spin asymmetry AL of W is given by [4],

where J is the polar angle of the out-going W+/- relative to the incoming polarized proton beam. STAR and PHENIX experiments at RHIC carried out the longitudinal single spin asymmetry measurements in the weak boson production in the polarized p+p collisions from 2011, 2012 and 2013 data sets, at center of mass energy of 510GeV, and beam polarization as high as +/- about 60%. Shown in Fig. 5 are the W single spin asymmetry AL vs lepton rapidity measured by the STAR (LH) and the PHENIX experiments (Center), also shown are the flavor-identified polarized u-bar and d-bar distributions from DSSV and NNPDF global fits (RH), see [4].

RHIC high energy W program has produced a substantial improvement of our understanding of the flavor-identified sea-quark polarization in the proton.

arXiv:1406.7122

Figure 5. Longitudinal single spin asymmetry vs lepton rapidity in W production from STAR(LH) and PHENIX ( Cent.) at RHIC; RH: flavor-identified sea quark helicity distributions extracted from global fits from NNPDF and DSSV collaborations.

3. Physics of TMD To access the 3-D parton distributions, we go beyond the 1-dimmensional description by including the parton internal transverse momentum kT relative to the nucleon momentum direction. As discussed above, there are 8 of them at leading twist. Several novel phenomena arise from the interactions involving these TMDs in transversely polarized nucleon scatterings.

In recent years, transverse spin physics has attracted tremendous attention from both experimental and theoretical communities, and as a result, significant progress has been made in understanding the spin dependent hadronic processes in QCD at high energy and the nucleon structure beyond the one-dimensional parton picture. Spin-dependent 3-D image of the parton distributions in the nucleon is under intense investigation and is the focus of the future EIC program, and also several novel transverse spin-dependent asymmetries have been confirmed experimentally in both SIDIS and p+p collisions, including Sivers and Collins transverse single spin asymmetries (TSSAs). One good example of the recent highlight is the study of the sign change in the Sivers transverse single spin asymmetry between SIDIS and Drell-Yan processes. Sivers TSSA were first measured in the Drell-Yan and W+/- production processes in hadronic interactions by the COMPASS experiment at CERN and the STAR experiment at RHIC, respectively. Although statistically limited, the data are consistent or in favor of the sign change scenario.

3.1 Physics of the transeverse single spin asymmetry in p+p collisions Contrary to the original expectation of very small TSSA, first experimental data from ZGS Do we understand this?

Kane, Pumplin, Repko, PRL, 1978 ! Early attempt: Do we understand this? showed surprisingly large asymmetry2 [5], orders of magnitudes larger than the naïve expectation AN ~ mq/pT. Such large TSSA have been persistent Kane,from thePumplin early, Repkofixed ,target PRL, 1978 energy all the way ps = 510GeV/c2 ! Early attempt: to the top collider energy at RHIC, AAAB/nicbVBNS8NAEJ34WetXVDx5WSyCp5pURS9C0YMeK9gPaGPZbDft0s0m7m6EEgr+FS8eFPHq7/Dmv3Hb5qCtDwYe780wM8+POVPacb6tufmFxaXl3Ep+dW19Y9Pe2q6pKJGEVknEI9nwsaKcCVrVTHPaiCXFoc9p3e9fjfz6I5WKReJOD2LqhbgrWMAI1kZq27st9SB1qobo4tR10DWtHZH7UtsuOEVnDDRL3IwUIEOlbX+1OhFJQio04VippuvE2kux1IxwOsy3EkVjTPq4S5uGChxS5aXj84fowCgdFETSlNBorP6eSHGo1CD0TWeIdU9NeyPxP6+Z6ODcS5mIE00FmSwKEo50hEZZoA6TlGg+MAQTycytiPSwxESbxPImBHf65VlSKxXd42Lp9qRQvsziyMEe7MMhuHAGZbiBClSBQArP8Apv1pP1Yr1bH5PWOSub2YE/sD5/ACCflEk= , see [1,6] . 2 AB(pT , ~s) Cross section: / At the leading+ twist,+... it was expected that the TSSA of high pT hadrons produced in the Cross section: transverselyAB(pT , ~s) polarized p+p collisions+ should+... be very small, as shown in Fig. 5. / " # AN =

AAACR3icbVBNS8MwAE3n15xfVY9egkMQxNGKoBdh6sWTTHAfsHYlzdIZlqYlSR2j9N958erNv+DFgyIezbaCuu1B4PE+SPL8mFGpLOvVKCwsLi2vFFdLa+sbm1vm9k5DRonApI4jFomWjyRhlJO6ooqRViwICn1Gmn7/euQ3H4mQNOL3ahgTN0Q9TgOKkdKSZ3YuvVt4AZ1AIJxCR9JeiDqpk8RIiGiQweNfrRsN+ETN0tng0dygZ5atijUGnCV2TsogR80zX3QdJyHhCjMkZdu2YuWmSCiKGclKTiJJjHAf9UhbU45CIt10vEMGD7TShUEk9OEKjtW/jRSFUg5DXydDpB7ktDcS53ntRAXnbkp5nCjC8eSiIGFQRXA0KuxSQbBiQ00QFlS/FeIHpCdVevqSHsGe/vIsaZxUbKti352Wq1f5HEWwB/bBIbDBGaiCG1ADdYDBE3gDH+DTeDbejS/jexItGHlnF/xDwfgBPKe0BQ== " + #

mq Asymmetry: AB(pT , ~s) AB(pT , ~s) = mq ↵s AZIMUTHALAsymmetry: ASYMMETRIES FOR HADRONAB(...pT , ~s) AB(pT , PHYSICAL~s) REVIEW D 83, 034021 (2011) = ↵s / pT which has been presented and discussed at length ina transverse momenta. Analogously, z and k  are the light- / pT series of papers (see, e.g., Refs. [39,42,43]). We will then cone momentum fraction and the transverse? momentum of (pred.) present the expression of the polarized cross section for theToo the observedsmall pion insideto theexplain jet with respect available to (w.r.t.) the data! A 0 AAAB/XicbVDLSgMxFM34rPU1PnZugkWom2FGBF1W3biSCvYB7ThkMpk2NMkMSUaoQ/FX3LhQxK3/4c6/MW1noa0HAodz7uHenDBlVGnX/bYWFpeWV1ZLa+X1jc2tbXtnt6mSTGLSwAlLZDtEijAqSENTzUg7lQTxkJFWOLga+60HIhVNxJ0epsTnqCdoTDHSRgrs/Yvg5j6vmkzkHI9gV1EO3cCuuI47AZwnXkEqoEA9sL+6UYIzToTGDCnV8dxU+zmSmmJGRuVupkiK8AD1SMdQgThRfj65fgSPjBLBOJHmCQ0n6u9EjrhSQx6aSY50X816Y/E/r5Pp+NzPqUgzTQSeLoozBnUCx1XAiEqCNRsagrCk5laI+0girE1hZVOCN/vledI8cTzX8W5PK7XLoo4SOACHoAo8cAZq4BrUQQNg8AiewSt4s56sF+vd+piOLlhFZg/8gfX5A4ZTk/Y= N process of interest, discussing in detail the different par- jet (parton c) direction of motion. ⇠ tonic contributions to the process; we will finally list the a=A;SA f^ x ; k contains all information on the   a=A;SA a a a a0 ð ? Þ azimuthal! asymmetries What that do can be Too measuredwe smallneed? and their to polarizationexplain state available of the initial parton data!a, which depends in physical content. In Sec. III we will present phenomeno- turnFigure on the 5. (experimentally Due to small fixed) quark parent mass, hadron mAqpolar-, very small TSSA was expected for the outgoing hadrons with large pT. logical results for the azimuthal asymmetries discussed in ization state and on the soft, nonperturbativeµ⌫↵ dynamics the kinematical configuration ofA theN RHIC experiments,i~sp at (p~encoded h inp~ theT eight) leading-twisti✏ polarized andphµ transverses⌫ p↵ph0 !different What c.m. energiesdo we and forneed? central-/ and forward-· momentum–dependent⇥ ) parton distribution functions, rapidity jet production. In particular, we will first present Need a phase,whichTo a willexplain spin be discussed theflip, in observed the enough following. largea=A;SA isvectors TSSA the s, several new ideas were proposed. Among them, here are results for the totally maximized effects, by taking all aa0 TMD functions saturated to natural positivity bounds andtwohelicity leading density matrixformalisms of parton a. Analogously,that can the generate po- large TSSA in SIDIS and p+p collisions: adding in sign all possible partonic contributions. This will larizationµ⌫↵ state of parton b inside the unpolarized hadron B AN i~sp (p~h p~T ) 1.i ✏ Siverspb=Bhµ function^s⌫ p↵p0 – a TMD function of the correlation of proton spin and parton transverse assess the! potential Vanish phenomenological without relevance of each parton’sis encoded into   transverse fb=B xb; k hb . motion: b b0 ð ? Þ effect. We will then consider/ more carefully· those⇥ effects ) The Mmomentum^ ’s are the pQCDk . leading-orderFor the helicitytransversely polarized proton, none-zero Sivers function implies c;d;a;b T involving the Sivers and Boer-Mulders distributions and scattering amplitudes for the hard partonic process ab cd. the Collins fragmentationNeed function, a phase, forA which direct phenomeno- a spin probe flip, enoughparton^  for parton’stransverse vectors momentum transverse! distribution motion, is not azimuthally symmetric, thus generating The D ; z; k  ’s are the soft leading-twist TMD logical parametrizations obtained by fitting combined data c c0 ð ? Þ fragmentationfinit functionse TSSA describing in the the fragmentation outgoing pro- hadron distribution [7]. for azimuthal asymmetries in SIDIS, Drell-Yan, and eþeÀ collisions are available. Section IV containsSpin-orbital our final re- cess2. ofcorrelation,Collins the scattered function (polarized) partonQCD – ac TMDinto quantum the finalfunction that interference correlates the transvers spin of a fragmenting quark marks and conclusions. leading pion inside the jet. ! Vanish without parton’s transverseAs already(transversity) said, we motion: will consider and as initial the particles transverseA, B, momentum distribution of the daughter hadrons [8]. two spin-1=2 hadrons (typically, two protons) with hadron II. FORMALISM AndB unpolarized their equivalents and hadron A in a pure in transverse the colinear spin state twist-3 approach, kT integrated quark-gluon-quark and tri- A In this section we present and summarize the expres-gluondenoted correlation by SA, with polarization functions (pseudo)vector [1,9P ].. sions of the polarized crossA direct section and probe of the measurable for parton’sEj and pj transverseare, respectively, the motion, energy and three- azimuthal asymmetries for the process A"B jet momentum of the observed jet.  X, where A and B are typically a pp or p!p pair.þ Unless otherwise stated, we will always work in the AB 15 Sinceþ most of the formalism has been already presentedNonehadronic-zero c.m. Sivers frame, with and hadron CollinsA moving alongfunctions the were first observed by the HERMES and the COMPASS in Refs. [39,42,43], we will shortly recall the main ingre- ^ Spin-orbital correlation,CollinsZcm direction; QCD we quantum willEffect define XZ cminterferenceas the in production Jets (Mid-Rapidity) dients of the approach, discussing more extensively onlyexperimentsplaneþ containing thein colliding polarized beamsð Þ and SIDIS the observed scattering jet, s in charged hadron productions, see recent reviews of relevant details specific to the process considered. with p > 0. We therefore have, neglecting all masses the historyð jÞXcm in [1]. Latest theoretical development show that these physics processes can be studies Within a generalized TMD parton model approach in- (see also Fig. 1): cluding spin and intrinsic parton motion effects, and as-in• p+pFirst coll measurementisions also, see of[1 0Collinss] and referenceseffect in � therein.+ � collisions suming factorization, the invariant differential cross section for the process A SA B jet  X can be written, at leading twist inð theÞ soft! TMDþ functions,þ as • � = 500 GeV follows:

A SA B jet  X • Multi-dimensional binning Ejd ð Þ ! þ þ dxadxb 2 3 2 2 d k a ? d pjdzd k  ¼ a;b;c;d;  16 xaxbs ? X f gZ �� − � d2k a=A;SA f^ x ;k b=B f^ x ;k b   a=A;SA a a   b=B b b  ? a a0 ð ? Þ b b0 ð ? Þ ^ ^ ^  M ; ; ; Mà ; ; ;  s^ t^ u^ D ; z;k  : (1)  c d a b c0 d a0 b0 ð þ þ Þ c c0 ð ? Þ In an LO pQCD approach the scattered parton c in the hard elementary process ab cd is identified with the observed fragmentation jet.! Let us summarize briefly the physical meaning of the terms in Eq. (1). Full details and technical aspects can be found in Refs. [39,42,43]. We sum over all allowed partonic processes contributing to the physical process observed.  stays for a sum over f g all partonic helicities,  1=2 1 for quark (gluon) FIG. 1 (color online). Kinematical configuration for the pro- partons, respectively. x ¼Æand k ðÆ are,Þ respectively, the a;b a;b Figurecess A6S. ALH:B Collinsjet  asymmetryX in the hadronic can c.m.be studied reference with hadron distribution inside a jet. RH: Collins asymmetry measured by STAR in initial parton light-cone momentum? fractions and intrinsic frame. ð Þ ! þ þ transverely polarized p+p collisions at center of mass envery 500GeV.Phys.Rev . D97, 032004 (2018)

034021-3 Comparison with 3.2 Sivers function and sign change Phys. Lett. B773, 300-306 (2017) Among the 8 TMD functions, the Sivers functionarXiv:1707.00913 has attracted particular interest. Sivers function describes the correlation of intrinsic parton transverse motion and the nucleon’s spin direction. This function can be accessed either in the SIDIS process or Drell-Yan in p+p collisions. Besides providing the 3-D parton structure of the nucleon, the most interesting property of Sivers function is its process dependence, and as a consequence of the QCD gauge invariance and factorization assumption, the sign of Sivers function in SIDID and Drell-Yan (and W+/-) is expected to be opposite [11]. This nontrivial process-dependence of Sivers function is a fundamental prediction of QCD. With the first confirmation of non-zero Sivers functions in SIDIS from COMPASSSivers(func4on(is(process(dependent( and HERMES, as well as several JLab experiments, it is extremely important to test this predictionSivers(func4on(is(process(dependent(. Recently, first experimental studies were carried out by the COMPASS experiment • In%• QatIn%C CERND:%QCD:% for Drell-Yan and by the STAR experiment at RHIC for the weak boson W+/- productions 2 – atSign%change%of%the%quark%– much Sign%change%of%the%quark% larger Q , seeSivers FigSivers%func:on%between%DIS%s. %func:on%between%DIS%7 and 8. and%and%Drell DrellRYan%process.%RYan%process.%

DIS:DIS: attractive attractive Drell-Yan:Drell-Yan: repulsive repulsive

As a result: As a result:

• Urgent%need%to%achieve%these%milestone%measurements!%• Urgent%need%to%achieve%these%milestone%measurements!% • Figure 7. It is predicted Sivers function changes sign between SIDIS and Drell-Yan processes. RH: First Drell-Yan TSSA Determine%the%• Determine%the%result fromSiversSivers COMPASS,%func:on%from%global%QCD%analysis.%%func:on%from%global%QCD%analysis.% although statistically limited, it is consistent with the prediction. QCD%Evolu:on%2018% Xuan%Li%(LANL)% 5% QCD%Evolu:on%2018% Xuan%Li%(LANL)% 5%

Figure 8. TSSA of W+/- production in p+p collisions measured by the STAR experiment at RHIC. 3.3 Future prospects In coming years, RHIC will continue pursuing spin physics program with the current STAR and the future sPHENIX experiments, with possible forward detector upgrade for both detectors [3]. Also, there are several new polarized proton beam or target experiments either under commissioning now or being proposed. The upcoming Fermilab SpinQuest/E1039 fixed target experiment will utilize a polarized proton target (NH3) to study Drell-Yan TSSA using the 120GeV proton beam from the Main Injector [12]. The kinematics of the dimuon spectrometer selects sea- quarks from the thus directly accessing sea quarks’ Sivers function in the x range from 0.1 to 0.3, see Fig. 9.

February16,2015 14:58 WSPC/CRC 9.75x6.5 1560064

A. Klein

February16,2015 14:58 WSPC/CRC 9.75x6.5 1560064

Sea￿est SpectrometerA. Klein

Targets: LH2, LD2, C, Fe, W 4 y tion arra Sta pe osco S od es 3 H tub Hadron Absorbers (stop beam, , 2 & rop s 1 P on rray tati e a S scop r odo be H ham muon identification) ift c on t Dr eV ir ne T lid ag K t Fig. 4. ProjectedSo g m statisticaler precisiongne and comparison with two calculations usin orb Ma Foc abs dron er Ha ump Magnets (focussing, momentum m d Bea determination) 4. Summary r rbe bso n a ) dro all Ha on w 4 tracking stations, consist of eam (Ir B Experiment E1039 will be the first measurement25 m of proton induced Drell-Yan with Hodoscopes apolarizedtarget.Throughthechosenkinematicswewillbeabletomeasurethe& D 2 , L ts LH 2 rge d ta Dri￿ Chambers (St. 1-3) or Sivers Single Spins Asymmetryoli for the sea quarks. Our measurement will allow us Prop. Tubes (St. 4) to determine the magnitude and sign of this asymmetry for 4 different bins in Fig. 4. Projected statistical precision and comparison with two calculations xtarget withFigure excellent 9. SpinQuest/E0139 statistical precision. experiment Because at Fermila theMass Siversb distributionwill asymmetrystudy seaquark fi￿ed is with related Sivers estimated function. LH: SpinQuest to the quarkspectrometer; orbital angular RH: Projected momentum, Drell-Yan such TSSA a measurement from 2-year 4.operation will Summary help of solve E0139 the at Fermilab. componentsExperiment E1039 will be the first measurement of proton induced Drell-Yan with nucleon spin puzzle. A nonzero asymmetry requires an orbitalapolarizedtarget.Throughthechosenkinematicswewillbeabletomeasurethe angular momentum Sivers Single Spin Asymmetry for the sea quarks. Our measurement will allow us Int. J. Mod. Phys. Conf. Ser. 2015.37. Downloaded from www.worldscientific.com At CERN, there is a new COMPASS++/AMBERWell fi￿ed: proposal to continue hadron physics program contribution from the sea quarks, and would be a “smoking gun”to determine for Lu the¯ =0,while magnitude and sign of this asymmetry for 4 different bins in using the M2 beamline at SPS. This new QCDDetectors facility &xtarget trackingwith will excellent tool̸ statisticalcover work precision. as a expected broad Because the range Sivers asymmetry of isphysics related topics azeroasymmetrywouldleadtoLu¯ =0,whichwouldmakethe“spinpuzzle”evento the quark orbital angular momentum, such a measurement will help solve the nucleon spin puzzle. A nonzero asymmetry requires an orbital angular momentum Int. J. Mod. Phys. Conf. Ser. 2015.37. Downloaded from www.worldscientific.com from nucleon structures, hadron spectroscopyDrell–Yan to dark cancontribution be matter selected from the sea withsearch quarks, and would [1 be3 a]. “smoking gun” for Lu¯ =0,while more puzzling. ̸ azeroasymmetrywouldleadto2 Lu¯ =0,whichwouldmakethe“spinpuzzle”even mass > 4.2 moreGeV puzzling./c by 60.217.68.152 on 07/25/18. Re-use and distribution is strictly not permitted, except for Open Access articles. by 60.217.68.152 on 07/25/18. Re-use and distribution is strictly not permitted, except for Open Access articles. References References Kei Nagai (Academia Sinica) Light-antiquark Flavor Asymmetry1. See forQNP example: 2018 HERMES @ Tsukuba Collab.14th (A. Nov., Airapetian 2018et 9 al. / 17) Phys. Rev. Lett. 94,012002 (2005); Phys. Rev. Lett. 103 (2009) 152002; COMPASS Collab. (V. Yu. Alexakhin et 1. See for example: HERMES Collab. (A. Airapetian et al.) Phys. Rev.al.) Phys. Lett. Rev.94 Lett.,01200294,202002(2005);JlabHallACollab(X.Qianet al.) Phys. Rev. Lett. 107,072003(2011). (2005); Phys. Rev. Lett. 103 (2009) 152002; COMPASS Collab. (V.2. See Yu. for example: Alexakhin S. D. Bass, Rev.et Mod. Phys. 77,1257(2005). 3. M. Deka et al.,arXiv:1312.4816(2013). al. ) Phys. Rev. Lett. 94,202002(2005);JlabHallACollab(X.Qian4. NMC Collabet (P. al. Amaudruz) Phys.et al.) Phys. Rev. Lett. 66,2712(1991). Rev. Lett. 107,072003(2011). 1560064-6 2. See for example: S. D. Bass, Rev. Mod. Phys. 77,1257(2005). 3. M. Deka et al.,arXiv:1312.4816(2013). 4. NMC Collab (P. Amaudruz et al.) Phys. Rev. Lett. 66,2712(1991).

1560064-6

Figure 10. LH: COMPASS++/AMBER, a proposed New QCD Facility at M2 beam line of the SPS at CERN; RH: Proposed Spin Phys ics Detector at NICA to study polarized nucleon structures.

At Dubna, Nuclotron based Ion Collider fAcility (NICA) is a flagship project of the Joint Institute for Nuclear Research which is expected to be in operation by 2021. Main tasks of “NICA Facility” are study of hot and dense baryonic matter, investigation the polarization phenomena and the nucleon spin structure. The Spin Physics Detector (SPD) will carry out nucleon spin structure studies with Drell-Yan, J/Psi and other observables [14].

4. Conclusions

Tremendous progress has been made in the last decades in studying the 3-D structures of the nucleon and their roles in strong interactions. Recent polarized proton-proton or hadron collision data have shed new light in our understanding of phenomena that have puzzled us for decades, such as the “”, the unexpected large TSSA in hadron production in transversely polarized p+p and DIS interactions, and the first experimental hints of process dependent sign change in Sivers function etc. Due to space limitation, only a few selected recent results are presented, and there are many interesting results that I just even couldn’t touch upon at all. I encourage interested readers to check out the review articles listed below and references therein.

Acknowledgment

The author thanks many experimental colleagues from RHIC, Fermilab, HERMES and COMPASS experiments and theory friends for valuable discussions. This work was supported in part by the Office of Nuclear Physics within the U.S. DOE Office of Science and Los Alamos National Laboratory LDRD program office.

References [1] J. Ralston and D. Soper, Nucl. Phys. B152, 109 (1979); R. L. Jeffe and X. Ji, Phys. Rev. Lett. 67, 552 (1991); M. Amselmino et al, arXiv:0907.3999 J.C. Peng and J.-W. Qiu, Prog. Part. Nucl. Phys. 76, 43 (2014). M.G. Perdekamp and F. Yuan, Annu. Rev. Nucl. Part. Sci. 65, 429 (2015). E. C. Aschenauer et al, RHIC-Spin Collaboration, arXiv:1501.01220 [2] EMC Collaboration, Phys. Lett. B 206 (2) 364 (1888) EMC Collaboration, Nucl. Phys. B 328, 1 (1989). R. L. Jaffe, Physics Today 48, 9, 24 (1995) T. P. Cheng and L-F Li, Phys. Rev. Lett. 62, 1441 (1989) S. Brodsky, J. Ellis, and M. Karliner, Phys. Lett. B 206, 309 (1988). [3] D. de Florian, R. Sassot, M. Stratmann and W. Vogelsang, Phys. Rev. Lett. 113, 012001 (2014); NNPDFPol Collaboration, Nucl. Phys. B 887, 276 (2014); STAR Collaboration, Phys. Rev. Lett. 115, 092002 (2015); PHENIX Collaboration, Phys. Rev. D 90, 012007 (2014); J. Blumlein and H. Bottcber, Nucl. Phys. B 841, 205 (2010); E. Leader, A. V. Sidorov and D. B. Stamenov, Phys. Rev. D 82, 114018 (2010); M. Hirai and S. Kumano, Nucl. Phys. B813, 106 (2009); C. Bourrely, F. Buccella and J. Soffer, Phys. Rev. D 83, 074008 (2011); P. Jimenez-Delgado, A. Accardi and W. Melnitchouk, Phys. Rev. D 89, 034025 (2014) E. C. Aschenauer et al, RHIC-spin collaboration, arXiv:1602.03922 [4] C. Bourrely and J. Soffer, Phys. Lett. B 314, 132 (1993); P. M. Nadolsky and C. P. Yuan, Nucl. Phys. B 666, 31 (2003); D. de Florian and W. Vogelsang, Phys. Rev. D 81, 094020 (2010); STAR Collaboration, Phys. Rev. Lett. 113, 072301 (2014); PHENIX Collaboration, Phys. Rev. D 93, 051103(R) (2016); PHENIX Collaboration, Phys. Rev. D 98, 032007 (2018). [5] G. L. Kane, J. Pumplin and W. Repko, Phys. Rev. Lett. 41, 1689 (1978); R.D. Klem et al, ZGS, Phys. Rev. Lett. 36, 929 (1976); E704, PLB261, 201 (1991), PLB264, 462 (1991); C. E. Allgower et al. (E925 Collaboration), PRD65, 092008 (2002); [6] STAR, Phy. Rev. Lett. 92, 171801 (2004); and more. [7] D. W. Sivers, Phys. Rev. D 41 83 (1990) [8] J. C. Collins, Nucl. Phys. B 396, 161 (1993) [9] A. V. Efremov and O. V. Teryaev, Sov. J. Nucl. Phys. 36, 140 (1982), Phys. Lett. B 150 383 (1985); J.W. Qiu and G. F. Sterman, Phys. Rev. Lett. 67, 2264 (1991), Nucl. Phys. B 378, 52 (1992) [10] F. Yuan, Phys. Rev. Lett. 100, 032003 (2008); T. C. Roger and P. J. Mulders, Phys. Rev. D 81, 094006 (2010); Umberto D’Alesio et al PRD 83 034021 (2011); A. Metz and J. Zhou, Phys. Lett. B 700, 11 (2011) L. Gamberg, B. B. Kang and A. Prokudin, Phys. Rev. Lett. 110, 232301 (2013); [11] J. C. Collins, Phys. Lett. B 536, 43 (2002); Z. B. Kang and J. W. Qiu, Phys. Rev. Lett.103, 172001 (2009) [12] SpinQuest/E1039 at Fermilab, A. Chen et al, arXiv:1901.09994, [13] COMPASS++/AMBER, arXiv:1808.00848 [14] NICA/SPD, arXiv: 1408.3959.