# Model Comparison of the Transverse Momentum Spectra of Charged Hadrons Produced inPbPb Collision at sNN=5.02 TeV

Authors: Rohit Gupta; Satyajit Jena
Journal: Advances in High Energy Physics
 (2022)

---

## Abstract

Transverse momentum,pT, spectra are of prime importance in order to extract crucial information about the evolution dynamics of the system of particles produced in the collider experiments. In this work, the transverse momentum spectra of charged hadrons produced in PbPb collision at 5.02 TeV have been analyzed using different distribution functions in order to gain strong insight into the information that can be extracted from the spectra. We have also discussed the applicability of the unified statistical framework on the spectra of charged hadron at 5.02 TeV


---

## Body

## 1. Introduction

The quest to develop the understanding of Quantum Chromodynamics (QCD) matter created during the very early universe has been the primary motivation behind several theoretical as well as experimental studies in particle physics. Whether the QCD phase transition is first order, second order, or simple crossover and what are the critical temperature and the order parameters of phase transition, etc., are some of the long-standing questions that intrigue the researchers to explore QCD phase diagram by utilizing data from different heavy-ion collider experiments like Relativistic Heavy Ion Collider (RHIC) at BNL and in Large Hadron Collider (LHC) at CERN. The QCD matter under discussion is popularly known as Quark-Gluon Plasma (QGP) which cooks up under extreme conditions of temperature and energy density. When the temperatureT and the energy density ε reach a critical value, quarks no longer remain confined inside their individual hadron and instead become free to move over a nuclear volume forming a deconfined QGP state. Such extreme condition has been achieved by colliding heavy ions moving at the ultrarelativistic speed in the experiments at RHIC and LHC and provided us the tool to investigate and to uncover the mystery of the fireball that filled our universe few microseconds after the Big Bang.The theoretical reasoning behind the formation of QGP comes from the asymptotic freedom due to the nature of the running coupling strength of constituent partons. The theory of quark gluon interaction, QCD, dictates that the coupling strength between quarks and gluons increases with a decrease in the energy scales (or increase in length scale). On the other hand, the length scale can be decreased by bringing partons closer to each other. As the length scale reduces, the coupling strength also decreases, and it reaches a point where quarks and gluons asymptotically appear to be free from their nucleonic volume changing its hadronic degrees of freedom to partonic degrees of freedom. Hence, a dense hot soup of quarks and gluons forms, which we call QGP. This is exactly what experiments like RHIC and LHC had achieved by colliding highly energetic nucleons travelling at ultrarelativistic speed. However, the time scale of such process is extremely small so it is not possible to directly probe and characterize it using present day technologies. In experiments, we only receive the information about kinematics of final state particles that are free streaming to the detectors. We utilize the kinematics of these final state particles to get information about the underlying process that led to the formation of these particles.The transverse momentum,pT, spectra of final state particle are a crucial kinematic observable that can be utilized as an effective probe to understand the thermodynamical properties of the system produced in the collision. It is also an essential observable to understand the dynamics of QGP and the quark-hadron phase transition. Several theoretical models [1–4] have been developed to characterize pT spectra and to extract different physical parameters that can be further used to enhance our understanding of the system produced during such collisions. Most of these models are developed based on the statistical and thermodynamical approach since it is almost impossible to apply perturbative field theory calculations because of the high QCD coupling strength at low pT (due to asymptotic freedom).In this work, we have performed a comparative study of different models onpT spectra of charged hadrons produced in PbPb collision at the recently published result at sNN=5.02 TeV [5] along with the 2.76 TeV [6] data measured by the ALICE experiment at CERN.


## 2. Models

As discussed in the previous section, it is difficult to apply perturbative QCD at low-pT region of spectra as the coupling strength becomes extremely high at a low momentum transfer scale. Hence, we rely on the phenomenological models. Among several such models, statistical thermal approaches are widely used to explain transverse momentum spectra. Initial work by Koppe [7, 8], Fermi [9, 10], and Hagedorn [11, 12] developed the fundamentals of the statistical description of particle spectra in high energy physics. These initial works were followed by several applications of statistical thermodynamics in studying the particle production mechanism using the momenta carried out by outgoing particles. In this section, we will discuss different models developed to explain the spectra, and we will present the results obtained by fitting the pT spectra data produced in 5.02 TeV collision. Some of the models discussed below require the information of mass of the constituent particle. Since a larger contribution of charged hadrons comes from pion (~70%), we have used the mass of pion while analyzing the spectra of charged hadrons.


### 2.1. Boltzmann Distribution

For a system that is considered to be of purely thermal origin, the most natural choice to explain the statistical distribution of the energy of constituent particles is the standard Boltzmann-Gibbs statistics. Boltzmann distribution has long been used to explain energy spectra [1, 13] of the classical statistical system. Because the temperature of the system during heavy-ion collision is extremely high and Fermi-Dirac and Bose-Einstein statistics approach Maxwell-Boltzmann statistics at high temperature, we can use Boltzmann distribution to explain the spectra of particles produced during the heavy-ion collision. In the case of Boltzmann distribution, we can write number density as [1, 13]
(1)d2N2πpTdpTdy=mTgV2π3exp−mTT.Here,mT is the transverse mass and is related to transverse momentum by mT=m2+pT2, g is the spin degeneracy factor, and V is volume of the system under consideration. In Figure 1, we have fitted the pT spectra of charged hadrons that produced PbPb collision at 5.02 TeV with the Boltzmann distribution function (Equation (1)). We observe that the data deviate significantly from the BG function at low and high pT regions of the spectra. This deviation finds its origin in the low number of particles in the collision. BG statistics apply to the system where the number of particles in the system should be of the order of the Avogadro number; however, only a few thousand particles are produced in heavy-ion collision. At the same time, the applicability of BG statistics is limited to the system where entropy of the ensemble is extensive as well as additive in nature. To tackle these fundamental issues, we need to find a solution for nonextensive formalism which can be applied to the system under discussion. Therefore, a generalization of BG statistics was proposed by Tsallis in his seminal work in Ref. [2] to tackle nonextensivity in the system. By doing so, it overcomes the limitation of BG statistics. In this context, we provide a brief discussion on Tsallis statistics and the corresponding fit result of PbPb collision at sNN=5.02 TeV.Figure 1
The transverse momentum spectra of charged hadrons at different centralities produced at collision energies of 5.02 TeV [5] fitted with Boltzmann (Equation (1)).


### 2.2. Tsallis Distribution

To provide a generalized formalism to explain nonextensive systems such as the system produced in the heavy-ion collision, Tsallis statistics, also known as nonextensive statistics, were proposed as a generalization of Boltzmann-Gibbs statistics in 1988 by C. Tsallis. In Tsallis statistics, a parameter′q′ has been introduced, which takes care of the nonextensivity in the system. This new parameter also acts as a scaling factor for the number of particles to make this statistic applicable to the system with a low number of particles compared to the Avogadro number. Here, normal exponential is replaced by q-exponential as
(2)expqx=1−q−1x−1/q−1,which in the limit q⟶1 gives us normal exponential.Further, entropy in the case of Tsallis statistics [2] is given as
(3)Sq=k1−∑ipiqq−1.The functional form of transverse momentum spectra in the case of Tsallis statistics can be obtained by replacing normal exponential in Equation (1) with q-exponential:
(4)12πpTd2NdpTdy=gVmT2π31+q−1mTT−q/q−1.This distribution function converges to Boltzmann distribution in the limitq⟶1. Also, Tsallis statistics have been proved to be thermodynamically consistent following the laws of thermodynamics [14]. This form of transverse momentum spectra has been shown in different works [14–18] to nicely explain the spectra in a limited pT region. Additionally, it has been observed that there are certain systems [19] with significant long-range interactions; the entropy in such system can be nonadditive and nonextensive, and the Tsallis formalism has been used to explain such system [20].Figure2 shows the fitting of Tsallis distribution to the pT spectra of charged particles produced in PbPb collision at sNN=5.02 TeV. It is evident from the figure that Tsallis distribution provides a much better fit to the spectra compared to the Boltzmann distribution. The detailed comparison of χ2/NDF values is given in Table 1; we can observe that there is a significant improvement in explaining the spectra of PbPb collision at sNN=5.02 TeV.Figure 2
The transverse momentum spectra of charged hadrons at different centralities produced at collision energies of 5.02 TeV [5] fitted with Tsallis (Equation (4)).Table 1
Table of parameter values obtained after fitting charged hadron spectra at 3 different centralities with different functions.
VariableMethodEnergy (TeV)0 to 5%40 to 50%70 to 80%TeffMeVBG2.76316.95±2.74316.3±3.58305.23±5.755.02319.747±0.67330.97±0.67338.522±0.68Tsallis2.76163±3.995137.52±3.839111.22±4.1135.02176.101±1.049148.433±0.775117.75±0.640Unified function2.76393.535±44.567311.339±77.432296.085±215.7355.02407.448±3.322369.097±18.529329.4±34.236TkinMeVBW2.76124.25±67.877163.18±1.700157.73±11.5915.02113.571±14.841157.637±4.268164.249±1.847TBW2.7676.983±5.93952.006±6.02621.143±20.1565.02—76.292±1.06244.363±1.163qTsallis2.761.0945±0.00271.1186±0.00291.138±0.00345.021.09736±0.000721.12023±0.000571.1416±0.00052TBW2.761.0135±0.01021.0548±0.00881.0967±0.01785.02—1.04678±0.002311.08609±0.00147q-Weibull2.761.0212±0.01411.0613±0.01661.0905±0.02185.021.00274±0.004091.04755±0.003461.08525±0.00359Unified function2.761.0478±0.00381.0935±0.00841.1325±0.03125.021.04808±0.000181.08476±0.001771.1315±0.00317χ2/NDFBG2.7625.34529.81924.09075.02396.783735.2721080.25Tsallis2.761.9941.12260.510615.0234.300336.368723.6161BW2.760.2960.6030.5255.022.0872112.045220.6522TBW2.760.6660.5150.3515.02—4.748031.74155q-Weibull2.760.110.0590.0185.021.416062.193111.15846Unified function2.760.1010.0520.0175.021.709141.964821.08821However, we observed that it does not explain the data at low and highpT part of spectra. This might be due to the various particle production mechanisms in the heavy-ion collision, affecting the nature of spectra. Therefore, the applicability of Tsallis statistics is only limited to soft pT regimes of the spectra (although it deviates narrowly at low pT) whereas for high pT, dominated by hard processes, it does not explain. As shown and discussed in Ref. [21–23], the Tsallis formalism becomes extremely complex to incorporate hard processes, and eventually, a modification is required to explain very high pT part of the particle spectra. On the contrary, the particle production is dominated by hard QCD processes at the high pT region of the spectra, which could be explained by a QCD-inspired power law form of function such as
(5)fpT=1NdNdpT=ApT1+pTp0−n.A unified model to explain the contribution of both soft and hard processes topT spectra is still an open problem. Before we discuss the unified formalism [24, 25], we would like to discuss several other approaches to explain pT spectra. Apart from the purely statistical models discussed above, several hydrodynamics-based models are also developed to explain the spectra considering the initial state fluctuation and the collective behaviour of the QCD matter created during the collision.


### 2.3. Blast-Wave (BW)

It is experimentally established that the system produced during the high energy collision has an azimuthal anisotropy because of the difference in flow velocities along different directions. This azimuthal anisotropy is a result of some initial state geometrical effects that arise during the collision. Thus, it can be argued that the outgoing particles must carry the imprint of such effects and can influence the nature of the spectra. To incorporate such effects on spectra, several models are proposed [3, 4]. The Blast-Wave model is a hydrodynamics-inspired model developed to include the flow properties along with the random thermal motion of particles in order to give a complete picture of the evolution dynamics of QGP, including the azimuthal anisotropy. Transverse momentum spectra in the case of the Blast-Wave model [3] are given as
(6)dNpTdpT∝∫0RrdrmTI0pTsinhρrTkin×K1mTcoshρrTkin.Here,mT is the transverse mass and I0 and K are modified Bessel functions. Also, ρr=tanh−1β and β is the transverse radial flow velocity which is parameterized in the form of a power law of the form βtr=βsr/Rn, where βs represent the surface flow velocity and n is the exponent of the flow profile. The average transverse flow velocity βT is given in terms of βs and exponent n as βT=βs×2/2+n [3]. Finally, r/R represents the radial position of the thermal source.The temperature extracted using the BW function has a different interpretation as compared to the temperature extracted from the Boltzmann-Gibbs function (Equation (1)) or Tsallis function (Equation (4)). Fitting with BW function gives us kinetic freeze-out temperature Tkin whereas fitting with BG or Tsallis distribution gives us effective temperature Teff. These two quantities are related as [26]
(7)Teff=Tkin+fβt.Here,fβt is a function of transverse flow velocity βt.We show the result of Blast-Wave fitting to thepT spectra of particles produced in 5.02 TeV collision in Figure 3. The fitting yields the kinetic temperature as 113.6 MeV and the average flow velocity close to 0.63c for most central collision. The χ2/NDF values are close to unity for central collisions, but they start to increase with the decrease in the nuclear overlap by approaching a value of 20 for most peripheral collisions. This formalism does not include the concept of nonextensive; therefore, a generalization of the Blast-Wave model to include a nonextensive system is discussed in the next section.Figure 3
The transverse momentum spectra of charged hadrons at different centralities produced at collision energies of 5.02 TeV [5] fitted with Blast-Wave function (Equation (6)).


### 2.4. Tsallis Blast-Wave (TBW)

The Tsallis Blast-Wave method is an extension of the Blast-Wave method which was introduced to take into account the effect of nonextensivity in the system along with the flow properties. As discussed elsewhere in the article, Tsallis formalism can be used to tackle the system with nonextensivity; the Boltzmann distribution used in the Blast-Wave model is replaced by the Tsallis distribution in order to get the Tsallis Blast-Wave function [3, 4]. Transverse momentum spectra in the case of TBW are given as
(8)dNpTdpT∝mT∫−YYcoshydy∫−ππdϕ∫0Rrdr×1+q−1TmTcoshycoshρ−pTsinhρcosϕ−1/q−1.Here,y represent rapidity and ρ is the flow profile as described in Blast-Wave (BW).After its introduction, the TBW method has been used in different works to fitpT spectra of particles and extract information related to the thermodynamics and hydrodynamics evolution of the system produced in the heavy-ion collision. For TBW fit, the value of n is set to one [3, 4].We have shown the fitting ofpT spectra with TBW function in Figure 4. The value of kinetic freeze-out temperature extracted from the TBW fit shows a declining trend with the increase in centrality, and the average flow velocity is around 0.4c. The TBW model gives a holistic view of particle production in high energy collision, including statistical and hydrodynamical description but only at the low pT part of the spectra. Although the result obtained by this method gives a better explanation at low pT, it does not take care of hardness in the spectra in the high pT region, which is dominated by the particles produced in hard scattering processes.Figure 4
The transverse momentum spectra of charged hadrons at different centralities produced at collision energies of 5.02 TeV [5] fitted with Tsallis Blast-Wave function (Equation (8)).We also tried to fit this data with a statistical model as proposed in Ref. [27]. The q-Weibull distribution is an extension of Weibull distribution, which was described by Swedish mathematician Waloddi Weibull in 1951. Weibull distribution is a continuous probability distribution and is given as
(9)Px,λ,k=kλxλk−1e−x/λk,x≥0,0,x<0,where λ is the scale parameter and k is the shape parameter of distribution, with k&λ>0.Weibull distribution finds its application in the system where dynamical evolution is driven by fragmentation and sequential branching [28, 29]. Since the evolution of the system in hadrons and heavy-ion collision is dominated by a perturbative QCD-based parton cascade model, we can apply q-Weibull distribution to study particle spectra.Incorporating the Tsallis framework in Weibull distribution gives us theq-Weibull distribution [27]:
(10)Pqx,q,λ,k=kλxλk−1eq−x/λk,where
(11)eq−x/λk=1−1−qxλk1/1−q.In Figure5, we have performed q-Weibull fit to the transverse momentum spectra data, and the fit shows good agreement with the data. Its relevance in particle spectra lies in its ability to fit spectra over a wide range of pT as compared to the other models described above with good accuracy. The models discussed so far mostly described the data at low pT quite well as they do not include the hard QCD processes in their construction. In the next section, we will discuss a unified formalism that describes both soft and hard processes in a consistent manner.Figure 5
The transverse momentum spectra of charged hadrons at different centralities produced at collision energies of 5.02 TeV [5] fitted with q-Weibull function (Equation (10)).


### 2.5. Unified Distribution Function

Pearson distribution is a generalized probability distribution, which, under different limits on its parameter, reduces to different distribution functions like exponential, Gaussian, Gamma distribution, and Student’st-distribution. It was introduced in 1895 by Karl Pearson in his seminal work [30]. It has been applied quite successfully in different fields such as geophysics, statistics, and financial marketing. It has been introduced in the context of particle production in heavy-ion collision for the first time in Ref. [24].The transverse momentum spectra in the case of unified distribution are given as(12)12πpTd2NdpTdy=B′1+pTp0−n1+q−1pTT−q/q−1.As it has been described in Ref. [24, 25], this form of pT distribution function nicely takes care of both the “hard” and “soft” parts of the spectra. Further, the unified formalism also provides a connection to the flow of particle as discussed in Ref. [24].Figure6 represents a good agreement between transverse momentum spectra and the unified distribution function in the form of Equation (12). Apart from providing the best fit to the transverse momentum spectra (as is evident from the χ2/NDF values provided in Table 1), unified distribution is also more physical in the sense that it is thermodynamically consistent [25], and it is also proved to be backward compatible with Tsallis distribution under limiting conditions.Figure 6
The transverse momentum spectra of charged hadrons at different centralities produced at collision energies of 5.02 TeV [5] fitted with unified distribution function.


## 2.1. Boltzmann Distribution

For a system that is considered to be of purely thermal origin, the most natural choice to explain the statistical distribution of the energy of constituent particles is the standard Boltzmann-Gibbs statistics. Boltzmann distribution has long been used to explain energy spectra [1, 13] of the classical statistical system. Because the temperature of the system during heavy-ion collision is extremely high and Fermi-Dirac and Bose-Einstein statistics approach Maxwell-Boltzmann statistics at high temperature, we can use Boltzmann distribution to explain the spectra of particles produced during the heavy-ion collision. In the case of Boltzmann distribution, we can write number density as [1, 13]
(1)d2N2πpTdpTdy=mTgV2π3exp−mTT.Here,mT is the transverse mass and is related to transverse momentum by mT=m2+pT2, g is the spin degeneracy factor, and V is volume of the system under consideration. In Figure 1, we have fitted the pT spectra of charged hadrons that produced PbPb collision at 5.02 TeV with the Boltzmann distribution function (Equation (1)). We observe that the data deviate significantly from the BG function at low and high pT regions of the spectra. This deviation finds its origin in the low number of particles in the collision. BG statistics apply to the system where the number of particles in the system should be of the order of the Avogadro number; however, only a few thousand particles are produced in heavy-ion collision. At the same time, the applicability of BG statistics is limited to the system where entropy of the ensemble is extensive as well as additive in nature. To tackle these fundamental issues, we need to find a solution for nonextensive formalism which can be applied to the system under discussion. Therefore, a generalization of BG statistics was proposed by Tsallis in his seminal work in Ref. [2] to tackle nonextensivity in the system. By doing so, it overcomes the limitation of BG statistics. In this context, we provide a brief discussion on Tsallis statistics and the corresponding fit result of PbPb collision at sNN=5.02 TeV.Figure 1
The transverse momentum spectra of charged hadrons at different centralities produced at collision energies of 5.02 TeV [5] fitted with Boltzmann (Equation (1)).


## 2.2. Tsallis Distribution

To provide a generalized formalism to explain nonextensive systems such as the system produced in the heavy-ion collision, Tsallis statistics, also known as nonextensive statistics, were proposed as a generalization of Boltzmann-Gibbs statistics in 1988 by C. Tsallis. In Tsallis statistics, a parameter′q′ has been introduced, which takes care of the nonextensivity in the system. This new parameter also acts as a scaling factor for the number of particles to make this statistic applicable to the system with a low number of particles compared to the Avogadro number. Here, normal exponential is replaced by q-exponential as
(2)expqx=1−q−1x−1/q−1,which in the limit q⟶1 gives us normal exponential.Further, entropy in the case of Tsallis statistics [2] is given as
(3)Sq=k1−∑ipiqq−1.The functional form of transverse momentum spectra in the case of Tsallis statistics can be obtained by replacing normal exponential in Equation (1) with q-exponential:
(4)12πpTd2NdpTdy=gVmT2π31+q−1mTT−q/q−1.This distribution function converges to Boltzmann distribution in the limitq⟶1. Also, Tsallis statistics have been proved to be thermodynamically consistent following the laws of thermodynamics [14]. This form of transverse momentum spectra has been shown in different works [14–18] to nicely explain the spectra in a limited pT region. Additionally, it has been observed that there are certain systems [19] with significant long-range interactions; the entropy in such system can be nonadditive and nonextensive, and the Tsallis formalism has been used to explain such system [20].Figure2 shows the fitting of Tsallis distribution to the pT spectra of charged particles produced in PbPb collision at sNN=5.02 TeV. It is evident from the figure that Tsallis distribution provides a much better fit to the spectra compared to the Boltzmann distribution. The detailed comparison of χ2/NDF values is given in Table 1; we can observe that there is a significant improvement in explaining the spectra of PbPb collision at sNN=5.02 TeV.Figure 2
The transverse momentum spectra of charged hadrons at different centralities produced at collision energies of 5.02 TeV [5] fitted with Tsallis (Equation (4)).Table 1
Table of parameter values obtained after fitting charged hadron spectra at 3 different centralities with different functions.
VariableMethodEnergy (TeV)0 to 5%40 to 50%70 to 80%TeffMeVBG2.76316.95±2.74316.3±3.58305.23±5.755.02319.747±0.67330.97±0.67338.522±0.68Tsallis2.76163±3.995137.52±3.839111.22±4.1135.02176.101±1.049148.433±0.775117.75±0.640Unified function2.76393.535±44.567311.339±77.432296.085±215.7355.02407.448±3.322369.097±18.529329.4±34.236TkinMeVBW2.76124.25±67.877163.18±1.700157.73±11.5915.02113.571±14.841157.637±4.268164.249±1.847TBW2.7676.983±5.93952.006±6.02621.143±20.1565.02—76.292±1.06244.363±1.163qTsallis2.761.0945±0.00271.1186±0.00291.138±0.00345.021.09736±0.000721.12023±0.000571.1416±0.00052TBW2.761.0135±0.01021.0548±0.00881.0967±0.01785.02—1.04678±0.002311.08609±0.00147q-Weibull2.761.0212±0.01411.0613±0.01661.0905±0.02185.021.00274±0.004091.04755±0.003461.08525±0.00359Unified function2.761.0478±0.00381.0935±0.00841.1325±0.03125.021.04808±0.000181.08476±0.001771.1315±0.00317χ2/NDFBG2.7625.34529.81924.09075.02396.783735.2721080.25Tsallis2.761.9941.12260.510615.0234.300336.368723.6161BW2.760.2960.6030.5255.022.0872112.045220.6522TBW2.760.6660.5150.3515.02—4.748031.74155q-Weibull2.760.110.0590.0185.021.416062.193111.15846Unified function2.760.1010.0520.0175.021.709141.964821.08821However, we observed that it does not explain the data at low and highpT part of spectra. This might be due to the various particle production mechanisms in the heavy-ion collision, affecting the nature of spectra. Therefore, the applicability of Tsallis statistics is only limited to soft pT regimes of the spectra (although it deviates narrowly at low pT) whereas for high pT, dominated by hard processes, it does not explain. As shown and discussed in Ref. [21–23], the Tsallis formalism becomes extremely complex to incorporate hard processes, and eventually, a modification is required to explain very high pT part of the particle spectra. On the contrary, the particle production is dominated by hard QCD processes at the high pT region of the spectra, which could be explained by a QCD-inspired power law form of function such as
(5)fpT=1NdNdpT=ApT1+pTp0−n.A unified model to explain the contribution of both soft and hard processes topT spectra is still an open problem. Before we discuss the unified formalism [24, 25], we would like to discuss several other approaches to explain pT spectra. Apart from the purely statistical models discussed above, several hydrodynamics-based models are also developed to explain the spectra considering the initial state fluctuation and the collective behaviour of the QCD matter created during the collision.


## 2.3. Blast-Wave (BW)

It is experimentally established that the system produced during the high energy collision has an azimuthal anisotropy because of the difference in flow velocities along different directions. This azimuthal anisotropy is a result of some initial state geometrical effects that arise during the collision. Thus, it can be argued that the outgoing particles must carry the imprint of such effects and can influence the nature of the spectra. To incorporate such effects on spectra, several models are proposed [3, 4]. The Blast-Wave model is a hydrodynamics-inspired model developed to include the flow properties along with the random thermal motion of particles in order to give a complete picture of the evolution dynamics of QGP, including the azimuthal anisotropy. Transverse momentum spectra in the case of the Blast-Wave model [3] are given as
(6)dNpTdpT∝∫0RrdrmTI0pTsinhρrTkin×K1mTcoshρrTkin.Here,mT is the transverse mass and I0 and K are modified Bessel functions. Also, ρr=tanh−1β and β is the transverse radial flow velocity which is parameterized in the form of a power law of the form βtr=βsr/Rn, where βs represent the surface flow velocity and n is the exponent of the flow profile. The average transverse flow velocity βT is given in terms of βs and exponent n as βT=βs×2/2+n [3]. Finally, r/R represents the radial position of the thermal source.The temperature extracted using the BW function has a different interpretation as compared to the temperature extracted from the Boltzmann-Gibbs function (Equation (1)) or Tsallis function (Equation (4)). Fitting with BW function gives us kinetic freeze-out temperature Tkin whereas fitting with BG or Tsallis distribution gives us effective temperature Teff. These two quantities are related as [26]
(7)Teff=Tkin+fβt.Here,fβt is a function of transverse flow velocity βt.We show the result of Blast-Wave fitting to thepT spectra of particles produced in 5.02 TeV collision in Figure 3. The fitting yields the kinetic temperature as 113.6 MeV and the average flow velocity close to 0.63c for most central collision. The χ2/NDF values are close to unity for central collisions, but they start to increase with the decrease in the nuclear overlap by approaching a value of 20 for most peripheral collisions. This formalism does not include the concept of nonextensive; therefore, a generalization of the Blast-Wave model to include a nonextensive system is discussed in the next section.Figure 3
The transverse momentum spectra of charged hadrons at different centralities produced at collision energies of 5.02 TeV [5] fitted with Blast-Wave function (Equation (6)).


## 2.4. Tsallis Blast-Wave (TBW)

The Tsallis Blast-Wave method is an extension of the Blast-Wave method which was introduced to take into account the effect of nonextensivity in the system along with the flow properties. As discussed elsewhere in the article, Tsallis formalism can be used to tackle the system with nonextensivity; the Boltzmann distribution used in the Blast-Wave model is replaced by the Tsallis distribution in order to get the Tsallis Blast-Wave function [3, 4]. Transverse momentum spectra in the case of TBW are given as
(8)dNpTdpT∝mT∫−YYcoshydy∫−ππdϕ∫0Rrdr×1+q−1TmTcoshycoshρ−pTsinhρcosϕ−1/q−1.Here,y represent rapidity and ρ is the flow profile as described in Blast-Wave (BW).After its introduction, the TBW method has been used in different works to fitpT spectra of particles and extract information related to the thermodynamics and hydrodynamics evolution of the system produced in the heavy-ion collision. For TBW fit, the value of n is set to one [3, 4].We have shown the fitting ofpT spectra with TBW function in Figure 4. The value of kinetic freeze-out temperature extracted from the TBW fit shows a declining trend with the increase in centrality, and the average flow velocity is around 0.4c. The TBW model gives a holistic view of particle production in high energy collision, including statistical and hydrodynamical description but only at the low pT part of the spectra. Although the result obtained by this method gives a better explanation at low pT, it does not take care of hardness in the spectra in the high pT region, which is dominated by the particles produced in hard scattering processes.Figure 4
The transverse momentum spectra of charged hadrons at different centralities produced at collision energies of 5.02 TeV [5] fitted with Tsallis Blast-Wave function (Equation (8)).We also tried to fit this data with a statistical model as proposed in Ref. [27]. The q-Weibull distribution is an extension of Weibull distribution, which was described by Swedish mathematician Waloddi Weibull in 1951. Weibull distribution is a continuous probability distribution and is given as
(9)Px,λ,k=kλxλk−1e−x/λk,x≥0,0,x<0,where λ is the scale parameter and k is the shape parameter of distribution, with k&λ>0.Weibull distribution finds its application in the system where dynamical evolution is driven by fragmentation and sequential branching [28, 29]. Since the evolution of the system in hadrons and heavy-ion collision is dominated by a perturbative QCD-based parton cascade model, we can apply q-Weibull distribution to study particle spectra.Incorporating the Tsallis framework in Weibull distribution gives us theq-Weibull distribution [27]:
(10)Pqx,q,λ,k=kλxλk−1eq−x/λk,where
(11)eq−x/λk=1−1−qxλk1/1−q.In Figure5, we have performed q-Weibull fit to the transverse momentum spectra data, and the fit shows good agreement with the data. Its relevance in particle spectra lies in its ability to fit spectra over a wide range of pT as compared to the other models described above with good accuracy. The models discussed so far mostly described the data at low pT quite well as they do not include the hard QCD processes in their construction. In the next section, we will discuss a unified formalism that describes both soft and hard processes in a consistent manner.Figure 5
The transverse momentum spectra of charged hadrons at different centralities produced at collision energies of 5.02 TeV [5] fitted with q-Weibull function (Equation (10)).


## 2.5. Unified Distribution Function

Pearson distribution is a generalized probability distribution, which, under different limits on its parameter, reduces to different distribution functions like exponential, Gaussian, Gamma distribution, and Student’st-distribution. It was introduced in 1895 by Karl Pearson in his seminal work [30]. It has been applied quite successfully in different fields such as geophysics, statistics, and financial marketing. It has been introduced in the context of particle production in heavy-ion collision for the first time in Ref. [24].The transverse momentum spectra in the case of unified distribution are given as(12)12πpTd2NdpTdy=B′1+pTp0−n1+q−1pTT−q/q−1.As it has been described in Ref. [24, 25], this form of pT distribution function nicely takes care of both the “hard” and “soft” parts of the spectra. Further, the unified formalism also provides a connection to the flow of particle as discussed in Ref. [24].Figure6 represents a good agreement between transverse momentum spectra and the unified distribution function in the form of Equation (12). Apart from providing the best fit to the transverse momentum spectra (as is evident from the χ2/NDF values provided in Table 1), unified distribution is also more physical in the sense that it is thermodynamically consistent [25], and it is also proved to be backward compatible with Tsallis distribution under limiting conditions.Figure 6
The transverse momentum spectra of charged hadrons at different centralities produced at collision energies of 5.02 TeV [5] fitted with unified distribution function.


## 3. Result

We have performed the fitting of transverse momentum spectra using differentpT distribution functions like the Boltzmann-Gibbs function, Tsallis distribution function, Blast-Wave, Tsallis Blast-Wave, and q-Weibull and unified framework. We have used the recently released pT spectra data of charged hadrons produced in PbPb collision at 5.02 TeV measured by ALICE experiment [5]. For comparison, we have also performed a similar analysis on the charged hadrons pT spectra data produced in 2.76 TeV PbPb collision [6].In the above section, we have provided the fitting plot for5.02 TeV pT spectra data. Fitting has been carried out for different collision centrality classes including 0‐5% (most central collision), 40‐50% (midcentral collision), and 70‐80% (peripheral collision). We have used the ROOT [31] data analysis framework to perform the fittings.In Tables1 and 2 along with Figures 7–10, best fit values of different parameters are presented along with the χ2/NDF value which is a measure of the goodness of fit.Table 2
The best fit value of exponent “n” and average transverse flow velocity obtained by fitting the charged hadron transverse momentum spectra using BW and TBW models.
Centralityn (BW)βT (BW)βT (TBW)2.76 TeV5.02 TeV2.76 TeV5.02 TeV2.76 TeV5.02 TeV0 to 5%1.2633±0.77370.9471±0.12360.5544±0.14420.6303±0.02910.4241±0.0048—5 to 10%1.3660±0.61370.9538±0.13220.534±0.10590.6316±0.03080.4244±0.00480.4181±0.001410 to 20%1.5316±1.64410.9740±0.14650.5047±0.25280.6302±0.03340.4247±0.00520.4192±0.001320 to 30%1.8038±0.55461.1591±0.17040.4642±0.07140.5862±0.03380.4254±0.00520.4186±0.001430 to 40%2.1756±0.46731.4738±0.20040.4199±0.04870.5225±0.03190.4263±0.00540.4199±0.001240 to 50%2.6499±0.46991.9766±0.08140.3771±0.03890.4479±0.00950.4281±0.00480.4189±0.001450 to 60%3.1212±0.49932.4192±0.07990.3446±0.03410.4033±0.00750.4305±0.00510.4200±0.001660 to 70%3.6574±0.55023.2153±0.08770.3149±0.03090.3424±0.00540.4330±0.00520.4209±0.001670 to 80%3.9250±0.59073.4730±0.07970.3041±0.03050.3301±0.00480.4355±0.00830.4257±0.0013Figure 7
Value ofχ2/NDF obtained for different functions fitted with pT spectra data of particles produced at 2.76PbPb collision.Figure 8
Value ofχ2/NDF obtained for different functions fitted with pT spectra data of particles produced at 5.02PbPb collision.Figure 9
Value of nonextensivity parameter “q” for different functions fitted with pT spectra data of particles produced at 2.76 TeV PbPb collision.Figure 10
Value of nonextensivity parameter “q” for different functions fitted with pT spectra data of particles produced at 5.02 TeV PbPb collision.Looking at theq values, we can predict the nonextensivity in the system, or, in other words, we can quantify how much a system deviates from the equilibrium. Large q values indicate a nonequilibrium system, and the system approaches equilibrium as the q value approaches unity.From Table1, q values from all four methods, Tsallis, TBW, q-Weibull, and unified function show a decreasing trend as we move from most peripheral to most central collision indicating that the system is highly nonequilibrium in the peripheral collision, and it moves toward equilibrium as the overlap region between two colliding nuclei increases. The trend in this result shows a similar pattern reported previously for different energies in Ref. [3, 4, 27].Table2 shows a decreasing trend in the value of average transverse flow velocity estimated by the BW model as we go from central to peripheral collision. The corresponding numerical values for TBW fit do not show any centrality dependence, and also, the values are lower as compared to the BW model. This difference can be attributed to the absence of the nonequilibrium description in the BW model as it demands the system to be in thermal equilibrium.Apart from the nonextensivity parameter, theq-Weibull fit also provides the value of two more free parameters λ and k. The best fit value of these parameters is presented in Table 3. The value of parameter k increases as we go from central to peripheral collision for both the energies. The values are also consistently higher in 2.76 TeV as compared to the 5.02 TeV. The parameter λ shows a reverse trend with the values decreasing from central to peripheral collision. The physics interpretation of these parameters is still an open problem, and a detailed study of q-Weibull function is required to have a better understanding of the model.Table 3
The best fit value of parametersk and λ obtained by fitting the charged hadron transverse momentum spectra using q-Weibull model.
Centralitykλ2.76 TeV5.02 TeV2.76 TeV5.02 TeV0 to 5%0.8183±0.04810.7666±0.01270.1953±0.01700.1974±0.00495 to 10%0.8218±0.04930.7741±0.01320.1949±0.01720.1988±0.004910 to 20%0.8297±0.05030.7788±0.01180.1940±0.01720.1972±0.004420 to 30%0.8407±0.05200.7998±0.01150.1910±0.01710.1968±0.004030 to 40%0.8545±0.05460.8035±0.01100.1865±0.01720.1894±0.003740 to 50%0.8684±0.05740.8295±0.01110.1791±0.01710.1853±0.003550 to 60%0.8816±0.06210.8492±0.01190.1703±0.01750.1772±0.003560 to 70%0.8971±0.06900.8717±0.01560.1611±0.01820.1678±0.003270 to 80%0.8909±0.07830.8662±0.01230.1489±0.01970.1536±0.0033Among all the different distribution functions, we have obtained the best fit for the unified distribution, which is also evident from theχ2/NDF values. One possible reason for the low χ2/NDF in unified distribution as compared to other distributions is that it also takes care of hard perturbative QCD processes, which affect the pT distribution at higher pT range as described in Ref. [24]. Another interesting result which is underdiscussed is related to the scaling properties of q parameter. Since the q parameter acts as a scaling factor to apply statistical mechanics at the low number of particles, the number of final state particles decreases with the increase in centrality, and hence, the corresponding scaling parameter should increase with centrality. We observe this in Figures 9 and 10 for all four different methods indicating that the q parameter can be interpreted as a scaling factor depending upon the number of final state particles in the system.


## 4. Conclusion

In this work, we have analyzed thepT spectra of charged particles produced in PbPb collision at sNN=5.02 TeV using six different formalisms which are based on different physics inputs; for example, Boltzmann distribution is a purely thermal model, Blast-Wave models also include collective flow effect, and the unified distribution also takes into account hard QCD processes. The effects such as nonextensivity and collective flow need to be taken into account in a model for explaining such data. From the variation of q-parameter with centrality, we also observe that the system is highly nonequilibrium at the peripheral collision and moves toward the equilibrium as we move toward the most central collision. Further, a very high value of χ2/NDF indicates that a purely thermal model is not a good explanation for the pT spectra. Apart from this, the scaling properties of the q parameter have also been established, indicating its relevance in small systems.In conclusion, we can say that a complete model to explain transverse momentum spectra must include the nonextensivity and collective flow on top of random thermal motion and should also be able to explain the effects on highpT range arising due to hard pQCD processes. The essence of all these physics inputs is present in the unified distribution, and it also provides the best fit to the transverse momentum spectra as is evident from the values of χ2/NDF presented in Figures 7 and 8.



---

*Source: 5482034-2022-06-14.xml*