# Dependence of Temperatures and Kinetic Freeze-Out Volume on Centrality in Au-Au and Pb-Pb Collisions at High Energy

Authors: Muhammad Waqas; Fu-Hu Liu; Zafar Wazir
Journal: Advances in High Energy Physics
 (2020)

---

## Abstract

Centrality-dependent double-differential transverse momentum spectra of negatively charged particles (π−, K−, and p¯) at the mid(pseudo)rapidity interval in nuclear collisions are analyzed by the standard distribution in terms of multicomponent. The experimental data measured in gold-gold (Au-Au) collisions by the PHENIX Collaboration at the Relativistic Heavy Ion Collider (RHIC) and in lead-lead (Pb-Pb) collisions by the ALICE Collaboration at the Large Hadron Collider (LHC) are studied. The effective temperature, initial temperature, kinetic freeze-out temperature, transverse flow velocity, and kinetic freeze-out volume are extracted from the fitting to transverse momentum spectra. We observed that the mentioned five quantities increase with the increase of event centrality due to the fact that the average transverse momentum increases with the increase of event centrality. This renders that larger momentum (energy) transfer and further multiple scattering had happened in central centrality.


---

## Body

## 1. Introduction

One of the most important questions in high-energy collisions is the identification of various phases of dense matter. It is expected to reach a deconfined state of matter (quarks and gluons) at high energy or density. This state of matter is called Quark-Gluon Plasma (QGP), which was obtained in the early universe shortly after the big bang prior to the condensation in hadrons. The characterization of phase transition in finite systems is a fascinating multidisciplinary topic which has been studied for decades [1, 2] within different phenomenological applications. The Relativistic Heavy Ion Collider (RHIC) and Large Hadron Collider (LHC) have been providing excellent tools to determine the phase structure of the strongly interacting Quantum Chromodynamics (QCD) matter [3–5] and to study the properties of QGP [6–8].Within the framework of statistical thermal models, it is assumed that the initial stage of collisions of nuclei at the RHIC and LHC [9–11] gives a tremendous amount of temperature, where a hot and dense “fireball” over an extended region for a very short period of time (almost 10−22 seconds) is formed. The fireball consists of QGP, and it instantly cools which results in the expansion of the system (the change of the volume or density of the system), and the partons recombine into a blizzard of hadronic matter [12]. After the hadronization of the fireball, the hadrons continuously interact with each other and the particle number changes. This process results in a decrease of temperature and at a certain value where the reaction process stops and the temperature at this point is called the “chemical freeze-out temperature” (Tch). At the stage of chemical freeze-out, the yield ratios of different types of particles remain invariant [13].However, the rescattering process still takes place which continues to build up the collective (hydrodynamical) expansion. Resultantly, the matter becomes dilute and the mean free path of the given hadrons for the elastic reaction processes becomes comparable with the size of the system. At this stage, the rescattering process stops, which results in the decoupling of hadrons from the rest of the system [14]. This stage is called as the kinetic or thermal freeze-out stage, and the temperature at this stage is called the kinetic or thermal freeze-out temperature (T0). After this stage, the particle’s energy/momentum spectrum is frozen in time and it is the least stage of the system evolution. Meanwhile, at this stage, the particle’s movement is also affected by the flow effect which should be excluded when one extracts T0. To describe the flow effect, one may use the transverse flow velocity βT.The above discussedT0 and βT can be extracted from transverse momentum (pT) spectra of particles. Also, from pT spectra, one can extract the initial temperature Ti according to the color string percolation model [15–17]. Generally, if the flow effect is not excluded in the temperature parameter, this type of temperature is called the effective temperature (T). At least, three types of temperatures, T, Ti, and T0, can be extracted from pT spectra. Although the yield ratios of different types of particles can be also obtained from the normalizations of pT spectra and then Tch can be also extracted from pT spectra, we mainly extract other three types of temperatures and βT from pT spectra in this paper due to their more pending situations.In addition, volume is also an important parameter in high-energy collisions. The volume occupied by the ejectiles when the mutual interactions become negligible, and the only force they feel is the columbic repulsive force, is known as the kinetic freeze-out volume (V). There are various freeze-out volumes at various freeze-out stages, but we are only focusing on the kinetic freeze-out volume V in the present work. As we know,Vgives the information of the coexistence of phase transition and is important in the extraction of multiplicity, microcanonical heat capacity, and its negative branch or shape of the caloric curves under the thermal constraints [18–22]. In this paper, the fifth quantity extracted from pT spectra is V. By way of parenthesis, the mean pT, i.e., pT, is also obtained.The study of three types of temperatures, transverse flow velocity, and kinetic freeze-out volume is a very wide, interesting, and of course huge project. However, in this paper, we will only analyze the centrality dependences of the five quantities in gold-gold (Au-Au) collisions at 200 GeV and in lead-lead (Pb-Pb) collisions at 2.76 TeV. Only thepT spectra of negatively charged pions (π−), kaons (K−), and antiprotons (p¯) are used in the extraction. These representational spectra are enough to extract the necessary centrality dependences.The remainder of this paper orderly consists of the method and formalism, results and discussion, and summary and conclusions which are presented in Sections2, 3, and 4, respectively.


## 2. The Method and Formalism

Soft excitation and hard scattering processes are the two generally considered processes for the particle production. The soft excitation process contributes in a narrowpT range which is less than 2–3GeV/c or up to 4–5GeV/c and is responsible for the production of most of the light flavored particles. The soft excitation process has various choices of formalisms including but are not limited to the Hagedorn thermal model (Statistical-Bootstrap model) [23], the standard distribution [24], the blast-wave model with Boltzmann-Gibbs statistics [25–27], the blast-wave model with Tsallis statistics [28–30], and the current thermodynamical-related models [31–34]. The main contributor to the produced particles is the soft excitation process.If necessary, for the hard excitation process, there is limited choice of formalisms [35–37] and it can be described by the theory of strong interaction. In fact, the contribution of the hard scattering process is parameterized to an inverse power law, i.e., the Hagedorn function [23]:
(1)fHpT=ApT1+pTp0−n,where p0 and n are free parameters and A is the normalized constant related to the free parameters. The inverse power law is obtained from the calculus of QCD [3–5] and has at least three revisions, which is out of focus of the present work and will not be discussed further.Different probability density functions can be used to describe the contributions of the soft excitation and hard scattering processes. Due to few fractions and being earlier than the kinetic freeze-out stage, the hard scattering process does not contribute largely toT0 and βT in general. In fact, the contribution of the hard scattering process can be neglected if we study the pT spectra in a narrow range, say pT<2–3GeV/c or extending to pT<4–5GeV/c, for which only the contribution of the soft excitation process is indeed needed. In our opinion, various distributions show similar behaviors in case of fitting the data with acceptable representations, which results pT in similar pTpT2 with different parameters.For the spectra contributed by the soft excitation process, we can choose the standard distribution, as it is very close in concept to the ideal gas model. The standard distribution is the combination of the Boltzmann, Fermi-Dirac, and Bose-Einstein distributions. The probability density function of the standard distribution in terms ofpT at midrapidity is generally as follows [24]:
(2)fSpT=CgV′2π2pTpT2+m02×exppT2+m02T+S−1,where C is the normalization constant, V′ is the fitted kinetic freeze-out volume (in terms of interaction volume) of the emission source at the stage of kinetic freeze-out as discussed above, g is the degeneracy factor for pion and kaon (or proton) and has the value of 3 (or 2), m0 is the rest mass of the considered particle, S=−1+1 is for bosons (fermions), and T is the effective temperature as discussed above.By considering the experimental rapidity rangeymin,ymax around the midrapidity, the accurate form of Equation (2) is as follows [24]:
(3)fSpT=CgV′2π2pTpT2+m02∫yminymaxcoshy×exppT2+m02coshyT+S−1dy.T and V′ in Equations (2) and (3) are free parameters. In most cases, the single-component standard distribution is not enough for the simultaneous description of the low-pT region. So we have to use a two-component standard distribution.In particular, in some cases, the multicomponent (l-component) standard distribution has to be used, which can be demonstrated as follows:
(4)fSpT=∑i=1lkiCigVi′2π2pTpT2+m02∫yminymaxcoshy×exppT2+m02coshyTi+S−1dy,where Ci is the normalization constant, ki represents the fraction contributed by the ith component, and Ti and Vi′ are free parameters denoted the kinetic freeze-out temperature and volume, respectively, corresponding to the ith component.In case of considering both the contributions of the soft excitation and hard scattering processes, we use the superposition in principle:(5)f0pT=kSfSpT+1−kSfHpT,where kS is the contribution ratio of the soft excitation process. The contribution ranges of the soft excitation and hard scattering processes described by Equation (5) are intersectant in the low-pT region.Another type of superposition which uses the usual step functionθx based on the Hagedorn model [23] is
(6)f0pT=A1θp1−pTfSpT+A2θpT−p1fHpT, where A1 and A2 are the normalized constants which synthesize A1fSp1=A2fHp1. The contribution ranges of the soft excitation and hard scattering processes described by Equation (6) are segregative at pT=p1.In the present work, we will study only the first component in Equation (5) or (6), because we do not study a very wide pT range. Meanwhile, we use the two-component standard distribution, i.e., l=2, in Equation (4) in which the first or second component has no particular priority. As probability density functions, Equations (1)–(6) are normalized to 1 naturally. When we compare the probability density functions with the experimental data which appear usually in other forms, a suitable transformation for the probability density function is certainly needed. Considering the treatment of normalization, the real fitted kinetic freeze-out volume should be Vi′/Ciki which will be simply used in the following section as the symbols, Vi or V=∑i=1lVi.


## 3. Results and Discussion

Figure1 presents the event centrality-dependent double-differential pT spectra, 1/2πpTd2N/dpTdy, of the identified particles (π−, K−, and p¯) produced in Au-Au collisions at sNN=200GeV in the midpseudorapidity interval of ∣η∣<0.35, where N and y stands for the number of particles and rapidity, respectively. The symbols are the experimental data measured by the PHENIX Collaboration [38]. The spectra of centralities 0–5%, 5–10%, 10–15%, 15–20%, 20–30%, 30–40%, 40–50%, 50–60%, 60–70%, 70–80%, and 80–92% in the three panels are multiplied by 20, 10, 5, 2.5, 1.5, 1, 1, 1, 1, 1, and 1, respectively. In addition, the spectra of centralities 15–20%, 20–30%, 30–40%, 40–50%, 50–60%, 60–70%, 70–80%, and 80–92% in panel (a) are remultiplied by 0.8, 0.7, 0.6, 0.45, 0.3, 0.24, 0.17, and 0.15, respectively. The curves are our fitting results by using the two-component standard distribution, Equation (4) with l=2. The substantially togethered event centralities, the values of free parameters (T1, T2, V1, and V2), χ2, and number of degree of freedom (ndof) are listed in Table 1. One can see the well-approximate description of the model results to the experimental data of the PHENIX Collaboration in special pT ranges in high-energy Au-Au collisions at the RHIC.Figure 1
Centrality-dependent1/2πpTd2N/dpTdy of (a) π−, (b) K−, and (c) p¯ produced in ∣η∣<0.35 in Au-Au collisions at 200 GeV. The symbols represent experimental data of the PHENIX Collaboration [38], while the curves are the results of our fits by using the two-component standard distribution, Equation (4) with l=2. The spectra of centralities 0–5%, 5–10%, 10–15%, 15–20%, 20–30%, 30–40%, 40–50%, 50–60%, 60–70%, 70–80%, and 80–92% in the three panels are multiplied by 20, 10, 5, 2.5, 1.5, 1, 1, 1, 1, 1, and 1, respectively. In addition, the spectra of centralities 15–20%, 20–30%, 30–40%, 40–50%, 50–60%, 60–70%, 70–80%, and 80–92% in panel (a) are remultiplied by 0.8, 0.7, 0.6, 0.45, 0.3, 0.24, 0.17, and 0.15, respectively.
(a)(b)(c)Table 1
Values of effective temperatures (T1 and T2), volumes (V1 and V2), χ2, and number of degree of freedom (ndof) corresponding to the curves in Figures 1 and 2.
FigureParticleCentralityT1 (GeV)T2 (GeV)V1 ((fm)3)V2 ((fm)3)χ2ndofFigure1Au-Au200 GeVπ−0–5%0.181±0.0140.268±0.006745±1001967±30945245–10%0.138±0.0130.270±0.006564±802054±174142410–15%0.179±0.0150.250±0.005430±622004±1532332415–20%0.167±0.0120.247±0.003354±611862±320632420–30%0.151±0.0090.245±0.004360±391854±216872430–40%0.125±0.0150.237±0.006195±302005±1892762440–50%0.148±0.0140.227±0.004106±251585±1713882450–60%0.098±0.0110.227±0.00359±111048±1531462460–70%0.125±0.0100.274±0.00637±11984±109712470–80%0.112±0.0140.209±0.00718±3296±372562480–92%0.068±0.0150.221±0.0068.0±1.0213±136924K−0–5%0.201±0.0090.271±0.008196±361081±17817125–10%0.188±0.0070.255±0.00832±61116±71181210–15%0.211±0.0130.250±0.00637±61030±102241215–20%0.135±0.0140.248±0.00835±5836±105141220–30%0.219±0.0150.241±0.00610±1733±66181230–40%0.156±0.0100.238±0.0074.0±0.3665±82131240–50%0.130±0.0110.235±0.00813±1.2253±39241250–60%0.170±0.0130.230±0.0061.6±0.2577±4891260–70%0.163±0.0150.225±0.0060.50±0.02544±50351270–80%0.155±0.0150.220±0.0070.30±0.01505±602151280–92%0.130±0.0100.218±0.0050.030±0.001460±1306612p¯0–5%0.280±0.0170.320±0.00336±7786±12017185–10%0.272±0.0150.317±0.00231±7775±70201810–15%0.265±0.0180.313±0.00616±3747±861621815–20%0.185±0.0160.313±0.00523±4710±80351820–30%0.211±0.0130.309±0.00613±1678±6541411830–40%0.247±0.0180.317±0.00710±1634±50111840–50%0.244±0.0150.308±0.0084.0±0.5597±37391850–60%0.237±0.0130.300±0.0041.6±0.3574±70431860–70%0.210±0.0130.285±0.0060.50±0.02476±45281870–80%0.170±0.0140.267±0.0050.30±0.01455±501101880–92%0.172±0.0100.218±0.0040.40±0.01406±40518Figure2Pb-Pb2.76 TeVπ−0–5%0.179±0.0100.340±0.0042517±2084099±252206375–10%0.172±0.0120.331±0.0042226±1353468±2011013710–20%0.165±0.0100.312±0.0051395±702813±176703720–30%0.161±0.0140.300±0.0061012±1032427±231943730–40%0.147±0.0130.280±0.005534±563502±1652443740–50%0.129±0.0120.289±0.005307±383600±1981103750–60%0.108±0.0100.281±0.005176±323380±2132063760–70%0.114±0.0150:290±0:00594±133150±1471253570–80%0.130±0.0100.287±0.00745±92930±2065023580–90%0.118±0.0140.287±0.00616±21945±10211935K−0–5%0.277±0.0100.361±0.007170±363442±2083325–10%0.270±0.0090.359±0.009152±103237±240153210–20%0.252±0.0090.354±0.00797±113190±18783220–30%0.265±0.0090.350±0.00653±73136±200223230–40%0.247±0.0080.395±0.00627±43014±198123240–50%0.210±0.0080.341±0.00622±42828±1901193250–60%0.213±0.0110.345±0.0058.1±1.02659±210263260–70%0.226±0.0090.335±0.0055.8±1.32132±176463270–80%0.225±0.0100.332±0.0062.6±0.31736±134713280–90%0.170±0.0090.320±0.0060.50±0.051150±1207532p¯0–5%0.426±0.0120.421±0.00755±61965±19549335–10%0.300±0.0100.425±0.00575±61782±1643733310–20%0.427±0.0100.397±0.35250±51748±140433320–30%0.405±0.0130.395±0.00631±41690±130523330–40%0.400±0.0120.382±0.00529±41534±158203340–50%0.352±0.0110.383±0.00615±2401±1351093350–60%0.331±0.0120.380±0.0056.0±0.51357±125673360–70%0.310±0.0120.369±0.0068.7±1.01069±90753370–80%0.321±0.0150.333±0.0041.0±0.3978±501043380–90%0.310±0.0140.328±0.0060.70±0.04708±6012233Figure2 is the same as Figure 1, but it gives the results for π−, K−, and p¯ in different centrality bins in Pb-Pb collisions at 2.76 TeV in the midrapidity interval ∣y∣<0.5. The experimental data of the ALICE Collaboration is represented by the symbols [39], whereas the spectra of Pb-Pb are scaled by the factor of 2n and n changes from 9 to 0 with the change of event centrality such as from 0–5% to 80–90%. The related parameters and the existing centralities are listed together in Table 1. One can see the well-approximate description of the model results to the experimental data of the ALICE Collaboration in special pT ranges in high-energy Pb-Pb collisions at the LHC.Figure 2
Same as Figure1 but for the spectra of π−, K−, and p¯ in y<0.5 in Pb-Pb collisions at sNN=2.76TeV. The symbols represent the measured data of the ALICE Collaboration [39], where the spectra are scaled by factors 2n and n changes from 9 to 0 as the event centrality changes from 0–5% to 80–90%.
(a)(b)(c)To study the change trend of parameters with centrality, Figure3 shows the dependence of effective temperature on centrality for the productions of π−, K−, and p¯ in different centrality bins in the Au-Au and Pb-Pb collisions at 200 GeV and 2.76 TeV, respectively. Panel (a) shows the result for Au-Au collisions, while panel (b) shows the result for Pb-Pb collisions. Different symbols represent different particles. One can see the clear decrease of effective temperature from the central to peripheral collisions. The reason behind it is as follows: the more violent collisions in central collisions where it can get a degree of higher excitation and also involve more number of participants in interactions, while they decrease from the central to peripheral collisions. The effective temperatures in collisions at the LHC in different centrality bins are higher than those at the RHIC due to more energy deposition in collisions at the LHC.Figure 3
Dependence of effective temperatureT in different centrality bins C in (a) Au-Au collisions at 200 GeV and (b) Pb-Pb collisions at 2.76 TeV. The obtained values corresponded to the identified particles extracted from the experimental pT spectra.
(a)(b)Figure4 is the same as Figure 3; however, it shows the result for the dependence of kinetic freeze-out volume V in events with different centralities C, where V=V1+V2 due to the additivity of volume. One can see that the kinetic freeze-out volume decreases from the central to peripheral collisions, as the number of participant nucleons decreases from the central to peripheral collisions depending on the interaction volume. Due to the large number of binary collisions by the rescattering of partons, the system with more participants reaches quickly to an equilibrium state, but because of the decreases in centrality and the decreases in the number of participants, the system goes slowly to an equilibrium state. The large volume and more number of participants in the central collisions may indicate the occupation of superdense hadronic matter, but of course, further and more complete information about the local energy density of superhadronic matter is needed to study the possible phase transition of QGP. The figure shows a volume differential scenario too. The heavier the particle is, the smaller the kinetic freeze-out volume has, which shows the early freeze-out of heavier particles as compared to the lighter particles and suggests different freeze-out surfaces for different particles. Such result can be found in the literature [40, 41].Figure 4
Same as Figure3, but for dependence of kinetic freeze-out volume V on centrality C.
(a)(b)The dependence of the mean transverse momentumpT in different centrality events is shown in Figure 5. The symbols represent pT for different particles obtained from the fitting function Equation (4) with l=2 over a pT range from 0 to 5GeV/c, where the parameter values are listed in Table 1. One can see that pT decreases from the central to peripheral Au-Au and Pb-Pb collisions for all particle species, and it may be caused due to decreasing the participant nucleons from the central to peripheral collisions and this result is similar to Ref. [42]. It is also important to notice that p¯ spectra exhibit a concave shape in the peripheral events, which is well described by the power-law parametrization as observed in Ref. [43], but this curvature decreases with the increasing centrality and it leads to an almost exponential dependence on pT for the most central collisions. Furthermore, pT for heavier particles is larger than that for lighter ones, and pT at LHC energy is slightly larger than that at RHIC energy. The increase of pT in central collisions and with the massive mass of the particle may indicate the collective radial flow, and the same behavior is observed at a few GeV [44] and more than 10 GeV [45].Figure 5
Same as Figure3, but for dependence of the mean transverse momentum pT on centrality C.
(a)(b)Figure6 is the same as Figure 5; however, it demonstrates the result for the initial temperature Ti, where Ti is obtained by the root-mean-square pT divided by 2, i.e., pT2/2, according to the color string percolation model [15–17]. The symbols are the representation of the results obtained from the fitting function Equation (4) with l=2 over a pT range in 0–5GeV/c and with the parameter values listed in Table 1. The mass differential temperature scenario is also observed. It is necessary to point out that the initial temperature obtained in this work is larger than the effective temperature which is in agreement with the order of time evolution of the interacting system.Figure 6
Same as Figure3, but for dependence of the initial temperature Ti on centrality C.
(a)(b)Figure7 shows the dependence of (a, b) effective temperature T on the rest mass m0 and (c, d) mean transverse momentum pT on the mean energy or mean moving mass m¯ in the rest frame of the emission source for (a, c) Au-Au collisions at 200 GeV and (b, d) Pb-Pb collisions at 2.76 TeV in different centrality classes. The symbols represent the derived quantities according to the free parameters listed in Table 1, where T=T1V1+T2V2/V1+V2, pT=∫pTfSpTdpT, and m¯=pT2+pz2+m02, where pz is the longitudinal momentum and its distribution can be obtained according to pT distribution if we assume isotropic emission in the source rest frame [46]. The lines are the results fitted for the values of derived quantities if we assume linear correlations are existent.Figure 7
Dependence of (a, b) effective temperatureT on the rest mass m0 and (c, d) mean transverse momentum pT on the mean energy or mean moving mass m¯ for (a, c) Au-Au collisions at 200 GeV and (b, d) Pb-Pb collisions at 2.76 TeV in different centrality classes. The lines are the results fitted for the values of derived quantities.
(a)(b)(c)(d)The intercepts in Figures7(a) and 7(b) are regarded as the kinetic freeze-out temperature T0, and the slopes in Figures 7(c) and 7(d) are regarded as the transverse flow velocity βT [47–50]. The dependences of (a) T0 on C and (b) βT on C, as well as the correlation between (c) βT and βT, are presented in Figure 8. One can see that T0 and βT decrease with the increase of C, and T0 increases with the increase of βT, which renders that central collisions display higher excitation and quicker expansion than peripheral collisions due to more energy deposition in central collisions.Figure 8
Dependence of (a) kinetic freeze-out temperatureT0 on centrality C and (b) transverse flow velocity βT on centrality C.
(a)(b)(c)It should be noted thatT0 shown in Figure 8 is larger than 160–170 MeV which is regarded as the chemical freeze-out temperature Tch of phase transition from hadronic matter to QGP. As expected, T0 should be less than or equal to Tch due to time evolution. The difference between our results and Tch is regarded as different “thermometers” (methods) used in the extraction of temperature. In our opinion, a unified “thermometer” (method) should be used in the extraction of temperature. Or, one may find a relation to convert one temperature to another one, as what one did between Celsius temperature and degree Fahrenheit.Meanwhile, there is a positive correlation betweenT0 and βT as shown in Figure 8(c). Some studies show a negative correlation between T0 and βT when one uses the blast-wave model [25–30]. At present, one could not decide which correlation is correct. In our opinion, for a given pT spectrum, T0 and βT are negatively correlative if one uses the blast-wave model which gives a large T0 to correspond to a small βT. However, for a set of pT spectra with varying centralities and energies, the situation is dissimilar. The present work uses an alternative method to extract T0 and βT and obtain a positive correlation.In addition, one can see some fluctuations in Figures3–6 and 8. These fluctuations have no particular physics meaning but reflect the statistical and/or systematical fluctuations in the data itself. Although there are fluctuations in the dependences of parameters on centrality, one can see the general decreasing trend of parameters when decreasing the centrality. The considered parameters have a similar trend due to their consistent meaning on the energy deposition which is reflected in terms of the excitation and expansion degree.Both the mean transverse momentumpT and the initial temperature Ti are obtained from pT spectra. The relation of pT and Ti is certainly positive correlation. Figure 9 shows this correlation. One can see that pT increases with the increase of Ti. This correlation is natural since Ti is defined by pT2/2. One can also see that pT2/2 increases with the increase of collision energy and the size of the system, but the latter has a very small effect, so we can neglect it.Figure 9
Dependence of the mean transverse momentumpT on the initial temperature Ti.From Figures4–9, one can see that the considered quantities in the most peripheral Pb-Pb collisions at 2.76 TeV overlap the most central Au-Au collisions with lower energy of 200 GeV, which may be the indication of the formation of similar systems in the most peripheral collisions at higher energies and in the most central collisions at lower energies, and it may support the hypothesis of the effective energy for the particle production [51–54].In the considered Au-Au collisions at 200 GeV and Pb-Pb collisions at 2.76 TeV, the decreasing trend of temperatures and kinetic freeze-out volume from the central to peripheral collisions renders that more energy deposition and then higher excitation and quicker expansion in central collisions. Some studies [9, 13, 55–57] which use other methods such as the blast-wave model show that the kinetic freeze-out temperature in central collisions is less than that in peripheral collisions, though this opposite result can be explained as longer freeze-out time in central collisions.Indeed, the kinetic freeze-out temperature and transverse flow velocity and other quantities are model dependent. We notice that the current blast-wave model uses a small or almost zero transverse flow velocity in peripheral collisions and obtains a larger kinetic freeze-out temperature in peripheral collisions compared to central collisions. If we use a large transverse flow velocity in peripheral collisions, we can obtain a smaller kinetic freeze-out temperature in peripheral collisions compared to central collisions [58].In addition, the kinetic freeze-out temperatureT0 and transverse flow velocity βT are also transverse momentum range dependent. In our opinion, to obtain the parameters as accurately as possible, we should use the transverse momentum range as accurately as possible. The transverse momentum range should not be too narrow or too wide. A too narrow transverse momentum range will exclude the contributions of some particles which should be included. A too wide transverse momentum range will include the contributions of some particles which should be excluded. In fact, model- and pT range-independent T0 and βT are ideal.Generally, the mean transverse momentumpT and the root-mean-square transverse momentum pT2 are model independent. Obtaining the initial temperature by Ti=pT2/2 is a suitable treatment which is regardless of the model, though it is from the color string percolation model [15–17]. It is expected that T0 and βT are related to pT, which results in model-independent T0 and βT.As what we did in our recent work [59], let T0≡kpT/2 and βT≡1−kpT/2m¯, where k is a parameter which can be approximately taken to be 0.3−0.01lnsNN (sNN is in units of GeV) [60], 1/2 is used due to both contributions of the projectile and target participants, and m¯ denotes the mean energy of the considered particles in the source rest frame. If the pT range is wide enough, T0 and βT are also pT range independent.If we useT0 by the above new definition instead of the intercept in the linear relation between T and m0, the mean T0 (∼0.10–0.12 GeV) obtained from Figure 9 is obviously less than that in the present work which is too large compared to others. Meanwhile, if we use βT by the above new definition instead of the slope in the linear relation between pT and m¯, the mean βT (∼0.23–0.28 c) obtained from Figure 9 is less than that in the present work. Regardless of size, the new definitions of T0 and βT are model independent.Before summary and conclusions, we would like to point out that this paper fits only the transverse momentum spectra measured from collisions with varying centralities by the two-component standard distribution. Some centrality dependences of related parameters are found. In our recent work [60–63], we have fitted the spectra measured from collisions with varying energies by the (two-component) standard distribution and/or Tsallis statistics. Some spectra are from small system size, and others are from large system size. The related parameters are found to depend also on energy and the larger nucleus in the projectile and target nuclei.In particular, with the increasing energy, the kinetic freeze-out temperature increases quickly from a few GeV to around 10 GeV and then slowly or slightly from around 10 GeV to more than 10 TeV. This implies that around 10 GeV is a special energy at which the interaction mechanism had changed. In fact, the collision system had undergone a baryon-dominated to meson-dominated final state [64]. This implies that the critical energy of phase transition from hadronic matter to QGP is possibly existent at around 10 GeV.The dependence on the larger nucleus in the projectile and target nuclei is consistent to the dependence on centrality. This implies possibly that there is a critical size from the small to large systems and from the peripheral to central collisions. The data measured by the NA61/SHINE Collaboration [65] show that the nucleon number in the projectile or target nucleus on the onset of deconfinement is ≈10. Meanwhile, the energy on the onset of deconfinement is ≈10 GeV. This double-10 signature is very interesting and should be studied further by various models and methods in the future.We have studied three types of temperatures, namely, the effective temperature, initial temperature, and kinetic freeze-out temperature, in this paper. Although the three types of temperatures are extracted from the transverse momentum spectra, they have different physics meanings. The effective temperature is obtained directly from the fit function, which describes together the degree of the thermal motion and flow effect at the stage of kinetic freeze-out. In the case of excluding the contribution of flow effect from the effective temperature, we expect to obtain the kinetic freeze-out temperature which describes only the thermal motion. The initial temperature in this paper is quoted directly from the color string percolation model [15–17], which is expected to describe the excitation degree of the initial state as what we did in our recent work [66].It is regretful that the chemical freeze-out temperature is not discussed in this paper, though it has wider applications and discussions in the literature [67–69]. The chemical freeze-out temperature describes the excitation degree of the collision system at the stage of chemical freeze-out. Generally, the chemical freeze-out temperature can be obtained from the ratio of particle yields and can be used to map the phase diagram with the chemical potential. In the extensive statistics and/or axiomatic/generic nonextensive statistics [67–69], one may discuss the chemical and/or kinetic freeze-out parameters systematically.


## 4. Conclusions

The main observations and conclusions are summarized as follows:(a)
The transverse momentum spectra ofπ−, K−, and p¯ at mid(pseudo)rapidity produced in different centrality events in Au-Au collisions at 200 GeV and Pb-Pb collisions at 2.76 TeV have been analyzed. The experimental data measured by the PHENIX and ALICE Collaborations are fitted by the two-component standard distribution in which the temperature concept is quite close to the ideal gas model(b)
The effective temperature, initial temperature, kinetic freeze-out temperature, transverse flow velocity, and mean transverse momentum increase with the increase of event centrality from the peripheral to central collisions, which results in a higher excitation degree and quicker expansion velocity in central collisions. The kinetic freeze-out volume increases with the increase of event centrality from the peripheral to central collisions due to more number of participant nucleons taking part in central collisions(c)
The mass-dependent differential effective temperature, initial temperature, kinetic freeze-out volume, and mean transverse momentum are observed. The kinetic freeze-out temperature and transverse flow velocity extracted in this paper do not show a mass-dependent differential scenario due to the reason of methodology. Many quantities are model or method dependent(d)
The formation of similar systems is possible in the most peripheral nucleus-nucleus collisions at high energy and in the most central nucleus-nucleus collisions at low energy. This observation confirms the hypothesis of the effective energy for the particle production. There are many similarities in high-energy collisions



---

*Source: 8198126-2020-05-27.xml*