Properties of Chemical and Kinetic FreezeOuts in HighEnergy Nuclear Collisions
View this Special IssueResearch Article  Open Access
RuiFang Si, HuiLing Li, FuHu Liu, "Comparing Standard Distribution and Its Tsallis Form of Transverse Momenta in High Energy Collisions", Advances in High Energy Physics, vol. 2018, Article ID 7895967, 12 pages, 2018. https://doi.org/10.1155/2018/7895967
Comparing Standard Distribution and Its Tsallis Form of Transverse Momenta in High Energy Collisions
Abstract
The experimental (simulated) transverse momentum spectra of negatively charged pions produced at midrapidity in central nucleusnucleus collisions at the HeavyIon Synchrotron (SIS), Relativistic HeavyIon Collider (RHIC), and Large Hadron Collider (LHC) energies obtained by different collaborations are selected by us to investigate, where a few simulated data are taken from the results of FOPI Collaboration which uses the IQMD transport code based on Quantum Molecular Dynamics. A twocomponent standard distribution and the Tsallis form of standard distribution are used to fit these data in the framework of a multisource thermal model. The excitation functions of main parameters in the two distributions are analyzed. In particular, the effective temperatures extracted from the twocomponent standard distribution and the Tsallis form of standard distribution are obtained, and the relation between the two types of effective temperatures is studied.
1. Introduction
High energy heavyion (nucleusnucleus) collisions are an important method to simulate and study the big bang in the early universe, properties of new matter created in extreme conditions, accompanying phenomena in the creation, and physics mechanisms of the creation. Some models based on the quantum chromodynamics (QCD) and/or thermal and statistical methods can be used to analyze the equation of state (EoS) at finite temperature and density, properties of chemical and kinetic freezeouts in collision process, distribution laws of different particles in final state, and universality of hadroproduction in different systems [1–5]. The properties of nuclear matter and its phase transition to quarkgluon plasma (QGP) at high temperature and density can be obtained. With the developments in the methodologies of experimental techniques and theoretical studies, the collision energy per nucleon pair in the centerofmass system increases from high energy range which has a few to several hundred GeV to ultrahigh energy range which has presently a few to over ten TeV.
The temperature and density described the EoS showing that the new matter created in high and ultrahigh energy ranges is not similar to the ideal gaslike state of quarks and gluons expected by early theoretical models. Instead, the effects of strong dynamical coupling, longrange interactions, local memory, and others appear in the interior of interacting system. The rapid evolution of interacting system and the indirect measurements of some observable quantities result in that one can use the statistical method to study the distribution properties of some observable quantities such as (pseudo) rapidity, (transverse) momentum, (transverse) energy, azimuthal angle, elliptic flow, multiplicity, and others of finalstate fragments and particles [1–5]. Thus, some quantitative or qualitative results related to the properties of interacting system and particle production can be observed.
As the quantities which can be early measured in experiments, that is, the socalled “the first day” measureable quantities, the rapidity and transverse momentum distributions attract wide attentions due to their carryovers on the information of longitudinal extension and transverse expansion of the emission source in interacting system. With the increasing collision energy, the rapidity distribution range extends from a few rapidity units to over ten rapidity units, and the transverse momentum distribution range increases from 0 until a few GeV/c to 0 until over hundred GeV/c. Different functions and methods are used by different researchers to describe rapidity and transverse momentum distributions as well as other distributions which can be measured in experiments [1–5]. Based on a multisource thermal model [6–9], the rapidity and transverse momentum distributions obtained in experiments at different collision energies are studied by us in terms of twocylinder, Rayleigh, Boltzmann, Tsallis, and other distributions. In particular, comparing with rapidity distribution, transverse momentum distribution contains more abundant information and attracts wider attentions. Although one has Monte Carlo and other indirect methods to describe transverse momentum distributions, analytical functions are more expected to use.
Because of the same transverse momentum distribution being described by different functions to obtain values of different parameters, possible relations existing among different parameters can be studied. In this paper, based on the multisource thermal model [6–9], the standard distribution (Boltzmann, FermiDirac, and BoseEinstein distributions) and its Tsallis form are used to describe the transverse momentum distribution of finalstate particles produced in high energy nucleusnucleus collisions. The excitation functions of effective temperatures obtained by the two distributions are extracted and the relation between the two effective temperatures is studied.
The rest part of this paper is structured as follows. A brief description of the model and method is presented in Section 2. Results on comparisons with experimental (simulated) data and discussion are given in Section 3. Finally, we summarize our main observations and conclusions in Section 4.
2. The Model and Method
According to the multisource model [6–9], a few emission sources of produced particles are assumed to form in interacting system due to different reaction mechanisms and/or data examples. For each emission source, the thermal model or other similar models and distributions can be used to perform calculation on the production of particles. The potential models include [10], but are not limited to, ideal gaslike model, ideal hydrodynamic model, and viscous hydrodynamic model. In these models, the relativistic effect has to be particularly considered, and the quantum effect can be usually neglected. If we study in detail the interacting system and finalstate particles, both the relativistic and quantum effects have to be considered.
In the middle stage of collision process, the interacting system and emission sources in it can be regarded as to stay at the hydrodynamic state. After the stage of chemical freezeout, in particular after the stage of kinetic freezeout, the interacting system and emission sources in it should stay at the gaslike state. Otherwise, it is difficult to understand the kinetic information of singular particle measured in experiments. What had happened during the phase transition from the liquidlike state at the middle stage to the gaslike state at the final stage and why is beyond the focus of the present work. We shall not discuss this issue here.
According to the ideal gas model with the relativistic and quantum effects, the particle spectra can be described by the standard distribution. The number of particles is [11]where is the degeneracy factor, is the volume, is the momentum, is the energy, is the rest mass, is the chemical potential, and is the effective temperature; , , and correspond to the Boltzmann, FermiDirac, and BoseEinstein statistics, respectively; corresponds to plus +, and corresponds to minus −. The invariant momentum distribution of particles isThe normalized probability density distribution of particle momenta can be written aswhere is the normalized constant in the standard probability density distribution of momenta. It is related to the selection of parameters.
The normalized joint probability density distribution of particle rapidities and transverse momenta iswhere corresponds to plus + and corresponds to minus −. The normalized probability density distribution of particle rapidities is then written to bewhere denotes the maximum transverse momentum. This rapidity distribution is only for an emission source. In the case of considering multiple sources, we have to consider sources distribution in the rapidity space [2–4, 12–16]. This issue is beyond the focus of the present work, and we shall not discuss it anymore. The normalized probability density distribution of particle transverse momenta is written to bewhere and denote the maximum and minimum rapidities, respectively.
It should be noted that, in the above formulas, although the same symbol is used to represent the normalized constants in different formulas, these constants may be different from each other. In the case of considering multisource emission, we have to use the multicomponent distribution to describe the transverse momentum distribution of finalstate particles. If emission sources are considered, we havewhere denotes the normalized constant for the th component in components, denotes the contribution fraction of the th component in finalstate distribution, and denotes the effective temperature corresponding to the th component. There are temperature fluctuations among different components. In the case of considering multisource emission, we have the effective temperature of interacting system as . Generally, two or three emission sources are enough to describe the experimental data obtained in soft excitation process. That is, 2 or 3 in most cases.
If we consider the Tsallis form of standard distribution, the number of particles is [11, 17]where is an entropy index which characterizes the departing degree of the interacting system from the equilibrium state. Generally, we have ; if , the system stays in the equilibrium state. is the effective temperature. Other symbols have the same meanings as (1). The invariant momentum distribution of particles isThe normalized probability density distribution of particle momenta is
The normalized joint probability density distribution of particle rapidities and transverse momenta isThen, the normalized probability density distribution of particle rapidities isThe normalized probability density distribution of particle transverse momenta isIn the above formulas, although the same symbol is used to represent the normalized constants in different formulas, these constants may be different from each other. As discussed in [17], the Tsallis form has at least four types of function representations, though we choose only one that contains after and the index . We do not need to consider a multisource for the Tsallis form due to it covering a two or threecomponent standard distribution, and the two or threecomponent standard distribution describes well the transverse momentum spectrum of particles produced in soft excitation process.
It should be noted again that the above multicomponent (two or threecomponent) standard distribution and the Tsallis form of standard distribution can describe only the transverse momentum spectrum of particles produced in soft excitation process. The transverse momentum spectrum produced in soft excitation process covers a narrow range. For the transverse momentum spectrum covering a wide range, we have to consider the contribution of hard scattering process. According to the QCD calculus [18–20], we have an inverse powerlawto describe the transverse momentum spectrum produced in hard scattering process, where and are free parameters and is the normalized constant which is related to the free parameters. It is obvious that a twocomponent function is needed for a wide transverse momentum spectrum. The first component is the multicomponent (two or threecomponent) standard distribution or Tsallis form which describes the soft process, and the second component is the inverse powerlaw which describes the hard process. The application of the inverse powerlaw is beyond the focus of the present work. We shall not discuss it anymore.
In the above discussions, to obtain chemical potential of a given particle, the chemical freezeout temperature of the emission source is needed to know first of all. In the case of assuming the same chemical freezeout moment, the emission source has the sole . According to [21, 22], there is a relation among , the yield and mass of the first particle, the yield and mass of the second particle, and the ratio . We havewhere denote fermion and boson, respectively. If the fermion and boson are not needed to distinguish each other, we have . This results in a simple expression for (15); that is, .
In the framework of a statistical thermal model of noninteracting gas particles with the assumption of standard MaxwellBoltzmann statistics, there is an empirical expression for the chemical freezeout temperature [23–26],where denotes the energy per nucleon pair in the centerofmass system. Both the units of and are in GeV. The limiting value of is GeV.
In the framework of a thermal model with standard distribution, the chemical potentials of some particles can be obtained from the ratios of negatively to positively charged particles. According to [27], we havewhere the symbol of a given particle is used for its yield for the purpose of simplicity. Further, the chemical potentials of the mentioned particles areEmpirically, the chemical potential for baryon is [23–26]which is also obtained in the framework of a statistical thermal model of noninteracting gas particles with the assumption of standard MaxwellBoltzmann statistics, where both the units of and are in GeV.
We would like to point out that (16) and (23) should be modified in the framework of generalized nonextensive statistics when we use the Tsallis form of standard distribution. At the same time, (17)–(22) should be generalized within an analysis with the Tsallis form. To modify (16)–(23) is beyond our focus and ability. We shall not discuss these modifications here. Instead, as an approximate treatment, we use and obtained within an analysis with the standard distribution as those within the Tsallis form. In fact, the absolute value of is very small, and its effect on the transverse momentum spectra can be neglected. Therefore, this approximate treatment is acceptable.
It should be noted once more that, as mentioned in the above discussions, what we extract from the multicomponent standard distribution or the Tsallis form of standard distribution is the effective temperature, but not the real temperature of emission source. Generally, the transverse momentum spectrum contains both the contributions of thermal motion and flow effect. The real temperature is only a reflection of purely thermal motion, and the flow effect should not be included in it. As for the methods to obtain the real temperature by disengaging the contributions of thermal motion and flow effect, we can use the blastwave model based on the Boltzmann distribution [28–30], the blastwave model based on the Tsallis distribution [31], the improved Tsallis distribution [32, 33], some alternative methods [21, 29, 34–36], and others [37–40]. These methods themselves are beyond the focus of the present work. We shall not discuss them anymore.
3. Results and Discussion
The transverse momentum spectra of negatively charged pions produced in midrapidity range in 2.24 and 2.52 GeV central goldgold (AuAu) collisions [41] measured (simulated) by the FOPI Collaboration at the HeavyIon Synchrotron (SIS), 11.5 [42], 62.4, 130, and 200 GeV central AuAu collisions [29] measured by the STAR Collaboration at the Relativistic HeavyIon Collider (RHIC), 22.5 GeV central coppercopper (CuCu) [43] and 200 GeV central AuAu collisions [27] measured by the PHENIX Collaboration at the RHIC, and 2.76 TeV central leadlead (PbPb) collisions [44] measured by the ALICE Collaboration at the Large Hadron Collider (LHC) are selected to investigate. Among them, the results of FOPI Collaboration are given in Figure 1 with the simulated data (the last eight circles) of the IQMD transport code [45] which is based on Quantum Molecular Dynamics [46]. To avoid confusion, most results of the STAR Collaboration are given in Figure 2, and the results corresponding to 11.5 GeV are given in Figure 3. The results of PHENIX and ALICE Collaborations are given in Figures 3 and 4, respectively. In each figure, the symbols represent the experimental (simulated) data scaled by different amounts in some cases. The collision energy and type, centrality and midrapidity ranges, and scaled amount if not 1 are marked in the panel. The dashed and solid curves denote the results fitted by the twocomponent standard distribution and the Tsallis form of standard distribution. The values of parameters, , and degree of freedom (dof) are listed in Table 1 ordered by the energy from low to high. In particular, is the average weighted by the fractions of different components, is obtained by (16) and (22), and the values of in (22) at different energies are obtained from [47]. As a preliminary result, the values of for the first and second standard distributions and the Tsallis form are assumed to be the same. In the fitting, the method of least square is used to obtain the best parameter values. One can see that the twocomponent standard distribution and the Tsallis form of standard distribution describe approximately the transverse momentum spectra of negatively charged pions produced in central nucleusnucleus collisions in the energy range from SIS to LHC.

To study the excitation functions of free parameters, that is, the dependence of free parameters on collision energy, the relations (), (), , and are presented in Figures 5–8, respectively. The symbols and error bars in the figures denote the values of free parameters and their errors. Both the values of free parameters and their errors are taken from Table 1. The lines in Figures 5 and 6 are obtained by the method of least square. These lines can be described by linear functions , where the slope and intercept are listed in Table 1 and the unit of is in GeV. One can see that the four effective temperatures , , , and increase linearly with increase of . In particular, the relation between and can be obtained to be due to Table 1, which shows a linear relation between and . With increase of , has a minimum at about 10 GeV, and increases primitively and saturates at about 10 GeV.
Our results show some interesting features. Actually, one could as well say that there is no difference in the particle production in central nucleusnucleus collisions from a few GeV to a few TeV. This in some sense echoes recent studies of Sarkisyan et al. [1, 48]. In addition, our recent study shows that the same or similar fits to be good for protonproton collisions [40], though the parameter values in protonproton collisions are closer to those in peripheral nucleusnucleus collisions when comparing with central nucleusnucleus collisions. This suggests universality in particle production, as it is obtained in recent and previous studies of Sarkisyan et al. [1, 12, 48–50], but now for transverse momentum distribution as well. On the other hand, the multiplicity and transverse momentum distributions observed in different data samples can be uniformly fitted by multicomponent Erlang distribution [8, 51, 52], which also show the universality in particle production. Indeed, the universality in particle production exists not only in mean multiplicity and pseudorapidity density but also in multiplicity and transverse momentum distributions in some conditions.
Our observation that has a minimum at about 10 GeV and increases primitively and saturates at about 10 GeV is in agreement with recent work of Cleymans [53] in which the energy region GeV for heavyion collisions is indicated to be an interesting one. In fact, in this energy region, the final state has the highest net baryon density, and a transition from a baryon dominated to a meson dominated final state takes place. At the same time, ratios of strange particles to mesons show obviously maxima in this energy region [53]. At a slightly smaller energy (about 6~8 GeV), other works show some extremes or saturation in excitation functions of parameters. These parameters include, but are not limited to, the specific reduced curvature of netproton rapidity distribution [54–56], chemical freezeout temperature [57, 58], mean transverse mass minus rest mass [57], yield ratios of positive kaons to pions [57–59], squared speedofsound [60], string tension in Schwinger mechanism [61], width and fraction of fragmentation source [62], and width ratios of experimental negative pion rapidity distribution to Landau hydrodynamic model prediction [58].
In the above analyses, for a not too wide transverse momentum spectrum, a standard distribution is usually not enough to describe the spectrum. Generally, we need a twocomponent standard distribution to describe the not too wide spectrum. It is expected that, in the case of studying a wider transverse momentum spectrum, we need a threecomponent standard distribution to describe the wider spectrum. If a set of experimental data is described by the two or threecomponent standard distribution, it is also described by the Tsallis form of standard distribution [11]. If the two or threecomponent standard distribution describes a temperature fluctuation between two or among three emission sources, the Tsallis form of standard distribution describes a degree of nonequilibrium. The degree of nonequilibrium is characterized by the entropy index . A larger corresponds to a farther nonequilibrium among different emission sources. One can see from Table 1 that in most cases, which renders an approximate equilibrium among different emission sources or the whole interacting system stays in an approximate equilibrium state.
Both the two or threecomponent standard distribution and the Tsallis form of standard distribution describe only the results of soft excitation process. For the soft process, the particle spectrum appears with the characteristics of thermal emission phenomenon. Although the standard distribution describes the characteristics of thermal emission, some nonthermal emissions also obey the standard distribution. Even if the Tsallis form has less connection with thermal emission, they are relative due to the standard distribution. In the case of studying a very wide transverse momentum spectrum, for example, for a width of more than 5 GeV/c, to consider only the contribution of soft process is not enough in description of experimental data. To describe a wider transverse momentum spectrum, we have to consider simultaneously the contribution of hard scattering process. As mentioned in Section 2, according to the QCD calculus [18–20], the hard process can be described by the inverse powerlaw. Because of the hard process having no connection with the thermal emission, it does not affect the extraction of temperature parameter. In the case of extracting only temperature parameter, a too wide transverse momentum spectrum is not needed.
In the above analyses, the temperature extracted by us is in fact the effective temperature . It is neither the temperature at the kinetic freezeout nor the temperature at the chemical freezeout of the emission source or interacting system. Generally, can be extracted from the transverse momentum spectra, and can be extracted from the ratios of different types of particles. However, the temperature extracted from the transverse momentum spectra is not surely due to the contribution of flow effect. How to get rid of the contribution of flow effect is a question that is worth discussing. In the blastwave model [28–31], the mean transverse flow velocity is introduced. Thus, and can be simultaneously obtained based on the analysis of transverse momentum spectrum by the model. In addition, by using the standard distribution and its Tsallis form to analyze the transverse momentum spectra of different particles, we can obtain the linear relation between and . The intercept in the linear relation is regarded as [21, 29, 34–36]. We can also obtain the linear relation between and mean moving mass (mean energy). The slope in the linear relation is regarded as [37–40].
Generally speaking, the twocomponent standard distribution and the Tsallis form of standard distribution are the same in essentials while differing in minor points in the behaviors in the figures. The standard distribution corresponds to the classical statistical system which has shortrange interactions in interior and nonmultifractal structure in boundary. Some extensive thermodynamic quantities such as energy, momentum, internal energy, and entropy are linearly related to the system size and particle number. These quantities obey simply additive property. The statistical method and the microscopic description of system are adaptive. The entropy function is a power tool to study the microscopic dynamics of system under the macroscopic condition by describing the occupation number of phase spaces of the system. The Tsallis form breaks through the limitation of classical statistics by using the entropy index . The complex system with longrange interactions, local memory effect, strong dynamic correlation, fractal or multifractal occupation in phase space, and others can be described by the Tsallis form. The Tsallis form also causes the classical extensive quantities not to obey the simple additive property. Instead, the coupled item appears in the quantities and the nonextensive statistical effects are formed in the transverse and longitudinal dynamics [63–73].
In the above discussions, one can see that the two or threecomponent standard distribution can be described by the Tsallis form of the standard distribution. It does not mean that the single standard distribution cannot be described by the Tsallis form. In fact, by using a lower temperature and an entropy index that is closer to 1, the Tsallis form describes well the single standard distribution. The standard distribution is successfully replaced by the Tsallis form due to changing from 1 to a value that is greater than 1. This means that the interacting system changes from the classical and extensive statistical system to the nonextensive system, which is an essential change of the system properties. However, in some cases, the same set of experimental data can be described by both the (two or threecomponent) standard distribution obeying the extensive statistics and the Tsallis form obeying the nonextensive statistics. This means that in these cases there is no obvious boundary to distinguish extensive system and nonextensive system for a given interacting system. We have to examine which property is the main factor. Or, the interacting system in the present energy range stays in a transition gradation from extensive system to nonextensive system.
4. Conclusions
We summarize here our main observations and conclusions.
(a) The transverse momentum spectra of negatively charged pions produced in central nucleusnucleus collisions measured (simulated) in midrapidity range by different collaborations at the SIS, RHIC, and LHC are studied by the twocomponent standard distribution and the Tsallis form of standard distribution which are fitted into the frame of multisource thermal model. The two distributions describe approximately the experimental (simulated) data.
(b) The excitation functions of related parameters are analyzed. The four effective temperatures , , , and increase linearly with increase of . In particular, the relation between and can be obtained to be which shows a linear relation between and . With increase of , has a minimum at about 10 GeV, and increases primitively and saturates at about 10 GeV.
(c) There is no difference in the particle production in central nucleusnucleus collisions from a few GeV to a few TeV. Combining with other works, one can say that the same or similar fits are good for protonproton collisions. This suggests universality in particle production, as it is already obtained in mean multiplicity, pseudorapidity density, and multiplicity distribution, but now for transverse momentum distribution as well.
(d) The energy of GeV for heavyion collisions is indicated to have the highest net baryon density and the maximum ratios of strange particles to mesons and to take place a transition from the final state which has mainly baryons to the final state which has mainly mesons [53]. At a slightly smaller energy (about 6~8 GeV), other works show some extremes or saturation in excitation functions of some parameters [54–62]. These extremes and saturation are related to the search of soft point of equation of state.
(e) To be closer to the classical situation, the two or threecomponent standard distribution has an advantage over the Tsallis form of standard distribution due to similar statistics for the classical situation and standard distribution. However, the Tsallis form of standard distribution uses less parameter than the two or threecomponent standard distribution. If the two or threecomponent standard distribution describes a temperature fluctuation between two or among three sources, the Tsallis form of standard distribution describes a degree of nonequilibrium.
(f) In the considered energy range, different emission sources stay in an approximate equilibrium state or the whole interacting system stays in an approximate equilibrium state. There is no obvious boundary to distinguish extensive system and nonextensive system for a given interacting system. The interacting system stays in a transition gradation from extensive system to nonextensive system. To obtain only the kinetic freezeout temperature, we would rather use the two or threecomponent standard distribution due to it being closer to the classical situation.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
Comments on the manuscript and relevant communications from Edward K. G. Sarkisyan and YaHui Chen are highly acknowledged. This work was supported by the National Natural Science Foundation of China under Grants nos. 11575103 and 11747319, the Shanxi Provincial Natural Science Foundation under Grant no. 201701D121005, and the Fund for Shanxi “1331 Project” Key Subjects Construction.
References
 E. K. G. Sarkisyan, A. N. Mishra, R. Sahoo, and A. S. Sakharov, “Multihadron production dynamics exploring the energy balance in hadronic and nuclear collisions,” Physical Review D, vol. 93, no. 5, Article ID 054046, 2016. View at: Publisher Site  Google Scholar
 I. Bautista, J. G. Milhano, C. Pajares, and J. D. de Deus, “Multiplicity in pp and AA collisions: The same power law from energymomentum constraints in string production,” Physics Letters B, vol. 715, no. 13, pp. 230–233, 2012. View at: Publisher Site  Google Scholar
 Z. J. Jiang, Y. Zhang, H. L. Zhang, and H. P. Deng, “A description of the pseudorapidity distributions in heavy ion collisions at RHIC and LHC energies,” Nuclear Physics A, vol. 941, pp. 188–200, 2015. View at: Publisher Site  Google Scholar
 G. Wolschin, “Ultraviolet energy dependence of particle production sources in relativistic heavyion collisions,” Physical Review C, vol. 91, no. 1, Article ID 014905, 2015. View at: Publisher Site  Google Scholar
 M. SanchisLozano and E. SarkisyanGrinbaum, “Ridge effect and threeparticle correlations,” Physical Review D, vol. 96, no. 7, Article ID 074012, 2017. View at: Publisher Site  Google Scholar
 F.H. Liu, “Particle production in AuAu collisions at RHIC energies,” Physics Letters B, vol. 583, no. 12, pp. 68–72, 2004. View at: Publisher Site  Google Scholar
 F.H. Liu, “Dependence of charged particle pseudorapidity distributions on centrality and energy in p(d)A collisions at high energies,” Physical Review C, vol. 78, no. 1, Article ID 014902, 2008. View at: Publisher Site  Google Scholar
 F.H. Liu, “Unified description of multiplicity distributions of finalstate particles produced in collisions at high energies,” Nuclear Physics A, vol. 810, no. 1–4, pp. 159–172, 2008. View at: Publisher Site  Google Scholar
 F.H. Liu, “Charged particle production in d+Au collisions at GeV,” Physical Review C, vol. 69, no. 6, Article ID 067901, 2004. View at: Google Scholar
 S. Abreu, S. V. Akkelin, J. Alam et al., “Heavyion collisions at the LHC Last call for predictions,” Journal of Physics G, vol. 35, no. 5, Article ID 054001, 2008. View at: Google Scholar
 J. Cleymans and D. Worku, “Relativistic thermodynamics: transverse momentum distributions in highenergy physics,” The European Physical Journal A, vol. 48, Article ID 160, 2012. View at: Publisher Site  Google Scholar
 A. N. Mishra, R. Sahoo, E. K. G. Sarkisyan, and A. S. Sakharov, “Effectiveenergy budget in multiparticle production in nuclear collisions,” The European Physical Journal C, vol. 74, Article ID 3147, 2014. View at: Publisher Site  Google Scholar
 G. Wolschin, “Particle production sources at LHC energies,” Journal of Physics G: Nuclear and Particle Physics, vol. 40, no. 4, Article ID 045104, 2013. View at: Publisher Site  Google Scholar
 I. Bautista, C. Pajares, J. G. Milhano, and J. D. De Deus, “Rapidity dependence of particle densities in pp and AA collisions,” Physical Review C, vol. 86, no. 3, Article ID 034909, 2012. View at: Publisher Site  Google Scholar
 Z. J. Jiang and H. L. Zhang, “Pseudorapidity distributions of charged particles produced in pp collisions at centerofmass energies from 23.6 GeV to 900 GeV,” Modern Physics Letters A, vol. 29, no. 27, Article ID 1450130, 2014. View at: Publisher Site  Google Scholar
 Z. J. Jiang, H. L. Zhang, J. Wang, and K. Ma, “The evolutiondominated hydrodynamic model and the pseudorapidity distributions in high energy physics,” Advances in High Energy Physics, vol. 2014, Article ID 248360, 10 pages, 2014. View at: Publisher Site  Google Scholar
 H. Zheng and L. Zhu, “Comparing the Tsallis distribution with and without thermodynamical description in p+p collisions,” Advances in High Energy Physics, vol. 2016, Article ID 9632126, 10 pages, 2016. View at: Publisher Site  Google Scholar
 R. Odorico, “Does a transverse energy trigger actually trigger on largeP_{T} jets?” Physics Letters B, vol. 118, no. 13, pp. 151–154, 1982. View at: Publisher Site  Google Scholar
 G. Arnison, A. Astbury, B. Aubert (UA1 Collaboration) et al., “Transverse momentum spectra for charged particles at the CERN protonantiproton collider,” Physics Letters B, vol. 118, no. 13, pp. 167–172, 1982. View at: Google Scholar
 T. Mizoguchi, M. Biyajima, and N. Suzuki, “Analyses of whole transverse momentum distributions in and collisions by using a modified version of Hagedorn’s formula,” International Journal of Modern Physics A, vol. 32, no. 11, Article ID 1750057, 2017. View at: Publisher Site  Google Scholar
 S. Takeuchi, K. Murase, T. Hirano, P. Huovinen, and Y. Nara, “Effects of hadronic rescattering on multistrange hadrons in highenergy nuclear collisions,” Physical Review C, vol. 92, no. 4, Article ID 044907, 2015. View at: Publisher Site  Google Scholar
 C. Gale, Y. Hidaka, S. Jeon et al., “Production and elliptic flow of dileptons and photons in a matrix model of the quarkgluon plasma,” Physical Review Letters, vol. 114, no. 7, Article ID 072301, 2015. View at: Publisher Site  Google Scholar
 A. Andronic, P. BraunMunzinger, and J. Stachel, “Thermal hadron production in relativistic nuclear collisions,” Acta Physica Polonica B, vol. 40, pp. 1005–1012, 2009. View at: Google Scholar
 A. Andronic, P. BraunMunzinger, and J. Stachel, “The Horn, the hadron mass spectrum and the QCD phase diagram: the statistical model of hadron production in central nucleusnucleus collisions,” Nuclear Physics A, vol. 834, pp. 237c–240c, 2010. View at: Publisher Site  Google Scholar
 A. Andronic, P. BraunMunzinger, and J. Stachel, “Hadron production in central nucleusnucleus collisions at chemical freezeout,” Nuclear Physics A, vol. 772, no. 34, pp. 167–199, 2006. View at: Publisher Site  Google Scholar
 J. Cleymans, H. Oeschler, K. Redlich, and S. Wheaton, “Comparison of chemical freezeout criteria in heavyion collisions,” Physical Review C, vol. 73, no. 3, Article ID 034905, 2006. View at: Publisher Site  Google Scholar
 S. S. Adler, S. Afanasiev, C. Aidala (PHENIX Collaboration) et al., “Identified charged particles spectra and yields in Au+Au collisions at 200 GeV,” Physical Review C, vol. 69, Article ID 034909, 2004. View at: Google Scholar
 E. Schnedermann, J. Sollfrank, and U. Heinz, “Thermal phenomenology of hadrons from 200A GeV S+S collisions,” Physical Review C, vol. 48, no. 5, pp. 2462–2475, 1993. View at: Publisher Site  Google Scholar
 B. I. Abelev, M. M. Aggarwal, Z. Ahammed (STAR Collaboration) et al., “Systematic measurements of identified particles spectra in pp, d+Au, and Au+Au collisions at the STAR detector,” Physical Review C, vol. 79, no. 3, Article ID 034909, 2009. View at: Google Scholar
 B. I. Abelev, M. M. Aggarwal, Z. Ahammed (STAR Collaboration) et al., “Identified particle production, azimuthal anisotropy, and interferometry measurements in Au+Au collisions at 9.2 GeV,” Physical Review C, vol. 81, no. 2, Article ID 024911, 2010. View at: Google Scholar
 Z. Tang, Y. Xu, L. Ruan, G. van Buren, F. Wang, and Z. Xu, “Spectra and radial flow in relativistic heavy ion collisions with Tsallis statistics in a blastwave description,” Physical Review C, vol. 79, no. 5, Article ID 051901, 2009. View at: Publisher Site  Google Scholar
 D. Thakur, S. Tripathy, P. Garg, R. Sahoo, and J. Cleymans, “Indication of a differential freezeout in protonproton and heavyion collisions at RHIC and LHC energies,” Advances in High Energy Physics, vol. 2016, Article ID 4149352, 13 pages, 2016. View at: Publisher Site  Google Scholar
 T. Bhattacharyya, J. Cleymans, A. Khuntia, P. Pareek, and R. Sahoo, “Radial flow in nonextensive thermodynamics and study of particle spectra at LHC in the limit of small (q1),” The European Physical Journal A, vol. 52, Article ID 30, 2016. View at: Publisher Site  Google Scholar
 H. Heiselberg and A.M. Levy, “Elliptic flow and HanburyBrownTwiss correlations in noncentral nuclear collisions,” Physical Review C, vol. 59, no. 5, pp. 2716–2727, 1999. View at: Publisher Site  Google Scholar
 U. W. Heinz, Concepts of heavyion physics, Lecture Notes for Lectures Presented at the 2nd CERNLatinAmerican School of HighEnergy Physics, 114 June 2003, San Miguel Regla, Mexico, arXiv:hepph/0407360, 2004.
 R. Russo, Measurement of D^{+} meson production in pPb collisions with the ALICE detector [Ph.D. thesis], Universita degli Studi di Torino, Torino, Italy, 2015.
 H.R. Wei, F.H. Liu, and R. A. Lacey, “Kinetic freezeout temperature and flow velocity extracted from transverse momentum spectra of finalstate light flavor particles produced in collisions at RHIC and LHC,” The European Physical Journal A, vol. 52, Article ID 102, 2016. View at: Publisher Site  Google Scholar
 H.R. Wei, F.H. Liu, and R. A. Lacey, “Disentangling random thermal motion of particles and collective expansion of source from transverse momentum spectra in high energy collisions,” Journal of Physics G: Nuclear and Particle Physics, vol. 43, no. 12, Article ID 125102, 2016. View at: Publisher Site  Google Scholar
 H.L. Lao, H.R. Wei, F.H. Liu, and R. A. Lacey, “An evidence of massdependent differential kinetic freezeout scenario observed in PbPb collisions at 2.76 TeV,” The European Physical Journal A, vol. 52, Article ID 203, 2016. View at: Publisher Site  Google Scholar
 H.L. Lao, F.H. Liu, B.C. Li, M.Y. Duan, and R. A. Lacey, Examining the model dependence of extracting the kinetic freezeout temperature and transverse flow velocity in small collision system, arXiv:1708.07749, 2017.
 W. Reisdorf, M. Stockmeier, A. Andronic (FOPI Collaboration) et al., “Systematics of pion emission in heavy ion collisions in the 1A GeV regime,” Nuclear Physics A, vol. 781, no. 34, pp. 459–508, 2007. View at: Publisher Site  Google Scholar
 S. Das (for the STAR Collaboration), “Centrality dependence of freezeout parameters from the beam energy scan at STAR,” Nuclear Physics A, vol. 904–905, pp. 891c–894c, 2013. View at: Publisher Site  Google Scholar
 J. T. Mitchell (for the PHENIX Collaboration), “The PHENIX potential in the search for the QCD critical point, Talk given at the 3rd International Workshop on the Critical Point and Onset of Deconfinement, Florence, Italy, July 2006, arXiv:nuclex/0701079, 2006”. View at: Google Scholar
 M. Floris (for the ALICE Collaboration), “Identified particles in pp and PbPb collisions at LHC energies with the ALICE detector,” Journal of Physics G: Nuclear and Particle Physics, vol. 38, no. 12, Article ID 124025, 2011. View at: Publisher Site  Google Scholar
 C. Hartnack, R. K. Puri, J. Aichelin et al., “Modelling the manybody dynamics of heavy ion collisions: Present status and future perspective,” The European Physical Journal A, vol. 1, pp. 151–169, 1998. View at: Google Scholar
 J. Aichelin, “‘Quantum’ molecular dynamics−a dynamical microscopic nbody approach to investigate fragment formation and the nuclear equation of state in heavy ion collisions,” Physics Reports, vol. 202, pp. 233–360, 1991. View at: Publisher Site  Google Scholar
 S. Das (for the STAR Collaboration), “Identified particle production and freezeout properties in heavyion collisions at RHIC beam energy scan program,” in Proceedings of the EPJ Web of Conference, vol. 90, Article ID 08007, 2015. View at: Google Scholar
 E. K. G. Sarkisyan, A. N. Mishra, R. Sahoo, and A. S. Sakharov, “Centrality dependence of midrapidity density from GeV to TeV heavyion collisions in the effectiveenergy universality picture of hadroproduction,” Physical Review D, vol. 94, no. 1, Article ID 011501, 2016. View at: Publisher Site  Google Scholar
 E. K. G. Sarkisyan and A. S. Sakharov, “Relating multihadron production in hadronic and nuclear collisions,” The European Physical Journal C, vol. 70, no. 3, pp. 533–541, 2010. View at: Publisher Site  Google Scholar
 E. K. G. Sarkisyan and A. S. Sakharov, “On similarities of bulk observables in nuclear and particle collisions,” in Proceedings of the AIP Conference, vol. 828, pp. 35–41, 2006. View at: Google Scholar
 F.H. Liu, Q.W. Lü, B.C. Li, and R. Bekmirzaev, “A description of the multiplicity distributions of nuclear fragments in hA and AA collisions at intermediate and high energies,” Chinese Journal of Physics, vol. 49, pp. 601–620, 2011. View at: Google Scholar
 F.H. Liu, Y.Q. Gao, T. Tian, and B.C. Li, “Unified description of transverse momentum spectrums contributed by soft and hard processes in highenergy nuclear collisions,” The European Physical Journal A, vol. 50, Article ID 94, 2014. View at: Publisher Site  Google Scholar
 J. Cleymans, The physics case for the 10 GeV energy region, arXiv:1711.02882 [hepph], 2017.
 Y. B. Ivanov and D. Blaschke, “Robustness of the baryonstopping signal for the onset of deconfinement in relativistic heavyion collisions,” Physical Review C, vol. 92, no. 2, Article ID 024916, 2015. View at: Publisher Site  Google Scholar
 Y. B. Ivanov, “Baryon stopping as a probe of deconfinement onset in relativistic heavyion collisions,” Physics Letters B, vol. 721, no. 13, pp. 123–130, 2013. View at: Publisher Site  Google Scholar
 Y. B. Ivanov, “Alternative scenarios of relativistic heavyion collisions. I. Baryon stopping,” Physical Review C, vol. 87, no. 6, Article ID 064904, 2013. View at: Publisher Site  Google Scholar
 L. Kumar (for the STAR Collaboration), “Identified hadron production from the RHIC beam energy scan,” Journal of Physics G: Nuclear and Particle Physics, vol. 38, no. 12, Article ID 124145, 2011. View at: Publisher Site  Google Scholar
 A. Rustamov, “The horn, kink and step, dale: from few GeV to few TeV,” Central European Journal of Physics, vol. 10, pp. 1267–1270, 2012. View at: Publisher Site  Google Scholar
 D. T. Larsen, “Energy dependence of hadron spectra and yields in p+p and ^{7}Be+^{9}Be collisions from the NA61/SHINE experiment at the CERN SPS,” Journal of Physics: Conference Series, vol. 668, no. 1, Article ID 012020, 2016. View at: Publisher Site  Google Scholar
 F.H. Liu, L.N. Gao, and R. A. Lacey, “Searching for Minimum in Dependence of Squared SpeedofSound on Collision Energy,” Advances in High Energy Physics, vol. 2016, Article ID 9467194, 9 pages, 2016. View at: Publisher Site  Google Scholar
 L.N. Gao, F.H. Liu, and R. A. Lacey, “Excitation functions of parameters in Erlang distribution, Schwinger mechanism, and Tsallis statistics in RHIC BES program,” The European Physical Journal A, vol. 52, Article ID 137, 2016. View at: Publisher Site  Google Scholar
 L.N. Gao, F.H. Liu, Y. Sun, Z. Sun, and R. A. Lacey, “Excitation functions of parameters extracted from threesource (net)proton rapidity distributions in AuAu and PbPb collisions over an energy range from AGS to RHIC,” The European Physical Journal A, vol. 53, Article ID 61, 2017. View at: Publisher Site  Google Scholar
 C. Tsallis, “Possible generalization of BoltzmannGibbs statistics,” Journal of Statistical Physics, vol. 52, no. 12, pp. 479–487, 1988. View at: Publisher Site  Google Scholar  MathSciNet
 C. Tsallis and E. P. Borges, “Nonextensive statistical mechanics Applications to nuclear and high energy physics,” in Proceedings of the 10th International Workshop on Multiparticle Production, N. Antoniou, Ed., 2003. View at: Google Scholar
 W. M. Alberico, P. Czerski, A. Lavagno, M. Nardi, and V. Somá, “Signals of nonextensive statistical mechanics in high energy nuclear collisions,” Physica A: Statistical Mechanics and its Applications, vol. 387, no. 23, pp. 467–475, 2008. View at: Publisher Site  Google Scholar
 C. Tsallis, “Nonadditive entropy: the concept and its use,” The European Physical Journal A, vol. 40, pp. 257–266, 2009. View at: Publisher Site  Google Scholar
 G. Wilk and Z. Włodarczyk, “Power laws in elementary and heavyion collisions : a story of fluctuations and nonextensivity?” The European Physical Journal A, vol. 40, pp. 299–312, 2009. View at: Publisher Site  Google Scholar
 W. M. Alberico and A. Lavagno, “Nonextensive statistical effects in highenergy collisions,” The European Physical Journal A, vol. 40, pp. 313–323, 2009. View at: Publisher Site  Google Scholar
 C. Tsallis, “An introduction to nonadditive entropies and a thermostatistical approach to inanimate and living matter,” Contemporary Physiscs, vol. 55, no. 3, pp. 179–197, 2014. View at: Publisher Site  Google Scholar
 L. Marques, J. Cleymans, and A. Deppman, “Description of highenergy pp collisions using Tsallis thermodynamics: transverse momentum and rapidity distributions,” Physical Review D, vol. 91, no. 5, Article ID 054025, 2015. View at: Publisher Site  Google Scholar
 C.Y. Wong, G. Wilk, L. J. L. Cirto, and C. Tsallis, “From QCDbased hardscattering to nonextensive statistical mechanical descriptions of transverse momentum spectra in highenergy and collisions,” Physical Review D, vol. 91, no. 11, Article ID 114027, 2015. View at: Publisher Site  Google Scholar
 M. Rybczyński, G. Wilk, and Z. Włodarczyk, “System size dependence of the logperiodic oscillations of transverse momentum spectra,” in Proceedings of the EPJ Web of Conference, vol. 90, 6 pages, 2015. View at: Publisher Site  Google Scholar
 H. Shababi and P. Pedram, “The minimal length uncertainty and the nonextensive thermodynamics,” International Journal of Theoretical Physics, vol. 55, no. 6, pp. 2813–2823, 2016. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2018 RuiFang Si et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. The publication of this article was funded by SCOAP^{3}.