Advances in High Energy Physics
Volume 2019 (2019), Article ID 4825790, 15 pages
https://doi.org/10.1155/2019/4825790

Review Article

Introduction to the Transverse-Momentum-Weighted Technique in the Twist-3 Collinear Factorization Approach

Hongxi Xing and Shinsuke Yoshida

Institute of Quantum Matter and School of Physics and Telecommunication Engineering, South China Normal University, Guangzhou 510006, China

Correspondence should be addressed to Hongxi Xing; hxing@m.scnu.edu.cn and Shinsuke Yoshida; shinyoshida85@gmail.com

Received 21 December 2018; Accepted 31 March 2019; Published 2 June 2019

Guest Editor: Alexei Prokudin

Copyright © 2019 Hongxi Xing and Shinsuke Yoshida. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. The publication of this article was funded by SCOAP³.

Abstract.

The twist-3 collinear factorization framework has drawn much attention in recent decades as a successful approach in describing the data for single spin asymmetries (SSAs). Many SSAs data have been experimentally accumulated in a variety of energies since the first measurement was done in the late 1970s and it is expected that the future experiments like Electron-Ion-Collider will provide us with more data. In order to perform a consistent and precise description of the data taken in different kinematic regimes, the scale evolution of the collinear twist-3 functions and the perturbative higher-order hard part coefficients are mandatory. In this paper, we introduce the techniques for next-to-leading order (NLO) calculation of transverse-momentum-weighted SSAs, which can be served as a useful tool to derive the QCD evolution equation for twist-3 functions and to verify the QCD collinear factorization for twist-3 observables at NLO, as well as obtain the finite NLO hard part coefficients.

1. Introduction

The large single transverse spin asymmetries (SSAs) have been a longstanding problem over 40 years since it was turned out that the conventional perturbative calculation based on the parton model picture failed to describe the large SSAs which were experimentally observed in pion and polarized hyperon productions [1, 2]. In the recent several years, two QCD factorization frameworks have been proposed to study phenomenologically the observed SSAs: the transverse momentum dependent factorization approach [3–15] and the twist-3 collinear factorization approach [16–38]. These two frameworks are shown to be equivalent in the common applied kinematic region [39–46].

The twist-3 collinear factorization framework is a natural extension of the conventional perturbative QCD framework and it could give a reasonable description of the large SSAs. Measurements of SSAs at Relativistic-Heavy-Ion-Collider (RHIC) [47–51] have greatly motivated the theoretical work on developing the twist-3 framework, because it is a unique applicable framework for single hadron productions in proton-proton collision. A series of important works have been done in the past a few decades and the SSAs for the hadron production were completed at leading order (LO) with respect to the QCD strong coupling constant [16–38]. Recent numerical simulations based on the complete LO result confirmed that the twist-3 approach gives a reasonable description of the SSA data provided by RHIC [52, 53].

Electron-Ion-Collider (EIC) is a next-generation hadron collider expected to provide more data in different kinematic regimes for SSAs. In order to extract the fundamental structure of the nucleon from the measurements at a future EIC, comprehensive and precise calculations for SSAs in transversely polarized lepton-proton collision are highly demanded. It is well known that nonperturbative functions in the perturbative QCD calculation, in general, receive logarithmic radiative corrections and the evolution equation with respect to this logarithmic scale is necessary for a systematic treatment of the cross sections in wide range of energies. Most famous example is the DGLAP evolution equation of the twist-2 parton distribution functions (PDFs). Correct description of the small scale violation of the structure function controlled by the DGLAP equation was an important success of the QCD phenomenology in the early days. The twist-3 function is expected to have similar logarithmic dependence and its evolution equation will play an important role in global fitting of the SSA data accumulated in different energies. Consistent description of the data will be a good evidence that the twist-3 framework, one of major fundamental developments in recent QCD phenomenology, is a feasible theory to solve the 40-year mystery in high energy physics. The evolution equations for the twist-3 functions have been derived in two different methods. The first method is a calculation of the higher-order corrections to the nonperturbative function itself [54–61]. Since the nonperturbative function has the operator definition, we can investigate the infrared singularity of the operator through higher-order perturbative calculation. We can read the evolution equation from the infrared structure of the function. This is a standard technique and the evolution equations have been derived for the twist-3 distribution functions for initial state proton [54–59] and the twist-3 fragmentation functions for final state hadron [60, 61]. The second method which we will review in this paper is a transverse-momentum-weighted technique for the SSAs [62–67]. Except for deriving the QCD evolution equation for twist-3 nonperturbative functions, the transverse-momentum-weighted technique can be also used as a tool to verify the twist-3 collinear factorization at higher orders in strong coupling constant . There is also phenomenological interest related to this technique. One can use the standard dimensional regularization method to derive the NLO hard part coefficient for transverse-momentum-weighted SSAs, which can be used for high precision extraction of twist-3 functions from the relevant experimental data. The recent measurement of the transverse-momentum-weighted SSAs at COMPASS [68] strongly motivates the phenomenological application of the results reviewed in this paper. We expect more data will be produced in future COMPASS and EIC measurements.

The rest of the paper is organized as follows. In Section 2 we present the notation and the calculation of transverse-momentum-weighted SSAs at leading order for semi-inclusive deep inelastic scattering (SIDIS). In Section 3 we present the detail of NLO calculation for both real and virtual corrections, we show the cancelation of soft divergence in the sum of real and virtual corrections, and the collinear divergences can be absorbed into the redefinition of twist-3 Qiu-Sterman function and unpolarized leading twist fragmentation function. In Section 4 we review the application of the transverse-momentum-weighted technique to other processes that have been done in recent years. We conclude our paper in Section 5.

2. Transverse-Momentum-Weighted SSA at Leading Order

In this paper, we take the process of SIDIS as an example to show the techniques of perturbative calculation for transverse-momentum-weighted differential cross section at twist-3. We start this section by specifying our notation and the kinematics of SIDIS and present the calculation for transverse-momentum-weighted SSA at leading order (LO).

2.1. Notation

We consider the scattering of an unpolarized lepton with momentum on a transversely polarized proton with momentum and transverse spin and observe the final state hadron production with momentum ,

We focus on one-photon exchange process with the momentum of the virtual photon given by and its invariant mass . We define all vectors in the so-called hadron frame. We define to be the momentum for the parton that fragments into the final state hadron. The conventional Lorentz invariant variables in SIDIS are defined as

For clear understanding, we start with the -integrated cross section at leading twist in unpolarized lepton-proton scattering

There is only one hard scale in this case; therefore the differential cross section shown above can be reliably computed by using the standard collinear factorization formalism. The LO contribution is given by scattering amplitude . It is trivial that the LO cross section is proportional to the unpolarized PDFs and the unpolarized fragmentation functions (FFs),

where is the LO Born cross section with being the QED coupling constant. The bare results at contain infrared divergences which represent the long-range interaction in hadronic collision process. These divergences are canceled by the renormalization of the PDFs and FFs and the DGLAP evolution equations are derived as the renormalization group equations. The final result at NLO can be written as the convolution of finite hard part coefficient and nonperturbative functions (PDFs and FFs) [63]

where represents convolution.

The concept of the transverse-momentum-weighted technique is mostly the same with the twist-2 case. Notice that direct -integration of the cross section for unpolarized lepton scattering off transversely polarized proton vanishes due to the linear dependence of . Realize that the SSA is characterized in terms of three vectors: the momentum of the final state hadron, the momentum, and the spin of the initial state proton, which can be combined as

where is a two-dimensional antisymmetric tensor with ; is an arbitrary vector which satisfies and . We introduce a weight factor and consider the following transverse-momentum-weighted differential cross section

which is well defined after -integration. Since the virtuality is the only hard scale after the -integration, one can safely use the collinear twist-3 factorization formalism, and the technique in performing NLO calculation will follow those used at leading twist. Same technique has been applied to Drell-Yan dilepton production in proton-proton collisions [62, 66] and can be extended to polarized electron-positron collisions.

We recall the cross section for SIDIS presented in [64],

where is the leptonic tensor. We focus on the metric contribution and the SSA generated by initial state twist-3 distribution functions of the transversely polarized proton. Then we can factorize the nonperturbative part by introducing the usual twist-2 unpolarized fragmentation function

The hadronic tensor describes a scattering of the virtual photon and the transversely polarized proton. We will make the subscript implicit in the rest part of this paper for simplicity.

2.2. Leading Order

We demonstrate how to derive the LO cross section for the transverse-momentum-weighted SSA based on the collinear twist-3 framework and show that the LO cross section is proportional to the first moment of the TMD Sivers function. The twist-3 calculation is well formulated in the diagrammatic method. We consider a set of the general diagrams shown in Figure 1 and extract twist-3 contributions from these diagrams. We start from the first diagram in Figure 1, which can be expressed as

where and are the momenta of the parton from initial state proton and that fragments to the final state observed hadron, respectively, . The hard part at LO is given by

We perform -integration before the collinear expansion,

Because the -integration gives factor, we can identify the leading term in the collinear expansion as a twist-3 contribution. We perform the collinear expansion of the hard part around

and substitute the above expansion into the hadronic tensor as shown in (10)

where . Now we turn to the second and the third diagrams in Figure 1. These two diagrams can be expressed as

where the hard parts are given by

The -integration gives and then the leading term in collinear expansion gives twist-3 contribution again,

For the matrix element, we have to separate the components of the gluon field into “longitudinal" and “transverse" part as

The longitudinal part gives the leading contribution. It is straightforward to derive the Ward-Takahashi identities (WTIs) for the hard parts,

Finally the hadronic tensor shown in (15) can be expressed as

Combining (14) and (20), we can obtain the result

We can find that this matrix element corresponds to term in the first moment of the TMD correlator

where represents the Wilson line, is the nucleon mass, and we used the fact that the first moment of the TMD Sivers function gives the Qiu-Sterman function . The nonlinear term in the field strength tensor and the higher-order terms in the Wilson lines which have to be added to (21) come from the more gluon-linked diagrams in Figure 2. Finally we can derive the LO cross section formula as

where . As we demonstrated here, the LO cross section for the weighted SSA is proportional to the first moment of the TMD function (also known as the kinematical twist-3 function). We can expect the NLO contribution gives the evolution equation to the twist-3 Qiu-Sterman function, in this case. This technique is quite general so that we can apply the same technique to other TMD functions. When we focus on the twist-3 fragmentation effect, we can derive the evolution equation for the first Moment of the Collins function. When we consider the double spin asymmetry and change the weight factor , we can investigate other TMD distribution functions like the Worm-Gear and the Pretzelosity.

Figure 1: A series of LO diagrams in the diagrammatic method.

Figure 2: The NLO virtual correction diagrams in SIDIS.

3. Transverse-Momentum-Weighted SSA at NLO

In this section, we review the calculation for transverse-momentum-weighted SSA at NLO including both real and virtual corrections.

3.1. Virtual Correction

We first consider the NLO contribution from the virtual correction which is given by the scattering amplitude with one gluon loop. When we adopt the dimensional regularization scheme, the gauge-invariance of the cross section is maintained for the loop diagram. Then we can derive the same WTI as shown in (19) which is the consequence of gauge-invariance of the diagrams. Taking advantage of this fact, we just need to calculate simple diagrams shown in Figure 2 for the virtual correction. The calculation for this kind of one-loop diagrams has been well established, which is exactly the same as the vertex correction at leading twist. We follow the conventional technique here. All ultraviolet divergences can be canceled by the renormalization of the QCD Lagrangian. Then we can set the ultraviolet and infrared divergences are the same with each other in dimensions regularization approach, , and identify all divergences as infrared. In this definition, we do not have to think about the first and the third amplitudes in Figure 2 because these are exactly zero in the mass case as we considered here. The hard partonic cross section with the second amplitude is given by

where in -dimension, we made a change for -dimensional calculation, and we used the fact that

We perform the basic -dimensional calculation for each integration,

The complex-conjugate diagram gives the same contribution. Then we can show the cross section for the NLO virtual correction as

this is exactly the same as the virtual correction at leading twist. The strategy in the virtual correction calculation presented here is different with that shown in [63], in which the authors did not use the WTI shown in (19) and directly calculated which has one more external gluon line with the momentum in Figure 2. The authors obtained a consistent result with (28), which demonstrated the validation of WTI through explicit calculation by including all virtual diagrams shown in Figures 3 and 4 in [63]. We would like to comment that the WTI reduces much calculational cost. The direct calculation of takes tremendous time as they contain significant amount of tensor reduction and integration. These two calculations should be conceptually the same with each other as long as we correctly keep track of all imaginary contributions. We confirmed in this section the consistency mathematically in -scattering case. The consistency check in a more general way will be a future task in the collinear twist-3 factorization approach.

3.2. Real Correction

We now complete the NLO calculation by adding the real emission contribution represented by partonic scattering process. The calculation for scattering diagrams has been well studied in -unintegrated case. We just have to repeat the same calculation but in -dimension. We adopt the conventional technique by separating the propagator into the principle value part and imaginary part [18–33],

and we focus on the pole contribution which is required to generate the phase space for SSA. The derivation of the cross section for the pole contribution has been well developed so far based on the diagrammatic method we reviewed in Section 2. Here we recall the result derived in [23–28] as

where the matrix element is given by

We can construct the gauge-invariant expression (30) before performing the -integration. There are three types of pole contributions in SIDIS which are, respectively, known as soft-gluon pole contribution (SGP, ), hard pole contribution (HP, or ), and another hard pole contribution (HP2, or ) [23–28, 64], and the corresponding hard parts are given by

Typical diagrams for each pole contribution are shown in Figure 3 (full diagrams can be found in [64]). We write down the explicit form of each diagram in Figure 3 in order to help readers to follow our calculation,

where is the sum of the polarization vector . We can show the WTI for the pole diagrams

and it gives a useful relation

For HP and HP2 contributions, we can use (35) and do not have to perform -derivative directly as in (30). However, we cannot use this relation for SGP contribution because it contains a delta function . In [69], the authors found a reduction formula for SGP contribution as

where is the -scattering cross section without the external gluon line with momentum . Substituting (31) into (30) and using ((35), (36)), we can derive the following result:

where we used the Mandelstam variables

with , . Then the cross section can be written as

where we used the symmetry of the -integration as

The authors of [63] found that the factor in the denominator is essential to derive the correct evolution function of . They calculated all SGP and HP contributions associated with . After that, the HP2 contribution and all contributions were calculated in [64]. The results of all hard cross sections are listed in [64]. We perform the -integration,

where is the solid angle and it can be integrated out as

We carry out the -expansion as follows:

Then the cross section can be derived as follows:

where the derivative term was converted to the nonderivative term through the partial integral,

Combine the real (see (28)) and virtual corrections (see (49)), we find that the double-pole term, which represents a soft-collinear divergence, cancels out. The remaining collinear divergences, represented by the single pole , can be eliminated by the redefinition of Qiu-Sterman function and fragmentation function,

where is well known splitting function

and we adopted the -scheme

This collinear singularity is consistent with that of explored in [54–59]. After the renormalization, the final result for the transverse-momentum-weighted SSAs at NLO is given by

Note that the cross section above does not include the contribution from the gluon fragmentation channel. From the requirement that the physical cross section does not depend on the factorization scale ,

we can derive the LO evolution equation for ,

This evolution equation based on the transverse-momentum-weighted technique was first discussed in Drell-Yan process [62], in which the authors succeeded in deriving the terms in the parenthesis in the evolution equation. After that, the authors of [63] pointed out that an extra term also contributes to the evolution equation. The HP2 pole contribution and all terms were obtained in [64] within the method of transverse momentum weighting.

Figure 3: Typical diagrams for soft-gluon pole (left), hard pole (middle), and another hard pole (right). The red barred propagator gives the pole term.

4. Application to Other Processes

A lot of works on the transverse-momentum-weighted SSA have been done in recent years. We briefly summarize all related work in this section. The evolution equation of the Qiu-Sterman function (57) is still missing the gluon mixing contribution associated with the 3-gluon distribution functions defined by

where and are the structure constants of group. Figure 4 shows typical diagrams which give the gluon mixing contribution in SIDIS and Drell-Yan.

Figure 4: The diagrams which give the gluon mixing contribution to the evolution equation of in SIDIS (left) and Drell-Yan (right).

The calculation technique for these diagrams was developed in [29–33]. There is only the soft-gluon pole () contribution due to the interchange symmetry of the external gluon lines and then the cross section is expressed by four independent functions , , , and . The gluon mixing term in the flavor singlet evolution was discussed in both SIDIS [65, 67] and Drell-Yan [66] and the cross section up to the finite term is derived as

where , with the center of mass energy and the rapidity , is given by a linear combination of the 3-gluon distribution functions, and is the antiquark PDF. In Drell-Yan process, we use the transverse momentum of the virtual photon for the weighted cross section. Using the condition as in (56), we can derive the gluon mixing term of . The explicit form of the evolution kernel and the finite terms are shown in [66, 67]. Adding the mixing term to (57), the evolution equation for the Qiu-Sterman function, the first moment of the TMD Sivers function, was completed at LO with respect to QCD coupling constant . As a by-product of the work on Drell-Yan, the NLO cross section related to the first momentum of the TMD Boer-Mulders function was also derived as

where is chiral-odd transversity distribution and the twist-3 function corresponds to the first moment of the Boer-Mulders function and its exact definition can be found in [66]. The evolution kernel and the finite terms of the cross section are shown in [66]. The authors in [67] also discussed the evolution equation for the first moment of the TMD Collins fragmentation function. The calculation of the twist-3 fragmentation contribution is different from the distribution case because the cross section only receives the nonpole contribution of the hard scattering. The calculation technique for the nonpole contribution has been well developed in [26, 34–36]. It is not difficult to extend it to the -dimensional calculation. The NLO cross section was derived in [67] as

where , , and are the twist-3 fragmentation functions for spin- hadron (different definition is shown in [70]). All evolution kernels and finite terms are shown in [67].

The NLO transverse-momentum-weighted cross section has been completed for single inclusive hadron production in SIDIS and Drell-Yan process in proton-proton collisions by a series of work presented here. These results are useful not only for the derivation of the evolution equations, but also for the verification of twist-3 collinear factorization feasibility, as well as for the global analysis of the experimental data. Measurement of the weighted SSAs just began very recently [68] and more data will be provided by future experiments. The analysis of the data based on the NLO result will lead to better understanding of the origin of the SSAs. There are still some TMD functions which have not been discussed yet; we hope the techniques presented in this paper can help extending the application of the transverse-momentum-weighted technique to resolve all these open questions.

The transverse-momentum-weighted technique has also been extended to study the transverse momentum broadening effect for semi-inclusive hadron production in lepton-nucleus scattering [71] and Drell-Yan dilepton production in proton-nucleus scattering [72]. In these studies, the QCD evolution equation for twist-4 quark-gluon correlation function was derived for the first time, and the twist-4 (double scattering) collinear factorization at NLO was confirmed through explicit calculations [73]. The finite NLO hard parts were obtained for the transverse momentum broadening effect, which can be used in global analysis of world data to extract precisely the medium properties characterized by the twist-4 matrix elements.

5. Summary

We reviewed the transverse-momentum-weighted technique as a useful tool to derive the scale evolution equation for the twist-3 collinear function which is expressed by the first moment of the TMD function. We first demonstrated the calculation of the LO cross section formula in a pedagogical way. Then we showed the basic techniques for the NLO calculation for both the virtual correction and real emission contributions. A lot of work have been done on the Qiu-Sterman function [62–65] and recently the application of this technique to other twist-3 functions was also discussed [66, 67]. There is still room of the application to many other TMD functions by considering appropriate twist-3 observables. We hope our review paper will provide basic knowledge needed to work on this subject.

In the end, we would like to point out the importance of the -weighted SSA from the phenomenological point of view. We introduced this observable as a tool to derive the scale evolution equation by focusing on the -term in the cross section. However, the finite term is also important when we evaluate the cross section in order to compare it with the experimental data. The COMPASS experiment reported the data of the -weighted SSA very recently [68]. We expect that the data will be accumulated in the future experiments at COMPASS, JLab, and EIC and then the exact NLO cross section including the finite contribution will play an important role in the analysis of those weighted SSA data.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This research is supported by NSFC of China under Project no. 11435004 and research startup funding at SCNU.

References

R. D. Klem, J. E. Bowers, H. W. Courant et al., “Measurement of asymmetries of inclusive pion production in proton-proton interactions at 6 and 11.8 GeV/c,” Physical Review Letters, vol. 36, no. 16, pp. 929–931, 1976.

G. Bunce, “Λ⁰ hyperon polarization in inclusive production by 300-GeV protons on beryllium,” Physical Review Letters, vol. 36, p. 1113, 1976.

X. D. Ji, J. P. Ma, and F. Yuan, “QCD factorization for semi-inclusive deep-inelastic scattering at low transverse momentum,” Physical Review D, vol. 71, no. 12, 2005.

X. D. Jia, J. P. Maab, and F. Yuana, “QCD factorization for spin-dependent cross sections in DIS and drell-yan processes at low transverse momentum,” Physics Letters B, vol. 597, no. 3-4, pp. 299–308, 2004.

S. J. Brodsky, D. S. Hwang, and Schmidt I, “Final-state interactions and single-spin asymmetries in semi-inclusive deep inelastic scattering,” Physics Letters B, vol. 530, pp. 99–107, 2002.

S. J. Brodsky, D. S. Hwang, and I. Schmidt, “Initial-state interactions and single-spin asymmetries in Drell-Yan processes,” Nuclear Physics B, vol. 642, no. 1-2, pp. 344–356, 2002.

P. Mulders and R. Tangerman, “The complete tree-level result up to order 1/Q for polarized deep-inelastic leptoproduction,” Nuclear Physics B, vol. 461, no. 1-2, pp. 197–237, 1996.

A. Bacchetta, M. Diehl, K. Goeke, A. Metz, P. J. Mulders, and M. Schlegel, “Semi-inclusive deep inelastic scattering at small transverse momentum,” Journal of High Energy Physics, vol. 2007, article no 093, 2007.

D. Boer, P. J. Mulders, and F. Pijlman, “Universality of T-odd effects in single spin and azimuthal asymmetries,” Nuclear Physics B, vol. 667, pp. 201–241, 2003.

D. Boer and P. J. Mulders, “Time-reversal odd distribution functions in leptoproduction,” Physical Review D: Particles, Fields, Gravitation and Cosmology, vol. 57, pp. 5780–5786, 1998.

M. Anselmino, M. Boglione, U. D'Alesio, S. Melis, F. Murgia, and A. Prokudin, “Sivers effect in Drell-Yan processes,” Physical Review D: Particles, Fields, Gravitation and Cosmology, vol. 79, Article ID 054010, 2009.

M. Anselmino, M. Boglione, U. D'Alesio et al., “Sivers effect for pion and kaon production in semi-inclusive deep inelastic scattering,” The European Physical Journal A, vol. 39, no. 1, pp. 89–100, 2009.

Z. B. .Kang and J. W. Qiu, “Testing the time-reversal modified universality of the sivers function,” Physical Review Letters, vol. 103, Article ID 172001, 2009.

Z. B. Kang and J. W. Qiu, “Single transverse spin asymmetry of dilepton production near Z⁰ pole,” Physical Review D: Particles, Fields, Gravitation and Cosmology, vol. 81, Article ID 054020, 2010.

J. Collins, A. Efremov, K. Goeke et al., “Sivers effect in the Drell-Yan process at BNL RHIC,” Physical Review D Covering Particles, Fields, Gravitation, and Cosmology, vol. 73, Article ID 094023, 2006.

A. V. Efremov and O. V. Teryaev, “On spin effects in quantum chromodynamics,” Soviet Journal of Nuclear Physics, vol. 36, no. 1, pp. 242–246, 1982.

A. V. Efremov and O. V. Teryaev, “QCD asymmetry and polarized hadron structure function measurement,” Physics Letters B, vol. 150, no. 5, pp. 383–386, 1985.

J.-W. Qiu and G. F. Sterman, “Single transverse spin asymmetries,” Physical Review Letters, vol. 67, no. 17, pp. 2264–2267, 1991.

J.-W. Qiu and G. F. Sterman, “Single transverse spin asymmetries in direct photon production,” Nuclear Physics B, vol. 378, pp. 52–78, 1992.

J.-W. Qiu and G. F. Sterman, “Single transverse-spin asymmetries in hadronic pion production,” Physical Review D: Particles, Fields, Gravitation and Cosmology, vol. 59, Article ID 014004, 1999.

C. Kouvaris, J. W. Qiu, W. Vogelsang, and F. Yuan, “Single transverse-spin asymmetry in high transverse momentum pion production in pp collisions,” Physical Review D: Particles, Fields, Gravitation and Cosmology, vol. 74, Article ID 114013, 2006.

Y. Kanazawa and Y. Koike, “Polarization in hadronic Λ hyperon production and chiral-odd twist-3 distribution,” Physical Review D Covering Particles, Fields, Gravitation, and Cosmology, vol. 64, Article ID 034019, 2001.

H. Eguchi, Y. Koike, and K. Tanaka, “Twist-3 formalism for single transverse spin asymmetry reexamined: semi-inclusive deep inelastic scattering,” Nuclear Physics B, vol. 763, no. 1-2, pp. 198–227, 2007.

Y. Koike and K. Tanaka, “Master formula for twist-3 soft-gluon-pole mechanism to single transverse-spin asymmetry,” Physics Letters B, vol. 646, no. 5, pp. 232–241, 2007.

Y. Koike and K. Tanaka, “Universal structure of twist-3 soft-gluon-pole cross sections for single transverse-spin asymmetry,” Physical Review D: Particles, Fields, Gravitation and Cosmology, vol. 76, no. 6, Article ID 011502, 2007.

Z. B. Kang, F. Yuan, and J. Zhou, “Twist-three fragmentation function contribution to the single spin asymmetry in pp collisions,” Physics Letters B, vol. 691, pp. 243–248, 2010.

Z. B. Kang, A. Metz, J. W. Qiu, and J. Zhou, “Exploring the structure of the proton through polarization observables in lpjet X,” Physical Review D: Particles, Fields, Gravitation and Cosmology, vol. 84, Article ID 034046, 2011.

L. Gamberg and Z. B. Kang, “Single transverse spin asymmetry of prompt photon production,” Physics Letters B, vol. 718, no. 1, pp. 181–188, 2012.

H. Beppu, Y. Koike, K. Tanaka, and S. Yoshida, “Contribution of twist-3 multigluon correlation functions to single spin asymmetry in semi-inclusive deep inelastic scattering,” Physical Review D: Particles, Fields, Gravitation and Cosmology, vol. 82, no. 5, Article ID 054005, 2010.

Y. Koike and S. Yoshida, “Three-gluon contribution to the single spin asymmetry in Drell-Yan and direct-photon processes,” Physical Review D: Particles, Fields, Gravitation and Cosmology, vol. 85, Article ID 034030, 2012.

Y. Koike and S. Yoshida, “Probing the three-gluon correlation functions by the single spin asymmetry in p↑p→DX,” Physical Review D: Particles, Fields, Gravitation and Cosmology, vol. 84, Article ID 014026, 2011.

Z. B. Kang and J. W. Qiu, “Single transverse-spin asymmetry for D-meson production in semi-inclusive deep inelastic scattering,” Physical Review D: Particles, Fields, Gravitation and Cosmology, vol. 78, Article ID 034005, 2008.

Z. B. Kang, J. W. Qiu, W. Vogelsang, and F. Yuan, “Accessing trigluon correlations in the nucleon via the single spin asymmetry in open charm production,” Physical Review D: Particles, Fields, Gravitation and Cosmology, vol. 78, no. 11, Article ID 114013, 9 pages, 2008.

F. Yuan and J. Zhou, “Collins function and the single transverse spin asymmetry,” Physical Review Letters, vol. 103, Article ID 052001, 2009.

A. Metz and D. Pitonyak, “Fragmentation contribution to the transverse single-spin asymmetry in proton–proton collisions,” Physics Letters B, vol. 723, no. 4-5, pp. 365–370, 2013.

K. Kanazawa and Y. Koike, “Contribution of the twist-3 fragmentation function to the single transverse-spin asymmetry in semi-inclusive deep inelastic scattering,” Physical Review D: Particles, Fields, Gravitation and Cosmology, vol. 88, no. 7, Article ID 074022, 15 pages, 2013.

Z.-T. Liang, A. Metz, D. Pitonyak, A. Schafer, Y.-K. Song, and J. Zhou, “Double spin asymmetry ALT in direct photon production,” Physics Letters B, vol. 712, no. 3, pp. 235–239, 2012.

A. Metz, D. Pitonyak, A. Schaefer, and J. Zhou, “Analysis of the double-spin asymmetry ALT in inelastic nucleon-nucleon collisions,” Physical Review D Covering Particles, Fields, Gravitation, and Cosmology, vol. 86, Article ID 114020, 2012.

X. Ji, J. W. Qiu, W. Vogelsang, and F. Yuan, “Unified picture for single transverse-spin asymmetries in hard-scattering processes,” Physical Review Letters, vol. 97, Article ID 082002, 2006.

X. Ji, J. W. Qiu, W. Vogelsang, and F. Yuan, “Single transverse-spin asymmetry in Drell-Yan production at large and moderate transverse momentum,” Physical Review D Covering Particles, Fields, Gravitation, and Cosmology, vol. 73, Article ID 094017, 2006.

X. Ji, J. W. Qiu, W. Vogelsang, and F. Yuand, “Single transverse-spin asymmetry in semi-inclusive deep inelastic scattering,” Physics Letters B, vol. 638, no. 2-3, pp. 178–186, 2006.

Y. Koike, W. Vogelsang, and F. Yuan, “On the relation between mechanisms for single-transverse-spin asymmetries,” Physics Letters B, vol. 659, no. 5, pp. 878–884, 2008.

A. Bacchetta, D. Boer, M. Diehl, and P. Mulders, “Matches and mismatches in the descriptions of semi-inclusive processes at low and high transverse momentum,” Journal of High Energy Physics, vol. 023, 2008.

D. Boer, Z. Kang, W. Vogelsang, and F. Yuan, “Test of the universality of naive-time-reversal-odd fragmentation functions,” Physical Review Letters, vol. 105, no. 20, 2010.

Z.-B. Kang, J.-W. Qiu, W. Vogelsang, and F. Yuan, “Observation concerning the process dependence of the sivers functions,” Physical Review D: Particles, Fields, Gravitation and Cosmology, vol. 83, Article ID 094001, 2011.

L. Gamberg and Z.-B. Kang, “Process dependent Sivers function and implication for single spin asymmetry in inclusive hadron production,” Physics Letters B, vol. 696, no. 1-2, pp. 109–118, 2011.

J. Adams, “Cross sections and transverse single-spin asymmetries in forward neutral-pion production from proton collisions at √s=200  GeV,” Physical Review Letters, vol. 92, Article ID 171801, 2004.

B. I. Abelev, M. M. Aggarwal, Z. Ahammed et al., “Forward neutral-pion transverse single-spin asymmetries in p+p collisions at √s=200  GeV,” Physical Review Letters, vol. 101, Article ID 222001, 2008.

S. Heppelmann, “Measurement of transverse single spin asymmetry A_N in eta mass region at large feynman X_F with the star forward pion detector,” Nuclear Experiment, 2009.

S. S. Adler, S. Afanasiev, and C. Aidala, “Measurement of transverse single-spin asymmetries for midrapidity production of neutral pions and charged hadrons in polarized collisions at GeV,” Physical Review Letters, vol. 95, Article ID 202001, 2005.

I. Arsene, I. G. Bearden, D. Beavis et al., “Single-transverse-spin asymmetries of identified charged hadrons in polarized pp collisions at √s=62.4  GeV,” Physical Review Letters, vol. 101, Article ID 042001, 2008.

K. Kanazawa, Y. Koike, A. Metz, and D. Pitonyak, “Towards an explanation of transverse single-spin asymmetries in proton-proton collisions: the role of fragmentation in collinear factorization,” Physical Review D: Particles, Fields, Gravitation and Cosmology, vol. 89, Article ID 111501, 2014.

L. Gamberg, Z. Kang, D. Pitonyak, and A. Prokudin, “Phenomenological constraints on A N in p ↑ p → π X from Lorentz invariance relations,” Physics Letters B, vol. 770, pp. 242–251, 2017.

Z. B. Kang and J. W. Qiu, “Evolution of twist-3 multiparton correlation functions relevant to single transverse-spin asymmetry,” Physical Review D Covering Particles, Fields, Gravitation, and Cosmology, vol. 79, Article ID 016003, 2009.

J. Zhou, F. Yuan, and Z. T. Liang, “QCD evolution of the transverse momentum dependent correlations,” Physical Review D Covering Particles, Fields, Gravitation, and Cosmology, vol. 79, Article ID 114022, 2009.

V. M. Braun, A. N. Manashov, and B. Pirnay, “Scale dependence of twist-three contributions to single spin asymmetries,” Physical Review D: Particles, Fields, Gravitation and Cosmology, vol. 80, Article ID 114002, 2009.

A. Schafer, J. Zhou, and D. Phys. Rev, “Note on the scale evolution of the Efremov-Teryaev-Qiu-Sterman function TF(x,x),” Physical Review D Covering Particles, Fields, Gravitation, and Cosmology, vol. 85, Article ID 117501, 2012.

J. P. Ma and Q. Wang, “Scale dependence of twist-3 quark–gluon operators for single spin asymmetries,” Physics Letters B, vol. 715, no. 1-3, pp. 157–163, 2012.

Z. B. Kang and J. W. Qiu, “QCD evolution of naive-time-reversal-odd parton distribution functions,” Physics Letters B, vol. 713, no. 3, pp. 273–276, 2012.

Z. B. Kang, “QCD evolution of naive-time-reversal-odd fragmentation functions,” Physical Review D Covering Particles, Fields, Gravitation, and Cosmology, vol. 83, Article ID 036006, 2011.

J. P. Ma and G. P. Zhang, “Evolution of chirality-odd twist-3 fragmentation functions,” Physics Letters B, vol. 772, pp. 559–566, 2017.

W. Vogelsang and F. Yuan, “Next-to-leading order calculation of the single transverse spin asymmetry in the Drell-Yan process,” Physical Review D Covering Particles, Fields, Gravitation, and Cosmology, vol. 79, Article ID 094010, 2009.

Z. B. Kang, I. Vitev, H. Xing, and D. Phys. Rev, “ransverse momentum-weighted Sivers asymmetry in semi-inclusive deep inelastic scattering at next-to-leading order,” Physical Review D: Particles, Fields, Gravitation and Cosmology, vol. 87, Article ID 034024, 2013.

S. Yoshida, “New pole contribution to Ph⊥-weighted single-transverse spin asymmetry in semi-inclusive deep inelastic scattering,” Physical Review D: Particles, Fields, Gravitation and Cosmology, vol. 93, Article ID 054048, 2016.

L. Dai, Z. Kang, A. Prokudin, and I. Vitev, “Next-to-leading order transverse momentum-weighted Sivers asymmetry in semi-inclusive deep inelastic scattering: The role of the three-gluon correlator,” Physical Review D: Particles, Fields, Gravitation and Cosmology, vol. 92, no. 11, Article ID 114024, 2015.

A. P. Chen, J. P. Ma, and G. P. Zhang, “One-loop corrections to single spin asymmetries at twist-3 in Drell-Yan processes,” Physical Review D Covering Particles, Fields, Gravitation, and Cosmology, vol. 95, Article ID 074005, 2017.

A. P. Chen, J. P. Ma, and G. P. Zhang, “One-loop corrections of single spin asymmetries in semi-inclusive DIS,” Physical Review D: Particles, Fields, Gravitation and Cosmology, vol. 97, Article ID 054003, 2018.

M. G. Alexeev, G. D. Alexeev, A. Amoroso et al., “Measurement of PT-weighted Sivers asymmetries in leptoproduction of hadrons,” Nuclear Physics B, vol. 940, pp. 34–53, 2019.

Y. Koike and K. Tanaka, “Master formula for twist-3 soft-gluon-pole mechanism to single transverse-spin asymmetry,” Physics Letters B, vol. 668, pp. 458-459, 2008.

K. Kanazawa, Y. Koike, A. Metz, D. Pitonyak, and M. Schlegel, “Operator constraints for twist-3 functions and Lorentz invariance properties of twist-3 observables,” Physical Review D Covering Particles, Fields, Gravitation, and Cosmology, vol. 5, Article ID 054024, 2016.

Z. B. Kang, E. Wang, X. N. Wang, and H. Xing, “Transverse momentum broadening in semi-inclusive deep inelastic scattering at next-to-leading order,” Physical Review D Covering Particles, Fields, Gravitation, and Cosmology, vol. 94, no. 11, Article ID 114024, 2016.

Z. B. Kang, J. W. Qiu, X. N. Wang, and H. Xing, “Next-to-leading order transverse momentum broadening for Drell-Yan production in p+A collisions,” Physical Review D: Particles, Fields, Gravitation and Cosmology, vol. 94, no. 7, Article ID 074038, 2016.

Z. B. Kang, E. Wang, X. N. Wang, and H. Xing, “Next-to-leading order QCD factorization for semi-inclusive deep inelastic scattering at twist 4,” Physical Review Letters, vol. 112, no. 10, Article ID 102001, 2014.