# The Unified Hydrodynamics and the Pseudorapidity Distributions in Heavy Ion Collisions at BNL-RHIC and CERN-LHC Energies

Authors: Z. J. Jiang; J. Wang; K. Ma; H. L. Zhang
Journal: Advances in High Energy Physics
 (2015)

---

## Abstract

The charged particles produced in nucleus-nucleus collisions are divided into two parts. One is from the hot and dense matter created in collisions. The other is from leading particles. The hot and dense matter is assumed to expand according to unified hydrodynamics and freezes out into charged particles from a space-like hypersurface with a fixed proper time ofτFO. The leading particles are conventionally taken as the particles which inherit the quantum numbers of colliding nucleons and carry off most of incident energy. The rapidity distributions of the charged particles from these two parts are formulated analytically, and a comparison is made between the theoretical results and the experimental measurements performed in Au-Au and Pb-Pb collisions at the respective BNL-RHIC and CERN-LHC energies. The theoretical results are well consistent with experimental data.


---

## Body

## 1. Introduction

Since the elementary work of Landau in 1953 [1], relativistic hydrodynamics has been applied to calculate a large number of variables developed in the context of nucleon or nucleus collisions at high energy. In particular, owing to the successful description of elliptic flow measured at the Relativistic Heavy-Ion Collider (RHIC) at Brookhaven National Laboratory (BNL) [2] and recently at the Large Hadron Collider (LHC) at CERN [3], the hydrodynamic research has entered into a more active phase. It has now been widely accepted as one of the best approaches for understanding the space-time evolution of the matter created in such collisions [4–10].Although at present there are powerful numerical approaches to deal with certain hydrodynamic problems, this will require a very large scale of calculation and skillful sophisticated techniques to avoid instabilities. On the contrary, analytical solutions, since their simple and transparent forms, usually providing us with an invaluable insight into the characteristics of matter created in collisions, are always our pursuit of the goal though usually at the price of some ideal assumptions.Due to the tremendous complexity of hydrodynamic equations, the progress in finding exact solutions is not going well. Up till now, most of this work mainly involves1+1 dimensional flows with simple equation of state [11–19]. The 3+1 dimensional hydrodynamics is less developed, and no general exact solutions are known so far.An important application of1+1 dimensional hydrodynamics is the analysis of the pseudorapidity distributions of the charged particles produced in nucleon or nucleus collisions. In this paper, by taking into account the effect of leading particles, we will discuss such distributions in the context of unified hydrodynamics [14]. The main points of this model are listed in Section 2. Its solution is then exploited in Section 3 to formulate the rapidity distributions of the charged particles frozen out from fluid at a space-like hypersurface with a fixed proper time τFO. In Section 4, the theoretical results are compared with the experimental data performed in nucleus-nucleus collisions at BNL-RHIC and CERN-LHC energies. Section 5 is about conclusions.


## 2. A Brief Description of the Unified Hydrodynamics

The1+1 expansion of a perfect fluid follows the equations:(1)∂T00∂t+∂T10∂z=0,∂T01∂t+∂T11∂z=0,where t is the time, z is the longitudinal coordinate along beam direction, and(2)Tμν=ε+puμuν-pgμνis the energy-momentum tensor, and gμν=diag1,-1,-1,-1,  uμ,  ε, and p are, respectively, the metric tensor, 4-velocity, energy density, and pressure of fluid. For a constant speed of sound, ε and p are related by the equation of state:(3)ε=gp,where 1/g=cs is the speed of sound. Investigations have shown that g changes very slowly with energy and centrality [20–23]. For a given incident energy, it can be well taken as a constant.Using (2) and (3), and noticing the light-cone components of the 4-velocity(4)u±=e±y,where y is the ordinary rapidity of fluid, (1) reads(5)e2y-12g+1∂+p+e2yg+1p∂+y+1-e-2y2g+1∂-p+e-2yg+1p∂-y+∂+p-∂-p=0,e2y+12g+1∂+p+e2yg+1p∂+y+1+e-2y2g+1∂-p-e-2yg+1p∂-y-∂+p-∂-p=0,where ∂+ and ∂- are the compact notation of partial derivatives with respect to light-cone coordinates z±=t±z=x0±x1=τe±ηST, τ=z+z- is the proper time, and ηST=1/2lnz+/z- is the space-time rapidity of fluid. The solutions of above equations are(6)g∂+lnp=-g+122∂+y-g2-12e-2y∂-y,g∂-lnp=g+122∂-y+g2-12e2y∂+y.The key ingredient of unified hydrodynamics is that it generalizes the relation betweeny and ηST by (7)2y=lnu+-lnu-=lnF+z+-lnF-z-,where F±z± are a priori arbitrary function obeying equation(8)F±F±′′=A22,where A is a constant. In case of(9)F±z±=z±.Equation (7) reduces to y=ηST, returning to the boost-invariant picture of Hwa-Bjorken. Otherwise, (7) describes the nonboost-invariant geometry of Landau. Accordingly, (7) unifies the Hwa-Bjorken and Landau hydrodynamics together. It paves a way between these two models.Substituting (7) into (6), we have(10)g∂+lnp=-g+124f+′f++g2-14f-′f+,g∂-lnp=-g+124f-′f-+g2-14f+′f-,where f±=F±/H and H is an arbitrary constant. From the above equations, we can get the entropy density of fluid:(11)sz+,z-=s0pp0g/g+1=s0exp-g+14l+2+l-2+g-12l+l-,where, by definition(12)l±z±=lnf±,and in terms of l±, we have from (7) and (8)(13)yz+,z-=12l+2-l-2,(14)z±=2h∫0l±eu2du,where h=H/A.


## 3. The Rapidity Distributions of the Charged Particles Frozen out from Fluid

By using (11), we can obtain the rapidity distributions of the charged particles frozen out from fluid or hot and dense matter created in collisions. To this end, we first evaluate the entropy distributions of the fluid on a space-like hypersurface with a fixed proper time τFO, from which the fluid will freeze out into the charged particles. Such distributions take the form(15)dSdy=suμdλμdyτFO=suμnμdλdyτFO,where nμ is the 4-dimensional unit vector normal to the hypersurface(16)nμnμ=n+n-=1.Also dλμ=dλnμ, and dλ is the space-like infinitesimal length element along hypersurface(17)dλ=dλμdλμ=-dz+dz-.Considering a hypersurface(18)ϕz+,z-=C,where C is a constant, we have(19)∂+ϕdz++∂-ϕdz-=0.The light-cone components of unit vector nμ can be expressed by ∂±ϕ as(20)n+=∂-ϕ∂+ϕ∂-ϕ,n-=∂+ϕ∂+ϕ∂-ϕ,dλ=-dz+dz-=dz-∂-ϕ∂+ϕ=dz+∂+ϕ∂-ϕ.Then, the expression on the right-hand side of (15) is(21)uμnμdλ=12u+n-+u-n+dλ=12u+∂+ϕ+u-∂-ϕdz-∂+ϕ=12u+∂+ϕ+u-∂-ϕdz+∂-ϕ.Furthermore, known from (14),(22)dz±=hl±expl±2dl±2.In terms of ∂±ϕ and l±, dl±2 can be written as(23)dl±2=±2∂∓ϕexpl∓2/l∓∂+ϕexpl+2/l++∂-ϕexpl-2/l-dy.Thus, (22) becomes(24)dz±=±2h∂∓ϕ∂+ϕl-exp-l-2+∂-ϕl+exp-l+2dy.Substituting it into (21) and then into (15), we have(25)dSdy=se1/2l+2+l-2∂+ϕey+∂-ϕe-y∂+ϕl-ey+∂-ϕl+e-yτFO.In the above equation, the right-hand side is evaluated at the hypersurface with proper timeτFO. Known from (18), this hypersurface can be taken as(26)ϕz+,z-=τFO2=z+z-=C,z+dz-+z-dz+=0.Comparing the above equation with (19), we get(27)∂±ϕ=z∓.Inserting it together with (11) into (25) and noticing the proportional relation between entropy and the number of the charged particles, we have(28)dNFluidb,sNN,ydy=Cb,sNNe-g-1l+-l-2/4z-ey+z+e-yz-l-ey+z+l+e-y,where Cb,sNN, independent of rapidity y, is an overall normalization constant. b is the impact parameter, and sNN is the center-of-mass energy per pair of nucleons.


## 4. Comparison with Experimental Measurements and the Rapidity Distributions of Leading Particles

From (28), we can get the rapidity distributions of the charged particles generated in collisions. Figure 1 presents such distributions for π+, π-, K+, K-, p, and p- resulting from 5% most central Au-Au collisions at sNN=200 GeV. The symbols are the experimental measurements given by BRAHMS Collaboration at BNL-RHIC [24–26]. The solid curves are the theoretical predictions from (28). From this figure, we can see that (28) coincides well with the experimental data of all the charged particles with the exception of p. For proton p, an obvious discrepancy appears in the rapidity interval of y=2.0~3.0. This might be caused by the effects of leading particles.Figure 1
The rapidity distributions of the specified charged particles in central Au-Au collisions atsNN=200 GeV. The scattered symbols are the experimental measurements [24–26]. The solid curves are the results from unified hydrodynamics of (28).The leading particles are believed to be formed outside the nucleus, that is, outside the colliding region [27, 28]. The generation of leading particles is therefore free from fluid evolution. Accordingly, their rapidity distributions are beyond the scope of hydrodynamics and should be treated separately.In our previous work [6], we once argued that the rapidity distribution of leading particles takes the Gaussian form:(29)dNLeadb,sNN,ydy=NLeadb,sNN2πσexp-y-y0b,sNN22σ2,where NLeadb,sNN, y0b,sNN, and σ are, respectively, the number of leading particles, central position, and width of distribution. This conclusion comes from the consideration that, for a given incident energy, different leading particles resulting from each time of nucleus-nucleus collisions have approximately the same amount of energy or rapidity. Then, the central limit theorem [29, 30] guarantees the rationality of above argument. Actually, as one can see from the shapes of the curves in Figure 1 any kind of charged particles resulting in collisions forms a good Gaussian rapidity distribution.Termy0b,sNN in (29) is the average position of leading particles. It should increase with incident energy and centrality. The value of σ relies on the relative energy or rapidity differences among leading particles. It should not, at least not apparently, depend on the incident energy, centrality, and even colliding system. The concrete values of y0 and σ can be determined by tuning the theoretical predictions to experimental data.By definition, leading particles mean the particles which carry on the quantum numbers of colliding nucleons and take away most of incident energy. Then, the number of leading particles is equal to that of participants. For nucleon-nucleon, such asp-p collisions, there are only two leading particles. They are separately in projectile and target fragmentation region. For an identical nucleus-nucleus collision, the number of leading particles is (30)NLeadb,sNN=NPartb,sNN2,where NPartb,sNN is the total number of participants equaling [31, 32](31)NPartb,sNN=∫nPartb,sNN,sd2s,where the integral variable s is the transverse coordinates in the almond-shaped colliding region with respect to the center of one nucleus. nPartb,sNN,s is the total number of participants in the flux tube with a unit bottom area, located at position s along beam direction. It takes the form(32)nPartb,sNN,s=TAs1-exp-σNNinsNNTBs-b+TBs-b1-exp-σNNinsNNTAs,where σNNinsNN, the inelastic nucleon-nucleon cross section, increases slowly with energy. For instance, for sNN=200, 130, 62.4 GeV and 2.76 TeV, σNNin=42, 41, 36, and 64±5 mb [33, 34], respectively. The subscripts A and B in the above equation represent the projectile and target nucleus. TAs or TBs-b is the nuclear thickness function with the value equaling the nucleon number in the flux as defined above. It is equal to(33)Ts=∫ρs,zdz,where(34)ρr=ρ01+expr-r0/ais the well-known Woods-Saxon distribution of nuclear density. a and r0, taking somewhat different values in different papers [31], are, respectively, the skin depth and radius of nucleus. In this paper, we take the value of a=0.54 fm and r0=1.12A1/3-0.86A-1/3 fm, where A is the mass number of nucleus. Here, ρ0 is another constant determined by condition(35)∫ρrdV=A.Known from the investigations in [6], (31) can give a correct result in different nucleus-nucleus collisions at different centrality and energy from BNL-RHIC to CERN-LHC scales.From rapidity distributions, we can get pseudorapidity distributions by relation [35](36)dNb,sNN,ηdη=1-m2mT2cosh2ydNb,sNN,ydy,where mT=m2+pT2 is the transverse mass and pT is the transverse momentum. The first factor on the right-hand side of the above equation is actually the Jacobian determinant. This transformation is closed by another relation:(37)y=12lnpT2cosh2η+m2+pTsinhηpT2cosh2η+m2-pTsinhη.Taking into account the contributions from both the freeze-out of fluid and leading particles, the rapidity distributions in (36) can be written as(38)dNb,sNN,ydy=dNFluidb,sNN,ydy+dNLeadb,sNN,ydy.Inserting the above equation or the sum of (28) and (29) into (36), we can get the pseudorapidity distributions of the charged particles. The results are shown in Figures 2, 3, 4, and 5, which are, respectively, for distributions in different centrality of Au-Au collisions at sNN=200, 130, and 62.4 GeV and in different centrality of Pb-Pb collisions at sNN=2.76 TeV. The solid dots are the experimental data [36, 37]. The dashed curves are the results obtained from unified hydrodynamics of (28). The dashed-dotted curves are the results got from leading particles of (29). The solid curves are the results acquired from (38), that is, the sums of the dashed and dashed-dotted curves. It can be seen that the theoretical results are in good agreement with experimental measurements.Figure 2
The pseudorapidity distributions of the charged particles produced in different centrality of Au-Au collisions atsNN=200 GeV. The solid dots are the experimental measurements [36]. The dashed curves are the results from unified hydrodynamics of (28). The dashed-dotted curves are the results from leading particles of (29). The solid curves are the sums of the dashed and dashed-dotted ones.Figure 3
The pseudorapidity distributions of the charged particles produced in different centrality of Au-Au collisions atsNN=130 GeV. The solid dots are the experimental measurements [36]. The dashed curves are the results from unified hydrodynamics of (28). The dashed-dotted curves are the results from leading particles of (29). The solid curves are the sums of the dashed and dashed-dotted ones.Figure 4
The pseudorapidity distributions of the charged particles produced in different centrality of Au-Au collisions atsNN=62.4 GeV. The solid dots are the experimental measurements [36]. The dashed curves are the results from unified hydrodynamics of (28). The dashed-dotted curves are the results from leading particles of (29). The solid curves are the sums of the dashed and dashed-dotted ones.Figure 5
The pseudorapidity distributions of charged particles produced in different centrality of Pb-Pb collisions atsNN=2.76 TeV. The solid dots are the experimental measurements [37]. The dashed curves are the results from unified hydrodynamics of (28). The dashed-dotted curves are the results from leading particles of (29). The solid curves are the sums of the dashed and dashed-dotted ones.In calculations, the parameterg in (28) takes the values of g=8.16 and 5.55 in the respective Au-Au and Pb-Pb collisions at all of the concerned incident energy and centrality [20]. The width parameter σ in (29) takes the value of 0.85 and 0.90 in Au-Au and Pb-Pb collisions, respectively, at different incident energy and centrality. As the analyses given above, σ is independent of incident energy and centrality and also not evidently dependent on colliding system. The center parameter y0 in (29) takes the values of 2.64–2.89, 2.63–2.83, 2.61–2.82, and 3.40–3.97 for centrality from small to large in Au-Au collisions at sNN=200, 130, and 62.4 GeV and in Pb-Pb collisions at sNN=2.76 TeV, respectively. As pointed out earlier, for a given colliding system, y0 increases with energy and centrality.


## 5. Conclusions

By generalizing the relation between ordinary rapidityy and space-time rapidity ηST, unified hydrodynamics integrates the features of Hwa-Bjorken and Landau two famous hydrodynamic models together. In case of linear equation of state, this hydrodynamic model can be solved analytically. The exact solutions can then be used to formulate the rapidity distributions of the charged particles frozen out from the fluid at the space-like hypersurface with a fixed proper time τFO. The unified hydrodynamics is successful in describing the rapidity distributions of π±, K±, and p-. However, it is not as good in describing the corresponding distributions of protons which might contain the leading particles.In our previous work [6], we once discussed the same distributions in the context of evolution-dominated hydrodynamics. This hydrodynamic model differs from the one used in this paper in two ways. (1) Different initial conditions: the former assumes that the fluid is initially at rest. Its motion is totally dominated by the following evolution. The latter, as mentioned above, employs the initial condition of (7). It plays a connection between Hwa-Bjorken and Landau hydrodynamics. (2) Different freeze-out conditions: the former assumes that the freeze-out of fluid occurs at the space-like hypersurface with a fixed temperature TFO. The latter, as stated above, takes such hypersurface as one with a fixed proper time τFO.As for leading particles, we argue that they possess the Gaussian rapidity distribution normalized to the number of participants. This is the same as that proposed in [6]. Here, for the purpose of completion and applications, we list out the most points of this proposition. It is interesting to notice that the investigations of the present paper once again show that, for a given colliding system, the central position y0 of Gaussian rapidity distribution increases with incident energy and centrality, while the width σ of the distribution is irrelevant to them and also almost independent of colliding system. These are consistent with the conclusions arrived at in [6].Comparing with experimental data carried out by BRAHMS and PHOBOS Collaboration at BNL-RHIC in different centrality of Au-Au collisions atsNN=200, 130, and 62.4 GeV and by ALICE Collaboration at CERN-LHC in different centrality of Pb-Pb collisions at sNN=2.76 TeV, we can see that, although the charged particles frozen out from fluid play a dominant role, leading particles are also essential in characterizing the measured distributions. Only after the total contributions from both unified hydrodynamics and leading particles are taken into account together can the experimental measurements be matched up properly.



---

*Source: 430606-2015-06-10.xml*