System size dependence of thermodynamic variables at kinetic freeze-out in high-energy collisions using the Tsallis distribution
Abstract
We use a thermodynamically consistent form of Tsallis distribution to study the dependence of various thermodynamic quantities on the system size in high-energy collisions. The charged hadron spectra obtained in +, +Pb, Xe+Xe, and Pb+Pb collisions at LHC are used to determine the energy density, pressure, particle density, entropy density, mean free path, Knudsen number, heat capacity, isothermal compressibility, expansion coefficient, and speed of sound at the kinetic freeze-out surface. These quantities are studied as a function of the system size. Notably, the rate of increase (or decrease) in these thermodynamic variables is found to be more rapid in small systems such as + and +Pb collisions than in large systems such as Xe+Xe and Pb+Pb collisions. This may be due to the small volume of the hadronic system in small collision systems at kinetic freeze-out. It is observed that high-multiplicity + collisions produce similar thermodynamic conditions as peripheral heavy-ion collisions at kinetic freeze-out.
I Introduction
High-energy collisions provide a unique opportunity to study the behavior of matter at extraordinarily high temperatures and densities. In the last few decades, considerable progress has been made in the study of the properties of a new state of matter made up of de-confined quarks and gluons called Quark-Gluon Plasma (QGP). This state is believed to have formed when normal baryonic matter undergoes a first-order or crossover phase transition at very high temperature and/or high baryon density [1, 2, 3]. Previously, it was believed that QGP is only formed in heavy-ion collisions; however, recent results show that high-multiplicity + collisions behave similarly to heavy-ion collisions, posing the question of whether QGP droplets are formed in such collisions [4, 5, 6, 7]. Consequently, there is a pressing need to gain deeper insights into high-multiplicity + collisions in order to understand the properties of the produced system.
To better understand the behavior of matter resulting from high-energy collisions, it is crucial to comprehensively understand its thermodynamic properties [8, 9, 10, 11, 12, 13, 14, 15]. In this study, we have explored various thermodynamic properties of the hadronic medium produced in both small collision systems (+ and +Pb) and large collision systems (Xe+Xe and Pb+Pb) at LHC energies.
We examine a scenario in which the fireball resulting from high-energy collisions cools down as it undergoes spacetime expansion. The final-state particles decouple from the system after Kinetic Freeze-Out (KFO), which represents the surface of the last elastic scattering of the produced hadrons. The production of final-state particles in abundance warrants the use of statistical methods to describe such systems. One such approach involves studying the transverse momentum () spectra of the produced particles in the framework of Tsallis statistics. Tsallis distribution was initially introduced as a generalization of the Boltzmann-Gibbs distribution and it accounts for deviations from the standard equilibrium statistical description [16]. Over time, Tsallis distribution has proven to be a valuable tool for describing the hadron spectra in high-energy collisions. The Tsallis distribution function is defined as [17]:
(1) |
where is the energy of the particle, is the entropy index, is the temperature, and is the chemical potential. is indicative of the degree of deviation from equilibrium, with representing complete equilibrium. Although various forms of the Tsallis distribution have been used to study the hadron spectra, ensuring thermodynamic consistency in their formulations poses a significan challenge [18, 19, 20, 21]. In the present study, we have employed Eq. (1), which has been demonstrated to maintain thermodynamic consistency [17]. The thermodynamic quantities of interest were derived using Eq. (1) in accordance with the following relations [22, 23, 24]:
(2) | ||||
(3) | ||||
(4) | ||||
(5) |
where is the degeneracy factor, is the Tsallis volume; and , , and correspond to entropy density, number density, and energy density, respectively. Tsallis volume, , should not be identified with the volume determined from femtoscopy [25]. Additionally, we have also studied thermodynamic response functions, such as heat capacity (), isothermal compressibility (), and expansion coefficient (), along with mean free path (), Knudsen number (), and squared speed of sound () in the produced hadronic medium. The primary objective of this work is to understand the interplay between collision energy and system size on various thermodynamic quantities, contrasting their behavior in small and large systems.
The remainder of this paper is organized as follows: In Section II, we discuss the Tsallis distribution that is used to fit the transverse momentum spectra of the charged hadrons. In Section III, we calculate the various thermodynamic quantities and discuss their dependence on charged particle multiplicity. Finally, in Section IV, we summarize our findings.
II Transverse momentum spectra
The -spectra of different particle species can be accurately described using the Tsallis distribution. In this framework, the invariant yield of a particle is defined as:
(6) |
This expression can be represented , transverse mass (, where is the particle mass), and rapidity () as follows:
(7) |
Assuming at LHC energies, the expression for the invariant yield at mid-rapidity () becomes
(8) |
In the present study, since we are studying charged hadron spectra in various collision systems, therefore, a transformation from rapidity () to pseudorapidity () is necessary. This is achieved using the following relation:
(9) |
which simplifies to:
(10) |
at mid-rapidity.
The charged hadron spectra from various collision systems were fitted with the sum of Tsallis distributions for pions, kaons, and protons. The resulting expression of the fit function is as follows [8]:
(11) |
where , and and are the degeneracy factors. Factor of accounts for the antiparticles. This notation is followed throughout the paper.
We have used Eq. (11) to fit the charged hadron spectra obtained in Pb+Pb collisions at = 2.76 and 5.02 TeV, Xe+Xe collisions at = 5.44 TeV, +Pb collisions at = 5.02 TeV, and + collisions at = 5.02 and 13 TeV, measured by the ALICE Collaboration [26, 27, 28]. The extracted fit parameters are used to calculate various thermodynamic variables. To achieve a fair comparison of these thermodynamic quantities across different collision systems, we have fitted the spectra measured within the same pseudorapidity window of 0.8. Furthermore, we have chosen to limit the fits to 5 GeV/c to investigate bulk properties, as high- particles are generated by hard processes. In a previous study ([8]), various thermodynamic quantities were calculated by fitting the spectra in the full range in Pb+Pb collisions at = 2.76 and 5.02 TeV using Tsallis distribution. It should be noted that the selection of the fitting range impacts the fit parameters and consequently influences the resulting thermodynamic variables [29].
Figure 1 shows the Tsallis fits to the charged hadron spectra in various centrality classes of Pb+Pb collisions at = 2.76 and 5.02 TeV using Eq. (11). It is evident that the Tsallis fits provide a more accurate description of the charged hadron spectra in peripheral Pb+Pb collisions compared to central collisions for both collision energies [8].
Similary, Fig. 2 presents the Tsallis distribution fits to the charged hadron spectra in different centrality classes of Xe+Xe collisions at = 5.44 TeV, and 0-100% centrality in +Pb collisions at = 5.02 TeV. As was the case in Pb+Pb collisions, Tsallis fits provide a more accurate description of the charged hadron spectra in peripheral Xe+Xe collisions compared to central collisions. We also observe that the Tsallis distribution provide a good description of the charged particle spectra in +Pb collisions.
Figure 3 shows Tsallis fits to the charged hadron spectra measured in different multiplicity classes of + collisions at = 5.02 and 13 TeV. Notably, Tsallis fits provide a better description of the charged hadron spectra in lower multiplicity classes in + collisions compared to higher multiplicity classes.
Figure 4 shows the dependence of fit parameters (, , and ) on charged particle multiplicity (), which is an indicator of the size of the hadronic system produced in high-energy collisions [27, 30, 28, 31]. The findings from Fig. 4 lead to the following conclusions:
-
•
The Tsallis volume, , increases linearly with increasing and hence is a good indicator of the system size. However, this increase appears to be slower for + collisions compared to Pb+Pb and Xe+Xe collisions. Similar trend was also observed in femtoscopy studies where small systems have smaller radii compared to large systems extracted with two-pion and three-pion cumulants [32]. This suggests that at the same , the volume of the hadronic system produced in + collisions is smaller than in Pb+Pb and Xe+Xe collisions. The single measurement in +Pb collisions falls within the trend observed in + collisions at the corresponding , suggesting a potential convergence of behavior in small collision systems.
-
•
The temperature, , increases with increasing . It is observed that rises sharply in Pb+Pb and Xe+Xe collisions compared to + collisions with increasing .
-
•
The entropy index, , for + collisions increases with increasing , suggesting a deviation from equilibrium in high-multiplicity + collisions. The -parameter for +Pb collisions aligns with the trend observed in + collisions at the corresponding . However, for Pb+Pb and Xe+Xe collisions, -parameter decreases with increasing , suggesting that the system approaches equilibrium more closely in central collisions than in peripheral collisions. It is speculated that the initial energy density induced hot spots caused by the Color Glass Condensate formalism may generate significant temperature fluctuations especially in small systems [33]. In high multiplicity + collisions this will cause an increase in the value of , contrary to central heavy-ion collisions. We also observe that increases with increasing center-of-mass energy in both small and large collision systems, indicating a tendency for the system to deviate further from equilibrium in collisions with higher energy.
The parameters obtained from the Tsallis fits are used to calculate various thermodynamic quantities, which are explored in the following section.
III Thermodynamic Variables
III.1 Energy
The Tsallis parameters, , , and are used to compute the energy density using the following equation:
(12) |
The left panel of Fig. 5 shows the energy density, , and the right panel shows the total energy, (= ) of the hadronic medium at kinetic freeze-out plotted as a function of in various collision systems.
It is observed that the energy density and total energy increase with increasing . Furthermore, we find that the energy density exhibits a steeper rise in small collision systems compared to large collision systems, leading to a discontinuity at 30. This difference might arise because the transverse overlap area in small systems is much smaller than that in large systems for the same event multiplicities [33, 34]. This is further corroborated by the fact that the discontinuity disappears when is plotted against and a linear scaling between the two is observed. This shows that high multiplicity + collisions are capable of producing energy densities that are comparable to the most peripheral heavy-ion collisions [35].
III.2 Pressure
Examining the pressure of hadrnoic medium produced in high energy collisions can provide valuable insights into its equation of state. It can be evaluated using the following equation:
(13) |
Figure 6 shows the pressure of the hadronic medium at kinetic freeze-out as a function of , exhibiting an increase with increasing . It is also observed that the pressure rises more rapidly in + and +Pb collisions than in Pb+Pb and Xe+Xe collisions, resulting to a discontinuity at 30. This can possibly be attributed to the higher initial densities in + and +Pb collisions, which may result in a larger pressure exerted against the surrounding environment [33, 36, 37].
Within Landau’s hydrodynamic model framework, the equation of state for an ultra-relativistic hadron gas is, , where is the square of the speed of sound [38]. Figure 7 shows the linear correlation between pressure and energy density. A linear fit to all the data points results in = 0.262 0.002. A more detailed discussion on is presented later.
III.3 Particle number
The particle number density can be evaluated from the expression:
(14) |
Figure 8 shows increase in the number density, , and total number of particles, (= ), with increasing . It is observed that the increase in particle density is more pronounced for + and +Pb collisions than for Pb+Pb and Xe+Xe collisions, which could be attributed to the smaller volume of the former collision systems. Moreover, a linear increase of with further strengthens the argument that particle number density is closely related to the volume of the system.
III.4 Entropy
The entropy plays a major role in understanding the evolution of the system produced in heavy-ion collisions. Entropy density can be evaluated using the following expression:
(15) |
In Fig. 9, we observe that the entropy density, , and total entropy, (= ), increase with increasing . This can be attributed to an increase in both the particle number and the volume of the system. Additionally, we note that the entropy density increases more rapidly in + and +Pb collisions than in Pb+Pb and Xe+Xe collisions, leading to a discontinuity at 30. It can be inferred from an almost linear relationship between and that the dissimilarity in the rate of increase may be ascribed to the volume of the hadronic system. We have also explicitly verified that the thermodynamic relation, , holds.
III.5 Mean free path and Knudsen number
The mean free path () of a system is defined as the average distance that a particle travels between two successive collisions in the system. Mathematically, the mean free path is expressed as [10]:
(16) |
where is the particle density, and represents the scattering cross-section. In the specific case of hard core hadron radius approximation, with fm, [10, 39, 40]. Furthermore, Knudsen number () is a dimensionless parameter, defined as the ratio of the mean free path () to the radius, (= ), of the hadronic system [14, 15]. A small value of signifies that the gas system is in the continuum region, allowing the application of hydrodynamics. Generally, when , the fluid-dynamical criterion is satisfied, indicating fluid-like behavior, whereas a value greater than unity indicates that the particles are mostly free-streaming [41, 42, 43].
Figure 10, shows both and decrease with increasing . The decrease in results from the increase in the particle number density with larger system sizes. Consequently, the hadrons experience more frequent collisions, leading to a decrease in the average distance between successive collisions. The decrease appears to be more rapid in + collisions because of the smaller volume of the hadronic system compared to large collision systems at the same charged particle multiplicity.
It was observed that, except for collisions with large , the Knudsen number () for the hadronic gas at KFO was greater than 1. This suggests that the mean free path () is comparable to the system size (), indicating that the hadronic system at KFO is not in the continuum phase. Instead, the particles are predominantly free-streaming.
III.6 Heat capacity
The heat capacity () of a system quantifies the amount of energy required to increase its temperature by one unit. Mathematically, this can be expressed as [10, 11]:
(17) |
Figure 11 shows the variation in as a function of . We observe an increase in with increasing . This could be a consequence of the increase in the particle number with . As the number of particles in a system increases, a higher amount of energy is required to increase its temperature by one unit.
We also observe a steeper increase in with in the case of + collisions, leading to a discontinuity at 30. This might be due to the higher value of the -parameter in + collisions. Further the system is away from equilibrium, as indicated by a higher -value, more is the change in its internal energy as temperature rises, resulting in a higher [10].
III.7 Isothermal compressiblity
The isothermal compressibility () of a system is a measure of the extent to which its volume changes in response to an external pressure. For an ideal fluid that is both incompressible and non-viscous, is . It can be evaluated from the expression [10, 12, 13]:
(18) |
Figure 12 shows isothermal compressibility as a function of . We observe that decreases with increasing , suggesting that the system approaches near ideal behavior with increasing multiplicity. This can be explained by the fact that higher multiplicity classes consist of a greater number of particles, thus requiring higher pressure to achieve a small change in volume [12].
III.8 Expansion coefficient
The expansion coefficient or isobaric expansivity () of a system is the measure of relative expansion per degree change in temperature at a constant pressure. For the hadronic medium at KFO, is calculated using the following expression [10]:
(19) |
Figure 13 shows as a function . The value of is observed to be negative in the studied collision systems. increases with increasing . This is a consequence of the increasing particle density with . For a denser system, the increase in temperature produces a small change in its volume, whereas for a less dense system, the change in volume is comparatively higher.
III.9 Speed of sound
The speed of sound in hadronic medium is closely related to its equation of state. The squared speed of sound is given as [11, 12]:
(20) |
For a non-interacting massless ideal gas, the expected value of is 1/3 times the speed of light squared [44]. Figure 14 shows plotted as a function of . shows a weak dependence with . We observe that for studied collision systems and energies, on an average = 0.266 0.001.
IV Summary
The charged hadron spectra measured by the ALICE Collaboration in Pb+Pb collisions at = 2.76 and 5.02 TeV, Xe+Xe collisions at = 5.44 TeV, +Pb collisions at = 5.02 TeV, and + collisions at = 5.02 and 13 TeV have been studied in the framework of Tsallis statistics. We find that the Tsallis fits provide a better description of charged hadron spectra in peripheral Pb+Pb and Xe+Xe collisions than in more central collisions. In addition, we also observe an improvement in the fitting for lower multiplicity classes in + collisions and a reasonably good description of the spectra in +Pb collisions. Furthermore, the dependence of the fit parameters (, , and ) on the was also noted. It is observed that:
-
•
Tsallis volume of the hadronic medium increases with increasing , but the growth rate of volume in + collisions is slower compared to Pb+Pb and Xe+Xe collisions, while +Pb collisions show a similar trend to + collisions. This suggests that small collision systems produce smaller hadronic systems at KFO for the same compared to large collision systems.
-
•
Temperature increases with increasing with a slow increase in small collision systems and a sharp rise in large collision systems.
-
•
The entropy index, , increases with increasing in collisions, suggesting a deviation from equilibrium in high multiplicity collisions. -parameter in +Pb collisions is observed to follow the same trend as collisions at the corresponding . However, for Pb+Pb and Xe+Xe collisions, decreases with increasing , indicating that the system is closer to equilibrium in central collisions.
The thermodynamic properties of the hadronic medium in high-energy collisions, focusing on the dependence of energy density, pressure, particle density, entropy density, mean free path, Knudsen number, heat capacity, isothermal compressibility, expansion coefficient, and speed of sound on are studied. We find that:
-
•
Energy density increases with increasing , with small collision systems having a higher energy density in the common region, possibly because of the smaller volumes of the hadronic medium compared to large collision systems.
-
•
Pressure increases with increasing , with small collision systems exhibiting a faster increase. This might be explained by the fact that less space for particles to move around in small collision systems results into a more rapid momentum exchange.
-
•
Number density increases with increasing , with smaller systems exhibiting a faster increase owing to the smaller volume of the produced hadronic system.
-
•
Entropy density increases with increasing , with smaller systems having a higher rate of increase.
-
•
The mean free path decreases with increasing due to the increased particle number density. Knudsen number was found to be much greater than unity, suggesting that produced hadronic medium is in the free streaming region.
-
•
Heat capacity increases with increasing . This is a consequence of increase in particle number with increasing . Increasing the temperature of a system by one unit will require more energy when the system contains a greater number of particles.
-
•
Isothermal compressibility decreases with increasing indicating near-ideal behavior with increasing multiplicity.
-
•
Expansion coefficient is negative for the studied collision systems and it increases with increasing . This is due to the higher particle density associated with greater , resulting in a smaller change in volume for the same amount of pressure in a relatively denser system.
-
•
Squared speed of sound shows a weak dependence with . We observe an average value of = 0.266 0.001 for the studied collisions systems and energies.
This study characterizes the thermodynamic properties of the produced hadronic system in different collisions across a range of LHC energies. We observe that high multiplicity + collisions produce a hadronic medium with thermodynamic properties similar to peripheral heavy-ion collisions. These observations may contribute to the efforts to understand the properties of the medium formed in high-multiplicity + collisions.
Acknowledgements.
CJ acknowledges the financial support from DAE-DST Project No. 3015/I/2021/Gen/RD-I/13283.References
- Ejiri [2008] S. Ejiri, Phys. Rev. D 78, 074507 (2008).
- Aoki et al. [2006] Y. Aoki, G. Endrődi, Z. Fodor, S. D. Katz, and K. K. Szabó, Nature 443, 675 (2006).
- Adams et al. [2005] J. Adams et al. (STAR Collaboration), Nucl. Phys. A 757, 102 (2005).
- Adam et al. [2017] J. Adam et al. (ALICE Collaboration), Nature Phys. 13, 535 (2017).
- Bjorken et al. [2013] J. D. Bjorken, S. J. Brodsky, and A. S. Goldhaber, Phys. Lett. B 726, 344 (2013).
- Khachatryan et al. [2017] V. Khachatryan et al. (CMS Collaboration), Phys. Lett. B 765, 193 (2017).
- Khuntia et al. [2019a] A. Khuntia, H. Sharma, S. K. Tiwari, R. Sahoo, and J. Cleymans, Eur. Phys. J. A 55 (2019a).
- Azmi et al. [2020] M. D. Azmi, T. Bhattacharyya, J. Cleymans, and M. Paradza, J. Phys. G: Nucl. Part. Phys 47, 045001 (2020).
- Deb et al. [2021] S. Deb, G. Sarwar, R. Sahoo, and J. e Alam, Eur. Phys. J. A 57 (2021).
- Sahu and Sahoo [2021a] D. Sahu and R. Sahoo, Physics 3, 207 (2021a).
- Khuntia et al. [2016] A. Khuntia, P. Sahoo, P. Garg, R. Sahoo, and J. Cleymans, Eur. Phys. J. A 52 (2016).
- Jain et al. [2023] S. Jain, R. Gupta, and S. Jena, Universe 9, 170 (2023).
- Khuntia et al. [2019b] A. Khuntia, S. K. Tiwari, P. Sharma, R. Sahoo, and T. K. Nayak, Phys. Rev. C 100, 014910 (2019b).
- Sarkar and Ghosh [2017] N. Sarkar and P. Ghosh, Phys. Rev. C 96, 044901 (2017).
- Sarkar et al. [2019] N. Sarkar, P. Deb, and P. Ghosh, arxiv:1905.06532v3 (2019).
- Tsallis [1988] C. Tsallis, J. Stat. Phys. 52, 479 (1988).
- Cleymans and Worku [2012a] J. Cleymans and D. Worku, J. Phys. G: Nucl. Part. Phys. 39, 025006 (2012a).
- Hui et al. [2018] J.-Q. Hui, Z.-J. Jiang, and D.-F. Xu, Adv. High Energy Phys. 2018, 7682325 (2018).
- Bíró et al. [2017] G. Bíró, G. Barnaföldi, T. Biró, K. Ürmössy, and Á. Takács, Entropy 19, 88 (2017).
- Si et al. [2018] R.-F. Si, H.-L. Li, and F.-H. Liu, Adv. High Energy Phys. 2018, 7895967 (2018).
- Lao et al. [2017] H.-L. Lao, F.-H. Liu, and R. A. Lacey, Eur. Phys. J. A 53 (2017).
- Cleymans and Worku [2012b] J. Cleymans and D. Worku, Eur. Phys. J. A 48 (2012b).
- Bhattacharyya et al. [2016] T. Bhattacharyya, J. Cleymans, and S. Mogliacci, Phys. Rev. D 94, 094026 (2016).
- Bhattacharyya et al. [2017] T. Bhattacharyya, J. Cleymans, P. Garg, P. Kumar, S. Mogliacci, R. Sahoo, and S. Tripathy, J. Phys.: Conf. Ser. 878, 012016 (2017).
- Cleymans and Paradza [2020] J. Cleymans and M. W. Paradza, Physics 2, 654 (2020).
- Acharya et al. [2018] S. Acharya et al. (ALICE Collaboration), J. High Energ. Phys. 2018, 13 (2018).
- Acharya et al. [2019a] S. Acharya et al. (ALICE Collaboration), Phys. Lett. B 788, 166 (2019a).
- Acharya et al. [2019b] S. Acharya et al. (ALICE Collaboration), Eur. Phys. J. C 79 (2019b).
- Patra et al. [2021] R. N. Patra, B. Mohanty, and T. K. Nayak, Eur. Phys. J. Plus 136 (2021).
- Abelev et al. [2013] B. Abelev et al. (ALICE Collaboration), Phys. Rev. C 88, 044910 (2013).
- Acharya et al. [2020] S. Acharya et al. (ALICE Collaboration), Phys. Lett. B 800, 135043 (2020).
- Grosse-Oetringhaus [2014] J. F. Grosse-Oetringhaus, Nuclear Physics A 931, 22 (2014), qUARK MATTER 2014.
- Liu et al. [2023] L. Liu, Z.-B. Yin, and L. Zheng, Chinese Phys. C 47, 024103 (2023).
- Loizides et al. [2018] C. Loizides, J. Kamin, and D. d’Enterria, Phys. Rev. C 97, 054910 (2018).
- Sahu and Sahoo [2021b] D. Sahu and R. Sahoo, J. Phys. G: Nucl. Part. Phys. 48, 125104 (2021b).
- Shuryak and Zahed [2013] E. Shuryak and I. Zahed, Phys. Rev. C 88, 044915 (2013).
- Castorina et al. [2020] P. Castorina, A. Iorio, D. Lanteri, H. Satz, and M. Spousta, Phys. Rev. C 101, 054902 (2020).
- Mohanty and Alam [2003] B. Mohanty and J.-e. Alam, Phys. Rev. C 68, 064903 (2003).
- Kadam and Mishra [2015] G. P. Kadam and H. Mishra, Phys. Rev. C 92, 035203 (2015).
- Bugaev et al. [2013] K. A. Bugaev, D. R. Oliinychenko, A. S. Sorin, and G. M. Zinovjev, Eur. Phys. J. A 49 (2013).
- Huovinen and Molnar [2009] P. Huovinen and D. Molnar, Phys. Rev. C 79, 014906 (2009).
- Bouras et al. [2010] I. Bouras, E. Molnár, H. Niemi, Z. Xu, A. El, O. Fochler, C. Greiner, and D. H. Rischke, Phys. Rev. C 82, 024910 (2010).
- Gombeaud et al. [2009] C. Gombeaud, T. Lappi, and J.-Y. Ollitrault, Phys. Rev. C 79, 054914 (2009).
- Ferrer and Hackebill [2023] E. Ferrer and A. Hackebill, Nuclear Physics A 1031, 122608 (2023).