# Pseudorapidity Distribution of Charged Particles and Square Speed of Sound Parameter inp

Authors: Ya-Qin Gao; Tian Tian; Li-Na Gao; Fu-Hu Liu
Journal: Advances in High Energy Physics
 (2014)

---

## Abstract

Pseudorapidity distributions of charged particles produced in proton-proton (p
-
p) or proton-antiproton (p
-
p
-) collisions over an energy range from 0.053 to 7 TeV are studied by using the four-component Landau hydrodynamic model. The results calculated by the model are in agreement with the experimental data of the UA5, PHOBOS, UA1, P238, CDF, ALICE, and CMS Collaborations which present orderly from low to high energies. According to the distribution widths of different components, the values and some features of square speed of sound parameter c
s
2 for “participant” and “spectator” quark components are obtained. It is shown that the speed of sound for “participant” quark components agrees approximately with that for “spectator” quark components in the error ranges. The present work is useful for studying nucleus-nucleus collisions in the related energy range.


---

## Body

## 1. Introduction

In the laboratory conditions, the only way for creating and studying the interacting systems with hadronic or partonic degrees of freedom at extremely high energy and density is to investigate heavy-ion (nucleus-nucleus) collisions at very high center-of-mass energies. The relativistic heavy ion collider (RHIC) in the United States and the large hadron collider (LHC) in Switzerland have been built, respectively [1, 2]. Particularly, the LHC was originally designed to accelerate two protons up to total energy of 14 TeV, and it will accelerate heavy ions to collide at center-of-mass energy (s
N
N) of 5.5 TeV. Such high energy collisions offer a new environment which leads to a new significant extension of the kinematic range in longitudinal rapidity and transverse momentum [3]. It also provides a new chance for us to understand systematically the particle statistical behavior, production process, interaction mechanism, and correlation phenomena.The multiplicity and (pseudo)rapidity distributions of final-state particles can often be used to test different theoretical models and wonderful ideas. The rapidity distributiond
N
ch
/
d
y and the pseudorapidity distribution d
N
ch
/
d
η of charged particles are very important quantities for the study of mainly particle production mechanism in high energy nucleon-nucleon and nucleus-nucleus collisions. The studies of  d
N
ch
/
d
y and d
N
ch
/
d
η  in proton-proton (p-p) or proton-antiproton (p-p
¯) collisions not only can provide reference and baseline for heavy-ion collisions, but also can reflect the main particles contribution to heavy-ion collisions. In p-p or p-p
¯ collisions at present colliders, the center-of-energies reach a region of ultrahigh energies. Meanwhile, the highest energy nuclear experiments can be performed. Most of  high-density experimental data can help us realize the creation of a new state of matter, namely, the quark gluon plasma (QGP) [4, 5], which is a thermalized system consisting of strong coupling quarks and gluons in a very small region. It is believed that the QGP was the main matter in our early universe shortly after the Big Bang.Many models have been introduced in the field of high energy collisions, for example, the hydrodynamics model [6], the thermodynamics model [7], the hadron resonance gas model [8], the fireball model [9, 10], and so forth. Generally speaking, a given model treats the different collisions by different considerations, and different models approach the same collisions in different methods. In our previous work [11], based on the participant-spectator model [12] and Landau hydrodynamic model [13–17], the pseudorapidity distributions of charged particles produced in nucleus-nucleus collisions at the RHIC and LHC have been analyzed, and the values of square speed of sound (c
s
2) have been obtained. We notice that the related data for p-p or p-p
¯ collisions are absent in our previous work [11].To make up the limitation of our previous work [11], in this work, we will use the same method to study the pseudorapidity distributions of charged particles produced in p-p or p-p
¯ collisions over an energy range from 0.053 to 7 TeV and to extract the values of speed of sound parameter for the “participant” and “spectator” quark components which correspond, respectively, to the participant and spectator nucleons in nucleus-nucleus collisions. The quoted experimental data are taken from the UA5 (p-p
¯ collisions at 0.053, 0.2, 0.546, and 0.9 TeV) [18], PHOBOS (p-p collisions at 0.2 and 0.41 TeV) [19], UA1 (p-p
¯ collisions at 0.54 TeV) [20], P238 (p-p
¯ collisions at 0.63 TeV) [21], CDF (p-p
¯ collisions at 0.63 and 1.8 TeV) [22], ALICE (p-p collisions at 0.9 and 2.36 TeV) [23], and CMS (p-p collisions at 0.9, 2.36, and 7 TeV) [24, 25] Collaborations.


## 2. The Model

The picture of participant-spectator model [12] is the base of many theoretical models such as the nuclear fireball model [26], the two-component model [27], the nuclear overlap model [28], and the multisource thermal model [29, 30]. In the framework of participant-spectator model [12], the projectile and target nuclei penetrating through each other in collisions. Then, a central participant region, a target spectator region, and a projectile spectator region are expected to form in the collisions. Further, we can divide the participant region into two parts: a target participant region and a projectile participant region. In the multisource thermal model [29, 30], these regions can extend to different cylinders in rapidity space. In addition, the contributions of leading nucleons can fall under the spectator regions. In high energy nucleus-nucleus collisions, both the participant and spectator regions consist of nucleons.In high energy hadron-hadron collisions, the situation is different. If we use the picture of participant-spectator model [12], both the participant and spectator regions will be quarks. Similarly, we have four source components: the target spectator quark component, the target participant quark component, the projectile participant quark component, and the projectile spectator quark component, which are presented orderly from low to high in the rapidity space. According to [31, 32], if we use the Landau hydrodynamic model [13–17] to describe the hadron-hadron collisions, only the central pseudorapidity region can be fitted. This means that the Landau hydrodynamic model [13–17] gives a narrow pseudorapidity distribution, and the contributions of  leading nucleons are not included in the model. To describe the pseudorapidity distribution in a wide region, we can apply the model for each quark component in the hadron-hadron collisions.The pseudorapidity distribution of charged particles produced in each quark component is given by a Gaussian function [17, 33]
(1)
d
N
ch
d
η
=
K
(
s
N
N
)
1
/
2
2
π
σ
η
2
exp
⁡
⁡
(
-
η
2
2
σ
η
2
)
with
(2)
σ
η
=
8
3
c
s
2
1
-
c
s
4
ln
⁡
⁡
(
s
N
N
2
m
p
)
,
where K is the normalization constant, s
N
N is in the units of GeV, c
s denotes the speed of sound, m
p denotes the proton mass in the units of GeV/c
2 with c
=
1, and σ
η is the pseudorapidity distribution width.Strictly speaking, there are some differences between the rapidity distribution and the pseudorapidity distribution. In fact, the differences between the two distributions at high energies can be neglected due to the small values [34]. If we need to distinguish the two distributions, we should calculate them severally or make a complete transformation [35]. In the present work, we do not distinguish the two distributions for the purpose of simpleness. Equation (1) can be rewritten as the function of c
s
2
(3)
d
N
ch
d
η
=
K
(
s
N
N
)
1
/
2
[
16
π
3
c
s
2
1
-
c
s
4
ln
⁡
⁡
(
s
N
N
2
m
p
)
]
-
1
/
2
v
·
exp
⁡
⁡
{
-
η
2
[
16
3
c
s
2
1
-
c
s
4
ln
⁡
⁡
(
s
N
N
2
m
p
)
]
-
1
}
due to (2). The final-state pseudorapidity distribution is a weighted sum of four equations (3) with different rapidity shifts and weights.Forp-p or p-p
¯ collisions, the rapidity shift (y
TS) of target spectator is opposite to the rapidity shift (y
PS) of projectile spectator, and the rapidity shift (y
TP) of target participant is opposite to the rapidity shift (y
PP) of projectile participant. The weight (k) of target participant is equal to that of projectile participant, and the weight [
(
1
-
2
k
)
/
2
] of target spectator is equal to that of projectile spectator.


## 3. Comparisons with Experimental Data

Figures1(a)–1(d) show the pseudorapidity distributions of charged particles produced in inelastic p-p
¯ collisions at s
N
N = 0.053, 0.2, 0.546, and 0.9 TeV, respectively. The symbols in the range of η
>
0 represent the experimental data of the UA5 Collaboration [18], and those in the range of  η
<
0 are symmetrical reflection at the midpseudorapidity. The curves are our calculated results with the four-component Landau hydrodynamic model. The values of parameters obtained by fitting the experimental data are given in Table 1 with the values of χ
2
/
dof (χ
2 per degree of freedom), where c
s
2
(
P
) and c
s
2
(
S
) denote the square speeds of sound for the participants and spectators, respectively. One can see that the calculated results are in agreement with the experimental data at the four energies. Both values of c
s
2 for the participants and spectators are approximately equal to each other in the error ranges, and they do not depend obviously on s
N
N in the considered energy range.Table 1
Values ofc
s
2, other parameters, and χ
2/dof corresponding to the curves in Figures 1–3.
Figure
Collisions
s
N
N  (TeV)
c
s
2  (P)
c
s
2  (S)
y
PP
y
PS
k
χ
2/dof
Figure1(a)
p
-
p
-
0.053
0.24
±
0.06
0.24
±
0.06
1.50
±
0.08
3.10
±
0.50
0.472
±
0.016
0.062
Figure1(b)
p
-
p
-
0.200
0.19
±
0.02
0.19
±
0.07
1.81
±
0.06
4.00
±
0.05
0.419
±
0.012
0.064
Figure1(c)
p
-
p
-
0.546
0.20
±
0.03
0.20
±
0.10
2.01
±
0.07
4.30
±
0.09
0.437
±
0.015
0.038
Figure1(d)
p
-
p
-
0.900
0.22
±
0.04
0.21
±
0.09
2.01
±
0.11
3.50
±
0.05
0.398
±
0.015
0.026
Figure2(a)
p
-
p
0.200
0.20
±
0.01
0.18
±
0.02
1.81
±
0.04
3.90
±
0.20
0.419
±
0.009
0.009
Figure2(b)
p
-
p
0.410
0.20
±
0.01
0.21
±
0.03
1.94
±
0.03
3.95
±
0.15
0.402
±
0.010
0.009
Figure2(c)
p
-
p
-
0.540
0.20
±
0.05
0.18
±
0.04
1.78
±
0.07
4.00
±
0.50
0.420
±
0.020
0.023
Figure2(d)
p
-
p
-
0.630
0.22
±
0.04
0.21
±
0.06
1.81
±
0.19
3.95
±
0.05
0.343
±
0.024
0.017
Figure2(e)
p
-
p
-
0.630
0.15
±
0.03
0.22
±
0.05
1.79
±
0.05
3.95
±
0.02
0.318
±
0.011
0.035
Figure2(f)
p
-
p
-
1.800
0.16
±
0.03
0.24
±
0.05
1.99
±
0.07
4.00
±
0.33
0.321
±
0.011
0.009
Figure3(a)
p
-
p
0.900
0.21
±
0.06
0.21
±
0.05
2.00
±
0.03
3.00
±
0.17
0.417
±
0.011
0.002
Figure3(b)
p
-
p
0.900
0.18
±
0.03
0.20
±
0.05
2.00
±
0.03
3.50
±
0.30
0.417
±
0.008
0.011
Figure3(c)
p
-
p
2.360
0.22
±
0.06
0.26
±
0.03
2.00
±
0.03
3.98
±
0.09
0.326
±
0.004
0.005
Figure3(d)
p
-
p
2.360
0.18
±
0.04
0.24
±
0.03
2.00
±
0.04
4.00
±
0.13
0.295
±
0.006
0.024
Figure3(e)
p
-
p
7.000
0.18
±
0.05
0.23
±
0.02
2.00
±
0.04
4.10
±
0.09
0.321
±
0.004
0.004The pseudorapidity distributions,d
N
ch
/
d
η, of charged particles produced in inelastic p-p
¯ collisions at s
N
N
=
0.053 (a), 0.2 (b), 0.546 (c), and 0.9 TeV (d). The symbols represent the experimental data of the UA5 Collaboration [18], and the curves are our calculated results.
(a)
(b)
(c)
(d)The pseudorapidity distributions of charged particles produced in inelasticp-p collisions at s
N
N
=
0.2 and 0.41 TeV measured by the PHOBOS Collaboration [19], p-p
¯ collisions at s
N
N
=
0.54 TeV measured by the UA1 Collaboration [20], p-p
¯ collisions at s
N
N
=
0.63 TeV measured by the P238 Collaboration with |
η
| ranging from 1.5 to 5.5 [21], and p-p
¯ collisions at s
N
N
=
0.63 and 1.8 TeV measured by the CDF Collaboration with |
η
|
=
0
–
3.5 [22] are presented in Figures 2(a)–2(f), respectively. The symbols and curves represent the same meanings as those in Figure 1. By fitting the experimental data, the obtained values of c
s
2 and χ
2
/
dof are given in Table 1. Once again the four-component Landau hydrodynamic model describes the experimental data well. Although the energy has a large change, the values of c
s
2 for the participants and spectators are approximately the same in the error ranges, and they have no obvious change with changing the energy.The pseudorapidity distributions of charged particles in inelasticp-p collisions at s
N
N
=
0.2 (a) and 0.41 TeV (b), as well as in inelastic p-p
¯ collisions at s
N
N
=
0.54 (c), 0.63 (d, e), and 1.8 TeV (f). The symbols represent the experimental data of the PHOBOS [19] (a, b), UA1 [20] (c), P238 [21] (d), and CDF Collaborations [22] (e, f), respectively. The curves are our calculated results.
(a)
(b)
(c)
(d)
(e)
(f)Figure3 gives the pseudorapidity distributions of charged particles produced in non-single-diffractive p-p collisions at s
N
N = 0.9, 2.36, and 7 TeV. The symbols in the left panel represent the experimental data obtained by the CMS Collaboration [24, 25], and those in the right panel stand for the experimental data of the ALICE Collaboration [23]. The curves are our calculated results. The obtained values of c
s
2 and χ
2
/
dof are given in Table 1. We notice that the model describes the experimental data at the LHC energies. Both values of c
s
2 for the participants and spectators seem to be the same in the error ranges, and they do not depend obviously on the energy.The pseudorapidity distributions of charged particles in non-single-diffractivep-p collisions at s
N
N
=
0.9 (a, b), 2.36 (c, d), and 7 TeV (e). The symbols represent the experimental data of the CMS [24, 25] (a, c, e) and ALICE Collaborations [23] (b, d), and the curves are our calculated results.
(a)
(b)
(c)
(d)
(e)To see clearly the dependences ofc
s
2
(
P
) and c
s
2
(
S
) on s
N
N, the values of c
s
2
(
P
) and c
s
2
(
S
) in Table 1 are displayed in Figure 4. Different symbols represent the results extracted from different experiments as marked in the figure. The open and closed symbols represent the values of c
s
2
(
P
) and c
s
2
(
S
), respectively. One can see that both the c
s
2
(
P
) and c
s
2
(
S
) fall mainly into a range consisting of two lines: c
s
2
=
0.12 and c
s
2
=
0.30. These results are in agreement with the hadron resonance gas model which results in, respectively, c
s
2
≈
0.23 and c
s
2
≈
0.12 for the existing regions of hadron resonances including and excluding pions at the temperature of ~85 MeV and c
s
2
≈
0.14
–
0.15 for the two existing regions at the temperature of ~190 MeV [36, 37]. The present results are also in agreement with the lattice quantum chromodynamics theory which gives c
s
2
=
0.12
–
0.16 and 0.31 at the temperatures of 125 and 400 MeV, respectively [38].Figure 4
Correlations betweenc
s
2 and s
N
N. The closed and open symbols represent, respectively, the values of c
s
2
(
P
) and c
s
2
(
S
) obtained in Figures 1–3 and listed in Table 1. The two lines described by c
s
2
=
0.12 and c
s
2
=
0.30 show a range for different c
s
2.Both values ofc
s
2 for the participants and spectators are approximately the same in the error ranges and they do not depend obviously on s
N
N over an energy range from 0.053 to 7 TeV. From Figure 4, one can also see that c
s
2
(
P
) seems to be less than c
s
2
(
S
) at TeV energies. Because the fitted errors are large for both c
s
2, it is hard to determine which one is less than the other.The correlations betweeny
PP and ln
⁡
⁡
s
N
N as well as y
PS and ln
⁡
s
N
N are presented in Figures 5(a) and 5(b), respectively, where s
N
N  is in the units of GeV. The symbols represent the values of rapidity shifts y
PP and y
PS which are obtained in Figures 1–3 and listed in Table 1. The lines are our fitted results which are described by
(4)
y
PP
=
(
0.071
±
0.010
)
ln
⁡
⁡
s
N
N
+
(
1.460
±
0.071
)
with χ
2
/
dof
=
0.514 and
(5)
y
PS
=
(
0.017
±
0.030
)
ln
⁡
⁡
s
N
N
+
(
3.852
±
0.199
)
with χ
2
/
dof
=
0.863, respectively. One can see the linear relation existing between the rapidity shifts and ln
⁡
⁡
s
N
N. With the increase of the logarithmic center-of-mass energy, both rapidity shifts increase slightly.Correlations betweeny
PP and ln
⁡
s
N
N (a), as well as y
PS and ln
⁡
s
N
N (b), where s
N
N is in the units of GeV. The symbols represent the values of rapidity shifts y
PP and y
PS obtained in Figures 1–3 and listed in Table 1 and the lines are our fitted results.
(a)
(b)The correlation betweenk and ln
⁡
⁡
s
N
N is shown in Figure 6, where s
N
N is in the units of GeV. The symbols represent the values of contribution ratio of the projectile (or target) participant quark component which are obtained in Figures 1–3 and listed in Table 1. The line is our fitted result described by
(6)
k
=
-
(
0.026
±
0.001
)
ln
⁡
⁡
s
N
N
+
(
0.547
±
0.011
)
with χ
2
/
dof
=
0.646. One can see the linear relation existing between the contribution ratio and ln
⁡
⁡
s
N
N. Particularly, the contribution ratio of participants decreases and that of spectators increases with the increase of the logarithmic center-of-mass energy.Figure 6
Correlations betweenk and ln
⁡
s
N
N, where s
N
N  is in the units of GeV. The symbols represent the values of contribution ratio k of projectile (or target) participant quark component obtained in Figures 1–3 and listed in Table 1, and the line is our fitted result.


## 4. Conclusions

To conclude, the pseudorapidity distributions of charged particles produced inp-p collisions at s
N
N
=
0.2, 0.41, 0.9, 2.36, and 7 TeV and in p
-
p
¯ collisions at s
N
N
=
0.053, 0.2, 0.54, 0.546, 0.63, 0.9, and 1.8 TeV have been studied in the present work. The calculated results of the four-component Landau hydrodynamical model satisfactorily describe the experimental data of the UA5, PHOBOS, UA1, P238, CDF, ALICE, and CMS Collaborations. In the model, we have considered the application of the participant-spectator model [12] in hadron-hadron collisions. Then, the interacting system of hadron-hadron collisions is divided into four components: the target and projectile “spectator” quark components and the target and projectile “participant” quark components. For each component, the Landau hydrodynamical model [13–17] is applied, which results in the pseudorapidity distribution being a Gaussian type.Based on the above successful descriptions of the pseudorapidity distributions of charged particles produced inp-p or p-p
¯ collisions over an energy range from 0.053 to 7 TeV, the square speeds of sound, c
s
2, for both the participants and spectators are obtained and found to be the same in the fitted error ranges. Both c
s
2 are approximately independent of s
N
N in the considered energy range and fall mainly into the range of 0.12–0.30 which is in agreement with hadron resonance gas model which gives c
s
2
≈
0.12
–
0.23 at the temperature of ~85 MeV and c
s
2
≈
0.14
–
0.15 at ~190 MeV [36, 37]. The values of c
s
2 obtained in the present work are also in agreement with the lattice quantum chromodynamics theory which gives c
s
2
≈
0.12
–
0.31 in the temperature range of  125–400 MeV [38].In the process of extracting the speed of sound, the rapidity shifts and contribution ratio of the participant and spectator quark components are naturally obtained. There are linear relations existing betweeny
PP and ln
⁡
⁡
s
N
N, y
PS and ln
⁡
⁡
s
N
N, as well as k and ln
⁡
⁡
s
N
N. Both rapidity shifts increase slightly with the increase of the logarithmic center-of-mass energy. The contribution ratio of participants decreases and that of the spectators increases with the increase of the logarithmic center-of-mass energy.


## 5. Discussions

The basis for Landau’s equation (essentially a Gaussian pseudorapidity distribution) is the applicability of Landau hydrodynamics, namely, the creation of an initial state of dense Lorentz contracted matter that expands later according to the equations of 2  +  1 dimensional relativistic hydrodynamics. It seems that such an initial state is not created in hadron-hadron collision at intermediate or not too high energy, and the evolution of the interacting system would not be covered by hydrodynamics. However, hadron-hadron collision at high or ultrahigh energy can be regarded as a small dense system comparing with nucleus-nucleus collisions. This small dense system can produce a multihadron final-state state and can be approximately described by Landau’s equation.The square speed of sound extracted from the present work is in agreement with that from nucleus-nucleus collisions at RHIC and LHC energies in our previous work [11]. The similar values of speed of sound in both hadron-hadron and nucleus-nucleus collisions provides the insight into some common laws and universal mechanisms in multihadron productions [31, 32]. We believe that there are more similar quantities in multihadron productions in the “elementary” and “complex” collisions.We would like to point out that the present work is the first one which extracts the square speed of sound in hadron-hadron collisions from the pseudorapidity distributions of charged particles over a very wide energy range. What we use is the four-component Landau hydrodynamic model which is based on the participant-spectator model [12] and the Landau hydrodynamic model [13–17]. Our work can provide references and baselines for nucleus-nucleus collisions at high energies. We hope to extract the square speed of sound in proton-proton and lead-lead collisions at the highest energy at the LHC in the near future.



---

*Source: 569079-2014-04-22.xml*