# Particle Production at CBM Energies in a Thermal Model Approach

Authors: A. Prakash; P. K. Srivastava; B. K. Singh
Journal: Advances in High Energy Physics
 (2014)

---

## Abstract

The compressed baryonic matter (CBM) experiment planned at the Facility for Antiproton
and Ion Research (FAIR) will provide a major scientific effort for exploring the properties
of strongly interacting matter in the high baryon density regime. One of the important
goals behind such experiment is to precisely determine the equation of state (EOS) for the
strongly interacting matter at extremely large baryon density. In this paper, we have used
some successful models for RHIC and LHC energies to predict different particle ratios and
the total multiplicity of various hadrons in the CBM energy range, that is, from 10 A GeV to
40 A GeV lab energies, which corresponds to 4.43 A GeV and 8.71 A GeV center-of-mass
energies. Our main emphasis is to estimate the strange particles enhancement as well as an
increase in the net baryon density at CBM experiment. We have also compared the model
results with the experimental data obtained at alternating gradient synchrotron (AGS) and
super proton synchrotron (SPS).


---

## Body

## 1. Introduction

Relativistic heavy ion collisions provide a unique way to create and investigate the QCD matter at extreme temperature (T) and/or density (ρ
B) [1, 2]. These experiments are expected to produce short lived bubbles of medium, in which the quarks and gluons are the active degrees of freedom [3]. However, the analysis of such collision events is complicated due to a transient nature of the reaction. In these reactions, the object of interest is the formation of a fireball which consists of a compressed and/or heated zone and it is of finite small size (order of few Fermi) and further it fades away in a very short duration of less than 1
0
-
22 seconds [4]. The only observation comes through the particles (most abundantly hadrons) produced and emitted from the collision zone into the detectors. The main features of particle multiplicity, that is, the variation of number of particles produced with respect to collision energy, the momentum spectra, and the relative abundances of different species, have led us in the understanding of the nature and properties of strong interactions [5, 6]. Several experimental programs were planned to explore the properties of strongly interacting matter and to search the possible existence of QGP, for example, alternating gradient synchrotron (AGS) and relativistic heavy-ion collider (RHIC) experiments at BNL [7–9], super proton synchrotron (SPS) and large hadron collider (LHC) at CERN [10–12], and so forth. The compressed baryonic matter (CBM) experiment at Facility for Antiproton and Ion Research (FAIR) machine will provide a similar hot, dense situation in the laboratory. The aim of the CBM experiment is to obtain information on (i) the properties of hadrons in dense or hot baryonic/hadronic matter, (ii) the restoration of chiral symmetry at high temperatures and high baryon densities, (iii) the deconfinement phase transition from hadronic to quark-gluon matter at high temperatures and/or high baryon densities, and (iv) the nuclear equation of state at high baryon densities [13, 14].These experimental efforts need reliable theoretical guidance to proceed in the right direction. Lattice QCD which is the first principle method to calculate the properties of QCD provides some quality information on the strongly interacting matter at finite temperature and zero (or very small) chemical potential (μ) [15–17]. However, lattice QCD approach is severely limited at finite chemical potential due to a notorious sign problem. Perturbative calculations are only useful at very high temperatures as the QCD coupling constant becomes small. In recent years, relativistic fluid dynamics and statistical mechanics emerged as two important theoretical tools to explore the properties of hot and dense QCD matter created in the heavy ion collision experiments. Standard statistical models or thermal models assume a system in global thermodynamical equilibrium. It means that the intensive variables like pressure (P), temperature (T), and chemical potential (μ) are constant throughout the volume of the system. The most popular thermal model is ideal hadron gas (IHG) model and its modified versions are known as excluded volume models. These models assume that the hadrons are produced according to their statistical weights from an equilibrated hadron gas at freeze-out. These models were used in the recent past in order to deduce the multiplicities and ratios of hadrons, and agreements with the experimental data were found to be excellent [18–21]. On the other hand, hydrodyanmical or transport models assume that the system created at initial time is in local thermodynamic equilibrium [22, 23]. It means that the thermodynamical variables such as pressure and temperature of the system vary with space-time point. However, the variation is so slow that, for any point, one can assume thermodynamic equilibrium in some neighborhood about that point. After this initial time, the system expands and cools off until the hadronic freeze-out occurs and produces the hadrons. These models again show good agreement with the data obtained from RHIC and LHC [23].In this paper, our main motivation is to show what one should expect in context to particle multiplicities and particle ratios at CBM energies, in view of equilibrated statistical as well as transport models. We further compare the model results with the available experimental data of AGS and SPS energies. The rest of the paper is organized as follows: in the second section, we will provide a brief description of a recently developed thermal model with excluded volume correction. Further Section3 will provide the results regarding the hadron multiplicities and their ratios and comparison among various models. We mainly emphasize the net baryon density created at freeze-out, the production of strange particles and the asymmetry in the particle-antiparticle production, and so forth because these issues are still unsettled in order to get information on the QCD phase transition. The last section will be devoted to the conclusion and future prospects drawn from this study.


## 2. EOS for a Hot and/or Dense Hadron Gas

Recently we have proposed [24–28] a thermodynamically consistent excluded volume model for the hot and dense hadron gas (HG). In this model, the grand canonical partition function for the HG with full quantum statistics and after suitably incorporating excluded volume correction is [24–28]
(1)
ln
⁡
Z
i
ex
=
g
i
6
π
2
T
∫
V
i
0
V
-
∑
j
N
j
V
j
0
‍
d
V
×
∫
0
∞
‍
k
4
d
k
k
2
+
m
i
2
1
[
exp
⁡
(
(
E
i
-
μ
i
)
/
T
)
+
1
]
,
where g
i is the degeneracy factor of ith species of baryons, E
i is the energy of the particle (E
i
=
k
2
+
m
i
2), V
i
0 is the eigenvolume of one baryon of ith species and ∑
j
‍
N
j
V
j
0 is the total occupied volume, and N
j represents total number of baryons of jth species.Now we can write (1) as follows:
(2)
ln
⁡
Z
i
ex
=
V
(
1
-
∑
j
‍
n
j
ex
V
j
0
)
I
i
λ
i
,
where I
i represents the integral:
(3)
I
i
=
g
i
6
π
2
T
∫
0
∞
‍
k
4
d
k
k
2
+
m
i
2
1
[
exp
⁡
(
E
i
/
T
)
+
λ
i
]
,
and λ
i
=
exp
⁡
(
μ
i
/
T
) is the fugacity of the particle, n
j
ex is the number density of jth type of baryons after excluded volume correction, and it can be obtained from (2) as follows:
(4)
n
i
ex
=
λ
i
V
(
∂
ln
⁡
Z
i
ex
∂
λ
i
)
T
,
V
.
This leads to a transcendental equation:
(5)
n
i
ex
=
(
1
-
R
)
I
i
λ
i
-
I
i
λ
i
2
∂
R
∂
λ
i
+
λ
i
2
(
1
-
R
)
I
i
′
,
where I
i
′ is the partial derivative of I
i with respect to λ
i and R
=
∑
i
‍
n
i
ex
V
i
0 is the fractional occupied volume depending on n
i
ex. We can write R in an operator equation as follows [24]:
(6)
R
=
R
1
+
Ω
^
R
,
where R
1
=
R
0
/
(
1
+
R
0
) with R
0
=
∑
‍
n
i
0
V
i
0
+
∑
‍
I
i
′
V
i
0
λ
i
2; n
i
0 is the density of pointlike baryons of ith species and the operator Ω
^ has the following form:
(7)
Ω
^
=
-
1
1
+
R
0
∑
i
‍
n
i
0
V
i
0
λ
i
∂
∂
λ
i
.
Using Neumann iteration method and retaining the series up to Ω
^
2 term, we get
(8)
R
=
R
1
+
Ω
^
R
1
+
Ω
^
2
R
1
.
Equation (8) can be solved numerically. Finally, we get the total pressure [24, 25] of the hadron gas as follows:
(9)
p
HG
ex
=
T
(
1
-
R
)
∑
i
‍
I
i
λ
i
+
∑
i
‍
P
i
meson
.In (9), the first term represents the pressure due to all types of baryons, where excluded volume correction is incorporated and the second term gives the total pressure due to all mesons in HG having a pointlike size. This makes it clear that we consider the hard-core repulsion arising between two baryons only. Essentially, we consider that the mesons can interpenetrate each other but baryons can not owing to their hard-core size. In this calculation, we have taken an equal volume V
0
=
4
π
r
3
/
3 for each type of baryon with a hard-core radius r
=
0.8 fm. We have taken all baryons and mesons and their resonances having masses up to 2
GeV
/
c
2 in our calculation for the HG pressure. We have also imposed the condition of strangeness conservation by putting ∑
i
‍
S
i
(
n
i
s
-
n
-
i
s
)
=
0, where S
i is the strangeness quantum number of the ith hadron and n
i
s
(
n
-
i
s
) is the strange (antistrange) hadron density, respectively. Using this constraint equation, we get the value of strange chemical potential in terms of μ
B.In relativistic heavy ion collisions, a hot, dense matter is formed and it is often referred to as a “fireball.” The physical variables of the fireball are the volumeV, energy density ϵ,and baryon density ρ
B, which are in fact related to T and μ
B of the fireball. When cooling or expansion of the fireball starts, it goes through two types of freeze-out stages, namely, chemical freeze-out and thermal freeze-out. When inelastic collisions between constituents of the fireball cease, we call this the chemical freeze-out stage. Later when elastic collisions are also stopped, this is called the thermal freeze-out. In our model, we suppose that at the time of chemical freeze-out all the hadrons are emitted from the entire freeze-out hypersurface. Now the main point is to correlate thermodynamical variables T and μ
B of the fireball to the experimental variable s
N
N. In our model, we vary the T and μ
B of our model and observe how the experimental data of particle ratio are fitted at a given s
N
N. Similarly, we repeat the exercise for all the data available ranging from SIS to RHIC energy and thus get a set of T and μ
B for a particular s
N
N. Now after such fitting procedure, we get a parameterization between s
N
N and μ
B or T as follows:
(10)
μ
B
=
a
1
+
b
s
N
N
,
T
=
c
-
d
μ
B
2
-
e
μ
B
4
,
where the parameters a, b, c, d, and e have been determined from the  best  fit:  a
=
1.482
±
0
.
003
GeV,  b
=
0.351
±
0
.
009
Ge
V
-
1,  c
=
0.163
±
0
.
002 GeV,  d
=
0.170
±
0
.
020
Ge
V
-
1, and e
=
0.015
±
0
.
010
Ge
V
-
3. The details of the fitting procedure in thermal models can be found in [18]. We have also done this fitting for various other hadron gas models in order to obtain the values of T and μ
B with respect to s
N
N and all this is shown in Table I of [24].


## 3. Results and Discussion

The extracted freeze-out temperature in statistical thermal models of HG generally increases monotonically with the collision energy. However, the corresponding net baryon density exhibits a more complicated behavior [36]. Figure 1 demonstrates the variation of ρ
B with s
N
N. To calculate ρ
B, we first use our parameterization given by (10) to get the values of μ
B and T for given s
N
N. At these T and μ
B, we have calculated the value of net baryon density ρ
B. In the present excluded volume model, the net baryon density increases with s
N
N (see Figure 1), reaches a maximum value near s
N
N
=
8-9 GeV, and then starts decreasing. The maximum freeze-out net baryon density is approximately half of the normal nuclear density, ρ
0
=
0.15
f
m
-
3. We have compared our model results with the net baryon density obtained by Begun et al. [37] using excluded volume model based on [38]. They have calculated the baryon density at various s
N
N using 0.5 fm as hard-core radius (r
b) for each baryon and all the mesons are treated as pointlike particles, that is, (r
m
=
0). On the other hand, we have taken the value of 0.8 fm as the hard-core radius for each baryon. The quantitative values of maximum freeze-out density achieved in heavy ion collisions are almost equal in both models. However, the position of the maximum in the net baryon density is somewhat shifted towards the lower energy side in the results obtained by Begun et al. [37]. The maximum value of ρ
B obtained in our model as well as in the excluded volume model of Begun et al. is also lower than the value obtained in [36] where an HRG model is used without any excluded volume correction. Obviously, the excluded volume effect shifts the net baryon density achieved at freeze-out towards the lower value.Figure 1
Variation of net baryon density (ρ
B) at freeze-out with respect to center-of-mass energy (s
N
N).Figure2 presents the variation of total multiplicity of π
+ with respect to s
N
N. In the CBM energy range, the fireball volume at freeze-out, extracted in the excluded volume model approach, appears almost constant [28] for all emitted particles. We have taken 5000 f
m
3 as the fireball volume in order to calculate the total multiplicity of hadrons. We compare our results with the experimental data obtained at AGS and SPS [30–35] at low energies. We have also shown the total multiplicity of π
+ obtained from transport models like HSD and UrQMD. Both HSD [29] and UrQMD [39, 40] models usually employ the concepts of string, quark, and diquark, (q
,
q
-
,
q
q
,
q
-
q
-) as well as the hadronic degrees of freedom. However, the numerical evaluations are quite different in HSD as compared to UrQMD. The UrQMD includes all baryonic resonances up to an invariant mass of 2 GeV as well as mesonic resonances up to 1.9 GeV [39, 40]. However, HSD incorporates only the baryon octet and decuplet states and N
*
(
1440
), N
*
(
1535
) as well as their antiparticles together with the 0
- and 1
- meson octets. Higher baryonic resonances are discarded as the resonance structure (above Δ peak) is not clearly seen experimentally even in the photoabsorption by light nuclei [41–43]. In contrast to the UrQMD at low energy baryon-baryon and meson-baryon collisions, HSD includes the direct (nonresonant) meson production. Our excluded volume model suitably describes the data up to 5 GeV. However, all three models yield larger multiplicity for π
+ relative to the data above 5 GeV, although the multiplicity of the produced π
+ in excluded volume model is closer to the data in comparison to other two models. In principle, the yield of each particle is mainly governed by the particle fugacity essentially determined from the chemical freeze-out parameters. It also depends on the size of the system (or volume V) in which some variations can always be thought to appear.Figure 2
Variation of total multiplicity of producedπ
+ with respect to s
N
N. Dash-dotted and dotted curves are the results obtained from HSD and UrQMD models, respectively [29]. Experimental data is taken from [30–35].Figure3 demonstrates the variation of K
+ multiplicity with respect to s
N
N. We have used the same freeze-out volume, that is, 5000
f
m
3, as that used in the π
+ multiplicity. We compare our model results with the results obtained from HSD and UrQMD simulations. We have also shown the experimental data from CERN-SPS and RHIC-AGS for comparison [30–35]. The results obtained from both HSD and UrQMD do not match with the data beyond the energy 5 GeV. The overall level of agreement with the excluded volume model results is quite good. However, HSD results also suitably describe the data below 5 GeV. It should be emphasized here that the thermal statistical model should be used with caution at lower energies since the number of produced particles is very small.Figure 3
Variation of total multiplicity of producedK
+ with respect to s
N
N. Dash-dotted and dotted curves are the results obtained from HSD and UrQMD models, respectively [29]. Experimental data is taken from [30–35].In Figure4, we present the variation of total multiplicity of K
- with s
N
N. The results obtained from our model are in excellent agreement with the data [30–35]. However, HSD and UrQMD results satisfy the experimental data only below 5 GeV. Above 5 GeV, the total multiplicity of K
- obtained from HSD and UrQMD is relatively small and does not agree with the data.Figure 4
Variation of total multiplicity of producedK
- with respect to s
N
N. Dash-dotted and dotted curves are the results obtained from HSD and UrQMD models, respectively [29]. Experimental data is taken from [30–35].Strangeness enhancement has been proposed as one of the early and important signals of QGP formation [45–49]. It has been argued that if a quark-gluon plasma is formed from compressed nuclear matter as it may happen in the nuclear fragmentation region and/or in the low energy “stopping regime,” then the abundance of s and s
- quark would be highly enhanced compared to that of light u or d quark [45, 46, 50, 51]. This is possibly due to the Pauli exclusion principle which strongly suppresses the further creation of light-quark pairs [52] in the medium having abundant u, d quarks. This asymmetry in the flavor composition generated by a baryon-rich QGP should result in a large production of K
+, K
-, Λ, Λ
-, and so forth [53]. This effect is even more evident in the case of the multistrange hyperons. A striking observation reported by NA49 collaboration is a pronounced and a sharp maximum in the excitation function of K
+
/
π
+ ratio at 30 AGeV [31]. This sharp maximum which is also known as “horn” is not seen in p+p collisions. As K
+ is by far the most abundant carrier of antistrangeness at SPS energies, it also provides a good measure of the total strangeness produced in the collision. The ratio K
+
/
π
+ represents the strangeness to entropy ratio. A sharp maximum in this quantity was predicted by the statistical model indicating the early stage as a consequence of the transition to a deconfined state [53]. A similar maximum at the same beam energy is also reported by the same collaboration for other strange particles like Λ’s and Ξ
- [32]. These observations confirm that this particular feature is not given by K
+ alone, but they represent the total strangeness content of the final state [54]. The measurement of the excitation function of strangeness production by NA49 collaboration has renewed a fresh stimulating discussion about the role of strangeness as a signature for the deconfinement phase transition.Figure5 demonstrates the variation of K
+
/
π
+ ratio with respect to s
N
N. We have compared our model results with the experimental data obtained from SPS experiment [30–35]. We have further compared them with the other results like HSD, UrQMD, and Boltzmann-Uehling-Uhlenbeck (BUU) transport [44] models. We find that both HSD and UrQMD models fail to give agreement with the experimental data. However, BUU model shows better results in comparison to HSD and UrQMD. Our model yields better results in agreement with the experimental data. The authors in [29] have suggested that the overestimation of π
+ multiplicity in thermal models gives theoretical curve lying below the experimental data. The similar fireball volume used for K
+ and π
+ in the excluded volume model can also contribute to such disagreement. Here, we use the same freeze-out volume for the emission of all kinds of hadrons from the fireball surface. CBM experiment will definitely provide an important insight in the understanding of strange particle production mechanism and more vitally address the mechanism for the existence of the “horn” like behavior in the lower energy region.Figure 5
Variation ofK
+
/
π
+ ratio with respect to s
N
N. The results obtained from HSD and UrQMD are taken from [29]. The results obtained from BUU transport model are extracted from [44]. Experimental data is taken from [30–35].In Figure6, we present the variation of K
-
/
π
- with respect to s
N
N. Our model suitably describes the data, but below 5 GeV we again notice a large disagreement. However, HSD and UrQMD both also give complete disagreement with the experimental data almost in the entire energy range, that is, from 3 to 9 GeV.Figure 6
Variation ofK
-
/
π
- ratio with respect to s
N
N. The results obtained from HSD and UrQMD are taken from [29]. Experimental data is taken from [30–35].A large production of anti-baryons over baryons in heavy ion collisions are considered as the signal of QGP formation [55–59]. In heavy ion collisions, the system has a nonzero baryon number density arising due to nuclear stopping. At small and moderate center-of-mass energies as existing in the case of CBM experiment, the nuclear stopping is large in comparison to the RHIC and/or LHC energies where nuclear transparency is found to dominate. Nuclear stopping leads to an asymmetry between the production of hadrons and antihadrons in the produced particles other than the initial finite baryon number. However, there is a possibility that hadronization can also generate additional particle-antiparticle asymmetry. This asymmetry can be measured by the ratio of yields of antihadrons to hadrons [60, 61]. Ratios of yield of antiprotons to protons (p
-
/
p) and that of antikaons to kaons (K
-
/
K
+) are the representatives of two such significant observables measuring the hadron-antihadron asymmetry in heavy ion collisions [60, 61]. The ratio p
-
/
p carries the information regarding baryon-antibaryon asymmetry and the ratio K
-
/
K
+ almost cancels the effect of strangeness production and indicates the asymmetry between charged mesons and their antiparticles generated in the hot, dense medium.In Figures7 and 8, we have shown the variations of p
-
/
p and K
-
/
K
+ with s
N
N, respectively, as obtained in our excluded volume model with the center-of-mass energy. In both cases, the production of antiparticle to particle is less at lowest CBM energy. However, the asymmetry in p
- and p production is larger as also observed between K
- and K
+. As the energy increases, the production of antiparticle over particle increases since there is an increase in nuclear transparency. However, it is important to state here that the excluded volume model does not agree with the experimental data for the net nucleon density at RHIC highest energy [24] also.Figure 7
Variation ofp
-
/
p with respect to s
N
N. Experimental data is taken from [30–35].Figure 8
Variation ofK
-
/
K
+ with respect to s
N
N. The results obtained from HSD and UrQMD are taken from [29]. Experimental data is taken from [30–35].In summary, we have calculated the net baryon density at freeze-out in the CBM energy range which comes out to be the maximum achievable density in heavy ion collisions. Further, we have calculated the total multiplicities of various produced hadrons, for example,π
+, K
+, K
- using a constant freeze-out volume which is equal to 5000
f
m
3 for all the hadrons. We have also calculated the ratio of K
+
/
π
+ and K
-
/
π
- in CBM energy range. Almost all the models fail to reproduce the “horn” in K
+
/
π
+ ratio. Furthermore, we have also calculated the particle to antiparticle ratio like p
-
/
p and K
-
/
K
+. In CBM energy range, these ratios increase rapidly. However, due to lack of experimental data, we do not have any meaningful comparison, and hence we cannot precisely understand the asymmetry between hadron and antihadron production in heavy ion collisions.In conclusion, we face many worthwhile questions regarding particle production at CBM energy regime. More vital questions are regarding the formation of a thermal fireball and the applicability of thermal equilibrium in the statistical system. But it is interesting to investigate what we should precisely expect. So we should precisely know what we expect from various equlibrium-nonequilibrium models. Our excluded volume model gives a better agreement with the available experimental data in comparison to other models. In the CBM energy range, most of the multiplicities as well as particle ratios show somewhat peculiar behavior. However, there is scarcity of experimental data. The data available suffer from poor statistics also. The upcoming CBM experiment having a high luminosity beam will provide a unique opportunity to perform systematic and comprehensive measurements, with better statistics of bulk and rare particle behavior. This will help us to understand the particle production mechanism and also possibly find the existence of QGP at CBM energy.



---

*Source: 983861-2014-03-04.xml*