arxiv: v3 [hep-ph] 19 Apr 2016

Transcription

1 Kinetic freeze-out temperature and flow velocity extracted from transverse momentum spectra of final-state light flavor particles produced in collisions at RHIC and LHC Hua-Rong Wei a, Fu-Hu Liu a,, and Roy A. Lacey b, a Institute of Theoretical Physics, Shanxi University, Taiyuan, Shanxi , China b Departments of Chemistry & Physics, Stony Brook University, Stony Brook, NY 794, USA arxiv: v3 [hep-ph] 9 Apr 06 Abstract: The transverse momentum spectra of final-state light flavor particles produced in proton-proton (p-p), copper-copper (Cu-Cu), gold-gold (Au-Au), lead-lead (Pb- Pb), and proton-lead (p-pb) collisions for different centralities at relativistic heavy ion collider (RHIC) and large hadron collider (LHC) energies are studied in the framework of a multisource thermal model. The experimental data measured by the STAR, CMS, and ALICE Collaborations are consistent with the results calculated by the multi-component Erlang distribution and Tsallis Statistics. The effective temperature and real temperature (kinetic freeze-out temperature) of interacting system at the stage of kinetic freeze-out, the mean transverse flow velocity and mean flow velocity of particles, and the relationships between them are extracted. The dependences of effective temperature and mean (transverse) momentum on rest mass, moving mass, centrality, and center-of-mass energy, and the dependences of kinetic freeze-out temperature and mean (transverse) flow velocity on centrality, center-of-mass energy, and system size are obtained. Keywords: Kinetic freeze-out temperature, Flow and transverse flow velocities, Transverse momentum spectrum, Erlang distribution, Tsallis statistics PACS:.38.Mh, 5.75.Dw, 4.0.Pa Introduction Since the relativistic heavy ion collider (RHIC) and large hadron collider (LHC) successfully run, the evolution process of collision system and properties of quark-gluon plasma (QGP) [] formed in a high-temperature and high-density extreme condition, attract more interests and studies. Many theoretical and experimental methods are used to study the high energy nucleus-nucleus collisions which result in multi-particle productions. The model analysis is a simple but effective method in the study of high energy nucleus-nucleus collisions. It can extract the information of interacting system and QGP by analyzing various spectra of final-state products with distribution laws of different models. For example, by using the blast-wave model [], thermal and statistical model fuhuliu@63.com; fuhuliu@sxu.edu.cn Roy.Lacey@Stonybrook.edu

2 [3], Landau hydrodynamic model [4 7], multisource thermal model [8 0], and so forth, one can describe transverse momentum (mass) spectrum, azimuthal distribution, and rapidity distribution of final-state products to extract temperature, non-equilibrium degree, longitudinal extension, and speed of sound of interacting system, and flow velocity and chemical potential of particles [ 4], especially at the state of kinetic freeze-out, which renders that model analyzes have made great contributions to study the properties of reaction system and QGP, as well as new physics. In high energy collisions, the interacting system at the kinetic freeze-out (the last stage of collisions) stays at a thermodynamic equilibrium state or local equilibrium state, when the particle emission process is influenced by not only the thermal motion but also the flow effect. In other words, it is interesting to study the temperature of interacting system and flow velocity of particles at the stage of kinetic freeze-out. The effective temperature extracted from the transverse momentum spectrum [5 5] includes thermal motion and flow effect of particles, where the thermal motion is actually the reflection of the real temperature of emission source. From dissecting the effective temperature, it is possible to obtain the real temperature of interacting system (kinetic freeze-out temperature) and mean (transverse) flow velocity of particles. The relationships between effective temperature, real temperature, flow velocity are not totally clear. Although the theories of studying kinetic freeze-out temperature are many, their results are different from each other in some cases. This indicates that studying more their relations is important and needful. In the present work, we use the Erlang distribution with one-, two-, or three-component and Tsallis statistics in the multisource thermal model [8 0] to describe the transverse momentum spectra of final-state particles produced in proton-proton(p-p), copper-copper (Cu-Cu), gold-gold (Au-Au), lead-lead (Pb-Pb), and proton-lead (p-pb) collisions with different centrality intervals over a s NN (center-of-mass energy per nucleon pair) or s (only for p-p collision in some cases) range fro. to 7 TeV. The Monte Carlo method is used to calculate the results and to see the statistical fluctuations in the process of calculation. The calculated results are compared with the experimental data of the STAR [5 9], CMS [0], and ALICE [ 5] Collaborations. From comparison and analysis, the kinetic freeze-out temperature of interacting system, mean transverse flow velocity and mean flow velocity of particles, and their relations are then extracted. The structure of the present work is as followings. The model and formulism are shortly described in section. Results and discussion are given in section 3. In section 4, we summarize our main observations and conclusions. The model and formulism The present work is based on the multisource thermal model [8 0], which assumes that many emission sources are formed in high energy collisions. Due to the existent of different interacting mechanisms in collisions and measurement of different event samples in experiments, these sources are classified into a few groups. Generally, we assume that sources in the same group stay at a local equilibrium state, which means they have the same excitation degree and a common temperature. The emission process of all the

3 sources in different groups results in the final-state spectrum, which can be described by a multi-component distribution law, because of a local equilibrium state corresponding to a singular distribution. Using different distributions to describe transverse momentum ( ) spectra of final-state particles can obtain different information. For example, a multicomponent Erlang distribution can directly give the mean transverse momentum of each group, while the Tsallis statistics can show the effective temperature of the whole interacting system which may have group-by-group fluctuations in different local thermal equilibriums. According to the model [0], particles generated from one emission source is assumed to obey an exponential function of transverse momentum distribution f ij (p tij ) = [ p tij exp p ] tij, () p tij where p tij is the transverse momentum contributed by the ith source in the jth group, and p tij is the mean value of p tij. All the sources in the jth group results in the folding of exponential functions f j ( ) = p m j [ T (m j )! p tij exp p ] T, () m j p tij wherem j isthesource number inthejthgroup, and denotesthetransverse momentum contributed by the m j sources. This is the Erlang distribution. The contribution of all l groups of sources can be expressed as l f E ( ) = k j f j ( ), (3) j= where k j denotes the relative weight contributed by the jth group and meets the normalization l j= k j =. This is the multi-component Erlang distribution. Then, using the inverse slope parameter p tij, we can obtain the mean transverse momentum of final-state particles to be l = k j m j p tij. (4) j= The Tsallis statistics is in fact the sum of contributions of two or three standard distributions. Although the Tsallis statistics does not give each local temperature like multi-component distribution which reveals fluctuations from a local equilibrium state to another one, it use an average temperature of the whole interacting system to describe the effect of local temperature fluctuations. So, the Tsallis statistics is widely used in high energy collisions [6 35]. According to the Tsallis statistics [6 3], we use the formalism of unit-density function of and rapidity (y) d N [ = C p T dydp +coshy T + q T ] q p q T +coshy, (5) 3

4 where N is the number of particles, C = gv/(π) is the normalization constant, g and V are degeneracy factor and volume respectively, is the rest mass of considered particle, T is the mean effective temperature over fluctuations in different groups, and q (q > ) is the factor (entropy index) to characterize the degree of non-equilibrium among different groups. With an integral for y in Eq. (5), the normalized Tsallis distribution is obtained and can be written as f T ( ) = N dn ymax [ = C T p T dp + coshy T y min + q T ] q p q T +coshy dy, (6) where C T denotes the normalization constant which results in 0 f T ( )d =, y min is the minimum rapidity, and y max is the maximum rapidity. Based on the above two distributions, we can use the Monte Carlo method to obtain a series of. Under the assumption of isotropic emission in the source rest frame, the space angle θ and azimuthal angel φ of particles satisfy the distributions of f θ (θ) = (/)sinθ and f φ (φ) = /(π) respectively. By the Monte Carlo method, a series of θ and φ are obtained. Correspondingly, the x-component, y-component, and (longitudinal) z-component of momentum are p x = cosφ, p y = sinφ, and p z = /tanθ, respectively. Then, the momentum p = /sinθ or p = p T +p z, the energy E = p +m 0, the Lorentz factor γ = / (p/e), and the moving mass m = γ, as well as their averages, p, γ, and m are acquired. Particularly, the values of mean transverse momentums according to analytical function and Monte Carlo method are almost the same. 3 Results and discussion Figure presents the transverse momentum spectra of various identified hadrons produced in p-p collision at center-of-mass energy (a) s = 0., (b) 0.9, (d).76, and (e) (f) 7 TeV, where N EV on the axis denotes the number of inelastic collisions events. The symbols represent the experimental data of (a) π +, K +, and p measured by the STAR Collaboration at midrapidity y < 0. [5], (b) (d) (e) π +, K +, and p measured by the CMS Collaboration in the range y < [0], Λ, φ, and Ξ + Ξ + measured by the ALICE Collaboration in the range y < 0.8 [], as well as (f) Ξ and Ω measured by the ALICE Collaboration in the range y < []. The errors include both the statistical and systematic errors. The solid and dashed curves are our results calculated by using the one- or two-component Erlang distribution and Tsallis statistics, respectively. The values of free parameters (m, p ti, k, m, and p ti ), normalization constant (N E0 ), and χ per degree of freedom (χ /dof) corresponding to the one- or two-component Erlang distribution are listed in Table, and the values of free parameters (T and q), normalization constant (N T0 ), and χ /dof corresponding to the Tsallis statistics are given in Table, where the normalization constants (N E0 and N T0 ) are used for comparisons between the normalized curves and data points. In particular, the value of χ /dof for Ξ + Ξ + in Figure is in fact the value of χ due to the number of data points being less than that of parameters. One can see that the one- or two-component Erlang distribution and Tsallis statistics describe the experimental data of the considered particles in p-p collision 4

5 at different energies. From Tables and, one can see that the numbers of sources in different groups are, 3, or 4. The effective temperature T increases and the non-equilibrium degree parameter q decreases with increase of rest mass, which reflects non-simultaneous productions of different types of particles, while T and q increase with increase of centerof-mass energy. The normalization constants N E0 and N T0 decrease with increase of rest mass, and increase with increase of center-of-mass energy. It is interesting to note that the product of T and q increases with increase of rest mass and center-of-mass energy. Figure shows the transverse momentum spectra of (a) KS 0, (b) Λ, Ξ, and (d) Ω producedincu-cucollisionsat s NN = 0.TeV.Thesymbolsrepresent theexperimental data of the STAR Collaboration in y < and different centrality (C) intervals of 0 0%, 0 30%, and 40 60% [6]. The error bars are combined statistical and systematic errors. The solid and dashed curves are our results calculated by using the one- or twocomponent Erlang distribution and Tsallis statistics, respectively. For clarity, the results for different intervals are scaled by different amounts shown in the panels. The values of free parameters, normalization constants, and χ /dof are displayed in Tables and. Obviously, the one- or two-component Erlang distribution and Tsallis statistics describe well the experimental data of the considered particles in 0. TeV Cu-Cu collisions with different centrality intervals. The numbers of sources in different groups are,, 3, or 4. The parameter T increases and the parameter q decreases with increases of rest mass and centrality. The product of T and q increases with increases of rest mass and centrality. The parameters N E0 and N T0 decrease with increases of rest mass and decrease of centrality. The spectra of (a) π + for y <, (b) KS 0 for y <, p for y <, (d) φ for y <, (e) Λ for y <, and (f) Ξ for y < 0.75 produced in Au-Au collisions at s NN = 0. TeV as a function of centrality are given in Figure 3. The experimental data were recorded by the STAR Collaboration [7 9], and scale factors for different centralities are applied to the spectra in the panels for clarity. The uncertainties on the data points for π +, KS, 0 p, and Λ are statistical and systematic combined. While the uncertainties for φ and Ξ are only statistical uncertainties and systematic uncertainties respectively. The results calculated by using the one-, two-, or three-component Erlang distribution and Tsallis statistics are shown in the solid and dashed curves, respectively. The values of free parameters (m, p ti, k, m, and p ti ), normalization constant, and χ /dof corresponding to the one- or two-component Erlang distribution in Figures 3(b) 3(f) are listed in Table. The values of free parameters (m, p ti, k, m, p ti, k, m 3, and p ti3 ), normalization constant, and χ /dof corresponding to the three-component Erlang distribution in Figure 3(a) are listed in Table 3. The values of free parameters, normalization constant, and χ /dof corresponding to the Tsallis statistics are given in Table. Once more, the two types of distributions are in good agreement with the experimental data of the considered particles in 0. TeV Au-Au collisions with different centrality intervals. The numbers of sources in different groups are,, 3, or 4. The parameter T increases and the parameter q decreases with increases of rest mass and centrality. The product of T and q increases with increases of rest mass and centrality. ThenormalizationconstantsN E0 andn T0 decreasewithincreaseofrestmassanddecrease of centrality. The spectra of (a) π +, (b) K +, p, and (d) φ produced in central (0 5%), semi- 5

6 central (50 60%), and peripheral (80 90%) Pb-Pb collisions at s NN =.76 TeV are displayed in Figure 4. The symbols represent the experimental data measured by the ALICE Collaboration at midrapidity y < [3, 4]. The uncertainties on the data points are combined statistical and systematic ones. The fitted results with the oneor two-component Erlang distribution and Tsallis statistics are plotted by the solid and dashed curves, respectively. The values of free parameters, normalization constants, and χ /dof are exhibited in Tables and. As can be seen, the two types of distribution laws are consistent with the experimental data of the considered particles in.76 TeV Pb-Pb collisions with different centrality classes. The numbers of sources in different groups are,, or 3. The effective temperature T increases with increases of rest mass and centrality. The parameter q decreases with increases of rest mass and centrality. The product of T and q increases with increases of rest mass and centrality. The parameters N E0 and N T0 decrease with increase of rest mass and decrease of centrality. Figure 5 exhibits the spectra of (a) π + + π, (b) K + + K, p + p, and (d) Λ + Λ produced in central (0 5%), semi-central (40 60%), and peripheral (80 00%) p-pb collisions at s NN = 5.0 TeV. The ALICE experimental data in 0 < y < are represented by different symbols [5]. The error bars are combined statistical and systematic errors. Our results analyzed by the two-component Erlang distribution and Tsallis statistics are given by the solid and dashed curves, respectively. The values of free parameters, normalization constants, and χ /dof are summarized in Tables and. Once again, the experimental data can be well described by the two types of fit functions for spectra in all centrality bins. The numbers of sources in different groups are, 3, or 4. The effective temperature T increases with increases of rest mass and centrality. With increases of rest mass and centrality, the parameters q, N E0, and N T0 decrease, and the product of T and q increases. In the above descriptions, it is easy to notice that the numbers of sources m j in the jth group are in the range m j 4. The m j is so small that we think these sources corresponding to a few partons which include sea and valent quarks and gluons in highenergy collisions. Generally, the transverse momentum spectrum is contributed by the sum of soft and hard parts. The soft excitation process is a strong interacting process where a few sea quarks andgluons taken part in, and the hard scattering process is a more violent collision among a few valent quarks in incident nucleons. In spectrum, the soft excitationandhardscatteringprocesses correspondtoanarrowlow- andawidehigh- regions respectively [36 38]. And in the describing by two-component Erlang distribution, the first and second components correspond to the soft and hard processes respectively. Although the low- region contributed by soft process is narrow, the contribution of soft excitation is main, which can be seen from the relative weight k > 50%. Particularly, when the region of spectrum is narrow enough, the two-component Erlang distribution actually is the (one-component) Erlang distribution for the contribution of the second component being neglected. On the contrary, when the high- region is very wide, the hard process can not be described by one component distribution, which means the two-component distribution would expand to the three-component Erlang distribution. The temperature parameter T extracted from the Tsallis distribution is actually the effective temperature of emission sources, which can not reflect the transverse excitation of interacting system at the stage of kinetic freeze-out. In order to obtain the real temper- 6

7 ature(kinetic freeze-out temperature) of emission sources and(transverse) flow velocity of final-state particles, we study the linear dependences of mean transverse momentum, mean momentum p, and effective temperatures T on particle rest mass and mean moving mass m. Figures 6 8 show the center-of-mass energy and centrality dependences of, p, and T on respectively, and Figures 9 show that on m respectively. One can see that the three quantities, p, and T increase with increase of and m. From p-p collision, the three quantities increase with increase of center-of-mass energy, and from Cu-Cu, Au-Au, Pb-Pb, and p-pb collisions, the three quantities increase with increase of centrality. To see clearly the dependences of, p, and T on and m, we fit the linear relationships which can be written as and = C 0 +k T, (7) p = C p +k p, (8) T = T 0 +k, (9) = C 0 + u T m, (0) p = C p + u m, () T = C T +k m, () respectively, where the units of temperature, momentum, velocity, and mass are GeV, GeV/c, c, and GeV/c, respectively, where c = is in natural units; The intercepts (C 0, C p, T 0, C 0, C p and C T ) have the same units as corresponding dependent variables; The slopes k T, k p, u T, and u are in the units of c, while k and k are in the units of c. The values of intercepts, slopes, and χ are given in Table 4. From Figures 6 and Table 4, we can see the intercepts and slopes in different linear correlations. In all cases, the intercepts with different center-of-mass energies (p-p collision) or centrality bins (Cu-Cu, Au-Au, Pb-Pb, or p-pb collisions) have the tendency of converging to one point, which means the intercepts are nearly equal to each other or do not change obviously with center-of-mass energy or centrality, while the slopes have the tendencies of increasing with center-of-mass energy and centrality. The intercepts and slopes obtained from, p, and T correlations are larger than those from m, p m, and T m correlations, but the changes about intercepts are much larger than those about slopes. Besides, the increments of slopes with center-ofmass energy or centrality in, p, and T correlations are larger than those in m, p m, and T m correlations, which renders slow changes in the latter cases. In function T = T 0 + k, the quantity T extracted directly from the distribution mentioned above, is the effective temperature which includes thermal motion and flow effect of particles. As the temperature corresponding to massless ( = 0) particle, T 0 is regarded as the source real temperature at the kinetic freeze-out (or the kinetic freezeout temperature of interacting system) [39 43]. The flow effect of particles is shown in quantity of k, where the slope k has the dimension of the square of velocity. At the 7

8 same time, in function T = C T m, although the intercept C +k T has the same dimension as temperature, it s values shown in Figure and Table 4 approximately equal to zero, which have no physics meaning. Incorrelations, m, p, and p m, the slopes (k T, u T, k p, and u ) have the dimension of velocity and are considered to be related to mean transverse flow velocity and mean flow velocity. It is interesting to find the relationships between the intercepts and slopes, especially relationships between the several flow velocities and k (or k ). Figure exhibits the correlations between intercepts (a) C p C 0, (b) C p C 0, and T 0 C 0 in different collision systems with different centrality classes at different energies, and corresponding fitting are executed. One can see that C p and C 0, C p and C 0 have explicit linear relations C p = (π/)c 0 and C p = (π/)c 0 with almost zero χ, which is duetotheassumptionofisotropicemissioninthesourcerest frame. InT 0 C 0 correlation, the T 0 increases with increase of C 0, and they are basically compatible with the linear relation T 0 = (0.9 ±0.073)C 0 with χ /dof = Other relations among intercepts do not show an obvious law and are not shown in the panels due to trivialness. At the same time, the correlations between different slopes (a) k p k T, (b) k k T, k T u T, (d) u u T, and (e) k u T in different collision systems with different centrality classes at different energies, as well as corresponding fitting are shown in Figure 3. Once again, based on the assumption of isotropic emission in the source rest frame, the k p and k T from p and, u and u T from p m and m also have the same explicit linear relations k p = (π/)k T and u = (π/) u T with almost zero χ. From Figures 3(b) and 3(d), one can see that k (or k ) increases with increase of k T (or u T ), and the parameters are in good agreement with the fitted power function relations k = [(/3)±0.08]kT with χ /dof =.546 and k = [(/)±0.70] u T with χ /dof =.47, respectively. The relationship between k T and u T respectively from and m correlations is fitted by the line k T = (3.0±0.84) u T (0.74±0.08) with χ /dof = Other relations among slopes do not show an obvious law and are not shown in the panels due to trivialness. In addition, we would like to point out that in the Monte Carlo calculation, some conservation laws (such as energy conservation and momentum conservation) and physics limitation (such as flow velocity < c) are used so that we can obtain reasonable values and relations. It is noticed that most of u T values (0.339 c) extracted from m are slightly less than c and k T ( c) from, while most of u values (3 0.80c) extracted from p m is more close to c and less than k p (46.44c) from. Considering overlarge k T and k p, we regard u T and u as the mean transverse flow velocity and mean flow velocity of particles respectively, and they obey the relation of u = (π/) u T. In addition, u T and k approximatively obey k = (/) u T. The mean T 0 of interacting system obtained from T does dot change in error range or approximates independent of energy, centrality, particle type, and system size. The mean T 0 corresponding to the Tsallis distribution is (0.07 ± 0.007) GeV, which is less than that [(0.37 ± 0.007)] GeV from the Boltzmann distribution in previous work [43], and is also less than that (0.77 GeV) from an exponential shape of transverse mass spectrum [4]. Our result is close to that (0.098 GeV) from the blast-wave model [44]. Comparing with chemical freeze-out temperature, the kinetic freeze-out temperature in 8

9 present work is obviously less than that (0.70 GeV) of theoretical critical point of the QCD (quantum chromodynamics) prediction [45 47], and that (6 GeV) from particle ratios in a thermal and statistical model [3]. In the above discussions, we present in fact a method to extract the kinetic freezeout temperature and (transverse) flow velocity from transverse momentum spectra in the multisource thermal model [8 0], in which the sources are described by different distribution laws. The method to extract the kinetic freeze-out temperature is also used in other literature [39 43], and the method to extract the (transverse) flow velocity is seldom discussed elsewhere. Although the method used in the present work is different from the blast-wave model [] which discusses radial flow, it offers anyhow another way to study the kinetic freeze-out temperature and (transverse) flow velocity. In our recent work [48], the multisource thermal model is revised with the blast-wave picture, where we assume that the fragments and particles produced by thermal reason in the sources are pushed away by a blast-wave. The blast-wave causes the final-state products to be effected by isotropic and anisotropic flows. 4 Conclusions From the above discussions, we obtain following conclusions. (a) The transverse momentum distributions of final-state particles produced in p-p, Cu-Cu, Au-Au, Pb-Pb, and p-pb collisions with different centrality bins over an energy range fro. to 7 TeV, can be described by one-, two-, or three-component Erlang distribution and Tsallis statistics in the framework of multisource thermal model. The calculated results are consistent with the experimental data of π ±, K ±, K 0 S, p, p, φ, Λ, Λ, Ξ, and Ω, etc. measured by the STAR, CMS, and ALICE Collaborations. (b) In most cases, the data of spectra are fitted by the two-component Erlang distribution, where the first component corresponding to a narrow low- region is contributed by the soft excitation process in which a few sea quarks and gluons taken part in and accounts for a larger proportion, and the second component for a wide high- region indicates hard scattering process which is a more violent collision among a few valent quarks in incident nucleons. The mean transverse momentum extracted from multi-component Erlang distribution increases with increases of center-of-mass energy, particle mass, event centrality, and system size. TheTsallis distributionusestwofreeparameterst andq todescribetheeffective temperature and the non-equilibrium degree of the interacting system respectively. The present work shows that q decreases with increases of particle mass and event centrality, and increases with increase of center-of-mass energy. In physics, q = corresponds to an equilibrium state and a large q corresponds to a state departing far from equilibrium. Our study indicates that a high center-of-mass energy results in interacting system deviating from the equilibrium state, and the interacting system of central collisions and heavier particles are closer to the equilibrium state. The extracted T increases with increases of center-of-mass energy, particle mass, event centrality, and system size. (d) The intercept T 0 in T correlation is regarded as the mean kinetic freezeout temperature of interacting system, which do not depend on center-of-mass energy, 9

10 event centrality, and system size in error range, which renders that different interacting systems at the stage of kinetic freeze-out stay at the same phase state. The average of T 0 corresponding to the Tsallis distribution is (0.07±0.007) GeV, which is less than the effective temperature, chemical freeze-out temperature, and expected critical temperature (30 65 MeV) of the QGP formation [49]. (e) The mean transverse flow velocity u T and mean flow velocity u of particles are obtained from m and p m correlations respectively. The present work shows that u T (0.339-c) and u (3-0.80c) have narrower range and are more close to c than those from and p, respectively. On the assumption of isotropic emission in the source rest frame, there is u = (π/) u T. Besides, u T and k from corresponding T m correlation approximatively meet k = (/) u T, which means u T and u can be obtained from correlation T m based on k = (/) u T and u = (π/) u T relationships. In addition, the mean transverse flow velocity and mean flow velocity have a slightly increase tendency with center-of-mass energy, event centrality, and system size. Conflict of Interests The authors declare that there is no conflict of interests regarding the publication of this paper. Acknowledgments This work was supported by the National Natural Science Foundation of China under Grant No and the US DOE under contract DE-FG0-87ER4033.A008. References [] D.H. Rischke, Prog. Part. Nucl. Phys. 5, 97 (004). [] E. Schnedermann, J. Sollfrank, U.W. Heinz, Phys. Rev. C, 48, 46 (993). [3] J. Stachel, A. Andronic, P. Braun-Munzinger, K. Redlich, J. Phys. Conf. Ser. 509, 009 (04). [4] L.D. Landau, Izvestiya akademii nauk: series fizicheskikh 7, 5 (953), in English- Translation: Collected Papers of L. D. Landau, D. Ter-Haarp, Ed., p. 569, Pergamon, Oxford, UK, 965. [5] I.M. Khalatnikov, Journal of Experimental and Theoretical Physics, 7, (954) (in Russian). [6] P.A. Steinberg, Nucl. Phys. A, 75, 43 (005). [7] L.-N. Gao, F.-H. Liu, Adv. High Energy Phys. 05, (05). 0

11 [8] F.-H. Liu, Nucl. Phys. A 80, 59 (008). [9] F.-H. Liu, J.-S. Li, Phys. Rev. C 78, (008). [0] F.-H. Liu, Y.-Q. Gao, T. Tian, B.-C. Li, Eur. Phys. J. A 50, 94 (04). [] F.-H. Liu, H.-R. Wei, R. A. Lacey, Eur. Phys. J. A, 5, 43 (05). [] L.-N. Gao, F.-H. Liu, Adv. High Energy Phys. 05, 8473 (05). [3] H.-R. Wei, Y.-H Chen, L.-N. Gao, F.-H. Liu, Adv. High Energy Phys. 04, 7863 (04). [4] H. Zhao, F.-H. Liu, Adv. High Energy Phys. 04, (04). [5] STAR Collaboration (B.I. Abelev et al.), Phys. Rev. C 79, (009). [6] STAR Collaboration (G. Agakishiev et al.), Phys. Rev. Lett. 08, 0730 (0). [7] STAR Collaboration (B.I. Abelev et al.), Phys. Rev. Lett. 97, 530 (006). [8] STAR Collaboration (B.I. Abelev et al.), Phys. Rev. Lett. 99, 30 (007). [9] STAR Collaboration (J. Adams et al.), Phys. Rev. Lett. 98, 0630 (007). [0] CMS Collaboration (S. Chatrchyan et al.), Eur. Phys. J. C 7, 64 (0). [] ALICE Collaboration (K. Aamodt et al.), Eur. Phys. J. C 7, 594 (0). [] ALICE Collaboration (M. Floris et al.), J. Phys. G 38, 405 (00). [3] ALICE Collaboration (B. Abelev et al.), Phys. Rev. C 88, (03). [4] ALICE Collaboration (B. Abelev et al.), Phys. Rev C 9, (05). [5] ALICE Collaboration (B. Abelev et al.), Phys. Lett. B 78, 5 (04). [6] C. Tsallis, J. Stat. Phys. 5, 479 (988). [7] C. Tsallis, Braz. J. Phys. 39, 337 (009). [8] W.M. Alberico, A. Lavagno, Eur. Phys. J. A 40, 33 (009). [9] G. Wilk, Z. W lodarczyk, Phys. Rev. C 79, (009). [30] G. Wilk, Z. W lodarczyk, in Proceeding of the 6th International Workshop on Relativistic Aspects of Nuclear Physics (RANP000), São Paulo, Brazil, October 7-0, 000, edited by T. Kodama, C.E. Aguiar, H.-T. Elze, F. Grassi, Y. Hama, G. Krein (World Sciencetific, Singapore, 00) pp , arxiv:hep-ph/0089v. [3] W.M. Alberico, P. Czerski, A. Lavagno, M. Nardi, V. Somá, arxiv:hep-ph/0507v.

12 [3] B.-C. Li, Y.-Z. Wang, F.-H. Liu, Phys. Lett. B 75, 35 (03). [33] T.S. Biró, Eur. Phys. J. A 40, 55 (009). [34] T.S. Biró, Physica A 39, 33 (03). [35] T.S. Biró, G. Purcsel, K. Urmossy, Eur. Phys. J. A 40, 35 (009). [36] J. Cleymans, D. Worku, Eur. Phys. J. A 48, 60 (0). [37] K. Urmossy, G.G. Barnaföldi, Sz Harangozó, T.S. Biró, Z. Xu, arxiv: [38] K. Urmossy, T.S. Biró, G.G. Barnaföldi, Z. Xu, arxiv: [39] U.W. Heinz, Concepts of heavy-ion physics, Lecture Notes for Lectures Presented at the nd CERN Latin-American School of High-Energy Physics, San Miguel Regla, Mexico, June -4, 003, arxiv:hep-ph/ [40] PHENIX Collaboration (S.S. Adler et al.), Phys. Rev. C 69, (004). [4] S. Takeuchi, K. Murase, T. Hirano, P. Huovinen, Y. Nara, Phys. Rev. C 9, (05). [4] R. Russo, Measurement of D + meson production in p-pb collisions with the ALICE detector, Ph.D. thesis, Universitaâ degli Studi di Torino, Italy, 05, arxiv: [43] H.-R. Wei, F.-H. Liu, R. A. Lacey, arxiv: [44] ALICE Collaboration (J. Adam et al.), arxiv: [45] Z. Fodor, S.D. Katz, JHEP 0404, 050 (004). [46] S. Datta, R.V. Gavai, S. Gupta, Nucl. Phys. A , 883c (03). [47] N. Xu for the STAR Collaboration, Nucl. Phys. A 93, (04). [48] F-H Liu, H-L Lao, Indian. J. Phys. 90, accepted (06). [49] M. Floris, Nucl. Phys. A 93, 03 (04).

13 Table. Values of free parameters, normalization constant, and χ /dof corresponding to one- or two-component Erlang distribution in Figures 5 except Figure 3(a). The value of χ /dof for Ξ + Ξ + in Figure is in fact the value of χ due to less data points. Figure Type m p ti (GeV/c ) k m p ti (GeV/c ) N E0 χ /dof Figure (a) π ± ± K + p 0.80 ± ± ± ± Figure (b) π ± ± ± ± K ± ± ± ± p 4 0 ± ± ± ± Figure Λ ± ± φ ± ± Ξ + Ξ ± ± (0.8) Figure (d) π ± ± ± ± K + p ± ± ± ± ± ± ± ± Figure (e) π ± ± ± ± K ± ± ± ± p ± ± ± ± Figure (f) Ξ ± ± ± (.30 ± 00) Ω ± ± ± (.480 ± 0.600) Figure (a) 0-0% 0.30 ± ± ± ± % 0.30 ± ± ± ± ± ± ± ± Figure (b) 0-0% ± ± ± ± % ± ± ± ± ± ± ± ± Figure 0-0% ± ± % ± ± ± ± Figure (d) 0-0% ± ± % ± ± ± ± Figure 3(b) ± ± ± ± % ± ± ± ± % ± ± ± ± ± ± ± ± % 0.30 ± ± ± ± Figure 3 0-% ± ± ± ± % ± ± ± ± % ± ± ± ± ± ± ± ± % ± ± ± ± Figure 3(d) ± ± % ± ± % ± ± % ± ± % 0.45 ± ± Figure 3(e) ± ± ± ± % ± ± ± ± % ± ± ± ± ± ± ± ± % ± ± ± ± Figure 3(f) ± ± % ± ± % ± ± ± ± % ± ± Figure 4(a) 0.34 ± ± ± ± % 0.94 ± ± ± ± % 0.70 ± ± ± ± Figure 4(b) ± ± ± ± % ± ± ± ± % ± ± ± ± Figure ± ± % ± ± ± ± % ± ± ± ± Figure 4(d) ± ± % 80 ± ± % 70 ± ± Figure 5(a) 0.75 ± ± ± ± ± ± ± ± % 0.75 ± ± ± ± Figure 5(b) ± ± ± ± ± ± ± ± % ± ± ± ± Figure ± ± ± ± ± ± ± ± % ± ± ± ± Figure 5(d) ± ± ± ± ± ± ± ± % ± ± ± ±

14 Table. Values of free parameters, normalization constants, and χ /dof corresponding to Tsallis distribution in Figures 5. The value of χ /dof for Ξ + Ξ + in Figure is the value of χ due to less data points. Figure Type T (GeV) q N T0 χ /dof Figure (a) π ± ± ± K + p 0.05 ± ± ± ± ± ± Figure (b) π ± ± ± K + p 0. ± ± ± ± ± ± Figure Λ 0.75 ± ± ± φ 0.74 ± ± ± Ξ + Ξ ± ± ± 0.00 (0.348) Figure (d) π ± ± ± K ± ± ± p 9 ± ± ± Figure (e) π ± ± ± K + p 0.7 ± ± ± ± ± ± Figure (f) Ξ 0.88 ± ± 0.00 (.389 ± 0.34) Ω 0.5 ± ± (.80 ± 0.3) Figure (a) 0-0% 0.3 ± ± ± % 0.8 ± ± ± ± ± ± Figure (b) 0-0% 0.7 ± ± ± % 0.98 ± ± ± ± ± ± Figure 0-0% 0.43 ± ± ± % 0.8 ± ± ± ± ± ± Figure (d) 0-0% ± ± ± % 0.33 ± ± ± ± ± ± Figure 3(a) 0-% 0. ± ± ± % 0.0 ± ± ± % 0.09 ± ± ± ± ± ± % 0.03 ± ± ± Figure 3(b) 0.66 ± ± ± % 0.63 ± ± ± % 0.6 ± ± ± ± ± ± % 7 ± ± ± Figure 3 0-% 0.5 ± ± ± % 0.0 ± ± ± % 0.03 ± ± ± ± ± ± % 0.60 ± ± ± Figure 3(d) 0.55 ± ± ± % 0.53 ± ± ± % 0.50 ± ± ± % 0. ± ± ± % 0.03 ± ± ± Figure 3(e) 0.88 ± ± ± % 0.78 ± ± ± % 0.73 ± ± ± ± ± ± % 0.50 ± ± ± Figure 3(f) ± ± ± % 0.35 ± ± ± % ± ± ± ± ± ± % 0.38 ± ± ± Figure 4(a) 0. ± ± ± % 0.09 ± ± ± % 0.07 ± ± ± Figure 4(b) 0.58 ± ± ± % 0.60 ± ± ± % 0.5 ± ± ± Figure ± ± ± % ± ± ± % 0.00 ± ± ± Figure 4(d) ± ± ± % ± ± ± % 0.56 ± ± ± Figure 5(a) ± ± ± ± ± ± % 0.07 ± ± ± Figure 5(b) 0.77 ± ± ± ± ± ± % ± ± ± Figure ± ± ± ± ± ± % 0.09 ± ± ± Figure 5(d) 0.30 ± ± ± ± ± ± % 0.00 ± ± ± Table 3. Values of free parameters, normalization constant, and χ /dof corresponding to the three-component Erlang distribution in Figure 3(a), where m = m = m 3 = which are not listed in the column. Figure Type p ti (GeV/c ) k p ti (GeV/c ) k p ti3 (GeV/c ) N E0 χ /dof Figure 3(a) 0-% 0.0 ± ± ± ± ± ± % 0. ± ± ± ± ± ± % 0.5 ± ± ± ± ± ± ± ± ± ± ± ± % 0.00 ± ± ± ± ± ±

15 Table 4. Values of intercepts, slopes, and χ corresponding to the lines in Figures 6. Figure Correlation Type Intercept Slope χ Figure 6(a) p-p 0. TeV 0.3 ± ± p-p 0.9 TeV 0.45 ± ± p-p.76 TeV ± ± p-p 7 TeV ± ± Figure 6(b) 0-0% 0.95 ± ± % ± ± ± ± Figure 6 0-% ± ± % ± ± % ± ± ± ± % 0.96 ± ± Figure 6(d) 0.37 ± ± % 0.36 ± ± % 0.35 ± ± Figure 6(e) ± ± ± ± % ± ± Figure 7(a) p p-p 0. TeV 06 ± ± p-p 0.9 TeV ± ± p-p.76 TeV 9 ± ± p-p 7 TeV 0.66 ± ± Figure 7(b) 0-0% ± ± % ± ± ± ± Figure 7 0-% 54 ± ± % 76 ± ± % 6 ± ± ± ± % ± ± Figure 7(d) 85 ± ± % 69 ± ± % 0 ± ± Figure 7(e) ± ± ± ± % 90 ± ± Figure 8(a) T p-p 0. TeV ± ± p-p 0.9 TeV ± ± p-p.76 TeV ± ± p-p 7 TeV ± ± Figure 8(b) 0-0% ± ± % ± ± ± ± Figure 8 0-% ± ± % ± ± % ± ± ± ± % ± ± Figure 8(d) ± ± % ± ± % ± ± Figure 8(e) ± ± ± ± % ± ± Figure 9(a) m p-p 0. TeV 0.67 ± ± p-p 0.9 TeV 0.07 ± ± p-p.76 TeV 0.48 ± ± p-p 7 TeV 0.60 ± ± Figure 9(b) 0-0% ± ± % 0.99 ± ± ± ± Figure 9 0-% 0.05 ± ± % 0.6 ± ± % 0.3 ± ± ± ± % ± ± Figure 9(d) 0.08 ± ± % 0.0 ± ± % ± ± Figure 9(e) 0.0 ± ± ± ± % 0.47 ± ± Figure 0(a) p m p-p 0. TeV 0.6 ± ± p-p 0.9 TeV 0.36 ± ± p-p.76 TeV 0.3 ± ± p-p 7 TeV 0.5 ± ± Figure 0(b) 0-0% ± ± % 0.3 ± ± ± ± Figure 0 0-% 0.65 ± ± % 0.83 ± ± % 0.78 ± ± ± ± % 0.0 ± ± Figure 0(d) 0.7 ± ± % 0.60 ± ± % ± ± Figure 0(e) 8 ± ± ± ± % 0.3 ± ± Figure (a) T m p-p 0. TeV ± ± p-p 0.9 TeV ± ± p-p.76 TeV ± ± p-p 7 TeV ± ± Figure (b) 0-0% 0.09 ± ± % ± ± ± ± Figure 0-% 0.00 ± ± % 0.00 ± ± % ± ± ± ± % ± ± Figure (d) ± ± % ± ± % 0.0 ± ± Figure (e) ± ± ± ± % 0.06 ± ±

16 (/π ) (d N/dyd ) ( - ) (a) p-p, 0. TeV, y <0. π + K + p (x0.) Erlang Tsallis (/N EV ) (d N/dyd ) ( - ) (b) p-p, 0.9 TeV, y < π + K + p Erlang Tsallis (/N EV ) (d N/dyd ) ( - ) p-p, 0.9 TeV, y <0.8 Λ φ Ξ +Ξ + Erlang Tsallis (/N EV ) (d N/dyd ) ( - ) (d) p-p,.76 TeV, y < π K + p Erlang Tsallis (/N EV ) (d N/dyd ) ( - ) (e) p-p, 7 TeV, y < π K + p Erlang Tsallis (/πn EV ) (d N/dyd ) ( - ) (f) p-p, 7 TeV, y < Ξ Ω Erlang Tsallis Figure. Transverse momentum spectra of various identified hadrons produced in p-p collisions at (a) s = 0., (b) 0.9, (d).76, and (e)(f) 7 TeV. The symbols represent the experimental data of (a) π +, K +, and p measured by the STAR Collaboration at midrapidity y < 0. [5], (b)(d)(e) π +, K +, and p measured by the CMS Collaboration in the range y < [], Λ, φ, and Ξ + Ξ + measured by the ALICE Collaboration in the range y < 0.8 [], as well as (f) Ξ and Ω measured by the ALICE Collaboration in the range y < [3]. The errors include statistical and systematic errors. The solid and dashed curves are our results with oneor two-component Erlang distribution and Tsallis statistics, respectively. For clarity, the results for p in Figure (a) are scaled by 0. shown in the panels. 6

17 (/π ) (d N/dyd ) ( - ) Cu-Cu, 0. TeV, y <, K 0 S 0-0% (x0 0 ) 0-30% (x0 - ) (x0-4 ) Erlang Tsallis (/π ) (d N/dyd ) ( - ) Cu-Cu, 0. TeV, y <, Λ 0 (a) 0 (b) Cu-Cu, 0. TeV, y <, Ξ Cu-Cu, 0. TeV, y <, Ω (/π ) (d N/dyd ) ( - ) (/π ) (d N/dyd ) ( - ) (d) Figure. Transverse momentum spectra of (a) K 0 s, (b) Λ, Ξ, and (d) Ω produced in Cu- Cu collisions at s NN = 0. TeV. The symbols represent the experimental data of the STAR Collaboration in y < and different centrality intervals of 0 0%, 0 30%, and 40 60% [6]. The error bars are combined statistical and systematic errors. The solid and dashed curves are our results calculated by using the one- or two-component Erlang distribution and Tsallis statistics, respectively. For clarity, the results for different intervals are scaled by different amounts shown in the panels. 7