A Universal Description of Pseudorapidity Distributions in Both Nucleus-Nucleus and pp Collisions at Currently Available Energies

Investigations have shown that the collective motion appears not only in nucleus-nucleus but also in p-p collisions. The best tool for depicting such collective motion is relativistic hydrodynamics. In this paper, the collective motion is assumed to obey the hydrodynamicmodel which integrates the features of Landau andHwa-Bjorken theory and is one of a very few analytically solvable models. The fluid is then supposed to freeze out into charged particles from a space-like hypersurface with a fixed time of tFO. The investigations of present paper show that this part of charged particles together with leading particles, which, by conventional definition, carry on the quantum numbers of colliding nucleons and take away the most part of incident energy, can give a proper universal description to the pseudorapidity distributions of charged particles measured in both nucleus-nucleus and p-p collisions at currently available energies.


Introduction
In recent years, especially with the operations of BNL-RHIC and after CERN-LHC, the nature of matter created in nucleus or particle collisions has been undergoing extensive and deep experimental and theoretical investigations.One of the most important achievements arrived at from these investigations is that the matter formed in collisions shows the characteristics of collective motion as a nearly perfect fluid with a small shear viscosity over entropy density .
The best approach to describing the spatiotemporal evolution of the fluid-like matter is the relativistic hydrodynamics, which was first put forward by Landau in his pioneering work in 1953 [32].Owing to the high degree of nonlinearity and interconnection of hydrodynamic equations, numerical simulations are, as usual, widely employed to deal with them especially for 2-or 3-dimensional expansions or situations including viscosity.In numerical simulations, besides a powerful calculation system, there is also a need of a sophisticated skill for avoiding instabilities in solving partial differential hydrodynamic equations.
From the time of Landau till now, the analytical solutions of hydrodynamics are mainly limited to 1 + 1 expansions for a perfect fluid with a simple equation of state [8][9][10][11][12][13][14][15][16][17].There are also a few analytical discussions about 2 + 1 or 3 + 1 flows [18,19].Along with the discovery of collective motion of matter created in high energy physics in recent years, the analytical work has entered into a stage of rapid developments and achieved a number of good results.For example, by generalizing the relation between ordinary rapidity and space-time rapidity, [9] integrates Landau and Hwa-Bjorken [33,34] two models together, proposing a model called unified hydrodynamics in this paper and presenting a set of exact solutions.By taking advantage of the traditional scheme of Khalatnikov potential [35], [10] solved analytically the hydrodynamic equations and gave a pack of simple exact solutions for a perfect fluid with a linear equation of state.By taking into account the work done by the fluid elements on each other, [11][12][13] generalized the Hwa-Bjorken model for an accelerationless system to the model for an accelerated one and obtained a class of exact analytical solutions of relativistic hydrodynamics.

A Brief Description of the Model
The key ingredients of the model are as follows.
(1) The evolution of the fluid-like matter created in collisions is dominated by the energy-momentum conservation: where  ] = ( 0 ,  1 ) = (, ),  is the time and  is the longitudinal component of coordinates along beam direction.
] is the energy-momentum tensor, which, for a perfect fluid, takes the form as where  ] =  ] = diag(1, −1),   , , and  are the metric tensor, 4-velocity, energy density, and pressure of the fluid, respectively.For a constant speed of sound,  and  are related by the equation of state: where 1/ √  =  s is the speed of sound.Equations ( 2) and ( 3) make (1) translate to where  + and  − are the compact notation of partial derivatives with respect to light-cone coordinates  ± =  ±  =  0 ±  1 =  ± S ,  = √ +  − is the proper time,  S = 1/2 ln( + / − ) is the space-time rapidity, and  is the ordinary rapidity of the fluid.
(2) Equation ( 4) is a set of complicated differential equations with high nonlinearity and coupling between variables  and .In order to solve it, the relation between  and  S is generalized to the form [9]: where  ± =  0 ±  1 =  ± are the light-cone components of 4-volicity, and  ± ( ± ) are the two arbitrary functions.In case of (5) reduces to  =  S , returning to the boost-invariant picture of Hwa-Bjorken.Otherwise, (5) describes the non-boostinvariant geometry of Landau.Accordingly, (5) unifies the Landau with Hwa-Bjorken hydrodynamics.It plays a role of bridge between these two models.By using ( 5), (4) becomes where the prime stands for the derivatives with respect to  + or  − ,  ± =  ± /, and  is an arbitrary constant.Completing the integration of the above equation, it becomes  ln  = − ( + 1) In order to keep the consistency of the above two equations, Δ ± ( ± ) should take Inserting them into (8), we can obtain the final solution [9]: Advances in High Energy Physics where  is the entropy density of the fluid and  0 is its initial scale., as mentioned above, takes as a constant in this paper, and ℎ = /, and  is an arbitrary constant fitting condition: where the double prime represents the second derivatives with respect to  + or  − .As  = 0, the above equation reduces to   ± = 0 which gives the solution of  ± ( ± ) =  ± + , returning to the Hwa-Bjorken case of ( 6) up to a constant .Hence, the values of  or ℎ = / describe the deviation of system from Hwa-Bjorken expansions.
(3) As the expansion of the fluid lasts to the time of  FO , the inelastic interactions between the particles in the fluid cease, and the ratios of different kinds of particles remain unchanged.At this moment, the collective motion of the fluid comes to an end, and the fluid decouples or freezes out into the detected particles from a space-like hypersurface with a fixed time of  FO .The rapidity distributions of entropy of the fluid at this hypersurface take the form: where   is the 4-dimensional unit vector normal to the hypersurface: d  = d  , and d is the space-like infinitesimal length element along hypersurface: where the minus sign accounts for the space-like characteristic of d.Considering that the number of charged particles is proportional to the amount of entropy; from solution ( 10) and ( 13), we can obtain the rapidity distributions of the produced charged particles [9]: where  is an overall normalization constant. stands for an arbitrary space-like hypersurface.
(4) The right-hand side of ( 16) is evaluated on the space-like hypersurface with the time equaling  FO .Such hypersurface can be therefore taken as where  is an arbitrary constant.This equation gives Thus, ( 16) turns into (5) In nucleus or particle collisions, apart from the freezeout of the fluid, leading particles also have certain contribution to the measured charged particles.Leading particles are believed to be formed outside the nucleus, that is, outside the colliding region [42,43].The motion and generation of leading particles are therefore free from hydrodynamic descriptions.As we have argued before that the rapidity distribution of leading particles takes the Gaussian form [15,16], where  Lead ,  0 , and  are the number of leading particles, central position, and width of distribution, respectively.This conclusion comes from the consideration that, for a given incident energy, different leading particles resulting from each time of collisions have approximately the same amount of energy or rapidity.Then, the central limit theorem [44,45] guarantees the plausibility of the above argument.Actually, experimental measurements have shown that any kind of charged particles presents a good Gaussian rapidity distribution [46][47][48]. 0 in (20) is the average position of leading particles.It should increase with incident energies and centrality cuts.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 energies, centrality cuts, and even colliding systems.The specific values of  0 and  can be determined by tuning the theoretical predictions to experimental data.
By definition, leading particles mean the particles inheriting the quantum numbers of colliding nucleons.Hence, the number of leading particles is equal to that of participants.For p-p collisions, there are only two leading particles.They are separately in the projectile and target fragmentation region.For an identical nucleus-nucleus collision, the number of leading particles where  Part is the total number of participants, which can be determined in theory by Glauber model [49,50].[36].The dashed curves are the results from unified hydrodynamics of (19).The dashed-dotted curves are the results from leading particles of (20).The solid curves are the sums of dashed and dashed-dotted curves.

Comparison with Experimental Measurements in Nucleus-Nucleus Collisions
Having rapidity distributions of ( 19) and ( 20), the pseudorapidity distribution measured in experiments can be expressed as where  T is the transverse momentum,  T = √ 2 +  2 T is the transverse mass, and is the total rapidity distribution from both the freeze-out of the fluid and leading particles.Substituting (24) or ( 19) and ( 20) into (22), we can get the pseudorapidity distributions of the charged particles.The results are shown in Figures 1 and 2, which are for different centrality Au-Au and Pb-Pb collisions at √ NN = 200 GeV and 2.76 TeV, respectively.The solid dots in the figures are the experimental measurements [36,37].The dashed curves are the results got from unified hydrodynamics of (19).The dashed-dotted curves are the results obtained from leading The solid dots are the experimental measurements [37].The dashed curves are the results from unified hydrodynamics of (19).The dashed-dotted curves are the results from leading particles of (20).The solid curves are the sums of dashed and dashed-dotted curves.
particles of (20).The solid curves are the sums of dashed and dashed-dotted curves.It can be seen that the combined contributions from both unified hydrodynamics and leading particles match up well with experimental data.
Experiments have shown that the overwhelming majority of charged particles produced in heavy ion collisions at high energy consists of pions, kaons, and protons with proportions of about 83%, 12%, and 5%, respectively [51].These proportions are not evidently dependent on colliding energies and systems.Furthermore, for a given incident energy, the transverse momentum  T changes very slowly with centrality cut.For a specific type of charged particle, it can be well taken to be a constant for centrality cuts from 0-55%, which we are interested in.In Au-Au collisions at √ NN = 200 GeV, this constant is about 0.45, 0.65, and 0.93 GeV/ for pions, kaons, and protons, respectively.In calculations, the mass and transverse momentum in ( 22) and ( 23) take the values of  = 0.22 GeV and  T = 0.50 GeV/ in Au-Au collisions at √ NN = 200 GeV, which are the mean values of those of pions, kaons, and protons.In Pb-Pb collisions at √ NN = 2.76 TeV, they take the values of  = 0.22 GeV and  T = 0.62 GeV/.Here,  T takes the mean value from experimental measurements [4].
The factor  in (3) takes the values of 6.6 and 4.5, respectively, for different centrality Au-Au and Pb-Pb collisions.That is, it decreases with increasing energies.It should change like this.Since  decreases with increasing temperature of the fluid [52], and the higher the energy, the higher the temperature of the created fluid.
The normalization constant  in (19) takes the values of (6.45, 6.07, 5.31, 4.56, 3.84, 3.12, 2.57, 2.15, 1.71, 1.39, 1.06) × 10 3 and (1.74, 1.44, 1.08, 0.74) × 10 4 for centrality cuts from smaller to larger in Au-Au and Pb-Pb collisions, respectively.The parameter ℎ in (11) takes the values of (5.41, 3.08, 2.53, 2.05, 1.01, 0.92, 0.83, 0.55, 0.47, 0.36, 0.22) × 10 −4 and (1.02, 0.78, 0.51, 0.29) × 10 −5 for centrality cuts from smaller to larger in Au-Au and Pb-Pb collisions, respectively.It can be seen that ℎ decreases with increasing energies and centrality cuts.It should vary in this way.Since the larger the energy and centrality cut, the more transparent the nucleus.The produced charged particles will then be located in a wider rapidity region.While, the region of rapidity distributions Advances in High Energy Physics of ( 19) is mainly determined by ℎ.The wider the region of rapidity distributions, the smaller the value of ℎ.This can be seen in Figure 3, which is for Au-Au collisions normalized to 0-3% centrality cut.The wider curve is the corresponding rapidity distribution in Figure 1, obtained by using ℎ = 5.41 × 10 −4 .The narrower one is for ℎ = 5.41 × 10 −2 .Hence, the region of rapidity distributions decreases with increasing ℎ.This forms the basis for the above fitted ℎ values.The central parameter  0 in (20) takes the values of 2.91, 2.93, 2.95, 2.98, 3.01, 3.05, 3.10, 3.17, 3.20, 3.25, and 3.30 and 3.61, 3.62, 3.64, and 3.65 for centrality cuts from smaller to larger in Au-Au and Pb-Pb collisions, respectively.As mentioned above, it increases with centrality cuts and energies.The width parameter  in (20) takes a constant of  = 0.85 for different centrality cuts in both Au-Au and Pb-Pb collisions.As stated above, it is irrelevant to centrality cuts, incident energies, and colliding systems.

Comparison with Experimental Measurements in p-p Collisions
Compared with nucleus-nucleus collisions, p-p collisions are a relatively simpler and smaller system and therefore produce the matter in a much tinier volume.If such matter expands in the same way as that generated in nucleus-nucleus collisions has been experiencing endless disputations.Such situation has changed since a series of findings in recent years [20][21][22][23][24][25][26][27][28][29][30][31].The investigations of [20] have shown that the hydrodynamic predictions of multiplicity productions in p-p (p) or e + -e − collisions at even several of GeV match up well with experimental data.The measurements from CMS Collaboration at CERN-LHC have shown that [21], just like in nucleus-nucleus collisions, the ridge structures, the signal of collective motion, also appear in pp collisions.References [22,23] further substantiated that these observed ridge structures can be well understood in the framework of hydrodynamics.In [24], the STAR Collaboration at BNL-RHIC has measured the HBT (Hanbury-Brown-Twiss) radii for p-p collisions as a function of multiplicity.The measured results are well favored by the predictions of hydrodynamic models made in [25][26][27].
In [28,29], Sarkisyan et al. utilized the combined model of Landau hydrodynamics plus combinations of the constituent quarks in participants to deal with the multihadron productions and transverse momentum distributions in pp (p) collisions.By employing the hydrodynamic solutions known as Gubser flow [18], [30] demonstrated the collective radial flow of high multiplicity p-p collisions.In our previous work [31], we used the evolution-dominated hydrodynamics together with the effects of leading particles to describe the pseudorapidity distributions of charged particles produced in p-p collisions.At present paper, from (19) and (20), we can also get such distributions.The results are shown in Figure 4 for p-p collisions at energies from 23.6 to 900 GeV.The solid dots are the experimental measurements [38][39][40].The dashed curves are the results from unified hydrodynamics of (19).The dashed-dotted curves are the results from leading particles of (20).The solid curves are the sums of dashed and dasheddotted curves.It can be seen that the theoretical results are in well consistent with experimental data.
In calculations, as in Au-Au or Pb-Pb collisions, the mass  in ( 22) and ( 23) takes the value of 0.22 GeV.The transverse momentum  T varies with energies according to relation [41]: where  T and √ are in units of GeV/ and GeV, respectively.The width parameter  in (20) takes the same value as that in nucleus-nucleus collisions; that is,  = 0.85.This certifies the above argument again that  is independent of energies and colliding objects.The factors , ℎ, and  0 take the values of 8.4-5.2,18.4-4.0× 10 −5 , and 1.91-3.00,respectively, for energies from lower to higher.Just as in nucleus-nucleus collisions,  and ℎ decrease with increasing energies, and,  The pseudorapidity distributions of produced charged particles in p-p collisions at energies from √ = 23.6 to 900 GeV.The solid dots are the experimental measurements [38][39][40].The dashed curves are the results from unified hydrodynamics of (19).The dashed-dotted curves are the results from leading particles of (20).The solid curves are the sums of dashed and dashed-dotted curves.
0 , on the contrary, increases with energies.Their variations against energies are presented in Figure 5.The solid dots represent the above fitted values.The solid curves are drawn from relations  = 11.3311− 0.5324 ln  + 0.0065ln 2 , ln ℎ = 23.3049− 3.7056 ln  + 0.0877ln 2 ,  0 = 0.3657 + 0.2709 ln  − 0.0055ln 2 , (26) where √ is in units of GeV.It can be seen that all the solid dots are well seated on the curves.The stars and circles in this figure are the corresponding predictions for p-p collisions at CERN-LHC energies of 2.36 and 7 TeV, respectively.From these predictions, we can get the pseudorapidity distributions of produced charged particles in these two situations, and the results are shown in Figure 6.The solid dots are the experimental data [41].The dashed curves are the predicted results from unified hydrodynamics of (19).The dasheddotted curves are the predicted results from leading particles of (20).The solid curves are the sums of dashed and dashed-dotted curves.It can be seen that the theoretical predictions agree well with the available measurements in the midpseudorapidity regions.
At the end of this section, we would like to point out that, at present, we have no knowledge about the viscosities in matter created in TeV energies especially for small system of p-p collisions.Here, just as what we have done in case of nucleus-nucleus collisions at RHIC energies, we ignore the viscous effects for the purpose of simplifications.

Conclusions
The charged particles in nucleus and particle collisions are supposed to have the same producing mechanism.That is, TeV.The solid dots are the experimental measurements [41].The dashed curves are the predicted results from unified hydrodynamics of (19).The dashed-dotted curves are the predicted results from leading particles of (20).The solid curves are the sums of dashed and dashed-dotted curves.

Figure 1 :
Figure1: The pseudorapidity distributions of produced charged particles in different centrality Au-Au collisions at √ NN = 200 GeV.The solid dots are the experimental measurements[36].The dashed curves are the results from unified hydrodynamics of(19).The dashed-dotted curves are the results from leading particles of(20).The solid curves are the sums of dashed and dashed-dotted curves.

Figure 2 :
Figure 2: The pseudorapidity distributions of produced charged particles in different centrality Pb-Pb collisions at √ NN = 2.76 TeV.The solid dots are the experimental measurements[37].The dashed curves are the results from unified hydrodynamics of(19).The dashed-dotted curves are the results from leading particles of(20).The solid curves are the sums of dashed and dashed-dotted curves.

Figure 3 :
Figure 3: The rapidity distributions of produced charged particles for two different ℎ values in Au-Au collisions normalized to 0-3% centrality cut.

Figure 4 :
Figure4: The pseudorapidity distributions of produced charged particles in p-p collisions at energies from √ = 23.6 to 900 GeV.The solid dots are the experimental measurements[38][39][40].The dashed curves are the results from unified hydrodynamics of(19).The dashed-dotted curves are the results from leading particles of(20).The solid curves are the sums of dashed and dashed-dotted curves.

Figure 6 :
Figure6: The pseudorapidity distributions of produced charged particles in p-p collisions at energies of √ = 2.36 and 7 TeV.The solid dots are the experimental measurements[41].The dashed curves are the predicted results from unified hydrodynamics of(19).The dashed-dotted curves are the predicted results from leading particles of(20).The solid curves are the sums of dashed and dashed-dotted curves.