Next Article in Journal
Quantum Communications Feasibility Tests over a UK-Ireland 224 km Undersea Link
Next Article in Special Issue
Some New Results Involving Past Tsallis Entropy of Order Statistics
Previous Article in Journal
Neural Adaptive H Sliding-Mode Control for Uncertain Nonlinear Systems with Disturbances Using Adaptive Dynamic Programming
Previous Article in Special Issue
On Complex Matrix-Variate Dirichlet Averages and Its Applications in Various Sub-Domains
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Rapidity and Energy Dependencies of Temperatures and Volume Extracted from Identified Charged Hadron Spectra in Proton–Proton Collisions at a Super Proton Synchrotron (SPS)

1
Department of Physics, Xinzhou Normal University, Xinzhou 034000, China
2
State Key Laboratory of Quantum Optics and Quantum Optics Devices, Institute of Theoretical Physics, Shanxi University, Taiyuan 030006, China
3
Laboratory of High Energy Physics, Physical-Technical Institute of Uzbekistan Academy of Sciences, Chingiz Aytmatov Str. 2b, Tashkent 100084, Uzbekistan
4
Department of Natural Sciences, National University of Science and Technology MISIS (NUST MISIS), Almalyk Branch, Almalyk 110105, Uzbekistan
*
Authors to whom correspondence should be addressed.
Entropy 2023, 25(12), 1571; https://doi.org/10.3390/e25121571
Submission received: 11 October 2023 / Revised: 16 November 2023 / Accepted: 20 November 2023 / Published: 22 November 2023

Abstract

:
The standard (Bose–Einstein/Fermi–Dirac, or Maxwell–Boltzmann) distribution from the relativistic ideal gas model is used to study the transverse momentum ( p T ) spectra of identified charged hadrons ( π , π + , K , K + , p ¯ , and p) with different rapidities produced in inelastic proton–proton ( p p ) collisions at a Super Proton Synchrotron (SPS). The experimental data measured using the NA61/SHINE Collaboration at the center-of-mass (c.m.) energies s = 6.3 , 7.7, 8.8, 12.3, and 17.3 GeV are fitted well with the distribution. It is shown that the effective temperature ( T e f f or T), kinetic freeze-out temperature ( T 0 ), and initial temperature ( T i ) decrease with the increase in rapidity and increase with the increase in c.m. energy. The kinetic freeze-out volume (V) extracted from the π , π + , K , K + , and p ¯ spectra decreases with the rapidity and increase with the c.m. energy. The opposite tendency of V, extracted from the p spectra, is observed to be increasing with the rapidity and decreasing with the c.m. energy due to the effect of leading protons.

1. Introduction

The existence of confinement and asymptotic freedom in Quantum Chromodynamics (QCD) has led to many conjectures about the thermodynamic and transport properties of hot and dense matter. Because of confinement, nuclear matter should be composed of low-energy hadrons, and it is considered a weakly interacting gas of hadrons. On the other hand, at very high energies, asymptotic freedom means that the interactions between quarks and gluons are very weak, and the nuclear matter is considered as a weakly coupling gas of quarks and gluons. There should be a phase transition between these two configurations, in which the degrees of freedom of hadrons disappear and Quark–Gluon Plasma (QGP) is formed, which is generated at a sufficiently high temperature or density [1,2,3,4,5,6]. QGP existed in the very early universe (a few microseconds after the Big Bang), and some forms of this matter may still exist in the core of neutron stars. Ultra-relativistic heavy-ion collisions have provided opportunities to systematically create and study different phases of bulk nuclear matter.
Several experiments performed at the Super Proton Synchrotron (SPS) [7,8], Relativistic Heavy Ion Collider (RHIC) [2,3,9,10,11,12,13,14,15], and Large Hadron Collider (LHC) [16,17,18,19] have reported abundant experimental data. The system of proton–proton ( p p ) collisions is usually used as a reference measurement for heavy ion collisions, as it has several valence quarks involved in the collisions. Collective flow is one of the characteristics of the thermal dense medium of this strongly interacting matter. The generated medium expands collectively such that the flow effect is expected to be distinguished from the thermal motion, which reflects the temperature. The heavy ion physics community has been fascinated by observing unexpected collective behavior in high-multiplicity p p collision events. It is therefore necessary and important to study p p collisions.
The transverse momentum ( p T ) spectra of identified charged hadrons produced in relativistic or high-energy collisions contain abundant information on the collision dynamics and the evolution properties of the system from the initial stage to the end of freeze-out phase [20]. Traditionally, it is believed that the flattening of the p T spectra with high multiplicity is a signal for the formation of a mixed phase of de-confined partons and hadrons. In the hydrodynamical model, the slopes of p T spectra are co-determined by the kinetic freeze-out temperature and the transverse expansion flow of the collision system [21]. The study of p T spectra can reveal information related to the effective temperature ( T e f f or T) of the system. A plateau-like region observed in the excitation function of T is considered a possible signal for the formation of mixed-phases, similar to the temperature dependence of entropy observed in the first-order phase transition. In addition, in order to understand the phase transition from QGP to hadronic matter, the transverse momentum density is often studied.
In the physical process of high-energy heavy ion collisions, at least four temperatures are often used, namely initial temperature ( T i ), chemical freeze-out temperature ( T c h ), kinetic (or thermal) freeze-out temperature ( T 0 ), and T. These temperatures correspond to different stages of collisions. The excitation degree of the interaction system at the initial stage is described by T i , at which hadrons undergo elastic and inelastic interactions in the hadronic medium. Due to the shortage of research methods, there is limited research on T i in the community, which should be based on the p T . With the decrease in temperature, the system begins to form hadronic matter and enters the chemical freeze-out stage. Under the condition of maintaining a certain degree of local dynamic equilibrium through quasi-elastic resonance scattering, the final stable hadronic yield has almost no change [22,23,24,25]. The T c h , and baryon chemical potential ( μ B ) at this stage can be obtained by using various thermodynamic models [3,26,27,28]. After the chemical freeze-out stage, the system further expands as the interactions become weak. Finally, the system enters the kinetic freeze-out stage as the elastic collisions between hadrons disappear.
In this paper, the p T spectra of identified charged hadrons ( π , π + , K , K + , p ¯ , and p) with different rapidities produced in inelastic p p collisions at the center-of-mass (c.m.) energies s = 6.3 , 7.7, 8.8, 12.3, and 17.3 GeV at the SPS [29] are studied, where the c.m. energy is also referred to as collision energy. Although the nonextensive distribution of the Tsallis statistics [30,31,32,33,34,35] has been widely used in recent years, the standard (Bose–Einstein/Fermi–Dirac, or Maxwell–Boltzmann) distribution from the relativistic ideal gas model is still used to extract T directly and then to obtain the average transverse momentum ( p T ), root-mean-square transverse momentum ( p T 2 ), T 0 , and T i indirectly.
The remainder of this paper is structured as follows. The formalism and method are described in Section 2. Results and discussion are provided in Section 3. In Section 4, we summarize our main observations and conclusions.

2. Formalism and Method

The particles produced in inelastic p p collisions are thought to be controlled by two main mechanisms or excitation degrees. The low- p T region, which is less than 1–2 GeV/c is dominated by the soft excitation process [36,37]. The high- p T region that is more than 1–2 GeV/c is governed by the hard scattering process [36,37]. The soft process corresponds to a low excitation degree, and the hard process implies a high excitation degree. The two-mechanism scheme is only one possible choice in understanding particle production. If the particles are distributed in a very wide p T region, one should consider the multiple mechanisms or excitation degrees. If the particles are distributed in a relatively narrow p T region, one may choose the single mechanism or excitation degree. In the two-mechanism scenario, it is currently believed that most light-flavor particles are produced in the soft process. The spectrum in the low- p T region shows exponential behavior, which can be fitted by the thermal distribution [38,39,40]. Heavy-flavor particles and some light-flavor particles are produced in the hard process. The spectrum in high- p T region shows inverse power-law behavior and can be fitted using the Hagedorn [41,42], Tsallis–Levy [31,32], or Tsallis–Pareto-type functions [32,33,34,35].
In this investigation, the light particle spectra in the low- p T region in inelastic p p collisions at the SPS are studied by using the most basic thermal distribution, the standard distribution, which comes from the relativistic ideal gas model. The invariant particle momentum (p) distribution described by the standard distribution can be given by [30]
E d 3 N d 3 p = 1 2 π p T d 2 N d y d p T = g V ( 2 π ) 3 E exp E μ T + S 1 ,
where N is the particle number, g is the degeneracy factor, V is the volume, μ is the chemical potential,
E = p 2 + m 0 2 = m T cosh y
is the energy,
m T = p T 2 + m 0 2
is the transverse mass,
y = 1 2 ln 1 + β z 1 β z = tanh 1 ( β z )
is the rapidity, β z is the longitudinal velocity, and S = 1 , 1, and 0 correspond to the Bose–Einstein, Fermi–Dirac, and Maxwell–Boltzmann statistics, respectively.
For the wide p T spectra, if a multi-component standard distribution
E d 3 N d 3 p = 1 2 π p T d 2 N d y d p T = i = 1 n g V i ( 2 π ) 3 E exp E μ T i + S 1
can be used in the fit, one may obtain multiple temperatures, that is, the temperature fluctuation. Here, n denotes the number of components. Let k i ( i = 1 , 2, …, n) denote the relative fraction of the i-th component, and V i and T i are the volume and temperature corresponding to the i-th component, respectively. Naturally, one has
V = i = 1 n V i , T = i = 1 n k i T i , i = 1 n k i = 1 .
Here, k i = V i / V .
Because of the temperature fluctuation, there are interactions among different subsystems or local sources due to the exchange of heat energy. This causes the couplings of entropy functions of various subsystems. The total entropy is then the sum of the entropies of subsystems plus the entropies of the couplings. The temperature fluctuation in the multi-component standard distribution is a way to explain the origin of the Tsallis distribution. Generally, the p T spectra, which can be fitted using the multi-component standard distribution, can also be fitted using the Tsallis distribution. Because of the influence of the entropy index (q), the temperature value extracted from the Tsallis distribution is smaller than that from the multi-component standard distribution. In fact, in the fit using the Tsallis distribution, increasing T and/or q can increase the particle yield in the high- p T region conveniently.
The data sample analyzed in the present work is in the low- p T region. This implies that the standard distribution can be used. In the standard distribution, the unit-density function of y and p T is written as
d 2 N d y d p T = g V ( 2 π ) 2 p T m T cosh y × exp m T cosh y μ T + S 1 .
Then, the density function of p T is
d N d p T = g V ( 2 π ) 2 p T m T y min y max cosh y × exp m T cosh y μ T + S 1 d y ,
where y min and y max are the minimum and maximum rapidities in the rapidity interval, respectively. The density function of y is
d N d y = g V ( 2 π ) 2 cosh y 0 p T max p T m T × exp m T cosh y μ T + S 1 d p T ,
where p T max is the maximum p T in the considered rapidity interval. Although p T max can be mathematically infinite, it is only large enough in physics due to the limitations of the conservation of energy and momentum.
No matter what the specific form of particle momentum distribution is used, the probability density function of p T is written in general as
f ( p T ) = 1 N d N d p T .
Naturally, f ( p T ) is normalized to 1. That is,
0 f ( p T ) d p T = 1 .
One has the average transverse momentum,
p T = 0 p T f ( p T ) d p T 0 f ( p T ) d p T = 0 p T f ( p T ) d p T ,
and the root-mean-square p T ,
p T 2 = 0 p T 2 f ( p T ) d p T 0 f ( p T ) d p T = 0 p T 2 f ( p T ) d p T .
In principle, there are three independent chemical potentials, baryon ( μ B ), electric charge or isospin ( μ I ), and strangeness ( μ S ), which are related to the three conserved charges. Although the chemical potential, μ π ( μ K or μ p ), of the pion (kaon or proton) can be written in terms of the above three chemical potentials [43,44,45,46,47,48,49], we obtained them by using an alternative method in the present work for more convenience.
Considering the yield ratio [ k j ( j = π , K, and p)] of negatively to positively charged hadrons ( j to j + ), the corresponding chemical potentials ( μ j and μ j + ), and the corresponding source temperature ( T j and T j + ), one has that the relationship between k j and μ j is [20,50,51,52,53]
k j j j + = exp μ j T j μ j + T j + = exp 2 μ j T j
if the conditions
T j = T j + = T j , μ j = μ j + = μ j
are satisfied. Here, j and j + also denote the yields of negative and positive hadrons respectively. k j can be obtained simply from the experimental data, and T j should be the chemical kinetic-freezing temperature T c h , which is slightly larger than or equal to the effective temperature T due to the short lifetime of the system formed in p p collisions. One has T j T in this work.
Further, one has
μ j = 1 2 T j ln k j .
Obviously, μ j is energy-dependent due to T j and k j being energy-dependent. Based on a collection of large amounts of experimental data, our previous work [52,53] presents the excitation functions of μ j in p p and central heavy ion collisions, which can be used for a direct extraction for this study. In particular, μ j decreases quickly with the increase in energy in p p collisions in the concerned SPS energy range. However, the tendency of μ π in central heavy ion collisions is opposite to that in p p collisions, though the tendency of μ K is similar, and that of μ p is also similar in the two collisions. The three μ j in both the collisions are close to 0 at around 100 GeV and above.
The chemical freeze-out temperature T c h in central heavy ion collisions is also energy-dependent [43,44,45,46,47,48,49], which shows a tendency for a rapid increase at a few GeV and then saturation at dozens of GeV and above. In view of the fact that the tendency of T c h has a parameterized excitation function with unanimity in the community, the present work does not study T c h parameter.
Generally, the kinetic freeze-out temperature T 0 has a tendency of a rapid increase at a few GeV, and then an ambiguous tendency (increase, decrease, or saturation) appears at dozens of GeV and above. It is worth studying the tendency of T 0 further. A thermal-related method shows that [54]
T 0 = p T 2 κ 0 ,
where κ 0 = 3.07 is a coefficient, and a value 2 is introduced by us because two participant partons (one from the projectile and the other from the target) are assumed to contribute to p T . This formula gives an approximate consistent tendency of T 0 as another thermal-related method [55], which shows T 0 to be proportional to p T and the coefficient to be energy-related, though the results from the two methods are not the same.
The initial temperature T i , which is comparable to the experimental data, is less studied in the community. According to the string percolation model [56,57,58], T i is expressed as
T i = p T 2 2 F ( ξ ) ,
where
F ( ξ ) = 1 exp ( ξ ) ξ
is the color-suppression factor related to the dimensionless percolation density parameter ξ . In p p collisions, F ( ξ ) 1 due to the low string overlap probability. As an initial quantity, T i should reflect the excitation degree of the system at the parton level. Correspondingly, the final quantity T 0 should also be extracted at the parton level. This is also the reason that the value of 2 is introduced by us in the denominator of the T 0 expression if one assumes that two participant partons are the energy sources in the formation of a particle.
The kinetic energy of a particle’s directional movement should not be reflected in the temperature parameters. The experimental data used in this paper were all measured in the forward-rapidity region. In order to remove the influence of directional motion, one can directly shift the forward rapidity and its interval to the mid-rapidity with the same interval width during the fitting process. In this paper, we integrate y from y min = 0.1 to y max = 0.1 in the fit to give a more accurate result, though y 0 and cosh y 1 near the mid-rapidity. The small difference (<1%) between the accurate and approximate calculations appears mainly in the normalization but not in the temperature parameter.
The method of least squares based on obtaining the minimum χ 2 is adopted to obtain the best parameters and their uncertainties. The treatment method is given in Appendix A.

3. Results and Discussion

Figure 1 and Figure 2 show the rapidity-dependent double differential p T spectra, d 2 N / d y d p T , of π and π + respectively, produced in inelastic p p collisions at the SPS. Panels (a)–(e) correspond to the results of s = 6.3 , 7.7, 8.8, 12.3, and 17.3 GeV, respectively. The symbols represent the experimental data at different y, with an interval width of 0.2 units, measured using the NA61/SHINE Collaboration [29], and the curves are our results fitted from the Bose–Einstein distribution. In order to see the fitting effect more clearly, the experimental data and fitting results at different rapidities are multiplied by different factors labeled in the panel for scaling. The values of related free parameters (T), the normalization constant (V), χ 2 , and the number of degrees of freedom (ndof) for the curves in Figure 1 and Figure 2 are listed in Table A1 in Appendix B. One can see that the fitting results with the Bose–Einstein distribution are in good agreement with the experimental data of π and π + spectra, measured using the NA61/SHINE Collaboration in p p collisions at different s and in different y intervals.
Similarly, Figure 1, Figure 2, Figure 3 and Figure 4 show the rapidity-dependent d 2 N / d y d p T of K and K + , respectively, produced in inelastic p p collisions at different s . The values of T, V, and χ 2 /ndof for the curves in Figure 3 and Figure 4 are listed in Table A2 in Appendix B. One can see that the fitting results from the Bose–Einstein distribution are in agreement with the experimental data of K and K + , measured by the NA61/SHINE Collaboration in p p collisions at different s and in different y intervals.
Similar to Figure 1, Figure 2, Figure 3 and Figure 4, Figure 5 and Figure 6 show the rapidity-dependent d 2 N / d y d p T of p ¯ and p, respectively, produced in inelastic p p collisions at different s . The experimental data of p ¯ at s = 6.3 GeV in Figure 5 are not available. The values of T, V, and χ 2 /ndof for the curves in Figure 5 and Figure 6 are listed in Table A3 in Appendix B. One can see that the p T spectra of p ¯ and p in p p collisions are shown to obey approximately the Fermi–Dirac distribution.
To show more intuitively the dependence of the free parameter T and derived quantities (the kinetic freeze-out temperature T 0 and initial temperature T i ) on rapidity, y, and c.m. energy, s , Figure 7, Figure 8, Figure 9 and Figure 10 show the relations of Ty, T 0 y, T i y, and Vy at different s , respectively, and Figure 11, Figure 12, Figure 13 and Figure 14 show the relations of T s , T 0 s , T i s , and V s at different y, respectively. Panels (a)–(f) correspond to the results from π , π + , K , K + , p ¯ , and p spectra, respectively. These figures show some changing trends of parameters.
In most cases, one can generally see that T, T 0 , and T i decrease (increase) with the increase in y ( s ). There is a tendency of saturation for the three temperatures at s = 7.7 GeV and above. Being the initial energy of a saturation effect, 7.7 GeV is a special energy at which the reaction products are proton-dominated and above which the products are meson-dominated. For π , π + , K , K + , and p ¯ spectra, the extracted V also decreases (increases) with the increase in y ( s ). However, for p spectra, the extracted V shows an opposite tendency, increasing (decreasing) with the increase in y ( s ).
There is an isospin and mass independence of T. This property is exactly that of T c h , which implies a single scenario of chemical freeze-out. However, although T 0 and T i are isospin-independent, they increase with the increase in mass. The mass dependence of T 0 is a reflection of a mass-dependent differential kinetic freeze-out scenario or multiple kinetic freeze-out scenarios. The mass dependence of T i means that the formation moments of different particles are different. With the increase in T 0 ( T i ), massive particles are emitted (formed) earlier. On average, this work shows that p ¯ ( p ) are emitted (formed) earlier than K , and K are emitted (formed) earlier than π , though the relaxation times for the emissions (formations) of different particles can overlap.
Except for V from the p spectra, the tendencies of other parameters from the p spectra, and the tendencies of parameters from the spectra of other particles are easy to understand. It is expected that the local system in the mid-rapidity region has more deposited energy than that in the forward region. Meanwhile, the collision system at a higher energy has more deposited energy than that at lower energy. This results in a higher excitation degree (then higher temperature) at the mid-rapidity and more produced particles (then larger volume) at a higher energy.
The V tendency from the p spectra is opposite to that from the spectra of other particles. The reason is that the pre-existing leading protons affect the p spectra. Because of the leading protons appearing in the forward region, the number of protons and then the volume of a proton source in the fixed interval are small at the mid-rapidity. At a higher energy, the leading protons appear in the more forward region, which leads to a smaller V in the fixed interval in the rapidity space. In the present work, the fixed interval is that Δ y = y max y min = 0.2 .
The values of V depend on particle mass and charge. Excluding the case of p, which contains pre-existing leading protons in the p p system, V decreases significantly with the increase in mass, and positive hadrons correspond to the larger V of the emission source. This is because the larger the mass, the more difficult it is to produce this particle. Meanwhile, there is an electromagnetic exclusion (attraction) between positive (negative) hadrons and pre-existing protons. This causes larger (smaller) V with an emission source of positive (negative) hadrons.
Generally, the effective temperature T is proportional to the mean transverse momentum p T . The present work shows that T π 0.351 p T π , T π + 0.348 p T π + , T K 0.284 p T K , T K + 0.293 p T K + , T p ¯ 0.234 p T p ¯ , and T p 0.240 p T p . Here, the type of a particle appears as the subscript label of the related quantity. The ratio of T / p T is approximately independent of a particle mass. This is consistent with the ratios of T 0 / p T and T i / p T 2 , which are independent of particle mass according to Equations (17) and (18).
As only a free parameter, T does not show an obvious dependence on particle type or mass. However, it is hard to extract exact information from T because it is not a real temperature, because it also contains the contribution of transverse flow. T 0 is smaller than T i due to the fact that T 0 is “measured” at the kinetic freeze-out stage (the final one), and T i is “measured” at the initial stage. From the initial stage to the final one, the system becomes colder and colder. This is indeed observed in the present work.
In the above discussions, although chemical potential μ runs through the entire process, it is an insensitive quantity in the fit and not a free parameter due to the fact that it depends on T c h and k j . Our previous work [52,53] shows that, from 6.3 to 17.3 GeV, μ π + , μ K + , and μ p are around 0.041–0.017, 0.110–0.042, and 0.510–0.180 GeV, respectively, which are directly used in this work. These results have excluded the contributions from resonance decays [59]. Although the resonance decays contribute considerably to the yields of negative and positive hadrons, they contribute to the yield ratios, and then, the chemical potentials are small [52,53].
Before the summary and conclusions, it should be pointed out that the data sets analyzed by us are in a narrow and low- p T range and obey the standard distribution. We believe that even if the narrow spectra are in a high- p T range, the standard distribution can be used, and a high temperature can be obtained. The success of this work reflects that the classical concept and distribution can still play a great role in the field of high-energy collisions, though the application is in a local region. In our opinion, when researchers search for novel theoretical models, they first need to take into account classical theories.
Although the topic has been extensively studied in many papers for the SPS, RHIC, and LHC heavy-ion collisions and outline the validity of a nonextensive statistical distribution [60,61,62,63,64,65,66], those investigations used the spectra in a wide p T range. It is unanimous that for the wide p T spectra, a two-, three-, or multi-component standard distribution is needed in the fit. Then, a temperature fluctuation can be observed from the multi-component standard distribution. At this point, the Tsallis distribution is needed. This is the relationship between the standard distribution and the Tsallis distribution in the fit process.
In addition, in comparison with Hanbury–Brown–Twiss (HBT) results [67], large values of volume are obtained in the present work. The reason is that different volumes are studied. Generally, the former describes the system size in the initial state of collisions, and the latter is a reflection of the size of an expanded fireball in the final state (at the kinetic freeze-out) of collisions. Obviously, the latter is much larger than the former. The values of the three temperatures obtained in the present work seem reasonable.

4. Summary and Conclusions

The main observations and conclusions are summarized here.
(a) The transverse momentum spectra of the identified charged hadrons ( π , π + , K , K + , p ¯ , and p) with different rapidities produced in proton–proton collisions at center-of-mass energies s = 6.3 , 7.7, 8.8, 12.3, and 17.3 GeV have been studied using the standard distribution. The fitted results are in agreement with the experimental data measured by the NA61/SHINE Collaboration at the SPS. The effective temperature T, kinetic freeze-out temperature T 0 , initial temperature T i , and kinetic freeze-out volume V are extracted. The present work shows that the standard distribution coming from the relativistic ideal gas model works well in some cases.
(b) In most cases, T, T 0 , and T i decrease with the increase in rapidity y and increase with the increase in s . There is a tendency of saturation for the three temperatures at s = 7.7 GeV and above. From a quick increase to a slow saturation in the three temperatures, the transition energy 7.7 GeV is the boundary for proton-dominated and meson-dominated final states. For the spectra of produced hadrons ( π , π + , K , K + , and p ¯ ), the extracted V also decreases with the increase in y and increases with the increase in s . For the spectra of p, the extracted V increases with the increase in y and decreases with the increase in s . This is opposite to other hadrons because p contains the pre-existing leading protons, which affect the result.
(c) The three temperatures do not show an obvious isospin dependence. However, V shows a significant isospin dependence. The reason for the isospin dependence of V is the electromagnetic interactions between positive (negative) hadrons and pre-existing protons. The exclusion (attraction) between positive (negative) hadrons and pre-existing protons causes larger (smaller) V of an emission source of positive (negative) hadrons. Compared with the three temperature types, V shows a larger mass dependence. The mass dependence of V is also a reflection of a mass-dependent differential kinetic freeze-out scenario or multiple kinetic freeze-out scenario.

Author Contributions

Conceptualization, F.-H.L. and K.K.O.; Methodology, F.-H.L. and K.K.O.; Software, P.-P.Y.; Validation, F.-H.L. and K.K.O.; Formal analysis, P.-P.Y.; Investigation, P.-P.Y.; Resources, P.-P.Y.; Data curation, P.-P.Y.; Writing—original draft, P.-P.Y.; Writing—review and editing, F.-H.L. and K.K.O.; Visualization, P.-P.Y.; Supervision, F.-H.L. and K.K.O.; Project administration, P.-P.Y. and F.-H.L.; Funding acquisition, P.-P.Y., F.-H.L. and K.K.O. All authors have read and agreed to the published version of the manuscript.

Funding

The work of P.-P.Y. was supported by the Shanxi Provincial Natural Science Foundation under Grant No. 202203021222308 and the Doctoral Scientific Research Foundations of Shanxi Province and Xinzhou Normal University. The work of F.-H.L. was supported by the National Natural Science Foundation of China under Grant No. 12147215, the Shanxi Provincial Natural Science Foundation under Grant No. 202103021224036, and the Fund for Shanxi “1331 Project” Key Subjects Construction. The work of K.K.O. was supported by the Agency of Innovative Development under the Ministry of Higher Education, Science and Innovations of the Republic of Uzbekistan within the fundamental project No. F3-20200929146 on the analysis of open data on heavy-ion collisions at RHIC and LHC.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

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

Conflicts of Interest

The authors declare no conflict of interest. The funding agencies have no role in the design of the study; in the collection, analysis, or interpretation of the data; in the writing of the manuscript; or in the decision to publish the results.

Appendix A. The Method to Obtain the Parameter and Its Uncertainty

Let
y i = f ( x i ) , i = 1 , 2 , , n
be the model value of the i-th fitting point. The physical quantities or parameters, λ ( λ 1 , λ 2 , , λ j ), can be obtained by fitting the experimental data, where j is the number of parameters, which includes the normalization constant. One has
χ 2 = i = 1 n [ f ( x i ) Y i ] 2 ( δ Y i ) 2 ,
where n is the number of fitting points, Y i represents the experimental value, and δ Y i represents the uncertainty of the experimental value, usually including statistical and systematic uncertainties.
Due to the small particle number of p T samples being studied in this paper, the parameter uncertainty is assumed to follow the Student’s distribution (shortened to the t-distribution) [68],
f ( λ ν ) = ( ν ν + λ 2 ) ν + 1 2 Γ ( ν + 1 2 ) μ π Γ ( ν 2 ) ,
where ν represents the ndof and Γ ( x ) represents the Gamma function. With the increase in ν , the t-distribution gradually approaches the normal or Gaussian distribution N ( 0 , 1 ) . When ν approaches 1, the t-distribution approaches the Cauchy distribution.
In the present fitting, a 0.5 % confidence interval is used to describe the parameter uncertainty. This means that there is a 0.5 % probability that the parameter will fall within ( λ j t σ , λ j + t σ ), where σ is the standard deviation of each parameter and t satisfies the equation
P = t σ t σ f ( λ ) d λ = 0.005 .
The standard deviation of each parameter can be calculated using
σ = ( J λ T J λ ) 1 s 2 ν ,
where
( J λ ) i j = [ f ( x i ) Y i ] λ j
is the Jacobian matrix and determinant of the model, J λ T is the transpose of the Jacobian matrix, the superscript 1 represents matrix inversion, and
s 2 = i = 1 n [ f ( x i ) Y i ] 2
is the variance. Then, the corresponding best parameter is given by
λ j [ λ j t σ j , j , λ j + t σ j , j ]
with an uncertainty of t σ j , j .

Appendix B. Parameter Tables Obtained in the Fitting Process

Table A1. Values of T, V, χ 2 , and ndof corresponding to the curves in Figure 1 and Figure 2 for π and π + produced in p p collisions, respectively, where the values of s are given together. The values of χ 2 are reserved as an integer or the first non-zero decimal (if less than 1). The symbol “−” indicates that the data are not available.
Table A1. Values of T, V, χ 2 , and ndof corresponding to the curves in Figure 1 and Figure 2 for π and π + produced in p p collisions, respectively, where the values of s are given together. The values of χ 2 are reserved as an integer or the first non-zero decimal (if less than 1). The symbol “−” indicates that the data are not available.
s (GeV)y π π +
T (GeV)V ( fm 3 ) χ 2 /ndofT (GeV)V ( fm 3 ) χ 2 /ndof
  0.1 0.118 ± 0.001 ( 8.119 ± 0.177 ) × 10 1 13 / 16
  0.3 0.123 ± 0.001 ( 7.755 ± 0.167 ) × 10 1 17 / 16 0.115 ± 0.002 ( 1.268 ± 0.072 ) × 10 2 0.09 / 6
  0.5 0.119 ± 0.001 ( 7.247 ± 0.157 ) × 10 1 15 / 16 0.108 ± 0.001 ( 9.805 ± 0.269 ) × 10 1 2 / 10
  0.7 0.118 ± 0.001 ( 6.514 ± 0.141 ) × 10 1 11 / 16 0.107 ± 0.001 ( 1.017 ± 0.028 ) × 10 2 1 / 10
  0.9 0.114 ± 0.001 ( 5.516 ± 0.121 ) × 10 1 11 / 16 0.098 ± 0.001 ( 1.007 ± 0.029 ) × 10 2 3 / 10
6.3 1.1 0.108 ± 0.001 ( 4.464 ± 0.098 ) × 10 1 16 / 16 0.100 ± 0.001 ( 7.106 ± 0.179 ) × 10 1 2 / 9
  1.3 0.097 ± 0.001 ( 3.466 ± 0.081 ) × 10 1 14 / 15 0.087 ± 0.001 ( 8.072 ± 0.468 ) × 10 1 1 / 6
  1.5 0.088 ± 0.001 ( 2.472 ± 0.061 ) × 10 1 9 / 14 0.088 ± 0.001 ( 5.903 ± 0.423 ) × 10 1 1 / 4
  1.7 0.081 ± 0.001 ( 3.957 ± 0.196 ) × 10 1 4 / 5
  1.9 0.070 ± 0.001 ( 2.736 ± 0.112 ) × 10 1 4 / 6
  2.1 0.066 ± 0.001 ( 2.036 ± 0.099 ) × 10 1 0.1 / 5
  0.1 0.123 ± 0.001 ( 9.450 ± 0.204 ) × 10 1 43 / 16 0.130 ± 0.003 ( 1.261 ± 0.137 ) × 10 2 0.5 / 3
  0.3 0.124 ± 0.001 ( 9.141 ± 0.197 ) × 10 1 33 / 16 0.121 ± 0.001 ( 1.423 ± 0.052 ) × 10 2 0.5 / 6
  0.5 0.122 ± 0.001 ( 8.475 ± 0.183 ) × 10 1 32 / 16 0.121 ± 0.001 ( 1.356 ± 0.037 ) × 10 2 2 / 9
  0.7 0.122 ± 0.001 ( 7.944 ± 0.172 ) × 10 1 19 / 16 0.119 ± 0.001 ( 1.253 ± 0.035 ) × 10 2 1 / 9
  0.9 0.119 ± 0.001 ( 6.865 ± 0.149 ) × 10 1 20 / 16 0.116 ± 0.001 ( 1.165 ± 0.032 ) × 10 2 2 / 9
7.7 1.1 0.115 ± 0.001 ( 5.715 ± 0.124 ) × 10 1 17 / 16 0.116 ± 0.001 ( 9.373 ± 0.251 ) × 10 1 0.5 / 8
  1.3 0.107 ± 0.001 ( 4.585 ± 0.104 ) × 10 1 14 / 15 0.108 ± 0.001 ( 8.131 ± 0.354 ) × 10 1 1 / 7
  1.5 0.101 ± 0.001 ( 3.526 ± 0.081 ) × 10 1 12 / 15 0.098 ± 0.001 ( 7.273 ± 0.343 ) × 10 1 1 / 7
  1.7 0.092 ± 0.001 ( 2.600 ± 0.066 ) × 10 1 12 / 13 0.086 ± 0.001 ( 5.933 ± 0.223 ) × 10 1 4 / 8
  1.9 0.075 ± 0.001 ( 3.963 ± 0.141 ) × 10 1 12 / 8
  2.1 0.069 ± 0.001 ( 2.440 ± 0.115 ) × 10 1 35 / 9
  0.1 0.124 ± 0.001 ( 1.015 ± 0.025 ) × 10 2 53 / 16 0.136 ± 0.003 ( 1.130 ± 0.104 ) × 10 2 0.2 / 4
  0.3 0.125 ± 0.001 ( 9.894 ± 0.219 ) × 10 1 70 / 16 0.124 ± 0.001 ( 1.460 ± 0.057 ) × 10 2 1 / 7
  0.5 0.124 ± 0.001 ( 9.394 ± 0.208 ) × 10 1 69 / 16 0.121 ± 0.001 ( 1.405 ± 0.044 ) × 10 2 1 / 8
  0.7 0.122 ± 0.001 ( 8.704 ± 0.193 ) × 10 1 43 / 16 0.123 ± 0.001 ( 1.318 ± 0.036 ) × 10 2 1 / 9
  0.9 0.126 ± 0.001 ( 7.596 ± 0.245 ) × 10 1 2 / 7 0.121 ± 0.001 ( 1.235 ± 0.035 ) × 10 2 1 / 9
8.8 1.1 0.119 ± 0.001 ( 6.808 ± 0.245 ) × 10 1 0.2 / 6 0.118 ± 0.001 ( 1.093 ± 0.031 ) × 10 2 0.4 / 9
  1.3 0.114 ± 0.001 ( 5.736 ± 0.229 ) × 10 1 1 / 5 0.113 ± 0.001 ( 9.364 ± 0.273 ) × 10 1 1 / 9
  1.5 0.102 ± 0.001 ( 4.687 ± 0.238 ) × 10 1 0.2 / 4 0.107 ± 0.001 ( 8.245 ± 0.270 ) × 10 1 0.2 / 8
  1.7 0.096 ± 0.001 ( 3.175 ± 0.076 ) × 10 1 11 / 15 0.096 ± 0.001 ( 6.446 ± 0.198 ) × 10 1 4 / 9
  1.9 0.086 ± 0.001 ( 2.337 ± 0.062 ) × 10 1 8 / 13 0.084 ± 0.001 ( 4.803 ± 0.152 ) × 10 1 4 / 9
  2.1 0.076 ± 0.001 ( 3.612 ± 0.137 ) × 10 1 3 / 7
  0.1 0.126 ± 0.001 ( 1.172 ± 0.025 ) × 10 2 88 / 16 0.132 ± 0.002 ( 1.351 ± 0.106 ) × 10 2 0.3 / 4
  0.3 0.126 ± 0.001 ( 1.149 ± 0.025 ) × 10 2 91 / 16 0.128 ± 0.002 ( 1.507 ± 0.076 ) × 10 2 0.4 / 6
  0.5 0.125 ± 0.001 ( 1.123 ± 0.024 ) × 10 2 69 / 16 0.131 ± 0.001 ( 1.379 ± 0.046 ) × 10 2 2 / 8
  0.7 0.124 ± 0.001 ( 1.034 ± 0.022 ) × 10 2 60 / 16 0.124 ± 0.001 ( 1.426 ± 0.046 ) × 10 2 0.4 / 7
  0.9 0.122 ± 0.001 ( 9.590 ± 0.206 ) × 10 1 51 / 16 0.121 ± 0.001 ( 1.334 ± 0.039 ) × 10 2 1 / 8
12.3 1.1 0.122 ± 0.001 ( 8.621 ± 0.185 ) × 10 1 43 / 16 0.122 ± 0.001 ( 1.208 ± 0.033 ) × 10 2 1 / 9
  1.3 0.118 ± 0.001 ( 7.327 ± 0.158 ) × 10 1 27 / 16 0.120 ± 0.001 ( 1.060 ± 0.026 ) × 10 2 1 / 11
  1.5 0.114 ± 0.001 ( 6.270 ± 0.136 ) × 10 1 19 / 16 0.118 ± 0.001 ( 9.376 ± 0.233 ) × 10 1 1 / 11
  1.7 0.109 ± 0.001 ( 5.063 ± 0.112 ) × 10 1 20 / 16 0.113 ± 0.001 ( 7.718 ± 0.020 ) × 10 1 2 / 11
  1.9 0.106 ± 0.001 ( 6.229 ± 0.016 ) × 10 1 3 / 11
  2.1 0.097 ± 0.001 ( 4.876 ± 0.014 ) × 10 1 4 / 10
  0.1 0.127 ± 0.001 ( 1.312 ± 0.028 ) × 10 2 89 / 16 0.137 ± 0.003 ( 1.234 ± 0.145 ) × 10 2 0.05 / 3
  0.3 0.129 ± 0.001 ( 1.293 ± 0.027 ) × 10 2 99 / 16 0.127 ± 0.002 ( 1.514 ± 0.127 ) × 10 2 0.1 / 4
  0.5 0.127 ± 0.001 ( 1.248 ± 0.026 ) × 10 2 93 / 16 0.128 ± 0.002 ( 1.469 ± 0.087 ) × 10 2 0.4 / 5
  0.7 0.127 ± 0.001 ( 1.176 ± 0.025 ) × 10 2 92 / 16 0.126 ± 0.002 ( 1.502 ± 0.064 ) × 10 2 1 / 6
  0.9 0.126 ± 0.001 ( 1.104 ± 0.023 ) × 10 2 87 / 16 0.125 ± 0.002 ( 1.442 ± 0.061 ) × 10 2 2 / 6
17.3 1.1 0.124 ± 0.001 ( 1.022 ± 0.021 ) × 10 2 66 / 16 0.125 ± 0.001 ( 1.342 ± 0.045 ) × 10 2 1 / 7
  1.3 0.122 ± 0.001 ( 9.378 ± 0.971 ) × 10 1 37 / 16 0.120 ± 0.001 ( 1.358 ± 0.045 ) × 10 2 1 / 7
  1.5 0.120 ± 0.001 ( 8.359 ± 0.177 ) × 10 1 27 / 16 0.120 ± 0.001 ( 1.189 ± 0.039 ) × 10 2 0.2 / 7
  1.7 0.116 ± 0.001 ( 7.247 ± 0.154 ) × 10 1 20 / 16 0.117 ± 0.001 ( 1.079 ± 0.035 ) × 10 2 0.4 / 7
  1.9 0.112 ± 0.001 ( 5.898 ± 0.126 ) × 10 1 16 / 16 0.120 ± 0.001 ( 9.162 ± 0.030 ) × 10 1 0.1 / 7
  2.1 0.106 ± 0.001 ( 4.619 ± 0.100 ) × 10 1 20 / 16 0.118 ± 0.001 ( 7.601 ± 0.024 ) × 10 1 0.2 / 7
Table A2. Values of T, V, χ 2 , and ndof corresponding to the curves in Figure 3 and Figure 4 for K and K + produced in p p collisions, respectively.
Table A2. Values of T, V, χ 2 , and ndof corresponding to the curves in Figure 3 and Figure 4 for K and K + produced in p p collisions, respectively.
s (GeV)y K K +
T (GeV)V ( fm 3 ) χ 2 /ndofT (GeV)V ( fm 3 ) χ 2 /ndof
  0.1 0.121 ± 0.001 ( 3.013 ± 0.186 ) × 10 0 2 / 8 0.110 ± 0.002 ( 9.184 ± 0.990 ) × 10 0 1 / 5
  0.3 0.116 ± 0.001 ( 2.824 ± 0.195 ) × 10 0 1 / 7 0.108 ± 0.002 ( 9.001 ± 0.863 ) × 10 0 3 / 6
  0.5 0.115 ± 0.002 ( 2.422 ± 0.284 ) × 10 0 2 / 4 0.098 ± 0.003 ( 6.165 ± 1.215 ) × 10 0 6 / 4
6.3 0.7 0.081 ± 0.001 ( 1.748 ± 0.247 ) × 10 0 3 / 4 0.100 ± 0.001 ( 5.302 ± 0.461 ) × 10 0 18 / 10
  0.9 0.094 ± 0.001 ( 1.190 ± 0.108 ) × 10 0 4 / 6 0.091 ± 0.001 ( 5.201 ± 0.359 ) × 10 0 5 / 10
  1.1 0.098 ± 0.002 ( 5.859 ± 0.722 ) × 10 1 4 / 4 0.099 ± 0.001 ( 2.675 ± 0.166 ) × 10 0 10 / 9
  1.3 0.096 ± 0.005 ( 4.229 ± 1.372 ) × 10 1 0.09 / 2 0.081 ± 0.001 ( 2.467 ± 0.243 ) × 10 0 3 / 6
  1.5 0.073 ± 0.003 ( 2.041 ± 0.654 ) × 10 1 0.7 / 2 0.078 ± 0.001 ( 1.935 ± 0.200 ) × 10 0 5 / 6
  0.1 0.130 ± 0.002 ( 4.775 ± 0.457 ) × 10 0 0.5 / 5 0.137 ± 0.002 ( 1.240 ± 0.087 ) × 10 1 1 / 7
  0.3 0.115 ± 0.002 ( 4.123 ± 0.335 ) × 10 0 2 / 6 0.138 ± 0.002 ( 1.211 ± 0.065 ) × 10 1 1 / 9
  0.5 0.136 ± 0.003 ( 3.434 ± 0.386 ) × 10 0 1 / 4 0.142 ± 0.001 ( 1.083 ± 0.084 ) × 10 1 1 / 6
  0.7 0.111 ± 0.002 ( 3.004 ± 0.244 ) × 10 0 1 / 6 0.146 ± 0.001 ( 9.481 ± 0.497 ) × 10 0 1 / 9
  0.9 0.120 ± 0.002 ( 2.526 ± 0.197 ) × 10 0 1 / 6 0.136 ± 0.001 ( 8.165 ± 0.439 ) × 10 0 1 / 9
7.7 1.1 0.116 ± 0.002 ( 1.703 ± 0.162 ) × 10 0 1 / 5 0.129 ± 0.001 ( 7.598 ± 0.423 ) × 10 0 1 / 9
  1.3 0.108 ± 0.002 ( 1.097 ± 0.108 ) × 10 0 1 / 5 0.130 ± 0.001 ( 6.466 ± 0.358 ) × 10 0 4 / 9
  1.5 0.099 ± 0.003 ( 6.563 ± 1.121 ) × 10 1 0.3 / 3 0.120 ± 0.001 ( 5.394 ± 0.317 ) × 10 0 2 / 9
  1.7 0.063 ± 0.002 ( 2.875 ± 0.807 ) × 10 1 0.4 / 2 0.110 ± 0.001 ( 4.168 ± 0.277 ) × 10 0 1 / 8
  1.9 0.095 ± 0.001 ( 3.140 ± 0.256 ) × 10 0 2 / 7
  2.1 0.078 ± 0.001 ( 2.172 ± 0.229 ) × 10 0 7 / 6
  0.1 0.124 ± 0.002 ( 5.460 ± 0.357 ) × 10 0 3 / 7 0.124 ± 0.001 ( 1.278 ± 0.072 ) × 10 1 3 / 9
  0.3 0.128 ± 0.002 ( 5.890 ± 0.431 ) × 10 0 0.4 / 6 0.126 ± 0.001 ( 1.259 ± 0.076 ) × 10 1 1 / 8
  0.5 0.133 ± 0.002 ( 4.015 ± 0.250 ) × 10 0 3 / 7 0.142 ± 0.001 ( 9.172 ± 0.468 ) × 10 0 2 / 9
  0.7 0.119 ± 0.002 ( 3.591 ± 0.245 ) × 10 0 1 / 7 0.139 ± 0.001 ( 8.748 ± 0.463 ) × 10 0 1 / 9
8.8 0.9 0.124 ± 0.002 ( 3.099 ± 0.239 ) × 10 0 0.4 / 6 0.140 ± 0.001 ( 8.885 ± 0.470 ) × 10 0 0.5 / 9
  1.1 0.134 ± 0.002 ( 2.426 ± 0.185 ) × 10 0 1 / 6 0.146 ± 0.002 ( 8.079 ± 0.418 ) × 10 0 1 / 9
  1.3 0.106 ± 0.002 ( 1.478 ± 0.144 ) × 10 0 1 / 5 0.146 ± 0.001 ( 6.816 ± 0.352 ) × 10 0 1 / 9
  1.5 0.089 ± 0.002 ( 9.181 ± 1.199 ) × 10 1 0.2 / 4 0.131 ± 0.001 ( 5.690 ± 0.313 ) × 10 0 1 / 9
  1.7 0.091 ± 0.002 ( 5.613 ± 0.733 ) × 10 1 2 / 4 0.128 ± 0.001 ( 4.577 ± 0.252 ) × 10 0 2 / 9
  1.9 0.080 ± 0.003 ( 3.513 ± 0.860 ) × 10 1 1 / 2 0.132 ± 0.001 ( 3.929 ± 0.228 ) × 10 0 5 / 8
  0.1 0.125 ± 0.002 ( 7.420 ± 0.487 ) × 10 0 1 / 7 0.129 ± 0.002 ( 1.404 ± 0.093 ) × 10 1 1 / 7
  0.3 0.124 ± 0.002 ( 7.399 ± 0.563 ) × 10 0 0.3 / 6 0.139 ± 0.002 ( 1.419 ± 0.075 ) × 10 1 2 / 9
  0.5 0.132 ± 0.002 ( 6.808 ± 0.504 ) × 10 0 0.4 / 6 0.139 ± 0.001 ( 1.398 ± 0.073 ) × 10 1 1 / 9
  0.7 0.132 ± 0.002 ( 6.064 ± 0.452 ) × 10 0 0.4 / 6 0.152 ± 0.002 ( 1.278 ± 0.070 ) × 10 1 0.1 / 8
12.3 0.9 0.115 ± 0.001 ( 4.823 ± 0.380 ) × 10 0 0.5 / 6 0.139 ± 0.002 ( 1.128 ± 0.072 ) × 10 1 1 / 7
  1.1 0.113 ± 0.001 ( 3.704 ± 0.262 ) × 10 0 1 / 7 0.140 ± 0.002 ( 1.005 ± 0.523 ) × 10 1 0.4 / 9
  1.3 0.111 ± 0.001 ( 2.861 ± 0.207 ) × 10 0 0.4 / 7 0.137 ± 0.001 ( 8.783 ± 0.463 ) × 10 0 1 / 9
  1.5 0.113 ± 0.002 ( 2.064 ± 0.196 ) × 10 0 0.08 / 5 0.134 ± 0.001 ( 7.117 ± 0.375 ) × 10 0 0.4 / 9
  1.7 0.094 ± 0.002 ( 1.293 ± 0.169 ) × 10 0 2 / 4 0.133 ± 0.001 ( 6.053 ± 0.354 ) × 10 0 0.3 / 8
  0.1 0.137 ± 0.002 ( 9.281 ± 0.516 ) × 10 0 0.4 / 8 0.128 ± 0.001 ( 1.480 ± 0.088 ) × 10 1 2 / 8
  0.3 0.137 ± 0.002 ( 9.562 ± 0.601 ) × 10 0 0.05 / 7 0.129 ± 0.002 ( 1.512 ± 0.090 ) × 10 1 1 / 8
  0.5 0.130 ± 0.002 ( 8.807 ± 0.511 ) × 10 0 0.3 / 8 0.133 ± 0.002 ( 1.566 ± 0.092 ) × 10 1 3 / 8
  0.7 0.135 ± 0.002 ( 8.375 ± 0.527 ) × 10 0 0.1 / 7 0.149 ± 0.002 ( 1.543 ± 0.086 ) × 10 1 1 / 8
  0.9 0.136 ± 0.001 ( 7.181 ± 0.381 ) × 10 0 1 / 9 0.143 ± 0.002 ( 1.331 ± 0.075 ) × 10 1 0.2 / 8
17.3 1.1 0.121 ± 0.001 ( 6.036 ± 0.346 ) × 10 0 1 / 44 0.139 ± 0.002 ( 1.149 ± 0.065 ) × 10 1 0.4 / 8
  1.3 0.115 ± 0.001 ( 4.485 ± 0.283 ) × 10 0 2 / 8 0.140 ± 0.002 ( 9.986 ± 0.633 ) × 10 0 0.1 / 7
  1.5 0.115 ± 0.001 ( 3.557 ± 0.225 ) × 10 0 0.5 / 8 0.133 ± 0.002 ( 8.618 ± 0.554 ) × 10 0 1 / 7
  1.7 0.102 ± 0.001 ( 2.318 ± 0.189 ) × 10 0 4 / 6 0.127 ± 0.002 ( 6.873 ± 0.449 ) × 10 0 0.4 / 7
  1.9 0.098 ± 0.002 ( 1.667 ± 0.167 ) × 10 0 0.1 / 5 0.124 ± 0.002 ( 5.261 ± 0.401 ) × 10 0 1 / 6
  2.1 0.076 ± 0.001 ( 0.833 ± 0.115 ) × 10 0 2 / 4 0.118 ± 0.002 ( 2.968 ± 0.205 ) × 10 0 6 / 7
Table A3. Values of T, V, χ 2 , and ndof corresponding to the curves in Figure 5 and Figure 6 for p ¯ and p produced in p p collisions, respectively.
Table A3. Values of T, V, χ 2 , and ndof corresponding to the curves in Figure 5 and Figure 6 for p ¯ and p produced in p p collisions, respectively.
s (GeV)y p ¯ p
T (GeV)V ( fm 3 ) χ 2 /ndofT (GeV)V ( fm 3 ) χ 2 /ndof
  0.1 0.104 ± 0.001 ( 1.662 ± 0.245 ) × 10 1 4 / 7
  0.3 0.095 ± 0.001 ( 1.993 ± 0.244 ) × 10 1 9 / 10
  0.5 0.091 ± 0.001 ( 2.176 ± 0.286 ) × 10 1 8 / 10
  0.7 0.097 ± 0.001 ( 2.359 ± 0.273 ) × 10 1 17 / 10
6.3 0.9 0.089 ± 0.001 ( 3.633 ± 0.520 ) × 10 1 4 / 9
  1.1 0.090 ± 0.001 ( 3.326 ± 0.630 ) × 10 1 3 / 6
  1.3 0.074 ± 0.001 ( 4.321 ± 0.963 ) × 10 1 1 / 6
  1.5 0.064 ± 0.001 ( 4.180 ± 1.066 ) × 10 1 2 / 6
  1.7 0.043 ± 0.001 ( 2.855 ± 1.225 ) × 10 1 0.3 / 6
  0.1 0.114 ± 0.002 ( 3.326 ± 0.458 ) × 10 1 9 / 7 0.135 ± 0.002 ( 1.270 ± 0.121 ) × 10 1 1 / 9
  0.3 0.109 ± 0.002 ( 3.236 ± 0.457 ) × 10 1 14 / 7 0.133 ± 0.002 ( 1.405 ± 0.136 ) × 10 1 1 / 9
  0.5 0.109 ± 0.002 ( 2.025 ± 0.332 ) × 10 1 10 / 6 0.133 ± 0.002 ( 1.704 ± 0.165 ) × 10 1 1 / 9
  0.7 0.111 ± 0.003 ( 1.433 ± 0.301 ) × 10 1 10 / 5 0.125 ± 0.002 ( 2.061 ± 0.211 ) × 10 1 0.3 / 9
7.7 0.9 0.065 ± 0.001 ( 6.943 ± 1.945 ) × 10 2 12 / 5 0.124 ± 0.002 ( 2.637 ± 0.272 ) × 10 1 0.2 / 9
  1.1 0.126 ± 0.002 ( 3.075 ± 0.302 ) × 10 1 2 / 9
  1.3 0.114 ± 0.001 ( 3.436 ± 0.384 ) × 10 1 1 / 9
  1.5 0.104 ± 0.001 ( 3.801 ± 0.505 ) × 10 1 1 / 8
  1.7 0.086 ± 0.001 ( 4.103 ± 0.783 ) × 10 1 1 / 6
  0.1 0.105 ± 0.003 ( 3.835 ± 1.138 ) × 10 1 6 / 3 0.135 ± 0.002 ( 1.318 ± 0.125 ) × 10 1 0.2 / 9
  0.3 0.112 ± 0.002 ( 3.713 ± 0.535 ) × 10 1 4 / 7 0.137 ± 0.002 ( 1.417 ± 0.134 ) × 10 1 0.3 / 9
  0.5 0.104 ± 0.002 ( 2.982 ± 0.647 ) × 10 1 4 / 5 0.133 ± 0.002 ( 1.632 ± 0.156 ) × 10 1 1 / 9
  0.7 0.131 ± 0.002 ( 1.682 ± 0.241 ) × 10 1 29 / 6 0.127 ± 0.002 ( 1.928 ± 0.194 ) × 10 1 0.1 / 9
8.8 0.9 0.094 ± 0.002 ( 8.812 ± 1.583 ) × 10 2 11 / 6 0.128 ± 0.002 ( 2.374 ± 0.236 ) × 10 1 0.08 / 9
  1.1 0.127 ± 0.002 ( 2.931 ± 0.293 ) × 10 1 0.08 / 9
  1.3 0.125 ± 0.002 ( 3.365 ± 0.342 ) × 10 1 1 / 9
  1.5 0.112 ± 0.001 ( 3.630 ± 0.411 ) × 10 1 0.2 / 9
  1.7 0.097 ± 0.001 ( 3.748 ± 0.584 ) × 10 1 1 / 7
  1.9 0.085 ± 0.001 ( 3.785 ± 0.730 ) × 10 1 1 / 6
  0.1 0.100 ± 0.001 ( 9.128 ± 1.170 ) × 10 1 11 / 9 0.145 ± 0.002 ( 7.934 ± 0.667 ) × 10 0 1 / 10
  0.3 0.101 ± 0.001 ( 9.005 ± 1.106 ) × 10 1 5 / 9 0.148 ± 0.002 ( 8.276 ± 0.678 ) × 10 0 1 / 10
  0.5 0.104 ± 0.002 ( 7.481 ± 0.905 ) × 10 1 2 / 9 0.141 ± 0.002 ( 9.202 ± 0.872 ) × 10 0 1 / 9
  0.7 0.112 ± 0.001 ( 6.565 ± 0.747 ) × 10 1 2 / 9 0.145 ± 0.002 ( 1.100 ± 0.101 ) × 10 1 0.2 / 9
12.3 0.9 0.103 ± 0.001 ( 4.409 ± 0.586 ) × 10 1 5 / 8 0.142 ± 0.002 ( 1.324 ± 0.105 ) × 10 1 1 / 11
  1.1 0.111 ± 0.001 ( 2.464 ± 0.303 ) × 10 1 4 / 8 0.135 ± 0.001 ( 1.683 ± 0.139 ) × 10 1 0.1 / 11
  1.3 0.103 ± 0.002 ( 1.558 ± 0.232 ) × 10 1 8 / 7 0.132 ± 0.001 ( 2.091 ± 0.177 ) × 10 1 0.3 / 11
  1.5 0.100 ± 0.002 ( 7.733 ± 1.344 ) × 10 2 11 / 6 0.127 ± 0.001 ( 2.468 ± 0.218 ) × 10 1 1 / 11
  1.7 0.119 ± 0.001 ( 2.891 ± 0.274 ) × 10 1 2 / 11
  1.9 0.104 ± 0.001 ( 3.098 ± 0.434 ) × 10 1 1 / 8
  0.1 0.126 ± 0.002 ( 1.719 ± 0.206 ) × 10 0 0.4 / 8 0.150 ± 0.002 ( 7.374 ± 0.755 ) × 10 0 0.2 / 8
  0.3 0.102 ± 0.001 ( 1.499 ± 0.174 ) × 10 0 11 / 10 0.150 ± 0.002 ( 7.430 ± 0.760 ) × 10 0 0.2 / 8
  0.5 0.126 ± 0.001 ( 1.544 ± 0.149 ) × 10 0 3 / 10 0.150 ± 0.002 ( 7.882 ± 0.810 ) × 10 0 0.08 / 8
  0.7 0.133 ± 0.002 ( 1.438 ± 0.144 ) × 10 0 1 / 9 0.143 ± 0.002 ( 9.096 ± 0.964 ) × 10 0 0.2 / 8
  0.9 0.121 ± 0.002 ( 1.086 ± 0.136 ) × 10 0 2 / 8 0.151 ± 0.002 ( 1.061 ± 0.109 ) × 10 1 0.5 / 8
17.3 1.1 0.127 ± 0.002 ( 0.805 ± 0.111 ) × 10 0 1 / 7 0.141 ± 0.002 ( 1.191 ± 0.149 ) × 10 1 0.09 / 7
  1.3 0.123 ± 0.002 ( 0.568 ± 0.066 ) × 10 0 4 / 8 0.138 ± 0.002 ( 1.555 ± 0.196 ) × 10 1 0.04 / 7
  1.5 0.161 ± 0.005 ( 0.410 ± 0.081 ) × 10 0 1 / 5 0.135 ± 0.002 ( 1.950 ± 0.248 ) × 10 1 0.02 / 7
  1.7 0.136 ± 0.005 ( 0.214 ± 0.064 ) × 10 0 2 / 4 0.136 ± 0.002 ( 2.515 ± 0.392 ) × 10 1 0.2 / 6
  1.9 0.071 ± 0.003 ( 0.057 ± 0.029 ) × 10 0 1 / 3 0.142 ± 0.003 ( 2.868 ± 0.491 ) × 10 1 0.001 / 5
  2.1 0.132 ± 0.003 ( 3.033 ± 0.531 ) × 10 1 0.4 / 5

References

  1. Arsene, I.; Bearden, I.G.; Beavis, D.; Besliu, C.; Budick, B.; Bggild, H.; Chasman, C.; Chasman, C.; Christensen, C.H.; Christiansen, P.; et al. Quark gluon plasma and color glass condensate at RHIC? The perspective from the BRAHMS experiment. Nucl. Phys. A 2005, 757, 1–27. [Google Scholar] [CrossRef]
  2. Adcox, K.; Adler, S.S.; Afanasiev, S.; Aidala, C.; Ajitanand, N.N.; Akiba, Y.; Al-Jamel, A.; Alexander, J.; Amirikas, R.; Aoki, K.; et al. Formation of dense partonic matter in relativistic nucleus-nucleus collisions at RHIC: Experimental evaluation by the PHENIX collaboration. Nucl. Phys. A 2005, 757, 184–283. [Google Scholar] [CrossRef]
  3. Adams, J.; Aggarwal, M.M.; Ahammed, Z.; Amonett, J.; Anderson, B.D.; Arkhipkin, D.; Averichev, G.S.; Badyal, S.K.; Bai, Y.; Balewski, J.; et al. Experimental and theoretical challenges in the search for the quark gluon plasma: The STAR collaboration’s critical assessment of the evidence from RHIC collisions. Nucl. Phys. A 2005, 757, 102–183. [Google Scholar] [CrossRef]
  4. Schukraft, J. Heavy ion physics with the ALICE experiment at the CERN LHC. Phil. Trans. R. Soc. Lond. A 2012, 370, 917–932. [Google Scholar]
  5. Braun-Munzinger, P.; Stachel, J. The quest for the quark-gluon plasma. Nature 2007, 448, 302–309. [Google Scholar] [CrossRef] [PubMed]
  6. Harris, J.W.; Müller, B. The search for the quark-gluon plasma. Ann. Rev. Nucl. Part. Sci. 1996, 46, 71–107. [Google Scholar] [CrossRef]
  7. Heinz, U.W.; Jacob, M. Evidence for a new state of matter: An assessment of the results from the CERN lead beam programme. arXiv 2000, arXiv:nucl-th/0002042v1. [Google Scholar]
  8. Podlaski, P. Results on system size dependence of strangeness production in the CERN SPS energy range from NA61/SHINE. EPJ Web Conf. 2023, 276, 03008. [Google Scholar] [CrossRef]
  9. Tannenbaum, M.J. Recent results in relativistic heavy ion collisions: From ‘a new state of matter’ to ‘the perfect fluid’. Rept. Prog. Phys. 2006, 69, 2005–2060. [Google Scholar] [CrossRef]
  10. Abdulhamid, M.I.; Aboona, B.E.; Adam, J.; Adamczyk, L.; Adams, J.R.; Aggarwal, I.; Aggarwal, M.M.; Ahammed, Z.; Anderson, D.M.; Aschenauer, E.C.; et al. Energy dependence of intermittency for charged hadrons in Au+Au collisions at RHIC. Phys. Lett. B 2023, 845, 138165. [Google Scholar] [CrossRef]
  11. Abdulhamid, M.I.; Aboona, B.E.; Adam, J.; Adamczyk, L.; Adams, J.R.; Agakishiev, G.; Aggarwal, I.; Aggarwal, M.M.; Ahammed, Z.; Aitbaev, A. Beam energy dependence of triton production and yield ratio (Nt × Np/ N d 2 ) in Au+Au collisions at RHIC. Phys. Rev. Lett. 2023, 130, 202301. [Google Scholar] [CrossRef] [PubMed]
  12. Aboona, B.E.; Adam, J.; Adamczyk, L.; Adams, J.R.; Aggarwal, I.; Aggarwal, M.M.; Ahammed, Z.; Anderson, D.M.; Aschenauer, E.C.; Atchison, J. Search for the chiral magnetic effect in Au+Au collisions at s N N = 27 GeV with the STAR forward event plane detectors. Phys. Lett. B 2023, 839, 137779. [Google Scholar] [CrossRef]
  13. Adam, J.; Adamczyk, L.; Adams, J.R.; Adkins, J.K.; Agakishiev, G.; Aggarwal, M.M.; Ahammed, Z.; Alekseev, I.; Anderson, D.M.; Aoyama, R.; et al. Bulk properties of the system formed in Au + Au collisions at s N N = 14.5 GeV at the BNL STAR detector. Phys. Rev. C 2020, 101, 024905. [Google Scholar] [CrossRef]
  14. Helmut, S. Limits of confinement: The first 15 years of ultra-relativistic heavy ion studies. Nucl. Phys. A 2003, 715, 3c–19c. [Google Scholar]
  15. Back, B.B.; Baker, M.D.; Ballintijn, M.; Barton, D.S.; Becker, B.; Betts, R.R.; Bickley, A.A.; Bindel, R.; Budzanowski, A.; Busza, W.; et al. The PHOBOS perspective on discoveries at RHIC. Nucl. Phys. A 2005, 757, 28–101. [Google Scholar] [CrossRef]
  16. Aamodt, K.; Abelev, B.; Quintana, A.A.; Adamová, D.; Adare, A.M.; Aggarwal, M.M.; Rinella, G.A.; Agocs, A.G.; Salazar, S.A.; Ahammed, Z.; et al. Elliptic flow of charged particles in Pb-Pb collisions at 2.76 TeV. Phys. Rev. Lett. 2010, 105, 252302. [Google Scholar] [CrossRef] [PubMed]
  17. Aamodt, K.; Abelev, B.; Quintana, A.A.; Adamová, D.; Adare, A.M.; Aggarwal, M.M.; Rinella, G.A.; Agocs, A.G.; Salazar, S.A.; Ahammed, Z.; et al. Charged-particle multiplicity density at mid-rapidity in central Pb-Pb collisions at s N N = 2.76 TeV. Phys. Rev. Lett. 2010, 105, 252301. [Google Scholar] [CrossRef]
  18. Aamodt, K.; Quintana, A.A.; Adamová, D.; Adare, A.M.; Aggarwal, M.M.; Rinella, G.A.; Agocs, A.G.; Salazar, S.A.; Ahammed, Z.; Ahmad, N.; et al. Centrality dependence of the charged-particle multiplicity density at mid-rapidity in Pb-Pb collisions at s N N = 2.76 TeV. Phys. Rev. Lett. 2011, 106, 032301. [Google Scholar] [CrossRef]
  19. Aamodt, K.; Abelev, B.; Quintana, A.A.; Adamová, D.; Adare, A.M.; Aggarwal, M.M.; Rinella, G.A.; Agocs, A.G.; Agostinelli, A.; Salazar, S.A.; et al. Higher harmonic anisotropic flow measurements of charged particles in Pb-Pb collisions at s N N = 2.76 TeV. Phys. Rev. Lett. 2011, 107, 032301. [Google Scholar] [CrossRef]
  20. Adler, S.S.; Afanasiev, S.; Aidala, C.; Ajitanand, N.N.; Akiba, Y.; Alexander, J.; Amirikas, R.; Aphecetche, L.; Aronson, S.H.; Averbeck, R.; et al. Identified charged particle spectra and yields in Au+Au collisions at s N N = 200 GeV. Phys. Rev. C 2004, 69, 034909. [Google Scholar] [CrossRef]
  21. Van Hove, L. Multiplicity dependence of pt spectrum as a possible signal for a phase transition in hadronic collisions. Phys. Lett. B 1982, 118, 138–140. [Google Scholar] [CrossRef]
  22. Andronic, A.; Braun-Munzinger, P.; Stachel, J. Hadron production in central nucleus-nucleus collisions at chemical freeze-out. Nucl. Phys. A 2006, 772, 167–199. [Google Scholar] [CrossRef]
  23. Cleymans, J.; Oeschler, H.; Redlich, K.; Wheaton, S. Comparison of chemical freeze-out criteria in heavy-ion collisions. Phys. Rev. C 2006, 73, 034905. [Google Scholar] [CrossRef]
  24. Andronic, A.; Braun-Munzinger, P.; Stachel, J. Thermal hadron production in relativistic nuclear collisions. Acta Phys. Pol. B 2009, 40, 1005–1012. [Google Scholar]
  25. Andronic, A.; Braun-Munzinger, P.; Stachel, J. The horn, the hadron mass spectrum and the QCD phase diagram: The statistical model of hadron production in central nucleus-nucleus collisions. Nucl. Phys. A 2010, 834, 237c–240c. [Google Scholar] [CrossRef]
  26. Redlich, K.; Cleymans, J.; Oeschler, H.; Tounsi, A. Particle production and equilibration in heavy ion collisions. Acta Phys. Polon. B 2002, 33, 1609–1628. [Google Scholar]
  27. Wheaton, S.; Cleymans, J.; Hauer, M. THERMUS: A thermal model package for ROOT. Comput. Phys. Commun. 2009, 180, 84–106. [Google Scholar] [CrossRef]
  28. Andronic, A.; Beutler, F.; Braun-Munzinger, P.; Redlich, K.; Stachel, J. Statistical hadronization of heavy flavor quarks in elementary collisions: Successes and failures. Phys. Lett. B 2009, 678, 350–354. [Google Scholar] [CrossRef]
  29. Aduszkiewicz, A.; Ali, Y.; Andronov, E.; Antićić, T.; Baatar, B.; Baszczyk, M.; Bhosale, S.; Blondel, A.; Bogomilov, M.; Brandin, A.; et al. Measurements of π±, K±, p and p ¯ spectra in proton–proton interactions at 20, 31, 40, 80 and 158 GeV/c with the NA61/SHINE spectrometer at the CERN SPS. Eur. Phys. J. C 2017, 77, 671. [Google Scholar] [CrossRef]
  30. Cleymans, J.; Worku, D. Relativistic thermodynamics: Transverse momentum distributions in high-energy physics. Eur. Phys. J. A 2012, 48, 160. [Google Scholar] [CrossRef]
  31. Abelev, B.I.; Adams, J.; Aggarwal, M.M.; Ahammed, Z.; Amonett, J.; Anderson, B.D.; Anderson, M.; Arkhipkin, D.; Averichev, G.S.; Bai, Y.; et al. Strange particle production in p + p collisions at s = 200 GeV. Phys. Rev. C 2007, 75, 064901. [Google Scholar] [CrossRef]
  32. Tsallis, C. Possible generalization of Boltzmann-Gibbs statistics. J. Stat. Phys. 1988, 52, 479–487. [Google Scholar] [CrossRef]
  33. Biro, T.S.; Purcsel, G.; Urmossy, K. Non-extensive approach to quark matter. Eur. Phys. J. A 2009, 40, 325–340. [Google Scholar] [CrossRef]
  34. Zheng, H.; Zhu, L.L.; Bonasera, A. Systematic analysis of hadron spectra in p + p collisions using Tsallis distributions. Phys. Rev. D 2015, 92, 074009. [Google Scholar] [CrossRef]
  35. Zheng, H.; Zhu, L.L. Can Tsallis distribution fit all the particle spectra produced at RHIC and LHC? Adv. High Energy Phys. 2015, 2015, 180491. [Google Scholar] [CrossRef]
  36. Wang, X.N.; Hwa, R.C. The effect of jet production on the multiplicity dependence of average transverse momentum. Phys. Rev. D 1989, 39, 187–194. [Google Scholar] [CrossRef] [PubMed]
  37. Sjöstrand, T.; van Zijl, M. A multiple-interaction model for the event structure in hadron collisions. Phys. Rev. D 1987, 36, 2019–2041. [Google Scholar] [CrossRef]
  38. Parvan, A.S. Non-extensive statistics effects in transverse momentum spectra of hadrons. arXiv 2015, arXiv:1502.01581. [Google Scholar]
  39. Rath, R.; Khuntia, A.; Sahoo, R.; Cleymans, J. Event multiplicity, transverse momentum and energy dependence of charged particle production, and system thermodynamics in pp collisions at the Large Hadron Collider. J. Phys. G 2020, 47, 055111. [Google Scholar] [CrossRef]
  40. Schnedermann, E.; Sollfrank, J.; Heinz, U. Thermal phenomenology of hadrons from 200A GeV S+S collisions. Phys. Rev. C 1993, 48, 2462–2475. [Google Scholar] [CrossRef]
  41. Hagedorn, R. Multiplicities, pT distributions and the expected hadron ⟶ quark-gluon phase transition. Riv. Nuovo Cimento 1983, 6, 1–50. [Google Scholar] [CrossRef]
  42. Adamczyk, L.; Adkins, J.K.; Agakishiev, G.; Aggarwal, M.M.; Ahammed, Z.; Ajitanand, N.N.; Alekseev, I.; Anderson, D.M.; Aoyama, R.; Aparin, A.; et al. Bulk properties of the medium produced in relativistic heavy-ion collisions from the beam energy scan program. Phys. Rev. C 2017, 96, 044904. [Google Scholar] [CrossRef]
  43. Braun-Munzinger, P.; Stachel, J.; Wessels, J.P.; Xu, N. Thermal equilibration and expansion in nucleus-nucleus collisions at the AGS. Phys. Lett. B 1995, 344, 43–48. [Google Scholar] [CrossRef]
  44. Andronic, A.; Braun-Munzinger, P.; Stachel, J. Thermal hadron production in relativistic nuclear collisions: The hadron mass spectrum, the horn, and the QCD phase transition. Phys. Lett. B 2009, 673, 142–145. [Google Scholar] [CrossRef]
  45. Abelev, B.I.; Aggarwal, M.M.; Ahammed, Z.; Anderson, B.D.; Arkhipkin, D.; Averichev, G.S.; Bai, Y.; Balewski, J.; Barannikova, O.; Barnby, L.S.; et al. Systematic measurements of identified particle spectra in pp, d+Au and Au+Au collisions from STAR. Phys. Rev. C 2009, 79, 034909. [Google Scholar] [CrossRef]
  46. Cleymans, J.; Oeschler, H.; Redlich, K. Influence of impact parameter on thermal description of relativistic heavy ion collisions at (1–2)A GeV. Phys. Rev. C 1999, 59, 1663–1673. [Google Scholar] [CrossRef]
  47. Braun-Munzinger, P.; Heppe, I.; Stachel, J. Chemical equilibration in Pb+Pb collisions at the SPS. Phys. Lett. B 1999, 465, 15–20. [Google Scholar] [CrossRef]
  48. Manninen, J.; Becattini, F. Chemical freeze-out in ultra-relativistic heavy ion collisions at s N N = 130 and 200 GeV. Phys. Rev. C 2008, 78, 054901. [Google Scholar] [CrossRef]
  49. Andronic, A.; Braun-Munzinger, P.; Redlich, K. Decoding the phase structure of QCD via particle production at high energy. Nature 2018, 561, 321–330. [Google Scholar] [CrossRef]
  50. Koch, P.; Rafelski, J.; Greiner, W. Strange hadron in hot nuclear matter. Phys. Lett. B 1983, 123, 321–330. [Google Scholar] [CrossRef]
  51. Braun-Munzinger, P.; Magestro, D.; Redlich, K.; Stachel, J. Hadron production in Au-Au collisions at RHIC. Phys. Lett. B 2001, 518, 41–46. [Google Scholar] [CrossRef]
  52. Lao, H.-L.; Gao, Y.-Q.; Liu, F.-H. Energy dependent chemical potentials of light particles and quarks from yield ratios of antiparticles to particles in high energy collisions. Universe 2019, 5, 152. [Google Scholar] [CrossRef]
  53. Lao, H.-L.; Gao, Y.-Q.; Liu, F.-H. Light particle and quark chemical potentials from negatively to positively charged particle yield ratios corrected by removing strong and weak decays. Adv. High Energy Phys. 2020, 2020, 5064737. [Google Scholar] [CrossRef]
  54. Gardim, F.G.; Giacalone, G.; Luzum, M.; Ollitrault, J.Y. Thermodynamics of hot strong-interaction matter from ultrarelativistic nuclear collisions. Nat. Phys. 2020, 16, 615–619. [Google Scholar] [CrossRef]
  55. Waqas, M.; Peng, G.-X.; Ajaz, M.; Haj Ismail, A.; Wazir, Z.; Li, L.-L. Extraction of different temperatures and kinetic freeze-out volume in high energy collisions. J. Phys. G 2022, 49, 095102. [Google Scholar] [CrossRef]
  56. Gutay, L.J.; Hirsch, A.S.; Pajares, C.; Scharenberg, R.P.; Srivastava, B.K. De-confinement in small systems: Clustering of color sources in high multiplicity p ¯ p collisions at s = 1.8 TeV. Int. J. Mod. Phys. E 2015, 24, 1550101. [Google Scholar] [CrossRef]
  57. Scharenberg, R.P.; Srivastava, B.K.; Pajares, C. Exploring the initial stage of high multiplicity proton–proton collisions by determining the initial temperature of the quark-gluon plasma. Phys. Rev. D 2019, 100, 114040. [Google Scholar] [CrossRef]
  58. Sahoo, P.; De, S.; Tiwari, S.K.; Sahoo, R. Energy and centrality dependent study of deconfinement phase transition in a color string percolation approach at RHIC energies. Eur. Phys. J. A 2018, 54, 136. [Google Scholar] [CrossRef]
  59. Yu, N.; Luo, X.F. Particle decay from statistical thermal model in high-energy nucleus-nucleus collisions. Eur. Phys. J. A 2019, 55, 26. [Google Scholar] [CrossRef]
  60. Biyajima, M.; Kaneyama, M.; Mizoguchi, T.; Wilk, G. Analyses of kt distributions at RHIC by means of some selected statistical and stochastic models. Eur. Phys. J. C 2005, 40, 243–250. [Google Scholar] [CrossRef]
  61. Cleymans, J.; Hamar, G.; Levai, P.; Wheaton, S. Near-thermal equilibrium with Tsallis distributions in heavy-ion collisions. J. Phys. G 2009, 36, 064018. [Google Scholar] [CrossRef]
  62. Shao, M.; Yi, L.; Tang, Z.B.; Chen, H.F.; Li, C.; Xu, Z.B. Examination of the species and beam energy dependence of particle spectra using Tsallis statistics. J. Phys. G 2010, 37, 085104. [Google Scholar] [CrossRef]
  63. Wong, C.-Y.; Wilk, G.; Cirto, L.J.L.; Tsallis, C. From QCD-based hard-scattering to nonextensive statistical mechanical descriptions of transverse momentum spectra in high-energy pp and p p ¯ collisions. Phys. Rev. D 2015, 91, 114027. [Google Scholar] [CrossRef]
  64. Hui, J.-Q.; Jiang, Z.-J.; Xu, D.-F. A description of the transverse momentum distributions of charged particles produced in heavy ion collisions at RHIC and LHC energies. Adv. High Energy Phys. 2018, 2018, 7682325. [Google Scholar] [CrossRef]
  65. Tripathy, S.; Tiwari, S.K.; Younus, M.; Sahoo, R. Elliptic flow in Pb+Pb collisions at s N N = 2.76 TeV at the LHC using Boltzmann transport equation with non-extensive statistics. Eur. Phys. J. A 2018, 54, 38. [Google Scholar] [CrossRef]
  66. Rybczyński, M.; Włodarczyk, Z. Tsallis statistics approach to the transverse momentum distributions in p-p collisions. Eur. Phys. J. C 2014, 74, 2785. [Google Scholar] [CrossRef]
  67. Zhang, S.; Ma, Y.G.; Chen, J.H.; Zhong, C. Beam energy dependence of Hanbury-Brown-Twiss radii from a blast-wave model. Adv. High Energy Phys. 2016, 2016, 9414239. [Google Scholar] [CrossRef]
  68. Hogg, R.V.; McKean, J.W.; Craig, A.T. Introduction to Mathematical Statistics, 8th ed.; Pearson: Boston, MA, USA, 2018. [Google Scholar]
Figure 1. The spectra of π produced in p p collisions at s = (a) 6.3, (b) 7.7, (c) 8.8, (d) 12.3, and (e) 17.3 GeV at different y with an interval width of 0.2. The symbols represent the experimental data measured by the NA61/SHINE Collaboration [29] and the curves are the fitting results from the Bose–Einstein distribution.
Figure 1. The spectra of π produced in p p collisions at s = (a) 6.3, (b) 7.7, (c) 8.8, (d) 12.3, and (e) 17.3 GeV at different y with an interval width of 0.2. The symbols represent the experimental data measured by the NA61/SHINE Collaboration [29] and the curves are the fitting results from the Bose–Einstein distribution.
Entropy 25 01571 g001
Figure 2. The spectra of π + produced in p p collisions at s = (a) 6.3, (b) 7.7, (c) 8.8, (d) 12.3, and (e) 17.3 GeV at different y. The symbols represent the experimental data measured by the NA61/SHINE Collaboration [29] and the curves are the fitting results from the Bose–Einstein distribution.
Figure 2. The spectra of π + produced in p p collisions at s = (a) 6.3, (b) 7.7, (c) 8.8, (d) 12.3, and (e) 17.3 GeV at different y. The symbols represent the experimental data measured by the NA61/SHINE Collaboration [29] and the curves are the fitting results from the Bose–Einstein distribution.
Entropy 25 01571 g002
Figure 3. The spectra of K produced in p p collisions at s = (a) 6.3, (b) 7.7, (c) 8.8, (d) 12.3, and (e) 17.3 GeV at different y. The symbols represent the experimental data measured by the NA61/SHINE Collaboration [29] and the curves are the fitting results from the Bose–Einstein distribution.
Figure 3. The spectra of K produced in p p collisions at s = (a) 6.3, (b) 7.7, (c) 8.8, (d) 12.3, and (e) 17.3 GeV at different y. The symbols represent the experimental data measured by the NA61/SHINE Collaboration [29] and the curves are the fitting results from the Bose–Einstein distribution.
Entropy 25 01571 g003
Figure 4. The spectra of K + produced in p p collisions at s = (a) 6.3, (b) 7.7, (c) 8.8, (d) 12.3, and (e) 17.3 GeV at different y. The symbols represent the experimental data measured by the NA61/SHINE Collaboration [29] and the curves are the fitting results from the Bose–Einstein distribution.
Figure 4. The spectra of K + produced in p p collisions at s = (a) 6.3, (b) 7.7, (c) 8.8, (d) 12.3, and (e) 17.3 GeV at different y. The symbols represent the experimental data measured by the NA61/SHINE Collaboration [29] and the curves are the fitting results from the Bose–Einstein distribution.
Entropy 25 01571 g004
Figure 5. The spectra of p ¯ produced in p p collisions at s = (a) 7.7, (b) 8.8, (c) 12.3, and (d) 17.3 at different y. The symbols represent the experimental data measured by the NA61/SHINE Collaboration [29] and the curves are the fitting results from the Fermi–Dirac distribution.
Figure 5. The spectra of p ¯ produced in p p collisions at s = (a) 7.7, (b) 8.8, (c) 12.3, and (d) 17.3 at different y. The symbols represent the experimental data measured by the NA61/SHINE Collaboration [29] and the curves are the fitting results from the Fermi–Dirac distribution.
Entropy 25 01571 g005
Figure 6. The spectra of p produced in p p collisions at s = (a) 6.3, (b) 7.7, (c) 8.8, (d) 12.3, and (e) 17.3 GeV at different y. The symbols represent the experimental data measured by the NA61/SHINE Collaboration [29] and the curves are the fitting results from the Fermi–Dirac distribution.
Figure 6. The spectra of p produced in p p collisions at s = (a) 6.3, (b) 7.7, (c) 8.8, (d) 12.3, and (e) 17.3 GeV at different y. The symbols represent the experimental data measured by the NA61/SHINE Collaboration [29] and the curves are the fitting results from the Fermi–Dirac distribution.
Entropy 25 01571 g006
Figure 7. Dependence of T on y at different s from the spectra of (a) π , (b) π + , (c) K , (d) K + , (e) p ¯ , and (f) p.
Figure 7. Dependence of T on y at different s from the spectra of (a) π , (b) π + , (c) K , (d) K + , (e) p ¯ , and (f) p.
Entropy 25 01571 g007
Figure 8. Dependence of T 0 on y at different s from the spectra of (a) π , (b) π + , (c) K , (d) K + , (e) p ¯ , and (f) p.
Figure 8. Dependence of T 0 on y at different s from the spectra of (a) π , (b) π + , (c) K , (d) K + , (e) p ¯ , and (f) p.
Entropy 25 01571 g008
Figure 9. Dependence of T i on y at different s from the spectra of (a) π , (b) π + , (c) K , (d) K + , (e) p ¯ , and (f) p.
Figure 9. Dependence of T i on y at different s from the spectra of (a) π , (b) π + , (c) K , (d) K + , (e) p ¯ , and (f) p.
Entropy 25 01571 g009
Figure 10. Dependence of V on y at different s from the spectra of (a) π , (b) π + , (c) K , (d) K + , (e) p ¯ , and (f) p.
Figure 10. Dependence of V on y at different s from the spectra of (a) π , (b) π + , (c) K , (d) K + , (e) p ¯ , and (f) p.
Entropy 25 01571 g010
Figure 11. Dependence of T on s at different y from the spectra of (a) π , (b) π + , (c) K , (d) K + , (e) p ¯ , and (f) p.
Figure 11. Dependence of T on s at different y from the spectra of (a) π , (b) π + , (c) K , (d) K + , (e) p ¯ , and (f) p.
Entropy 25 01571 g011
Figure 12. Dependence of T 0 on s at different y from the spectra of (a) π , (b) π + , (c) K , (d) K + , (e) p ¯ , and (f) p.
Figure 12. Dependence of T 0 on s at different y from the spectra of (a) π , (b) π + , (c) K , (d) K + , (e) p ¯ , and (f) p.
Entropy 25 01571 g012
Figure 13. Dependence of T i on s at different y from the spectra of (a) π , (b) π + , (c) K , (d) K + , (e) p ¯ , and (f) p.
Figure 13. Dependence of T i on s at different y from the spectra of (a) π , (b) π + , (c) K , (d) K + , (e) p ¯ , and (f) p.
Entropy 25 01571 g013
Figure 14. Dependence of V on s at different y from the spectra of (a) π , (b) π + , (c) K , (d) K + , (e) p ¯ , and (f) p.
Figure 14. Dependence of V on s at different y from the spectra of (a) π , (b) π + , (c) K , (d) K + , (e) p ¯ , and (f) p.
Entropy 25 01571 g014
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Yang, P.-P.; Liu, F.-H.; Olimov, K.K. Rapidity and Energy Dependencies of Temperatures and Volume Extracted from Identified Charged Hadron Spectra in Proton–Proton Collisions at a Super Proton Synchrotron (SPS). Entropy 2023, 25, 1571. https://doi.org/10.3390/e25121571

AMA Style

Yang P-P, Liu F-H, Olimov KK. Rapidity and Energy Dependencies of Temperatures and Volume Extracted from Identified Charged Hadron Spectra in Proton–Proton Collisions at a Super Proton Synchrotron (SPS). Entropy. 2023; 25(12):1571. https://doi.org/10.3390/e25121571

Chicago/Turabian Style

Yang, Pei-Pin, Fu-Hu Liu, and Khusniddin K. Olimov. 2023. "Rapidity and Energy Dependencies of Temperatures and Volume Extracted from Identified Charged Hadron Spectra in Proton–Proton Collisions at a Super Proton Synchrotron (SPS)" Entropy 25, no. 12: 1571. https://doi.org/10.3390/e25121571

APA Style

Yang, P. -P., Liu, F. -H., & Olimov, K. K. (2023). Rapidity and Energy Dependencies of Temperatures and Volume Extracted from Identified Charged Hadron Spectra in Proton–Proton Collisions at a Super Proton Synchrotron (SPS). Entropy, 25(12), 1571. https://doi.org/10.3390/e25121571

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop