Properties of Chemical and Kinetic FreezeOuts in HighEnergy Nuclear Collisions
View this Special IssueResearch Article  Open Access
Z. J. Jiang, Dongfang Xu, Yan Huang, "A Description of Pseudorapidity Distributions of Charged Particles Produced in Au+Au Collisions at RHIC Energies", Advances in High Energy Physics, vol. 2018, Article ID 1369098, 8 pages, 2018. https://doi.org/10.1155/2018/1369098
A Description of Pseudorapidity Distributions of Charged Particles Produced in Au+Au Collisions at RHIC Energies
Abstract
In heavy ion collisions, charged particles come from two parts: the hot and dense matter and the leading particles. In this paper, the hot and dense matter is assumed to expand according to the hydrodynamic model including phase transition and decouples into particles via the prescription of CooperFrye. The leading particles are as usual supposed to have Gaussian rapidity distributions with the number equaling that of participants. The investigations of this paper show that, unlike low energy situations, the leading particles are essential in describing the pseudorapidity distributions of charged particles produced in high energy heavy ion collisions. This might be due to the different transparencies of nuclei at different energies.
1. Introduction
The BNL Relativistic Heavy Ion Collider (RHIC) accelerates nuclei up to the centerofmass energies from a dozen GeV to 200 GeV per nucleon. In the past decade, the measurements from such collisions have triggered an extensive research for the properties of matter at extreme conditions of very high temperature and energy densities [1–33]. One of the most important achievements from such research is the discovery that the matter created in nucleusnucleus collisions at RHIC energies is in the state of strongly coupled quarkgluon plasma (sQGP) exhibiting a clear collective behavior nearly like a perfect fluid with very low viscosity [10–33].
The best approach for describing the spacetime evolution of fluidlike sQGP is the relativistic hydrodynamics. However, owing to the formidable complexities of hydrodynamic equations, the most analytical work so far is mainly limited to 1+1 expansion for a perfect fluid with simple equation of state, which can be found as an important application in the analysis of the pseudorapidity distributions of charged particles in high energy physics. In this paper, combining the effects of leading particles, we will discuss such distributions in the framework of 1+1 hydrodynamic model including phase transition [10].
2. A Brief Introduction to the Model
Here, for the purpose of completeness and applications, we will list the key ingredients of the hydrodynamic model [10].
The expansion of fluid is subject to the conservation of energy and momentum. This is reflected in continuity equationwhere is the energymomentum tensor. For a perfect fluidwhereis the 4velocity of fluid and is its rapidity. and in (2) are the energy density and pressure of fluid, which meet the thermodynamical relationswhere and are the temperature and entropy density of fluid, respectively. To close (1), another relation, namely, the equation of stateis needed, where is the sound speed of fluid, which takes different values in sQGP and in hadronic phase.
Project (1) to the direction of and the direction perpendicular to , respectively. This leads to equationsEquation (6) is the continuity equation for entropy conservation. Equation (7) means the existence of a scalar function satisfying relationsFrom and Legendre transformation, Khalatnikov potential is introduced via relationIn terms of , the variables and can be expressed aswhere is the initial temperature of fluid and . Through the above equations, the coordinate base of is transformed to that of , and (6) is translated into the socalled telegraphy equation
Along with the expansions of matter created in collisions, it becomes cooler and cooler. As its temperature drops from the initial to the critical , phase transition occurs. The matter transforms from sQGP state to hadronic state. The produced hadrons are initially in the violent and frequent collisions. The major part of these collisions is inelastic. Hence, the abundances of identified hadrons are changing. Furthermore, the mean free paths of these primary hadrons are very short. The movement of them is still like that of a fluid meeting (11) with only difference being the value of . In sQGP, , which is the sound speed of a massless perfect fluid, being the maximum of . In the hadronic state, . At the point of phase transition, that is, as , is discontinuous.
The solution of (11) for the sector of sQGP is [10]where is a constant determined by tuning the theoretical results to experimental data and is the 0thorder modified Bessel function of the first kind, and
In the sector of hadrons, the solution of (11) is [10]where
3. The Pseudorapidity Distributions of Charged Particles
(1) The Invariant Multiplicity Distributions of Charged Particles Frozen Out from sQGP. From Khalatnikov potential , the rapidity distributions of charged particles frozen out from fluidlike sQGP read [34]where is the area of overlap region of collisions, being the function of impact parameter or centrality cuts. Inserting (10) into the above equation, the part in the round brackets becomes
With the expansions of hadronic matter, it continues becoming cooler. According to the prescription of CooperFrye [34], as the temperature drops to the freezeout temperature , the inelastic collisions among hadrons cease. The yields of identified hadrons remain unchanged becoming the measured results in experiments. The invariant multiplicity distributions of charged particles equal [10, 15, 34]where is the transverse mass of produced charged particle with rest mass . in (18) is the baryochemical potential. For Fermi charged particles, in the denominator of (18), and for Bosons, . The meaning of (18) is evident. It is the convolution of with the energy of the charged particles in the state with temperature .
The integral interval of in (18) is . The integrand is evaluated with . At this moment, the fluid freezes out into the charged particles. Replacing in (17) by of (14), it becomeswherewhere is the 1storder modified Bessel function of the first kind.
(2) The Invariant Multiplicity Distributions of Leading Particles. Investigations have shown that the leading particles are formed outside the overlap region of collisions [35, 36]. The generation and movement of them are therefore beyond the scope of hydrodynamic description and should be treated separately.
In our previous work [24–26], we once argued that the rapidity distributions of leading particles take the Gaussian formwhere and are, respectively, the central position and width of distributions. It can be expected that should increase with increasing energies and centrality cuts, while should not, at least not apparently, depend on the energies, centrality cuts, and even colliding systems. The specific values of them can be determined by comparing the theoretical results with experimental data. in (21) represents the number of leading particles, which, for an identical nucleusnucleus collision, equals half of the number of participants.
The investigations have shown that [37], for certain rapidity, the invariant multiplicity distributions of leading particles possess the formwhere is a constant. Then, as a function of rapidity, the invariant multiplicity distributions of leading particles can be written aswhich is normalized to .
(3) The Pseudorapidity Distributions of Charged Particles. Writing invariant multiplicity distributions in terms of pseudorapidity, we havewhereTo fulfill the transformation of (24), another relationis in order.
Substituting (18) and (23) into (24) and carrying out the integration of , we can get the pseudorapidity distributions of charged particles produced in high energy heavy ion collisions. Figures 1 and 2 show such distributions in Au+Au collisions at and 62.4 GeV, respectively. The solid dots in the figures are the experimental measurements [5]. The dashed, dasheddotted, and dotted curves are, respectively, the contribution from pions, kaons, and protons got from the hydrodynamic result of (18). The dottedstar curves are the components of leading particles obtained from (23). The solid curves are the sums of the four types of curves. for each curve is listed in Table 1. It can be seen that the combined contribution from both hydrodynamics and leading particles matches up well with experimental data.

Experiments have shown that the overwhelming majority of charged particles produced in Au+Au collisions at GeV consists of pions, kaons, and protons with proportions of about 84%, 12%, and 4%, respectively [38], which are roughly independent of energies, centrality cuts, and colliding systems. In calculations, the ratios of these three kinds of particles take about the same as these values. in (15) takes the values of and 0.42 for and 62.4 GeV from the investigations of [15, 39–41]. in (15) takes the wellrecognized value of MeV. The freezeout temperature takes the values of MeV from the studies of [6], which also shows that the baryochemical potential in (18) is about equal to 20 and 50 MeV for and 62.4 GeV, respectively. For the most central collisions at these two different energies, in (15) takes the values of and 0.68 GeV referring to those given in [15]. This allows us to determine the constant in (16) to be , , and for pions, kaons, and protons, respectively. Keeping unchanged, is fixed for the rest centrality cuts by making theoretical results fit in with experimental data. The results are listed in Table 1. It can be seen that decreases slowly with increasing centralities especially in the first four cuts. Table 1 also lists the central position in (21). As addressed above, it increases with increasing energies and centralities. The width parameter in (21) takes the value of a constant of , being independent of energies and centrality cuts. The parameter in (23) takes the value of for the two different energies.
Figure 3 shows the pseudorapidity distributions of charged particles produced in Au+Au collisions at GeV. for each curve is listed in Table 1. The meanings of different types of curves are the same as those in Figures 1 and 2. It can be seen that, in the absence of leading particles, the hydrodynamics alone can give a good description to the experimental observations. This is different from Figures 1 and 2, where leading particles are essential in fitting experimental data. This difference might be caused by the different transparencies of nuclei in different energies. As the analyses given in [42], in central Au+Au collisions at GeV, the leading particles are located at about . This position is far away from the midrapidity region, where, relative to the low yields of charged particles frozen out from sQGP, the effect of leading particles is evident which should be considered separately. On the contrary, in case of Au+Au collisions at GeV, . This position is so close to the midrapidity region that the effect of leading particles is hidden by the large yields of charged particles generated from the freezeout of sQGP. Therefore, there is no need to consider the contribution of leading particles separately.
In drawing Figure 3, takes the values as those listed in Table 1. and MeV. The other parameters, such as , , and , are the same as those used in drawing Figures 1 and 2.
4. Conclusions
By taking into consideration the effect of leading particles, the hydrodynamic model incorporating the phase transition is used to analyze the pseudorapidity distributions of charged particles produced in Au+Au collisions at RHIC energies.
The hydrodynamic model contains rich information about transport coefficients of sQGP, such as the sound speed in sQGP, the sound speed in hadronic phase, the phase transition temperature , the chemical freezeout temperature , the baryochemical potential , and the initial temperature . With the exception of , the other five coefficients take the values either from the wellknown theoretical results or from experimental measurements. As for , there are no widely accepted results so far. In our calculations, in the most central Au+Au collisions at and 62.4 GeV takes the value referring to those given by other investigations, which enables us to ascertain the constant in (16). In the rest centrality cuts and in Au+Au collisions at GeV, is determined by maintaining unchanged and comparing the theoretical results with experimental data.
The leading particles, by conventional definition, are the particles carrying on the quantum numbers of colliding nucleons and taking away the most part of incident energy. They are separately in projectile and target fragmentation region. The present investigations show that the importance of leading particles in describing the pseudorapidity distributions of charged particles produced in heavy ion collisions is related to the incident energy. At high energy, owing to the high transparency of nuclei, the contribution of leading particles is evident and indispensable, while, at low energy, as a result of poor transparency of nuclei the effect of leading particles is integrated with the results of freezeout of sQCD. It does not need to be dealt with separately.
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this paper.
Acknowledgments
This work is supported by the Shanghai Key Lab of Modern Optical System.
References
 C. Adler and STAR Collaboration, “Measurement of inclusive antiprotons from Au+Au collisions at = 130 GeV,” Physical Review Letters, vol. 87, Article ID 262302, 2001. View at: Google Scholar
 K. Adcox and PHENIX Collaboration, “Centrality dependence of , , , and production from = 130 GeV Au+Au collisions at RHIC,” Physical Review Letters, vol. 88, Article ID 242301, 2002. View at: Google Scholar
 S. S. Adler and PHENIX collaboration, “Elliptic flow of identified hadrons in Au+Au collisions at = 200 GeV,” Physical Review Letters, vol. 91, Article ID 182301, 2003. View at: Google Scholar
 K. Adcox and PHENIX Collaboration, “Single identified hadron spectra from = 130 GeV Au+Au collisions,” Physical Review C, vol. 69, Article ID 024904, 2004. View at: Google Scholar
 B. Alver and PHOBOS Collaboration, “Chargedparticle multiplicity and pseudorapidity distributions measured with the PHOBOS detector in Au+Au, Cu+Cu, d+Au, and p+p collisions at ultrarelativistic energies,” Physical Review C, vol. 83, Article ID 024913, 2011. View at: Google Scholar
 B. I. Abelev and STAR Collaboration, “Systematic measurements of identified particle spectra in p + p, d + Au, and Au+Au collisions at the STAR detector,” Physical Review C, vol. 79, Article ID 034909, 2009. View at: Google Scholar
 A. Adare and PHENIX Collaboration, “Spectra and ratios of identified particles in Au+Au and d+Au collisions at = 200 GeV,” Physical Review C, vol. 88, Article ID 024906, 2013. View at: Google Scholar
 L. C. Arsene and BRAHMS collaboration, “Rapidity and centrality dependence of particle production for identified hadrons in Cu+Cu collisions at = 200 GeV,” Physical Review C, vol. 94, Article ID 014907, 2016. View at: Google Scholar
 L. Adamczyk and STAR Collaboration, “Centrality dependence of identified particle elliptic flow in relativistic heavy ion collisions at = 7.7–62.4 GeV,” Physical Review C, vol. 93, Article ID 014907, 2016. View at: Google Scholar
 N. Suzuki, “Onedimensional hydrodynamical model including phase transition,” Physical Review C nuclear physics, vol. 81, no. 4, Article ID 044911, 2010. View at: Publisher Site  Google Scholar
 A. Bialas, R. A. Janik, and R. Peschanski, “Unified description of Bjorken and Landau 1+1 hydrodynamics,” Physical Review C nuclear physics, vol. 76, no. 5, Article ID 054901, 2007. View at: Publisher Site  Google Scholar
 G. Beuf, R. Peschanski, and E. N. Saridakis, “Entropy flow of a perfect fluid in (1 + 1) hydrodynamics,” Physical Review C: Nuclear Physics, vol. 78, no. 6, Article ID 064909, 2008. View at: Publisher Site  Google Scholar
 T. Csörgo, M. I. Nagy, and M. Csanád, “New family of simple solutions of relativistic perfect fluid hydrodynamics,” Physics Letters B, vol. 663, no. 4, pp. 306–311, 2008. View at: Publisher Site  Google Scholar
 M. Csanád, M. I. Nagy, and S. Lökös, “Exact solutions of relativistic perfect fluid hydrodynamics for a QCD Equation of State,” The European Physical Journal A, vol. 48, no. 11, pp. 173–178, 2012. View at: Publisher Site  Google Scholar
 T. Mizoguchi, H. Miyazawa, and M. Biyajima, “A potential including the Heaviside function in the 1 + 1 dimensional hydrodynamics by Landau : IIts basic properties and application to data at RHIC energies,” The European Physical Journal A, vol. 40, no. 1, pp. 99–108, 2009. View at: Publisher Site  Google Scholar
 T. R. Taylor and C. Vafa, “Analytic solution for relativistic transverse flow at the softest point,” Physics Letters B, vol. 474, no. 12, pp. 21–26, 2000. View at: Publisher Site  Google Scholar
 T. S. Bíró, “Generating new solutions for relativistic transverse flow at the softest point,” Physics Letters B, vol. 487, no. 12, pp. 133–139, 2000. View at: Publisher Site  Google Scholar
 C.Y. Wong, “Landau hydrodynamics reexamined,” Physical Review C: Nuclear Physics, vol. 78, no. 5, Article ID 054902, 2008. View at: Publisher Site  Google Scholar
 A. Bialas and R. Peschanski, “Asymmetric (1+1)dimensional hydrodynamics in highenergy collisions,” Physical Review C: Nuclear Physics, vol. 83, no. 5, Article ID 054905, 2011. View at: Publisher Site  Google Scholar
 E. K. G. Sarkisyan and A. S. Sakharov, “Relating multihadron production in hadronic and nuclear collisions,” The European Physical Journal C, vol. 70, no. 3, pp. 533–541, 2010. View at: Publisher Site  Google Scholar
 C. Gale, S. Jeon, and B. Schenke, “Hydrodynamic modeling of heavyion collisions,” International Journal of Modern Physics A, vol. 28, no. 11, Article ID 1340011, 2013. View at: Publisher Site  Google Scholar
 U. Heinz and R. Snellings, “Collective flow and viscosity in relativistic heavyion collisions,” Annual Review of Nuclear and Particle Science, vol. 63, pp. 123–151, 2013. View at: Publisher Site  Google Scholar
 A. N. Mishra, R. Sahoo, E. K. G. Sarkisyan, and A. S. Sakharov, “Effectiveenergy budget in multiparticle production in nuclear collisions,” The European Physical Journal C, vol. 74, article 3147, 2014. View at: Google Scholar
 Z. J. Jiang, Q. G. Li, and H. L. Zhang, “Revised Landau hydrodynamic model and the pseudorapidity distributions of charged particles produced in nucleusnucleus collisions at maximum energy at the BNL Relativistic Heavy Ion Collider,” Physical Review C: Nuclear Physics, vol. 87, no. 4, Article ID 044902, 2013. View at: Publisher Site  Google Scholar
 Z. J. Jiang, Y. Zhang, H. L. Zhang, and H. P. Deng, “A description of the pseudorapidity distributions in heavy ion collisions at RHIC and LHC energies,” Nuclear Physics A, vol. 941, pp. 188–200, 2015. View at: Publisher Site  Google Scholar
 Z. W. Wang, Z. J. Jiang, and Y. S. Zhang, “The investigations of pseudorapidity distributions of final multiplicity in Au+Au collisions at high energy,” Journal of University of Shanghai for Science and Technology, vol. 31, pp. 322–326, 2009. View at: Google Scholar
 J. Qian and U. Heinz, “Hydrodynamic flow amplitude correlations in eventbyevent fluctuating heavyion collisions,” Physical Review C: Nuclear Physics, vol. 94, no. 2, Article ID 024910, 2016. View at: Publisher Site  Google Scholar
 E. K. G. Sarkisyan, A. N. Mishra, R. Sahoo, and A. S. Sakharov, “Multihadron production dynamics exploring the energy balance in hadronic and nuclear collisions,” Physical Review D: Particles, Fields, Gravitation and Cosmology, vol. 93, no. 7, Article ID 079904, 2016. View at: Google Scholar
 E. K. G. Sarkisyan, A. N. Mishra, R. Sahoo, and A. S. Sakharov, “Centrality dependence of midrapidity density from GeV to TeV heavyion collisions in the effectiveenergy universality picture of hadroproduction,” Physical Review D: Particles, Fields, Gravitation and Cosmology, vol. 94, no. 1, Article ID 011501, 2016. View at: Publisher Site  Google Scholar
 J. NoronhaHostler, M. Luzum, and J.Y. Ollitrault, “Hydrodynamic predictions for 5.02 TeV PbPb collisions,” Physical Review C nuclear physics, vol. 93, no. 3, Article ID 034912, 2016. View at: Publisher Site  Google Scholar
 H. Niemi, K. J. Eskola, and R. Paatelainen, “Eventbyevent fluctuations in a perturbative QCD + saturation + hydrodynamics model: Determining QCD matter shear viscosity in ultrarelativistic heavyion collisions,” Physical Review C: Nuclear Physics, vol. 93, no. 2, Article ID 024907, 2016. View at: Publisher Site  Google Scholar
 F. G. Gardim, F. Grassi, M. Luzum, and J. NoronhaHostler, “Hydrodynamic predictions for mixed harmonic correlations in 200 GeV Au+Au collisions,” Physical Review C: Nuclear Physics, vol. 95, no. 3, Article ID 034901, 2017. View at: Publisher Site  Google Scholar
 M. Alqahtani, M. Nopoush, and M. Strickland, “Quasiparticle anisotropic hydrodynamics for central collisions,” Physical Review C: Nuclear Physics, vol. 95, no. 3, Article ID 034906, 2017. View at: Publisher Site  Google Scholar
 F. Cooper and G. Frye, “Landau's hydrodynamic model of particle production and electronpositron annihilation into hadrons,” Physical Review D: Particles, Fields, Gravitation and Cosmology, vol. 11, pp. 192–213, 1975. View at: Publisher Site  Google Scholar
 A. Berera, M. Strikman, W. Toothacker, W. Walker, and J. Whitmore, “The limiting curve of leading particles from hadronnucleus collisions at infinite A,” Physics Letters B, vol. 403, no. 12, pp. 1–7, 1997. View at: Publisher Site  Google Scholar
 J. J. Ryan, “Proceedings of the Annual Meeting of the Division of Particles and Fields of the APS,” World Scientific, Singapore, 1993. View at: Google Scholar
 A. E. Brenner, D. C. Carey, J. E. Elias et al., “Experimental study of singleparticle inclusive hadron scattering and associated multiplicities,” Physical Review D: Particles, Fields, Gravitation and Cosmology, vol. 26, no. 7, pp. 1497–1553, 1982. View at: Publisher Site  Google Scholar
 S. S. Adler and PHENIX collaboration, “Identified charged particle spectra and yields in Au+Au collisions at = 200 GeV,” Physical Review C, vol. 69, Article ID 034909, 2004. View at: Google Scholar
 A. Adare and PHENIX Collaboration, “Scaling properties of azimuthal anisotropy in Au+Au and Cu+Cu collisions at = 200 GeV,” Physical Review Letters, vol. 98, Article ID 162301, 2007. View at: Google Scholar
 L. N. Gao, Y. H. Chen, H. R. Wei, and F. H. Liu, “Speed of sound parameter from RHIC and LHC heavyion data,” Advances in High Energy Physics, vol. 2013, Article ID 450247, 8 pages, 2013. View at: Publisher Site  Google Scholar
 S. Borsányi, G. Endrodi, Z. Fodor et al., “The QCD equation of state with dynamical quarks,” Journal of High Energy Physics, vol. 2010, no. 11, article 77, pp. 1–31, 2010. View at: Publisher Site  Google Scholar
 Z.J. Jiang, H.L. Zhang, J. Wang, K. Ma, and L.M. Cai, “The evolutiondominated hydrodynamics and the pseudorapidity distributions in nucleusnucleus collisions at low energies at the BNL relativistic heavy ion collider,” Chinese Journal of Physics, vol. 52, no. 6, pp. 1676–1685, 2014. View at: Google Scholar
Copyright
Copyright © 2018 Z. J. Jiang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. The publication of this article was funded by SCOAP^{3}.