Advances in High Energy Physics

Volume 2019, Article ID 4825790, 15 pages

https://doi.org/10.1155/2019/4825790

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

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; nc.ude.uncs.m@gnixh and Shinsuke Yoshida; moc.liamg@58adihsoynihs

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^{3}.

#### 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 asFor clear understanding, we start with the -integrated cross section at leading twist in unpolarized lepton-proton scatteringThere 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 aswhere 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 sectionwhich 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 functionThe 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 aswhere 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 byWe 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 aswhere the hard parts are given byThe -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 asThe 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 asCombining (14) and (20), we can obtain the resultWe can find that this matrix element corresponds to term in the first moment of the TMD correlatorwhere 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 aswhere . 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.