arxiv: v3 [nucl-ex] 4 Oct 2016

Transcription

1 Disentangling random thermal motion of articles and collective exansion of source from transverse momentum sectra in high energ collisions Hua-Rong Wei a, Fu-Hu Liu a,1, and Ro A. Lace b,2 a Institute of Theoretical Phsics, Shanxi Universit, Taiuan, Shanxi , China b Deartments of Chemistr & Phsics, Ston Brook Universit, Ston Brook, NY 11794, USA arxiv: v3 [nucl-ex] 4 Oct 2016 Abstract: In the framework of a multisource thermal model, we describe exerimental results of the transverse momentum sectra of final-state light flavour articles roduced in gold-gold (Au-Au), coer-coer (Cu-Cu), lead-lead (Pb-Pb), roton-lead (-Pb), and roton-roton (-) collisions at various energies, measured b the PHENIX, STAR, ALICE, and CMS Collaborations, b using the Tsallisstandard (Tsallis form of Fermi-Dirac or Bose-Einstein), Tsallis, and two- or three-comonent standard distributions which can be in fact regarded as different tes of thermometers or thermometric scales and seedometers. A central arameter in the three distributions is the effective temerature which contains information on the kinetic freeze-out temerature of the emitting source and reflects the effects of random thermal motion of articles as well as collective exansion of the source. To disentangle both effects, we extract the kinetic freeze-out temerature from the intercet of the effective temerature (T) curve as a function of article s rest mass (m 0 ) when lotting T versus m 0, and the mean transverse flow velocit from the sloe of the mean transverse momentum ( ) curve as a function of mean moving mass (m) when lotting versus m. Kewords: Kinetic freeze-out temerature, Transverse flow velocit, Transverse momentum sectrum, Tsallis-standard distribution, Tsallis distribution, Standard distribution PACS: q, Ag, Ld 1 Introduction Comaring with fixed target exeriments, the Relativistic Heav Ion Collider (RHIC) in the USA and the Large Hadron Collider (LHC) in Switzerland attract more studies due to their exhibitions on the evolution rocess of interacting sstem in collisions at higher energies. The exeriments at RHIC and LHC are more comlex and difficult for some limited conditions, which render ver limited measurable quantities. It is exected that more interacting information can be extracted from these limited measurable quantities. For examle, b analzing the transverse momentum sectra of final-state articles, one can obtain the kinetic freeze-out temerature of interacting sstems (emission sources) and the transverse flow velocit of roduced articles. These two quantities reflect the excited degree of emission sources and hdrodnamic exansion icture of interacting sstems at the stage of kinetic freeze-out, when hadrons are no longer interactive and their momenta do not change [1]. In high energ collisions, one can use different distribution laws to describe the transverse momentum sectra of final-state articles. In the framework of a multisource thermal or statistical model, we can use the standard (Boltzmann, Fermi-Dirac, or Bose-Einstein) distribution [2], the multi-comonent standard distribution, the Tsallis distribution (statistics), the Tsallis-standard (Tsallis forms of the standard) distribution [3], the Erlang distribution, the multi-comonent Erlang distribution [4], the Lév distribution [5, 6], the blast-wave function [7], the ower law, and so forth, to describe the article transverse momentum sectrum contributed b a given emission source. From the distributions mentioned above, one can directl extract the effective temerature, which is actuall not the real temerature (kinetic freeze-out temerature) of the emission source. The real temerature of the emission source should reflect urel thermal motion of articles in the source, while the effective temerature extracted from the transverse momentum sectrum includes both thermal motion and flow effect of articles, and is greater than the 1 fuhuliu@163.com; fuhuliu@sxu.edu.cn 2 Ro.Lace@Stonbrook.edu 1

2 real temerature. Onl if the effect of flow is excluded in the extraction of source temerature, can we obtain the kinetic freeze-out temerature. In this aer, within the framework of the multisource thermal model [8 10], we use the Tsallisstandard distribution, Tsallis distribution, and two- or three-comonent standard distribution [3] to describe the transverse momentum sectra of final-state light flavour articles roduced in gold-gold (Au-Au), coer-coer (Cu-Cu), lead-lead (Pb-Pb), roton-lead (-Pb), and roton-roton (-) collisions with different centralit intervals over a center-of-mass energ ( s NN ) range from 0.2 to 7 TeV. Although we can obtain the analtical exressions for the mentioned distributions, the Monte Carlo method is used to see the statistical fluctuations in the rocess of calculation. The results from the Monte Carlo method are comared with the exerimental data of the PHENIX [11], STAR [12], ALICE [13 18], and CMS Collaborations [19]. The kinetic freeze-out temerature of interacting sstem and transverse flow velocit of final-state articles are then extracted from the comarisons and analses. 2 The model and formulism In this aer, the transverse momentum sectra are described in the framework of the multisource thermal model [8 10], which assumes that man emission sources are formed in high energ collisions and are searated into a few grous resulting from different interacting mechanisms in the collisions as well as different event samles in exeriment measurements. The sources in the same grou have the same excitation degree and sta at a common local equilibrium state, which can be described b using different distribution laws. The final-state distribution is attributed to all sources in different grous, which results in a multi-temerature emission rocess if we use the standard distribution. This also means that the transverse momentum sectra can be described b a multi-comonent standard distribution which can be fitted b the Tsallis distribution. In fact, we can adot one-comonent Tsallis-standard (T-S) distribution (Tsallis form of the Boltzmann, Fermi-Dirac, or Bose-Einstein distribution), one-comonent Tsallis distribution, and two- or three-comonent standard distribution (Boltzmann, Fermi-Dirac, or Bose-Einstein distribution), to describe the transverse momentum sectra of final-state light flavour articles. According to refs. [3, 20], we use the uniform exressions of the Tsallis-standard distribution and Tsallis distribution which have more than one version [20 26], resectivel. Based on the invariant article momentum distribution, as well as the unit-densit function of transverse momentum ( ) and raidit(), the Tsallis-standard transverse momentum distribution and the Tsallis transverse momentum distribution are derived [3, 20]. In the resent work, the formalism of the Tsallis-standard transverse momentum distribution is adoted to be [3, 20] f T-S( ) = 1 N dn max = C T-S0 2 T d +m2 0 cosh T min {[ 1± q T-S 1 T T-S )] ± 1 } ( 2T +m20 cosh µ q 1 1 T-S +S d, (1) where C T-S0 is the normalization constant gv/(2π) 2 which results from 0 f T-S( )d = 1; g, V, N, µ, and m 0 are degenerac factor, volume, article number, chemical otential, and article s rest mass, resectivel; min is the minimum raidit and max is the maximum raidit; The + and in the ± sign are for 2 T +m2 0 cosh > µ and 2 T +m2 0cosh µ, resectivel; S has values 0, +1, and 1, which denote the Boltzmann, Fermi-Dirac, and Bose-Einstein distributions, resectivel; T T-S is the effective temerature of emission sources and q T-S is the entro index or nonequilibrium degree factor. In the calculation, since the effect of chemical otential can be ignored (i.e. µ = 0) in collisions at RHIC and LHC energies, we have 2 T +m2 0cosh > µ and take + in the ± sign. The final Tsallisstandard distribution is given b max {[ cosh ] 1 } 2T +m20 cosh q 1 1 T-S +S d. f T-S( ) = 1 dn = C T-S0 2 T N d +m q T-S 1 T min T T-S (2) Further, considering S = 0 and the ower index q T-S/(q T-S 1) and other limitations, the Tsallis transverse momentum distribution is obtained and can be written as [20] f T ( ) = 1 dn max [ = C T0 2 T N d +m2 0 cosh 1+ q ] q T T 1 2 T T min T +m2 0 cosh q T 1 d, (3) T 2

3 where C T0 is the normalization constant which gives 0 f T ( )d = 1, T T is the effective temerature, and q T is the entro index or nonequilibrium degree factor. Sometimes the uer index q T /(q T 1) is relaced b 1/(q T 1) due to q T being ver close to 1 and alication of mean field aroximation. The latter obtains a smaller q T. Particularl, in the resent work, we also use the two- or three-comonent standard distribution, which is different from the distributions introduced above. The standard Boltzmann, Fermi-Dirac, and Bose-Einstein distributions for the ith comonent (grou) can be uniforml shown as f i ( ) = 1 N dn max [ = C i0 2 T d +m2 0 cosh ex T min ( 2 T +m 2 0 cosh T i ) +S] 1 d, (4) where C i0 is the normalization constant which gives f 0 i ( )d = 1 and T i is the effective temerature for the ith comonent. In final state, the sectrum is contributed b the l comonents of the distribution; that is f S ( ) = 1 dn l = k i f i ( ), (5) N d where k i is the relative weight contributed b the ith comonent. This is the multi-comonent standard distribution. Considering the relative contribution of each comonent, we have the mean effective temerature to be l T S = k i T i, (6) i=1 where T S reflects the mean excitation degree for different comonents and can be used to describe the effective temerature of emission source. We would like to oint out that Eq. (6) is not a simle additive treatment for different temeratures, but an average weighted b different k i, where l i=1 k i = 1. It is exected that we can obtain the relation between the effective temerature T (T T-S, T T, or T S ) and the article s rest mass m 0. A linear fitting can obtain the intercet T 0 (T T-S0, T T0, or T S0 ) in linear relation T m 0. As the temerature corresonds to massless (m 0 = 0) article, T 0 is regarded as the source real temerature, or the kinetic freeze-out temerature of interacting sstem. According to refs. [11, 27 30], we have the relation between T and T 0, i=1 T = T 0 +am 0. (7) According to refs. [27, 28], the sloe a can be given b v 2 0 /2 and v 0 is the (average and transverse) radial flow velocit which is valid for low onl. Considering different distribution laws, v 0 can be v T-S0, v T0, or v S0, and T 0 can be T T-S0, T T0, or T S0, corresonding to Tsallis-standard, Tsallis, or standard distribution, resectivel. In the above discussions, different distribution laws are in fact regarded as different tes of thermometers or thermometric scale and seedometers. Since a = v 2 0/2 in Eq. (7) is onl valid for low region [27, 28], whereas the effective temerature is extracted in the resent work for fits in range which goes beond 2 GeV/c, we give u to extract radial flow v 0 from the sloe a in Eq. (7). Although other works [11, 29, 30] also regard the intercet T 0 in Eq. (7) as the kinetic freeze-out temerature, different relations between radial flow velocit v 0 and sloe a are used, and in some cases the relation is undetermined. In view of uncertain relations between radial flow velocit and sloe in wide range, we give u to extract radial flow velocit from the sloe, and extract onl the kinetic freeze-out temerature from T 0 in Eq. (7). We need an alternative method to extract mean (transverse) flow velocit. In our ver recent work [31], the linear relations between T and m 0 [Eq. (7)], mean transverse momentum and m 0, mean momentum and m 0, T and mean moving mass m, and m, as well as and m are studied. From the analses on dimension and quantit, we regard the intercet in T m 0 relation [Eq. (7)] as the kinetic freeze-out temerature and the sloe in m (or m) relation as the mean transverse flow velocit u T (or mean flow velocit u ). Particularl, for m and m relations, we have = 0 + u T m (8) and = 0 + u m, (9) 3

4 1/(2 )d 2 N/(dd ) ( -2 ) 1/(2 )d 2 N/(dd ) ( -2 ) 1/(2 )d 2 N/(dd ) ( -2 ) %(. ) Au-Au s NN = 0.2 TeV Tsallis-standard %( ) (e) 1/(2 )d 2 N/(dd ) ( -2 ) 1/(2 )d 2 N/(dd ) ( -2 ) 1/(2 )d 2 N/(dd ) ( -2 ) Tsallis standard %(. ) %( ) %( ) %( ) Figure 1. Transverse momentum sectra for π +, π, K +, K, (e), and (f) roduced in Au-Au collisions with different centralit bins at s NN = 0.2 TeV. The data are measured b the PHENIX Collaboration at midraidit [11] and are scaled verticall as quoted in the figure. The dotted, dashed, and solid curves are our results calculated b using the Tsallis-standard, Tsallis, and two-comonent standard distributions, resectivel. (f) where 0 and 0 denote the mean transverse momentum and mean momentum of massless article, resectivel. We would like to oint out that, in the case of describing meanwhile the same exerimental data, the three distributions result in nearl the same (or ), the same m (or m) relation, and the same u T (or u ). From the above analses, we see that the linear relations between T and m 0, and m, as well as and m [Eqs. (7) (9)] can be resectivel obtained from the same set of arameter values which are extracted from the same set of exerimental data. This means that the extraction rocesses of T 0, u T, and u are indeendent, though their values are entangled due to the same set of arameters. In articular, v 0 and the sloe in Eq. (7) are related to (transverse) flow velocit. However, there is no obvious and exact relation between them [11, 27 30]. Anhow, we think that the intercet in Eq. (7), the sloe in Eq. (8), and the sloe in Eq. (9) can rovide a set of alternative methods to extract the kinetic freeze-out temerature, transverse flow velocit, and flow velocit, resectivel [31]. Thus, we can use this set of alternative methods in the resent work to extract searatel the kinetic freeze-out temerature and transverse flow velocit according to transverse momentum sectra. This set of alternative methods is also used to judge different kinetic freeze-out scenarios in our another recent work [32] in which an evidence of mass-deendent differential kinetic freeze-out scenario is observed, while the single and double kinetic freeze-out scenarios are eliminated. 3 Results and discussion Figure 1 resents the centralit deendence of sectra for π +, π, K +, K, 4

5 1/(2 )d 2 N/(dd ) ( -2 ) 1/(2 )d 2 N/(dd ) ( -2 ) % ( x ) 20-30% ( x ) 40-60% ( x ) Cu-Cu s NN = 0.2 TeV K 0 S Cu-Cu /(2 )d 2 N/(dd ) ( -2 ) 1/(2 )d 2 N/(dd ) ( -2 ) Tsallis-standard Tsallis standard Cu-Cu Cu-Cu Figure 2. Transverse momentum sectra of KS 0, Λ, Ξ, and Ω roduced in Cu-Cu collisions at s NN = 0.2 TeV in three centralit intervals. The smbols reresent the exerimental data recorded b the STAR Collaboration in the raidi range < 0.5 [12]. The dotted, dashed, and solid curves are our results calculated b using the Tsallis-standard, Tsallis, and two- or three-comonent standard distributions, resectivel. (e), and (f) roduced in Au-Au collisions at center-of-mass energ s NN = 0.2 TeV. The data measured b the PHENIX Collaboration at midraidit [11] are reresented b different smbols which corresond to different centralit (C) bins from 0 5% to 80 92% scaled b different amounts as shown in the anels for clarit. The error bars are statistical onl. The dotted, dashed, and solid curves are our results calculated b using the Tsallis-standard, Tsallis, and two-comonent standard distributions, resectivel. The values of free arameters (T T-S, q T-S, T T, and q T ), normalization constants (N T-S0 and N T0 ) for comarisons between curves and data, and χ 2 er degree of freedom (χ 2 /dof) for Tsallisstandard and Tsallis distributions are given in Table 1. The values of free arameters (T 1, k 1, and T 2 ), mean effective temerature (T S ), normalization constant (N S0 ) for comarisons between curves and data, and χ 2 /dof for two-comonent standard distribution are listed in Table 2. One can see that all three tes of distribution laws are consistent with the exerimental data. The effective temeratures T T-S, T T, and T S increase with the increase of centralit or article mass, and T T-S T T < T S for a given set of data. We would like to oint out that the increase of centralit and the decrease of centralit ercentage have the same meaning. Figure 2 shows sectra of KS 0, Λ, Ξ, and Ω roduced in Cu-Cu collisions at snn = 0.2 TeV in different centralit intervals of 0 10%, 20 30%, and 40 60%. The smbols reresent the exerimental data recorded b the STAR Collaboration in the raidi range < 0.5[12]. For clarit, the results for different C intervals are scaled b different amounts shown in the anels. The uncertainties on the data oints include both statistical and sstematic errors. The dotted, dashed, and solid curves are our results based on the Tsallis-standard, Tsallis, and two- or three-comonent standard transverse momentum distributions, resectivel. The values of free arameters (T T-S, q T-S, T T, q T, T 1, k 1, T 2, k 2, and T 3 ), T S, normalization constants (N T-S0, N T0, and N S0 ), and χ 2 /dof are dislaed in Tables 1 and 2. Obviousl, the exerimental data can be described b the three tes of fit functions for in all centralit bins. The effective temeratures T T-S, T T, and T S increase with the increase of centralit or article mass, and T T-S T T < T S for a given set of data. The sectra of π +, K +,, and φ roduced in central (0 5%), semi-central (50 60%), and eriheral (80 90%) Pb-Pb collisions at s NN = 2.76 TeV are dislaed in Fig. 3, where N ev denotes the number of events. The data oints are measured b the ALICE Collaboration at midraidit ( < 0.5) [13, 14]. The error bars combined both statistical and sstematic uncertainties. The fitting 5

6 1/(2 N ev )d 2 N/(dd ) ( -2 ) 1/(2 N ev )d 2 N/(dd ) ( -2 ) % 50-60% 80-90% Pb-Pb s NN = 2.76 TeV /(2 N ev )d 2 N/(dd ) ( -2 ) d 2 N/(dd ) ( -1 ) Tsallis-standard Tsallis standard Figure 3. Transverse momentum sectra of π +, K +,, and φ roduced in Pb-Pb collisions at s NN = 2.76TeV in three centralitintervals. The smbols reresentthe exerimental datameasured b the ALICE Collaboration at mid-raidit ( < 0.5) [13, 14]. The fitting results with three tes of distributions are lotted b the curves. 1/(2 N ev )d 2 N/(dd ) ( -2 ) 1/(2 N ev )d 2 N/(dd ) ( -2 ) 1/(2 N ev )d 2 N/(dd ) ( -2 ) Pb s NN = 5.02 TeV 0<< % % % (e) K 0 S 1/(2 N ev )d 2 N/(dd ) ( -2 ) 1/(2 N ev )d 2 N/(dd ) ( -2 ) 10-6 Tsallis-standard Tsallis standard Figure 4. Transverse momentum sectra of π ±, +, K ±, Λ + Λ, and (e) KS 0 in -Pb collisions in different centralit classes at s NN = 5.02 TeV. The exerimental data reresented b the smbols are erformed b the ALICE Collaboration in 0 < < 0.5 [15]. Our fitting results are exhibited b the curves. 6

7 (1/N ev )d 2 N/(dd ) ( -1 ) (1/N ev )d 2 N/(dd ) ( -1 ) 1/(2 N inel. ev. )d 2 N/(dd ) ( -1 ) s =0.9 TeV - s =2.76 TeV - s =7 TeV Tsallis-standard Tsallis standard (e) (1/N ev )d 2 N/(dd ) ( -1 ) (1/N ev )d 2 N/(dd ) ( -1 ) (1/N inel )d 2 N/(dd ) ( -1 ) s =0.9 TeV - s =7 TeV - s =7 TeV Figure 5. Transverse momentum sectra of various identified hadrons roduced in - collisions. The smbols in Figs. 5, 5, and 5 reresent the exerimental data for π +, K +, and measured b the CMS Collaboration [19] in the range < 1 at s = 0.9,2.76,7 TeV, resectivel. The data oints of the ALICE Collaboration measured in < 0.8 at s = 0.9 TeV for Λ, φ, and Ξ + Ξ + [16], in < 0.5 at s = 7 TeV for Ξ and Ω [17], and in < 0.5 at s = 7 TeV for π ±, K ±, and + [18] are shown in Figs. 5, 5(e), and 5(f), resectivel. The fitting curves are our results. (f) results with Tsallis-standard, Tsallis, and two- or three-comonent standard distributions are lotted b the dotted, dashed, and solid curves, resectivel. The values of free arameters, T S, normalization constants, and χ 2 /dof are shown in Tables 1 and 2. These fitting functions describe the exerimental data. Once again, the effective temeratures T T-S, T T, and T S increase with increase of centralit or article mass, and T T-S T T < T S for a given set of data. Figure 4 exhibits sectra of π ±, +, K ±, Λ + Λ, and (e) KS 0 roduced in - Pb collisions in different centralit classes (0 5%, 40 60%, and %) at s NN = 5.02 TeV. The exerimental data reresented b different smbols are erformed b the ALICE Collaboration in the raidit interval 0 < < 0.5 [15]. The statistical and sstematic uncertainties are combined in the error bars. Our results from the analses of the Tsallis-standard, Tsallis, and two- or three-comonent standard distributions are exhibited b the dotted, dashed, and solid curves, resectivel. The values of free arameters, T S, normalization constants, and χ 2 /dof are summarized in Tables 1 and 2. As can be seen, the fitting results are in agreement with the exerimental data. The effective temeratures T T-S, T T, and T S increase with the increase of centralit or article mass, and T T-S T T < T S for a given set of data. Figure 5 gives sectra of various identified hadrons roduced in - collisions at different energies, where N inel denotes the number of inelastic events. The smbols in Figs. 5, 5, and 5 reresent the exerimental data for π +, K +, and measured b the CMS Collaboration [19] in the range < 1 at s = 0.9, 2.76, and 7 TeV, resectivel. The data oints of the ALICE Collaboration measured in < 0.8 at s = 0.9 TeV for Λ, φ, and Ξ + Ξ + [16], in < 0.5 at s = 7 TeV for Ξ and Ω [17], and in < 0.5 at s = 7 TeV for π ±, K ±, and + [18] are shown in Figs. 5, 5(e), and 7

8 T T-S T T Au-Au s NN = 0.2 TeV 0-5% 15-20% 30-40% 50-60% 70-80% 80-92% T S (e) m 0 (GeV/c 2 ) (f) m 0 (GeV/c 2 ) Figure 6. Particle mass and centralit deendences of T T-S, T T, and T S for ositive [left:,, and (e)] and negative articles [right:,, and (f)] in Au-Au collisions at s NN = 0.2 TeV. The smbols reresent the arameter values extracted from Fig. 1 and listed in Tables 1 and 2. The lines reresent linear fits to the results from each centralit bin as a function of mass using Eq. (7). 5(f), resectivel. The uncertainties corresonding to combined statistics and sstematics are shown as error bars. The fitting results (dotted, dashed, and solid curves) using mentioned functional forms (Tsallis-standard, Tsallis, and two- or three-comonent standard distributions) are suerimosed. The values of free arameters, T S, normalization constants, and χ 2 /dof are exhibited in Tables 1 and 2. Once more, the Tsallis-standard, Tsallis, and two- or three-comonent standard distributions can describe the exerimental data of the considered articles. The effective temeratures T T-S, T T, and T S increase with the increase of article mass, and T T-S T T < T S for a given set of data. To see clearl the deendences of effective temeratures T T-S, T T, and T S on centralit and article mass, the related values listed in Tables 1 and 2, which are extracted from figures 1, 2, 3, and 4, are analzed in Figs. 6, 7, 8 (left anel), and 8 (right anel), resectivel. Different smbols corresond to different centralit intervals shown in the anels. The lines are our fitted results. At the same time, Fig. 9 shows the deendences of effective temeratures on article mass for - collisions at the LHC energies. Different smbols corresond to different effective temeratures which are listed in Tables 1 and 2 and extracted from Fig. 5. The lines are our fitted results. Conclusions obtained from Figs. 1 5 and Tables 1 and 2 can be clearl seen from Figs In Figs. 6 9, the lines are worthwhile for further investigation. In fact, according to Eq. (7), the intercets at m 0 = 0 are real temeratures (kinetic freeze-out temeratures) T T-S0, T T0, and T S0, which are extracted from Tsallis-standard, Tsallis, and two- or three-comonent standard distributions, resectivel. We show the deendences of kinetic freeze-out temeratures on centralit C and energ snn in Fig. 10 for different sizes of collisions, where the related values for - collisions are comared with central (0 5%) nucleus-nucleus collisions. In most cases, T T-S0 T T0 < T S0. The kinetic freeze-out temeratures decrease slightl with the decrease of centralit (or with increase of centralit ercentage), and do not show an obvious change with energ and size for the central, semi-central, and eriheral collisions. In the considered energ range from 0.2 to 7 TeV, the mean kinetic freeze-out temeratures (... ) extracted from different distribution laws for the mentioned three centralities are marked in the anels b the secific values and dashed lines. We see that T T-S0 T T0 < T S0, and T 0 central > T 0 semi-central T 0 eriheral. The indeendences of kinetic freeze-out temeratures on energ and size render that different interacting sstems at the stage of kinetic freeze-out reached the same excitation degree of hadronic matter in the considered energ range. We noticed that the value of kinetic freeze-out temerature T S0 [(156±6) MeV] in central collisions 8

9 T T-S T T T S Cu-Cu s NN = 0.2 TeV K 0 S 0-10% 20-30% 40-60% m 0 (GeV/c 2 ) Figure 7. Same as Fig. 6, but for deendences of T T-S, T T, and T S on article mass and centralit in Cu-Cu collisions at s NN = 0.2 TeV. The smbols reresent the arameter values extracted from Fig. 2. T T-S Pb-Pb s NN = 2.76 TeV -Pb s NN = 5.02 TeV 0< T T-S T T % 50-60% 80-90% T T % 40-60% % T S (e) m 0 (GeV/c 2 ) T S (f) K 0 S m 0 (GeV/c 2 ) () ( ) Figure 8. Same as Fig. 6, but for article mass and centralit deendences of T T-S, T T, and T S for articles in Pb-Pb collisions at s NN = 2.76 TeV [left:,, and (e)], and in -Pb collisions at snn = 5.02 TeV [right:,, and (f)]. The smbols reresent the arameter values extracted from Figs. 3 and 4. 9

10 T T s = 0.9 TeV T T-S T T T S s = 7 TeV m 0 (GeV/c 2 ) m 0 (GeV/c 2 ) T T s = 2.76 TeV s = 7 TeV m 0 (GeV/c 2 ) m 0 (GeV/c 2 ) () Figure 9. Same as Fig. 6, but for deendences of T T-S, T T, and T S on article mass in - collisions at s = 0.9 TeV, s = 2.76 TeV, s = 7 TeV for ositive articles, and s = 7 TeV for charged articles. The smbols reresent the arameter values extracted from Fig. 5. T T-S0 T T0 T S (e) Cu-Cu,. 0.2 TeV Au-Au, 0.2 TeV ( +, +, ) Au-Au, 0.2 TeV ( -, -, ) Pb-Pb, TeV -Pb, 5.02 TeV -,. 0.9 TeV -, 2.76 TeV -, 7 TeV C (%) T 0 T 0 T T T-S0, - T T0, - T S0, - T T-S0, Cu-Cu, 0-10% T T0, Cu-Cu, 0-10% T S0, Cu-Cu, 0-10% T T-S0, Au-Au, 0-5% T T0, Au-Au, 0-5% T S0, Au-Au, 0-5% T T-S0, Pb-Pb, 0-5% central collisions T T-S0, Cu-Cu, 40-60% T T0, Cu-Cu, 40-60% T S0, Cu-Cu, 40-60% T T-S0, Au-Au, 50-60% T T0, Au-Au, 50-60% T S0, Au-Au, 50-60% T T-S0, Pb-Pb, 50-60% T T0, Pb-Pb, 50-60% semi-central collisions T T-S0, Au-Au, 80-92% T T0, Au-Au, 80-92% T S0, Au-Au, 80-92% T T-S0, Pb-Pb, 80-90% T T0, Pb-Pb, 80-90% T S0, Pb-Pb, 80-90% T T-S0, -Pb, % T T0, -Pb, % T S0, -Pb, % eriheral collisions T T0, Pb-Pb, 0-5% T S0, Pb-Pb, 0-5% T T-S0, -Pb, 0-5% T T0, -Pb, 0-5% T S0, -Pb, 0-5% T S0, Pb-Pb, 50-60% T T-S0, -Pb, 40-60% T T0, -Pb, 40-60% T S0, -Pb, 40-60% (f) s NN (TeV) Figure10. Centralit,energ,andsizedeendencesofkineticfreeze-outtemeraturesT T-S0, T T0, andt S0. The left anel [,, and (e)] is mainl for different kinetic freeze-out temeratures and centralities, and the right anel [,, and (f)] is mainl for three centralities and different energies. The smbols reresent our fitted results from Figs The results for - collisions are lotted for comarisons with central (0 5%) nucleus-nucleus collisions. The dashed lines marked with secific values in the right anel reresent the three averages T T-S0, T T0, and T S0 for three centralities, where T T-S0 T T0 < T S0. 10

11 [Fig. 10] is less than that of chemical freeze-out temerature (170 MeV) used in some theoretical estimations for the critical oint of QCD (quantum chromodnamics) [33 35]. Our result is also less than that of kinetic freeze-out temerature (177 MeV) extracted from an exonential shae of transverse mass sectrum [11], and equal to that of chemical freeze-out temerature (156 MeV) extracted from article ratios in a thermal and statistical model [36]. The latter item means that the kinetic freeze-out temerature obtained in the resent work is greater than that obtained from the thermal and statistical model [36], or the kinetic and chemical freeze-outs had haened at almost the same time. We are inclined to the latter. These results are not exactl consistent with each other even if we consider the kinetic freeze-out temerature being less than the chemical freeze-out temerature. It is undoubted that more analses are needed in the future. To extract the transverse flow velocit, we analze m relations in Fig. 11, where the related arameters listed in Tables 1 and 2 are used to calculate and m, where the three distributions result in the same relations. Different smbols corresond to different centralit intervals shown in the anels. The lines are our fitted results. From the sloes, we can obtain the mean transverse flow velocities of articles roduced in collisions with different centralities and at different energies. The related results are shown in Figs. 12 and 12 resectivel, where the results for - collisions are shown at C = 2.5% for comarisons. One can see that the mean transverse flow velocit in central collisions is greater than that in eriheral collisions. The relative difference is about 10%. The mean transverse flow velocit u T [(0.443±0.071)c] for central collisions obtained in the resent work is comarable with the radial flow velocit [( )c] [37, 38] of the blast-wave model [7, 37, 38], while the value of u T [(0.403±0.079)c] for eriheral collisions obtained in the resent work is far from the radial flow velocit ( 0) [37] of the blast-wavemodel [7, 37, 38]. In fact, u T contains isotroic radial flow and anisotroic collective flow (directed flow and others). In central collisions, the effect of anisotroic collective flow is small [39], u T contains mainl the isotroic radial flow. In eriheral collisions, the isotroic radial flow is nearl zero [37], u T contains mainl the anisotroic collective flow [39]. According to Fig. 12, u T in central collisions (which contains mainl the isotroic radial flow) is greater than that in eriheral collisions (which contains mainl the anisotroic collective flow). Because of u T containing the isotroic radial flow and anisotroic collective flow together, it is different from the original blast-wave model [7, 37, 38] in which onl the isotroic radial flow is contained. Our exlanation on flow comonents is comarable with the imroved blast-wave model [40, 41, 42] in which the third arameter is introduced to describe the anisotroic transverse flow generated in noncentral collisions [40] and the fourth arameter is introduced to take into account the anisotroic shae of the source in coordinate sace [41]. Although the deendence of ellitic flow on transverse momentum is described, the velocit of anisotroic transverse flow is not available in the imroved blast-wave model [40, 41, 42]. As an indirect measurement of model indeendence, the resent work is hoed to cause a further discussion on the flow velocit. The kinetic freeze-out temerature and transverse flow velocit are extracted from the transverse momentum sectra, though the methodolog used in the resent work is different from that of the blastwave model with the Boltzmann-Gibbs distribution or the Tsallis statistics [7, 37, 38], or the quark recombination model with the Cooer-Fre formula and flow [43]. In the latter two models, the kinetic freeze-out temerature and radial flow velocit can be obtained from certain formula descrition. In the resent work, we obtain the effective temerature from the formula of transverse momentum distribution, and kinetic freeze-out temerature from the intercet in linear T m 0 relation [Eq. (7)] [11, 27, 29 32] and the transverse flow velocit from the sloe in linear m relation [Eq. (8)] [32]. The method used in the resent work needs two stes to obtain the kinetic freeze-out temerature and transverse flow velocit, though the rocess of calculation is simle. Although the toic discussed b us was intensivel discussed in the literature and similar fitting and discussion were carried out in the literature man times, we have used an alternative methodolog [11, 27, 29 32] to extract the kinetic freeze-out temerature and transverse flow velocit, and the analsis is more extensive than anthing done before. This methodolog is different from that of the original blast-wave model [7, 37, 38] or the quark recombination model [43]. In addition, differing from the original blastwave model which studies the kinetic freeze-out temerature and isotroic radial flow, the resent work focuses on the kinetic freeze-out temerature and transverse flow which contains together the isotroic radial flow and anisotroic collective flow. We would like to emhasize that our recent work [31] investigated articularl the methodolog based on the analses of transverse momentum sectra, in which the linear relations between T and m 0, 11

12 < > < > < > Au-Au s 0-5% NN = 0.2 TeV 15-20% 30-40% +, K +, ) 50-60% 70-80% 80-92% m (GeV/c 2 ) Cu-Cu s NN = 0.2 TeV % 20-30% 40-60% m (GeV/c 2 ) -Pb s NN = 5.02 TeV 0-5% 50-60% % (e) m (GeV/c 2 ) < > < > < > Au-Au s 0-5% NN = 0.2 TeV 15-20% 30-40% ( -, K -, ) 50-60% 70-80% 80-92% m (GeV/c 2 ) Pb-Pb s NN = 2.76 TeV 0-5% 50-60% 80-90% m (GeV/c 2 ) TeV 2.76 TeV s = 7 TeV (f) m (GeV/c 2 ) Figure 11. Deendences of on m in 0.2 TeV Au-Au collisions for ositive articles, 0.2 TeV Au-Au collisions for negative articles, 0.2 TeV Cu-Cu collisions, 2.76 TeV Pb-Pb collisions, and (e) 5.02 TeV -Pb collisions with different centralities, as well as (f) - collisions at different energies. The smbols are our calculated results due to the arameter values listed in Tables 1 and 2. The lines are our fitted results. 0.6 Au-Au, 0.2 TeV ( +, +, ) Au-Au, 0.2 TeV ( -, -, ) Cu-Cu, 0.2 TeV Pb-Pb, 2.76 TeV -Pb, 5.02 TeV -, 0.9 TeV -, 2.76 TeV -, 7 TeV central collisions semi-central collisions eriheral collisions <u T > C (%) s NN (TeV) Figure 12. Deendences of u T on centralit and energ for different collisions marked in the anel. The smbols are the sloes in Fig. 11. The dotted, dashed, and dotted-dashed lines in Fig. 12 are the mean u T over energ for central, semi-central, and eriheral collisions, resectivel. 12

13 and m 0, and m 0, T and m, and m, as well as and m are studied. At the same time, the relations among some intercets and sloes are studied. The resent work focuses on an alication of the above methodolog. Different kinetic freeze-out temeratures are obtained due to different tes of thermometers or thermometric scales and the same transverse flow velocit is obtained due to the same transverse momentum sectrum. Our descritions on different thermometers or thermometric scales can be comared with those in thermodnamics. In addition, the other difference is that in our recent work [31], the Tsallis distribution and a one- or two-comonent Erlang distribution is used, while in the resent work, the Tsallis-standard, Tsallis, and two- or three-comonent standard distributions are used. The differences in the descrition of article sectra with the three different distributions are marginal and onl visible in a few cases at high momentum. This is also reflected in the similar values of the χ 2 /dof as listed in Tables 1 and 2 of the resent work. In our oinion, among different tes of thermometers or thermometric scales, the standard distribution and its multi-comonent distribution can be used to extract the standard temerature which is closest to the thermodnamic temerature in thermal hsics. The other thermometers can be regarded as non-standard thermometers which can be corrected if necessar. 4 Conclusions We summarize here our main observations and conclusions. The transverse momentum sectra of final-state light flavour articles (π ±, K ±, KS 0,,, φ, Λ, Λ, Ξ, and Ω, etc.), roduced in Au-Au, Cu-Cu, Pb-Pb, -Pb, and - collisions over an energ range from 0.2 to 7 TeV, are described b the Tsallis-standard, Tsallis, and two- or three-comonent standard distributions which are used in the multisource thermal model as different tes of thermometers or thermometric scales and seedometers. The characteristic of multile temeratures in multicomonent standard distribution can be described b the Tsallis-standard and Tsallis distributions. The calculated results are in agreement with the exerimental data recorded b the PHENIX, STAR, ALICE, and CMS Collaborations. Based on the three tes of thermometers or thermometric scales and seedometers, the effective temeratures (T T-S, T T, and T S ) and kinetic freeze-out temeratures (T T-S0, T T0, and T S0 ) of interacting sstem at the stage of kinetic freeze-out, and the mean transverse flow velocit of final-state articles are successivel extracted from the transverse momentum sectra. In the extraction of the kinetic freeze-out temerature and transverse flow velocit from the transverse momentum sectra, the methodolog used in the resent work is different from that of the original blast-wave model with the Boltzmann-Gibbs distribution or the Tsallis statistics, or the quark recombination model with the Cooer-Fre formula and flow. The effect of transverse flow is not directl considered in formulas describing the transverse momentum sectra. Instead, we obtain the effective temerature from the formulas, the kinetic freeze-out temerature from the intercet b lotting the effective temerature versus the article s rest mass, and the transverse flow velocit from the sloe b lotting the mean transverse momentum versus the mean moving mass. Our treatment disentangles naturall the random thermal motion of articles and collective exansion of sources. The resent work shows that the effective temeratures T T-S, T T, and T S increase with the increase of centralit or article s rest mass, and satisf the relationshi T-S T T < T S for a given set of data. In most cases, we have T T-S0 T T0 < T S0. The kinetic freeze-out temeratures decrease slightl with the decrease of centralit, and do not show an obvious change with energ and size for the central, semi-central, and eriheral collisions. The mean kinetic freeze-out temeratures derived from the three tes of thermometers show T T-S0 T T0 < T S0, and T 0 central > T 0 semi-central T 0 eriheral. The indeendences of kinetic freeze-out temeratures on energ and size reveal that different interacting sstems at the stage of kinetic freeze-out reached the same excitation degree of hadronic matter in the considered energ range. Generall, the effective temerature is greater than kinetic freeze-out temerature, since the former containing both the thermal motion and flow effect while the latter containing onl the thermal motion. The resent work also shows that the central collisions have a higher kinetic freeze-out temerature than that of eriheral collisions. This is consistent with those of the chemical freeze-out temerature and effective temerature. Based on the three temeratures we conclude that the excitation degree of central collisions is higher than that of eriheral collisions at both the chemical and kinetic freeze-out stages. In addition, at the stage of kinetic freeze-out, the transverse flow velocit in central collisions is larger than that in eriheral collisions, indicating a larger squeeze and exansion in central collisions 13

14 which also show a higher temerature and excitation degree. Thus, the kinetic freeze-out temeratures, which describe the random thermal motions of articles, and the transverse flow velocities, which describe the collective exansions of sources, aint a consistent icture. (e) Theresentworkis analicationofourrecentwork[31]inwhich analternativemethod isinvestigated to extract searatel the kinetic freeze-out temerature and flow velocit. If we regard the standard distribution and its multi-comonent distribution as the most exact and standard thermometer due to the being closest to the thermodnamic temerature in thermal hsics, we obtain a kinetic freeze-out temerature [(156 ± 6) MeV in central collisions] which is equal to the chemical freeze-out temerature extracted from article ratios [36]. This indicates that the two freeze-outs had haened at almost the same time even at the RHIC and LHC energies. We would like to treat other two thermometers used in the resent work as non-standard thermometers, which can be corrected if necessar. Another alication [32] of our alternative method shows an evidence of mass-deendent differential kinetic freeze-out scenario, while the single and double kinetic freeze-out scenarios are eliminated. (f) In central collisions, the transverse flow mainl contains the isotroic radial flow, since the effect of anisotroic collective flow is small and can be neglected in the extraction of transverse flow which contains mainl the isotroic radial flow. In contrast, the transverse flow of eriheral collisions includes mainl the anisotroic collective flow, since the radial flow is nearl zero. In this work, the transverse flow contains both the isotroic radial flow and anisotroic collective flow, which outerforms the original blast-wave model since it onl includes the radial flow. We also found that the radial flow in central collisions is greater than anisotroic collective flow in eriheral collisions. Acknowledgments This work was suorted b the National Natural Science Foundation of China under Grant No and the US DOE under contract DE-FG02-87ER40331.A008. References [1] Guta N 2008 A Stud of fluctuations in multilicit distributions at ultrarelativistic heav ion interactions, Ph.D. Thesis, Universit of Jammu, India, htts:// thesis.df. [2] Meng C-R 2009 Chin. Phs. Lett [3] Liu F-H, Gao Y-Q and Wei H R 2014 Adv. High Energ Phs [4] Gao Y-Q, Tian C-X, Duan M-Y, Li B-C and Liu F-H 2012 Pramana J. Phs [5] Wilk G and Wldarczk Z 2000 Phs. Rev. Lett [6] Adams J et al. (STAR Collaboration) 2005 Phs. Rev. C [7] Schnedermann E, Sollfrank J and Heinz U 1993 Phs. Rev. C [8] Liu F-H 2004 Phs. Lett. B [9] Liu F-H 2008 Phs. Rev. C [10] Liu F-H 2008 Nucl. Phs. A [11] Adler S S et al. (PHENIX Collaboration) 2004 Phs. Rev. C [12] Agakishiev G et al. (STAR Collaboration) 2012 Phs. Rev. Lett [13] Abelev B et al. (ALICE Collaboration) 2013 Phs. Rev. C [14] Abelev B et al. (ALICE Collaboration) 2015 Phs. Rev. C [15] Abelev B et al. (ALICE Collaboration) 2014 Phs. Lett. B

15 [16] Aamodt K et al. (ALICE Collaboration) 2011 Eur. Phs. J. C [17] Floris M et al. (ALICE Collaboration) 2011 J. Phs. G [18] Adam J et al. (ALICE Collaboration) 2015 Eur. Phs. J. C [19] Chatrchan S et al. (CMS Collaboration) 2012 Eur. Phs. J. C [20] Clemans J and Worku D 2012 Eur. Phs. J. A [21] Biró T S 2009 Eur. Phs. J. A [22] Biró T S, Purcsel G and Ürmöss K 2009 Eur. Phs. J. A [23] Tsallis C 1988 J. Stat. Phs [24] Tsallis C 2009 Braz. J. Phs [25] Biró T S 2001 Phsica A [26] Biró T S, Barnaföldi G G, Ván P and Ürmöss K 2014 arxiv: [he-h]. [27] Takeuchi S, Murase K, Hirano T, Huovinen P and Nara Y 2015 Phs. Rev. C [28] Heiselberg H and Lev A M 1999 Phs. Rev. C [29] Heinz U W 2004 arxiv:he-h/ [30] Russo R 2015 Measurement of D + meson roduction in -Pb collisions with the ALICE detector, Ph.D. Thesis, Universitaâ degli Studi di Torino, Ital, arxiv: [nucl-ex]. [31] Wei H-R, Liu F-H and Lace R A 2016 Eur. Phs. J. A [32] Lao H-L, Wei H-R, Liu F-H and Lace R A 2016 Eur. Phs. J. A [33] Fodor Z and Katz S D 2004 JHEP [34] Datta S, Gavai R V and Guta S 2013 Nucl. Phs. A c. [35] Xu N (for the STAR Collaboration) 2014 Nucl. Phs. A [36] Stachel J, Andronic A, Braun-Munzinger P and Redlich K 2014 J. Phs. Conf. Ser [37] Tang Z B, Xu YC, Ruan L J, van Buren G, Wang F Q and Xu Z B 2009Phs. Rev.C79, (R). [38] Shao M, Yi L, Tang Z B, Chen H F, Li C and Xu Z B 2010 J. Phs. G [39] Shi S S (for the STAR Collaboration) 2016 J. Phs. Conf. Ser [40] Huovinen P, Kolb P F, Heinz U, Ruuskanen P V and Voloshin S 2001 Phs. Lett. B [41] Adler C et al. (STAR Collaboration) 2001 Phs. Rev. Lett [42] Sun X, Masui H, Poskanzer A M and Schmah A 2015 Phs. Rev. C [43] Ürmöss K and Biró T S 2010 Phs. Lett. B

16 Table 1. Values of free arameters, normalization constants, and χ 2 /dof corresonding to Tsallis-standard and Tsallis distributions in Figs The values of χ 2 /dof for Ξ + Ξ + in Fig. 5 are in fact the values of χ 2 values due to less data oints. The temeratures are in the units of GeV. 16 Figure Te T T-S q T-S N T-S0 χ 2 /dof T T q T N T0 χ 2 /dof 1 0-5% ± ± ± ± ± ± % ± ± ± ± ± ± % ± ± ± ± ± ± % ± ± ± ± ± ± % ± ± ± ± ± ± % ± ± ± ± ± ± % ± ± ± ± ± ± % ± ± ± ± ± ± % ± ± ± ± ± ± % ± ± ± ± ± ± % ± ± ± ± ± ± % ± ± ± ± ± ± % ± ± ± ± ± ± % ± ± ± ± ± ± % ± ± ± ± ± ± % ± ± ± ± ± ± % ± ± ± ± ± ± % ± ± ± ± ± ± % ± ± ± ± ± ± % ± ± ± ± ± ± % ± ± ± ± ± ± % ± ± ± ± ± ± % ± ± ± ± ± ± % ± ± ± ± ± ± (e) 0-5% ± ± ± ± ± ± % ± ± ± ± ± ± % ± ± ± ± ± ± % ± ± ± ± ± ± % ± ± ± ± ± ± % ± ± ± ± ± ± (f) 0-5% ± ± ± ± ± ± % ± ± ± ± ± ± % ± ± ± ± ± ± % ± ± ± ± ± ± % ± ± ± ± ± ± % ± ± ± ± ± ± % ± ± ± ± ± ± % ± ± ± ± ± ± % ± ± ± ± ± ± % ± ± ± ± ± ± % ± ± ± ± ± ± % ± ± ± ± ± ± % ± ± ± ± ± ± % ± ± ± ± ± ± % ± ± ± ± ± ± % ± ± ± ± ± ± % ± ± ± ± ± ± % ± ± ± ± ± ±