{'fraction_non_alphanumeric': 0.0488766651202863, 'fraction_numerical': 0.016303267280800584, 'mean_word_length': 4.390069655295589, 'pattern_counts': {'":': 0, '<': 0, '<?xml version=': 0, '>': 0, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 5, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Energy dependent kinetic freeze-out temperature and transverse flow velocity in high energy collisions

Transverse momentum spectra of negative and positive pions produced at mid-(pseudo)rapidity in inelastic or non-single-diffractive proton-proton collisions and in central nucleus-nucleus collisions over an energy range from a few GeV to above 10 TeV are analyzed by a (two-component) blast-wave model with Boltzmann-Gibbs statistics and with Tsallis statistics respectively. The model results are in similarly well agreement with the experimental data measured by a few productive collaborations who work at the Heavy Ion Synchrotron (SIS), Super Proton Synchrotron (SPS), Relativistic Heavy Ion Collider (RHIC), and Large Hadron Collider (LHC), respectively. The energy dependent kinetic freeze-out temperature and transverse flow velocity are obtained and analyzed. Both the quantities have quick increase from the SIS to SPS, and slight increase or approximate invariability from the top RHIC to LHC. Around the energy bridge from the SPS to RHIC, the considered quantities in proton-proton collisions obtained by the blast-wave model with Boltzmann-Gibbs statistics show more complex energy dependent behavior comparing with the results in other three cases.


Introduction
The initial state, chemical freeze-out, and kinetic freeze-out are three main stages undergone by the interacting system in particle-particle, particle-nucleus, and nucleus-nucleus collisions at high energies. As the basic element (proton and neutron) in proton-nucleus and nucleus-nucleus collisions, proton-proton collisions display some similarities to proton-nucleus and nucleusnucleus collisions. In particular, proton-proton collisions and gold-gold (copper-copper, lead-lead) collisions have been studied in large-scale experiments at the Heavy Ion Synchrotron (SIS), Super Proton Synchrotron (SPS), Relativistic Heavy Ion Collider (RHIC), and Large Hadron Collider (LHC). A few productive collaborations have been reporting abundant data on the particle ratios, (pseudo)rapidity spectra, transverse momentum (mass) spectra, invariant mass spectra, anisotropic flows, nuclear modification factor, and so on. Some useful information related to particle production and system evolution can be extracted from these data.
"Temperature is surely one of the central concepts in thermodynamics and statistical mechanics" [1]. The temperature in subatomic physics describes the excitation degree of the interacting system. Due to the system evolution, the temperature is expected to change at different stages. Undoubtedly, the kinetic freeze-out temperature (the temperature at the kinetic freeze-out stage) is less than or equal to the chemical freeze-out temperature (the temperature at the chemical freeze-out stage) due to the fact that the kinetic freeze-out happens posteriorly or simultaneously compared with the chemical freeze-out. The temperature at the initial state is the highest due to the largest compression and density during the process of system evolution. After the initial state, the interacting system undergoes transverse expansion and longitudinal extension. At the same time, some particles are emitted during the evolution process, though most particles are emitted at the kinetic freezeout. The kinetic freeze-out temperature and concomitant transverse flow velocity can be extracted from the transverse momentum spectra of identified particles.
In order to extract the kinetic freeze-out temperature and transverse flow velocity and to study their dependence on energy, one can use different models to analyze the transverse momentum spectra. These models include, but are not limited to, the blast-wave model with Boltzmann-Gibbs statistics [2,3] or with Tsallis statistics [4][5][6] and the alternative method [7][8][9][10][11] with standard distribution or with Tsallis distribution, where the standard distribution denotes the Boltzmann, Fermi-Dirac, and Bose-Einstein distributions. The blast-wave model can obtain simultaneously the kinetic freeze-out temperature and transverse flow velocity from a single formula. The alternative method needs a few steps to obtain the two quantities, where the intercept in the linear relation of effective temperature against rest mass is regarded as the kinetic freeze-out temperature [7][8][9][10][11], and the slope in the linear relation of mean transverse momentum against mean moving mass (mean energy) is regarded as the transverse flow velocity [9][10][11]. The effective temperature is in fact the temperature parameter in the standard distribution in which the contribution of flow effect is not excluded. Comparatively, the blast-wave model is more convenient than the alternative method.
As we know, the energy dependent kinetic freeze-out temperature and transverse flow velocity in high energy collisions is an open question at present. We think that it is hard to conclude whether or not there is an increase, decrease, or invariant trend with the increase of energy. In fact, the situations on the energy dependence of the kinetic freeze-out temperature and transverse flow velocity in the literature are contradictory. Some of them show an increase trend [12][13][14][15][16][17][18][19][20], some of them show a decrease trend [15][16][17][18][19][20], and others show an invariant trend [12,13] in the energy dependent kinetic freezeout temperature and/or transverse flow velocity, with the increase of energy from the RHIC to LHC. Our recent work [21] shows an increase or invariant trend. As for the energy range from the SIS to SPS, the trend is always incremental in different articles in the literature. Indeed, It is necessary to do more studies on this topic and to conclude some affirmative conclusions.
In this paper, the blast-wave model with Boltzmann-Gibbs statistics [2,3] and with Tsallis statistics [4] is used to extract the kinetic freeze-out temperature and transverse flow velocity by describing the transverse momentum spectra of negative and positive pions produced at mid-(pseudo)rapidity in inelastic (INEL) or nonsingle-diffractive (NSD) proton-proton (pp or p-p) collisions and in central gold-gold (Au-Au) [copper-copper (Cu-Cu), lead-lead (Pb-Pb)] collisions over an energy range from the SIS to LHC. The model results are compared with the experimental data measured at the SIS, SPS, RHIC, and LHC by the FOPI [22], NA61/SHINE [23], PHENIX [24,25], STAR [26][27][28], ALICE [29,30], and CMS [31,32] Collaborations. In some cases, the two-component blast-wave model is used due to the fact that the single model cannot fit the data very well. This paper is structured as follows. The formalism and method are described in Section 2. Results and discussion are given in Section 3. In Section 4, we summarize our main observations and conclusions.

Formalism and method
Generally, the spectra in the high transverse momentum region are contributed by the hard scattering process which is described by the quantum chromodynam-ics (QCD) calculus [33][34][35] or the Hagedorn function [36,37] which is in fact an inverse power law which has at least, in our opinion, three revisions. In refs. [38], [39][40][41][42][43], and [44], the three revisions of the Hagedorn function (inverse power law) are given respectively, though specific nomenclatures on the revised Hagedorn functions are not mentioned in some cases. We shall not discuss the Hagedorn function and its revisions due to the fact that the hard scattering process has no contribution to the kinetic freeze-out temperature and transverse flow velocity. Instead, the soft excitation process contributes the spectra in low transverse momentum region in which the kinetic freeze-out temperature and transverse flow velocity can be extracted.
In order to extract the kinetic freeze-out temperature and transverse flow velocity, we should analyze the spectra in low transverse momentum region. Although one can choose various functions to describe the mentioned spectra, the blast-wave model with Boltzmann-Gibbs statistics [2,3] and with Tsallis statistics [4][5][6] is a convenient consideration. For the spectra in high transverse momentum region, the mentioned model is also a possible choice with larger parameter values.
According to refs. [2,3], in the blast-wave model with Boltzmann-Gibbs statistics, the first component with the kinetic freeze-out temperature, T 1 , and transverse flow velocity, β T 1 , results in the probability density function of transverse momenta, p T , to be where C 1 is the normalized constant, m T = p 2 T + m 2 0 is the transverse mass, m 0 is the rest mass, r/R is the relative radial position, I 0 and K 1 are the modified Bessel functions of the first and second kinds respectively, ρ 1 = tanh −1 [β 1 (r)] is the boost angle, β 1 (r) = β S1 (r/R) n0 is a self-similar flow profile, β S1 is the flow velocity on the surface, and n 0 = 2 is used in the original form [2]. As a mean of β 1 (r), R 0 rβ 1 (r)dr = 2β S1 /(n 0 + 2). According to ref. [4][5][6], in the blast-wave model with Tsallis statistics, the first component with T 1 and β T 1 results in the probability density function of p T to be where q is an entropy index that characterizes the degree of non-equilibrium, φ denotes the azimuthal angle, and n 0 = 1 is used in the original form [4]. It should be noted that f (p T , T 1 , β T 1 ), T 1 , β T 1 , and C 1 in Eq. (2) are different from those in Eq. (1), though the same symbols are used. The second component has the same form as the first one, but with the kinetic freeze-out temperature, T 2 , and transverse flow velocity, β T 2 . The two-component blast-wave model can be structured as where k denotes the contribution fraction (rate) of the first component. According to Hagedorn's model [36], we may also use the usual step function to structure the two-component blast-wave model. That is where A 1 and A 2 are constants which result in the two components to be equal to each other at p T = p 1 , (3) and (4) can be used to extract the kinetic freeze-out temperature, T 0 , and transverse flow velocity, β T , in the two-component blast-wave model. The first component is used to extract T 0 and β T due to its contribution in low p T region by the soft process, and the second component is not used to extract T 0 and β T due to its contribution in high p T region by the hard process. We have and In the case of using Eq. (3) to get the parameter values of two components, k is directly given by Eq. (3). In the case of using Eq. (4) to get the parameter values of two components, k is expressed by due to the fact that Eq. (4) is the probability density function which results naturally in the normalization. There are little differences between values of T 0 (β T ) extracted from Eqs. (3) and (4). Generally, T 1 (β T 1 ) from Eq. (3) is less than that from Eq. (4) by ≤ 5%, and T 2 (β T 2 ) from Eq. (3) is larger than that from Eq. (4) by ≤ 5%. We are inclined to use Eq. (3) due to its convenient fit in obtaining a smooth curve, though there is a decussate region between the contributions of the two components in Eq. (3). Oppositely, with Eq. (4) it is hard to get a smooth curve at p T = p 1 , though there is no entanglement between the two components in Eq. (4).
To describe the system and statistics clearly, we divide the subject investigated in the present work into four cases. The case I is the case which describes pp collisions by the blast-wave model with Boltzmann-Gibbs statistics. The case II is the case which describes pp collisions by the blast-wave model with Tsallis statistics. The case III is the case which describes nucleus-nucleus collisions by the blast-wave model with Boltzmann-Gibbs statistics. And the case IV is the case which describes nucleus-nucleus collisions by the blast-wave model with Tsallis statistics.

Results and discussion
The transverse momentum spectra of π − and π + produced at mid-(pseudo)rapidity in INEL or NSD pp collisions at high energies are presented in Fig.1, where different mid-(pseudo)rapidity (y or η) intervals and collision energies ( √ s) are marked in the panels. The closed (open) symbols presented in panels (a)-(e) represent the data of π − (π + ) measured by the NA61/SHINE [23], PHENIX [24], STAR [26], ALICE [29], and CMS [31,32] Collaborations, respectively, where in panel (a) only the spectra of π − are available, and panel (c) is for NSD events and the other panels are for INEL events. The solid (dashed) curves are our results calculated by Eqs.
(1) and (3) for π − (π + ) (case I), and the dotted (dotdashed) curves are our results calculated by Eqs. (2) and (3) for π − (π + ) (case II). Different forms of the spectra are used due to different Collaborations, where N , E, p, σ, and N EV denote the particle number, energy, momentum, cross-section, and event number, respectively. In some cases, different amounts marked in the panels are used to scale the data for clarity. Panels (a ′ )-(e ′ ) show the ratios of data to fit obtained from Eqs. (1) and (3) for the case I, corresponding to the collisions shown in the p T spectra; and panels (a ′′ )-(e ′′ ) show the ratios of data to fit obtained from Eqs. (2) and (3) for the case II. The values of free parameters (T 1 , T 2 , β T 1 , β T 2 , and q if available), normalization constant (N 0 ), χ 2 , and degrees of freedom (dof) corresponding to the curves in Fig. 1 are listed in Tables 1 and 2 for the cases I and II respectively. One can see that Eq. (3) describes similarly well the p T spectra at mid-(pseudo)rapidity in INEL or NSD pp collisions over an energy range from a few GeV to above 10 TeV. The main parameters (T 0 and β T ) show some laws in the considered energy range, which will be discussed later.   Transverse momentum spectra of π − and π + produced at mid-(pseudo)rapidity in INEL or NSD pp collisions at high energies, where the mid-(pseudo)rapidity intervals and energies are marked in the panels. Panels (a)-(e) represent the data measured by the NA61/SHINE [23], PHENIX [24], STAR [26], ALICE [29], and CMS [31,32] Collaborations, respectively, by various symbols, where in panel (a) only the spectra of π − are available, and panel (c) is for NSD events and other panels are for INEL events. In some cases, different amounts marked in the panels are used to scale the data for clarity. The solid (dashed) curves are our results calculated by Eqs. (1) and (3) for π − (π + ) (case I), and the dotted (dot-dashed) curves are our results calculated by Eqs. (2) and (3) for π − (π + ) (case II). Panels (a ′ )-(e ′ ) and (a ′′ )-(e ′′ ) are for ratios of the data to fit obtained from Eqs. (1) and (3) for the case I and from Eqs. (2) and (3) for the case II respectively, corresponding to the collisions shown in the pT spectra.  Table 1. Values of free parameters (T1, T2, βT 1, and βT 2), normalization constant (N0), and χ 2 /dof corresponding to the solid and dashed curves in Fig. 1 in which different data are measured in different mid-(pseudo)rapidity intervals at different energies by different Collaborations. In some cases, there is no the second component.  Figure 2 is the same as Fig. 1, but it shows the results at mid-rapidity in central Au-Au (Cu-Cu, Pb-Pb) collisions at high energies. The symbols in panels (a)-(f) represent the data measured by the FOPI [22], STAR [27], STAR [28], PHENIX [25], STAR [26], and ALICE [30] Collaborations, respectively, with different mid-rapidity intervals, centrality intervals, and collision energies. The solid (dashed) curves are our results calculated by Eqs. (1) and (3) for π − (π + ) (case III), and the dotted (dot-dashed) curves are our results calculated by Eqs. (2) and (3) for π − (π + ) (case IV). Panels (a ′ )-(f ′ ) and (a ′′ )-(f ′′ ) show the ratios of data to fit. The values of various parameters, χ 2 , and dof corresponding to the curves in Fig. 2 are listed in Tables 3 and 4 for the cases III and IV respectively. One can see that Eqs. (2) and (3)  in center-of-mass system. The energy dependent T 0 and β T are presented in the upper and lower panels in Fig. 3 respectively, where the left and right panels correspond to the results obtained from Eqs. (1) and (3) and from Eqs. (2) and (3) respectively. The closed and open symbols represent the parameter values corresponding to π − and π + respectively, which are listed in Tables 1 and 2 for the cases I and II respectively, for INEL or NSD pp collisions, and in Tables 3 and 4 for the cases III and IV respectively, for central Au-Au (Cu-Cu, Pb-Pb) collisions. One can see that the difference between the results of π − and π + can be neglected. In the energy dependent T 0 and β T in INEL or NSD pp collisions obtained by the blast-wave model with Boltzmann-Gibbs statistics (case I), there are a hill at √ s ≈ 10 GeV, a drop at dozens of GeV, and then an increase from dozes of GeV to above 10 TeV. The energy dependent T 0 and β T in other three cases show a simple structure. That is, there is a quick increase from the SIS to SPS. Then, a slight increase or similar invariability appears from the RHIC to LHC. Special structure is not observed around the energy bridge from the SPS to RHIC in other three cases, though lower T 0 and β T are obtained in the cases II and IV in which the Tsallis statistics is used.
It should be noted that the fit of the first two points in Fig. 1(a) is to high in some cases. It means T 0 (β T ) for pp collisions at the SPS could be higher than what is shown in Fig. 3 in the case of improving the fit. Contrarily, the fit of the first point in Fig. 2(d) is too low. It means T 0 (β T ) for central Cu-Cu collisions at 22.5 GeV could be lower than that what is shown in Fig. 3. We give up to improve further the fit of the first one or two points due to the trend of whole fit being worse. Alternatively, we may use a multi-component model to improve further the fit. The changed amounts in T 0 (β T ) caused by the further improvement are small.
Since T 0 (β T ) is strongly related to mean transverse momentum, p T , or root-mean-square transverse mo-mentum, p 2 T , we plot p T against √ s ( √ s N N ) and p 2 T against √ s ( √ s N N ) in the upper and lower panels in Fig. 4 respectively, to check if the discontinuous hill in the case I still exists, where the left and right panels correspond to the results obtained from Eqs. (1) and (3) and from Eqs. (2) and (3) respectively. One can see that the hill in Fig. 3 is still existent in Fig. 4 for the case I. Similar trend as what is shown in Fig. 3 is observed in Fig. 4 for the other three cases.
It should be noted that the above calculations in Fig.  4 are done according to the functions of curves in the p T range from 0 to 5 GeV/c, though the experimental p T ranges are different and much shorter. If we reduce the p T range, the trend does not change obviously due to small fraction of particles with high p T . We give up to use the experimental data to calculate p T and p 2 T due to the fact that some data are incomplete in the considered p T range, in particular in the very-low-p T region.
Comparing with the case I, the different situation of energy dependence of T 0 (β T ) in the other three cases is caused by larger collision systems for central nucleusnucleus collisions and/or one more parameter (q) in the Tsallis statistics. In central nucleus-nucleus collisions, some statistical and/or dynamical fluctuations existed in small collision systems are smoothed because of the successional or cascade nucleon-nucleon collisions and the multiple scattering of secondary particles. In the Tsallis statistics, some fluctuations existed in the Boltzmann-Gibbs statistics are smoothed due to the introduction of q which causes more entanglements in the selections of T 0 and β T . These factors cause the energy dependent T 0 (β T ) in the other three cases to be different from that in the case I.
In any case, the special structure around 10 GeV in the energy dependent T 0 (β T ) in the case I or the different trend which turns a corner in other three cases indicates that this energy is special as pointed out by Cleymans [45]. The present work shows that, in most cases, the energy dependent T 0 (β T ) increases quickly before 10 GeV and then saturates or increases slowly after 10 GeV with the increase of collision energy, which also indicates this special energy (11 GeV more specifically [45]). Some special properties result in this special energy. The final state has the highest net baryon density at this energy. Meanwhile, there is a transition from a baryon-dominated to a meson-dominated final state, and the ratios of strange particles to mesons show clear and pronounced maximums [45].
At 11 GeV, the chemical freeze-out temperature in central nucleus-nucleus collisions is T ch ≈ 151 MeV [45]. The present work shows that, at this energy, the kinetic freeze-out temperature in INEL or NSD pp collisions is T 0 ≈ 115 MeV and that in central nucleus-nucleus collisions is T 0 ≈ 120 MeV, obtained from Eqs. (1) and (3). The differences between T ch and T 0 and be- Table 2. Same as Table 1, but one more parameter (q) is added and the values correspond to the dotted and dot-dashed curves in Fig. 1.   [22], STAR [27], STAR [28], PHENIX [25], STAR [26], and ALICE [30] Collaborations, respectively, by various symbols. The solid (dashed) curves are our results calculated by Eqs. (1) and (3) for π − (π + ) (case III), and the dotted (dot-dashed) curves are our results calculated by Eqs.  Table 3. Same as Table 1, but for the solid and dashed curves in Fig. 2 in which the results at mid-rapidity in central nucleus-nucleus collisions are presented.
Collab. and Type   Table 4. Same as Table 1, but one more parameter (q) is added and the values correspond to the dotted and dot-dashed curves in Fig. 2 in which the results at mid-rapidity in central nucleus-nucleus collisions are presented.
Collab. and Type  The closed and open symbols represent the parameter values corresponding to π − and π + respectively, which are listed in Tables 1-4, where √ sNN for nucleus-nucleus collisions is omitted for concision.
tween the two T 0 render that the chemical freeze-out happens obviously earlier than the kinetic freeze-out due to higher temperature at the chemical freeze-out. The energy deposition in central nucleus-nucleus collisions is larger than that in INEL or NSD pp collisions due to higher temperature in central nucleus-nucleus collisions. The values of T 0 (β T ) obtained from Eqs. (2) and (3) are systematically lower than those obtained from Eqs.
According to an ideal fluid consideration, the time evolution of temperature follows T f = T i (τ i /τ f ) 1/3 , where T i (= 300 MeV) and τ i (= 1 fm) are the initial temperature and proper time respectively [46,47], and Generally, in the low-p T region, T 0 (β T ) obtained only from the pion spectra is less than that obtained averagely by weighting the yields of various types of light particles. The treatment in the present work on the extraction of T 0 (β T ) is approximate. However, the fraction ratio of the pion yield to total yield in high energy collisions is larger (∼85%). This approximate treatment is acceptable. The energy dependent T 0 (β T ) obtained from the pion spectra is similar to that obtained from the weighted average of the spectra of various light particles.
From the RHIC to LHC, there are three trends for the change of T 0 . That is the incremental [14], decreas-  ing [15][16][17][18][19][20], or invariant trend [12,13] in the energy dependent T 0 . However, there is only one trend for the change of β T . That is the incremental trend [12][13][14][15][16][17][18][19][20] in the energy dependent β T . The decreasing or invariant trend is not observed in the energy dependent β T . The present work shows that, in the four cases (two collisions and two statistics), both T 0 and β T do not decrease from the RHIC to LHC. The energy deposition at the LHC is more than that at the RHIC.
The inconsistent results on the trend of T 0 (β T ) are in fact caused by different models or methods. Similar to different types of thermometers in thermodynamics, the standard method extracting nuclear temperature is needed in high energy collisions. In our opinion, the standard method should be related to the Boltzmann-Gibbs statistics and the standard distribution (the Boltzmann, Fermi-Dirac, and Bose-Einstein distributions) which are the foundation of the ideal gas model in thermodynamics. The blast-wave model with Boltzmann-Gibbs statistics [2,3] and the alternative method [5][6][7][8][9][10][11] with standard distribution are suitable to be candidates, though other models and methods also fit the data very well.
In the above discussions, we have used the twocomponent blast-wave model in some cases. In our very recent work [48], a superposition of the blast-wave model and an inverse power-law has been particularly used in these cases in pp collisions. It is noted that T 0 (β T ) is determined by the first component due to the fact that the spectra in low p T range is mainly contributed by the first component. From the point of view of extraction of T 0 (β T ), one can consider only the contribution of the first component. Even the second component can be given up in the fit process and the corresponding p T region can be left there.
In the last part of this section, we would like to discuss further some issues. The description of blastwave model came from Ref. [2] assumes the local thermal equilibrium based on a hydrodynamical framework. From Cooper-Frye formula, one can derive Eq. (1) for Boltzmann-Gibbs statistics. Due to the very short interacting time, the collision system is possibly not in the thermal equilibrium with an unified temperature, but in a few local thermal equilibriums with different temperatures. Naturally and in the simplest way, we can use the two-component blast-wave model, which basically assumes that the both low-p T and high-p T particles severally obey the above hydrodynamical assumption.
Although the system of pp collisions at low energy is probably not in thermal equilibrium or local thermal equilibriums due to low multiplicity, the present work treats pp collisions at high energies in which the multiplicities in most cases are not too low. In particular, small collision systems such as pp and p-nucleus collisions show abundant collective behaviors [49] which are similar to those in nucleus-nucleus collisions. These similarities to nucleus-nucleus collisions render that the idea of local thermal equilibrium for pp collisions at high energies may be the truth. The application of the (twocomponent) blast-wave model in pp collisions at high energies is acceptable.
In addition, we would like to emphasize that, comparing with central nucleus-nucleus collisions or with the blast-wave model with Tsallis statistics, we have observed the inconsistent trend for the energy dependence of T 0 (β T ), i.e. the hill, in pp collisions analyzed by the blast-wave model with Boltzmann-Gibbs statistics. This implies that the successional or cascade nucleon-nucleon collisions and multiple scattering of secondary particles in central nucleus-nucleus collisions play important roles in the collision process. Due to these successional collisions and multiple scattering, some statistical and/or dynamical fluctuations in pp collisions are smoothed obviously in central nucleusnucleus collisions. In the blast-wave model with Tsallis statistics, the extra entropy index q smoothes naturally these fluctuations.

Conclusions
To conclude, the transverse momentum spectra of π − and π + produced at mid-(pseudo)rapidity in INEL or NSD pp collisions and in central Au-Au (Cu-Cu, Pb-Pb) collisions over an energy range from the SIS to LHC have been analyzed by the (two-component) blastwave model with Boltzmann-Gibbs statistics and with Tsallis statistics. The model results are in similarly good agreement with the experimental data of FOPI, NA61/SHINE, STAR, PHENIX, ALICE, and CMS Col-laborations. The values of related parameters, the kinetic freeze-out temperature T 0 and the transverse flow velocity β T , are extracted from the fit process and the energy dependent parameters are obtained.
In INEL or NSD pp collisions, and in the analysis by the blast-wave model with Boltzmann-Gibbs statistics, both the energy dependent T 0 and β T show a complex structure. There is a hill at √ s ≈ 10 GeV, a drop at dozens of GeV, and an increase from dozes of GeV to above 10 TeV. In central Au-Au (Cu-Cu, Pb-Pb) collisions, or in the analysis by the blast-wave model with Tsallis statistics, both the energy dependent T 0 and β T show simple structure. Form the SIS to LHC, there is a quick increase before 10 GeV, and then a slight increase or saturation after 10 GeV. The differences among the energy dependent parameters in different cases are caused by different collision systems and statistics. Large system smoothes the fluctuations due to the successional collisions and multiple scattering. The Tsallis statistics smoothes the fluctuations due to the entropy index.
Some special properties result in the special energy of 10 GeV (11 GeV more specifically [45]). At this special energy, not only the final state has the highest net baryon density, but also the transition from a baryon-dominated to a meson-dominated final state has happened. At the same time, the ratios of strange particles to mesons show clear and pronounced maximums at this energy. Regardless of the baryon-dominated and meson-dominated final state, they possibly undergo a parton-dominated intermediate state at higher energy. To search for the critical energy at which the parton-dominated intermediate state appears initially is a long-term target, though the critical energy is possibly around the energy bridge from the SPS to RHIC.

Data Availability
All data are quoted from the mentioned references. As a phenomenological work, this paper does not report new data.

Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this paper.