Particle Production at CBM Energies in a Thermal Model Approach
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).
Relativistic heavy ion collisions provide a unique way to create and investigate the QCD matter at extreme temperature () and/or density () [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 . 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 seconds . 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 (), temperature (), 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 .
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 Section 3 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] where is the degeneracy factor of th species of baryons, is the energy of the particle (), is the eigenvolume of one baryon of th species and is the total occupied volume, and represents total number of baryons of th species.
Now we can write (1) as follows: where represents the integral: and is the fugacity of the particle, is the number density of th type of baryons after excluded volume correction, and it can be obtained from (2) as follows: This leads to a transcendental equation: where is the partial derivative of with respect to and is the fractional occupied volume depending on . We can write in an operator equation as follows : where with ; is the density of pointlike baryons of th species and the operator has the following form: Using Neumann iteration method and retaining the series up to term, we get Equation (8) can be solved numerically. Finally, we get the total pressure [24, 25] of the hadron gas as follows:
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 for each type of baryon with a hard-core radius fm. We have taken all baryons and mesons and their resonances having masses up to in our calculation for the HG pressure. We have also imposed the condition of strangeness conservation by putting , where is the strangeness quantum number of the th hadron and is the strange (antistrange) hadron density, respectively. Using this constraint equation, we get the value of strange chemical potential in terms of .
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 volume , energy density ,and baryon density , which are in fact related to and 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 and of the fireball to the experimental variable . In our model, we vary the and of our model and observe how the experimental data of particle ratio are fitted at a given . Similarly, we repeat the exercise for all the data available ranging from SIS to RHIC energy and thus get a set of and for a particular . Now after such fitting procedure, we get a parameterization between and or as follows: where the parameters , , , , and have been determined from the best fit: , , GeV, , and . The details of the fitting procedure in thermal models can be found in . We have also done this fitting for various other hadron gas models in order to obtain the values of and with respect to and all this is shown in Table I of .
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 . Figure 1 demonstrates the variation of with . To calculate , we first use our parameterization given by (10) to get the values of and for given . At these and , we have calculated the value of net baryon density . In the present excluded volume model, the net baryon density increases with (see Figure 1), reaches a maximum value near -9 GeV, and then starts decreasing. The maximum freeze-out net baryon density is approximately half of the normal nuclear density, . We have compared our model results with the net baryon density obtained by Begun et al.  using excluded volume model based on . They have calculated the baryon density at various using fm as hard-core radius () for each baryon and all the mesons are treated as pointlike particles, that is, (). On the other hand, we have taken the value of 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. . The maximum value of obtained in our model as well as in the excluded volume model of Begun et al. is also lower than the value obtained in  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 2 presents the variation of total multiplicity of with respect to . In the CBM energy range, the fireball volume at freeze-out, extracted in the excluded volume model approach, appears almost constant  for all emitted particles. We have taken 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  and UrQMD [39, 40] models usually employ the concepts of string, quark, and diquark, () 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 GeV as well as mesonic resonances up to GeV [39, 40]. However, HSD incorporates only the baryon octet and decuplet states and , as well as their antiparticles together with the and 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 GeV. However, all three models yield larger multiplicity for relative to the data above 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 ) in which some variations can always be thought to appear.
Figure 3 demonstrates the variation of multiplicity with respect to . We have used the same freeze-out volume, that is, , 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 GeV. The overall level of agreement with the excluded volume model results is quite good. However, HSD results also suitably describe the data below 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.
In Figure 4, we present the variation of total multiplicity of with . 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 GeV. Above GeV, the total multiplicity of obtained from HSD and UrQMD is relatively small and does not agree with the data.
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 and quark would be highly enhanced compared to that of light or quark [45, 46, 50, 51]. This is possibly due to the Pauli exclusion principle which strongly suppresses the further creation of light-quark pairs  in the medium having abundant , quarks. This asymmetry in the flavor composition generated by a baryon-rich QGP should result in a large production of , , , , and so forth . 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 ratio at 30 AGeV . This sharp maximum which is also known as “horn” is not seen in p+p collisions. As 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 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 . A similar maximum at the same beam energy is also reported by the same collaboration for other strange particles like ’s and . These observations confirm that this particular feature is not given by alone, but they represent the total strangeness content of the final state . 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.
Figure 5 demonstrates the variation of ratio with respect to . 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  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  have suggested that the overestimation of multiplicity in thermal models gives theoretical curve lying below the experimental data. The similar fireball volume used for 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.
In Figure 6, we present the variation of with respect to . Our model suitably describes the data, but below 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 to GeV.
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 () and that of antikaons to kaons () are the representatives of two such significant observables measuring the hadron-antihadron asymmetry in heavy ion collisions [60, 61]. The ratio carries the information regarding baryon-antibaryon asymmetry and the ratio almost cancels the effect of strangeness production and indicates the asymmetry between charged mesons and their antiparticles generated in the hot, dense medium.
In Figures 7 and 8, we have shown the variations of and with , 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 and production is larger as also observed between and . 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  also.
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, , , using a constant freeze-out volume which is equal to for all the hadrons. We have also calculated the ratio of and in CBM energy range. Almost all the models fail to reproduce the “horn” in ratio. Furthermore, we have also calculated the particle to antiparticle ratio like and . 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.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
The authors gratefully acknowledge Professor C. P. Singh for many stimulating and inspiring discussions on the paper. One of the authors, P. K. Srivastava, is grateful to the University Grants Commission (UGC), New Delhi, for providing a research fellowship. A. Prakash and B. K. Singh are grateful to the Department of Science and Technology (DST), Government of India, New Delhi, for providing a research grant.
B. Friman, C. Hohne, J. Knoll et al., Eds., The CBM Physics Book, vol. 814 of Lecture Notes in Physics, Springer, Berlin, Germany, 2011.
L. Adamczyk, J. K. Adkins, G. Agakishiev et al., “Observation of an energy-dependent difference in elliptic flow between particles and anti-particles in relativistic heavy ion collisions,” Physical Review Letters, vol. 110, no. 14, Article ID 142301, 7 pages, 2013.View at: Google Scholar
A. Laszlo, PoS(CPOD07), 054, 2007.
F. Becattini, J. Cleymans, A. Keranen, E. Suhonen, and K. Redlich, “Features of particle multiplicities and strangeness production in central heavy ion collisions between 1.7A and 15A GeV/c,” Physical Review C, vol. 64, no. 2, Article ID 024901, 9 pages, 2001.View at: Publisher Site | Google Scholar
C. Alt, T. Anticic, B. Baatar et al., “Pion and kaon production in central Pb+Pb collisions at 20A and 30A GeV: evidence for the onset of deconfinement,” Physical Review C, vol. 77, no. 2, Article ID 024903, 10 pages, 2008.View at: Google Scholar
C. Alt, T. Anticic, B. Baatar et al., “Energy dependence of Λ and Ξ production in central Pb+Pb collisions at 20A, 30A, 40A, 80A, and 158A GeV measured at the CERN Super Proton Synchrotron,” Physical Review C, vol. 78, no. 3, Article ID 034918, 15 pages, 2008.View at: Publisher Site | Google Scholar
S. Scherer, S. A. Bass, M. Bleicher et al., “Critical review of quark gluon signatures,” Progress in Particle and Nuclear Physics, vol. 42, pp. 279–293, 1999.View at: Google Scholar
V. Friese, PoS(CPOD09), 005, 2009.