Analyzing Transverse Momentum Spectra of Pions, Kaons and Protons in p–p, p–A and A–A Collisions via the Blast-Wave Model with Fluctuations

The transverse momentum spectra of different types of particles, π±, K±, p and p¯, produced at mid-(pseudo)rapidity in different centrality lead–lead (Pb–Pb) collisions at 2.76 TeV; proton–lead (p–Pb) collisions at 5.02 TeV; xenon–xenon (Xe–Xe) collisions at 5.44 TeV; and proton–proton (p–p) collisions at 0.9, 2.76, 5.02, 7 and 13 TeV, were analyzed by the blast-wave model with fluctuations. With the experimental data measured by the ALICE and CMS Collaborations at the Large Hadron Collider (LHC), the kinetic freeze-out temperature, transverse flow velocity and proper time were extracted from fitting the transverse momentum spectra. In nucleus–nucleus (A–A) and proton–nucleus (p–A) collisions, the three parameters decrease with the decrease of event centrality from central to peripheral, indicating higher degrees of excitation, quicker expansion velocities and longer evolution times for central collisions. In p–p collisions, the kinetic freeze-out temperature is nearly invariant with the increase of energy, though the transverse flow velocity and proper time increase slightly, in the considered energy range.


Introduction
In high-energy collisions, one of the most important questions is the identification of various phases of dense matter. Quark-gluon plasma (QGP) [1][2][3], which is considered a new state of matter, was produced in the early universe shortly after the Big Bang, prior to the condensation in hadrons. High energy nucleus-nucleus (heavy ion) collisions at the large hadron collider (LHC) [4][5][6][7][8] provide another excellent environment with high temperatures and high density under which QGP are expected to form and to live for a longer lifetime than that at the relativistic heavy ion collider (RHIC) [9]. Presently, the LHC has performed four different types of collisions, proton-proton (p-p), proton-lead (p-Pb), lead-lead (Pb-Pb) and xenon-xenon (Xe-Xe) collisions, at different collision energies. The former two with low multiplicity are not expected to form QGP due to small systems and less energy deposition, though the deconfinement of quarks and gluons may appear. The latter two are expected to form QGP due to large systems and high energy. It is possible that the former two with high multiplicity will also form QGP due to violent collisions.
Kinetic freeze-out temperature, transverse flow velocity and proper time are three important parameters with which to characterize the thermal properties of different stages during high-energy p-p, proton-nucleus (p-A) and nucleus-nucleus (A-A) collisions. The proper time reflects the time elapsed from the initial collisions to the stage of kinetic freezeout. A large proper time represents that the interacting system has a long life. Temperature is an important concept in the physics of high-energy collisions. Usually, papers [6,[10][11][12][13] use four types of temperature: initial temperature, chemical freeze-out temperature, kinetic freeze-out temperature and effective temperature. Firstly, the initial temperature describes the excitation degree of the interacting system at the initial stage of collisions. Secondly, the chemical freeze-out temperature describes the excitation degree of the interacting system at the stage of chemical equilibrium, where the chemical components (relative fractions) of particles are no longer changed. Thirdly, the kinetic freeze-out temperature describes the excitation degree of the interacting system at the stage of kinetic and thermal equilibrium in which the (transverse) momentum spectra of particles are no longer changed. Fourthly, the effective temperature is not a real temperature; it is related to particle mass and describes the sum of the excitation degree of the interacting system and the effect of transverse flow at the stage of kinetic freeze-out, where the transverse flow resulted from the impact and squeeze reflects the hydrodynamic expansion of the interacting system.
Generally, since the initial stage of collisions happens earlier than (or alongside) the stage of chemical freeze-out, the initial temperature is larger than (or equal to) the chemical freeze-out temperature. The chemical freeze-out temperature is equal to or larger than the kinetic freeze-out temperature, as the chemical equilibrium comes about at the same time as or before the kinetic equilibrium. It is usually obtained from the particle ratios [14][15][16]. The effective temperature is larger than the kinetic freeze-out temperature due to mass and flow effects [17,18]. It can be extracted from the transverse momentum spectra by using some distribution laws, such as the standard (Boltzmann, Fermi-Dirac and Bose-Einstein) distributions, the Tsallis distribution and so forth. The kinetic freeze-out temperature can be extracted from the transverse momentum spectra using models such as the blast-wave model with Boltzmann-Gibbs statistics, the blast-wave model with Tsallis statistics and the blast-wave model with fluctuations. Usually, the initial, chemical freeze-out and kinetic freeze-out temperatures are real temperatures.
We were very interested in studying the transverse excitation, the dynamic expansion characteristic and the lifetime of the interacting system, as their relation is important for mapping the phase diagram, and with them we could obtain some information about the proper time. In order to extract some quantities and study their dependences on event centrality and collision energy, we may use some models to analyze the particle spectra. These models include, but are not limited to, the blast-wave model with Boltzmann-Gibbs statistics [19,20], the blast-wave model with Tsallis statistics [21][22][23] and the alternative method [24][25][26][27][28] based on the standard distribution (including Boltzmann, Fermi-Dirac and Bose-Einstein distributions) or Tsallis distribution. The blast-wave model is a traditional and current method; it has very widespread applicability. The model makes the simple assumption that particles are locally thermalized in a hard-sphere, uniform-density source at a kinetic freeze-out temperature and are moving with a common collective manner in a transverse radial flow velocity field. In the alternative method, the kinetic freeze-out temperature is regarded as the intercept in the linear relation between effective temperature and rest mass, and the transverse flow velocity is regarded as the slope in the linear relation between mean transverse momenta and mean energy. Like the alternative method, too, the traditional blast-wave model can only extract the kinetic freeze-out temperature and transverse flow velocity. By inheritance and through development, the blast-wave model with fluctuations can extract the kinetic freeze-out temperature, transverse flow velocity and proper time simultaneously.
To understand the thermal properties of different stages of high energy collisions, we applied the blast-wave model with fluctuations to study the transverse momentum (p T ) spectra of different particles produced in A-A, p-A and p-p collisions at the LHC. The kinetic freeze-out temperature (T 0 ), transverse flow velocity (β T ) and uncorrected proper time (τ) were extracted from Pb-Pb collisions at 2.76 TeV; p-Pb collisions at 5.02 TeV; Xe-Xe collisions at 5.44 TeV; and p-p collisions at 0.9, 2.76, 5.02, 7 and 13 TeV. The cited experimental data were measured by the ALICE [29][30][31][32] and CMS Collaborations [33,34]. The changing trends of related parameters with event centrality and collision energy were then obtained and analyzed.
The paper is organized as follows. The formalism and method are briefly described in Section 2. Results and discussion are given in Section 3. Finally, we summarize our main observations and conclusions in Section 4.

Formalism and Method
High-energy collisions are complex processes in which many emission sources are formed. Sources with the same excitation degree may form a local equilibrium state which can be described by the standard distribution. For different equilibrium states which have different excitation degrees, different temperature parameters may be used. At the same time, we neglect the quantum effect and chemical potential due to their small influences on the p T spectra in high-energy collisions.
According to the invariant phase-space source emission distribution of Schnedermann, Sollfrank and Heinz (SSH) [19], Tomášik, Wiedemann and Heinz [35] further developed the blast-wave model with fluctuations. In this model the invariant momentum distribution is calculated by integrating over the space-time coordinates of the source function. According to [35,36], the blast-wave model with fluctuations results in the invariant momentum distribution as follows: where τ is the proper time considered by us as an uncorrected quantity, m T = p 2 T + m 2 0 is the transverse mass, m 0 is the rest mass, µ 0 is chemical potential-which is chosen to be 0 at high energy, T 0 is the kinetic freeze-out temperature of emission source, I 0 is the modified Bessel functions of the first kind, is the transverse flow rapidity, β(r) = β S (r/R) n 0 is a self-similar flow profile, β S is the flow velocity on the surface of the thermal source, r and R are, respectively, the radial position and its maximum in the thermal system, η s is the source rapidity, n 0 is a free parameter which is customarily chosen to be 2 [19] due to the quadratic profile resembling the solutions of hydrodynamics closest [37] and G(r) and H(η s ) are, respectively, the transverse and longitudinal source distributions. Generally, transverse flow velocity is In the case of n 0 = 2 as used in [19], we have β T = 0.5β S [38]. n 0 is not a sensitive quantity. It does not matter if n 0 = 1 or n 0 = 2. When we choose n 0 = 1, the results are similar. To be compatible with [35,36], we used n 0 = 2.
For a single source emission, one has η s = 0, G(r) = 1 and H(η s ) = 1. In the natural system of units, whereh = 1, Equation (1) is simplified as follows: The blast-wave model with fluctuations results in the united density function of y and p T as follows: According to the general expression of invariant momentum, the normalized factor in Equation (1) is τ/4π 2 = gV/(2π) 3 . For pseudoscalar particles π + , π − , K + and K − , the degeneracy factor is g = 1. For p and p, we have g = 2 due to their spin being 1/2. For sphere volume, V = (4/3)πR 3 , where R is the maximum range of light particles at the kinetic freeze-out. Then, we obtain the real or corrected proper time τ 0 = (1.5τ) 1/3 for emissions of π + , π − , K + and K − ; and τ 0 = (0.75τ) 1/3 for emissions of p andp. Figure 1 presents the transverse momentum spectra, (1/N EV )(2π p T ) −1 d 2 N/(dydp T ), of (a) π + , (b) π − , (c) K + , (d) K − , (e) p and (f)p produced in Pb-Pb collisions in different centralities at √ s NN = 2.76 TeV in the rapidity range |y| < 0.5. N EV which is usually omitted on the vertical axis denotes the number of events. The symbols represent the experimental data measured by the ALICE Collaboration in different centralities [29] and the curves are our fitting results by using the blast-wave model with fluctuations, i.e., Equation (4). In the fitting, the method of least squares was used to get the minimized χ 2 . The substantially near event centralities, the values of free parameters' kinetic freeze-out temperature T 0 , the transverse flow velocity β T and proper time τ; and the χ 2 degrees of freedom (dof) and corrected proper time τ 0 corresponding to the fit are listed in Table 1. In some cases, the χ 2 is large, which indicates that the fitting was qualitative and approximately acceptable, and the large dispersion between the curve and data exists. In most cases, one can see the good approximate descriptions of the model results for the experimental data of the ALICE Collaboration in the p T spectra of the identified particles produced in different centralities for Pb-Pb collisions at √ s NN = 2.76 TeV. Table 1. Values of parameters (T 0 , β T and τ), χ 2 , dof and τ 0 corresponding to the fits of the blast-wave model with fluctuations. The collision types, centrality and particle types are listed for clarity.

Figure Centrality
Particle
To see the dispersions of the curve from the data, Figure 1 continued (Figure 2) presents the ratios of data/fit. The panels a*-f* correspond to Figure 1a-f, respectively. The different closed (open) symbols represent the data/fit values corresponding to different centralities marked in Figure 2. Indeed, in some p T regions, the ratios are large. This means that the dispersions between the curve and data are large, and the fits are only qualitative and approximate in some cases.   Table 1.
One can see the good approximation the model produced of the experimental data of the ALICE Collaboration in the p T spectra of identified particles produced in different centrality p-Pb collisions at √ s NN = 5.02 TeV. Figure 4 is the same as Figure 2, but it presents the ratios of data/fit, in which panels a*-c* correspond to Figure 3a-c, respectively. The different closed (open) symbols represent the data/fit values corresponding to different centralities marked in the figure. One can see again that, in some p T regions, the ratios are large, which means large dispersions between the curve and data. In some cases, the fits are only qualitative and approximate.   Figures 1 and 3, but it shows the spectra, (1/N EV )d 2 N/(dydp T ), of (a) (π + + π − )/2, (b) (K + + K − )/2 and (c) (p +p)/2 produced in Xe-Xe collisions at √ s NN = 5.44 TeV in |y| < 0.5. The symbols represent the experimental data measured by the ALICE Collaboration [31]. The related parameters and the centralities are listed together in Table 1. One can see the good approximation the model produced of the experimental data of the ALICE Collaboration in the p T spectra of particles produced in different centrality Xe-Xe collisions at √ s NN = 5.44 TeV. Figure 6 is the same as Figures 2 and 4, but it presents the ratios of data/fit in which panels a*-c* correspond to Figure 5a-  represent the experimental data measured by the ALICE Collaboration in centralities 0-5%, 5-10%, 10-20%, 20-30%, 30-40%, 40-50%, 50-60%, 60-70% and 70-90% [31]. The spectra were scaled by the factor of 2 n ; n changed from 8 to 0 with the event centrality changing from 0-5% to 70-90%. Figure 7 is the same as Figure 5, but it shows the spectra of (a) π + , (b) π − , (c) K + , (d) K − , (e) p and (f)p produced in p-p collisions in |y| < 1. The symbols represent the experimental data measured by the CMS Collaboration at √ s = 0.9, 2.76, 7 and 13 TeV [33,34]. The energy and related parameters are listed in Table 2. One can see the good approximation the model produced of the experimental data of the CMS Collaboration in the p T spectra of identified particles produced in p-p collisions at different energies.  Table 2. Values of parameters (T 0 , β T and τ), χ 2 , dof and τ 0 corresponding to the fits of blast-wave model with fluctuations. The collision type, energy type and particle type are listed for clarity.   Figure 8 is the same as Figure 6, but it presents the ratios of data/fit, in which panels a*-f* correspond to Figure 7a Figure 5, but showing the spectra of (a) π + , (b) π − , (c) K + , (d) K − , (e) p and (f) p, produced in p-p collisions in |y| < 1. The symbols represent the experimental data measured by the CMS Collaboration at √ s = 0.9, 2.76, 7 and 13 TeV [33,34]; the spectra for the same particle in different centralities were multiplied by different amounts shown in the panels for clarity. The spectra were scaled by a factor of 2 n ; n changed from 0 to 3 with the energy changing from 0.9 to 2.76 TeV. Figure 9 is the same as Figure 1, but it shows the spectra of (a) (π + + π − )/2, (b) (K + + K − )/2 and (c) (p +p)/2 produced in p-p collisions in the pseudorapidity range |η| < 0.8. The symbols represent the experimental data measured by the ALICE Collaboration at √ s = 5.02 TeV [34]. The energy and related parameters are listed in Table 2. One can see the good approximation the model produced of the experimental data of the ALICE Collaboration in the p T spectra of identified particles produced in p-p collisions at √ s = 5.02 TeV. Figure 10 is the same as Figure 2, but it presents the ratios of data/fit corresponding to Figure 9. The different closed symbols represent the data/fit values corresponding to different particles marked in the figure. The same conclusions obtained in the above discussions can be obtained from Figure 10.   Figure 1, but showing the spectra of (π + + π − )/2, (K + + K − )/2 and (p +p)/2 produced in p-p collisions in the pseudorapidity range |η| < 0.8. The symbols represent the experimental data measured by the ALICE Collaboration at √ s = 5.02 TeV [32]. From Figures 1, 3, 5, 7 and 9, one can see that we have used the species-dependent parameters to fit the spectra of pions, kaons and protons. This is not the usual way to use the blast-wave model, which fits various spectra simultaneously. We have examined the simultaneous fit of the model for various spectra and know that narrow and different p T ranges have to be used for different particles. We do not think that the simultaneous fitting of various spectra can be better in wide p T ranges. Instead, we may use the individual fit for different spectra and obtain better fits. If needed, we may use the parameters averaged by weighting the yields of different particles to give the simultaneous fit. In fact, the simultaneous fit was not needed in this case, though what we show by the end of this section is a weighted average of proper times for different particles. As will be seen, the individual fit reveals and confirms the mass-dependent differential kinetic freezeout scenario.
From the above discussions we know that, according to the blast-wave model with fluctuations, we can obtain not only the kinetic freeze-out temperature and transverse flow velocity but also the proper time. The advantage of this model is that the information about proper time can be obtained. This is very important for us to understand the specific evolution of high-energy collisions. The proper time reflects the lifetime of the interacting system. We may discuss the change rule of the proper time by this model. This is an important innovation of the present work. However, the disadvantage of the model is also obvious. In some cases, one component model cannot fit the data well. In fact, it is only applicable to the low-p T region, but not to the high-p T region. This problem is not only a disadvantage in the blast-wave model with fluctuations, but also a disadvantage in all thermal models.
To study the change trend of parameters with centrality and energy, Figure 11 shows the dependencies of kinetic freeze-out temperature T 0 on (a)-(c) centrality C and (d) energy √ s for the production of different particles in (a) Pb-Pb, (b) p-Pb, (c) Xe-Xe and (d) p-p collisions. Different symbols represent different meanings shown in the panels. One can see that, from peripheral to central collisions, T 0 slightly increases. The reasons are the more violent interactions in central collisions, where there are higher degrees of excitation, and more participant nucleons being involved in the collisions. T 0 also increases with the increase of particle mass. This is an evidence of a mass-dependent differential kinetic freeze-out scenario or a multiple kinetic freeze-out scenario [12,39]. In addition, with the increase of energy, T 0 is nearly invariant at the LHC energies, which implies the nearly saturated density in QGP phase. In p-p collisions, the dependence of T 0 on energy is similar to that in peripheral nuclear (A-A and p-A) collisions.  Figure 11. Dependencies of kinetic freeze-out temperature T 0 on (a-c) centrality C and (d) energy √ s for the productions of (a) π + , π − , K + , K − , p andp in different centrality Pb-Pb collisions at 2.76 TeV; (b) (π + + π − )/2, (K + + K − )/2 and (p +p)/2 in different centrality p-Pb collisions at 5.02 TeV; (c) (π + + π − )/2, (K + + K − )/2 and (p +p)/2 in different centrality Xe-Xe collisions at 5.44 TeV; and (d) π + , π − , K + , K − , p andp at 0.9, 2.76, 7 and 13 TeV, along with (π + + π − )/2, (K + + K − )/2 and (p +p)/2 at 5.02 TeV. The red closed and black open symbols in panel (a) represent positively and negatively charged particles respectively, which are quoted from Table 1. The blue closed symbols in panels (b,c) represent positively plus negatively charged particles divided by 2, which are quoted from Table 1. The red closed (black open) and blue closed symbols in panel (d) represent positively (negatively) and positively plus negatively charged particles divided by 2, which are quoted from Table 2. Figure 12 is the same as Figure 11, but it shows the dependence of transverse flow velocity β T on (a)-(c) centrality C and (d) energy √ s. Different symbols represent different β T of different particles. One can see in general the slight increase of β T from peripheral to central collisions. With the increase of energy, β T increases slightly. The present work confirms our previous work [27,40,41] which used the intercept-slope method and obtained a slight larger β T in central collisions than in peripheral collisions. In our opinion, the flow is produced in the inner core of the interacting system. Even for peripheral or protonproton collisions, there is non-zero flow velocity. From Figure 12, one can also see the mass-dependent β T in most cases. A heavy particle corresponds to a small β T due to its large inertia. The differences in β T for different particles decrease with the increase of √ s, where β T is large at high √ s. This reflects that we can neglect the mass effect in a strong flow field.  In the fit process, the free parameters T 0 and β T were not fixed. We looked for the most appropriate parameters according to the values of χ 2 in various cases. That is, we obtained the results without bias. As the two parameters were sensitive to the results, they had far less freedom to vary. However, it was inevitable that there would be a correlation between the parameters. The values of T 0 and β T are not independent of centrality classes of up to 60% centrality in nuclear collisions. The same T 0 and β T indicate that the given interaction system has the same excitation degree, expansion characteristic and matter form from 0 to 60% centrality. The same matter can be QGP from the considered central to semi-central collisions. Figure 13 is the same as Figure 11, but it shows the dependencies of proper time τ 0 on (a)-(c) centrality C and (d) energy √ s. Different symbols represent different τ 0 for different particles. One can see the increase of τ 0 from peripheral to central collisions in panels (a-c), as the number of participant nucleons increases from peripheral to central collisions. Due to the large number of binary collisions by the re-scattering of partons, the system with more participants reached equilibrium slowly. However, the small number of participant nucleons led to the system to go quickly to equilibrium.  Figure 13. Same as Figure 11, but showing the dependencies of proper time τ 0 on (a-c) centrality C and (d) energy √ s. Figure 13 also proves the mass-dependent differential scenario from τ 0 . The heavier the particle is, the smaller τ 0 is, which shows the early freeze-out of heavier particles as compared to the lighter particles. This also suggests that different particles have different freeze-out surfaces [39,42]. The result that pions correspond to a much larger τ 0 than protons means that the protons stop interacting while pions are still interacting. Due to protons having larger m 0 than pions, protons are left behind during the system evolution from the origin of collisions in the radial direction, which is the behavior of hydrodynamics [43]. In Figure 13d, one can see, with the increase of energy, that τ 0 increases slightly or is nearly invariant. This is consistent with Figures 11d and 12d.
To study further the dependence of proper time on centrality C and energy √ s, Figure 14a shows the dependencies of the average proper time τ 0 and the average rootmean-square proper time τ 2 0 on C. The red closed and black open symbols represent τ 0 and τ 2 0 , respectively, which are from Pb-Pb, p-Pb and Xe-Xe collisions and are also listed in Table 3 with the impact parameter b. These averages were obtained by different particle weights due to different contribution fractions of π ± , K ± , p andp. One can see the  and blue open symbols represent τ 0 and τ 2 0 , respectively, from p-p collisions. The red and blue closed symbols represent τ 0 and τ 2 0 , respectively, from Pb-Pb, p-Pb and Xe-Xe collisions, which were calculated from Table 3 in the range of b = 0-8 fm by different weights (∝ b) due to different contribution fractions. The cases of large b were not included to reduce the influence of the cold nuclear effect. The results for p-p collisions corresponding to red and blue open symbols were calculated by different particle weights. The values of proper times used in Figure 14b are also listed in Table 4. In p-p collisions, one can see that with an increase of energy, τ 0 and τ 2 0 are slightly increscent or nearly invariant. In nuclear collisions, at similar energies, the considered proper times for Pb-Pb collisions are larger than those for p-Pb collisions, and the latter are larger than Xe-Xe collisions. Generally, nuclear collisions correspond to larger proper times than p-p collisions.   Table 3. In panel (b), the red and blue open symbols represent τ 0 and τ 2 0 , respectively, from p-p collisions, which are listed in Table 4. The red and blue closed symbols represent τ 0 and τ 2 0 respectively from Pb-Pb, p-Pb and Xe-Xe collisions, which are also listed in Table 4.  Figure 14a shows that the order of τ 0 and τ 2 0 is Pb-Pb > p-Pb > Xe-Xe in central collisions. However, τ 0 and τ 2 0 in peripheral p-Pb collisions are larger than those in peripheral Pb-Pb and Xe-Xe collisions. This is caused by the cold nuclear effect in A-A collisions. In order to avoid the cold nuclear effect, we considered only the weighted average of the events from central to semi-central Pb-Pb and Xe-Xe collisions. This was the case for A-A collisions in Figure 14b, for which one can see the order of τ 0 and τ 2 0 was Pb-Pb > p-Pb > Xe-Xe > p-p. This order is concordant with that of the maximum size between the projectile and target, and that of the other size if the maximum size is the same. This is direct evidence for the statement that the maximum size of the nucleus determines the lifetime of the interacting system, and the other size also affects it. The larger the maximum size, the longer the system lifetime. If the maximum sizes are the same, the larger other size will determine the longer system lifetime. To see the common property of the symmetric and asymmetric collisions, Figure 15 shows the dependencies of τ 0 and τ 2 0 on impact parameter b. The symbols have the same meaning as Figure 14a. One can see that the large differences in Figure 14a disappear in Figure 15. No matter the symmetric and asymmetric collisions, τ 0 and τ 2 0 decrease with increases of impact parameter b. In the plot of τ 0 and τ 2 0 versus impact parameter, there is no obvious difference between symmetric and asymmetric collisions. We may say that, compared with the centrality-dependent proper times, the results of impact parameter-dependent proper times are more concordant among different collisions.

Summary and Conclusions
We summarize here our main observations and conclusions. (a) We analyzed the transverse momentum spectra of π ± , K ± , p andp produced in different centrality Pb-Pb collisions at √ s NN = 2.76 TeV, p-Pb collisions at √ s NN = 5.02 TeV and Xe-Xe collisions at √ s NN = 5.44 TeV. The experimental data measured by the ALICE Collaborations were approximately fitted by the blast-wave model with fluctuations. Meanwhile, the transverse momentum spectra of π ± , K ± , p andp produced in p-p collisions at √ s NN = 0.9, 2.76, 5.02, 7 and 13 TeV were analyzed. The experimental data measured by the ALICE and CMS Collaborations were approximately fitted by the model. The kinetic freeze-out temperature, transverse flow velocity and proper time were extracted. (b) The kinetic freeze-out temperature and transverse flow velocity increase or do not vary with increases of event centrality from peripheral to central collisions, indicating higher excitation degree and greater expansion velocity of the system in central collisions. The average proper time increases with increasing event centrality from peripheral to central collisions, indicating longer lifetime of the system in central collisions. This means that in central collisions, the system needs longer average proper time to reach equilibrium. The heavier particles correspond to shorter proper time, showing the early freeze-out of heavier particles as compared to the lighter particles. From kinetic freeze-out temperature, transverse flow velocity and proper time, the mass-dependent differential kinetic freeze-out scenario or multiple kinetic freeze-out scenario is confirmed.

Data Availability Statement:
The data used to support the findings of this study are included within the article and are cited at relevant places within the text as references.