Single Transverse-Spin Asymmetries in Drell-Yan Processes

The asymmetry in the angular distribution of Drell--Yan dilepton pairs generated in collisions where just one nucleon is transversely polarised has been examined in the literature with a variety of results, mainly differing by factors of two. We re-evaluate the asymmetry via twist-3 contributions in collinear factorisation. In order to allow complete and in-depth comparison with existing calculations, we supply all calculational details.


Introduction
In this paper we focus our attention on the study of the cross-section for Drell-Yan [1] processes in which an unpolarised proton and a transversely polarised antiproton annihilate to produce a dilepton pair. In particular, we shall calculate the single-spin asymmetry (SSA) defined as where S T is the polarisation of the antiproton in the plane orthogonal to its motion. The study of this asymmetry is of particular interest for a variety of reasons. Firstly, a number of experiments propose to measure it. Secondly, given the small number of quarks involved in such processes, SSAs are ideal observables to further test factorisation mechanisms in pQCD. Lastly, SSAs are related to particular parton distribution functions that describe the inner structure of hadrons. a e-mail: philip.ratcliffe@uninsubria.it The motivation to review this specific topic is the diversity of the results found in literature. The first calculation was performed in [2], where the authors obtained the following result: where f q is the twist-two (unpolarised) distribution for a quark of flavour q, the momentum fraction of the antiquark is defined by x ′ = Q 2 /(xs), s = (P + P ′ ) 2 is the CM frame energy squared and Q 2 the momentum squared of the virtual photon. The correlator T F (x, x) is a twist-three quark-quark-gluon spin-dependent matrix element (given e.g. in [2]), which will be described in detail later. The derivative term in the numerator of (2) was later doubted and the result given in [3,4] is The question of the derivative was further addressed in [5] without clear conclusions while in [6] its absence was argued owing to the absence of so-called double soft-gluon poles in this calculation. Result (3) was obtained by expanding the hadronic tensor in terms of twist-3 correlators depending on the hadron transverse spin. A parallel approach may be found in [7], in which factorisation is in terms of transverse momentum dependent quark-quark correlators. The final result is, however, in contrast with the previous two results. 2 A further attempt was made in [8], in which result (3) was found to be multiplied by a factor two. The SSA was recently recalculated in [9], where the authors noted that transverse-momentum dependence must be taken into account not only in the hadronic tensor, but also in the leptonic tensor; indeed, p T flows from the hadronic to the leptonic part via the virtual photon. The result found there was half of (3). In [10] the calculation was performed using partonic states and the result obtained coincides with that of [9]. The motivation given there for the difference with respect to the earlier papers appears dissimilar in form although the interplay between the transverse-momentum dependence in both the hadronic and leptonic tensors is central.
In order to attempt to clarify the situation, we have repeated the calculation taking into account the observations of [9]; we obtain, however, a result which is two times (3). In this paper we shall define all the partonic distribution functions and quantities used and show all calculational details.
The rest of the paper is organised as follows. In the next section we provide the definitions of the kinematical variables used and we define the quark-quark and quark-quark-gluon correlators with respect to the associated parton distributions. In section 3 we describe the Collins-Soper (CS) reference frame [11], which we use in the calculation, and the transverse-spin vector in these coordinates. We then focus our attention on the question of gauge invariance of the correlators and suppression of gauge-link operators. In the final two sections we examine the collinear expansion and show all calculational details involved in reaching the final result.

Notation and Definitions
In this section we review the notation and definitions used in this paper. Let us start by noting that in the ultrarelativistic regime the momentum carried by a quark is essentially in the hadron direction, or rather, in the beam direction. Therefore, given a pair of vectors: taking the third component along the beam direction, we may define and their light-cone projections The specific value of Λ will determine the chosen hadron reference frame.
Consequently, a generic vector may be written as P µ = (p + , p − , P T ), with P T = (p 1 , p 2 ) while a generic scalar product is P ·P ′ = P + P ′− + P − P ′+ − P T ·P ′ T . The power of this notation lies in the fact that, by choosing the beam direction as p + , as in a light-cone coordinate system, this component will dominate in scalar products and we may easily isolate effects induced by transverse spin and momentum. With these definitions, we may write the hadron momentum as This relation is true not only in the hadron rest frame (Λ = M/ √ 2, with M the hadron mass), but also in the infinite-momentum frame (Λ = P + → ∞).
In the second part of this section we define one of the fundamental objects used in the study of hadronic scattering processes: the quark-quark correlator. Let us start by noting that the transition matrices that describe the passage from an hadronic state with definite momentum and spin |P, S to a generic state |X are of the form X| ψ i |P, S . As a consequence, to obtain a transition probability, we need to multiply by the complex conjugate and sum over all intermediate states. It is therefore natural to define the quark-quark correlator Correlators written in this form are not measurable in physical processes and our interest will therefore be to connect this object with observables such as the hadron momentum and spin. Firstly, we note that the correlator is a 4×4 Dirac matrix, a Lorentz scalar and depends on three independent physical quantities [P, S, p], where P is the hadron momentum, S its spin and p the quark momentum. In parallel, there are only three independent matrices, they are [½, γ µ , γ 5 ]. Let us now multiply this set of vectors and matrices by each other respecting the following properties: -Invariance under parity. Parity changes the signs of momenta, but not of axial vectors such as spin. As a consequence, helicity, defined as λ = P·S |P| , flips sign under such transformations. For parity to be conserved, the following must be valid: φ(P, S, p) = γ 0 φ(P , −S,p) γ 0 .
-Invariance under charge conjugation and time reversal. This condition implies that: where C = iγ 0 γ 2 is the charge-conjugation operator. This is obvious since the transformation does not alter the spatial part of the spin operator and helicity does not change sign. -The correlator must be Hermitian, which implies: -The expansion must be linear in spin. This condition arises from the fact that Lorentz invariance implies that the hadronic tensor must be linear in the spinors u(P, S) andū(P, S). As a consequence, the tensors constructed with this requirement are spin independent or linear in spin: The set of combinations of vectors based on these requirements is as follows: Note that we need not consider all possible products between p, P and S, but only those linearly dependent on the spin S. This is because the coefficients A i multiplying each structure in the complete expansion will depend on p·P and p 2 .
Let us now construct the most general expansion of the correlator over the basis obtained by multiplying the elements of the set of vectors and matrices. We must pay attention to the condition mentioned earlier. The maximum number of products is limited by the fact that the product of a quantity with itself does not generate any new Dirac structure. In fact, we have ½ 2 = ½, forbidden by the request of linearity in spin. As a consequence, if r is the number of basic elements multiplying each other, we obtain r = 1 : ½, / p, / P , / Sγ 5 , p·Sγ 5 ; (14a) r = 2 : / p / P , / p/ Sγ 5 , p·S / pγ 5 , / P / Sγ 5 , p·S / P γ 5 ; (14b) The method produces 12 independent structures; not all, however, fulfill the requirements of time-reversal invariance and hermiticity. To see which satisfy these conditions it is convenient to find out which of these terms include a double or triple slash. To do this we use the relations and in which the only interesting terms are the first as the others have already been taken into account (recall that a·b is proportional to ½).
We now have the correlator expansion: In parallel, the correlator is a matrix in Dirac space; it may therefore be expanded over an orthonormal basis of γ matrices to obtain: where the letters S, P , V , A and T indicate the types of currents: scalar, pseudoscalar, vector, axial and tensor respectively. Comparing expressions (16) and (17), we obtain We may simplify these expressions, noting in this case that P µ ∼ P + , p µ = xP µ , M S µ = λP µ +M S µ T and we may neglect terms in M 2 /P + . As a consequence, neglecting quark transverse momentum, we can rewrite these equations in terms of the proton momentum and transverse spin only. Let us define new coefficients D i 4 depending only on the quark momentum fraction x. Keeping only the dominant terms, we redefine: As a consequence, it is possible to rewrite the correlator in terms of the structure functions f (x), ∆f (x) and ∆ T f (x). We obtain: Thus far we have worked at the twist-two level. Alternatively, following [12], it is possible to project the Dirac spinors onto 'good' and 'bad' components and show that they are connected by gluonic fields whenever the quark involved in the process is off-shell. In particular, defining the Hermitian projectors È ± = 1 The following relation is also true: We observe that the relation between '+' and '−' components depends on a term containing the covariant derivative, and thus on the Dirac equation. Therefore, if the quark is on-shell, it is possible to rewrite ψ − in terms of ψ + and the quark mass, but if the quark is off-shell there are no relations between 'good' and 'bad' components and the addition of new degrees of freedom that are not eliminable by a different choice of gauge is necessary. These are given by the interaction of the transverse gluon field with the quark. In order to isolate these dependences and create general twist-3 distributions, it is necessary to insert the n components into the Sudakov decomposition made earlier. Using relation (7) together with the expansion in (20), with convenient redefinitions, we obtain the following set of twist-2 to 4 distribution functions [13]: 1 1 Note that in [13] the normalisation S 2 = −M 2 is adopted.
x y P, S Figure 1 The Born diagram for quark-quark-gluon correlation matrix.
Wishing now to also display the dependence of the twist-3 correlators on the gluonic fields, we expand over independent distribution functions. Let us start by redefining the quark-quark-gluon correlation matrix In general this matrix will be associated with the Born diagram shown in Fig. 1.
In this case x and y represent the momentum fractions carried by the quarks. More precisely, x is the momentum fraction carried by the quark on the left of the diagram and y on the right. Consequently, x − y is the momentum fraction carried by the gluon. From now on an integral is implied over the quark and gluon momenta with the restrictions δ(xP −p) and δ(yP −p ′ ), where p and p ′ are the quark momenta.
We observe that this matrix contains a covariant derivative and thus separates into two parts. The first contains the operatorψ∂ α ψ and the second the gluonic fieldψA α ψ. As already seen for the twist-two distribution, in order to obtain a twist-τ distribution it is only necessary to consider τ partons with dynamically independent polarisations. If this independence does not hold, new higher-order contributions will be generated. In order to respect this independence therefore, we project out the "+" components from such a correlation matrix and, for the gauge choice A + = 0, only transverse gluon polarisations survive.
To perform this projection, we insert / n and restrict the index α to be transverse. Consequently, we may define the distributions G(x, y) andG(x, y) for the vector and axial-vector projections of φ α Dij : where ǫ αβ T ≡ ǫ +−αβ and the dots denote higher-twist contributions. Note with our normalisation G(x, y) and G(x, y) have dimensions of a mass.
The origin of these two possible choices of expansion over momentum and spin vectors lies in the possible gluon polarisation degrees of freedom. In fact, in the first definition the gluon polarisation is orthogonal to the hadron spin and momentum, in the second it is aligned with the spin. Demanding hermiticity of the operators, we find that G(x, y) is antisymmetric under the interchange x ↔ y whereasG(x, y) is symmetric. Note moreover that invariance under time reversal requires both to be real. Analogously, considering also the tensor interaction, we may define: × P, S|ψ(0)σ µν γ 5 iD α T (ηn)/ nψ(λn) |P, S ≡ g µα p ν + g να p µ S·n H(x, y) The relations between definitions (25a-c) and (20) have been studied extensively in the literature, in particular in [3]. We have the following: and M e(x) = 1 x dy E(x, y) − E(y, x) .

Reference frame and transverse momentum
In this section we focus our attention on the choice of the reference frame. As the polarisation is defined in the LAB frame transversely to the beam direction, it is natural to choose the Collins-Soper reference frame. With this choice we can isolate transverse-momentum effects, not only in the hadronic tensor, but also in the leptonic tensor. This frame is defined by the following: Let us give a geometric interpretation of this frame: ξ represents the photon momentum fraction in the beam direction. Thus, subtracting ξq µ from the proton momentum, we find all the dependence on the transverse momentum of the struck quark. As for the basis vectors, Z indicates the direction of the beam in the absence of transverse momentum, X the transverse-momentum direction, T the photon direction and Y , again transverse, completes the set of orthonormal vectors. A simplified representation of this frame is displayed in Fig. 2.
In this frame we can also define the angles we shall use in the calculation of the SSA. In particular, θ represents the angle between the Z axis and the direction of the outgoing leptons, φ 1 is the angle between the lepton plane and the transverse-spin direction, and φ is the angle between the lepton plane and the X axis.
As we have seen, the polarisation of the hadron is defined in the LAB frame, we thus now wish to show how the spin vector transforms on moving from the LAB to the CS frame. Let us start by noting that this transformation is possible with two Lorentz boosts in succession. The first is in the beam direction in order to set Q 3 to zero, the second is in the Q T direction to cancel this component too and leave only the temporal component Q. After this transformation it is clear that Q 2 corresponds to the mass of the dilepton pair. We represent the spin vector before the boost as: Applying the two boosts we obtain where we may however neglect terms in Q 2 T . We thus see that in this frame the spin vector acquires a temporal component. However, it is easy to show that, after contraction with the leptonic tensor, this has no physical consequence.

Gauge invariance in twist-three correlators
We now analyse the question of gauge invariance. The gauge choice A + = 0, is made for two reasons: firstly, examining expression (8), we note that z is a spacetime variable. Therefore, in order to connect the two different space-time points, we must insert a gauge-link operator into the correlator in the following way where the gauge link is The choice A + = 0 sets this operator to ½, as the dominant direction of z is the "+" direction. Secondly, in this gauge it is simple to transform a correlator of the form into a gauge-invariant expression. This may be achieved by noting that and thus, after integration by parts, we arrive at where φ α F is obtained by replacing A α with F +α in Eq. (35). Using this simple rule, we may then make the substitutions G A (x, y) → G F (x, y) andG A (x, y) → G F (x, y), where G F etc. are defined by the same replacements in (25a) and (25b). 7 Finally, in order to compare clearly with our calculation, we note that the matrix element T F appearing in the asymmetry reads × P S|ψ(0)γ + ǫ µν T S T ν gF + µ (ηn)ψ(λn) |P S . (38) Thus, from the structure of Eqs. (25a), (37) and (38), we see that the function G F is identical to T F , having the same dependence on transverse spin, with the substitution A α → F +α .

General form of the cross-section
We begin this section by giving the standard expression for the Drell-Yan cross-section: The factorisation theorem provides the possibility to divide the hard-scattering part from the soft part. Thus, defining φ 1 and φ 2 , the correlators for the proton and the antiproton, we may write the hadronic tensor in the following general form: × Tr φ 1 γ µφ 2 γ ν . (40) For the case in which we consider an extra exchanged gluon, we must include the gluon propagator in the trace, as we shall in the next section. Note that this general form contains δ 4 (p 1 + p 2 − q), which includes both the longitudinal and transverse components.
We now perform the collinear expansion of the crosssection. To do this, we expand in a neighbourhood of momentum transverse to the direction of the hadron motion. Calling this momentum p T , we note that in the CS frame its only non-zero components are along the X and Y directions and it is thus of the form αX σ + βŶ σ . Understanding the derivative as a gradient, we have: x ȳ x Figure 3 The twist-3 contribution to the Drell-Yan process (to which must also be added the Hermitian conjugate).
We therefore now need to differentiate the hadronic tensor with respect to the transverse momentum: which, in light-cone coordinates, is equivalent to Integrating the first member by parts gives the second with the opposite sign and so they cancel exactly. Thus, no linear contributions in p T to the DY process survive and therefore the derivative found in [2] is absent.

Calculation of the single-spin asymmetry
In this section we show all details of our calculation. We start by noting that, via the optical theorem, we need an imaginary part to obtain a non-zero SSA. Such a contribution can only come from diagrams that involve a single gluon exchange between the polarised hadron and the quark coming from the other hadron. The diagram relevant here is shown in Fig. 3. We now write an analytic expression for the hadronic tensor using the factorisation theorem. For the soft parts we write the product of two correlators defined as in (8) and (35).
For the hard part we see on-shell propagation of the quark before the photon production and must thus insert a quark propagator before the first vertex. Considering the diagram in Fig. 3, we may write: 8 The relevant contribution to the denominator here is 2x(y − x)P + 1 P − 2 =x(y − x)s. In the numerator we need only take the '+' component since the gluon comes from the proton and, using the fact that the LAB and CM frames correspond, (P + 1 ) 2 = (P − 2 ) 2 and therefore P + 1 = s 2p + . The propagator term then becomes Similarly, for the Hermitian conjugate of the diagram in Fig. 3 the propagator becomes the complex conjugate and the role of x and y are interchanged. We thus have: Only these two terms contribute to the SSA owing to the pole in the denominator at x = y, or simply when the gluon carries zero momentum. This source of imaginary part can be regularised via Cauchy's theorem and used for the SSA calculation [14,15].
Let us now examine the complete expression for the hadronic part. We must evaluate the integral The first term comes from the diagram with no gluon exchange and contributes to the denominator of the SSA. In order to evaluate it we make the expansion Recalling that / P 1 = P + 1 γ − , / P 2 = P − 2 γ + , / p = γ − and P + 1 P − 2 = s/2, we finally obtain: Note that the only the first term is symmetric and thus alone survives contraction with the leptonic tensor.
We now evaluate the second term of (47), containing the twist-3 correlator. For a polarised hadron the full expansion for the correlator becomes (recall that α is a transverse index), while for the unpolarised hadron we only have the vector term f (x) / P 1 . Using again / P 1 = P + 1 γ − , / P 2 = P − 2 γ + and / p + = γ − , we then have Evaluating traces and collecting the Levi-Civita tensor, we obtain the final expression Finally, we evaluate the third and last term of (47). Following the same procedure and summing this with other results, we obtain the complete expression: With our choice of gauge and using (37), this is equivalent to the gauge-invariant expression: 9 where only the first two terms are symmetric and thus survive contraction with the leptonic part. We now regularise the pole via Cauchy's theorem and take the imaginary part to obtain the following expression: Note that the δ-function sets the term inside square brackets to zero. Adding the first term of (47), performing the sum over colours and reinstating the factor e 2 /N c , we finally obtain: The final operation to perform is contraction with the leptonic tensor, following formula (41). We start by evaluating the first term in the collinear expansion (41), i.e. with p T identically zero: Inserting now Eq. (57) into the second term of formula (41), recalling that the derivative term vanishes and noting that only the Z {µ X ν} and Z {µ Y ν} terms in L µν survive, we find that this contribution is identical to the second term of Eq. (58) and thus implies a factor two. With this final observation, to first order in the collinear expansion, we have: which leads to the following SSA: sin 2θ sin φ 1 1 + cos 2 θ a e 2 a f a (x) G F a (x, x) a e 2 a f a (x) f a (x) .

Conclusions
As in [9,10], we find it crucial in our calculation to perform the collinear expansion with due respect to all possible ensuing transverse-momentum dependence. However, the extra contribution generated is of the same sign and magnitude as the naïve part, thus leading to an overall factor of two with respect to the older calculations and four for the later. Unfortunately, the level of detail provided in [8] renders full in-depth comparison impossible. We note too that, while an explanation is proposed in [9], the level of detail is again not sufficient to permit comparison. With regard to [10], the origin of the extra piece there does not appear to be the same as that of [9]. Indeed, the authors appear to claim that in earlier papers the error lies in incorrect expansion of the hadronic tensor whereas the authors of [9] point to failure to expand the leptonic tensor.
In conclusion, we have repeated the calculation of the transverse SSA in DY, including the contribution arising from collinear expansion of the leptonic tensor, as suggested in [9], and find a factor of four with respect to the most recent calculations presented in [9] and [10]; we thus find ourselves in agreement with the results of [8]. Unfortunately however, with the published detail available, we cannot point clearly to possible sources of error in the other papers.