Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
Next Article in Journal
String Interactions as a Source of Collective Behaviour
Next Article in Special Issue
Conserved Charge Fluctuations from RHIC BES and FXT
Previous Article in Journal
Reconstructing Modified and Alternative Theories of Gravity
Previous Article in Special Issue
Polarization Measurements of p and 3He Beams at RHIC and Future EIC Using the Polarized Atomic Hydrogen Gas Jet Target
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Scaling Behaviour of dN/dy in High-Energy Collisions

1
HUN-REN Wigner RCP, P.O. Box 49, H-1525 Budapest, Hungary
2
MATE Institute of Technology KRC, H-3200 Gyöngyös, Hungary
*
Author to whom correspondence should be addressed.
Universe 2024, 10(1), 45; https://doi.org/10.3390/universe10010045
Submission received: 30 November 2023 / Revised: 9 January 2024 / Accepted: 12 January 2024 / Published: 17 January 2024
(This article belongs to the Special Issue Multiparticle Dynamics)

Abstract

:
From a recently found family of analytic, finite and accelerating 1+1-dimensional solutions to perfect fluid relativistic hydrodynamics, we derive simple and powerful formulae to describe the rapidity and pseudorapidity density distributions. By introducing a new scaling function, we notice that the rapidity distribution data of the different experiments all collapse into a single curve. This data-collapsing (or -scaling) behaviour in the rapidity distributions suggests that high-energy p + p collisions may be described as collective systems.

1. Introduction

It is experimentally known that, when very fast particles collide, a large number of new particles are created. The study of the process of high-energy collisions based on thermodynamic methods was first proposed by Fermi [1,2]. However, Fermi’s theory proved to be incomplete, which, in short, was due to his incorrect view of the expansion of the complex system [3]. As Landau has shown, the expansion of the system can be studied using relativistic hydrodynamics [4].
In ref. [4], Landau has discussed in detail that, at the moment that two nucleons collide, a large number of particles are created and the mean free path of the collisions between the particles is small compared to the size of the system; thus, the medium reaches statistical equilibrium. The expansion of the medium should be studied in a hydrodynamic framework as a perfect fluid, since the mean free path during the expansion process remains small compared to the size of the system. Due to the presence of huge energy density, particles are continuously created and absorbed in the medium, so the system can hardly be characterised by the number of particles [3]. As the system expands, the interaction between particles decreases, so the mean free path increases. When the interaction becomes sufficiently small, the number of particles can be considered a physical property. When the mean free path becomes comparable to the size of the system, the latter disintegrates into separate particles [3].
The Hwa–Bjorken solution to relativistic hydrodynamics is suitable for describing 1+1-dimensional explosive perfect fluids [5,6]. The solution utilizes a Hubble-type velocity field and one of its main shortcomings is that it assumes zero relativistic acceleration. Consequently, it predicts a rapidity plateau, or a flat rapidity distribution, which does not describe any of the high-energy experimental data at RHIC or LHC energies (except, perhaps, in a narrow interval near midrapidity, where the distribution is approximately constant near its maximum). However, one of the great advantages of this Hwa–Bjorken solution is that it uses simple formulae, and, from this solution, Bjorken derived his famous formula for estimating the initial energy density of the hot and dense matter created in heavy-ion collisions [6]. This solution has previously been generalised to the case of an accelerating velocity field [7,8], but these generalisations only satisfy the equations of hydrodynamics under strong constraints. In one of these new exact solutions [7,8], the rate of acceleration may have an arbitrary constant value, but only for a super-hard, non-realistic equation of state, where the speed of sound and the speed of light are equal. Recently, Csörgő, Kasza, Csanád and Jiang published a 1+1-dimensional generalisation of the Hwa–Bjorken solution, where the rate of acceleration of the velocity field may have an arbitrary constant value, even for a realistic equation of state [9], referred to as the CKCJ solution.
This CKCJ family of perfect fluid solutions and its applications have been published previously. The calculation of the pseudorapidity distribution and the fitting of the resulting formula to experimental data were published in refs. [9,10,11]. The correction to Bjorken’s initial energy density estimate was presented in detail in ref. [12]. Based on measured data, we have determined the lifetime parameter of heavy-ion collisions for s N N = 130 and 200 GeV A u + A u reactions and published the result in refs. [11,13]. We calculated the initial energy density in s N N = 130 and 200 GeV A u + A u collisions, which were also published in [11]. The most recent result is the evaluation of the thermal photon spectrum, which we compared with experimental data in ref. [14].
Similar efforts were carried out by Csanád, Májer and Vargyas: they evaluated the direct photon spectrum [15,16] and hadronic observables [17] from a 1+3-dimensional analytic, but boost-invariant, solution to relativistic hydrodynamics for perfect fluids. Recently, a 1+1-dimensional perturbative solution to viscous hydrodynamics was used to derive an analytical formula for the rapidity and pseudorapidity distributions, analytically providing the correction of the bulk viscosity effect [18]. The shear and bulk viscous hydrodynamic deformation of the rapidity distributions were also investigated by Hirano and Monnai, using the Color Glass Condensate approach for relativistic heavy-ion collisions [19,20,21].
The justification for the use of a perfect fluid solution has been established in some of our previous works: we have found the property of asymptotic perfection for both nonrelativistic [22] and relativistic [23] Hubble-type viscous solutions. This asymptotic perfection means that viscous solutions become asymptotically equal with corresponding perfect fluid solutions, which implies that the results of the final state measurements of hadronic observables can be described not only by viscous solutions but also by the corresponding perfect fluid solutions.
In this manuscript, using a new formula computed from the solution reported in ref. [9], we find a data-collapsing or -scaling behaviour in experimental rapidity distribution data, which is one of the manifestations of a beautiful hydrodynamic scaling behaviour. This phenomenon was observed by Carruthers and Duong-van at intermediate energies for p + p π + + X and p + p π + X reactions in refs. [24,25]. They introduced a Gaussian scaling function for the rapidity densities from Landau hydrodynamics. In this paper, we derive the same scaling function from a different hydrodynamic solution, and we show that the data-collapsing behaviour in rapidity densities holds for high-energy collisions. Our theoretical insight is successfully cross-checked against experimental data, as shown explicitly in Section 4.

2. Recapitulation of the CKCJ Solution

In this section, we briefly recapitulate the 1+1-dimensional, analytical, accelerating family of CKCJ solutions to relativistic perfect fluid hydrodynamics from ref. [9]. The solutions of this family solve the equations of hydrodynamics at zero chemical potential ( μ = 0 ) under the following equation of state:
ε = κ 0 p ,
where ε stands for the energy density, while the pressure is denoted by p and κ 0 stands for their average (temperature-independent) constant of proportionality. This quantity κ 0 is directly related to the speed of sound in the medium ( κ 0 = c s 2 ), as recently clarified in ref. [22]. For the four-velocity ( u μ ), the CKCJ generalisation of the Hwa–Bjorken flow that allows the description of accelerating fluids reads as follows [9]:
u μ = cosh Ω , sinh Ω ,
where Ω is the fluid rapidity, which is a function of the coordinate rapidity η z , but it is independent from the longitudinal proper time τ (i.e., Ω Ω ( η z ) ) [9]. The new family of CKCJ solutions [9] can be summarised as follows:
η z ( H ) = Ω ( H ) H ,
Ω ( H ) = λ λ 1 κ 0 λ arctan κ 0 λ λ 1 tanh H ,
σ ( τ , H ) = σ 0 τ 0 τ λ V ( s ) 1 + κ 0 1 λ 1 sinh 2 ( H ) λ 2 ,
T ( τ , H ) = T 0 τ 0 τ λ κ 0 1 V ( s ) 1 + κ 0 1 λ 1 sinh 2 ( H ) λ 2 κ 0 ,
s ( τ , H ) = τ 0 τ λ 1 sinh ( H ) 1 + κ 0 1 λ 1 sinh 2 ( H ) λ / 2 ,
where T is the temperature and σ stands for the entropy density. The trajectory equation of the fluid element is defined by the scaling variable s, satisfying the scale equation u μ μ s = 0 . The function V ( s ) is a non-negative function of the scaling variable s. The parameter λ is an integration constant, but it has an important physical meaning, since it determines the acceleration of the velocity field: for λ > 1 , we talk about accelerating expansion, while, for 0 < λ < 1 , we consider decelerating expansion [9].
The above equations can be considered as a parametric family of solutions in the sense that the coordinate-rapidity dependence of the thermodynamic quantities (T, σ , p, ε ) and the elements of the velocity field ( Ω , v z ) are given through the parameter H = Ω ( η z ) η z [9].
It is important to note that the above equations define a finite family of solutions in terms of the coordinate rapidity, and the range of validity of the solutions is affected by the values of λ and κ 0 . However, in the accelerationless limiting case λ 1 , the width of the range of validity goes to infinity [9]. Thus, for low accelerations, this family of solutions can be applied in the physical domain of high-energy collisions. If the scaling function V ( s ) is set to V ( s ) = 1 , then the Hwa–Bjorken solution can be recovered in the λ 1 limit:
Ω = η z ,
σ ( τ ) = σ 0 τ 0 τ λ ,
T ( τ ) = T 0 τ 0 τ λ κ 0 .

3. Evaluation of the Rapidity Density

The Cooper–Frye formula can be used to calculate the rapidity distribution:
d N d y = 1 2 π d Σ μ p μ exp p μ u μ T F τ f , η z ,
where d Σ μ is the normal vector of the freeze-out hypersurface, while p μ is the four-momentum of the detected particles. The integral is performed on the freeze-out hypersurface, which is the set of points in spacetime τ , η z where the hadronic medium is freezing out. This integral is evaluated at vanishing chemical potential; thus, the fugacity yields a trivial factor of unity in the Boltzmann integral, where T F is the temperature of the freeze-out hypersurface, which is not a constant, since the points on the freeze-out hypersurface have different temperatures. The proper time associated with the freeze-out is denoted by τ f . For the temperature, we use Equation (6), and we consider the simplest possible case, i.e., when the scaling function is chosen to be unity: V ( s ) = 1 .
The new solution was found in 1+1-dimensional spacetime, so it was necessary to embed the formula for the rapidity distribution in 1+3-dimensional spacetime. To accomplish this, we first assumed that the temperature of the medium is homogeneous in the transverse plane, and, thus, the rapidity density is independent of the transverse coordinates. As a second step, we introduced the transverse mass m T instead of the particle mass m. In such a case, the integration of the distribution on the surface perpendicular to the beam direction can be performed trivially, but the integration on momentum space must also be performed. Turning to the transverse mass as an integration variable, the rapidity distribution embedded in the three-dimensional space can be obtained by integrating the double-differential spectrum over the transverse mass. We used saddle point approximation in the calculations, which led to the following result:
d N d y d N d y y = 0 cosh α ( κ 0 ) 2 1 y α exp m T f cosh α ( κ 0 ) y α 1 ,
where α ( κ 0 ) = ( 2 λ κ 0 ) / ( λ κ ) and α ( 1 ) = α . The kinetic freeze-out temperature is denoted by T f and y stands for the rapidity.
In Equation (12), the contribution of radial flow is not present, but, from the p T -spectra, it is known that the radial flow effect is not negligible. While the embedding of the rapidity distribution in three spatial dimensions is an efficient tool, it makes the assumption that the temperature is homogeneous in the transverse plane. For this reason, the effect of radial flow does not appear, so we artificially include it in the equations. We introduce the effective temperature T eff , which is equal to the sum of the kinetic freeze-out temperature T f and the contribution of the radial flow. In the formula for the rapidity distribution, we simply use the notation T f T eff to include the radial flow contribution into the distribution [9,10]:
d N d y d N d y y = 0 cosh α ( κ 0 ) 2 1 y α exp m T eff cosh α ( κ 0 ) y α 1 .
The midrapidity density can be expressed by the following formulae:
d N d y y = 0 = 1 λ ( 2 λ 1 ) V x V p ( 2 π ) 3 exp m T eff ,
V x = R 2 π τ f ,
V p = ( 2 π m T eff ) 3 / 2 .
Here, R 2 π is the finite size of the transverse plane, V x is the volume of the coordinate space and V p is the volume of the momentum space. While the midrapidity density depends on several parameters of the model, for simplicity, we now consider it only as a normalisation constant, which will greatly facilitate the fitting of experimental data.

3.1. The Approximate Formula for the Rapidity Distribution

We now summarise, in a few lines, that further approximations can be made on Equation (13) of the rapidity distribution in the region where | y | α = 2 + 1 / ( λ 1 ) . Obviously, this region widens rapidly if one approaches the accelerationless limiting case ( λ 1 ). In this approximation, the rapidity distribution becomes a Gaussian distribution [11]:
d N d y N 2 π Δ y 2 1 / 2 exp y 2 2 Δ y 2 .
In this approximation, the two new parameters are the average multiplicity N and Δ y , which characterises the width of the distribution. Both quantities can be expressed in terms of the parameters of the rapidity density embedded in the three-dimensional space [11]:
1 Δ y 2 = λ 1 2 1 + 1 + 1 κ 0 1 2 + m T eff ,
N = 2 π Δ y 2 1 / 2 d N d y y = 0 .
Such a Gaussian rapidity distribution is by no means a unique result that is characteristic only to the CKCJ solution. Other calculations and model results that are based on Landau hydrodynamics usually predict a similar result; see, e.g., refs. [3,24,25,26,27,28,29,30,31,32,33]. It is shown in ref. [34] that Landau hydrodynamics can describe both the measured rapidity distributions and the results of the EPOS 1.99 event generator. Recently, the Landau hydrodynamical model was used in combination with an ANN (Artifical Neural Network) simulation model to estimate the d N / d y distributions in central A u + A u collisions for two future experimental facilites: FAIR and NICA [35].
Despite the similarity of the Gaussian rapidity distributions in the Landau and in the CKCJ solutions, there is an important technical difference: the Landau solution is an implicit solution, where the time and coordinate variables are given in terms of the temperature and the fluid rapidity, while the CKCJ solution is a parametric, but explicit, solution, where the observables are expressed in terms of the acceleration parameter λ , the average speed of sound κ 0 , the effective temperature (slope parameter) at midrapidity T eff and the midrapidity density d N / d y y = 0 —see Equations (13) and (17). This becomes advantageous when analysing the scaling properties of the rapidity distribution, as it gives a beautiful example of hydrodynmamical scaling.
Using Equation (17), the experimental data can be described by fitting two parameters ( N and Δ y ), while the original formula (Formula (13)) has four parameters to be fitted ( κ 0 , λ , T eff , d N / d y y = 0 ). This result is a wonderful manifestation of hydrodynamic scaling behaviour: two different fits of the data, resulting in different parameter sets, can lead to the same curve, provided that the values of the two relevant combinations of the four original parameters, (18) and (19), remain the same. From Equation (17) of the rapidity distribution, it is also clear that the rapidity density can be normalised by both the midrapidity value of the distribution and the mean multiplicity.

3.2. The Approximate Formula for the Pseudorapidity Distribution

The pseudorapidity distribution can be calculated from the double-differential invariant momentum spectrum and can be written as the product of two factors (as detailed in ref. [11]). One factor is the rapidity density, as written in Equation (13). The other factor is the Jacobian determinant J , calculated for a mean rapidity-dependent transverse momentum: p T ( y ) [11]. The results in refs. [36,37] have shown that the effect of transverse flow goes through the rapidity dependence of the average transverse momentum dependence of J . The rapidity-dependent average transverse mass was evaluated by saddle-point approximation, which limits our investigation to the nearly exponential ( p T < 1 2 GeV) kinematic limit. In the leading order, the rapidity dependence of the average transverse mass has Lorentzian shape [9,11]:
m T ( y ) = p T ( y ) 2 + m 2 = m + T eff 1 + b y 2 ,
where b = α ( κ 0 ) / ( 2 α 2 ) controls the width of the Lorentz distribution. This result was mentioned first just after Equation (36) in ref. [9] and it has been discussed, in more detail, that the CKCJ solution has a realistic feature: a rapidity-dependent decrease in the average transverse momentum. According to Equation (20), the physical meaning of the effective temperature becomes clear: it is related to the thermally averaged transverse mass at midrapidity. A Lorentzian shape for the rapidity-dependence of the average transverse momentum has been observed by the EHS/NA22 experiment in slightly asymmetric hadron–proton collisions [36].
The pseudorapidity distribution and the partial results of its derivation were described in detail in ref. [11], so we leave it aside for now. In this manuscript, we focus on an approximate formula that holds in the limiting case λ 1 . Note that, in this limit, b 0 , so the average rapidity-dependent transverse momentum can be approximated by its central value: p T ( y ) p T ( 0 ) = p T = T eff 2 + 2 m T eff . In this case, the Jacobi determinant and the pseudorapidity distribution can be approximated using the following equations [11]:
J cosh ( η p ) D 2 + cosh 2 ( η p ) ,
d N d η p N 2 π Δ y 2 cosh ( η p ) D 2 + cosh 2 ( η p ) exp y 2 η p 2 Δ y 2 ,
where η p is the pseudorapidity and D = m / p T stands for the depthness parameter. In Equation (22), the y ( η p ) rapidity can be approximated as follows [11]:
y ( η p ) tanh 1 cosh ( η p ) D 2 + cosh 2 ( η p ) tanh η p .

4. The Scaling of dN/dy Data

The Gaussian density is expressed by Equation (17) and is determined by the normalisation factor and the width of the distribution. These parameters thus carry the characteristics of the different reactions, such as the collision energy, the size of the colliding system and the centrality of the collision. The normalisation factor and the width of the distribution can be scaled out from Equation (17), and, with that, we obtain a curve that is independent from the characteristics of the reactions. We introduce the variable x, which is the ratio of the rapidity y to the width Δ y :
x = y Δ y .
Then, we divide the rapidity distribution by the normalisation factor, introducing a scale-independent distribution:
d N d y y = 0 1 d N d y = exp x 2 2 ,
where the expression of the normalisation factor is given by Equation (19), which is valid in the | y | α range. Our idea is that, if we scale back the rapidity distribution data from different reactions with the normalisation factor and Δ y as above, then, regardless of the reaction, all the data series will follow the curve of Equation (25). However, in most cases, we do not have d N / d y measurements available, so, in order to check our conjecture, we have transformed the available pseudorapidity density data. For this reason, first we used Equation (22) to determine the parameters D, N and Δ y from the pseudorapidity distribution data. Thus, for each dataset, we obtained different values of D, N and Δ y . Using the D parameter, we extracted the Jacobian determinant to transform the pseudorapidity distribution data points into rapidity distribution data points. Then the mean multiplicity N and the width parameter Δ y were used to rescale the rapidity distributions. Using this method, we checked the manifestation of the scaling behaviour on 11 datasets: PHOBOS 0–30% A u + A u at s N N = 20 , 62.4, 130 [38], CMS p + p at s = 7 , 8 and 13 TeV and 0–80% X e + X e at s N N = 5.44 TeV [39,40,41,42] and ALICE P b + P b at s N N = 5.02 TeV for central, mid-central and peripheral collisions [43]. We note that the PHOBOS data required a special approach, as PHOBOS did not publish the statistical and systematic errors of the pseudorapidity distributions separately. For this reason, we used the best estimate for the statistical errors in the PHOBOS data, described in detail in ref. [11].
A comparison of these data with Equation (25) is shown in Figure 1. The vertical errors of the data points correspond to the quadrature of the statistical errors associated with the pseudorapidity distributions and the uncertainty in the Jacobian determinants. The horizontal error of the data points is due to the uncertainty in Δ y . Figure 1 convincingly illustrates that our conjecture was right, i.e., all the data points measured in the different reactions are clustered (within error) on the blue curve, which is the manifestation of hydrodynamic scaling behaviour in the experimental data.
For the pseudorapidity distributions, we could not find a scale function similar to Equation (25) to describe the collapse of the experimental datasets of d N / d η p , as we saw for the rapidity distributions in Figure 1. This is due to the presence of the depth parameter D, which controls the dip of the pseudorapidity distribution around midrapidity. This can be easily seen by evaluating Equation (22) for the case of zero pseudorapidity:
d N d η p η p = 0 = d N d y y = 0 1 1 + D 2 ,
where we used Equation (19) for the midrapidity density. Thus, if this dip is close to zero ( D 0 ), we recover the formula of the rapidity distribution and the scaling behaviour is restored. This phenomenon is illustrated in Figure 2.

5. Conclusions

In this manuscript, we introduced a new approximate formula to describe the rapidity distributions, based on the solution published in ref. [9]. By rescaling the new formula, we introduced a new scale function to show the data collapsing in the rapidity distribution measurements. This data-collapsing (or scaling) behaviour was tested on 11 datasets from three different experiments (CMS, PHOBOS, ALICE). It is important to emphasise that these datasets are not only from nucleus–nucleus collisions, but also from proton–proton collisions, which are also described by the curve of the scale function. From this result, we can conclude that p + p collisions can also be treated as collective systems, as far as their rapidity or pseudorapidity distributions are concerned. Let us mention that the p + p data can be well-described with low speed of sound values ( c s 2 = 0.1 ), which is typical for fluids. This suggests that a strongly interacting quark–gluon plasma may be formed in p + p collisions at RHIC and at LHC energies. Although p + p collisions at RHIC energies have not been investigated in this manuscript, it is expected that the rapidity distributions of these systems do not violate the scaling behaviour presented in Section 4. It has been shown, in ref. [37], that p T spectra and HBT radii measured in s = 200 GeV p + p collisions at RHIC can be described by the Buda–Lund hydrodynamical model.
Let us emphasise that the hydrodynamical description of small systems (hadron–hadron collisions) is not a new idea. The EHS/NA22 collaboration previously estimated some hydrodynamic model parameters (radial flow, kinetic freeze-out temperature, geometric radius of the particle source, mean lifetime) via invariant momentum distributions of π particles and Bose–Einstein correlation functions in ( π + / K + ) p reactions at s 22 GeV, where the mean multiplicity of charged particles was less then 10 [36].

Author Contributions

Writing—original draft preparation, G.K.; writing—review and editing, T.C. and G.K.; conceptualisation, T.C.; supervision, T.C.; visualisation, G.K. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by NKFIH (National Office of Research, Innovation and Development, Hungary) K-133046 and K-147557, and the KKP-2023 Research Excellence Programme of MATE, Hungary. The APC was funded by NKFIH K-133046.

Data Availability Statement

No new data were created or analyzed in this study. Data sharing is not applicable to this article.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Fermi, E. High-energy nuclear events. Prog. Theor. Phys. 1950, 5, 570–583. [Google Scholar] [CrossRef]
  2. Fermi, E. Angular Distribution of the Pions Produced in High Energy Nuclear Collisions. Phys. Rev. 1951, 81, 683–687. [Google Scholar] [CrossRef]
  3. Belenkij, S.Z.; Landau, L.D. Hydrodynamic theory of multiple production of particles. Nuovo Cim. Suppl. 1956, 3S10, 15. [Google Scholar] [CrossRef]
  4. Landau, L.D. On Multiple Production of Particles during Collisions of Fast Particles. Izv. Akad. Nauk Ser. Fiz. 1953, 15, 51. [Google Scholar] [CrossRef]
  5. Hwa, R.C. Statistical Description of Hadron Constituents as a Basis for the Fluid Model of High-Energy Collisions. Phys. Rev. D 1974, 10, 2260. [Google Scholar] [CrossRef]
  6. Bjorken, J.D. Highly Relativistic Nucleus-Nucleus Collisions: The Central Rapidity Region. Phys. Rev. D 1983, 27, 140–151. [Google Scholar] [CrossRef]
  7. Csörgő, T.; Nagy, M.I.; Csanád, M. A New family of simple solutions of perfect fluid hydrodynamics. Phys. Lett. B 2008, 663, 306–311. [Google Scholar] [CrossRef]
  8. Nagy, M.I.; Csörgő, T.; Csanád, M. Detailed description of accelerating, simple solutions of relativistic perfect fluid hydrodynamics. Phys. Rev. C 2008, 77, 024908. [Google Scholar] [CrossRef]
  9. Csörgő, T.; Kasza, G.; Csanád, M.; Jiang, Z.F. New exact solutions of relativistic hydrodynamics for longitudinally expanding fireballs. Universe 2018, 4, 69. [Google Scholar] [CrossRef]
  10. Csörgő, T.; Kasza, G.; Csanád, M.; Jiang, Z.F. A new and finite family of solutions of hydrodynamics. Part I: Fits to pseudorapidity distributions. Acta Phys. Polon. B 2019, 50, 27–35. [Google Scholar] [CrossRef]
  11. Kasza, G.; Csörgő, T. Lifetime estimations from RHIC Au+Au data. Int. J. Mod. Phys. A 2019, 34, 1950147. [Google Scholar] [CrossRef]
  12. Kasza, G.; Csörgő, T. A new and finite family of solutions of hydrodynamics: Part II: Advanced estimate of initial energy densities. In Proceedings of the 13th Workshop on Particle Correlations and Femtoscopy, Cracow, Poland, 22–26 May 2018. [Google Scholar] [CrossRef]
  13. Csörgő, T.; Kasza, G. A new and finite family of solutions of hydrodynamics: Part III: Advanced estimate of the life-time parameter. In Proceedings of the 13th Workshop on Particle Correlations and Femtoscopy, Cracow, Poland, 22–26 May 2018. [Google Scholar] [CrossRef]
  14. Kasza, G. Describing the thermal radiation in Au + Au collisions at s N N = 200 GeV by an analytic solution of relativistic hydrodynamics. arXiv 2023, arXiv:2311.03568. [Google Scholar]
  15. Csanád, M.; Májer, I. Initial temperature and EoS of quark matter from direct photons. Phys. Part. Nucl. Lett. 2011, 8, 1013–1015. [Google Scholar] [CrossRef]
  16. Csanád, M.; Májer, I. Equation of state and initial temperature of quark gluon plasma at RHIC. Cent. Eur. J. Phys. 2012, 10, 850–857. [Google Scholar] [CrossRef]
  17. Csanád, M.; Vargyas, M. Observables from a solution of 1 + 3 dimensional relativistic hydrodynamics. Eur. Phys. J. A 2010, 44, 473–478. [Google Scholar] [CrossRef]
  18. Jiang, Z.F.; Yang, C.B.; Ding, C.; Wu, X.Y. Pseudo-rapidity distribution from a perturbative solution of viscous hydrodynamics for heavy ion collisions at RHIC and LHC. Chin. Phys. C 2018, 42, 123103. [Google Scholar] [CrossRef]
  19. Monnai, A.; Hirano, T. Longitudinal Viscous Hydrodynamic Evolution for the Shattered Colour Glass Condensate. Phys. Lett. B 2011, 703, 583–587. [Google Scholar] [CrossRef]
  20. Monnai, A.; Hirano, T. Viscous Hydrodynamic Evolution with Non-Boost Invariant Flow for the Color Glass Condensate. J. Phys. G 2011, 38, 124168. [Google Scholar] [CrossRef]
  21. Monnai, A.; Hirano, T. Viscous hydrodynamic deformation in rapidity distributions of the color glass condensate. In Proceedings of the Seventh Workshop on Particle Correlations and Femtoscopy (WPCF2011), Tokyo, Japan, 20–24 September 2011. [Google Scholar] [CrossRef]
  22. Kasza, G.; Csernai, L.P.; Csörgő, T. New, Spherical Solutions of Non-Relativistic, Dissipative Hydrodynamics. Entropy 2022, 24, 514. [Google Scholar] [CrossRef]
  23. Csörgő, T.; Kasza, G. New, multipole solutions of relativistic, viscous hydrodynamics. arXiv 2020, arXiv:2003.08859. [Google Scholar] [CrossRef]
  24. Carruthers, P.; Duong-Van, M. New scaling law based on the hydrodynamical model of particle production. Phys. Lett. B 1972, 41, 597–601. [Google Scholar] [CrossRef]
  25. Carruthers, P.; Doung-van, M. Rapidity and angular distributions of charged secondaries according to the hydrodynamical model of particle production. Phys. Rev. D 1973, 8, 859–874. [Google Scholar] [CrossRef]
  26. Shuryak, E.V. Multiparticle production in high energy particle collisions. Yad. Fiz. 1972, 16, 395–405. [Google Scholar]
  27. Zhirov, O.V.; Shuryak, E.V. Multiple Production of Hadrons and Predictions of the Landau Theory. Yad. Fiz. 1975, 21, 861–867. [Google Scholar]
  28. Steinberg, P. Landau hydrodynamics and RHIC phenomena. Acta Phys. Hung. A 2005, 24, 51–57. [Google Scholar] [CrossRef]
  29. Wong, C.Y. Landau Hydrodynamics Revisited. Phys. Rev. C 2008, 78, 054902. [Google Scholar] [CrossRef]
  30. Wong, C.Y.; Sen, A.; Gerhard, J.; Torrieri, G.; Read, K. Analytical Solutions of Landau (1 + 1)-Dimensional Hydrodynamics. Phys. Rev. C 2014, 90, 064907. [Google Scholar] [CrossRef]
  31. Sen, A.; Gerhard, J.; Torrieri, G.; Read, K.; Wong, C.Y. Hydrodynamics from Landau initial conditions. J. Phys. Conf. Ser. 2015, 630, 012042. [Google Scholar] [CrossRef]
  32. Gazdzicki, M.; Gorenstein, M.; Seyboth, P. Onset of deconfinement in nucleus-nucleus collisions: Review for pedestrians and experts. Acta Phys. Polon. B 2011, 42, 307–351. [Google Scholar] [CrossRef]
  33. Jiang, Z.J.; Li, Q.G.; Zhang, H.L. The revised Landau hydrodynamic model and the pseudorapidity distributions of produced charged particles in high energy heavy ion collisions. J. Phys. G 2013, 40, 025101. [Google Scholar] [CrossRef]
  34. Rahman, R.M.A.; Tawfik, A.N.; El-Bakry, M.Y.; Habashy, D.M.; Hanafy, M. Rapidity distribution within Landau hydrodynamical model and EPOS event-generator at AGS, SPS, and RHIC energies. Phys. Scr. 2022, 97, 065305. [Google Scholar] [CrossRef]
  35. Habashy, D.M.; El-Bakry, M.Y.; Tawfik, A.N.; Abdel Rahman, R.M.; Hanafy, M. Particles multiplicity based on rapidity in Landau and artificial neural network (ANN) models. Int. J. Mod. Phys. A 2022, 37, 2250002. [Google Scholar] [CrossRef]
  36. Agababyan, N.M. et al. [EHS/NA22 Collaboration] Estimation of hydrodynamical model parameters from the invariant spectrum and the Bose-Einstein correlations of π mesons produced in (π+/K+)p interactions at 250 GeV/c. Phys. Lett. B 1998, 422, 359–368. [Google Scholar] [CrossRef]
  37. Csörgő, T.; Csanád, M.; Lörstad, B.; Ster, A. Analysis of identified particle yields and Bose-Einstein (HBT) correlations in p+p collisions at RHIC. Acta Phys. Hung. A 2005, 24, 139–144. [Google Scholar] [CrossRef]
  38. Alver, B. et al. [PHOBOS Collaboration] Charged-particle multiplicity and pseudorapidity distributions measured with the PHOBOS detector in Au+Au, Cu+Cu, d+Au, and p+p collisions at ultrarelativistic energies. Phys. Rev. C 2011, 83, 024913. [Google Scholar] [CrossRef]
  39. Khachatryan, V. et al. [CMS Collaboration] Transverse-momentum and pseudorapidity distributions of charged hadrons in pp collisions at s = 7 TeV. Phys. Rev. Lett. 2010, 105, 022002. [Google Scholar] [CrossRef] [PubMed]
  40. Chatrchyan, S. et al. [CMS and TOTEM Collaborations] Measurement of pseudorapidity distributions of charged particles in proton-proton collisions at s = 8 TeV by the CMS and TOTEM experiments. Eur. Phys. J. C 2014, 74, 3053. [Google Scholar] [CrossRef]
  41. Sirunyan, A.M. et al. [CMS Collaboration] Measurement of charged pion, kaon, and proton production in proton-proton collisions at s = 13 TeV. Phys. Rev. D 2017, 96, 112003. [Google Scholar] [CrossRef]
  42. Sirunyan, A.M. et al. [CMS Collaboration] Pseudorapidity distributions of charged hadrons in xenon-xenon collisions at s N N = 5.44 TeV. Phys. Lett. B 2019, 799, 135049. [Google Scholar] [CrossRef]
  43. Adam, J. et al. [ALICE Collaboration] Centrality dependence of the pseudorapidity density distribution for charged particles in Pb-Pb collisions at s N N = 5.02 TeV. Phys. Lett. B 2017, 772, 567–577. [Google Scholar] [CrossRef]
Figure 1. Manifestation of hydrodynamic scaling behaviour on experimental data: by transforming the pseudorapidity distributions measured in different reactions, we calculated the data points associated with the rapidity densities, which we rescaled by the normalisation factor. The width of the distributions were also scaled by introducing the variable x. All of the resulting data points lie within error on the scale-independent blue curve obtained from our theoretical calculations. The datasets were taken from CMS [39,40,41,42], ALICE [43] and PHOBOS [38].
Figure 1. Manifestation of hydrodynamic scaling behaviour on experimental data: by transforming the pseudorapidity distributions measured in different reactions, we calculated the data points associated with the rapidity densities, which we rescaled by the normalisation factor. The width of the distributions were also scaled by introducing the variable x. All of the resulting data points lie within error on the scale-independent blue curve obtained from our theoretical calculations. The datasets were taken from CMS [39,40,41,42], ALICE [43] and PHOBOS [38].
Universe 10 00045 g001
Figure 2. The pseudorapidity distribution (divided by d N / d y y = 0 ) for different D values. The width parameter Δ y is fixed to 1. It can be clearly seen that, as D approaches 0 more closely, midrapidity dip in the pseudorapidity distribution becomes smaller, and, for D = 0 , we recover Equation (17) for the shape of the rapidity distribution.
Figure 2. The pseudorapidity distribution (divided by d N / d y y = 0 ) for different D values. The width parameter Δ y is fixed to 1. It can be clearly seen that, as D approaches 0 more closely, midrapidity dip in the pseudorapidity distribution becomes smaller, and, for D = 0 , we recover Equation (17) for the shape of the rapidity distribution.
Universe 10 00045 g002
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Kasza, G.; Csörgő, T. Scaling Behaviour of dN/dy in High-Energy Collisions. Universe 2024, 10, 45. https://doi.org/10.3390/universe10010045

AMA Style

Kasza G, Csörgő T. Scaling Behaviour of dN/dy in High-Energy Collisions. Universe. 2024; 10(1):45. https://doi.org/10.3390/universe10010045

Chicago/Turabian Style

Kasza, Gábor, and Tamás Csörgő. 2024. "Scaling Behaviour of dN/dy in High-Energy Collisions" Universe 10, no. 1: 45. https://doi.org/10.3390/universe10010045

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop