# Next-to-Leading Order Differential Cross Sections forJ/ψ, ψ(2S), and Υ Production in Proton-Proton Collisions at a Fixed-Target Experiment Using the LHC Beams

Authors: Yu Feng; Jian-Xiong Wang
Journal: Advances in High Energy Physics
 (2015)

---

## Abstract

Using nonrelativistic QCD (NRQCD) factorization, we calculate the yields forJ/ψ, ψ(2S), and Υ(1S) hadroproduction at s=72 GeV and 115 GeV including the next-to-leading order QCD corrections. Both these center-of-mass energies correspond to those obtained with 7 TeV and 2.76 TeV nucleon beam impinging a fixed target. We study the cross section integrated in pt as a function of the (center-of-mass) rapidity as well as the pt differential cross section in the central rapidity region. Using different NLO fit results of the NRQCD long-distance matrix elements, we evaluate a theoretical uncertainty which is certainly much larger than the projected experimental uncertainties with the expected 20 fb−1 to be collected per year with AFTER@LHC for pp collision at the center of mass energy s≃115 GeV.


---

## Body

## 1. Introduction

Nonrelativistic quantum chromodynamics (NRQCD) [1] is the most systematic factorization scheme to describe the decay and production of heavy quarkonia. It allows one to organize the theoretical calculations as double expansions in both the coupling constant αs and the heavy-quark relative velocity v. In the past few years, significant progress has been made in next-to-leading order (NLO) QCD calculations based on NRQCD. Calculations and fits of NRQCD long-distance matrix elements (LDMEs) for both the J/ψ yield and polarization in hadroproduction have been carried out [2–6] as well as for Υ hadroproduction [7, 8]. Using these LMDEs, one can in principle predict the transverse momentum pt differential cross section at any energies. In addition, in a recent study [9], we have discussed the implication of these fits on the energy dependence of the cross sections integrated in pt.In this paper, we predict these differential cross sections for the kinematics of a fixed-target experiment using the LHC beams (AFTER@LHC) [10]. In practice, 7 TeV protons on targets yield to a c.m.s energy close to 115 GeV and 72 GeV for 2.76 TeV nucleons (as in the case of a Pb beam). This corresponds to a range very seldom explored so far, significantly higher than that at CERN-SPS and not far from BNL-RHIC. With the typical luminosity of the fixed-target mode, which allows for yearly luminosities as large as 20 fb−1 in pp collision at s≃ 115 GeV, AFTER@LHC is expected to be a quarkonium and heavy-flavor observatory [10, 11]. In general, the opportunities of a fixed-target experiment using the LHC beam for spin and heavy-ion physics are discussed in [10, 12–14]. With the calculation at s=72 GeV, which is supposed to be a baseline rate where nuclear effects would be added, we confirm that charmonium yields can easily reach 109 per year and 106 for bottomonium at s≃ 115 GeV.


## 2. Next-to-Leading Order Calculation

Following the NRQCD factorization formalism [1], the cross section for quarkonium hadroproduction H can be expressed as (1)dσpp⟶H+X=∑i,j,n∫dx1dx2GpiGpjdσ^ij⟶QQ¯nXOnH,where p is either a proton or an antiproton, Gpi(j) is the parton distribution function (PDF) of p, the indices i,j run over all possible partonic species, and n denotes the color, spin, and angular momentum states of the intermediate QQ¯ pair. For ψ and Υ, namely, the 3S1 quarkonium sates, their leading CO states of relative order O(v4) are 1S0[8], 3S1[8], and 3PJ[8]. Along with the CS transition 3S1[1], we call the total CS + CO contributions as direct production. The short-distance coefficient (SDC) dσ^ will be calculated perturbatively, while the long-distance matrix elements (LDMEs) 〈OnH〉 are governed by nonperturbative QCD effects.Now let us take a look at the parton level processes relevant to this work. As it is well known, the CO contributions to hadroproduction appear atαs2 [15] and their Born contributions are (2)q+q¯⟶QQ¯S318,g+g⟶QQ¯S108,P3J=0,28,where q(q¯) denotes the light quarks (antiquarks).Up toαs3, QCD corrections include real and virtual corrections. One inevitably encounters ultra-violet (UV), infrared (IR), and Coulomb divergences when dealing with the virtual corrections. UV divergences from self-energy and triangle diagrams are canceled upon the renormalization procedure. For the real emission corrections, three kinds of processes should be considered: (3)g+g⟶QQ¯S311,S108,S318,P3J=0,28+g,g+qq¯⟶QQ¯S108,S318,P3J=0,28+qq¯,q+q¯⟶QQ¯S108,S318,P3J=0,1,28+g,some of which involve IR singularities in phase space integration and we adopt the two-cutoff phase space slicing method [16] to isolate these singularities by introducing two small cutoffs, δs and δc. For technical details, we refer readers to [17, 18].One has to note that in (3), the S31[1] production in gg fusion is not really correction. Strictly speaking, it is only the Born-order contribution for hadroproduction with a jet. In fact, all the real emission processes in (3) will be taken as Born-order contributions of quarkonium-jet production.As regards to thept dependent differential cross section, and the QCD NLO corrections in this case are up to αs4, which involves the real emission processes (4)g+g⟶QQ¯n+g+g,g+g⟶QQ¯n+q+q¯,g+qq¯⟶QQ¯n+g+qq¯,q+q¯⟶QQ¯n+g+g,q+q¯⟶QQ¯n+q+q¯,q+q¯⟶QQ¯n+q′+q¯′,q+q⟶QQ¯n+q+q,q+q′⟶QQ¯n+q+q′,where q, q′ denote light quarks with different flavors and (QQ¯)n can be either S31[1], S10[8], S31[8], or P3J[8]. One can find the details of such computations at this order in [18, 19] and some examples in [2, 3, 6–8].All of these calculations are made with the newly updated Feynman Diagram Calculation package [20].


## 3. Constrains on the LDMEs

The color-singlet (CS) LDMEs are estimated from the wave functions at the origin by〈OH(S31[1])〉=3Nc/2π|RH(0)|2, where the wave functions are obtained via potential model calculation [21]. This gives |RJ/ψ(0)|2=0.81GeV3, |Rψ(2S)(0)|2=0.53GeV3, and |RΥ(1S)(0)|2=6.5 GeV3. In the following, we will refer to this contribution as the CSM results when performed separately.The color-octet (CO) LDMEs can only be extracted from data. As for now, SDC are known up to NLO accuracy and the fits of LDMEs can be thus performed at NLO. However, different results are obtained when different dataset is used. We made a selection of these fits in order to assess the theoretical uncertainty induced by the LDMEs. We briefly discuss below these different fit results.In theJ/ψ case, seven groups of LDMEs [2, 5, 6, 22–25] are collected in Table 1. They are extracted by fitting the data of hadroproduction yield [2] or combined with polarization [5, 6] on pp collisions. The first one [22] was based on a wider set of data including ep and γγ system with pt>1 GeV. In [5, 6], the data with pt<7 GeV are excluded in their fit. The fit in [23, 24] took the ηc measurement (pt≥6 GeV) into consideration. Only one of them is used [24] since their results are almost the same. The last one incorporates the leading-power fragmentation corrections together with the QCD NLO corrections, which results in a different SDC and may result in different LDMEs. In [2], Ma et al. fit the data with pt>7 GeV by two linear combinations of LDMEs: (5)M0,r0J/ψ=OJ/ψS108+r0mc2OJ/ψP308,M1,r1J/ψ=OJ/ψS318+r1mc2OJ/ψP308,from which we extract the value of LDMEs by restricting 〈OJ/ψ(S10[8])〉 and 〈OJ/ψ(S31[8])〉 to be positive to get a loose constraint on the 〈OJ/ψ(P30[8])〉 range, from which we choose the center value in order to obtain the three LDMEs (Ma et al. (2011) in Table 1).Table 1
The values of LDMEs forJ/ψ hadroproduction (in units of GeV3).
References
O
J
/
ψ
(
3
S
1
[
1
]
)
O
J
/
ψ
(
1
S
0
[
8
]
)
O
J
/
ψ
(
3
S
1
[
8
]
)
O
J
/
ψ
(
3
P
0
[
8
]
)
/
m
Q
2
Butenschoen and Kniehl (2011) [22]
1.32
3.0
×
1
0
-
2
1.7
×
1
0
-
3
-
4.0
×
1
0
-
3
Chao et al. (2012) [5]
1.16
8.9
×
1
0
-
2
3.0
×
1
0
-
3
5.6
×
1
0
-
3
Ma et al. (2011) [2]
1.16
3.9
×
1
0
-
2
5.6
×
1
0
-
3
8.9
×
1
0
-
3
Gong et al. (2013) [6]
1.16
9.7
×
1
0
-
2
-
4.6
×
1
0
-
3
-
9.5
×
1
0
-
3
Zhang et al. (2014) [23]
0.24~0.90
(0.4~1.1)×10-2
1.0
×
1
0
-
2
1.7
×
1
0
-
2
Han et al. (2015) [24]
1.16
0.7
×
1
0
-
2
1.0
×
1
0
-
2
1.7
×
1
0
-
2
Bodwin et al. (2014) [25]
0
9.9
×
1
0
-
2
1.1
×
1
0
-
2
4.9
×
1
0
-
3As regards theψ(2S), only two NLO analyses have been done (see [2, 6]), both of which excluded the data with pt<7 GeV in their fit. To extract the LDME values from the fit results of Ma et al., the same method is used as for the J/ψ. For Υ(1S), we use three groups of LDMEs [8, 26, 27]. Both of them have exactly accounted for the direct production and the feed-down contributions. In the fit of [26], only the data in pt>15 GeV region are used, while in [8, 27] the region is pt>8 GeV. They all describe the high pt yield data at Tevatron and LHC very well. We gather the LDMEs of ψ(2S) and Υ(1S) in Table 2.Table 2
The values of LDMEs forψ(2S) and Υ(1S) hadroproduction (in units of GeV3).
H
References
O
H
(
3
S
1
[
1
]
)
O
H
(
1
S
0
[
8
]
)
O
H
(
3
S
1
[
8
]
)
O
H
(
3
P
0
[
8
]
)
/
m
Q
2
ψ
(
2
S
)
Gong et al. (2013) [6]
0.76
-
1.2
×
1
0
-
4
3.4
×
1
0
-
3
4.2
×
1
0
-
3
Ma et al. (2011) [2]
0.76
1.4
×
1
0
-
2
2.0
×
1
0
-
3
1.6
×
1
0
-
3
Υ
(
1
S
)
Gong et al. (2014) [8]
9.28
11.2
×
1
0
-
2
-
4.1
×
1
0
-
3
-
6.7
×
1
0
-
3
Han et al. (2014) [26]
9.28
3.5
×
1
0
-
3
5.8
×
1
0
-
2
3.6
×
1
0
-
2
Feng et al. (2015) [27]
9.28
13.6
×
1
0
-
2
6.1
×
1
0
-
3
-
9.3
×
1
0
-
3


## 4. Numerical Results

The differential cross sections as a function of (center-of-mass) rapidity and transverse momentum are considered in this section. In both cases, the CTEQ6M parton distribution functions [28] and the corresponding two-loop QCD coupling constants αs are used. The charm quark mass is set to be mc=1.5 GeV, while the bottom quark mass is taken to be mb=4.75 GeV. The renormalization and factorization scales are chosen as μr=μf=2mQ for rapidity distribution plots, while for the plots of transverse momentum distribution they are μr=μf=μT, with μT=(2mQ)2+pt2. The NRQCD scale is chosen as μΛ=mQ. It is important to note that different choices of these scales may be adopted which can bring some additional uncertainties in our predictions. We show these uncertainties from the scales combined to those from the quark masses for the plot of the rapidity distribution of the yield (and only for the CO contributions). The scale dependence is estimated by varying μr, μf, by a factor of 1/2 and 2 with respect to their central values and quark masses varying 0.1 GeV up and down for J/ψ, as well as 0.25 GeV for Υ. Branching ratios are taken from PDG [29], which give B[J/ψ→μμ] = 0.0596, B[ψ(2S)→μμ] = 0.0079, and B[Υ(1S)→μμ] = 0.0248, respectively. The two phase space cutoffs δs=103 and δc=δs/50 are chosen and the insensitivity of the results on different choices for these cutoffs has been checked.


### 4.1.dσ/dy up to αs3

First, we study thept-integrated cross section (where the whole pt region is integrated) as a function of rapidity. The QCD NLO corrections are up to αs3 in this case. In Figures 1 and 2, we show the rapidity distribution of direct J/ψ, ψ(2S), and Υ(1S) production cross section at center of mass energy s=72 GeV and 115 GeV, respectively. We first discuss the “unnormalized” CO channel contributions shown in Figure 1, where the CO LDMEs are set to unity for all three production channels. For ψ(2S), the CSM is different from J/ψ only by a factor, we therefore do not perform it separately. Obviously, the CSM results (red lines) for both J/ψ and Υ(1S) are small compared with the CO channels. The dominant CO channel for J/ψ is P3J[8] transition, while for Υ(1S) it is S10[8]. Besides, these “unnormalized” contributions for J/ψ display a clear hierarchy, but for Υ(1S), little difference between S31[8] and P3J[8] contributions shows up.Figure 1
The “unnormalized” CO channel contributions for directJ/ψ (a) and Υ(1S) (b) hadroproduction at the c.m.s energy 72 GeV (dot lines) and 115 GeV (dashed lines), respectively. The CO LDMEs for all the channels are set to unity.
(a)
(b)Figure 2
Rapidity distribution of differential cross section for directJ/ψ (a), ψ(2S) (b), and Υ(1S) (c) hadroproduction at the center of mass energy s=72 GeV and s=115 GeV, respectively. The lines are the uncertainty from LDMEs values while the color areas are scales and masses uncertainties.
(a)
(b)
(c)Adopting the LDMEs in Tables1 and 2, we present the rapidity distribution of cross section for various cases in Figure 2. The lines are the central values with different groups of LDMEs, while the colored areas are the uncertainties from scales and quark masses. Only the boundary lines are shown with scales and mass uncertainties. For the J/ψ, six groups of NRQCD results are shown as a band, the boundaries of which has a distance within factor 10. The values of the cross sections are roughly in the region of 104~105 pb. The CSM results are systematically below the full NRQCD band, again by a factor 10. Without a surprise, the CSM seems to be negligible for total NRQCD results. However, some words of cautious are in order here. In fact, as we have discussed in [9], the LO CSM contribution explains the data very well from the RHIC to LHC energies, while the CO LDMEs extracted from pt-differential NLO correction would lead to pt-integrated cross sections overshooting the data. Only the fit from Butenschoen and Kniehl [22], including rather low pt data, provides an acceptable description of the pt-integrated cross section. In other words, most of the predictions in Figure 2 may overshoot the data. One should indeed stress that most of the fits which we used are based on large pt data, while this rapidity distribution bears on the whole pt region. We suppose the band from the LDMEs of Butenschoen and Kniehl [22]; namely, the lowest band (red dashed line) would probably give the best prediction of the J/ψ yields, although their LDMEs clearly face difficulty to describe the polarization data. To be complete, let us mention that, following an observations also done in [9], the CSM yield may underestimate the measurements below RHIC energy.As regards theψ(2S), two groups of LDMEs lead to a consistent predictions which give a cross section around 103 pb at both s=72 GeV and 115 GeV. Accounting for the uncertainties of the scales and quark masses, the cross sections reach 104 pb in the central rapidity region. Nevertheless, these results overestimated the data as discussed in [9].In theΥ(1S) case, two curves are close and a third one is slightly different. Yet, their difference is only in pb units. At RHIC energies and below [9], this reproduces quite well the Υ(1S) data. This should thus also be the case at the energies considered here.


### 4.2.dσ/dpt up to αs4

Now let us discuss the differential cross sections in the transverse momentumpt. In Figure 3, the pt distributions of direct J/ψ, ψ(2S), and Υ(1S) hadroproduction are presented. For J/ψ and ψ(2S), the yields are dominated by the CO contributions, which is larger than that of the CSM at NLO by at least one order of magnitude. The various groups of LDMEs predict J/ψ and ψ(2S) differential yields which are much less spread than for the pt integrated yields; this is expected since the fits are based on a similar distribution but at different energies. Only the one from [25] (the light blue dot-dashed line) seems to depart from the other ones, being from 2 to 10 times larger in J/ψ case. This may be understood by the fact that the fits in [25] have a different SDC compared with others, which would be the source of the difference.Figure 3
Transverse momentum distribution of differential cross section with the rapidityy=0 for direct J/ψ, ψ(2S), and Υ(1S) hadroproduction from (a) to (c), respectively.
(a)
(b)
(c)ForΥ(1S), the red dashed and blue dot-dashed lines are almost parallel with a tiny difference, while the green dot line is obviously lower at low pt region and crosses the other ones as pt increases. This explains the behavior of dσ/dy in Figure 2 with a visible difference between the green curve and the other two.


## 4.1.dσ/dy up to αs3

First, we study thept-integrated cross section (where the whole pt region is integrated) as a function of rapidity. The QCD NLO corrections are up to αs3 in this case. In Figures 1 and 2, we show the rapidity distribution of direct J/ψ, ψ(2S), and Υ(1S) production cross section at center of mass energy s=72 GeV and 115 GeV, respectively. We first discuss the “unnormalized” CO channel contributions shown in Figure 1, where the CO LDMEs are set to unity for all three production channels. For ψ(2S), the CSM is different from J/ψ only by a factor, we therefore do not perform it separately. Obviously, the CSM results (red lines) for both J/ψ and Υ(1S) are small compared with the CO channels. The dominant CO channel for J/ψ is P3J[8] transition, while for Υ(1S) it is S10[8]. Besides, these “unnormalized” contributions for J/ψ display a clear hierarchy, but for Υ(1S), little difference between S31[8] and P3J[8] contributions shows up.Figure 1
The “unnormalized” CO channel contributions for directJ/ψ (a) and Υ(1S) (b) hadroproduction at the c.m.s energy 72 GeV (dot lines) and 115 GeV (dashed lines), respectively. The CO LDMEs for all the channels are set to unity.
(a)
(b)Figure 2
Rapidity distribution of differential cross section for directJ/ψ (a), ψ(2S) (b), and Υ(1S) (c) hadroproduction at the center of mass energy s=72 GeV and s=115 GeV, respectively. The lines are the uncertainty from LDMEs values while the color areas are scales and masses uncertainties.
(a)
(b)
(c)Adopting the LDMEs in Tables1 and 2, we present the rapidity distribution of cross section for various cases in Figure 2. The lines are the central values with different groups of LDMEs, while the colored areas are the uncertainties from scales and quark masses. Only the boundary lines are shown with scales and mass uncertainties. For the J/ψ, six groups of NRQCD results are shown as a band, the boundaries of which has a distance within factor 10. The values of the cross sections are roughly in the region of 104~105 pb. The CSM results are systematically below the full NRQCD band, again by a factor 10. Without a surprise, the CSM seems to be negligible for total NRQCD results. However, some words of cautious are in order here. In fact, as we have discussed in [9], the LO CSM contribution explains the data very well from the RHIC to LHC energies, while the CO LDMEs extracted from pt-differential NLO correction would lead to pt-integrated cross sections overshooting the data. Only the fit from Butenschoen and Kniehl [22], including rather low pt data, provides an acceptable description of the pt-integrated cross section. In other words, most of the predictions in Figure 2 may overshoot the data. One should indeed stress that most of the fits which we used are based on large pt data, while this rapidity distribution bears on the whole pt region. We suppose the band from the LDMEs of Butenschoen and Kniehl [22]; namely, the lowest band (red dashed line) would probably give the best prediction of the J/ψ yields, although their LDMEs clearly face difficulty to describe the polarization data. To be complete, let us mention that, following an observations also done in [9], the CSM yield may underestimate the measurements below RHIC energy.As regards theψ(2S), two groups of LDMEs lead to a consistent predictions which give a cross section around 103 pb at both s=72 GeV and 115 GeV. Accounting for the uncertainties of the scales and quark masses, the cross sections reach 104 pb in the central rapidity region. Nevertheless, these results overestimated the data as discussed in [9].In theΥ(1S) case, two curves are close and a third one is slightly different. Yet, their difference is only in pb units. At RHIC energies and below [9], this reproduces quite well the Υ(1S) data. This should thus also be the case at the energies considered here.


## 4.2.dσ/dpt up to αs4

Now let us discuss the differential cross sections in the transverse momentumpt. In Figure 3, the pt distributions of direct J/ψ, ψ(2S), and Υ(1S) hadroproduction are presented. For J/ψ and ψ(2S), the yields are dominated by the CO contributions, which is larger than that of the CSM at NLO by at least one order of magnitude. The various groups of LDMEs predict J/ψ and ψ(2S) differential yields which are much less spread than for the pt integrated yields; this is expected since the fits are based on a similar distribution but at different energies. Only the one from [25] (the light blue dot-dashed line) seems to depart from the other ones, being from 2 to 10 times larger in J/ψ case. This may be understood by the fact that the fits in [25] have a different SDC compared with others, which would be the source of the difference.Figure 3
Transverse momentum distribution of differential cross section with the rapidityy=0 for direct J/ψ, ψ(2S), and Υ(1S) hadroproduction from (a) to (c), respectively.
(a)
(b)
(c)ForΥ(1S), the red dashed and blue dot-dashed lines are almost parallel with a tiny difference, while the green dot line is obviously lower at low pt region and crosses the other ones as pt increases. This explains the behavior of dσ/dy in Figure 2 with a visible difference between the green curve and the other two.


## 5. Summary

We evaluated the NLO QCD corrections for the directJ/ψ, ψ(2S), and Υ(1S) production at fixed-target LHC energies. We studied the cross section integrated in pt as a function of the rapidity as well as the pt differential cross section in the central rapidity region, including QCD corrections up to αs3 and αs4 contributions, respectively. To perform a reliable prediction, various sets of NRQCD long-distance matrix elements obtained from different fitting methods are considered as well as the uncertainties from the scales and quark masses. With the typical luminosity of the fixed-target mode, which allows for yearly luminosities as large as 20 fb−1 with 7 TeV proton beams, our predictions confirm that charmonium yields can easily reach 109 per year and 106 for bottomonium at the center of mass energy s≃115 GeV in pp collision.



---

*Source: 726393-2015-10-05.xml*