Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
License: CC BY 4.0
arXiv:2401.13629v1 [hep-ph] 24 Jan 2024

System size dependence of thermodynamic variables at kinetic freeze-out in high-energy collisions using the Tsallis distribution

Rishabh Sharma rishabhsharma@students.iisertirupati.ac.in    Krishan Gopal krishangopal@students.iisertirupati.ac.in    Sharang Rav Sharma sharang.rav@students.iisertirupati.ac.in    Chitrasen Jena cjena@iisertirupati.ac.in Department of Physics, Indian Institute of Science Education and Research (IISER) Tirupati, Tirupati 517507, Andhra Pradesh, India
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 p𝑝pitalic_p+p𝑝pitalic_p, p𝑝pitalic_p+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 p𝑝pitalic_p+p𝑝pitalic_p and p𝑝pitalic_p+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 p𝑝pitalic_p+p𝑝pitalic_p collisions produce similar thermodynamic conditions as peripheral heavy-ion collisions at kinetic freeze-out.

Suggested keywords
preprint: APS/123-QED

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 p𝑝pitalic_p+p𝑝pitalic_p 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 p𝑝pitalic_p+p𝑝pitalic_p 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 (p𝑝pitalic_p+p𝑝pitalic_p and p𝑝pitalic_p+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 (pTsubscript𝑝𝑇p_{T}italic_p start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT) 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]:

f(E,q,T,μ)=[1+(q1)EμT]1q1,𝑓𝐸𝑞𝑇𝜇superscriptdelimited-[]1𝑞1𝐸𝜇𝑇1𝑞1f(E,q,T,\mu)=\left[1+(q-1)\frac{E-\mu}{T}\right]^{\frac{-1}{q-1}},italic_f ( italic_E , italic_q , italic_T , italic_μ ) = [ 1 + ( italic_q - 1 ) divide start_ARG italic_E - italic_μ end_ARG start_ARG italic_T end_ARG ] start_POSTSUPERSCRIPT divide start_ARG - 1 end_ARG start_ARG italic_q - 1 end_ARG end_POSTSUPERSCRIPT , (1)

where E𝐸Eitalic_E is the energy of the particle, q𝑞qitalic_q is the entropy index, T𝑇Titalic_T is the temperature, and μ𝜇\muitalic_μ is the chemical potential. q𝑞qitalic_q is indicative of the degree of deviation from equilibrium, with q=1𝑞1q=1italic_q = 1 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]:

Entropy: SEntropy: 𝑆\displaystyle\text{Entropy: }SEntropy: italic_S =sV=gVd3p(2π)3(ffq1qf),absent𝑠𝑉𝑔𝑉superscript𝑑3𝑝superscript2𝜋3𝑓superscript𝑓𝑞1𝑞𝑓\displaystyle=sV=-gV\int\frac{d^{3}p}{(2\pi)^{3}}\left(\frac{f-f^{q}}{1-q}-f% \right),= italic_s italic_V = - italic_g italic_V ∫ divide start_ARG italic_d start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT italic_p end_ARG start_ARG ( 2 italic_π ) start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT end_ARG ( divide start_ARG italic_f - italic_f start_POSTSUPERSCRIPT italic_q end_POSTSUPERSCRIPT end_ARG start_ARG 1 - italic_q end_ARG - italic_f ) , (2)
Number: NNumber: 𝑁\displaystyle\text{Number: }NNumber: italic_N =nV=gVd3p(2π)3fq,absent𝑛𝑉𝑔𝑉superscript𝑑3𝑝superscript2𝜋3superscript𝑓𝑞\displaystyle=nV=gV\int\frac{d^{3}p}{(2\pi)^{3}}f^{q},= italic_n italic_V = italic_g italic_V ∫ divide start_ARG italic_d start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT italic_p end_ARG start_ARG ( 2 italic_π ) start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT end_ARG italic_f start_POSTSUPERSCRIPT italic_q end_POSTSUPERSCRIPT , (3)
Energy: EEnergy: 𝐸\displaystyle\text{Energy: }EEnergy: italic_E =ϵV=gVd3p(2π)3Efq,absentitalic-ϵ𝑉𝑔𝑉superscript𝑑3𝑝superscript2𝜋3𝐸superscript𝑓𝑞\displaystyle=\epsilon V=gV\int\frac{d^{3}p}{(2\pi)^{3}}Ef^{q},= italic_ϵ italic_V = italic_g italic_V ∫ divide start_ARG italic_d start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT italic_p end_ARG start_ARG ( 2 italic_π ) start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT end_ARG italic_E italic_f start_POSTSUPERSCRIPT italic_q end_POSTSUPERSCRIPT , (4)
Pressure: PPressure: 𝑃\displaystyle\text{Pressure: }PPressure: italic_P =gd3p(2π)3p23Efq,absent𝑔superscript𝑑3𝑝superscript2𝜋3superscript𝑝23𝐸superscript𝑓𝑞\displaystyle=g\int\frac{d^{3}p}{(2\pi)^{3}}\frac{p^{2}}{3E}f^{q},= italic_g ∫ divide start_ARG italic_d start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT italic_p end_ARG start_ARG ( 2 italic_π ) start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT end_ARG divide start_ARG italic_p start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG 3 italic_E end_ARG italic_f start_POSTSUPERSCRIPT italic_q end_POSTSUPERSCRIPT , (5)

where g𝑔gitalic_g is the degeneracy factor, V𝑉Vitalic_V is the Tsallis volume; and s𝑠sitalic_s, n𝑛nitalic_n, and ϵitalic-ϵ\epsilonitalic_ϵ correspond to entropy density, number density, and energy density, respectively. Tsallis volume, V𝑉Vitalic_V, should not be identified with the volume determined from femtoscopy [25]. Additionally, we have also studied thermodynamic response functions, such as heat capacity (CVsubscript𝐶𝑉C_{V}italic_C start_POSTSUBSCRIPT italic_V end_POSTSUBSCRIPT), isothermal compressibility (κTsubscript𝜅𝑇\kappa_{T}italic_κ start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT), and expansion coefficient (α𝛼\alphaitalic_α), along with mean free path (λ𝜆\lambdaitalic_λ), Knudsen number (Kn𝐾𝑛Knitalic_K italic_n), and squared speed of sound (cs2superscriptsubscript𝑐𝑠2c_{s}^{2}italic_c start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT) 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 pTsubscript𝑝𝑇p_{T}italic_p start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT-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:

Ed3Nd3p=gV(2π)3E[1+(q1)EμT]qq1.𝐸superscript𝑑3𝑁superscript𝑑3𝑝𝑔𝑉superscript2𝜋3𝐸superscriptdelimited-[]1𝑞1𝐸𝜇𝑇𝑞𝑞1E\frac{d^{3}N}{d^{3}p}=\frac{gV}{(2\pi)^{3}}E\left[1+(q-1)\frac{E-\mu}{T}% \right]^{\frac{-q}{q-1}}.italic_E divide start_ARG italic_d start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT italic_N end_ARG start_ARG italic_d start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT italic_p end_ARG = divide start_ARG italic_g italic_V end_ARG start_ARG ( 2 italic_π ) start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT end_ARG italic_E [ 1 + ( italic_q - 1 ) divide start_ARG italic_E - italic_μ end_ARG start_ARG italic_T end_ARG ] start_POSTSUPERSCRIPT divide start_ARG - italic_q end_ARG start_ARG italic_q - 1 end_ARG end_POSTSUPERSCRIPT . (6)

This expression can be represented pTsubscript𝑝𝑇p_{T}italic_p start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT, transverse mass (mT=pT2+m2subscript𝑚𝑇superscriptsubscript𝑝𝑇2superscript𝑚2m_{T}=\sqrt{p_{T}^{2}+m^{2}}italic_m start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT = square-root start_ARG italic_p start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + italic_m start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG, where m𝑚mitalic_m is the particle mass), and rapidity (y𝑦yitalic_y) as follows:

d2NpTdpTdy=gV(2π)2mTcosh(y)superscript𝑑2𝑁subscript𝑝𝑇𝑑subscript𝑝𝑇𝑑𝑦𝑔𝑉superscript2𝜋2subscript𝑚𝑇𝑦\displaystyle\frac{d^{2}N}{p_{T}dp_{T}dy}=\frac{gV}{(2\pi)^{2}}m_{T}\cosh(y)divide start_ARG italic_d start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_N end_ARG start_ARG italic_p start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT italic_d italic_p start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT italic_d italic_y end_ARG = divide start_ARG italic_g italic_V end_ARG start_ARG ( 2 italic_π ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG italic_m start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT roman_cosh ( italic_y )
×[1+(q1)mTcosh(y)μT]qq1.absentsuperscriptdelimited-[]1𝑞1subscript𝑚𝑇𝑦𝜇𝑇𝑞𝑞1\displaystyle\ \times\left[1+(q-1)\frac{m_{T}\cosh(y)-\mu}{T}\right]^{\frac{-q% }{q-1}}.× [ 1 + ( italic_q - 1 ) divide start_ARG italic_m start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT roman_cosh ( italic_y ) - italic_μ end_ARG start_ARG italic_T end_ARG ] start_POSTSUPERSCRIPT divide start_ARG - italic_q end_ARG start_ARG italic_q - 1 end_ARG end_POSTSUPERSCRIPT . (7)

Assuming μ=0𝜇0\mu=0italic_μ = 0 at LHC energies, the expression for the invariant yield at mid-rapidity (y=0𝑦0y=0italic_y = 0) becomes

d2NdpTdy=gV(2π)2mTpT[1+(q1)mTT]qq1.superscript𝑑2𝑁𝑑subscript𝑝𝑇𝑑𝑦𝑔𝑉superscript2𝜋2subscript𝑚𝑇subscript𝑝𝑇superscriptdelimited-[]1𝑞1subscript𝑚𝑇𝑇𝑞𝑞1\frac{d^{2}N}{dp_{T}dy}=\frac{gV}{(2\pi)^{2}}m_{T}p_{T}\left[1+(q-1)\frac{m_{T% }}{T}\right]^{\frac{-q}{q-1}}.divide start_ARG italic_d start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_N end_ARG start_ARG italic_d italic_p start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT italic_d italic_y end_ARG = divide start_ARG italic_g italic_V end_ARG start_ARG ( 2 italic_π ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG italic_m start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT italic_p start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT [ 1 + ( italic_q - 1 ) divide start_ARG italic_m start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT end_ARG start_ARG italic_T end_ARG ] start_POSTSUPERSCRIPT divide start_ARG - italic_q end_ARG start_ARG italic_q - 1 end_ARG end_POSTSUPERSCRIPT . (8)

In the present study, since we are studying charged hadron spectra in various collision systems, therefore, a transformation from rapidity (y𝑦yitalic_y) to pseudorapidity (η𝜂\etaitalic_η) is necessary. This is achieved using the following relation:

dNdpTdη=dNdpTdy(1m2mT2cosh2(y)),𝑑𝑁𝑑subscript𝑝𝑇𝑑𝜂𝑑𝑁𝑑subscript𝑝𝑇𝑑𝑦1superscript𝑚2superscriptsubscript𝑚𝑇2superscript2𝑦\frac{dN}{dp_{T}d\eta}=\frac{dN}{dp_{T}dy}\left(\sqrt{1-\frac{m^{2}}{m_{T}^{2}% \cosh^{2}(y)}}\right),divide start_ARG italic_d italic_N end_ARG start_ARG italic_d italic_p start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT italic_d italic_η end_ARG = divide start_ARG italic_d italic_N end_ARG start_ARG italic_d italic_p start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT italic_d italic_y end_ARG ( square-root start_ARG 1 - divide start_ARG italic_m start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG italic_m start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT roman_cosh start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( italic_y ) end_ARG end_ARG ) , (9)

which simplifies to:

dNdpTdη=dNdpTdy(pTmT),𝑑𝑁𝑑subscript𝑝𝑇𝑑𝜂𝑑𝑁𝑑subscript𝑝𝑇𝑑𝑦subscript𝑝𝑇subscript𝑚𝑇\frac{dN}{dp_{T}d\eta}=\frac{dN}{dp_{T}dy}\left(\frac{p_{T}}{m_{T}}\right),divide start_ARG italic_d italic_N end_ARG start_ARG italic_d italic_p start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT italic_d italic_η end_ARG = divide start_ARG italic_d italic_N end_ARG start_ARG italic_d italic_p start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT italic_d italic_y end_ARG ( divide start_ARG italic_p start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT end_ARG start_ARG italic_m start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT end_ARG ) , (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]:

d2NdpTdη=2V(2π)2pT2i=13gi[1+(q1)mT,iT]qq1,superscript𝑑2𝑁𝑑subscript𝑝𝑇𝑑𝜂2𝑉superscript2𝜋2superscriptsubscript𝑝𝑇2superscriptsubscript𝑖13subscript𝑔𝑖superscriptdelimited-[]1𝑞1subscript𝑚𝑇𝑖𝑇𝑞𝑞1\frac{d^{2}N}{dp_{T}d\eta}=2\frac{V}{(2\pi)^{2}}p_{T}^{2}\sum_{i=1}^{3}g_{i}% \left[1+(q-1)\frac{m_{T,i}}{T}\right]^{\frac{-q}{q-1}},divide start_ARG italic_d start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_N end_ARG start_ARG italic_d italic_p start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT italic_d italic_η end_ARG = 2 divide start_ARG italic_V end_ARG start_ARG ( 2 italic_π ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG italic_p start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ∑ start_POSTSUBSCRIPT italic_i = 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT italic_g start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT [ 1 + ( italic_q - 1 ) divide start_ARG italic_m start_POSTSUBSCRIPT italic_T , italic_i end_POSTSUBSCRIPT end_ARG start_ARG italic_T end_ARG ] start_POSTSUPERSCRIPT divide start_ARG - italic_q end_ARG start_ARG italic_q - 1 end_ARG end_POSTSUPERSCRIPT , (11)

where i=π+,K+,p𝑖superscript𝜋superscript𝐾𝑝i=\pi^{+},K^{+},pitalic_i = italic_π start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT , italic_K start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT , italic_p, and gπ=gK=1subscript𝑔𝜋subscript𝑔𝐾1g_{\pi}=g_{K}=1italic_g start_POSTSUBSCRIPT italic_π end_POSTSUBSCRIPT = italic_g start_POSTSUBSCRIPT italic_K end_POSTSUBSCRIPT = 1 and gp=2subscript𝑔𝑝2g_{p}=2italic_g start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT = 2 are the degeneracy factors. Factor of 2222 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 sNNsubscript𝑠𝑁𝑁\sqrt{s_{NN}}square-root start_ARG italic_s start_POSTSUBSCRIPT italic_N italic_N end_POSTSUBSCRIPT end_ARG = 2.76 and 5.02 TeV, Xe+Xe collisions at sNNsubscript𝑠𝑁𝑁\sqrt{s_{NN}}square-root start_ARG italic_s start_POSTSUBSCRIPT italic_N italic_N end_POSTSUBSCRIPT end_ARG = 5.44 TeV, p𝑝pitalic_p+Pb collisions at sNNsubscript𝑠𝑁𝑁\sqrt{s_{NN}}square-root start_ARG italic_s start_POSTSUBSCRIPT italic_N italic_N end_POSTSUBSCRIPT end_ARG = 5.02 TeV, and p𝑝pitalic_p+p𝑝pitalic_p collisions at s𝑠\sqrt{s}square-root start_ARG italic_s end_ARG = 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 |η|<𝜂absent|\eta|<| italic_η | < 0.8. Furthermore, we have chosen to limit the fits to pT<subscript𝑝𝑇absentp_{T}<italic_p start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT < 5 GeV/c to investigate bulk properties, as high-pTsubscript𝑝𝑇p_{T}italic_p start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT particles are generated by hard processes. In a previous study ([8]), various thermodynamic quantities were calculated by fitting the spectra in the full pTsubscript𝑝𝑇p_{T}italic_p start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT range in Pb+Pb collisions at sNNsubscript𝑠𝑁𝑁\sqrt{s_{NN}}square-root start_ARG italic_s start_POSTSUBSCRIPT italic_N italic_N end_POSTSUBSCRIPT end_ARG = 2.76 and 5.02 TeV using Tsallis distribution. It should be noted that the selection of the pTsubscript𝑝𝑇p_{T}italic_p start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT fitting range impacts the fit parameters and consequently influences the resulting thermodynamic variables [29].

Refer to caption
Refer to caption
Figure 1: Tsallis fits to the pTsubscript𝑝𝑇p_{T}italic_p start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT-spectra of charged hadrons in different centrality classes of Pb+Pb collisions at sNNsubscript𝑠𝑁𝑁\sqrt{s_{NN}}square-root start_ARG italic_s start_POSTSUBSCRIPT italic_N italic_N end_POSTSUBSCRIPT end_ARG = 2.76 TeV (left) and 5.02 TeV (right). The lower panel of each plot shows the ratio of the data and fit function.

Figure 1 shows the Tsallis fits to the charged hadron spectra in various centrality classes of Pb+Pb collisions at sNNsubscript𝑠𝑁𝑁\sqrt{s_{NN}}square-root start_ARG italic_s start_POSTSUBSCRIPT italic_N italic_N end_POSTSUBSCRIPT end_ARG = 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].

Refer to caption
Refer to caption
Figure 2: Tsallis fits to the pTsubscript𝑝𝑇p_{T}italic_p start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT-spectra of charged hadrons in different collision centralities of Xe+Xe collisions (left) and 0-100% centrality of p𝑝pitalic_p+Pb collisions (right) at sNNsubscript𝑠𝑁𝑁\sqrt{s_{NN}}square-root start_ARG italic_s start_POSTSUBSCRIPT italic_N italic_N end_POSTSUBSCRIPT end_ARG = 5.02 TeV. The lower panel of each plot shows the ratio of the data and fit function.

Similary, Fig. 2 presents the Tsallis distribution fits to the charged hadron spectra in different centrality classes of Xe+Xe collisions at sNNsubscript𝑠𝑁𝑁\sqrt{s_{NN}}square-root start_ARG italic_s start_POSTSUBSCRIPT italic_N italic_N end_POSTSUBSCRIPT end_ARG = 5.44 TeV, and 0-100% centrality in p𝑝pitalic_p+Pb collisions at sNNsubscript𝑠𝑁𝑁\sqrt{s_{NN}}square-root start_ARG italic_s start_POSTSUBSCRIPT italic_N italic_N end_POSTSUBSCRIPT end_ARG = 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 p𝑝pitalic_p+Pb collisions.

Refer to caption
Refer to caption
Figure 3: Tsallis fits to the pTsubscript𝑝𝑇p_{T}italic_p start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT-spectra of charged hadrons in different multiplicity classes of p𝑝pitalic_p+p𝑝pitalic_p collisions at sNNsubscript𝑠𝑁𝑁\sqrt{s_{NN}}square-root start_ARG italic_s start_POSTSUBSCRIPT italic_N italic_N end_POSTSUBSCRIPT end_ARG = 5.02 TeV (left) and 13 TeV (right). The lower panel of each plot shows the ratio of the data and fit function.

Figure 3 shows Tsallis fits to the charged hadron spectra measured in different multiplicity classes of p𝑝pitalic_p+p𝑝pitalic_p collisions at s𝑠\sqrt{s}square-root start_ARG italic_s end_ARG = 5.02 and 13 TeV. Notably, Tsallis fits provide a better description of the charged hadron spectra in lower multiplicity classes in p𝑝pitalic_p+p𝑝pitalic_p collisions compared to higher multiplicity classes.

Refer to caption
Refer to caption
Refer to caption
Figure 4: Tsallis parameters (volume, temperature, and q𝑞qitalic_q) as a function of dNch/dηdelimited-⟨⟩𝑑subscript𝑁𝑐𝑑𝜂\langle dN_{ch}/d\eta\rangle⟨ italic_d italic_N start_POSTSUBSCRIPT italic_c italic_h end_POSTSUBSCRIPT / italic_d italic_η ⟩, obtained by fitting the charged hadron pTsubscript𝑝𝑇p_{T}italic_p start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT-spectra in Pb+Pb collisions at sNNsubscript𝑠𝑁𝑁\sqrt{s_{NN}}square-root start_ARG italic_s start_POSTSUBSCRIPT italic_N italic_N end_POSTSUBSCRIPT end_ARG = 2.76 and 5.02 TeV, Xe+Xe collisions at sNNsubscript𝑠𝑁𝑁\sqrt{s_{NN}}square-root start_ARG italic_s start_POSTSUBSCRIPT italic_N italic_N end_POSTSUBSCRIPT end_ARG = 5.44 TeV, p𝑝pitalic_p+Pb collisions at sNNsubscript𝑠𝑁𝑁\sqrt{s_{NN}}square-root start_ARG italic_s start_POSTSUBSCRIPT italic_N italic_N end_POSTSUBSCRIPT end_ARG = 5.02 TeV, and p𝑝pitalic_p+p𝑝pitalic_p collisions at s𝑠\sqrt{s}square-root start_ARG italic_s end_ARG = 5.02 and 13 TeV using Eq. (11).

Figure 4 shows the dependence of fit parameters (V𝑉Vitalic_V, T𝑇Titalic_T, and q𝑞qitalic_q) on charged particle multiplicity (dNch/dηdelimited-⟨⟩𝑑subscript𝑁𝑐𝑑𝜂\langle dN_{ch}/d\eta\rangle⟨ italic_d italic_N start_POSTSUBSCRIPT italic_c italic_h end_POSTSUBSCRIPT / italic_d italic_η ⟩), 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, V𝑉Vitalic_V, increases linearly with increasing dNch/dηdelimited-⟨⟩𝑑subscript𝑁𝑐𝑑𝜂\langle dN_{ch}/d\eta\rangle⟨ italic_d italic_N start_POSTSUBSCRIPT italic_c italic_h end_POSTSUBSCRIPT / italic_d italic_η ⟩ and hence is a good indicator of the system size. However, this increase appears to be slower for p𝑝pitalic_p+p𝑝pitalic_p 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 dNch/dηdelimited-⟨⟩𝑑subscript𝑁𝑐𝑑𝜂\langle dN_{ch}/d\eta\rangle⟨ italic_d italic_N start_POSTSUBSCRIPT italic_c italic_h end_POSTSUBSCRIPT / italic_d italic_η ⟩, the volume of the hadronic system produced in p𝑝pitalic_p+p𝑝pitalic_p collisions is smaller than in Pb+Pb and Xe+Xe collisions. The single measurement in p𝑝pitalic_p+Pb collisions falls within the trend observed in p𝑝pitalic_p+p𝑝pitalic_p collisions at the corresponding dNch/dηdelimited-⟨⟩𝑑subscript𝑁𝑐𝑑𝜂\langle dN_{ch}/d\eta\rangle⟨ italic_d italic_N start_POSTSUBSCRIPT italic_c italic_h end_POSTSUBSCRIPT / italic_d italic_η ⟩, suggesting a potential convergence of behavior in small collision systems.

  • The temperature, T𝑇Titalic_T, increases with increasing dNch/dηdelimited-⟨⟩𝑑subscript𝑁𝑐𝑑𝜂\langle dN_{ch}/d\eta\rangle⟨ italic_d italic_N start_POSTSUBSCRIPT italic_c italic_h end_POSTSUBSCRIPT / italic_d italic_η ⟩. It is observed that T𝑇Titalic_T rises sharply in Pb+Pb and Xe+Xe collisions compared to p𝑝pitalic_p+p𝑝pitalic_p collisions with increasing dNch/dηdelimited-⟨⟩𝑑subscript𝑁𝑐𝑑𝜂\langle dN_{ch}/d\eta\rangle⟨ italic_d italic_N start_POSTSUBSCRIPT italic_c italic_h end_POSTSUBSCRIPT / italic_d italic_η ⟩.

  • The entropy index, q𝑞qitalic_q, for p𝑝pitalic_p+p𝑝pitalic_p collisions increases with increasing dNch/dηdelimited-⟨⟩𝑑subscript𝑁𝑐𝑑𝜂\langle dN_{ch}/d\eta\rangle⟨ italic_d italic_N start_POSTSUBSCRIPT italic_c italic_h end_POSTSUBSCRIPT / italic_d italic_η ⟩, suggesting a deviation from equilibrium in high-multiplicity p𝑝pitalic_p+p𝑝pitalic_p collisions. The q𝑞qitalic_q-parameter for p𝑝pitalic_p+Pb collisions aligns with the trend observed in p𝑝pitalic_p+p𝑝pitalic_p collisions at the corresponding dNch/dηdelimited-⟨⟩𝑑subscript𝑁𝑐𝑑𝜂\langle dN_{ch}/d\eta\rangle⟨ italic_d italic_N start_POSTSUBSCRIPT italic_c italic_h end_POSTSUBSCRIPT / italic_d italic_η ⟩. However, for Pb+Pb and Xe+Xe collisions, q𝑞qitalic_q-parameter decreases with increasing dNch/dηdelimited-⟨⟩𝑑subscript𝑁𝑐𝑑𝜂\langle dN_{ch}/d\eta\rangle⟨ italic_d italic_N start_POSTSUBSCRIPT italic_c italic_h end_POSTSUBSCRIPT / italic_d italic_η ⟩, 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 p𝑝pitalic_p+p𝑝pitalic_p collisions this will cause an increase in the value of q𝑞qitalic_q, contrary to central heavy-ion collisions. We also observe that q𝑞qitalic_q 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, V𝑉Vitalic_V, T𝑇Titalic_T, and q𝑞qitalic_q are used to compute the energy density using the following equation:

ϵ=2i=13gid3p(2π)3Ei(1+(q1)EiT)qq1,italic-ϵ2superscriptsubscript𝑖13subscript𝑔𝑖superscript𝑑3𝑝superscript2𝜋3subscript𝐸𝑖superscript1𝑞1subscript𝐸𝑖𝑇𝑞𝑞1\epsilon=2\sum_{i=1}^{3}g_{i}\int\frac{d^{3}p}{(2\pi)^{3}}E_{i}\left(1+(q-1)% \frac{E_{i}}{T}\right)^{\frac{-q}{q-1}},italic_ϵ = 2 ∑ start_POSTSUBSCRIPT italic_i = 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT italic_g start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ∫ divide start_ARG italic_d start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT italic_p end_ARG start_ARG ( 2 italic_π ) start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT end_ARG italic_E start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ( 1 + ( italic_q - 1 ) divide start_ARG italic_E start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT end_ARG start_ARG italic_T end_ARG ) start_POSTSUPERSCRIPT divide start_ARG - italic_q end_ARG start_ARG italic_q - 1 end_ARG end_POSTSUPERSCRIPT , (12)

The left panel of Fig. 5 shows the energy density, ϵitalic-ϵ\epsilonitalic_ϵ, and the right panel shows the total energy, E𝐸Eitalic_E (= ϵVitalic-ϵ𝑉\epsilon Vitalic_ϵ italic_V) of the hadronic medium at kinetic freeze-out plotted as a function of dNch/dηdelimited-⟨⟩𝑑subscript𝑁𝑐𝑑𝜂\langle dN_{ch}/d\eta\rangle⟨ italic_d italic_N start_POSTSUBSCRIPT italic_c italic_h end_POSTSUBSCRIPT / italic_d italic_η ⟩ in various collision systems.

It is observed that the energy density and total energy increase with increasing dNch/dηdelimited-⟨⟩𝑑subscript𝑁𝑐𝑑𝜂\langle dN_{ch}/d\eta\rangle⟨ italic_d italic_N start_POSTSUBSCRIPT italic_c italic_h end_POSTSUBSCRIPT / italic_d italic_η ⟩. 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 dNch/dηdelimited-⟨⟩𝑑subscript𝑁𝑐𝑑𝜂\langle dN_{ch}/d\eta\rangle⟨ italic_d italic_N start_POSTSUBSCRIPT italic_c italic_h end_POSTSUBSCRIPT / italic_d italic_η ⟩ \approx 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 E𝐸Eitalic_E is plotted against dNch/dηdelimited-⟨⟩𝑑subscript𝑁𝑐𝑑𝜂\langle dN_{ch}/d\eta\rangle⟨ italic_d italic_N start_POSTSUBSCRIPT italic_c italic_h end_POSTSUBSCRIPT / italic_d italic_η ⟩ and a linear scaling between the two is observed. This shows that high multiplicity p𝑝pitalic_p+p𝑝pitalic_p collisions are capable of producing energy densities that are comparable to the most peripheral heavy-ion collisions [35].

Refer to caption
Refer to caption
Figure 5: Energy density (left panel) and total energy (right panel) at kinetic freeze-out as a function of dNch/dηdelimited-⟨⟩𝑑subscript𝑁𝑐𝑑𝜂\langle dN_{ch}/d\eta\rangle⟨ italic_d italic_N start_POSTSUBSCRIPT italic_c italic_h end_POSTSUBSCRIPT / italic_d italic_η ⟩ in Pb+Pb collisions at sNNsubscript𝑠𝑁𝑁\sqrt{s_{NN}}square-root start_ARG italic_s start_POSTSUBSCRIPT italic_N italic_N end_POSTSUBSCRIPT end_ARG = 2.76 and 5.02 TeV, Xe+Xe collisions at sNNsubscript𝑠𝑁𝑁\sqrt{s_{NN}}square-root start_ARG italic_s start_POSTSUBSCRIPT italic_N italic_N end_POSTSUBSCRIPT end_ARG = 5.44 TeV, p𝑝pitalic_p+Pb collisions at sNNsubscript𝑠𝑁𝑁\sqrt{s_{NN}}square-root start_ARG italic_s start_POSTSUBSCRIPT italic_N italic_N end_POSTSUBSCRIPT end_ARG = 5.02 TeV, and p𝑝pitalic_p+p𝑝pitalic_p collisions at s𝑠\sqrt{s}square-root start_ARG italic_s end_ARG = 5.02 and 13 TeV calculated using Eq (12). The expression adNch/dη+b𝑎delimited-⟨⟩𝑑subscript𝑁𝑐𝑑𝜂𝑏a\langle dN_{ch}/d\eta\rangle+bitalic_a ⟨ italic_d italic_N start_POSTSUBSCRIPT italic_c italic_h end_POSTSUBSCRIPT / italic_d italic_η ⟩ + italic_b is fitted to the total energy and we find a=1.463±0.026𝑎plus-or-minus1.4630.026a=1.463\pm 0.026italic_a = 1.463 ± 0.026 GeV and b=1.163±0.298𝑏plus-or-minus1.1630.298b=-1.163\pm 0.298italic_b = - 1.163 ± 0.298 GeV with χ2/ndfsuperscript𝜒2𝑛𝑑𝑓\chi^{2}/ndfitalic_χ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT / italic_n italic_d italic_f = 0.207.

III.2 Pressure

Refer to caption
Figure 6: Pressure at kinetic freeze-out as a function of dNch/dηdelimited-⟨⟩𝑑subscript𝑁𝑐𝑑𝜂\langle dN_{ch}/d\eta\rangle⟨ italic_d italic_N start_POSTSUBSCRIPT italic_c italic_h end_POSTSUBSCRIPT / italic_d italic_η ⟩ in Pb+Pb collisions at sNNsubscript𝑠𝑁𝑁\sqrt{s_{NN}}square-root start_ARG italic_s start_POSTSUBSCRIPT italic_N italic_N end_POSTSUBSCRIPT end_ARG = 2.76 and 5.02 TeV, Xe+Xe collisions at sNNsubscript𝑠𝑁𝑁\sqrt{s_{NN}}square-root start_ARG italic_s start_POSTSUBSCRIPT italic_N italic_N end_POSTSUBSCRIPT end_ARG = 5.44 TeV, p𝑝pitalic_p+Pb collisions at sNNsubscript𝑠𝑁𝑁\sqrt{s_{NN}}square-root start_ARG italic_s start_POSTSUBSCRIPT italic_N italic_N end_POSTSUBSCRIPT end_ARG = 5.02 TeV, and p𝑝pitalic_p+p𝑝pitalic_p collisions at s𝑠\sqrt{s}square-root start_ARG italic_s end_ARG = 5.02 and 13 TeV calculated using Eq (13).
Refer to caption
Figure 7: P(ϵ)𝑃italic-ϵP(\epsilon)italic_P ( italic_ϵ ) equation of state of charged hadrons in Pb+Pb collisions at sNNsubscript𝑠𝑁𝑁\sqrt{s_{NN}}square-root start_ARG italic_s start_POSTSUBSCRIPT italic_N italic_N end_POSTSUBSCRIPT end_ARG = 2.76 and 5.02 TeV, Xe+Xe collisions at sNNsubscript𝑠𝑁𝑁\sqrt{s_{NN}}square-root start_ARG italic_s start_POSTSUBSCRIPT italic_N italic_N end_POSTSUBSCRIPT end_ARG = 5.44 TeV, p𝑝pitalic_p+Pb collisions at sNNsubscript𝑠𝑁𝑁\sqrt{s_{NN}}square-root start_ARG italic_s start_POSTSUBSCRIPT italic_N italic_N end_POSTSUBSCRIPT end_ARG = 5.02 TeV, and p𝑝pitalic_p+p𝑝pitalic_p collisions at s𝑠\sqrt{s}square-root start_ARG italic_s end_ARG = 5.02 and 13 TeV at kinetic freeze-out. The expression cs2ϵsuperscriptsubscript𝑐𝑠2italic-ϵc_{s}^{2}\epsilonitalic_c start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_ϵ is fitted to the data and we find cs2=0.262±0.002superscriptsubscript𝑐𝑠2plus-or-minus0.2620.002c_{s}^{2}=0.262\pm 0.002italic_c start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT = 0.262 ± 0.002 with χ2/ndfsuperscript𝜒2𝑛𝑑𝑓\chi^{2}/ndfitalic_χ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT / italic_n italic_d italic_f = 0.038.

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:

P=2i=13gid3p(2π)3p23Ei(1+(q1)EiT)qq1.𝑃2superscriptsubscript𝑖13subscript𝑔𝑖superscript𝑑3𝑝superscript2𝜋3superscript𝑝23subscript𝐸𝑖superscript1𝑞1subscript𝐸𝑖𝑇𝑞𝑞1P=2\sum_{i=1}^{3}g_{i}\int\frac{d^{3}p}{(2\pi)^{3}}\frac{p^{2}}{3E_{i}}\left(1% +(q-1)\frac{E_{i}}{T}\right)^{\frac{-q}{q-1}}.italic_P = 2 ∑ start_POSTSUBSCRIPT italic_i = 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT italic_g start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ∫ divide start_ARG italic_d start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT italic_p end_ARG start_ARG ( 2 italic_π ) start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT end_ARG divide start_ARG italic_p start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG 3 italic_E start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT end_ARG ( 1 + ( italic_q - 1 ) divide start_ARG italic_E start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT end_ARG start_ARG italic_T end_ARG ) start_POSTSUPERSCRIPT divide start_ARG - italic_q end_ARG start_ARG italic_q - 1 end_ARG end_POSTSUPERSCRIPT . (13)

Figure 6 shows the pressure of the hadronic medium at kinetic freeze-out as a function of dNch/dηdelimited-⟨⟩𝑑subscript𝑁𝑐𝑑𝜂\langle dN_{ch}/d\eta\rangle⟨ italic_d italic_N start_POSTSUBSCRIPT italic_c italic_h end_POSTSUBSCRIPT / italic_d italic_η ⟩, exhibiting an increase with increasing dNch/dηdelimited-⟨⟩𝑑subscript𝑁𝑐𝑑𝜂\langle dN_{ch}/d\eta\rangle⟨ italic_d italic_N start_POSTSUBSCRIPT italic_c italic_h end_POSTSUBSCRIPT / italic_d italic_η ⟩. It is also observed that the pressure rises more rapidly in p𝑝pitalic_p+p𝑝pitalic_p and p𝑝pitalic_p+Pb collisions than in Pb+Pb and Xe+Xe collisions, resulting to a discontinuity at dNch/dηdelimited-⟨⟩𝑑subscript𝑁𝑐𝑑𝜂\langle dN_{ch}/d\eta\rangle⟨ italic_d italic_N start_POSTSUBSCRIPT italic_c italic_h end_POSTSUBSCRIPT / italic_d italic_η ⟩ \approx 30. This can possibly be attributed to the higher initial densities in p𝑝pitalic_p+p𝑝pitalic_p and p𝑝pitalic_p+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, P=cs2ϵ𝑃superscriptsubscript𝑐𝑠2italic-ϵP=c_{s}^{2}\epsilonitalic_P = italic_c start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_ϵ, where cs2superscriptsubscript𝑐𝑠2c_{s}^{2}italic_c start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT 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 cs2superscriptsubscript𝑐𝑠2c_{s}^{2}italic_c start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT = 0.262 ±plus-or-minus\pm± 0.002. A more detailed discussion on cs2superscriptsubscript𝑐𝑠2c_{s}^{2}italic_c start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT is presented later.

III.3 Particle number

Refer to caption
Refer to caption
Figure 8: Number density (left panel) and total number of particles (right panel) at kinetic freeze-out as a function of dNch/dηdelimited-⟨⟩𝑑subscript𝑁𝑐𝑑𝜂\langle dN_{ch}/d\eta\rangle⟨ italic_d italic_N start_POSTSUBSCRIPT italic_c italic_h end_POSTSUBSCRIPT / italic_d italic_η ⟩ in Pb+Pb collisions at sNNsubscript𝑠𝑁𝑁\sqrt{s_{NN}}square-root start_ARG italic_s start_POSTSUBSCRIPT italic_N italic_N end_POSTSUBSCRIPT end_ARG = 2.76 and 5.02 TeV, Xe+Xe collisions at sNNsubscript𝑠𝑁𝑁\sqrt{s_{NN}}square-root start_ARG italic_s start_POSTSUBSCRIPT italic_N italic_N end_POSTSUBSCRIPT end_ARG = 5.44 TeV, p𝑝pitalic_p+Pb collisions at sNNsubscript𝑠𝑁𝑁\sqrt{s_{NN}}square-root start_ARG italic_s start_POSTSUBSCRIPT italic_N italic_N end_POSTSUBSCRIPT end_ARG = 5.02 TeV, and p𝑝pitalic_p+p𝑝pitalic_p collisions at s𝑠\sqrt{s}square-root start_ARG italic_s end_ARG = 5.02 and 13 TeV calculated using Eq (14). The expression adNch/dη𝑎delimited-⟨⟩𝑑subscript𝑁𝑐𝑑𝜂a\langle dN_{ch}/d\eta\rangleitalic_a ⟨ italic_d italic_N start_POSTSUBSCRIPT italic_c italic_h end_POSTSUBSCRIPT / italic_d italic_η ⟩ is fitted to the total number of particles and we find a=1.904±0.030𝑎plus-or-minus1.9040.030a=1.904\pm 0.030italic_a = 1.904 ± 0.030 and b=0.035±0.448𝑏plus-or-minus0.0350.448b=0.035\pm 0.448italic_b = 0.035 ± 0.448 with χ2/ndfsuperscript𝜒2𝑛𝑑𝑓\chi^{2}/ndfitalic_χ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT / italic_n italic_d italic_f = 0.043.

The particle number density can be evaluated from the expression:

n=2i=13gid3p(2π)3(1+(q1)EiT)qq1,𝑛2superscriptsubscript𝑖13subscript𝑔𝑖superscript𝑑3𝑝superscript2𝜋3superscript1𝑞1subscript𝐸𝑖𝑇𝑞𝑞1n=2\sum_{i=1}^{3}g_{i}\int\frac{d^{3}p}{(2\pi)^{3}}\left(1+(q-1)\frac{E_{i}}{T% }\right)^{\frac{-q}{q-1}},italic_n = 2 ∑ start_POSTSUBSCRIPT italic_i = 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT italic_g start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ∫ divide start_ARG italic_d start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT italic_p end_ARG start_ARG ( 2 italic_π ) start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT end_ARG ( 1 + ( italic_q - 1 ) divide start_ARG italic_E start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT end_ARG start_ARG italic_T end_ARG ) start_POSTSUPERSCRIPT divide start_ARG - italic_q end_ARG start_ARG italic_q - 1 end_ARG end_POSTSUPERSCRIPT , (14)

Figure 8 shows increase in the number density, n𝑛nitalic_n, and total number of particles, N𝑁Nitalic_N (= nV𝑛𝑉nVitalic_n italic_V), with increasing dNch/dηdelimited-⟨⟩𝑑subscript𝑁𝑐𝑑𝜂\langle dN_{ch}/d\eta\rangle⟨ italic_d italic_N start_POSTSUBSCRIPT italic_c italic_h end_POSTSUBSCRIPT / italic_d italic_η ⟩. It is observed that the increase in particle density is more pronounced for p𝑝pitalic_p+p𝑝pitalic_p and p𝑝pitalic_p+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 N𝑁Nitalic_N with dNch/dηdelimited-⟨⟩𝑑subscript𝑁𝑐𝑑𝜂\langle dN_{ch}/d\eta\rangle⟨ italic_d italic_N start_POSTSUBSCRIPT italic_c italic_h end_POSTSUBSCRIPT / italic_d italic_η ⟩ 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:

s=2i=13gid3p(2π)3[EiT(1+(q1)EiT)qq1\displaystyle s=2\sum_{i=1}^{3}g_{i}\int\frac{d^{3}p}{(2\pi)^{3}}\left[\frac{E% _{i}}{T}\left(1+(q-1)\frac{E_{i}}{T}\right)^{\frac{-q}{q-1}}\right.italic_s = 2 ∑ start_POSTSUBSCRIPT italic_i = 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT italic_g start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ∫ divide start_ARG italic_d start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT italic_p end_ARG start_ARG ( 2 italic_π ) start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT end_ARG [ divide start_ARG italic_E start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT end_ARG start_ARG italic_T end_ARG ( 1 + ( italic_q - 1 ) divide start_ARG italic_E start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT end_ARG start_ARG italic_T end_ARG ) start_POSTSUPERSCRIPT divide start_ARG - italic_q end_ARG start_ARG italic_q - 1 end_ARG end_POSTSUPERSCRIPT
+(1+(q1)EiT)1q1].\displaystyle\left.+\left(1+(q-1)\frac{E_{i}}{T}\right)^{\frac{-1}{q-1}}\right].+ ( 1 + ( italic_q - 1 ) divide start_ARG italic_E start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT end_ARG start_ARG italic_T end_ARG ) start_POSTSUPERSCRIPT divide start_ARG - 1 end_ARG start_ARG italic_q - 1 end_ARG end_POSTSUPERSCRIPT ] . (15)
Refer to caption
Refer to caption
Figure 9: Entropy density (left panel) and total entropy (right panel) at kinetic freeze-out as a function of dNch/dηdelimited-⟨⟩𝑑subscript𝑁𝑐𝑑𝜂\langle dN_{ch}/d\eta\rangle⟨ italic_d italic_N start_POSTSUBSCRIPT italic_c italic_h end_POSTSUBSCRIPT / italic_d italic_η ⟩ in Pb+Pb collisions at sNNsubscript𝑠𝑁𝑁\sqrt{s_{NN}}square-root start_ARG italic_s start_POSTSUBSCRIPT italic_N italic_N end_POSTSUBSCRIPT end_ARG = 2.76 and 5.02 TeV, Xe+Xe collisions at sNNsubscript𝑠𝑁𝑁\sqrt{s_{NN}}square-root start_ARG italic_s start_POSTSUBSCRIPT italic_N italic_N end_POSTSUBSCRIPT end_ARG = 5.44 TeV, p𝑝pitalic_p+Pb collisions at sNNsubscript𝑠𝑁𝑁\sqrt{s_{NN}}square-root start_ARG italic_s start_POSTSUBSCRIPT italic_N italic_N end_POSTSUBSCRIPT end_ARG = 5.02 TeV, and p𝑝pitalic_p+p𝑝pitalic_p collisions at s𝑠\sqrt{s}square-root start_ARG italic_s end_ARG = 5.02 and 13 TeV calculated using Eq (III.4). The expression adNch/dη+b𝑎delimited-⟨⟩𝑑subscript𝑁𝑐𝑑𝜂𝑏a\langle dN_{ch}/d\eta\rangle+bitalic_a ⟨ italic_d italic_N start_POSTSUBSCRIPT italic_c italic_h end_POSTSUBSCRIPT / italic_d italic_η ⟩ + italic_b is fitted to total entropy and we find a=20.489±0.368𝑎plus-or-minus20.4890.368a=20.489\pm 0.368italic_a = 20.489 ± 0.368 and b=15.012±4.939𝑏plus-or-minus15.0124.939b=15.012\pm 4.939italic_b = 15.012 ± 4.939 with χ2/ndfsuperscript𝜒2𝑛𝑑𝑓\chi^{2}/ndfitalic_χ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT / italic_n italic_d italic_f = 1.702.

In Fig. 9, we observe that the entropy density, s𝑠sitalic_s, and total entropy, S𝑆Sitalic_S (= sV𝑠𝑉sVitalic_s italic_V), increase with increasing dNch/dηdelimited-⟨⟩𝑑subscript𝑁𝑐𝑑𝜂\langle dN_{ch}/d\eta\rangle⟨ italic_d italic_N start_POSTSUBSCRIPT italic_c italic_h end_POSTSUBSCRIPT / italic_d italic_η ⟩. 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 p𝑝pitalic_p+p𝑝pitalic_p and p𝑝pitalic_p+Pb collisions than in Pb+Pb and Xe+Xe collisions, leading to a discontinuity at dNch/dηdelimited-⟨⟩𝑑subscript𝑁𝑐𝑑𝜂\langle dN_{ch}/d\eta\rangle⟨ italic_d italic_N start_POSTSUBSCRIPT italic_c italic_h end_POSTSUBSCRIPT / italic_d italic_η ⟩ \approx 30. It can be inferred from an almost linear relationship between S𝑆Sitalic_S and dNch/dηdelimited-⟨⟩𝑑subscript𝑁𝑐𝑑𝜂\langle dN_{ch}/d\eta\rangle⟨ italic_d italic_N start_POSTSUBSCRIPT italic_c italic_h end_POSTSUBSCRIPT / italic_d italic_η ⟩ 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, ϵ+P=Tsitalic-ϵ𝑃𝑇𝑠\epsilon+P=Tsitalic_ϵ + italic_P = italic_T italic_s, holds.

III.5 Mean free path and Knudsen number

Refer to caption
Refer to caption
Figure 10: Mean free path (left panel) and Knudsen number (right panel) at kinetic freeze-out as a function of dNch/dηdelimited-⟨⟩𝑑subscript𝑁𝑐𝑑𝜂\langle dN_{ch}/d\eta\rangle⟨ italic_d italic_N start_POSTSUBSCRIPT italic_c italic_h end_POSTSUBSCRIPT / italic_d italic_η ⟩ in Pb+Pb collisions at sNNsubscript𝑠𝑁𝑁\sqrt{s_{NN}}square-root start_ARG italic_s start_POSTSUBSCRIPT italic_N italic_N end_POSTSUBSCRIPT end_ARG = 2.76 and 5.02 TeV, Xe+Xe collisions at sNNsubscript𝑠𝑁𝑁\sqrt{s_{NN}}square-root start_ARG italic_s start_POSTSUBSCRIPT italic_N italic_N end_POSTSUBSCRIPT end_ARG = 5.44 TeV, p𝑝pitalic_p+Pb collisions at sNNsubscript𝑠𝑁𝑁\sqrt{s_{NN}}square-root start_ARG italic_s start_POSTSUBSCRIPT italic_N italic_N end_POSTSUBSCRIPT end_ARG = 5.02 TeV, and p𝑝pitalic_p+p𝑝pitalic_p collisions at s𝑠\sqrt{s}square-root start_ARG italic_s end_ARG = 5.02 and 13 TeV. Mean free path is evaluated using Eq (16) and the definition of Knudsen number is given in the text.

The mean free path (λ𝜆\lambdaitalic_λ) 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]:

λ=1nσ,𝜆1𝑛𝜎\lambda=\frac{1}{n\sigma},italic_λ = divide start_ARG 1 end_ARG start_ARG italic_n italic_σ end_ARG , (16)

where n𝑛nitalic_n is the particle density, and σ𝜎\sigmaitalic_σ represents the scattering cross-section. In the specific case of hard core hadron radius approximation, with rh=0.3subscript𝑟0.3r_{h}=0.3italic_r start_POSTSUBSCRIPT italic_h end_POSTSUBSCRIPT = 0.3 fm, σ=4πrh2𝜎4𝜋superscriptsubscript𝑟2\sigma=4\pi r_{h}^{2}italic_σ = 4 italic_π italic_r start_POSTSUBSCRIPT italic_h end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT [10, 39, 40]. Furthermore, Knudsen number (Kn𝐾𝑛Knitalic_K italic_n) is a dimensionless parameter, defined as the ratio of the mean free path (λ𝜆\lambdaitalic_λ) to the radius, R𝑅Ritalic_R (= (3V/4π)1/3superscript3𝑉4𝜋13(3V/4\pi)^{1/3}( 3 italic_V / 4 italic_π ) start_POSTSUPERSCRIPT 1 / 3 end_POSTSUPERSCRIPT), of the hadronic system [14, 15]. A small value of Kn𝐾𝑛Knitalic_K italic_n signifies that the gas system is in the continuum region, allowing the application of hydrodynamics. Generally, when Kn<0.5𝐾𝑛0.5Kn<0.5italic_K italic_n < 0.5, 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 λ𝜆\lambdaitalic_λ and Kn𝐾𝑛Knitalic_K italic_n decrease with increasing dNch/dηdelimited-⟨⟩𝑑subscript𝑁𝑐𝑑𝜂\langle dN_{ch}/d\eta\rangle⟨ italic_d italic_N start_POSTSUBSCRIPT italic_c italic_h end_POSTSUBSCRIPT / italic_d italic_η ⟩. The decrease in λ𝜆\lambdaitalic_λ 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 p𝑝pitalic_p+p𝑝pitalic_p 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 dNch/dηdelimited-⟨⟩𝑑subscript𝑁𝑐𝑑𝜂\langle dN_{ch}/d\eta\rangle⟨ italic_d italic_N start_POSTSUBSCRIPT italic_c italic_h end_POSTSUBSCRIPT / italic_d italic_η ⟩, the Knudsen number (Kn𝐾𝑛Knitalic_K italic_n) for the hadronic gas at KFO was greater than 1. This suggests that the mean free path (λ𝜆\lambdaitalic_λ) is comparable to the system size (R𝑅Ritalic_R), 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 (CVsubscript𝐶𝑉C_{V}italic_C start_POSTSUBSCRIPT italic_V end_POSTSUBSCRIPT) of a system quantifies the amount of energy required to increase its temperature by one unit. Mathematically, this can be expressed as [10, 11]:

CV=ϵT|Vsubscript𝐶𝑉evaluated-atitalic-ϵ𝑇𝑉C_{V}=\left.\frac{\partial\epsilon}{\partial T}\right|_{V}italic_C start_POSTSUBSCRIPT italic_V end_POSTSUBSCRIPT = divide start_ARG ∂ italic_ϵ end_ARG start_ARG ∂ italic_T end_ARG | start_POSTSUBSCRIPT italic_V end_POSTSUBSCRIPT (17)

Figure 11 shows the variation in CVsubscript𝐶𝑉C_{V}italic_C start_POSTSUBSCRIPT italic_V end_POSTSUBSCRIPT as a function of dNch/dηdelimited-⟨⟩𝑑subscript𝑁𝑐𝑑𝜂\langle dN_{ch}/d\eta\rangle⟨ italic_d italic_N start_POSTSUBSCRIPT italic_c italic_h end_POSTSUBSCRIPT / italic_d italic_η ⟩. We observe an increase in CVsubscript𝐶𝑉C_{V}italic_C start_POSTSUBSCRIPT italic_V end_POSTSUBSCRIPT with increasing dNch/dηdelimited-⟨⟩𝑑subscript𝑁𝑐𝑑𝜂\langle dN_{ch}/d\eta\rangle⟨ italic_d italic_N start_POSTSUBSCRIPT italic_c italic_h end_POSTSUBSCRIPT / italic_d italic_η ⟩. This could be a consequence of the increase in the particle number with dNch/dηdelimited-⟨⟩𝑑subscript𝑁𝑐𝑑𝜂\langle dN_{ch}/d\eta\rangle⟨ italic_d italic_N start_POSTSUBSCRIPT italic_c italic_h end_POSTSUBSCRIPT / italic_d italic_η ⟩. 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 CVsubscript𝐶𝑉C_{V}italic_C start_POSTSUBSCRIPT italic_V end_POSTSUBSCRIPT with dNch/dηdelimited-⟨⟩𝑑subscript𝑁𝑐𝑑𝜂\langle dN_{ch}/d\eta\rangle⟨ italic_d italic_N start_POSTSUBSCRIPT italic_c italic_h end_POSTSUBSCRIPT / italic_d italic_η ⟩ in the case of p𝑝pitalic_p+p𝑝pitalic_p collisions, leading to a discontinuity at dNch/dηdelimited-⟨⟩𝑑subscript𝑁𝑐𝑑𝜂\langle dN_{ch}/d\eta\rangle⟨ italic_d italic_N start_POSTSUBSCRIPT italic_c italic_h end_POSTSUBSCRIPT / italic_d italic_η ⟩ \approx 30. This might be due to the higher value of the q𝑞qitalic_q-parameter in p𝑝pitalic_p+p𝑝pitalic_p collisions. Further the system is away from equilibrium, as indicated by a higher q𝑞qitalic_q-value, more is the change in its internal energy as temperature rises, resulting in a higher CVsubscript𝐶𝑉C_{V}italic_C start_POSTSUBSCRIPT italic_V end_POSTSUBSCRIPT [10].

Refer to caption
Figure 11: Heat capacity at kinetic freeze-out as a function of dNch/dηdelimited-⟨⟩𝑑subscript𝑁𝑐𝑑𝜂\langle dN_{ch}/d\eta\rangle⟨ italic_d italic_N start_POSTSUBSCRIPT italic_c italic_h end_POSTSUBSCRIPT / italic_d italic_η ⟩ in Pb+Pb collisions at sNNsubscript𝑠𝑁𝑁\sqrt{s_{NN}}square-root start_ARG italic_s start_POSTSUBSCRIPT italic_N italic_N end_POSTSUBSCRIPT end_ARG = 2.76 and 5.02 TeV, Xe+Xe collisions at sNNsubscript𝑠𝑁𝑁\sqrt{s_{NN}}square-root start_ARG italic_s start_POSTSUBSCRIPT italic_N italic_N end_POSTSUBSCRIPT end_ARG = 5.44 TeV, p𝑝pitalic_p+Pb collisions at sNNsubscript𝑠𝑁𝑁\sqrt{s_{NN}}square-root start_ARG italic_s start_POSTSUBSCRIPT italic_N italic_N end_POSTSUBSCRIPT end_ARG = 5.02 TeV, and p𝑝pitalic_p+p𝑝pitalic_p collisions at s𝑠\sqrt{s}square-root start_ARG italic_s end_ARG = 5.02 and 13 TeV calculated using Eq (17).

III.7 Isothermal compressiblity

Refer to caption
Figure 12: Isothermal compressibility at kinetic freeze-out as a function of dNch/dηdelimited-⟨⟩𝑑subscript𝑁𝑐𝑑𝜂\langle dN_{ch}/d\eta\rangle⟨ italic_d italic_N start_POSTSUBSCRIPT italic_c italic_h end_POSTSUBSCRIPT / italic_d italic_η ⟩ in Pb+Pb collisions at sNNsubscript𝑠𝑁𝑁\sqrt{s_{NN}}square-root start_ARG italic_s start_POSTSUBSCRIPT italic_N italic_N end_POSTSUBSCRIPT end_ARG = 2.76 and 5.02 TeV, Xe+Xe collisions at sNNsubscript𝑠𝑁𝑁\sqrt{s_{NN}}square-root start_ARG italic_s start_POSTSUBSCRIPT italic_N italic_N end_POSTSUBSCRIPT end_ARG = 5.44 TeV, p𝑝pitalic_p+Pb collisions at sNNsubscript𝑠𝑁𝑁\sqrt{s_{NN}}square-root start_ARG italic_s start_POSTSUBSCRIPT italic_N italic_N end_POSTSUBSCRIPT end_ARG = 5.02 TeV, and p𝑝pitalic_p+p𝑝pitalic_p collisions at s𝑠\sqrt{s}square-root start_ARG italic_s end_ARG = 5.02 and 13 TeV calculated using Eq. (18).

The isothermal compressibility (κTsubscript𝜅𝑇\kappa_{T}italic_κ start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT) 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, κTsubscript𝜅𝑇\kappa_{T}italic_κ start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT is 00. It can be evaluated from the expression [10, 12, 13]:

κT=n/μn2.subscript𝜅𝑇𝑛𝜇superscript𝑛2\kappa_{T}=\frac{\partial n/\partial\mu}{n^{2}}.italic_κ start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT = divide start_ARG ∂ italic_n / ∂ italic_μ end_ARG start_ARG italic_n start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG . (18)

Figure 12 shows isothermal compressibility as a function of dNch/dηdelimited-⟨⟩𝑑subscript𝑁𝑐𝑑𝜂\langle dN_{ch}/d\eta\rangle⟨ italic_d italic_N start_POSTSUBSCRIPT italic_c italic_h end_POSTSUBSCRIPT / italic_d italic_η ⟩. We observe that κTsubscript𝜅𝑇\kappa_{T}italic_κ start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT decreases with increasing dNch/dηdelimited-⟨⟩𝑑subscript𝑁𝑐𝑑𝜂\langle dN_{ch}/d\eta\rangle⟨ italic_d italic_N start_POSTSUBSCRIPT italic_c italic_h end_POSTSUBSCRIPT / italic_d italic_η ⟩, 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

Refer to caption
Figure 13: Expansion coefficient at kinetic freeze-out as a function of dNch/dηdelimited-⟨⟩𝑑subscript𝑁𝑐𝑑𝜂\langle dN_{ch}/d\eta\rangle⟨ italic_d italic_N start_POSTSUBSCRIPT italic_c italic_h end_POSTSUBSCRIPT / italic_d italic_η ⟩ in Pb+Pb collisions at sNNsubscript𝑠𝑁𝑁\sqrt{s_{NN}}square-root start_ARG italic_s start_POSTSUBSCRIPT italic_N italic_N end_POSTSUBSCRIPT end_ARG = 2.76 and 5.02 TeV, Xe+Xe collisions at sNNsubscript𝑠𝑁𝑁\sqrt{s_{NN}}square-root start_ARG italic_s start_POSTSUBSCRIPT italic_N italic_N end_POSTSUBSCRIPT end_ARG = 5.44 TeV, p𝑝pitalic_p+Pb collisions at sNNsubscript𝑠𝑁𝑁\sqrt{s_{NN}}square-root start_ARG italic_s start_POSTSUBSCRIPT italic_N italic_N end_POSTSUBSCRIPT end_ARG = 5.02 TeV, and p𝑝pitalic_p+p𝑝pitalic_p collisions at s𝑠\sqrt{s}square-root start_ARG italic_s end_ARG = 5.02 and 13 TeV calculated using Eq. (19).

The expansion coefficient or isobaric expansivity (α𝛼\alphaitalic_α) of a system is the measure of relative expansion per degree change in temperature at a constant pressure. For the hadronic medium at KFO, α𝛼\alphaitalic_α is calculated using the following expression [10]:

α=n/Tn.𝛼𝑛𝑇𝑛\alpha=-\frac{\partial n/\partial T}{n}.italic_α = - divide start_ARG ∂ italic_n / ∂ italic_T end_ARG start_ARG italic_n end_ARG . (19)

Figure 13 shows α𝛼\alphaitalic_α as a function dNch/dηdelimited-⟨⟩𝑑subscript𝑁𝑐𝑑𝜂\langle dN_{ch}/d\eta\rangle⟨ italic_d italic_N start_POSTSUBSCRIPT italic_c italic_h end_POSTSUBSCRIPT / italic_d italic_η ⟩. The value of α𝛼\alphaitalic_α is observed to be negative in the studied collision systems. α𝛼\alphaitalic_α increases with increasing dNch/dηdelimited-⟨⟩𝑑subscript𝑁𝑐𝑑𝜂\langle dN_{ch}/d\eta\rangle⟨ italic_d italic_N start_POSTSUBSCRIPT italic_c italic_h end_POSTSUBSCRIPT / italic_d italic_η ⟩. This is a consequence of the increasing particle density with dNch/dηdelimited-⟨⟩𝑑subscript𝑁𝑐𝑑𝜂\langle dN_{ch}/d\eta\rangle⟨ italic_d italic_N start_POSTSUBSCRIPT italic_c italic_h end_POSTSUBSCRIPT / italic_d italic_η ⟩. 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

Refer to caption
Figure 14: Squared speed of sound at kinetic freeze-out as a function of dNch/dηdelimited-⟨⟩𝑑subscript𝑁𝑐𝑑𝜂\langle dN_{ch}/d\eta\rangle⟨ italic_d italic_N start_POSTSUBSCRIPT italic_c italic_h end_POSTSUBSCRIPT / italic_d italic_η ⟩ in Pb+Pb collisions at sNNsubscript𝑠𝑁𝑁\sqrt{s_{NN}}square-root start_ARG italic_s start_POSTSUBSCRIPT italic_N italic_N end_POSTSUBSCRIPT end_ARG = 2.76 and 5.02 TeV, Xe+Xe collisions at sNNsubscript𝑠𝑁𝑁\sqrt{s_{NN}}square-root start_ARG italic_s start_POSTSUBSCRIPT italic_N italic_N end_POSTSUBSCRIPT end_ARG = 5.44 TeV, p𝑝pitalic_p+Pb collisions at sNNsubscript𝑠𝑁𝑁\sqrt{s_{NN}}square-root start_ARG italic_s start_POSTSUBSCRIPT italic_N italic_N end_POSTSUBSCRIPT end_ARG = 5.02 TeV, and p𝑝pitalic_p+p𝑝pitalic_p collisions at s𝑠\sqrt{s}square-root start_ARG italic_s end_ARG = 5.02 and 13 TeV calculated using Eq. (20). The dotted line represents the value of cs2superscriptsubscript𝑐𝑠2c_{s}^{2}italic_c start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT for a non-interacting massless ideal gas.

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]:

cs2=(Pϵ)s=(P/Tϵ/T).superscriptsubscript𝑐𝑠2subscript𝑃italic-ϵ𝑠𝑃𝑇italic-ϵ𝑇c_{s}^{2}=\left(\frac{\partial P}{\partial\epsilon}\right)_{s}=\left(\frac{% \partial P/\partial T}{\partial\epsilon/\partial T}\right).italic_c start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT = ( divide start_ARG ∂ italic_P end_ARG start_ARG ∂ italic_ϵ end_ARG ) start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT = ( divide start_ARG ∂ italic_P / ∂ italic_T end_ARG start_ARG ∂ italic_ϵ / ∂ italic_T end_ARG ) . (20)

For a non-interacting massless ideal gas, the expected value of cs2superscriptsubscript𝑐𝑠2c_{s}^{2}italic_c start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT is 1/3 times the speed of light squared [44]. Figure 14 shows cs2superscriptsubscript𝑐𝑠2c_{s}^{2}italic_c start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT plotted as a function of dNch/dηdelimited-⟨⟩𝑑subscript𝑁𝑐𝑑𝜂\langle dN_{ch}/d\eta\rangle⟨ italic_d italic_N start_POSTSUBSCRIPT italic_c italic_h end_POSTSUBSCRIPT / italic_d italic_η ⟩. cs2superscriptsubscript𝑐𝑠2c_{s}^{2}italic_c start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT shows a weak dependence with dNch/dηdelimited-⟨⟩𝑑subscript𝑁𝑐𝑑𝜂\langle dN_{ch}/d\eta\rangle⟨ italic_d italic_N start_POSTSUBSCRIPT italic_c italic_h end_POSTSUBSCRIPT / italic_d italic_η ⟩. We observe that for studied collision systems and energies, on an average cs2superscriptsubscript𝑐𝑠2c_{s}^{2}italic_c start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT = 0.266 ±plus-or-minus\pm± 0.001.

IV Summary

The charged hadron spectra measured by the ALICE Collaboration in Pb+Pb collisions at sNNsubscript𝑠𝑁𝑁\sqrt{s_{NN}}square-root start_ARG italic_s start_POSTSUBSCRIPT italic_N italic_N end_POSTSUBSCRIPT end_ARG = 2.76 and 5.02 TeV, Xe+Xe collisions at sNNsubscript𝑠𝑁𝑁\sqrt{s_{NN}}square-root start_ARG italic_s start_POSTSUBSCRIPT italic_N italic_N end_POSTSUBSCRIPT end_ARG = 5.44 TeV, p𝑝pitalic_p+Pb collisions at sNNsubscript𝑠𝑁𝑁\sqrt{s_{NN}}square-root start_ARG italic_s start_POSTSUBSCRIPT italic_N italic_N end_POSTSUBSCRIPT end_ARG = 5.02 TeV, and p𝑝pitalic_p+p𝑝pitalic_p collisions at s𝑠\sqrt{s}square-root start_ARG italic_s end_ARG = 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 p𝑝pitalic_p+p𝑝pitalic_p collisions and a reasonably good description of the spectra in p𝑝pitalic_p+Pb collisions. Furthermore, the dependence of the fit parameters (V𝑉Vitalic_V, T𝑇Titalic_T, and q𝑞qitalic_q) on the dNch/dηdelimited-⟨⟩𝑑subscript𝑁𝑐𝑑𝜂\langle dN_{ch}/d\eta\rangle⟨ italic_d italic_N start_POSTSUBSCRIPT italic_c italic_h end_POSTSUBSCRIPT / italic_d italic_η ⟩ was also noted. It is observed that:

  • Tsallis volume of the hadronic medium increases with increasing dNch/dηdelimited-⟨⟩𝑑subscript𝑁𝑐𝑑𝜂\langle dN_{ch}/d\eta\rangle⟨ italic_d italic_N start_POSTSUBSCRIPT italic_c italic_h end_POSTSUBSCRIPT / italic_d italic_η ⟩, but the growth rate of volume in p𝑝pitalic_p+p𝑝pitalic_p collisions is slower compared to Pb+Pb and Xe+Xe collisions, while p𝑝pitalic_p+Pb collisions show a similar trend to p𝑝pitalic_p+p𝑝pitalic_p collisions. This suggests that small collision systems produce smaller hadronic systems at KFO for the same dNch/dηdelimited-⟨⟩𝑑subscript𝑁𝑐𝑑𝜂\langle dN_{ch}/d\eta\rangle⟨ italic_d italic_N start_POSTSUBSCRIPT italic_c italic_h end_POSTSUBSCRIPT / italic_d italic_η ⟩ compared to large collision systems.

  • Temperature increases with increasing dNch/dηdelimited-⟨⟩𝑑subscript𝑁𝑐𝑑𝜂\langle dN_{ch}/d\eta\rangle⟨ italic_d italic_N start_POSTSUBSCRIPT italic_c italic_h end_POSTSUBSCRIPT / italic_d italic_η ⟩ with a slow increase in small collision systems and a sharp rise in large collision systems.

  • The entropy index, q𝑞qitalic_q, increases with increasing dNch/dηdelimited-⟨⟩𝑑subscript𝑁𝑐𝑑𝜂\langle dN_{ch}/d\eta\rangle⟨ italic_d italic_N start_POSTSUBSCRIPT italic_c italic_h end_POSTSUBSCRIPT / italic_d italic_η ⟩ in p+p𝑝𝑝p+pitalic_p + italic_p collisions, suggesting a deviation from equilibrium in high multiplicity collisions. q𝑞qitalic_q-parameter in p𝑝pitalic_p+Pb collisions is observed to follow the same trend as p+p𝑝𝑝p+pitalic_p + italic_p collisions at the corresponding dNch/dηdelimited-⟨⟩𝑑subscript𝑁𝑐𝑑𝜂\langle dN_{ch}/d\eta\rangle⟨ italic_d italic_N start_POSTSUBSCRIPT italic_c italic_h end_POSTSUBSCRIPT / italic_d italic_η ⟩. However, for Pb+Pb and Xe+Xe collisions, q𝑞qitalic_q decreases with increasing dNch/dηdelimited-⟨⟩𝑑subscript𝑁𝑐𝑑𝜂\langle dN_{ch}/d\eta\rangle⟨ italic_d italic_N start_POSTSUBSCRIPT italic_c italic_h end_POSTSUBSCRIPT / italic_d italic_η ⟩, 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 dNch/dηdelimited-⟨⟩𝑑subscript𝑁𝑐𝑑𝜂\langle dN_{ch}/d\eta\rangle⟨ italic_d italic_N start_POSTSUBSCRIPT italic_c italic_h end_POSTSUBSCRIPT / italic_d italic_η ⟩ are studied. We find that:

  • Energy density increases with increasing dNch/dηdelimited-⟨⟩𝑑subscript𝑁𝑐𝑑𝜂\langle dN_{ch}/d\eta\rangle⟨ italic_d italic_N start_POSTSUBSCRIPT italic_c italic_h end_POSTSUBSCRIPT / italic_d italic_η ⟩, with small collision systems having a higher energy density in the common dNch/dηdelimited-⟨⟩𝑑subscript𝑁𝑐𝑑𝜂\langle dN_{ch}/d\eta\rangle⟨ italic_d italic_N start_POSTSUBSCRIPT italic_c italic_h end_POSTSUBSCRIPT / italic_d italic_η ⟩ region, possibly because of the smaller volumes of the hadronic medium compared to large collision systems.

  • Pressure increases with increasing dNch/dηdelimited-⟨⟩𝑑subscript𝑁𝑐𝑑𝜂\langle dN_{ch}/d\eta\rangle⟨ italic_d italic_N start_POSTSUBSCRIPT italic_c italic_h end_POSTSUBSCRIPT / italic_d italic_η ⟩, 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 dNch/dηdelimited-⟨⟩𝑑subscript𝑁𝑐𝑑𝜂\langle dN_{ch}/d\eta\rangle⟨ italic_d italic_N start_POSTSUBSCRIPT italic_c italic_h end_POSTSUBSCRIPT / italic_d italic_η ⟩, with smaller systems exhibiting a faster increase owing to the smaller volume of the produced hadronic system.

  • Entropy density increases with increasing dNch/dηdelimited-⟨⟩𝑑subscript𝑁𝑐𝑑𝜂\langle dN_{ch}/d\eta\rangle⟨ italic_d italic_N start_POSTSUBSCRIPT italic_c italic_h end_POSTSUBSCRIPT / italic_d italic_η ⟩, with smaller systems having a higher rate of increase.

  • The mean free path decreases with increasing dNch/dηdelimited-⟨⟩𝑑subscript𝑁𝑐𝑑𝜂\langle dN_{ch}/d\eta\rangle⟨ italic_d italic_N start_POSTSUBSCRIPT italic_c italic_h end_POSTSUBSCRIPT / italic_d italic_η ⟩ 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 dNch/dηdelimited-⟨⟩𝑑subscript𝑁𝑐𝑑𝜂\langle dN_{ch}/d\eta\rangle⟨ italic_d italic_N start_POSTSUBSCRIPT italic_c italic_h end_POSTSUBSCRIPT / italic_d italic_η ⟩. This is a consequence of increase in particle number with increasing dNch/dηdelimited-⟨⟩𝑑subscript𝑁𝑐𝑑𝜂\langle dN_{ch}/d\eta\rangle⟨ italic_d italic_N start_POSTSUBSCRIPT italic_c italic_h end_POSTSUBSCRIPT / italic_d italic_η ⟩. 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 dNch/dηdelimited-⟨⟩𝑑subscript𝑁𝑐𝑑𝜂\langle dN_{ch}/d\eta\rangle⟨ italic_d italic_N start_POSTSUBSCRIPT italic_c italic_h end_POSTSUBSCRIPT / italic_d italic_η ⟩ indicating near-ideal behavior with increasing multiplicity.

  • Expansion coefficient is negative for the studied collision systems and it increases with increasing dNch/dηdelimited-⟨⟩𝑑subscript𝑁𝑐𝑑𝜂\langle dN_{ch}/d\eta\rangle⟨ italic_d italic_N start_POSTSUBSCRIPT italic_c italic_h end_POSTSUBSCRIPT / italic_d italic_η ⟩. This is due to the higher particle density associated with greater dNch/dηdelimited-⟨⟩𝑑subscript𝑁𝑐𝑑𝜂\langle dN_{ch}/d\eta\rangle⟨ italic_d italic_N start_POSTSUBSCRIPT italic_c italic_h end_POSTSUBSCRIPT / italic_d italic_η ⟩, 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 dNch/dηdelimited-⟨⟩𝑑subscript𝑁𝑐𝑑𝜂\langle dN_{ch}/d\eta\rangle⟨ italic_d italic_N start_POSTSUBSCRIPT italic_c italic_h end_POSTSUBSCRIPT / italic_d italic_η ⟩. We observe an average value of cs2superscriptsubscript𝑐𝑠2c_{s}^{2}italic_c start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT = 0.266 ±plus-or-minus\pm± 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 p𝑝pitalic_p+p𝑝pitalic_p 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 p𝑝pitalic_p+p𝑝pitalic_p collisions.

Acknowledgements.
CJ acknowledges the financial support from DAE-DST Project No. 3015/I/2021/Gen/RD-I/13283.

References