Research Article  Open Access
The EvolutionDominated Hydrodynamic Model and the Pseudorapidity Distributions in High Energy Physics
Abstract
By taking into account the effects of leading particles, we discuss the pseudorapidity distributions of the charged particles produced in high energy heavy ion collisions in the context of evolutiondominated hydrodynamic model. The leading particles are supposed to have a Gaussian rapidity distribution normalized to the number of participants. A comparison is made between the theoretical results and the experimental measurements performed by BRAHMS and PHOBOS Collaboration at BNLRHIC in AuAu and CuCu collisions at GeV and by ALICE Collaboration at CERNLHC in PbPb collisions at TeV.
1. Introduction
Along with the successful description of elliptic flow and multiplicity production in heavy ion collisions [1–4], relativistic hydrodynamics has now been widely accepted as one of the most important tools for understanding the spacetime evolution of the matter created in collisions [5–11]. With the specified initial conditions, the equation of state, and the freezeout conditions, the motion of fluid relies only on the local energymomentum conservation and the assumption of local thermal equilibrium. From this point of view, hydrodynamics is simple and powerful. However, on the other hand, the initial conditions, the equation of state, and the freezeout conditions of fluid are not well known. Worse still is that the partial differential equations of relativistic hydrodynamics are highly nonlinear and coupled. It is a very hard thing to solve them analytically. From this point of view, hydrodynamics is tremendously complicated. This is the reason why the progress in finding exact hydrodynamic solutions is not going well. Up till now, most of this work is only limited in dimensional flows for the perfect fluid with the simple equation of state [12–23]. The dimensional hydrodynamics is less developed, and no general exact solutions are known so far.
In the present paper, by using the evolutiondominated hydrodynamics [12] and taking into account the contribution from leading particles, we will discuss the pseudorapidity distributions of the charged particles produced in heavy ion collisions. We will first give a brief introduction to the evolutiondominated hydrodynamics in Section 2. The obtained solutions are then used in Section 3 to formulate the rapidity distributions of charged particles produced in heavy ion collisions. Then, in Section 4, a comparison is made between the theoretical results and experiment measurements carried out by BRAHMS and PHOBOS Collaboration at BNLRHIC in AuAu and CuCu collisions at GeV [24–27] and by ALICE Collaboration at CERNLHC in PbPb collisions at TeV [28]. The last section is traditionally about conclusions.
2. A Brief Introduction to EvolutionDominated Hydrodynamics
Here, for the purpose of completion and applications, we will give a brief introduction to the evolutiondominated hydrodynamics [12].
The motion of a perfect fluid obeys the equation where is the 4vector of spacetime and is the energymomentum tensor, is the metric tensor, and is the 4vector of fluid velocity. Noticing that where is the ordinary rapidity of fluid, the expansion of liquid obeys equations where is the lightcone coordinates, is the proper time, and is the spacetime rapidity of fluid.
In case of vanishing chemical potential, where and are, respectively, the temperature and entropy density of fluid. From the above two equations, we can get relation where where is an arbitrary initial temperature scale. Equation (6) means the existence of a potential satisfying where are the lightcone variables of fluid velocity. In this way, (6) is automatically fulfilled.
Equation (4) is a complicated, nonlinear, and coupled one. In order to solve it, one introduces Khalatnikov potential where are the functions of implicitly defined by in (8) and can be expressed by Khalatnikov potential as Following from this relation, one can get From (4), we can get the equation of as where where is the speed of sound. Equation (12) is now a linear secondorder partial differential equation, which works for any form of .
Experimental investigations have shown that the speed of sound is a constant of about or , which is almost independent of interaction energy and system [29–32]. In this case, we can take Equation (12) becomes The Green’s function of meets equation It has solution where is the Heaviside step function and is the 0th order modified Bessel function of the first kind. Thus, we arrive at where stands for the distributions of sources of hydrodynamic flow.
In heavy ion collisions at high energy, owing to the violent compression of collision system along beam direction, the initial pressure gradient of created matter in this direction is very large. By contrast, the effect of initial flow of sources is negligible. The motion of liquid is mainly dominated by the following evolution. The typical example reflecting such fact is the Landau hydrodynamic model [14, 15], where the fluid is assumed initially at rest. In this evolutiondominated case, the source function may take the form as [12, 15, 33] where is a constant. Inserting it into (18), we finally obtain the solution
3. The Rapidity Distributions in High Energy Heavy Ion Collisions
As an application of the Khalatnikov potential of (20), we will now derive the rapidity distribution of the charged particles produced in high energy heavy ion collisions. To this end, we first evaluate the entropy distribution at freezeout temperature as a function of rapidity .
The entropy distribution at freezeout temperature is defined as the amount of entropy flowing through the hypersurface with a fixed temperature in a unit rapidity interval. It has the form as [12] where is the 4dimensional unit vector of the hypersurface is the spacelike slab element along hypersurface with fixed temperature , which is defined as meeting where the minus sign accounts for the spacelike characteristic of .
In the base, the fixedtemperature hypersurface can be conveniently defined by The tangential vector of the hypersurface is where the primes represent the derivatives with regard to . According to definitions, we have Owing to (25), the above equation turns into This equation together with (22) gives Equation (25) translates (23) into Making use of (28), we obtain Using (29) and (30), (21) reads Furthermore, known from (10), Deduced from (11),These two equations make (31) become For evolutiondominated hydrodynamics, substituting (20) into the above equation, we acquire where is the 1st order modified Bessel function of the first kind. is related to the initial temperature of fluid and is therefore dependent on the incident energy and collision centrality. Since should not vary very much with centrality cuts, we can expect that would approximately maintain a constant or at least depend weakly on collision centrality for a given incident energy.
As the entropy is proportional to the number of produced charged particles, we obtain the rapidity distribution where , independent of rapidity , is an overall normalization constant. is the impact parameter, and is the centerofmass energy per pair of nucleons.
4. Comparison with Experimental Measurements and the Rapidity Distributions of Leading Particles
Figure 1 shows the rapidity distributions for , , , , , and produced in central AuAu collisions at GeV. The scattered symbols are the experimental data [24–26]. The solid curves are the theoretical results from (36). In calculations, the parameter takes the value of . We can see from this figure that, except for proton , (36) fits well with experimental measurements. For proton , experimental data show an evident uplift in the rapidity interval between and 3.0. This may result from parts of leading particles, which are free from the description of hydrodynamics. Hence, in order to match up with experimental data, we should take these leading particles into account separately.
Considering that, for a given incident energy, the leading particles in each time of nucleusnucleus collisions have approximately the same amount of energy, then, according to the central limit theorem [34, 35], the leading particles should follow the Gaussian rapidity distribution. That is, where and are, respectively, the central position and width of distribution. In fact, as is known to all, the rapidity distribution of any charged particles produced in heavy ion collisions can be well represented by Gaussian form ([24–26]; also confer the shapes of the curves in Figure 1). It is obvious that should increase with incident energy and centrality cut. However, should not apparently depend on them. This is due to the fact that the relative energy differences among leading particles should not be too much for different incident energies and centrality cuts. in (37) is the number of leading particles. It is a function of energy and centrality.
It is well known that, in nucleonnucleon, such as , collisions, there are two leading particles. One is in projectile fragmentation region, and the other is in target fragmentation region. Then, in nucleusnucleus collisions, the leading particles should be those nucleons which participate in collisions, the socalled participants, which locate separately at projectile and target fragmentation regions. For collisions between two identical nuclei, each nucleus should have about the same number of participants. Hence, the number of leading particles appearing in projectile or target fragmentation region should be where is the number of total participants in two nuclei, which can be evaluated by formula [36, 37] where is the coordinates in the overlap region measured from the center of one nucleus. The integrand in above equation: where is the inelastic nucleonnucleon crosssection. It increases slowly with energy. For example, for GeV, mb [38], and, for TeV, mb [39]. The subscripts and in the above equation denote the projectile and target nucleus, respectively. is the thickness function defined as where is the WoodsSaxon distribution of nuclear density. and are, respectively, the skin depth and radius of nucleus. In this paper, they take the values of fm and fm [36], where is the mass number of nucleus.
Tables 1 and 2 show the mean numbers of total participants in different centrality AuAu and CuCu collisions at GeV and PbPb collisions at TeV. The numbers with and without errors are those given by experiments [27, 28] and (39), respectively. Due to space constraints, Table 1 only shows the numbers in the first nine centrality cuts. It can be seen that both sets of numbers coincide well.


Having the rapidity distributions of (36) and (37), the pseudorapidity distribution measured in experiments can be expressed as [40] where is the transverse momentum, is the transverse mass, and is the total rapidity distribution from both fluid evolution and leading particles.
Substituting (45) into (43), we can get the pseudorapidity distributions of charged particles. Figures 2, 3, and 4 show such distributions in different centrality AuAu and CuCu collisions at GeV and PbPb collisions at 2.76 TeV, respectively. The solid dots in figures are the experimental measurements [27, 28]. The dashed curves are the results got from evolutiondominated hydrodynamics of (36). The dotted curves are the results obtained from leading particles of (37). The solid curves are the results achieved from (45), that is, the sums of dashed and dotted curves. It can be seen that the theoretical results are well consistent with experimental measurements.
In calculations, the parameter in (36) takes the values of 2.80 in the first three centrality cuts, 2.98 in the following six ones, and 3.17 in the last two ones in AuAu collisions. In CuCu collisions, takes the value of 2.95 in the first three centrality cuts, 3.15 in the following six ones, and 3.53 in the last three ones. In PbPb collisions, takes the value of 5.85 for the first two centrality cuts and 6.04 for the second two ones. It can be seen that increases with incident energy and centrality cut. The width parameter in (37) takes a constant of 0.85 for all three kinds of collision systems in different centrality cuts. As the analyses given above, is independent of incident energy and centrality cut. What is more, it is also independent of collision system. The center parameter in (37) takes the values as listed in Tables 1 and 2. As stated early, increases with energy and centrality cut. From Table 1, we can see that, for a given incident energy and centrality cut, decreases with increasing nucleus size. This can be understood if we notice the fact that the larger the nucleus size, the more the collisions that the participants will undergo. Hence, the final leading particles will lose more energy or have smaller . The fitting value of in the top 3% most central AuAu collisions is in accordance with the experimental observation shown in Figure 1, which indicates that the leading particles are mainly in the range between and 3. Experimental investigations also have shown that [41], in the top 5% most central AuAu collisions at GeV, the rapidity loss of participants is up to ; then the leading particles should locate at Seeing the smaller centrality cut considered in our analysis, our above fitting result is also consistent with this measurement.
5. Conclusions
The charged particles produced in heavy ion collisions are divided into two parts. One is from the hot and dense matter created in collisions. The other is from leading particles.
Compared with the effect of pressure gradient, the effect of initial flow of the hot and dense matter is negligible. The motion of this matter is mainly governed by the evolution of fluid. This thus guarantees the rationality of evolutiondominated hydrodynamics. With the scheme of Khalatnikov potential, this theoretical model can be solved exactly, and the rapidity distribution of charged particles can be expressed in a simple analytical form in terms of 0th and 1st order modified Bessel function of the first kind with only two parameters and . takes the value from experiments. is fixed by fitting with experimental data.
For leading particles, we assume that the rapidity distribution of them possesses the Gaussian form with the normalization constant being equal to the number of participants, which can be figured out in theory. This assumption is based on the consideration that, for a given incident energy, the leading particles have about the same energy and coincides with the fact that any kind of charged particles takes on well the Gaussian form of rapidity distribution. It is interesting to notice that the width of Gaussian rapidity distribution is irrelevant to the incident energy, centrality cut, and collision system. The fitting values of , the central positions of Gaussian rapidity distribution, are in good accordance with experimental data.
Comparing with the experimental measurements made by BRAHMS and PHOBOS Collaboration at BNLRHIC in AuAu and CuCu collisions at GeV and by ALICE Collaboration at CERNLHC in PbPb collisions at TeV, we can see that the total contributions from both evolutiondominated hydrodynamics and leading particles are well consistent with experimental data.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
This work is partly supported by the Transformation Project of Science and Technology of Shanghai Baoshan District with Grant no. CXY201225 and by the Shanghai Leading Academic Discipline Project with Grant no. XTKX 2012.
References
 J.Y. Ollitrault, “Anisotropy as a signature of transverse collective flow,” Physical Review D, vol. 46, no. 1, pp. 229–245, 1992. View at: Publisher Site  Google Scholar
 S. S. Adler, S. Afanasiev, C. Aidala et al., “Elliptic flow of identified hadrons in Au+Au collisions at $\sqrt{{s}_{\text{N}\text{N}}}=200$ GeV,” Physical Review Letters, vol. 91, Article ID 182301, 2003. View at: Publisher Site  Google Scholar
 K. Aamodt, B. I. Abelev, M. Reicher et al., “Higher harmonic anisotropic flow measurements of charged particles in PbPb collisions at $\sqrt{{s}_{\text{N}\text{N}}}=2.76$ TeV,” Physical Review Letters, vol. 107, Article ID 032301, 2011. View at: Publisher Site  Google Scholar
 P. A. Steinberg, “Bulk dynamics in heavy ion collisions,” Nuclear Physics A, vol. 752, pp. 423–432, 2005. View at: Publisher Site  Google Scholar
 C. Y. Wong, “Landau hydrodynamics reexamined,” Physical Review C, vol. 78, Article ID 054902, 2008. View at: Publisher Site  Google Scholar
 Z. J. Jiang, Q. G. Li, and H. L. Zhang, “The revised Landau hydrodynamic model and the pseudorapidity distributions of produced charged particles in high energy heavy ion collisions,” Journal of Physics G: Nuclear and Particle Physics, vol. 40, Article ID 025101, 2013. 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
 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
 E. K. G. Sarkisyan and A. S. Sakharov, “Relating multihadron production in hadronic and nuclear collisions,” European Physical Journal C, vol. 70, no. 3, pp. 533–541, 2010. View at: Publisher Site  Google Scholar
 H. Song, S. A. Bass, U. Heinz, T. Hirano, and C. Shen, “200 A GeV Au+Au collisions serve a nearly perfect quarkgluon liquid,” Physical Review Letters, vol. 106, no. 19, Article ID 192301, 2011. View at: Publisher Site  Google Scholar
 R. Rvblewski and W. Florkowski, “Highly anisotropic hydrodynamics in $3+1$ spacetime dimensions,” Physical Review C, vol. 85, Article ID 064901, 2012. 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
 I. M. Khalatnikov, “Some questions of the relativistic hydrodynamics,” Journal of Experimental and Theoretical Physics, vol. 27, pp. 529–541, 1954 (Russian). View at: Google Scholar
 L. D. Landau, “On the multiple production of particles in fast particle collisions,” Izvestiya Akademii Nauk SSSR, vol. 17, pp. 51–64, 1953 (Russian). View at: Google Scholar
 S. Z. Belen'kji and L. D. Landau, “Hydrodynamic theory of multiple production of particles,” Il Nuovo Cimento Series 10, vol. 3, no. 1, pp. 15–31, 1956. View at: Publisher Site  Google Scholar
 R. C. Hwa, “Statistical description of hadron constituents as a basis for the fluid model of highenergy collisions,” Physical Review D, vol. 10, no. 7, pp. 2260–2268, 1974. View at: Publisher Site  Google Scholar
 J. D. Bjorken, “Highly relativistic nucleusnucleus collisions: the central rapidity region,” Physical Review D, vol. 27, no. 1, pp. 140–151, 1983. View at: Publisher Site  Google Scholar
 T. Csörgő, N. I. Nagy, and M. Csanád, “New family of simple solutions of relativistic perfect fluid hydrodynamics,” Physics Letters B, vol. 663, pp. 306–311, 2008. View at: Publisher Site  Google Scholar
 M. I. Nagy, T. Csörgő, and M. Csanád, “Detailed description of accelerating, simple solutions of relativistic perfect fluid hydrodynamics,” Physical Review C, vol. 77, Article ID 024908, 2008. View at: Publisher Site  Google Scholar
 M. Csanád, M. I. Nagy, and T. Csörgő, “Similar final states from different initial states using new exact solutions of relativistic hydrodynamics,” European Physical Journal ST, vol. 155, pp. 19–26, 2008. 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
 M. S. Borshch and V. I. Zhdanov, “Exact solutions of the equations of relativistic hydrodynamics representing potential flows,” Symmetry, Integrability and Geometry: Methods and Applications, vol. 3, article 116, 11 pages, 2007. View at: Google Scholar  Zentralblatt MATH
 M. Csanád, M. I. Nagy, and S. Lökös, “Exact solutions of relativistic perfect fluid hydrodynamics for a QCD Equation of State,” European Physical Journal A, vol. 48, p. 173, 2012. View at: Publisher Site  Google Scholar
 M. Murray, “Scanning the phases of QCD with BRAHMS,” Journal of Physics G: Nuclear and Particle Physics, vol. 30, pp. S667–S674, 2004. View at: Publisher Site  Google Scholar
 M. Murray, “Flavor dynamics,” Journal of Physics G: Nuclear and Particle Physics, vol. 35, Article ID 044015, 2008. View at: Publisher Site  Google Scholar
 I. G. Bearden, D. Beavis, C. Besliu et al., “Charged meson rapidity distributions in central Au+Au collisions at $\sqrt{{s}_{\text{N}\text{N}}}=200$ GeV,” Physical Review Letters, vol. 94, Article ID 162301, 2005. View at: Publisher Site  Google Scholar
 B. Alver, B. Back, M. Baker et al., “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: Publisher Site  Google Scholar
 E. Abbas, B. Abelev, J. Adam et al., “Centrality dependence of the pseudorapidity density distribution for charged particles in PbPb collisions at $\sqrt{{s}_{\text{N}\text{N}}}=2.76$ TeV,” Physics Letters B, vol. 726, p. 610, 2013. View at: Publisher Site  Google Scholar
 A. Adare, S. Afanasiev, C. Aidala et al., “Scaling properties of azimuthal anisotropy in Au+Au and Cu+Cu collisions at $\sqrt{{s}_{\text{N}\text{N}}}=200$ GeV,” Physical Review Letters, vol. 98, p. 162301, 2007. View at: Publisher Site  Google Scholar
 N. Armesto, N. Borghini, S. Jeon et al., “Heavyion collisions at the LHCLast call for predictions,” Journal of Physics G: Nuclear and Particle Physics, vol. 35, Article ID 054001, 2008. 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,” European Physical Journal A, vol. 40, no. 1, pp. 99–108, 2009. View at: Publisher Site  Google Scholar
 S. Borsányi, G. Endrődi, Z. Fodor et al., “The QCD equation of state with dynamical quarks,” Journal of High Energy Physics, vol. 77, p. 1, 2010. View at: Google Scholar
 S. Amai, H. Fukuda, C. Iso, and M. Sato, “Hydrodynamical treatment of multiple meson production in high energy nucleonnucleus collisions,” Progress of Theoretical Physics, vol. 17, pp. 241–287, 1957. View at: Publisher Site  Google Scholar
 T. B. Li, The Mathematical Processing of Experiments, Science Press, Beijing, China, 1980, (Chinese).
 J. Voit, The Statistical Mechanics of Financial Markets, Springer, Berlin, Germany, 2005.
 Z.J. Jiang, “The numbers of participants and nucleonnucleon collisions in highenergy heavyion collisions,” Acta Physica Sinica, vol. 56, no. 9, pp. 5191–5195, 2007 (Chinese). View at: Google Scholar
 Z. J. Jiang, Y. F. Sun, and Q. G. Li, “The energy and centrality dependences of the pseudorapidity distributions of the charged particles in Au+Au collisions,” International Journal of Modern Physics E, vol. 21, Article ID 1250002, 2012. View at: Publisher Site  Google Scholar
 B. B. Back, M. D. Baker, M. Ballintijn et al., “The PHOBOS perspective on discoveries at RHIC,” Nuclear Physics A, vol. 757, pp. 28–101, 2005. View at: Publisher Site  Google Scholar
 K. Aamodt, A. Abrahantes Quintana, D. Adamová et al., “Centrality dependence of the chargedparticle multiplicity density at midrapidity in PbPb collisions at $\sqrt{{s}_{\text{N}\text{N}}}=2.76$ TeV,” Physical Review Letters, vol. 106, Article ID 032301, 2011. View at: Publisher Site  Google Scholar
 C. Y. Wong, Introduction To High Energy Heavy Ion Collisions, Press of Harbin Technology University, Harbin, China, 2002, (Chinese), English edition: World Scientific, Singapore, 1994.
 I. G. Bearden, D. Beavis, C. Besliu et al., “Nuclear stopping in Au+Au collisions at $\sqrt{{s}_{\text{N}\text{N}}}=200$ GeV,” Physical Review Letters, vol. 93, Article ID 102301, 2004. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2014 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}.