The $\sin (2\phi-\phi_S)$ azimuthal asymmetry in the pion induced   Drell-Yan process within TMD factorization

Hui Li; Xiaoyu Wang; Zhun Lu

arXiv:1907.07095·hep-ph·March 18, 2020

The $\sin (2\phi-\phi_S)$ azimuthal asymmetry in the pion induced Drell-Yan process within TMD factorization

Hui Li, Xiaoyu Wang, Zhun Lu

PDF

TL;DR

This paper predicts the $ ext{sin}(2 heta- heta_S)$ asymmetry in pion-induced Drell-Yan processes using TMD factorization, modeling distributions with light-cone wavefunctions and comparing results with COMPASS data.

Contribution

It introduces a detailed TMD-based prediction of the single transverse-spin asymmetry in pion-induced Drell-Yan, incorporating nonperturbative evolution effects and model distributions.

Findings

01

Prediction of asymmetry consistent with COMPASS data in sign and magnitude.

02

Demonstration of the impact of TMD evolution on asymmetry estimates.

03

Numerical estimates of asymmetry as functions of kinematic variables.

Abstract

We investigate the single transverse-spin asymmetry with a $sin (2 ϕ - ϕ_{S})$ modulation in the pion-induced Drell-Yan process within the theoretical framework of the transverse momentum dependent factorization. The asymmetry is contributed by the convolution of the Boer-Mulders function and the transversity. We adopt the model results for the distributions of the pion meson from the light-cone wavefunction approach and the available parametrization for the distributions of the proton to numerically estimate the $sin (2 ϕ - ϕ_{S})$ asymmetry in $π^{-} p$ Drell-Yan at the kinematics of the COMPASS at CERN. To implement the TMD evolution formalism of parton distribution functions, we apply the recently extracted nonperturbative Sudakov form factor associated with the distribution functions of the proton and the pion. It is found that our prediction on the single transverse-spin…

Equations96

\frac{\partial ln F ~ ( x , b ; μ , ζ _{F} )}{\partial ζ _{F}} = \tilde{K} (b; μ),

\frac{\partial ln F ~ ( x , b ; μ , ζ _{F} )}{\partial ζ _{F}} = \tilde{K} (b; μ),

\frac{d K ~}{d ln μ} = - γ_{K} (α_{s} (μ)),

\frac{d K ~}{d ln μ} = - γ_{K} (α_{s} (μ)),

\frac{d ln F ~ ( x , b ; μ , ζ _{F} )}{d ln μ} = γ_{F} (α_{s} (μ); \frac{ζ _{F}^{2}}{μ ^{2}}),

\tilde{F} (x, b, Q) = F \times e^{- S (Q, b)} \times \tilde{F} (x, b, μ_{i}),

\tilde{F} (x, b, Q) = F \times e^{- S (Q, b)} \times \tilde{F} (x, b, μ_{i}),

S (Q, b) = S_{P} (Q, b) + S_{NP} (Q, b) .

S (Q, b) = S_{P} (Q, b) + S_{NP} (Q, b) .

S_{P} (Q, b) = \int_{μ_{b}^{2}}^{Q^{2}} \frac{d μ ˉ ^{2}}{μ ˉ ^{2}} [A (α_{s} (\overset{μ}{ˉ})) ln \frac{Q ^{2}}{μ ˉ ^{2}} + B (α_{s} (\overset{μ}{ˉ}))],

S_{P} (Q, b) = \int_{μ_{b}^{2}}^{Q^{2}} \frac{d μ ˉ ^{2}}{μ ˉ ^{2}} [A (α_{s} (\overset{μ}{ˉ})) ln \frac{Q ^{2}}{μ ˉ ^{2}} + B (α_{s} (\overset{μ}{ˉ}))],

A = n = 1 \sum \infty A^{(n)} (\frac{α _{s}}{π})^{n},

A = n = 1 \sum \infty A^{(n)} (\frac{α _{s}}{π})^{n},

B = n = 1 \sum \infty B^{(n)} (\frac{α _{s}}{π})^{n} .

A^{(1)}

A^{(1)}

A^{(2)}

B^{(1)}

S_{NP} (Q; b) = g_{2} (b) ln Q / Q_{0} + g_{1} (b) .

S_{NP} (Q; b) = g_{2} (b) ln Q / Q_{0} + g_{1} (b) .

S_{NP} = g_{1} b^{2} + g_{2} ln \frac{b}{b _{*}} ln \frac{Q}{Q _{0}} + g_{3} b^{2} ((x_{0} / x_{1})^{λ} + (x_{0} / x_{2})^{λ}) .

S_{NP} = g_{1} b^{2} + g_{2} ln \frac{b}{b _{*}} ln \frac{Q}{Q _{0}} + g_{3} b^{2} ((x_{0} / x_{1})^{λ} + (x_{0} / x_{2})^{λ}) .

S_{NP}^{q} (Q, b) + S_{NP}^{\overset{q}{ˉ}} (Q, b) = S_{NP} (Q, b),

S_{NP}^{q} (Q, b) + S_{NP}^{\overset{q}{ˉ}} (Q, b) = S_{NP} (Q, b),

S_{NP}^{f_{1, q / p}} (Q, b) = \frac{g _{1}}{2} b^{2} + \frac{g _{2}}{2} ln \frac{b}{b _{*}} ln \frac{Q}{Q _{0}} .

S_{NP}^{f_{1, q / p}} (Q, b) = \frac{g _{1}}{2} b^{2} + \frac{g _{2}}{2} ln \frac{b}{b _{*}} ln \frac{Q}{Q _{0}} .

S_{NP}^{f_{1, q / π}} = g_{1}^{π} b^{2} + g_{2}^{π} ln \frac{b}{b _{*}} ln \frac{Q}{Q _{0}},

S_{NP}^{f_{1, q / π}} = g_{1}^{π} b^{2} + g_{2}^{π} ln \frac{b}{b _{*}} ln \frac{Q}{Q _{0}},

\tilde{F}_{q / H} (x, b; μ) = i \sum C_{q \leftarrow i} \otimes F_{i / H} (x, μ) .

\tilde{F}_{q / H} (x, b; μ) = i \sum C_{q \leftarrow i} \otimes F_{i / H} (x, μ) .

C_{q \leftarrow i} \otimes F_{i / H} (x, μ) \equiv \int_{x}^{1} \frac{d ξ}{ξ} C_{q \leftarrow i} (x / ξ, b; μ) F_{i / H} (ξ, μ),

C_{q \leftarrow i} \otimes F_{i / H} (x, μ) \equiv \int_{x}^{1} \frac{d ξ}{ξ} C_{q \leftarrow i} (x / ξ, b; μ) F_{i / H} (ξ, μ),

\tilde{F}_{q / H} (x, b; Q) = e^{- \frac{1}{2} S_{P} (Q, b_{*}) - S_{NP}^{F_{q / H}} (Q, b)} F (α_{s} (Q)) i \sum C_{q \leftarrow i}^{F} \otimes F_{i / H} (x, μ_{b}),

\tilde{F}_{q / H} (x, b; Q) = e^{- \frac{1}{2} S_{P} (Q, b_{*}) - S_{NP}^{F_{q / H}} (Q, b)} F (α_{s} (Q)) i \sum C_{q \leftarrow i}^{F} \otimes F_{i / H} (x, μ_{b}),

S_{P}^{q} (Q, b_{*}) = S_{P}^{\overset{q}{ˉ}} (Q, b_{*}) = S_{P} (Q, b_{*}) /2.

S_{P}^{q} (Q, b_{*}) = S_{P}^{\overset{q}{ˉ}} (Q, b_{*}) = S_{P} (Q, b_{*}) /2.

\tilde{f}_{1, q / p} (x, b; Q)

\tilde{f}_{1, q / p} (x, b; Q)

\tilde{f}_{1, q / π} (x, b; Q)

f_{1, q / p} (x, k_{⊥}; Q) = \int_{0}^{\infty} \frac{d bb}{2 π} J_{0} (k_{⊥} b) e^{- \frac{1}{2} S_{P} (Q, b_{*}) - S_{NP}^{f_{1, q / p}} (Q, b)} f_{1, q / p} (x, μ_{b}),

f_{1, q / p} (x, k_{⊥}; Q) = \int_{0}^{\infty} \frac{d bb}{2 π} J_{0} (k_{⊥} b) e^{- \frac{1}{2} S_{P} (Q, b_{*}) - S_{NP}^{f_{1, q / p}} (Q, b)} f_{1, q / p} (x, μ_{b}),

f_{1, q / π} (x, k_{⊥}; Q) = \int_{0}^{\infty} \frac{d bb}{2 π} J_{0} (k_{⊥} b) e^{- \frac{1}{2} S_{P} (Q, b_{*}) - S_{NP}^{f_{1, q / π}} (Q, b)} f_{1, q / π} (x, μ_{b}),

h_{1, q / p} (x, b; Q) = e^{- \frac{1}{2} S_{P} (Q, b_{*}) - S_{NP}^{f_{1, q / p}} (Q, b)} h_{1, q / p} (x, μ_{b}),

h_{1, q / p} (x, b; Q) = e^{- \frac{1}{2} S_{P} (Q, b_{*}) - S_{NP}^{f_{1, q / p}} (Q, b)} h_{1, q / p} (x, μ_{b}),

h_{1, q / p} (x, k_{⊥}; Q) = \int_{0}^{\infty} \frac{d bb}{2 π} J_{0} (k_{⊥} b) e^{- \frac{1}{2} S_{P} (Q, b_{*}) - S_{NP}^{f_{1, q / p}} (Q, b)} h_{1, q / p} (x, μ_{b}),

h_{1, q / π}^{α ⊥} (x, b; μ) = (\frac{- i b _{⊥}^{α}}{2}) T_{q / π, F}^{(σ)} (x, x; μ),

h_{1, q / π}^{α ⊥} (x, b; μ) = (\frac{- i b _{⊥}^{α}}{2}) T_{q / π, F}^{(σ)} (x, x; μ),

\tilde{h}_{1, q / π}^{⊥ α (DY)} (x, b; μ) = \int d^{2} k_{⊥} e^{- i k_{⊥} \cdot b_{⊥}} \frac{k _{⊥}^{α}}{M _{π}} h_{1, q / π}^{⊥ (DY)} (x, k_{⊥}^{2}; μ) .

\tilde{h}_{1, q / π}^{⊥ α (DY)} (x, b; μ) = \int d^{2} k_{⊥} e^{- i k_{⊥} \cdot b_{⊥}} \frac{k _{⊥}^{α}}{M _{π}} h_{1, q / π}^{⊥ (DY)} (x, k_{⊥}^{2}; μ) .

T_{q / π, F}^{(σ)} (x, x; μ) = \int d^{2} k_{⊥} \frac{k _{⊥}^{2}}{M _{π}} h_{1, q / π}^{⊥} (x, k_{⊥}^{2}; μ) = 2 M_{π} h_{1, q / π}^{⊥ (1)} .

T_{q / π, F}^{(σ)} (x, x; μ) = \int d^{2} k_{⊥} \frac{k _{⊥}^{2}}{M _{π}} h_{1, q / π}^{⊥} (x, k_{⊥}^{2}; μ) = 2 M_{π} h_{1, q / π}^{⊥ (1)} .

h_{1, q / π}^{α ⊥} (x, b; Q) = (\frac{- i b _{⊥}^{α}}{2}) e^{- \frac{1}{2} S_{P} (Q, b_{*}) - S_{NP}^{f_{1, q / π}} (Q, b)} T_{q / π, F}^{(σ)} (x, x; μ_{b}) .

h_{1, q / π}^{α ⊥} (x, b; Q) = (\frac{- i b _{⊥}^{α}}{2}) e^{- \frac{1}{2} S_{P} (Q, b_{*}) - S_{NP}^{f_{1, q / π}} (Q, b)} T_{q / π, F}^{(σ)} (x, x; μ_{b}) .

\frac{k _{⊥}}{M _{π}} h_{1, q / π}^{⊥} (x, k_{⊥}; Q) = \int_{0}^{\infty} d b (\frac{b ^{2}}{2 π}) J_{1} (k_{⊥} b) e^{- \frac{1}{2} S_{P} (Q, b_{*}) - S_{NP}^{f_{1, q / π}} (Q, b)} h_{1, q / π}^{⊥ (1)} (x; μ_{b}) .

\frac{k _{⊥}}{M _{π}} h_{1, q / π}^{⊥} (x, k_{⊥}; Q) = \int_{0}^{\infty} d b (\frac{b ^{2}}{2 π}) J_{1} (k_{⊥} b) e^{- \frac{1}{2} S_{P} (Q, b_{*}) - S_{NP}^{f_{1, q / π}} (Q, b)} h_{1, q / π}^{⊥ (1)} (x; μ_{b}) .

π^{-} (P_{π}) + p^{↑} (P_{p}) ⟶ γ^{*} (q) + X ⟶ l^{+} (ℓ) + l^{-} (ℓ^{'}) + X,

π^{-} (P_{π}) + p^{↑} (P_{p}) ⟶ γ^{*} (q) + X ⟶ l^{+} (ℓ) + l^{-} (ℓ^{'}) + X,

s = (P_{π} + P_{p})^{2}, x_{π} = \frac{Q ^{2}}{2 P _{π} \cdot q}, x_{p} = \frac{Q ^{2}}{2 P _{p} \cdot q},

s = (P_{π} + P_{p})^{2}, x_{π} = \frac{Q ^{2}}{2 P _{π} \cdot q}, x_{p} = \frac{Q ^{2}}{2 P _{p} \cdot q},

x_{F} = 2 q_{L} / s = x_{π} - x_{p}, τ = Q^{2} / s = x_{π} x_{p}, y = \frac{1}{2} ln \frac{q ^{+}}{q ^{-}} = \frac{1}{2} ln \frac{x _{π}}{x _{p}},

Peer Reviews

No public reviews on file for this paper yet. If you reviewed it on a platform where reviews are public (OpenReview, ICLR, NeurIPS, ICML), you can paste yours below so the community can read it here.

Videos

No videos yet. Explain this paper in a talk, walkthrough, or lecture? Add one.

Full text

The $\sin(2\phi-\phi_{S})$ azimuthal asymmetry in the pion induced Drell-Yan process within TMD factorization

Hui Li

School of Physics, Southeast University, Nanjing 211189, China

Xiaoyu Wang

[email protected]

School of Physics and Engineering, Zhengzhou University, Zhengzhou, Henan 450001, China

Zhun Lu

[email protected]

School of Physics, Southeast University, Nanjing 211189, China

Abstract

We investigate the single transverse-spin asymmetry with a $\sin(2\phi-\phi_{S})$ modulation in the pion-induced Drell-Yan process within the theoretical framework of the transverse momentum dependent factorization. The asymmetry is contributed by the convolution of the Boer-Mulders function and the transversity. We adopt the model results for the distributions of the pion meson from the light-cone wavefunction approach and the available parametrization for the distributions of the proton to numerically estimate the $\sin(2\phi-\phi_{S})$ asymmetry in $\pi^{-}p$ Drell-Yan at the kinematics of the COMPASS at CERN. To implement the TMD evolution formalism of parton distribution functions, we apply the recently extracted nonperturbative Sudakov form factor associated with the distribution functions of the proton and the pion. It is found that our prediction on the single transverse-spin dependent asymmetry ${\sin(2\phi-\phi_{S})}$ as functions of $x_{p}$ , $x_{\pi}$ , $x_{F}$ and $q_{\perp}$ is qualitatively consistent with the recent COMPASS measurement in both sign and magnitude.

I Introduction

The Boer-Mulders function, denoted by $h_{1}^{\perp}$ , is one of the eight transverse momentum dependent (TMD) parton distribution functions (PDFs) describing the partonic structure of hadrons at leading-twist level. It represents the transversely polarization asymmetry of quarks inside an unpolarized hadron Boer:1997nt ; Boer:1999mm arising from the correlation between the quark spin and the quark transverse momentum, thereby it manifests the novel structure of hadrons. However, the very existence of the Boer-Mulders function was not so obvious. Similar to its chiral-even partner–the Sivers function $f_{1T}^{\perp}$ Sivers:1989cc , the Boer-Mulder function was initially thought to vanish under the constraint of (naive) time reversal invariance of QCD Collins:1992kk . The situation was changed after explicit model calculations Brodsky:2002cx ; Brodsky:2002rv ; Boer:2002ju incorporating gluon exchange between the struck quark and the spectator shows that the T-odd distributions can actually survive. The crucial ingredient in the argument is the Wilson lines (or the gauge links) appearing in the full gange-invariant definition of TMDs Collins:2002kn ; Ji:2002aa . The presence of the Wilson lines also indicates that the T-odd distributions, such as Sivers function and the Boer-Mulders function are process dependent, i.e., they change sign Collins:2002kn ; Brodsky:2002rv ; Boer:2002ju between the semi-inclusive deeply inelastic scattering and Drell-Yan process, a vital prediction which needs verification by future experimental measurement. In the last decades, the Boer-Mulders function of the proton as well as that of the pion has been studied intensively in literature Boer:2002ju ; Gamberg:2003ey ; Yuan:2003wk ; Pobylitsa:2003ty ; Bacchetta:2003rz ; Lu:2004au ; Lu:2006ew ; Gamberg:2007wm ; Burkardt:2007xm ; Bacchetta:2008af ; Zhang:2008nu ; Meissner:2008ay ; Courtoy:2009pc ; Gamberg:2009uk ; Lu:2009ip ; Barone:2009hw ; Barone:2010gk ; Pasquini:2010af ; Lu:2012hh ; Pasquini:2014ppa ; Lu:2016pdp ; Wang:2017onm by models and phenomenological analysis.

As the Boer-Mulders function is a chiral-odd distribution function, it has to be coupled with another chirlal-odd distribution/fragmentation function to survive in a high energy scattering process. A promising process for accessing the Boer-Mulders function is the unpolarized Drell-Yan process, which displays an azimuthal dependence of the final-state dilepton with the $\cos 2\phi$ modulation. As proposed by Boer Boer:1999mm , the coupling of two Boer-Mulder function from each incident hadrons can generate such asymmetry. However, recently studies based on higher order perturbative QCD Peng:2015spa ; Lambertsen:2016wgj ; Chang:2018pvk show that gluon radiation in hard scattering can also give rise to the $\cos 2\phi$ asymmetry substantially, making the extraction of the Boer-Mulders function rather difficult. In the unpolarized semi-inclusive deep inelastic scattering, the combination of the Boer-Mulders function and the Collins fragmentation function $H_{1}^{\perp}$ can lead to a similar $\cos 2\phi_{h}$ azimuthal asymmetry of the final state spin-0 hadron. But this asymmetry is contaminated by the so-call Cahn effectcahn ; Barone:2005kt ; Barone:2008tn , which is a higher-twist kinematical effect due to the transverse motion of the unpolarized quarks. A cleaner process for accessing the Boer-Mulders function is the single transversely polarized Drell-Yan. In this process, the convolution of the Boer-Mulders function and the transversity distribution $h_{1}$ can give rise to a $\sin(2\phi-\phi_{S})$ asymmetry Boer:1999mm ; Arnold:2008kf with $\phi_{S}$ the azimuthal angle of target transverse spin. This makes the transversity function an ideal probe in analysing the information of the Boer-Mulders function from single transversely polarized Drell-Yan because of less contribution from the background. Recently, the first measurement on the $\sin(2\phi-\phi_{S})$ asymmetry has been performed by the COMPASS Aghasyan:2017jop , which adopted a pion beam to collide on the transversely polarized nucleon target. Although no clear tendency is observed on the $\sin(2\phi-\phi_{S})$ asymmetry due to relatively large statistical uncertainties, it indeed indicates negative sign and substantial size.

In this work, we will study the estimate the $\sin(2\phi-\phi_{S})$ asymmetry of the pion-nucleon Drell-Yan process by considering the convolution $h_{1}^{\perp}\otimes h_{1}$ . The main purpose is to investigate the feasibility of accessing the Boer-Mulders function from single polarized Drell-Yan. The theoretical tool we adopt in this study is the TMD factorization Collins:1981uk ; Collins:1984kg ; Collins:2011zzd ; Ji:2004xq which is applicable in the region the transverse momentum of the dilepton $q_{\perp}$ is much smaller than the hard scale $Q$ . The TMD factorization has been widely applied to various high energy processes, such as the semi-inclusive deep inelastic scatering (SIDIS) Collins:1981uk ; Collins:2011zzd ; Ji:2004wu ; Aybat:2011zv ; Collins:2012uy ; Echevarria:2012pw , $e^{+}e^{-}$ annihilation Collins:2011zzd ; Pitonyak:2013dsu ; Boer:2008fr , Drell-Yan Collins:2011zzd ; Arnold:2008kf and W/Z production in hadron collision Collins:2011zzd ; Collins:1984kg ; Lambertsen:2016wgj . The TMD factorization can be also extended to the moderate $q_{\perp}$ region where an equivalence Ji:2006ub ; Ji:2006vf between the TMD factorization and the twist-3 collinear factorization is found. From the perspective of TMD factorization, the physical observables in the region $q_{\perp}\ll Q$ can be expressed as the convolution of the factors related to hard scattering and the well-defined TMD distributions or fragmentation functions (collectively called as TMDs). One of the main features of the TMD formalism is that it provides a systematic approach to deal with the evolution of the TMDs. In this formalism, the energy evolution (or the scale dependence) of TMDs are governed by the so-called Collins-Soper equation Collins:1981uk ; Collins:1984kg ; Collins:2011zzd ; Idilbi:2004vb . The solution of the evolution equation shows the changes of TMDs from a initial scale to another scale may be determined by an exponential form of the Sudakov-like form factor Collins:1984kg ; Collins:2011zzd ; Aybat:2011zv ; Collins:1999dz , which can be separated to the perturbative part and nonperturbative part. The former one is perturbatively calculable, while the later one is usually by phenomenological extraction from experimental data. In this paper, we wil consider the evolution of both the the pion Boer-Mulders function and the proton transversity to estimate the $\sin(2\phi-\phi_{S})$ asymmetry at the kinematics of COMPASS and compare the results with recent COMPASS measurement.

The rest of the paper is organized as follows. In Sec. II, we provide a detailed review on the TMD evolution formalism for the unpolarized and polarized TMDs involved in the calculation. particularly, we will present our choice on the nonperturbative Sudakov form factors associated with the TMDs. In Sec. III, we derive the theoretical expression of ${\sin(2\phi-\phi_{S})}$ asymmetry in the pion-nucleon Drell-Yan within the framework of TMD factorization. In Sec. IV, we estimate the asymmetry at the COMPASS kinematics using a model result of the pion Boer-Mulders function and a parametrization for proton transversity as input. We also provide some discussion based on our numerical result. We summarize the paper in Sec. V.

II The TMD evolution of distribution functions

In this section, we review the evolution formalism of the unpolarized distribution function $f_{1}$ , the Boer-Mulders function $h_{1}^{\perp}$ of the pion as well as the transversity function $h_{1}$ of the proton, within the TMD factorization.

TMD evolution is usually performed in the coordinate $b_{\perp}$ -space, where $b_{\perp}$ is conjugated to $k_{\perp}$ in the transverse momentum space via Fourier transformation Collins:1984kg ; Collins:2011zzd . One of the main advantages of $b_{\perp}$ -space is that the cross section can be expressed as the production of $b_{\perp}$ -dependent functions instead of the complicate convolution of functions in $k_{\perp}$ -space. In the TMD factorization based on different schemes (such as the CS-81 Collins:1981uk , JMY Ji:2004xq ; Ji:2004wu and Collins-11 schemes Collins:2011zzd ), the TMD distribution functions $\tilde{F}(x,b;\mu,\zeta_{F})$ in $b_{\perp}$ space depend on two energy scales. One is the renormalization scale $\mu$ which is related to the corresponding collinear PDFs, the other is the energy scale $\zeta_{F}$ serving as a cutoff to regularize the light-cone singularity in the operator definition of the TMD distributions. The two energy dependencies are encoded in different evolution equations. For the $\zeta_{F}$ dependence of the TMD distributions, it is determined by the Collins-Soper (CS) equation Collins:1981uk ( $b=|\bm{b}_{\perp}|$ ):

[TABLE]

while the $\mu$ dependence is derived from the renormalization group equation as

[TABLE]

with $\tilde{K}$ the CS evolution kernel, and $\gamma_{K}$ and $\gamma_{F}$ the anomalous dimensions. Solving those equations, one can obtain the general solution for the energy dependence of $\tilde{F}$ :

[TABLE]

where $\mathcal{F}$ is the factor related to the hard scattering, $S(Q,b)$ is the Sudakov form factor. Hereafter, we will set $\mu=\sqrt{\zeta_{F}}=Q$ , and express $\tilde{F}(x,b;\mu=Q,\zeta_{F}=Q^{2})$ as $\tilde{F}(x,b;Q)$ for simplicity. Eq. (4) demonstrates that the distribution $\tilde{F}$ at an arbitrary scale $Q$ can be determined by the same distribution at an initial scale $\mu_{i}$ through the evolution encoded by the exponential form $\mathrm{exp}(-S(Q,b))$ .

Although Eq. (4) provides the general structure for the evolution of TMD distributions in $b$ space, it is only possible to calculate the $b$ dependence of $\mathcal{F}$ perturbatively in the small $b$ region. In the large $b$ region, the $b$ -dependence of $\mathcal{F}$ turns to be nonperturbative. A convenient way to take into account the evolution behavior of $\tilde{F}(x,b;Q)$ in the large $b$ region is to include a nonperturbative Sudakov-like form factor $S_{\mathrm{NP}}$ . The latter one is usually given in a parameterized form, which can be obtained by fitting from the experimental data. To allow a smooth transition of $b$ from perturbative region to nonperturbative region as well as to avoid the hitting on the Landau pole, one can set a parameter $b_{\mathrm{max}}$ to be the boundary between the two different regions. The typical value of $b_{\mathrm{max}}$ is chosen around $1\ \mathrm{GeV}^{-1}$ to guarantee that $b_{\ast}$ is always in the perturbative region. A $b$ -dependent function $b_{\ast}(b)$ may be also introduced to have the property $b_{\ast}\approx b$ at small $b$ value and $b_{\ast}\approx b_{\mathrm{max}}$ at large $b$ value. There are several different choices on the form of $b_{\ast}(b)$ Collins:1984kg ; Collins:2016hqq ; Bacchetta:2017gcc . In this work we choose it as $b_{\ast}=b/\sqrt{1+b^{2}/b_{\rm max}^{2}}\ ,~{}b_{\rm max}<1/\Lambda_{\mathrm{QCD}}$ Collins:1984kg ; Kang:2015msa .

Combining the perturbative part and the nonperturbative part, one has the complete result for the Sudakov form factor appearing in Eq. (4)

[TABLE]

with the boundary of the two parts set by the $b_{\mathrm{max}}$ . The perturbative part $S_{\mathrm{P}}(Q,b)$ has been studied Echevarria:2014xaa ; Kang:2011mr ; Aybat:2011ge ; Echevarria:2012pw ; Echevarria:2014rua in details and has the following form:

[TABLE]

which is the same for different kinds of distribution functions, namely, $S_{P}$ is spin-independent. In addition, the coefficients $A$ and $B$ in Eq.(6) can be expanded as the series of $\alpha_{s}/{\pi}$ :

[TABLE]

In this work, we will take $A^{(n)}$ up to $A^{(2)}$ and $B^{(n)}$ up to $B^{(1)}$ in the accuracy of next-to-leading-logarithmic (NLL) order Collins:1984kg ; Landry:2002ix ; Qiu:2000ga ; Kang:2011mr ; Aybat:2011zv ; Echevarria:2012pw :

[TABLE]

A general form of the nonperturbative part of the Sudakov form factor $S_{\rm NP}(Q;b)$ was suggested in Ref. Collins:1984kg :

[TABLE]

Here, $g_{i}(b)$ are the functions of the impact parameter $b$ . Particularly, $g_{2}(b)$ contains the information on the large $b$ behavior of the evolution kernel $\tilde{K}$ , while $g_{1}(b)$ contains information about the intrinsic nonperturbative transverse motion of bound partons, i.e., it depends on the type of the hadron and quark flavor. It might also depend on the momentum fraction of the partons $x$ Su:2014wpa . It is also worth pointing out that $g_{2}(b)$ is universal for different types of TMDs and does not depend on the particular process, which is one of the important predictions of QCD factorization theorems involving TMDs Collins:2011zzd ; Aybat:2011zv ; Echevarria:2014xaa ; Kang:2015msa .

For $S_{\mathrm{NP}}$ associated with the $pp$ collision, a parametrization that can describe the SIDIS and Drell-Yan data with $Q$ values ranging from a few to ten GeV has been proposed in Ref. Su:2014wpa

[TABLE]

The parameters $g_{i}$ are fitted from the nucleon-nucleon Drell-Yan process data Ito:1980ev ; Antreasyan:1981uv ; Moreno:1990sf ; Affolder:1999jh ; Abbott:1999wk ; Abazov:2007ac ; Aaltonen:2012fi at the initial scale of $Q^{2}_{0}=2.4\ \mathrm{GeV}^{2}$ yielding $g_{1}=0.212,\ g_{2}=0.84,\ g_{3}=0$ . Since the nonperturbative form factor $S_{\mathrm{NP}}$ for quarks from the one proton and antiquarks from another proton satisfies Prokudin:2015ysa

[TABLE]

$S_{\mathrm{NP}}$ associated with a single TMD distribution function can be expressed as

[TABLE]

In our calculation of the pion-proton Drell-Yan process, we will adopt the above form factor for the unpolarized TMD distributions of the proton.

For the nonperturbative form factors of the pion distribution function, we adopt the parametrization proposed in Ref. Wang:2017zym

[TABLE]

which has the same form as that for the proton (in the case $g_{3}=0$ ). After fitting to the $\pi^{-}N$ Drell-Yan data Conway:1989fs , the values of the parameters $g_{1}^{\pi}$ and $g_{2}^{\pi}$ are obtained at the initial energy scale $Q^{2}_{0}=2.4\ \mathrm{GeV}^{2}$ as $g^{\pi}_{1}=0.082$ and $g^{\pi}_{2}=0.394$ . In the fit we also chose $b_{\mathrm{max}}=1.5\ \mathrm{GeV}^{-1}$ , in consistence with the choice in Ref .Su:2014wpa . We note that a form of $S^{f_{1,q/\pi}}_{\mathrm{NP}}$ motivated by the NJL model was given in Ref. Ceccopieri:2018nop .

Besides the Sudakov form factor in Eq. (4), another important element in Eq. (4) is the TMD distribution function at a fixed scale $\tilde{F}(x,b,\mu)$ . In the small $b$ region $1/Q\ll b\ll 1/\Lambda$ , the TMD distributions at a fixed scale $\mu$ can be expressed as the convolution of the perturbatively calculable hard coefficients $C$ and the corresponding collinear counterparts, which could be the collinear PDFs or the multiparton correlation functions Collins:1981uk ; Bacchetta:2013pqa

[TABLE]

The convolution $\otimes$ regarding the momentum fraction of $x$ is given by

[TABLE]

and $F_{i/H}(\xi,\mu)$ is the corresponding collinear counterpart of the TMD distribution of flavor $i$ in hadron $H$ at the energy scale $\mu$ , which could be a dynamic scale related to $b_{*}$ by $\mu_{b}=c_{0}/b_{*}$ , with $c_{0}=2e^{-\gamma_{E}}$ and the Euler Constant $\gamma_{E}\approx 0.577$ Collins:1981uk . $\sum_{i}$ is the sum of both quark and antiquark flavors.

It is straightforward to rewrite the scale-dependent TMD distribution function $\tilde{F}$ of the proton and the pion in $b$ space

[TABLE]

The factor of $\frac{1}{2}$ in front of $S_{\mathrm{P}}$ comes from the fact that $S_{\mathrm{P}}$ of quarks and antiquarks satisfies the relation Prokudin:2015ysa

[TABLE]

The hard coefficients $C^{F}_{q\leftarrow i}$ and $\mathcal{F}$ for $f_{1}$ and $h_{1}$ have been calculated up to next-to-leading order (NLO), while those for the Boer-Mulders function still remain in the leading order (LO). For consistency, in this work we will adopt the LO results of the $C$ coefficients for $f_{1}$ , $h_{1}^{\perp}$ and $h_{1}$ . That is, we take $\mathcal{F}=1$ and $C^{F}_{q\leftarrow i}=\delta_{qi}\delta(1-x)$ for $F=f_{1},h_{1}$ and $h_{1}^{\perp}$ . We also note that a calculation in Ref. Kang:2015msa shows that the NLO $C$ -coefficient for $h_{1}$ vanishes.

With all the ingredients above, we can obtain the unpolarized distribution function of the proton and pion in $b$ space as

[TABLE]

The distribution function in the transverse momentum space can be obtained by performing the Fourier transformation on the $\tilde{f}_{1,q/H}(x,b;Q)$

[TABLE]

where $J_{0}$ is the Bessel function of the first kind, and $k_{\perp}=|\bm{k}_{\perp}|$ .

Similar to the unpolarized distribution function, the transversity distribution of the proton in $b$ -space and $k_{\perp}$ space can be obtained as Kang:2015msa

[TABLE]

where the factors and coefficients related to the hard scattering are adopted at LO and the corresponding collinear distribution is the integrated transversity $h_{1}(x)$ . The nonperturbative Sudakov form factor associated with the proton transversity distribution is also assumed to be the same as that for unpolarized distribution function Kang:2015msa .

According to Eq. (17), in the small $b$ region, we can also express the Boer-Mulders function of the pion beam at a fixed energy scale $\mu$ in terms of the perturbatively calculable coefficients and the corresponding collinear correlation function

[TABLE]

where the hard coefficients are calculated up to LO, and the Boer-Mulders function in the $b$ space is defined as

[TABLE]

The collinear function $T^{(\sigma)}_{q/\pi,F}(x,x;\mu)$ is the chiral-odd twist-3 quark-gluon-quark correlation function, which is related to the first transverse moment of the Boer-Mulders function $h_{1,q/\pi}^{\perp(1)}$ by

[TABLE]

As for the nonperturbative part of the Sudakov form factor associated with the Boer-Mulders function, the information still remains unknown. In a practical calculation, we assume that it is the same as $S_{\mathrm{NP}}^{f_{1,q/\pi}}$ , i.e., $S_{\mathrm{NP}}^{h_{1,q/\pi}^{\perp}}=S_{\mathrm{NP}}^{f_{1,q/\pi}}$ . Therefore, we can obtain the Boer-Mulders function of the pion in $b$ -space as

[TABLE]

After performing the Fourier transformation back to the transverse momentum space, one can get the Boer-Mulders function as

[TABLE]

III Formalism of the ${\sin(2\phi-\phi_{S})}$ asymmetry in Drell-Yan process

In this section, we present the formalism of the ${\sin(2\phi-\phi_{S})}$ asymmetry in Drell-Yan process within TMD factorization following the procedure in Ref. Collins:2011zzd . We take into account the TMD evolution effects to obtain the theoretical expression of the ${\sin(2\phi-\phi_{S})}$ asymmetry, which arises from the convolution of Boer-Mulders function of the pion beam and the transversity distribution function of the proton target at the leading twist.

The process we study is the pion-induced Drell-Yan process

[TABLE]

where $P_{\pi}$ , $P_{p}$ and $q$ stand for the four-momenta of the incoming $\pi^{-}$ meson, the proton target and the virtual photon, respectively, $Q^{2}=q^{2}$ is the invariant mass square of the lepton pair, and $\uparrow$ denotes the transverse polarization of the target. We adopt the following kinematical variables to express the experimental observables Collins:1984kg ; Gautheron:2010wva

[TABLE]

where $s$ is the total center-of-mass (c.m.) energy squared; $x_{\pi}$ and $x_{p}$ are the Bjorken variables of the pion and proton, respectively; $q_{L}$ is the longitudinal momentum of the virtual photon in the c.m. frame of the incident hadrons; $x_{F}$ is the Feynman $x$ variable; and $y$ is the rapidity of the lepton pair. Thus, $x_{\pi}$ and $x_{p}$ can be expressed as functions of $x_{F}$ , $\tau$ and of $y$ , $\tau$

[TABLE]

In leading twist, the differential cross section in $\pi p$ Drell-Yan for a transversely polarized target has the following general form Gautheron:2010wva

[TABLE]

Here, $\phi_{S}$ represents the azimuthal angle of the target polarisation vector $S_{T}$ in the target rest frame, $\phi$ and $\theta$ denote the azimuthal and polar angles of the lepton momentum in the Collins-Soper frame, $\hat{\sigma}_{U}=F_{U}^{1}(1+\cos^{2}\theta)$ , with $F_{U}^{1}$ the unpolarized structure function. The symbol $D_{[f(\theta)]}$ denotes the depolarization factor that depends on $\theta$ only, and in LO it is simplified to $\sin^{2}\theta/(1+\cos^{2}\theta)$ . Furthermore, $A_{P}^{f[\phi,\phi_{S}]}$ denotes the azimuthal asymmetry with a modulation of $f[\phi,\phi_{S}]$ , where $P=U$ or $T$ denotes the polarization of the target proton ( $U$ for unpolarized, while $T$ for transversely polarized). The asymmetry $A_{P}^{f[\phi,\phi_{S}]}$ can be written as the ratio between the corresponding structure function $F_{P}^{f[\phi,\phi_{S}]}$ and the unpolarized structure function. In this work, we focus on the ${\sin(2\phi-\phi_{S})}$ asymmetry:

[TABLE]

The denominator can expressed as the convolution of the unpolarized distribution functions from each hadron

[TABLE]

while the numerator ( $\bm{h}=\hat{\bm{q}}\equiv\bm{q}_{\perp}/{|\bm{q}_{\perp}|}$ ) Boer:1999mm ; Arnold:2008kf

[TABLE]

is the convolution of the pion Boer-Mulders distribution and the proton transversity distribution. The convolution of TMDs in the transverse momentum space is defined through the following notation

[TABLE]

with $N_{c}=3$ being the the number of colors, $\bm{q}_{\perp},\bm{k}_{a\perp}$ , and $\bm{k}_{b\perp}$ denoting the transverse momenta of the lepton pair, quark and antiquark in the initial hadrons. Finally, $\omega(\bm{k}_{a\perp},\bm{k}_{b\perp})$ is an arbitrary function of $\bm{k}_{a\perp}$ and $\bm{k}_{b\perp}$ .

In general, it is more convenient to study the structure function first in the $b$ -space, in which the convolution of the TMD distributions can be resolved to the product of $b$ -dependent TMDs. The physical observables can be obtained through a Fourier transformation from the $b$ -space to the $k_{\perp}$ -space. Using the property of the following Fourier transformation

[TABLE]

One can obtain the spin-dependent structure function $F_{T}^{\sin(2\phi-\phi_{S})}$ as

[TABLE]

where, we have used Eqs. (25), (27) and (28). The unpolarized structure function can be expressed in a similar way:

[TABLE]

where the expression of the unpolarized distribution function in Eq. (21) is included and the definition of the unpolarized distribution function in $b$ -space is

[TABLE]

IV Numerical calculation

Using the framework set up above, in this section we present the numerical calculation of the $\sin(2\phi-\phi_{S})$ azimuthal asymmetry in the pion-induced transversely polarized Drell-Yan process. We estimate the asymmetry at the kinematics of the COMPASS Drell-Yan program and compare it with the recent COMPASS measurement. To do this we need to know the corresponding distribution functions of the pion meson, as well as those of the proton target. For the former one, as there is no extraction on the Boer-Mulders function of the pion meson, we apply a model result based on the light-cone wave function of the pion meson from Ref. Wang:2017onm at the model scale $\mu_{0}^{2}=0.25\mathrm{GeV}^{2}$ :

[TABLE]

where the values of the parameters are Xiao:2003wf ; Wang:2017onm ,

[TABLE]

The corresponding collinear twist-3 distribution $T^{\sigma}_{q,F}(x,x,\mu_{0})$ at the model scale can be obtained by using Eq. (28). For consistency, we apply the unpolarized distribution function of the pion meson $f_{1\pi}(x)$ using the same model.

For the collinear distributions of the proton, we resort to existing parametrizations, i.e., we adopt the NLO set of the CT10 parametrization Lai:2010vv (central PDF set) for the unpolarized distribution function $f_{1}(x)$ of the proton, and we choose the transversity distribution extracted from SIDIS data Kang:2015msa via the same TMD evolution formalism:

[TABLE]

where $g_{1}^{q}$ is helicity distribution function deFlorian:2009vb .

We apply the QCDNUM evolution package Botje:2010ay to perform the evolution of $f_{1,q/\pi}$ from the model scale $\mu_{0}$ to another energy. As for the energy evolution of the twist-3 collinear correlation function $T^{(\sigma)}_{q,F}$ , the evolution effect has been studied in Refs. Kang:2012em ; Kang:2008ey ; Vogelsang:2009pj ; Zhou:2008mz ; Braun:2009mi . For simplicity, we only consider the homogenous term in the evolution kernel

[TABLE]

with $\Delta_{T}\,P_{qq}$ being the evolution kernel for the transversity distribution function $h_{1}(x)$ . We customize the original code of QCDNUM to include the approximate kernel in Eq. (46). Similarly, we also include the kernel in Eq. (47) to solve the DGLAP evolution equations for the transversity distribution function of proton.

The COMPASS Collaboration at CERN has reported the first measurement of the transverse-spin-dependent azimuthal asymmetries in the Drell-Yan process Aghasyan:2017jop in which a $\pi^{-}$ beam with $P_{\pi}=\ 190\ \mathrm{GeV}$ collides on a polarized $\mathrm{NH}_{3}$ target Gautheron:2010wva ; Aghasyan:2017jop (which can serve as a transversely polarized nucleon target). The covered kinematical ranges are as follows

[TABLE]

In Fig. 1, we plot our numerical result of the $\sin(2\phi-\phi_{S})$ azimuthal asymmetry as functions of $x_{N},x_{\pi},x_{F}$ and $q_{\perp}$ in the pion-induced Drell-Yan process based on the TMD factorization formalism described in Eqs. (35), (41), and (40) at the kinematics of COMPASS. To make the TMD factorization valid in the kinematic region $q_{\perp}\ll Q$ , the integration over the transverse momentum $q_{\perp}$ is performed in the region of $0<q_{\perp}<2~{}\mathrm{GeV}$ , which is the same as the cut in Ref. Sun:2013hua . The upper panels of Fig. 1 show the asymmetries as functions of $x_{p}$ (left panel) and $x_{\pi}$ (right panel); and the lower panels depict the $x_{F}$ -dependent and $q_{\perp}$ -dependent asymmetries, respectively. In the figure we also show the experimental data measured by the COMPASS Collaboration Aghasyan:2017jop for comparison.

As shown in Fig. 1, in all the cases the $\sin(2\phi-\phi_{S})$ azimuthal asymmetry in the $\pi^{-}p$ Drell-Yan from our calculation is negative, in agreement with most of the data from COMPASS. Our estimate also shows that the asymmetry changes slightly with the change of $x_{N}$ , $x_{\pi}$ and $x_{F}$ , and the magnitude of the $x_{N}$ -, $x_{\pi}$ - and $x_{F}$ -dependent asymmetries is around 0.05 to 0.10. For the $q_{\perp}$ asymmetry, we find that its magnitude is about 0.05 to 0.15 and moderately increases with increasing $q_{\perp}$ in the region $q_{T}<2\mathrm{GeV}$ . Our numerical estimates show that the $A_{UT}^{\sin(2\phi-\phi_{S})}$ is sizable at the the kinematics of COMPASS and is qualitatively consistent with the COMPASS measurement after considered the uncertainties of the data. Our study demonstrate that, with the current knowledge on the distributions of the proton, it is promising to apply the evolution formalism of TMD distributions to study the SSA contributed by the chiral-odd distributions at the kinematics of COMPASS. Our calculation also indicates that the proton transversity distribution may be served as a probe to access the pion Boer-Mulders function as well as the corresponding nonperturbative Sudakov form factor in the context of the current formalism on the transversely polarized $\pi^{-}p$ Drell-Yan process.

V Conclusion

In this work, we applied the TMD factorization to study the $\sin(2\phi-\phi_{S})$ azimuthal asymmetry in the single transversely polarized $\pi^{-}p$ Drell-Yan process that is accessible at COMPASS. The asymmetry arises from the coupling of the Boer-Mulders function of the pion beam and the transversity distribution of the proton target. We took into account the TMD evolution of the asymmetry by including the Sudakov form factor for the TMD distributions of the pion and proton. The hard coefficients associated with the corresponding collinear functions are kept in the leading-order accuracy. For the transversity distribution of the proton used in the study, we employed a recent parametrization for which the same TMD evolution effect is considered. For the distributions of the pion meson, we chose the result from a model calculation incorporating the light-cone wave function approach. As the nonperturbative Sudakov form factor associated with the pion Boer-Mulders function is still unknown, we assume that it is the same as that of the unpolarized distribution function. The latter one has been extracted from the unpolarized $\pi N$ Drell-Yan data. We then calculated the $\sin(2\phi-\phi_{S})$ azimuthal asymmetry in the $\pi^{-}p$ Drell-Yan process at the kinematics of COMPASS. Our analysis demonstrated that, within the framework of TMD evolution, the $\sin(2\phi-\phi_{S})$ asymmetry at COMPASS can be qualitatively described (sign and magnitude) by the current analysis on the TMD distributions of the pion and the proton. Furthermore, our study may provide a framework to access the Boer-Mulders function of the pion and the corresponding nonperturbative Sudakov form factor through transversely polarized $\pi p$ data.

Acknowledgements

This work is partially supported by the NSFC (China) grants 11575043, 11847217 and 11120101004. X. Wang is supported by the China Postdoctoral Science Foundation under Grant No. 2018M640680.

Bibliography85

The reference list from the paper itself. Each links out to its DOI / PubMed record.

1(1) D. Boer and P. J. Mulders, Phys. Rev. D 57 , 5780 (1998).
2(2) D. Boer, Phys. Rev. D 60 , 014012 (1999).
3(3) D. W. Sivers, Phys. Rev. D 41 , 83 (1990).
4(4) J. C. Collins, Nucl. Phys. B 396 , 161 (1993).
5(5) S. J. Brodsky, D. S. Hwang and I. Schmidt, Phys. Lett. B 530 , 99 (2002).
6(6) S. J. Brodsky, D. S. Hwang and I. Schmidt, Nucl. Phys. B 642 , 344 (2002).
7(7) D. Boer, S. J. Brodsky and D. S. Hwang, Phys. Rev. D 67 , 054003 (2003) [hep-ph/0211110].
8(8) J. C. Collins, Phys. Lett. B 536 , 43 (2002) [hep-ph/0204004].

TL;DR

Contribution

Findings

Abstract

Peer Reviews

Videos

The sin⁡(2ϕ−ϕS)\sin(2\phi-\phi_{S})sin(2ϕ−ϕS​) azimuthal asymmetry in the pion induced Drell-Yan process within TMD factorization

Abstract

I Introduction

II The TMD evolution of distribution functions

III Formalism of the sin⁡(2ϕ−ϕS){\sin(2\phi-\phi_{S})}sin(2ϕ−ϕS​) asymmetry in Drell-Yan process

IV Numerical calculation

V Conclusion

Acknowledgements

The $\sin(2\phi-\phi_{S})$ azimuthal asymmetry in the pion induced Drell-Yan process within TMD factorization

III Formalism of the ${\sin(2\phi-\phi_{S})}$ asymmetry in Drell-Yan process