# OnJ/ψ and Υ Transverse Momentum Distributions in High Energy Collisions

Authors: Bao-Chun Li; Ting Bai; Yuan-Yuan Guo; Fu-Hu Liu
Journal: Advances in High Energy Physics
 (2017)

---

## Abstract

The transverse momentum distributions of final-state particles are very important for high energy collision physics. In this work, we investigateJ/ψ and Υ meson distributions in the framework of a particle-production source, where Tsallis statistics are consistently incorporated. The results are in good agreement with the experimental data in p-p and p-Pb collisions at LHC energies. The temperature of the emission source and the nonequilibrium degree of the collision system are extracted.


---

## Body

## 1. Introduction

The investigation of nuclear matter at high energy densities is the main purpose of Relativistic Heavy Ion Collider (RHIC) and Large Hadron Collider (LHC) [1–7]. As a new matter state, quark-gluon plasma (QGP) is a thermalized system, which consists of strongly coupled quarks and gluons in a limited region. The suppression of J/ψ meson with respect to proton-proton (p-p) collisions is regarded as a distinctive signature of the QGP formation and brings valuable insight into properties of the nuclear matter. In proton-nucleus (p-A) collisions, the prompt J/ψ meson suppression has also been observed at large rapidity [8]. The low-viscosity QGP may be created in the small system, 3He + Au collisions [9]. The heavy quarkonium production can be suppressed by the suppression cold-nuclear-matter (CNM) effects, such as nuclear absorption, nuclear shadowing (antishadowing), and parton energy loss.The transverse momentumpT spectra of identified particles produced in the collisions are a vital research for physicists. Now, different models have been suggested to describe the pT distributions of the final-state particles in high energy collisions [10–13], such as Boltzmann distribution, Rayleigh distribution, Erlang distribution, the multisource thermal model, and Tsallis statistics. Different phenomenological models of initial coherent multiple interactions and particle transport have been proposed to discuss the particle production in high energy collisions. Tsallis statistics can deal with nonequilibrated complex systems in condensed-matter research [14]. It is developed to describe the particle production in recent years [15–24].In our previous work [12], the temperature information was understood indirectly by an excitation degree. We have obtained the emission source location dependence of the exciting degree specifically. In this paper, the temperature of the emission source is given directly by combining a picture of the particle-production source with Tsallis statistics. We discuss the transverse momentum distributions of J/ψ in p-p collisions at sNN=7 TeV, sNN=8 TeV, and sNN=13 TeV and p-Pb collisions at sNN=5 TeV. And the Υ distributions in p-p collisions at sNN=8 TeV are also taken into account for comparison.


## 2. Tsallis Statistics in an Emission Source

According to the multisource thermal model [12] and the nuclear geometry picture, at the initial stage of the collision, two cylinder-shaped groups are formed along the beam direction. In the laboratory reference frame, it is assumed that the projectile cylinder is at the positive y space and the target cylinder is at the negative y space. The cylinders are not a real shape and are understood to be a rapidity range of the emission source. The projectile and target cylinders can be regarded as an emission source with a rapidity width. The observed particles are emitted from the emission source.With Tsallis statistics’ success in dealing with nonequilibrated complex systems in condensed-matter research [14], it has been used to understand the particle production in high energy physics. In order to describe the transverse momentum spectra in high energy collisions, several versions of Tsallis distribution are proposed [15–24]. Recently, an improved form of the Tsallis distribution was proposed [18–20] and can meet the thermodynamic consistency. The meson number [16, 17] is given by(1)N=gV∫d3p2π31+q-1E-μT1/1-q,where g, V, p, E, and μ are the degeneracy factor, the volume, the momentum, the energy, and the chemical potential, respectively. The parameter T is the temperature of the emission source and the parameter q is the nonequilibrium degree. Generally, q is greater than 1 and is close to 1. The corresponding momentum distribution is given by(2)Ed3Nd3p=gVE12π31+q-1E-μT1/1-q.Then, the transverse momentum distribution is(3)d2NpTdydpT=gVmTcosh⁡y2π21+q-1mTcosh⁡y-μT1/1-q.When the chemical potential is neglected, at midrapidity y=0, the pT distribution is(4)d2NpTdydpT=gVmT2π21+q-1mTT1/1-q.It is worth noting that the distribution function of pT is only the distribution of mesons emitted from an emission point at y=0 in the emission source, not the final-state distribution due to the nonzero rapidity width of the emission source. By summing the contributions of all emission points, the transverse momentum pT distribution is rewritten as (5)dNpTdpT=c∫yminymaxcosh⁡ydymT1+q-1mTcoshyT1/1-q,where c=gV/2π2 is a normalize constant and ymax and ymin are the maximum and minimum values of the observed rapidity.


## 3. Transverse Momentum Spectra and Discussions

Figure1 shows the double-differential cross section d2σ/dpTdy of J/ψ mesons in p-p collisions at sNN=7 TeV. Figures 1(a), 1(b), 1(c), and 1(d) present prompt J/ψ with no polarisation, J/ψ from b with no polarisation, prompt J/ψ with full transverse polarisation, and prompt J/ψ with full longitudinal polarisation, respectively. The experimental data in y bins are from the LHCb Collaboration [25]. The solid lines indicate our model results, which are in good agreement with the experimental data in all rapidity ranges. The parameters T and q taken for the calculation are listed in Table 1. In different rapidity ranges, the values of the temperature T are different and decrease with increasing the rapidity bins in all four figures of Figure 1. The closer the emission source is to the center y=0, the larger the excitation degree is. The values of q do not change regularly with the rapidity bins.Table 1
Values ofT and q used in Figures 1-2.
Figure
y bins
T (GeV)
q
χ
2
/
d
o
f
Figure1(a)
2.0
<
y
<
2.5
0.845
1.018
0.926
2.5
<
y
<
3.0
0.833
1.020
1.052
3.0
<
y
<
3.5
0.801
1.024
1.078
3.5
<
y
<
4.0
0.789
1.021
1.127
4.0
<
y
<
4.5
0.755
1.015
1.074
Figure1(b)
2.0
<
y
<
2.5
0.863
1.034
1.244
2.5
<
y
<
3.0
0.850
1.033
1.310
3.0
<
y
<
3.5
0.838
1.032
1.317
3.5
<
y
<
4.0
0.783
1.028
1.293
4.0
<
y
<
4.5
0.725
1.026
1.280
Figure1(c)
2.0
<
y
<
2.5
0.950
1.002
1.302
2.5
<
y
<
3.0
0.900
1.011
1.208
3.0
<
y
<
3.5
0.875
1.012
1.226
3.5
<
y
<
4.0
0.867
1.003
1.251
4.0
<
y
<
4.5
0.826
1.001
1.295
Figure1(d)
2.0
<
y
<
2.5
0.845
1.020
1.150
2.5
<
y
<
3.0
0.840
1.017
1.174
3.0
<
y
<
3.5
0.836
1.016
1.127
3.5
<
y
<
4.0
0.829
1.010
1.203
4.0
<
y
<
4.5
0.802
1.001
1.192
Figure2(a)
2.0
<
y
<
2.5
0.882
1.013
1.239
2.5
<
y
<
3.0
0.880
1.013
1.282
3.0
<
y
<
3.5
0.876
1.011
1.255
3.5
<
y
<
4.0
0.873
1.008
1.291
4.0
<
y
<
4.5
0.866
1.003
1.385
Figure2(b)
2.0
<
y
<
2.5
0.885
1.038
1.103
2.5
<
y
<
3.0
0.880
1.036
1.055
3.0
<
y
<
3.5
0.874
1.031
1.346
3.5
<
y
<
4.0
0.830
1.030
1.174
4.0
<
y
<
4.5
0.805
1.029
1.230Figure 1
Double-differential cross sectiond2σ/dpTdy as a function of pT in bins of y for (a) prompt J/ψ with no polarisation, (b) J/ψ from b with no polarisation, (c) prompt J/ψ with full transverse polarisation, and (d) prompt J/ψ with full longitudinal polarisation in p-p collisions at sNN=7 TeV. The symbols indicate the experimental data of the LHCb Collaboration [15]. The solid lines indicate our theoretical calculations.
(a)
(b)
(c)
(d)Figure2 presents the double-differential cross section d2σ/dpTdy of J/ψ mesons in p-p collisions at sNN=8 TeV. Figures 2(a) and 2(b) are prompt J/ψ and J/ψ from b, respectively. Figure 3 presents the double-differential cross section BiSd2σiS/dpTdy of Υ mesons in the same collision, where i= 1, 2, and 3 correspond to Υ1S, Υ2S, and Υ3S mesons, respectively, and B is the dimuon branching fraction. The experimental data are from [26]. The symbols and lines represent the same meanings as those in Figure 1. The model results are also in agreement with the experiment data. The parameters T and q for the calculation are listed in Tables 1 and 2. The freeze-out temperature T also decreases with increasing the rapidity bins for both J/ψ and Υ. In the same y ranges, the T values are smaller for Υ than that for J/ψ. For J/ψ, the T values are generally larger than that at sNN=7 TeV. q for J/ψ and for Υ are 1.003–1.038 and 1.170–1.127, respectively. At pT>12 GeV, the agreement is not very good in Figures 1-2 and the χ2/dof values are in the range 0.926–1.385.Table 2
Values ofT and q used in Figures 3–5.
Figure
y bins
T
q
χ
2
/
d
o
f
Figure3(a)
2.0
<
y
<
2.5
0.744
1.138
0.152
2.5
<
y
<
3.0
0.734
1.132
0.124
3.0
<
y
<
3.5
0.728
1.130
0.160
3.5
<
y
<
4.0
0.670
1.132
0.128
4.0
<
y
<
4.5
0.610
1.127
0.216
Figure3(b)
2.0
<
y
<
2.5
0.795
1.170
0.280
2.5
<
y
<
3.0
0.790
1.168
0.356
3.0
<
y
<
3.5
0.788
1.165
0.209
3.5
<
y
<
4.0
0.784
1.151
0.375
4.0
<
y
<
4.5
0.770
1.133
0.481
Figure3(c)
2.0
<
y
<
2.5
0.703
1.169
0.415
2.5
<
y
<
3.0
0.688
1.168
0.225
3.0
<
y
<
3.5
0.684
1.167
0.233
3.5
<
y
<
4.0
0.680
1.157
0.410
4.0
<
y
<
4.5
0.605
1.131
0.504
Figure4(a)
2.0
<
y
<
2.5
0.962
1.035
0.429
2.5
<
y
<
3.0
0.940
1.034
0.418
3.0
<
y
<
3.5
0.928
1.033
0.371
3.5
<
y
<
4.0
0.902
1.032
0.273
4.0
<
y
<
4.5
0.882
1.031
0.350
Figure4(b)
2.0
<
y
<
2.5
0.988
1.065
0.182
2.5
<
y
<
3.0
0.964
1.063
0.254
3.0
<
y
<
3.5
0.936
1.060
0.261
3.5
<
y
<
4.0
0.910
1.057
0.270
4.0
<
y
<
4.5
0.856
1.054
0.315
Figure5(a)
1.5
<
y
<
2.0
0.792
1.040
0.436
2.0
<
y
<
2.5
0.786
1.037
0.304
2.5
<
y
<
3.0
0.784
1.036
0.391
3.0
<
y
<
3.5
0.782
1.033
0.452
3.5
<
y
<
4.0
0.780
1.025
0.526
Figure5(b)
1.5
<
y
<
2.0
0.843
1.053
0.170
2.0
<
y
<
2.5
0.840
1.050
0.225
2.5
<
y
<
3.0
0.838
1.048
0.303
3.0
<
y
<
3.5
0.836
1.041
0.518
3.5
<
y
<
4.0
0.822
1.026
0.531Figure 2
Double-differential cross sectiond2σ/dpTdy as a function of pT in bins of y for (a) prompt J/ψ and (b) J/ψ from b in p-p collisions at sNN= 8 TeV. The symbols indicate the experimental data of the LHCb Collaboration [25]. The solid lines indicate our theoretical calculations.
(a)
(b)Figure 3
Double-differential cross sectiond2σ/dpTdy as a function of pT in bins of y for (a) Υ1S, (b) Υ2S, and (c) Υ3S mesons in p-p collision at sNN= 8 TeV. The symbols indicate the experimental data of the LHCb Collaboration [25]. The solid lines indicate our theoretical calculations.
(a)
(b)
(c)Figure4 presents the double-differential cross section d2σ/dpTdy of J/ψ mesons produced in p-p collisions at sNN=13 TeV. The experimental data are from [27]. The symbols and lines represent the same meanings as those in Figure 1. The results are also in agreement with the experiment data. The T and q values taken for the calculation are listed in Table 2. The temperature T decreases with increasing the rapidity bins. As the emission source draws closer to the center, the excitation degree gets larger. The T values are larger than that at sNN=8 TeV in the same y range. The q behavior is similar to that in Figures 1 and 2.Figure 4
Same as Figure2, but for sNN=13 TeV. The experimental data are taken from the LHCb Collaboration [26].
(a)
(b)For comparison, Figure5 presents the double-differential cross section d2σ/dpTdy of J/ψ mesons in p-Pb collisions at sNN=5 TeV. The experimental data are from [8]. The symbols and lines represent the same meanings as those in Figure 1. The results are also in agreement with the experiment data. The parameters T and q taken for the calculation are listed in Table 2. The temperature T of the emission source decreases with increasing the rapidity bins. The q values are between 1.031 and 1.065. In the calculation, the least-square-fitting testing provides a statistical indication of the most probable values of the two parameters. The values of χ2/dof corresponding to the curves in Figures 1–5 are given in Tables 1 and 2. By the comparison between the model results and experimental data, we obtain the temperature of the emission source. For J/ψ particles in the same y bins, the temperature increases with increasing the collision energy. In p-p collision at sNN=8 TeV, the temperature is smaller for Υ than that for J/ψ in the same y ranges.Figure 5
Same as Figure2, but for p-Pb collisions at sNN=5 TeV. The experimental data are taken from the LHCb Collaboration [27].
(a)
(b)


## 4. Conclusions

Final-state particle production in high energy collisions has attracted much attention since the attempt has been made to understand the properties of strongly coupled QGP. Thermal-statistical models have been successful in describing particle yields in various systems at different energies. The emission source temperature is very important for understanding the matter evolution inp-p collisions at high energy. In the rapidity space, different final-state particles are emitted from different positions due to stronger longitudinal flow. In our previous work [12], we have studied the transverse momentum spectra of strange particles produced in A-A collisions at sNN=62.4 and 200 GeV in the improved fireball model. The temperature T of the emission source was characterized indirectly by the excitation degree, which varies with location in the cylinder. In the present work, we can directly extract the specific temperature by Tsallis statistics.In this work, we have embedded consistently the Tsallis statistics into the picture of the multisource thermal model for descriptions of the transverse momentum spectra ofJ/ψ and Υ mesons produced in pp and p-Pb collisions over an energy range from 5 to 13 TeV, with different rapidity intervals. In most cases, the model results are in good agreement with the experimental data of the LHCb Collaboration. The agreement is not very good at pT>12 GeV in Figures 1-2. In other collisions, the maximum χ2/dof value is 0.531 and the minimum value is 0.124. By comparing with the experimental data, the temperature of the emission source and the nonequilibrium degree of the collision system are extracted. In p-p and p-Pb collisions, the temperature decreases with increasing the rapidity bins. The closer the emission source is to the y center, the larger the excitation degree is. For J/ψ particles in the same y bins, the temperature T increases with increasing the collision energy. So, the excitation degree of the emission source also increases with increasing the collision energy. In p-p collision at sNN=8 TeV, the temperature is smaller for Υ than that for J/ψ in the same y ranges. The parameter q does not change significantly, which means the collision system is not very unstable.In summary, the transverse momentum spectra ofJ/Ψ and Υ mesons produced in p-p and p-Pb collisions at high energies have been studied by combining a picture of the particle-production source with Tsallis statistics. The model results are compared with the experimental data of the LHCb Collaboration. Our investigations show the improved model is successful in the description of the transverse momentum distributions. At the same time, it can offer the temperature of the emission source and the information about the nonequilibrium degree in the collisions.



---

*Source: 9383540-2017-05-18.xml*