# Indication of a Differential Freeze-Out in Proton-Proton and Heavy-Ion Collisions at RHIC and LHC Energies

Authors: Dhananjaya Thakur; Sushanta Tripathy; Prakhar Garg; Raghunath Sahoo; Jean Cleymans
Journal: Advances in High Energy Physics
 (2016)

---

## Abstract

The experimental data from the RHIC and LHC experiments of invariantpT spectra for most peripheral A+A and p+p collisions are analyzed with Tsallis distributions in different approaches. The information about the freeze-out surface in terms of freeze-out volume, temperature, chemical potential, and radial flow velocity for π+, K+, and p and their antiparticles is obtained. Furthermore, these parameters are studied as a function of the mass of the particles. A mass dependent differential freeze-out is observed which does not seem to distinguish between particles and their antiparticles. Furthermore, a mass-hierarchy in the radial flow is observed, meaning heavier particles suffer lower radial flow. Tsallis distribution function at finite chemical potential is used to study the mass dependence of chemical potential. The peripheral heavy-ion and proton-proton collisions at the same energies seem to be equivalent in terms of the extracted thermodynamic parameters.


---

## Body

## 1. Introduction

High-energy heavy-ion collisions provide a unique opportunity to study the nuclear matter under extreme conditions, that is, at high temperature and/or density. Due to high multiplicities produced inA+A and p+p collisions, the statistical models are more suitable to describe the particle production mechanism. Such a statistical description of transverse momentum pT of final state particles produced in high-energy collisions has been proposed to follow a thermalized Boltzmann type of distribution as given by [1](1)Ed3σd3p≃Cexp⁡-pTTexp⁡.To account for the high-pT tail, a power-law in pT has been proposed [2–4], which empirically accounts for the possible QCD contributions. Hagedorn proposed a combination of both the aspects, which describes the experimental data over a wide pT range [5] and is given by(2)Ed3σd3p=C1+pTp0-n⟶exp⁡-npTp0forpT⟶0,p0pTnforpT⟶∞,where C, p0, and n are fitting parameters. This becomes a purely exponential function for small pT and a purely power-law function for large pT values. A finite degree of deviation from the equilibrium statistical description of identified particle pT spectra has already been observed by experiments at RHIC [6, 7] and LHC [8–11]. Contrary to a thermalized system, where 〈pT〉 is associated with the temperature of the hadronizing matter, one fails to make such a connection in case of systems which are far from thermal equilibrium. In the latter systems, the temperature fluctuates either event by event or within the same event [12]. This creates room for possible description of the pT spectra in high-energy hadronic and nuclear collisions, using the nonextensive Tsallis statistics [13–15]. A thermodynamically consistent nonextensive distribution function is given by [16](3)fmT=Cq1+q-1mTT-1/q-1.Here, mT is the transverse mass and q is called the nonextensive parameter, a measure of degree of deviation from equilibrium. Equations (2) and (3) are related through the following transformations for large values of pT:(4)n=1q-1,p0=Tq-1.In the limitq→1, one recovers the standard Boltzmann-Gibbs distribution (see (1)) from the Tsallis distribution (see (6)). Here the effective kinetic freeze-out temperature (Texp⁡) as obtained from the inverse slope of the pT-spectra using (1) is related to the Tsallis temperature T by(5)Texp⁡≃a+b·T,where a=0.31-0.654q+0.354q2 and b=27.35-55q+29.07q2 and they are obtained numerically for distributions with same mean transverse momentum, pT, as discussed in [17] for p+p collisions in the low-pT (0.15<pT<0.6 GeV/c) regime.Tsallis statistics is used widely to explain the particle spectra in high-energy collisions [12, 22–27] starting from elementary e++e-, hadronic, and heavy-ion collisions [6, 28–44]. The produced particles from the collisions carry the information about collision dynamics and the subsequent space-time evolution till the occurrence of the final freeze-out. The evolution of the partonic system created in high-energy experiments is generally believed to be best described by hydrodynamics of an almost ideal fluid [45]. This approach gives a fair description of data on the transverse momentum spectra of hadrons, which are treated as one of the important tools to understand the production dynamics of particles in high-energy collisions. The systematic analysis with the help of an appropriate model or approach guides us to understand various thermodynamical as well as hydrodynamical properties of the fireball at different stages of its evolution. The integrated yields of various hadronic species at different center of mass energies are used in the present work. The corresponding freeze-out parameters for each hadronic species at the time of their freeze-out can be obtained from the analysis of their respective transverse momentum distributions. Different forms of the invariant yields using Tsallis distribution are available in the literature [24, 46, 47]. In the present work we have used all of these forms to study the temperature T, chemical potential μ, radial flow β, volume V, and nonextensive parameter q. It should be mentioned here that the parameter V is not necessarily related to the volume one obtains from HBT like experimental measurements. Further, we study the mass dependence of these parameters, which are obtained by analyzing invariant transverse momentum spectra. For the present analysis, we have used the data of p+p and A+A collisions of different experiments at RHIC and LHC. We observe a clear mass dependence of the above parameters, and the behaviour is found to be consistent from most peripheral A+A collisions to p+p collisions. The obtained thermodynamic parameters in p+p collisions are similar to those extracted for most peripheral A+A collisions at the same center of mass energies. This indicates a thermodynamical similarity between both the systems at a given collision energy.In the heavy-ion collision, the interaction volume of fireball decreases from most central to most peripheral collisions. So the number of participant nucleons also decreases from most central collisions to most peripheral collisions depending on the interaction volume. The system having more participants will quickly reach the equilibrium because of large number of binary collisions by rescattering of partons/hadrons as can be the case in centralA+A collisions. But in case of peripheral collisions due to smaller number of participants the system will be away from equilibrium for a while as compared to central collisions. Such a nonequilibrium system is better described by Tsallis nonextensive statistics, giving information about the various thermodynamic parameters of the system.The paper is organized as follows. In Section2, we present three forms of Tsallis distribution functions. Firstly, we discuss invariant yields with and without chemical potential in the Tsallis function. Then we show the Tsallis form of invariant yields with radial flow which is introduced analytically in one of our recent works [24]. In Section 3, results and discussions are made. Lastly, we conclude our findings in Section 4.


## 2. Nonextensivity and Transverse Momentum Spectra

In Sections2 and 3, we discuss the transverse momentum spectra of identified particles (π+, K+, and p and their antiparticles) produced in RHIC and LHC experiments using different forms of invariant yields using Tsallis nonextensive statistics.


### 2.1. Nonextensive Statistics without Radial Flow

The Tsallis distribution function at mid-rapidity, with finite chemical potential and without radial flow [46], is given by(6)1pTd2NdpTdyy=0=gVmT2π21+q-1mT-μT-q/q-1,where mT is the transverse mass of a particle given by pT2+m2, g is the degeneracy, and μ is the chemical potential of the system. In view of higher center of mass energies, where μ≃0, the transverse momentum distribution function [47] becomes(7)1pTd2NdpTdyy=0=gVmT2π21+q-1mTT-q/q-1.


### 2.2. Nonextensive Statistics with Radial Flow

The value of the nonextensive parameterq for high-energy collisions is 1≤q≤1.2 [48]. To study the order of deviation of the pT-spectra from an equilibrium Boltzmann distribution, the Tsallis distribution function has been expanded in a Taylor series in view of (q-1) being very small, after successful inclusion of radial flow in a relativistic scenario. The details of the method are described in [24]. The functional form of the distribution up to first order in (q-1) is given by (8)12πpTdNdpTdy=gV2π32TrI0sK1r-sI1sK0r-q-1Tr2I0sK0r+K2r+4q-1TrsI1sK1r-q-1Ts2K0rI0s+I2s+q-14Tr3I0sK3r+3K1r-3q-12Tr2sK2r+K0rI1s+3q-12Ts2rI0s+I2sK1r-q-14Ts3I3s+3I1sK0r,where r≡γmT/T,s≡γβpT/T.Ins and Knr are the modified Bessel functions of the first and second kind. There are four parameters involved, namely, V, T, β, and q, where V is the volume, T is the Tsallis temperature, β is the radial flow velocity, and q is the Tsallis nonextensive parameter.We use (8) to fit the particle spectra of identified particles in heavy-ion collisions to study the radial flow parameter.


## 2.1. Nonextensive Statistics without Radial Flow

The Tsallis distribution function at mid-rapidity, with finite chemical potential and without radial flow [46], is given by(6)1pTd2NdpTdyy=0=gVmT2π21+q-1mT-μT-q/q-1,where mT is the transverse mass of a particle given by pT2+m2, g is the degeneracy, and μ is the chemical potential of the system. In view of higher center of mass energies, where μ≃0, the transverse momentum distribution function [47] becomes(7)1pTd2NdpTdyy=0=gVmT2π21+q-1mTT-q/q-1.


## 2.2. Nonextensive Statistics with Radial Flow

The value of the nonextensive parameterq for high-energy collisions is 1≤q≤1.2 [48]. To study the order of deviation of the pT-spectra from an equilibrium Boltzmann distribution, the Tsallis distribution function has been expanded in a Taylor series in view of (q-1) being very small, after successful inclusion of radial flow in a relativistic scenario. The details of the method are described in [24]. The functional form of the distribution up to first order in (q-1) is given by (8)12πpTdNdpTdy=gV2π32TrI0sK1r-sI1sK0r-q-1Tr2I0sK0r+K2r+4q-1TrsI1sK1r-q-1Ts2K0rI0s+I2s+q-14Tr3I0sK3r+3K1r-3q-12Tr2sK2r+K0rI1s+3q-12Ts2rI0s+I2sK1r-q-14Ts3I3s+3I1sK0r,where r≡γmT/T,s≡γβpT/T.Ins and Knr are the modified Bessel functions of the first and second kind. There are four parameters involved, namely, V, T, β, and q, where V is the volume, T is the Tsallis temperature, β is the radial flow velocity, and q is the Tsallis nonextensive parameter.We use (8) to fit the particle spectra of identified particles in heavy-ion collisions to study the radial flow parameter.


## 3. Results and Discussion

It is expected that the number of binary collisions in a system with smaller number of participant nucleons is quite low. Hence, the probability of mutual interaction (resulting in momentum transfer) between the system quanta partons or hadrons becomes less for the systems with small value of participating nucleons. This makes the system stay away from a possible thermal equilibrium. On the other hand, an appreciable increment in the number of binary collisions is observed in a system possessing a large number of participant nucleons and consequently the system can reach quickly its thermal equilibrium or be in the close vicinity of it. It has been shown in [49] that the nonextensivity parameter (q) is close to 1 in central Au + Au collisions and increases towards peripheral collisions. A higher value of q reflects that the system is away from thermal equilibrium. These results suggest that the degree of nonequilibrium is higher in peripheral collisions. Such systems are best described by Tsallis distributions, for example, the most peripheral A+A and p+p systems. We perform the fitting of pT spectra for pT<3 GeV/c around mid-rapidity at sNN = 200 GeV and sNN = 2.76 TeV. Keeping all the parameters free, we try to fit the spectra with (6), (7), and (8) independently for different particles like π+, K+, and p and their antiparticle for Pb + Pb and Au + Au most peripheral collisions using TMinuit class available in the ROOT library [50] to get a convergent solution. The same procedure is followed for p+p collisions as well.We have discussed the present work in three sections. Now, (7) is used to analyze invariant pT spectra of identified particles (π+, K+, and p and their antiparticle) for most peripheral Au + Au collisions at sNN = 200 GeV and most peripheral Pb + Pb collisions at sNN = 2.76 TeV along with p+p collisions at both energies. The fitting of pT spectra of identified particles for Au + Au collisions at 200 GeV and Pb + Pb collisions at 2.76 TeV are shown in Figures 11 and 12 for most peripheral collisions. Also, Figures 13 and 14 show the fitting of invariant spectra to p+p collisions at 200 GeV and 2.76 TeV, respectively, for identified particles. The fitting is performed up to the maximum pT of 3 GeV/c for A+A collisions and 2.5 GeV/c for p+p collisions although the K+pT spectra in Au + Au collision are fitted only up to 2.0 GeV/c due to unavailability of data [18]. The fittings of Tsallis function with μ=0 seem to be better for peripheral heavy-ion collisions. χ2/ndf obtained by minimization in TMinuit are shown in Figures 1 and 2.Figure 1
χ 2 / ndf values are shown for most peripheral Au + Au and p+p collisions at sNN = 200 GeV using (7) as a fit function.Figure 2
χ 2 / ndf values are shown for most peripheral Pb + Pb and p+p collisions at sNN = 2.76 TeV using (7) as a fit function.A systematic study of the extracted parameters like volume (V), temperature (T), and Tsallis nonextensive parameter (q) is made with identified particle masses. The values of extracted parameters are summarized in Table 1. We have calculated the radius of the fireball at freeze-out and plotted it as a function of the mass of the identified particles. Figures 3 and 4 show the variation of common mass dependent parameters at sNN = 200 GeV and sNN = 2.76 TeV. We have then compared the parameters obtained for peripheral heavy-ion collisions with that of p+p collisions at the same center of mass energies to understand the thermodynamic properties of the produced systems.Table 1
χ 2 / n d f and different extracted parameters from Tsallis distribution fit (see (7)) to the pT spectra of (a) π+, K+, and p for most peripheral Au + Au collisions and p+p collisions at sNN = 200 GeV; (b) π++π-/2, K++K-/2, and p+p-/2 for most peripheral Pb + Pb collisions and π++π-/2, K++K-/2 and p+p-/2 for p+p collisions at sNN = 2.76 TeV.
(a) Au + Au 200 GeV (peripheral) Particle χ 2 / n d f V (fm)3 T (GeV) q π + 0.178 2242.35 ± 18.23 0.093 ± 0.003 1.113 ± 0.004 K + 0.842 223.358 ± 0.010 0.085 ± 0.011 1.120 ± 0.002 p 0.157 45.62 ± 2.45 0.205 ± 0.012 1.025 ± 0.029  p + p 200 GeV  π+ 0.073 5293.1 ± 23.3 0.088 ± 0.004 1.103 ± 0.005 K+ 0.297 2443.7 ± 33.7 0.066 ± 0.002 1.161 ± 0.009 p 1.055 360.5 ± 9.6 0.113 ± 0.004 1.070 ± 0.000  (b) Pb + Pb 2.76 TeV (peripheral) Particle χ 2 / n d f V (fm)3 T (GeV) q π++π-2 1.626 6134.30 ± 12.61 0.082 ± 0.001 1.156 ± 0.002  K++K-2 0.211 629.69 ± 14.80 0.133 ± 0.003 1.120 ± 0.004  p+p-2 0.193 39.623 ± 1.052 0.238 ± 0.012 1.052 ± 0.009   p + p 2.76 TeV  π++π-2 0.985 1102.46 ± 6.82 0.070 ± 0.001 1.163 ± 0.001  K++K-2 0.094 127.13 ± 1.33 0.096 ± 0.005 1.144 ± 0.005  p+p-2 0.006 291.39 ± 10.62 0.069 ± 0.003 1.133 ± 0.002Figure 3
Mass dependence of the thermodynamic parameters using (7) for most peripheral Au + Au and p+p collisions at sNN = 200 GeV.Figure 4
Mass dependence of thermodynamic parameters using (7) for most peripheral Pb + Pb and p+p collisions at sNN = 2.76 TeV.It is found that the volume parameter decreases with increase in particle mass. It indicates that different particles have different freeze-out surfaces depending on their mass. This is observed for bothp+p and A+A collisions and thus indicates a mass dependent differential freeze-out scenario in high-energy collisions. In view of the above freeze-out picture, taking T as the freeze-out temperature, we observe from Figures 3 and 4 that the freeze-out temperature slightly increases with particle mass. This observation goes in line with the earlier findings of [51], which uses a slightly different approach. It clearly indicates that the heavier particles have higher freeze-out temperature as compared to lighter particles. This observation goes in line with an intuitive expectation of mass dependent particle freeze-out or differential freeze-out. In literature, one finds various freeze-out scenarios like single freeze-out [26], strange and nonstrange particles having two different freeze-outs [52, 53], and so forth. However, our findings from the analysis of high-energy particle spectra with nonextensive statistics reveal a consistent differential freeze-out scenario with particles and their counter parts, the antiparticles freezing out at the same time (characterized by freeze-out temperature, freeze-out volume) from the produced fireballs. The nonextensive measure, that is, the q-parameter obtained from our studies, is consistent with its value, 1<q<1.2, for both p+p and A+A collisions. This goes in line with the expectations of the q-values in high-energy collisions [48].Also, Figures3 and 4 indicate a similar mass ordering for p+p systems, which indicates that the systems formed in most peripheral A+A collisions and p+p collisions are of similar thermodynamic nature. However, it should be noted that the obtained parameters seem to be unphysical; that is, volume or the radius in p+p collisions seems to be higher than the peripheral A+A collisions at the same energy. This problem is observed to be circumvented, when one adds the radial flow to the distribution. This is seen in the subsequent discussions and in Table 3.Secondly, we perform the fitting of invariantpT spectra of identified particles using (6) which contains an additional parameter as chemical potential (μ). The same fitting procedure is used to fit Au + Au most peripheral collisions at sNN = 200 GeV for identified particles as well as their antiparticles. The Tsallis form of invariant yields including μ is used in Figure 15. It fits up to pT~3 GeV/c for most peripheral collisions. The goodness of fit can be seen from χ2/ndf values shown in Figure 5.Figure 5
χ 2 / ndf values are shown for most peripheral Au + Au collisions at sNN = 200 GeV using (6) as a fit function.The mass dependence of extracted parameters is shown in Figure6. Also, the values of the parameters are tabulated in Table 2. In Figure 6, it is shown that all the parameters follow the same mass dependent trend as is discussed in the previous cases. Besides, an additional parameter, chemical potential (μ), is shown.Table 2
χ 2 / n d f and different extracted parameters from Tsallis distribution fit (see (6)) to the pT spectra of π±, K±, p, and p- for most peripheral Au + Au collisions at sNN = 200 GeV.
Au + Au 200 GeV (peripheral) Particle χ 2 / ndf V (fm)3 T (GeV) q µ (GeV) π+ 0.191 2125.93± 56.06 0.094± 0.004 1.113± 0.004 0.0050± 0.0003 π- 0.325 1895.78± 275.34 0.097± 0.004 1.110± 0.004 0.0080± 0.0003 K+ 1.415 63.64± 12.60 0.089± 0.021 1.128± 0.025 0.086± 0.009 K- 1.309 113.82± 68.83 0.090± 0.012 1.126± 0.023 0.085± 0.003 p 0.170 10.01± 3.71 0.213± 0.029 1.026± 0.024 0.310± 0.190 p- 0.111 11.62± 3.48 0.183± 0.015 1.051± 0.014 0.437± 0.103Table 3
χ 2 / n d fand different extracted parameters from Tsallis distribution fit (see (8)) to thepTspectra of (a)π+,K+, andpfor most peripheral Au + Au collisions andp+pcollisions at sNN = 200 GeV; (b)π++π-/2, K++K-/2, andp+p-/2for most peripheral Pb + Pb collisions and π++π-/2, K++K-/2, andp+p-/2 forp+p collisions at sNN = 2.76 TeV.
(a) Au + Au 200 GeV (peripheral) Particle χ 2 / n d f V (fm)3 T (GeV) β q π + 16.558 302.680 ± 7.591 0.152±0.001 0.475±0.013 1.021±0.004 K + 2.2461 30.005±1.676 0.154±0.044 0.363±0.008 1.020±0.021 p 0.1825 17.731±1.952 0.247±0.000 0.147±0.018 1.004±0.041  p + p200 GeV  π + 1.370 60.260±0.506 0.123±0.015 0.455±0.015 1.025±0.015 K + 0.996 7.914±0.878 0.160±0.004 0.400±0.013 1.027±0.038 p 3.144 25.970±0.667 0.194±0.017 0.258±0.019 1.010±0.010  (b) Pb + Pb 2.76 TeV (peripheral) Particle χ 2 / n d f V (fm)3 T (GeV) β q π + + π - 2 5.657 444.714±7.931 0.185±0.008 0.480±0.028 1.027±0.006  K + + K - 2 1.025 149.038±11.465 0.200±0.133 0.460±0.021 1.023±0.010  p + p - 2 0.258 14.507±1.107 0.299±0.090 0.260±0.055 1.013±0.024   p + p2.76 TeV  π + + π - 2 18.290 191.039 ± 15.890 0.116 ± 0.003 0.601 ± 0.013 1.038±0.004  K + + K - 2 8.006 10.164±1.097 0.197±0.005 0.433 ± 0.016 1.023±0.005  p + p - 2 7.596 1.839±0.256 0.224 ± 0.006 0.310 ± 0.014 1.014±0.004Figure 6
Mass dependence of extracted thermodynamic parameters using (6) for different particles and their corresponding antiparticles for most peripheral Au + Au collisions at sNN = 200 GeV.It is to be mentioned here that the chemical potential is extracted from the invariant yields at the thermal freeze-out, where all the elastic processes cease. Hence the obtained chemical potential here is called thermal chemical potential. At this stage one can introduce chemical potentials for all particle species (μi). Following [54], the relationship between chemical potential at kinetic freeze-out and chemical potential at chemical freeze-out is given by(9)μth=μchTthTch+m1-TthTch,where (Tch,Tth) and (μch,μth) are temperature and chemical potential at chemical and thermal freeze-out, respectively. Equation (9) clearly shows a mass dependency of thermal chemical potential.The extracted parameters shown in Figure6 show a consistent picture for particles and antiparticles in case of most peripheral collisions. It is evident from Figure 6 that the μ-parameter increases with particle mass. It shows that particles and antiparticles have same freeze-out conditions. It is worth mentioning here that the parameters T and μ extracted from identified particle spectra using (6) are related by [17](10)T=T0+q-1μ,where T=T0 for μ=0. Our extracted parameters using (6) (i.e., T,μ, and q) and (7) (i.e., T≡T0) seem to obey the above relation.Finally, the invariantpT spectra of identified particles are fitted with (8) for most peripheral collisions of Au + Au at sNN = 200 GeV and Pb + Pb at sNN = 2.76 TeV. Besides, the fitting of identified particles for p+p collisions is done for the same energies. The fitting results for Au + Au collisions at sNN = 200 GeV and the results for Pb + Pb collisions at sNN = 2.76 TeV are shown in Figures 16 and 17, respectively. Also, the fitting results for p+p collisions at sNN = 200 GeV and sNN = 2.76 TeV are shown in Figures 18 and 19, respectively.The fitting is performed up topT~3 GeV/c in A+A collisions and pT~2.5 GeV/c in p+p collisions; only a slight deviation at higher pT is observed in pion spectra. The values of various parameters obtained from the fits are provided in Table 3. It is observed from the values of χ2/ndf (Table 3) that the performance of the present approach is quite satisfactory, except in case of pion, which can be seen from Figures 7 and 8.Figure 7
χ 2 / ndf values are shown for most peripheral Au + Au collision and p+p collisions at sNN = 200 GeV using (8) as a fit function.Figure 8
χ 2 / n d f values are shown for most peripheral Pb + Pb collisions and p+p collisions at sNN = 2.76 TeV using (8) as a fit function.Then we make a mass dependent study of these parameters. Figures9 and 10 represent the behaviour of the extracted parameters. The parameters like volume, temperature, and q show a similar mass dependent behaviour as that of the parameters obtained from the fitting using (6) and (7). Here, through (8), we have introduced an additional radial flow parameter (β). Figures 9 and 10 show the decrease of radial flow with increase of particle mass, which is a signature of hydrodynamic evolution of the system created in high-energy collisions. The change of temperature parameter with the use of different forms of Tsallis distribution could be observed in [55]. In our analysis we have also observed that not only the temperature but also the other parameters are changing with the use of different forms of Tsallis distribution. For example, it is observed that, with the inclusion of radial flow, the value of the q-parameter becomes less, as compared to that observed in the other forms of Tsallis distribution, while describing the pT-spectra. This may hint that the nonextensivity is shared by the dynamics of the system [56–58].Figure 9
Mass dependence of extracted thermodynamic parameters of Tsallis distribution including radial flow using (8) for most peripheral Au + Au collisions and p+p collisions at sNN = 200 GeV.Figure 10
Mass dependence of extracted thermodynamic parameters of Tsallis distribution including radial flow using (8) for most peripheral Pb + Pb collisions and p+p collisions at sNN = 2.76 TeV.Figure 11
Fitting of invariantpT spectra with Tsallis distribution using (7) for π+ [18], K+ [18], and p [18] in most peripheral Au + Au collisions at sNN = 200 GeV. Symbols represent the data points and line is the fit function.Figure 12
Fitting of invariantpT spectra with Tsallis distribution using (7) for π++π-/2 [19, 20], K++K-/2 [19, 20], and p+p¯/2 [19, 20] in most peripheral Pb + Pb collisions at sNN = 2.76 TeV. Symbols represent the data points and line is the fit function.Figure 13
Fitting of invariantpT spectra with Tsallis distribution using (7) for π+ [21], K+ [21], and p [21] in p+p collisions at sNN = 200 GeV. Symbols represent the data points and line is the fit function.Figure 14
Fitting of invariantpT spectra with Tsallis distribution using (7) for π++π-/2 [19, 20], K++K-/2 [19, 20], and p+p¯/2 [19, 20] in p+p collisions at sNN = 2.76 TeV. Symbols represent the data points and line is the fit function.Figure 15
Fitting of invariantpT spectra with Tsallis distribution including chemical potential using (6) for π+ [18], K+ [18], and p [18] in most peripheral Au + Au collisions at sNN = 200 GeV. Symbols represent the data points and the line is the fit function.Figure 16
Fitting of invariantpT spectra with Tsallis distribution including radial flow using (8) for π+ [18], K+ [18], and p [18] in most peripheral Au + Au collisions at sNN = 200 GeV.Figure 17
Fitting of invariantpT spectra with Tsallis distribution including radial flow using (8) for π++π-/2 [19, 20], K++K-/2 [19, 20], and p+p¯/2 [19, 20] in most peripheral Pb + Pb collisions at sNN = 2.76 TeV; symbols represent the data points and the line is the fit function.Figure 18
Fitting of invariantpT spectra with Tsallis distribution including radial flow using (8) for π+ [21], K+ [21], and p [21] in p+p collisions at sNN = 200 GeV. Symbols represent the data points and the line is the fit function.Figure 19
Fitting of invariantpT spectra with Tsallis distribution including radial flow using (8) for π++π-/2 [19, 20], K++K-/2 [19, 20], and p+p¯/2 [19, 20] in p+p collisions at sNN = 2.76 TeV. Symbols represent the data points and the line is the fit function.


## 4. Conclusion

We present a systematic study of the parameters extracted from transverse momentum spectra of produced identified particles in Au + Au collisions atsNN = 200 GeV and Pb + Pb collisions at sNN = 2.76 TeV. Furthermore, the same analysis is done for p+p collisions at both, RHIC and LHC energies. In the present work different forms of invariant yields using Tsallis distributions are used. Parameters obtained from three different forms in most peripheral Au + Au collisions at sNN = 200 GeV and in most peripheral Pb + Pb collisions at sNN = 2.76 TeV follow a similar trend. Similar behaviour has been observed in p+p collisions for both the energies. The chemical potential plays an important role in particle production mechanism. We observe that the particles produced at early stage have higher thermal chemical potential as compared to particles produced at the later stage. The radial flow extracted from the particle spectra of peripheral A+A collisions is similar to p+p collisions at the same center of mass energy. In the discussed low-pT region, the radial flow decreases with increasing mass, which goes in line with hydrodynamic description of the evolution of a fireball created in high-energy collisions. Therefore, the present analysis gives a systematic information of hydrodynamical and thermodynamical evolution of the system, while making a direct comparison with hadronic and peripheral nuclear collisions. We observe a mass dependent differential freeze-out scenario in high-energy hadronic and nuclear collisions, which is also supported by similar observations by other authors [59].



---

*Source: 4149352-2016-11-15.xml*