# Transverse Momentum Distributions of Hadrons Produced in Pb-Pb Collisions at LHC EnergysNN=2.76 TeV

Authors: Saeed Uddin; Inam-ul Bashir; Riyaz Ahmed Bhat
Journal: Advances in High Energy Physics
 (2015)

---

## Abstract

The transverse momentum spectra of several types of hadrons,p, p̅, K+, K-, Ks0, Λ, Ω, Ω̅, Ξ-, and Ξ̅ produced in most central Pb-Pb collisions at LHC energy sNN=2.76 TeV have been studied at midrapidity (|y|<0.5) using an earlier proposed unified statistical thermal freeze-out model. The calculated results are found to be in good agreement with the experimental data measured by the ALICE experiment at LHC. The model calculation fits provide the thermal freeze-out conditions in terms of the temperature and collective flow effect parameters for different particle species. Interestingly the model parameter fits to the experimental data reveal stronger collective flow in the system and lesser freeze-out temperatures of the different particle species as compared to Au-Au collisions at RHIC. The strong increase of the collective flow appears to be a consequence of the increasing particle density at LHC. The model used incorporates a longitudinal as well as transverse hydrodynamic flow. The chemical potential has been assumed to be nearly equal to zero for the bulk of the matter owing to high degree of nuclear transparency effect at such collision energies. The contributions from heavier decay resonances are also taken into account.


---

## Body

## 1. Introduction

The study of identified particle spectra in heavy-ion collisions at ultrarelativistic energies is an important tool to investigate the properties of the strongly interacting system created in such collisions. The study also helps us to learn about the final state distribution of baryon numbers among various particle species at the thermochemical freeze-out after the collision which is initially carried by the nucleons only [1].Within the framework of the statistical model, it is assumed that a hot and dense fireball is formed over an extended region for a brief period of time (~ a few fm/c) after the initial collision and it undergoes collective expansion leading to a decrease in its temperature and finally to the hadronization. After the hot fireball formed in such collisions, which initially has a very high density of partons (i.e., quarks and gluons), hadronizes, the hadrons keep rescattering with each other and continue to build up collective expansion. Consequently, the matter becomes dilute and the average distance between hadrons exceeds the range of the strong interactions. At this point of time, all scattering processes stop and the hadrons decouple; that is, a freeze out occurs [2].The hadronic abundance freeze-out (i.e., the chemical freeze-out) occurs earlier when the rates for inelastic processes, in which secondary hadrons are produced or the hadrons change their identity (via strangeness exchange process, etc.), become too small to keep up with the expansion. Since the corresponding inelastic cross sections are only a small fraction of the total cross section at lower (thermal) energies, the inelastic processes stop well before theelastic ones, leading to an earlier chemical freeze-out for the hadron abundances. Finally at a later stage the hadrons completely decouple from each other such that even the elastic processes also come to a stop. Consequently, the momentum spectra get frozen in time and a thermal (or hydrodynamical) freeze-out occurs. Thus chemicalfreeze-outprecedes thermal or kinetic freeze-out [3].Transverse momentum (pT) distributions of identified hadrons are the most common tools used to study the dynamics of high energy collisions.The identified particle spectra provide information, about both the temperature of the system and the collective flow at the time of thermal freeze-out. Collective flow depends on the internal pressure gradients created in the collision and is addressed by hydrodynamic models [4–6]. These effects are species-dependent. The produced hadrons are believed to carry information about the collision dynamics and the subsequent space-time evolution of the system.Hence an accurate measurement of the transverse momentum distributions of identified hadrons along with the rapidity spectra is essential for the understanding of the dynamics and the properties of the created matter up to the final thermal or hydrodynamical freeze-out in case of collective flow [7].It has been shown earlier [2] that this model successfully simultaneously explains the rapidity and transverse momentum distributions of hadrons and their ratios in Au-Au collisions at highest RHIC energy of sNN=200 GeV. In this paper, we briefly describe the model and use it to reproduce the transverse momentum distributions of hadrons produced in Pb-Pb collisions at sNN=2.76 TeV.


## 2. Model

In order to obtain the particle spectra in the overall rest frame of the hadronic fireball in our model we first define the invariant cross-section for given hadronic specie in the local rest frame of a hadronic fluid element. Since the invariant cross section will have the same value in all Lorentz frames [8], we can thus write(1)Ed3Nd3p=E′d3Nd3p′,where E(E′) is the energy of the particle and p(p′) is the momentum. The primed quantities on the RHS refer to the invariant spectra of given hadronic specie in the rest frame of the local hadronic fluid element, while the unprimed quantities on the LHS refer to the invariant spectra of the same hadronic specie in the overall rest frame of the hadronic fireball. The occupation number distribution of the hadrons in the momentum space follows the distribution function:(2)E′d3Nd3p′~E′e((E′-μ)/T)±1,where (+) sign and (−) sign are for fermions and bosons, respectively, and μ is the chemical potential of the given hadronic specie. For the temperatures under consideration and the large masses of hadrons it is safe to work with Boltzman distribution.In recent works [7, 9] it has been clearly shown that there is a strong evidence of increasing baryon chemical potential, μB, along the collision axis in the RHIC experiments. In view of this fact we write the expression for chemical potential as μB=a+by02 [7, 9, 10], where y0 is the rapidity of the expanding hadronic fluid element. Here a and b are the two model parameters which can be fixed by fitting the experimental data. In the model the value of a essentially defines the baryon chemical potential in the central region of the bulk hadronic matter formed, while b determines the rate of increase of baryon chemical potential along the (longitudinal) rapidity axis with y0. In case of very high degree of nuclear transparency the values of a and b will tend to vanish for the bulk of the matter. Further it is assumed that [7] the rapidity of the expanding hadronic fluid element y0αz or y0=ξz, where z is the longitudinal coordinate of the hadronic fluid element and ξ is a proportionality constant. The above conditions also ensure that, under the transformation z→-z, we will have y0→-y0, thereby preserving the symmetry of the hadronic fluid flow about z=0 along the rapidity axis in the centre of mass frame of the colliding nuclei. This leads to an expression for the longitudinal velocity component of the hadronic fluid element:(3)βzz=1-2exp⁡⁡2ξz+1=tanh⁡y0.The transverse velocity component of the hadronic fireball, βT is assumed to vary with the transverse coordinate r in accordance with the blast wave model as βT(r)=βTs(r/R)n [11], where n is an index which fixes the profile of βT(r) in the transverse direction and βTs is the hadronic fluid surface transverse expansion velocity and is fixed in the model by using the parameterization βTs=βT01-βz2 [7]. This relation is also required to ensure that the net velocity β of any fluid element must satisfy β=βT2+βz2<1. We also parameterize R, that is, the transverse radius of fireball as R=r0exp⁡(-z2/σ2) where σ fixes the width of the matter distribution in the transverse direction [7, 9] and z, as described above, is the longitudinal coordinate of hadronic fluid element. This is required as the colliding nuclei when passing through each other may still feel some drag thus resulting only in a partial transparency. Consequently, the collision axis will be populated by an extended hadronic matter rapidly moving away from each other with its transverse size decreasing rapidly following a Gaussian distribution along the z-axis.In our analysis, the contributions of various heavier hadronic resonances [10, 12] which decay after the thermal freeze-out of the hadronic matter has occurred are also taken into account. The invariant spectrum of a given decay product of a given parent hadron in the local rest frame of a hadronic fluid element is written as [7, 10, 12](4)E′d3Ndecayd3p′=12p′mhp*∫E-E+dEhEhd3Nhd3ph,where the subscript h stands for the decaying (parent) hadron. The two body decay kinematics gives the product hadron’s momentum and energy in the “rest frame of the decaying hadron” as p*=E*2-m21/2 and E*=(mh2-mj2+m2)/2mh, where mj indicates the mass of the other decay hadron produced along with the first one. The limits of integration are E±={mh/m2}{E′E*±p′p*}. The E′(Eh) and p′(ph) are, respectively, the product (decaying parent) hadron’s energy and momentum in the local rest frame of the hadronic fluid element. A Boltzmann type distribution for the massive decaying hadron in the local rest frame of the hadronic fluid element leads to the following final expression for the invariant cross section of the product hadron:(5)E′d3Nd3p′=12p′mhp*λhghe-αθE′E*×αθE′E*sinh⁡⁡αθp′p*-p′p*cosh⁡⁡αθp′p*+T2sinh⁡⁡αθp′p*αθE′E*sinh⁡⁡αθp′p*-p′p*cosh⁡⁡αθp′p*,where α and θ are given by mh/m2 and 1/T, respectively.


## 3. Results and Discussions

We employ the minimumx2/DoF method to fit the experimental data. We find that the model calculations results (shown by solid curves in all the cases) fit the experimental data quite well (shown by filled circles in all the cases). The experimental data are taken from the ALICE Collaboration for Pb-Pb collisions at sNN=2.76 TeV [13–15]. We have shown the (statistical + systematic) errors in all the cases.Over a fairly largepT range the hydrodynamical calculations show an approximate exponential behavior, whereas the tails of measured spectra show a significant deviation in the slope beyond 5 GeV at LHC. At RHIC this transition from exponential behavior takes place at pT≳3 GeV. The fraction of hadrons with very large pT (≥3 GeV at RHIC and ≥5 GeV at LHC) is however small. We have considered the (maximum) pT range up to 5 GeV in the present analysis. It is because that the statistical hydrodynamic calculations cannot describe the hadron spectra at such large transverse momenta. The hadrons detected in this region are essentially formed by the partons which are result of the hard processes. These originate from the direct fragmentation of high-energy partons of the colliding beams and therefore are not able to thermalize through the process of multiple collisions [16]. We therefore turn to softer hadrons which are assumed to be reasonably thermalized and form the bulk of the secondary matter produced.The applicability region of hydrodynamics at LHC is therefore predicted to be forpT≤4-5 GeV depending on the particle’s mass. This range is wider than at RHIC [17]. The transverse momentum distributions are found to be sensitive to the values of the thermal/kinetic freeze-out temperature T and the transverse flow parameter βT0, whereas it is found to be insensitive to the change in the values of σ in our model. In our analysis we have therefore fixed the value of the parameter σ=5.0. This value essentially determines the size of the hadronic matter distributed along the z-axis and has a strong effect on the shape of the rapidity spectra of the particles. In our earlier analysis [2] of the RHIC data the value of σ turned out to be nearly 4.2. However, it is expected to be large at the LHC energy. The insensitivity of the transverse momentum distribution to the parameter σ has been tested and it is found that the minimum x2/DoF for protons varies only from 0.611 to 0.614 if σ is varied from 4.0 to 6.0. We have taken the values of a and b both to be zero for all the hadrons under the assumption of a baryon symmetric matter expected to be formed under the condition of a high degree of nuclear transparency in the nucleus-nucleus collisions at LHC energy that is, an ideal Bjorken picture. Unlike the previous works we have in our present analysis treated the index parameter n as a free parameter. The values of the parameters, T, βT0 and n, at freeze-out are determined by obtaining a best fit to a given hadron’s transverse momentum spectrum. The value of ξ=1 is fixed for all the hadrons studied in this paper. The theoretical fits for the transverse momentum spectra of all the hadrons have been normalized at the first data point (i.e., at the lowest pT) to facilitate a proper comparison with the experimental data set.In Figure1, we have shown the transverse momentum spectra of protons and antiprotons. The values of the thermal/kinetic freeze-out temperature T, the transverse flow parameter βT0, and the index parameter n for protons as well as antiprotons are found to be same, that is, 102 MeV, 0.88, and 1.40, respectively, with a minimum x2/DoF of 0.61 for protons and 0.55 for antiprotons. The same values of the freeze-out parameters for protons and antiprotons indicate a simultaneous freeze-out of these particles in the dense hadronic medium.Figure 1
Transverse momentum spectra of protons (a) and antiprotons (b) for the centrality class (0–5)%.
(a)
(b)The transverse momentum spectra forK+ and K- shown in Figure 2 gives the value of (T, βT0, and n) as (103 MeV, 0.89, and 1.80) for Kaons and (105 MeV, 0.88, and 1.80) for anti-Kaons. The minimum x2/DoF for both the two cases turns out to be 0.34. The almost similar freeze-out parameters obtained for protons, antiprotons, Kaons, and anti-Kaons indicate a near simultaneous freeze-out of these particles.Figure 2
Transverse momentum spectra ofK+ (a) and K- (b) for the centrality class (0–5)%.
(a)
(b)The transverse momentum spectrum of neutral Kaon, that is,Ks0, is shown in Figure 3. The values of T, βT0, and n obtained from the spectra of Ks0 are, respectively, 125 MeV, 0.84, and 1.61 with the minimum x2/DoF=1.70. The Ks0 shows a larger thermal freeze-out temperature than the charged Kaons indicating its earlier freeze-out than K±.Figure 3
Transverse momentum spectra ofKs0 for the centrality class (0–5)%.The transverse momentum spectra of hyperons (i.e.,Λ, Ξ-, Ξ--, Ω and Ω-) are shown in Figures 4, 5, and 6. The spectrum of Λ gives the values of T, βT0, and n as 127 MeV, 0.84, and 1.06, respectively, with a minimum x2/DoF=0.52. These values for Ξ- are found to be 133 MeV, 0.81, and 0.90 while for Ξ-- these are 149 MeV, 0.80, and 1.25, respectively. The parameters for Ω and Ω- are (155 MeV, 0.77, and 1.22) and (154 MeV, 0.77, and 1.23). The minimum x2/DoF for Ξ- and Ξ-- are 0.38 and 0.50, whereas for Ω and Ω- the minimum x2/DoF are 0.10 and 0.20, respectively. The relatively smaller values of minimum x2/DoF for Ωs are due to larger experimental error bars.Figure 4
Transverse momentum spectra of lambdaΛ for the centrality class (0–5)%.Figure 5
Transverse momentum spectra ofΞ- (a) and Ξ-- (b) for the centrality class (0–10)%.
(a)
(b)Figure 6
Transverse momentum spectra ofΩ (a) and Ω- (b) for the centrality class (0–10)%.
(a)
(b)The values of the thermal/kinetic freeze-out temperatureT, the transverse flow parameter βT0, and the index parameter n for all the hadrons studied in this paper are again presented in Table 1 to facilitate a proper comparison.Table 1
Freeze-out parameters of various hadrons obtained from their transverse   momentum spectra.
Particle
T (MeV)
β
T
0
n
x
2
/
DoF
p
102
0.88
1.40
0.61
p
-
102
0.88
1.40
0.55
K+
103
0.89
1.80
0.34
K
−
105
0.88
1.80
0.34
K
s
0
125
0.84
1.61
1.70
Λ
127
0.84
1.06
0.52
Ξ
-
133
0.81
0.90
0.38
Ξ
-
-
149
0.80
1.25
0.48
Ω
155
0.77
1.22
0.10
Ω
-
154
0.77
1.23
0.20It is evident from Table1 that the lighter particles, that is, (anti)protons and Kaons, exhibit a lower thermal freeze-out temperature and a higher surface transverse expansion velocity compared to the heavy (multi)strange hyperons. The reason for this can be attributed to an early freeze-out for the massive particles (hyperons) when the thermal temperature is high and the collective flow is in the early stage of development and consequently βT0 is small. The early freeze-out of these particles is due to their smaller cross-section with the hadronic matter. This can also be understood in terms of the mean free path, λ, of a particle in a thermal environment which is given by λ=1/νρ, where ν is the mean thermal cross section of the particle with the surrounding matter having density ρ. The transverse flow velocity profile index n, using the best fit method for the hadrons considered here, lies in the range 1.61–1.80 for the lightest particles considered here (i.e., Kaons), for (anti)protons it is 1.40 and for the heavier multistrange hyperon it lies in the range 0.90–1.25. Hence in general n appears to be smaller for the heavier particles. Another observation is that, for the particles with lower n, the pT spectra show a shoulder-like shape at low transverse momenta.A comparison with a similar fit to the RHIC data [2, 6] shows that for the most central collisions the flow velocity increases significantly at LHC reaching almost 0.9 and that the kinetic freeze-out temperature drops below the one at RHIC. For RHIC [2] the values of T and βT0 were found to be in the range 163–188 MeV and 0.58–0.67, respectively.We have also compared the results of some other model calculations like VISH2+1, HKM, and Krakow models [13, 18] with our model results in Figure 7 for the transverse momentum spectra of (Ξ-+Ξ--)/2 and (Ω+Ω-)/2. The Krakow model results are available up to nearly 3 GeV for (Ξ-+Ξ--)/2 and nearly 2.5 GeV for (Ω+Ω-)/2.Figure 7
Transverse momentum spectra for(Ξ-+Ξ--)/2 (a) and (Ω+Ω-)/2 (b) are compared to hydrodynamic models EPOS, VISH2+1, and KRAKOW for the (0–10)% most central Pb-Pb collisions at sNN=2.76 TeV.
(a)
(b)In the first case the KRAKOW model results are seen to fit the experimental data reasonably well but provides a bad fit for the second case, that is,(Ω+Ω-)/2. The VISH2+1 and HKM describe the shape of these spectra somewhat better; however, the VISH2+1 overestimates the yield of (Ω+Ω-)/2 at larger pT. Also it is found [18] that Krakow model overestimates the yield of protons up to 1.75 GeV. In contrast, our model successfully reproduces the shape of the transverse momentum distributions of these hadrons in a wider pT range up to 5 GeV.


## 4. Conclusion

The transverse momentum spectra of the hadronsp, p-, K+, K-, Ks0, Λ, Ω, Ω-, Ξ- and Ξ-- are fitted quite well by using our model. The assumption of vanishing chemical potential at midrapidity shows the effects of almost complete transparency in Pb-Pb collisions at LHC energy of 2.76 TeV. We also observe an earlier freeze-out of hyperons as compared to lighter mass particles, that is, Kaons and protons. The protons, antiprotons, Kaons, and anti-Kaons have similar freeze-out conditions, which indicate their near simultaneous freeze-out from the dense hadronic medium. The larger values of βT0 at the LHC energy, as compared to those at RHIC, indicates a stronger flow effect present in the system at LHC. The spectra are compared with the predictions from some other hydrodynamic models also and it is found that a better fit is obtained by using our model covering a wider pT range up to 5 GeV.



---

*Source: 154853-2015-03-25.xml*