arxiv: v1 [nucl-th] 14 Feb 2018

Transcription

1 Chemical freezeout arameters within generic nonextensive statistics Abdel Nasser TAWFIK Egytian Center for Theoretical Physics (ECTP), Modern University for Technology and Information (MTI), 57 Cairo, Egyt World Laboratory for Cosmology And Particle Physics (WLCAPP), 57 Cairo, Egyt arxiv: v [nuclth] 4 Feb 208 Hayam Yassin and Eman R. Abo Elyazeed Physics Deartment, Faculty of Women for Arts, Science and Education, Ain Shams University, 577 Cairo, Egyt Abstract The article roduction in relativistic heavyion collisions seems to be created in a dynamically disordered system which can be best described by an extended exonential entroy. In distinguishing between the alicability of this and BoltzmannGibbs (BG) in generating various articleratios, generic (non)extensive statistics (GNS) is introduced to the hadron resonance gas model. Accordingly, the degree of (non)extensivity is determined by the ossible modifications in the hase sace. Both BG extensivity and Tsallis nonextensivity are included as very secial cases defined by secific values of the equivalence classes (c,d). We found that the article ratios at energies ranging between 3.8 and 2760 GeV are best reroduced by nonextensive statistics, where c and d range between 0.9 and. The resent wor aims at illustrating that the roosed aroach is well caable to manifest the statistical nature of the system on interest. We don t aim at highlighting deeer hysical insights. In other words, while the resulting nonextensivity is neither BG nor Tsallis, the freezeout arameters are found very comatible with BG and accordingly with the wellnown freezeout hasediagram, which is in an excellent agreement with recent lattice calculations. We conclude that the article roduction is nonextensive but should not necessarily be accomanied by a radical change in the intensive or extensive thermodynamic quantities, such as internal energy and temerature. Only, the two critical exonents defining the equivalence classes (c, d) are the hysical arameters characterizing the (non)extensivity. Keywords: Quantum statistical mechanics, article roduction in relativistic collisions, robability theory PACS Nos: d, Dw, Cw Introduction A comrehensive descrition for various henomenon in different highenergy exeriments [, 2] can be given by extensive [such as BoltzmannGibbs (BG)] and nonextensive (such as Tsallis) statistics. In both statistical aroaches, the thermodynamic consistency should be guaranteed. In roving that the fireballs or heavy resonances lead to a bootstra aroach, i.e. further fireballs consist of smaller fireballs and so on, extensive statistics was utilized by Hagedorn [3]. Assuming that the characteristic distributionfunction gets variations from a ossible symmetrical change, etc., nonextensive concet was first introduced for the article roductions [4]. The imlementation of nonextensive Tsallis statistics was first introduced in Refs. [5, 6, 7, 8], with a clear emhasize that the hase sace lays an essential role. It has been argued that substituting the Boltzmann factor by qexonential function with q > leads to a good agreement with the exerimental results at high energies, esecially the transverse momentum sectra. Is is widely acceted that the statistical fits of the the transverse momentum sectra characterize the inetic freezeout, which is conjectured to tae lace later (in the sense that the system cools down) after the chemical freezeout []. Recently, Tawfi exlained that this rocedure is not necessarily fully incororating the nonextensivity in the Corresonding author: a.tawfi@eng.mti.edu.eg

2 article roductions, even at the inetic freezeout [9, 0,, 2, 3]. The longrange fluctuations, the correlations, and the interactions besides the ossible modifications in the hase sace of the article roduction are not roerly incororated through Tsallis algebra. In longrange interactions, both thermodynamic and longtime limits do not commute. Therefore, generic nonextensive statistics (GNS) was introduced in Ref. [9, 0, ], in which the hase sace becomes resonsible in determining the degree of (non)extensivity. It was shown that the lattice thermodynamics is well reroduced when the roosed GNS become characterized by extensive critical exonents (, ), while the heavyion article ratios are only reroduced when the roosed GNS become nonextensive critical exonents, e.g. neither 0 nor. The latter differs from Tsallis [9, 0, ]. Nonextensive statistics becomes the relevant aroach for nonequilibrium stationary states. While zeroth law of thermodynamics in equilibrium introduces the temerature, a socalled hysical temerature was roosed when utilizing Tsallisalgebra, see for instance [4, 5, 6, 7]. It was concluded that if the inverse Lagrange multilier associated with constrained internal energy is regarded as the temerature, both Tsallis and Clausius entroies become identical. This temerature is believed to differ from the hysical one. Based on this assumtion, the hysical temerature was conjectured to be extraolated to the real temerature, as it is assumed that this quantity is identical to the freezeout temerature. More details shall be elaborated in section 2.3. It was concluded that, the assumtion on inverse Lagrange multilier[4] follows the same ath that nonextensivity. It is only fulfilled through modifications in the thermodynamic intensive and extensive quantities, such as internal energy and temerature. In the resent wor, we introduce a GNS aroach, which assumes that the hase sace determines the degree of (non)extensivity deending on two critical exonents defining equivalence classes (c, d), while the intensive and extensive thermodynamic quantities shouldn t necessarily to be modified. To summarize, the(c, d)temerature seems not distinguished from the T arameter of Tsallis statistics. In the resent aer, generalized thermostatistics (GTT) is alied to investigate the chemical freezout. The main line of the critical remars on Tsallis tye statistics is related to the interretation of the resulting temerature. The temerature is a well elaborated concet in GTT. Characterizing the statistical nature of the article roduction is the main focus of the resent wor. As introduced, the roosed GNS aroach is able to find out either the article roduction is to be described by extensive or nonextensive statistics. This basically differs from an ad hoc imlementation of BG and Tsallis aroaches. In Ref. [9], one of the authors (AT) confronted calculations from the hadron resonance gas (HRG) model with GNS to the corresonding lattice QCD thermodynamics and to various article ratios measured at 7.7 and 200 GeV. He concluded that the lattice QCD simulations are best reroduced when the scaling exonents (c, d) get the extensive values, i.e. both are unity, while the article roduction seems to have nonextensive classes, i.e. both differ from unity. In other words, the roosed GNS aroach seemed to suggest that the lattice calculations are extensive while the article ratios are stemming from nonextensive rocesses. Let us first emhasize that this conclusion is fully correct, as the lattice QCD simulations are actually based on extensivity (additivity) of subsystems (quars and glouns roceeded and communicated by CPUs). On the other hand, the nature of the article roductions in the relativistic heavyion collisions might be different. The resentworis devotedto a detailed analysisofthe articleroduction in a wide rangeofbeam energies. Accordingly, the statistical aroach whether extensive or nonextensive, which well reroduces the chemical freezeout arameters, shall be defined. For the sae of comleteness, we highlight that the lattice freezeout arameters are found comatible to the ones deduced from extensive fits of HRG calculations to various article ratios measured in a wide range of beam energies. This is a solid confirmation that the roosed statistical aroach, GNS, wors well for the article roductions. The resent aer is organized as follows. The roosed aroaches shall be shortly discussed in section 2. A reminder to GNS shall be given in section 2.. Sections 2.2 and 2.3 shall be devoted to the characterizations of (c,d) nonoextensiveentroy and the extraolation of T q to T ch, resectively. GNS in HRG model shall be detailed in section 2.4. In section 3., the statistical fits of various articleratios at different beamenergies shall be outlined. The freezeout arameters, the temeratures and the baryon chemical otentials, deduced from the statistical fits of GNS grandcanonical artition function to various articleratios measured in heavyion exeriments shall be determined and comared with the results deduced from extensive HRG calculations. The resulting equivalence classes (c, d) shall be discussed in 3.2. Section 4 outlines the conclusions and the final remars. 2

3 2 Aroaches To exlore the hase diagram of hot and dense nuclearmatters which can be estimated at the beamenergies 3.8, 4.2, 5, 6.4, 7.7, 8, 9,.5, 2, 7, 9.6, 27, 39, 62.4, 30, 200, and 2760 GeV lay an essential role. It was assumed that this can be accomlished from the statistical fit of thermal models, such as HRG, in which GNS aroach is imlemented [9, 0, ] to various articleratios measured in relativistic heavyion collisions. A detailed descritions for the generic nonextensive statistics can be found at [9, 0, ]. In the section that follows we give a short reminder to the generic nonextensive statistics. 2. Reminder to generic nonextensive statistics Shannon entroy (S ) is uniquely characterized by four axioms, Khinchin axioms:. SK, continuity: for any n N, S () is continuous with resect to n, i.e. the entroy deending on, contentiously. Accordingly, function g() becomes continuous, as well. 2. SK2, maximality: for given n N and for (,, n ) n, S (,, n ) reaches maximum at i = /n, where i =,,n. Consequently, for equidistribution (), the entroy becomes maximum and the function g() become a concave function. 3. SK3, exandability: S (,, n,0) = S (,, n ), i.e. adding a nonzero robability state does not change the entroy so that g(0) = SK4, generalized additivity: if ij 0, i = m i j= ij, and n i= i =, where i =,,n and j =,,m i, then S (,, nmn ) = S (,, n ) n i= is ( i i, im i i ). This means that the entroy of a system, which is divided into two subsystems A and B, is given as S A lus the exectation value of S B conditional on A. Fulfilling the four axioms characterizes(c, d) extensive entroy, where the equivalence classes (, ) recall BG statistics. Violating the fourth axiom, individually, characterizes (c, d) nonextensive entroy. The latter can be Tsallis, where (q,0), or GNS, where (c,d) tae other values. We recall again that, the scaling exonents c and d are corresonding to two indeendent asymtotic roerties of the entroy functional (s), namely [8, 9] (a) d z c = lim x 0 s(zx) s(x), () = lim x 0 s(x a ) x ac s(x), (2) with a and z are arbitrary variables ranging between 0 and, but not affecting the (non)extensivity. 2.2 Characterizations of (c,d) nonextensive entroy From formed rocesses, on which the hase sace deends, the extensive entroy can be comuted, where the hase sace linearly increases with the increase in the system size (N). On the other hand, for deformed rocesses, one needs a general reresentation of the entroy functional (s). Generalized families of classreresentation can be achieved through twotime differentiable, monotonically increasing functions. For instance, a generalization for Shannon entroy was roosed as H() = W i= h( i), where W gives the number of states of the hasesace rocesses [8]. From the asymtotic roerties of the traceform entroy functionals, following reresentations have been roosed [8]. Greresentation: G() = W i= g( i), where and is a generalized logarithm, see Eq. (6). g( i ) i ( i ), (3) 3

4 2. Sreresentation: From i 0 dylog(y) = ilog( i ) i, the second term in righthand side can be absorbed in α( i ), when the normalization constrains of the maximumentroyrincile is imlemented. α is a Lagrangian multilier. Accordingly, h( i ) = i log( i ) and s( i ) = i Both reresentations should give equivalent (c,d) entroy classes 0 dx(x). (4) s (c,d,r) (x) = r c A d ex(a)γ(d,a clog(x)) rx, (5) where A = cdr/( r( c)). Corresondingly, both exonential and logarithm should be generalized. The earlier shall be introduced in Eq. (2), while the latter reads ( (c,d,r) (x) = r [ x c r( c) ) ] d log(x). (6) dr Now we are left with the fifth roerty of (c,d) nonextensive entroy. This should be characterized by vanishing entroy roerty, 0 = S (c,d,r) () = b(c,d,r) 0 dx (c,d,r) (x), (7) W s (c,d,r) (b (c,d,r) i ), (8) i= where b (c,d,r) otimizes the asymtotic extensivity. When scaling s (c,d,r), even the reresentation given in Eq. (3) fulfils vanishingentroy roerty, as well. To summarize, besides SK, SK2, and SK3, the (c,d) nonextensive entroy is characterized by: traceform of the entroy functional, S() = W i= s( i), and vanishing entroy roerty, Eqs. (7)(8). 2.3 Extraolation of T q to T ch An ab initio assumtion was made that the Tsallis freezeout temerature can be extraolated to the BGtemerature (T ch ) deduced from extensive thermalstatistical models [20, 2], T q = T ch (q ), (9) where q is Tsallis nonextensive arameter, i.e. q >. is conjectured to deend on the energy transfer between a source and the surrounding. It was found that at T ch = 92±5 MeV, = (950±0) MeV [22]. As mentioned in section, The temerature deduced from the statistical fits of the transverse momentum distributions was interreted as the inetic freezeout temerature. To remain within the scoe of the resent aer, we leave this oint without any further discussion. In a classical ideal gas, the Tsallis entroy, which describes an exact entroy for microcanonical distribution of such gas, was roosed to relate Tsallis temerature, the qtemerature, (T q ) to the BG freezeout temerature (T ch ) [23] ( ) Sq T q = T ch ex, (0) C where C is the heat caacity of reservoir system and S q is Tsallisentroy. To summarize, the qtemerature, which as discussed in section is now as hysical temerature, is about 2 3 times smaller than T ch now as the temerature [24]. Accordingly, the resulting freezeout hasediagram largely differs from the one drawn from the extensive statistical fits of various measured 4

5 articlesratios and recent lattice QCD simulations, esecially the freezeout temerature [24]. Desite the wellnow signroblem in the lattice QCD calculations at finite chemical otential (µ b ), where the Monte Carlo techniques liely fails, both deconfinement and freezeout temeratures are conjectured to be comatible with each other, esecially at small chemical otentials [25]. This range of µ b is to be related to RHIC and LHC energies, where recise measurements for different articleratios are available[26, 27]. What is obtained sofar is that at vanishing µ b and at the corresonding T q, the resulting HRG thermodynamics, even when alying Tsallisalgebra, considerably differs from the one based on the firstrincile lattice calculations. This can be understood from the fact that the latter assume extensivity and additivity besides an overall thermal equilibrium. In light of this, it is aarent that the Tsallistye nonextensivity shouldn t be utilized in order to reroduce of the lattice calculations. In this regards, we first remind with the remarable success of HRG with BG statistics in reroducing the firstrincile lattice thermodynamics [28, 29, 30, 3, 32, 33, 34]. This is another argumentation of why the roosed GNS aroach is the roer one. It reroduces the lattice calculations only when the equivalence classes (c, d) get extensive values, i.e. (, ). The success in characterizing the statistical nature of the lattice calculations manifests the fact that the roosed aroach (GNS) is generic and its alicability is much wider than that of BG and Tsallis. As discussed, comrehensive wors have been conducted in order to bring hysical temerature closer to the the temerature, for instance Eqs. (9)(0). These wors assume that the nonextensivity in the article roduction can be characterized through a radical change in intensive or extensive thermodynamic quantities, such as the internal energy and the temerature. The resulting T q, which is lower than T ch, are extraolated. One of authors (AT) has first introduced GNS aroach to the statistical fits of article ratios at 200 and 2760 GeV [9]. The resent wor covers a wider range of beam energies. Accordingly, the two critical exonents defining the equivalence classes (c, d) can be determined. While both intensive or extensive thermodynamic quantities remain almost unchanged, (c, d) unambiguously characterize the nonextensive nature of the article roduction. In a future wor, we shall study the energy deendence of (c,d). 2.4 GNS in HRG model As introduced in Ref. [9, 0, ], both extensive and nonextensive statistical roerties can be determined by two scaling exonents defining equivalence classes (c, d) for both correlated and uncorrelated systems, for instance. In its thermodynamic limit, the statistical nature of the system of interest can be characterized, unambiguously. The remarable success of the thermal models (extensive) [, 2] in describing the article roduction and that of the Tsallis nonextensive algebra [5, 6] simly mean that the cost which should be aid, namely the radical reduction in the freezeout temerature, was not necessarily. Such an ad hoc imlementation of statistical aroaches is not the best way to decide whether the article ratios or transverse momentum distributions or any other henomenological sectra are extensive or nonextensive. A direct confrontation with GNS is this [35]. Here, the system is not enforced to be biased towards either extensive or nonextensive statistical descrition. With a generalized exonential function, Eq. (2), ln Z(T) = ±V N M B i g i (2) 3 0 ln[±ε c,d,r (x i )] d 3, () where V is the fireball volume and x i = [µ i E i ()]/T with E i () = 2 m 2 i being the disersion relation of ith state (article) and ± reresent fermions and bosons, resectively. The exonential function ε c,d,r (x i ) is generalized as [8, 9] { ε (c,d,r) (x) = ex d ( [W B( x/r) /d) W (B)] }, (2) c a where W is Lambert Wfunction which has real solutions at = 0 with d 0 and at = with d < 0, and B = ( c)r [ ] ( c)r ( c)r ex, (3) ( c)r where c and d are two critical exonents defining equivalence classes for all interacting and noninteracting systems. They give estimations for two scaling functions with two asymtotic roerties. r is almost a free arameter. It is assumed not affecting the (c, d)class. 5

6 for d > 0, r < /( c), for d = 0, r = /( c), and for d < 0, r > /( c). Alternatively, a articular function for r, namely r = ( ccd) was roosed [8, 9]. From the artition function, Eq. (), the thermodynamical roerties such as ressure () and number density (n) [36] can be derived, = N M B i= N M B n = ± i= g i T g i ln[±ε c,d,r (x i )] 2 d, (4) ] ε c,d,r (x i ) W 0 [B( xi r ) d [ ]) ( c)[±ε c,d,r (x i )](r x i ) (W 2 d. (5) 0 B( xi r ) d It should be remared that in Eq. (), the logarithmic function should also be generalized. This was given in Eq. (6), when assuming Greresentation. In the resent calculations, the ossible decay channels of heavy resonances are taen into consideration. The th article final number density (identical exressions for other thermodynamic quantities, such as ressure, entroy and energy density, can be deduced) is given as = n n final N M B l b l n l, (6) where b l is the effective branching ratio of lth hadron resonance into the th article of interest. It is noteworthy mentioning that the resonance decays might be seen as correlations and interactions among the roduced articles and thus the introduction of nonextensivity to the article roduction in final state becomes eligible. 3 Results 3. Chemical freezeout hasediagram The article ratios measured by the E866, NA44, NA49, and NA57 exeriments at the Suerroton Synchrotron (SPS) [37, 38, 39, 40] at 3.8, 4.2, 5, 6.4, 8, 9, 2, and 7 GeV and by the STAR exeriment at the Relativistic Heavy Ion Collider (RHIC) at 7.7,.5, 9.6, 27, 39, 30, and 200 GeV shall be comared with the HRG calculations. Together with the ALICE measurements, we construct a data set covering beam energies from 3.8 u to 2760 GeV. For the same set of article ratios, we fit our calculations based on GNS, section 2.4, Eq. (5), to the exerimental results. The free arameters are T ch and µ b, besides the scaling exonents (c, d). There values are deduced where chemical freezeout conditions, such as s/t 3 = 7 [, 4, 42, 43, 44, 45, 46, 47] is fulfilled and minimum χ 2 is obtained, simultaneously. The results are shown in Figs., 2, and 3. Results on the statistical fits for various articleratios measured by ALICE at 2760 GeV (LHC) [left anel (a)] and by STAR exeriment at 200 GeV (RHIC) [right anel (b)] are deicted in Fig.. The resulting T ch and µ b are summarized in Tab. and shall be deicted in Fig. 4. It is obvious that they are very comatible with the ones deduced from BG statistics [26]. But, the resulting exonents (c,d) for first and second nonextensive roerty, resectively, differ from BG (, ) and Tsallis (q, 0). Sofar, we conclude that the resulting T ch and µ b seem very comatible with the ones assuming extensive statistics (BG). Does this alone mean that the article roductions in heavyion collisions is based on extensive rocesses? The answer is aarently no [9, 0, ]. A ey measure would be the reroduction of firstrincile calculations, the lattice QCD simulations on thermodynamics, for instance. While HRG with BG statistics excellently models the lattice thermodynamics below the critical temerature and describes well the article roductions in heavyion collisions [28, 29, 30, 3, 32, 33, 34], the Tsallis nonextensivity doesn t do this job. Only, it is conjectured to characterize the inetic freezeout and fitted well with the 6

7 measured transversemomentum distribution [5, 6]. How to unify these together, esecially, when the Tsallis temerature is remarably smaller than the freezeout temerature? First, let us recall that the Tsallis temerature fails to reroduce the lattice thermodynamics [20, 2]. The alication of the GNS to lattice QCD thermodynamics and article ratios which measured at different energies of RHIC was already done by one of the authors (AT) in Ref. [48]. Second, the answer has been delivered, coule years ago. In revious sections, we discussed on this that, the Tsallis resulting temerature should be extraolated to the freezeout one, T ch. The resent wor rooses another answer. Figure comares exerimental results on various articleratios (circles), /, K /K, /, /, /, /, K /, K /, /, /, /, /, and / measured by ALICE and STARexeriment at energies 2760 GeV and 200 GeV, left and right anel, resectively, are statistically fitted by means of the HRG model in which GNS is imlemented. In our HRG calculations, all ossible decay channels yielding the articles of interest and the related branching ratios are taen into account, Eq. (6). For the decay channels with notyetmeasured robabilities, the rules given in Ref. [40] have been alied. But no finitesize correction was alied [49]. We sum u contributions from all hadron resonances listed in recent article data grou comilation with masses 2 GeV. This refers to 388 different states of mesons and baryons besides their antiarticles. For further details, interested readers can consult Ref. []. The number density can be derived from the artition function and accordingly the article ratios, Eq. (5), can be determined. Other GNS fits are illustrated in Fig. 2 (RHIC) and Fig. 3 (SPS). In Fig. 4, the freezeout arameters, T ch and µ b (closed circles) as deduced from the statistical fits of the genericnonextensive HRG are comared with the ones from the extensive HRG (dashed line). At a finite chemical otential, the freezeout temerature in both cases is determined at constant s/t 3, where s is the entroy density []. Aarently, both aroaches are very comatible with each other. In other words, there is almost no difference between the resulting hysical temerature and the temerature (T ch ). We also comare with other freezeout arameters (symbols) determined in different henomenologies; Andronic et al. [40], Tawfi et al. [26, 27], and UrQMD [50]. With the latter we mean various simulations by hybrid UrQMD version 3.4 at varying chemical otentials. At each value of the chemical otential, the simulated article ratios are fitted by means of extensive HRG. In other words, UrQMD simulations in this case are taen as exerimental results. To this end, it was shown that the UrQMD simulations agree well with measured articleratios. Further details can be taen from Ref. [50]. The resulting freezeout temerature and chemical otential are almost the same when extensive and GNS fits are alied. As discussed in earlier sections, it was widely believed that the nonextensivity in the article roduction should be accomanied by a radical change in intensive or extensive thermodynamic quantities, such as the internal energy and the temerature. Therefore, the resulting temerature, for instance, was taen as a socalled hysical one, which afterwards should be extraolated to the freezeout temerature. The shortcuts of these assumtions shall be discussed in a forthcoming wor. To remain within the scoe of the resent word, we reort on nonextensivestatistical fit of various articleratios, where the intensive and extensive thermodynamic roerties of the strongly interacting system remain almost unchanged, while the nonextensitivity is defined by the equivalence classes (c, d). Their values are found close to unity, Tab.. This refers to a generic nonextensivity, which is aarently not of Tsallistye. The latter is characterized by (q,0) with q >. It should be remared that the excellent fits of various articleratios to the GNS statistical aroach shouldn t be interreted due to adding extra arameters(c, d). This should be related to rightly imlementing generic statistical aroach. Should the statistical nature of roducing various articleratios were really Tsallis nonextensivity, the added extra arameters should be (q, 0). Should this were extensive, (, ) should be resulted in. On the other hand, the success of the GNS statistical aroach in describing the lattice thermodynamics at (, ) confirms the caability of the roosed GNS to determine whether the system of interest has an extensive or a nonextensive statistical nature and aarently allows its imlementation in other ambiguous rocesses. 3.2 Equivalence classes (c, d) Table summarizes the freezeout arameters (T ch and µ b ) and the scaling exonents (c, d) as determined from the statistical fits of HRG with the roosed generic nonextensivity, GNS, to various exerimental results. Both freezeout arameters are grahically illustrated in Fig. 4 and comared with other results, section 3.. 7

8 The resulting classes (c, d) differ from unity referring to LambertW exonentials characterizing the entroic equivalence classes. The generalized logarithm and exonential functions were given in Eq. (6) and (2), resectively. This means that the article roductions has a nonbg statistical nature. Over the whole energyrange, the entroic equivalence classes can be exressed as LambertW exonentials, which is generic. Both resulting classes are ositive but less than unity, excet at LHC energy. We observe a general trend that the exonent c is closer to unity than d. At s 20 GeV, d reaches minimum, while c nearly maximum values. The differences between the two exonents becomes great at lowsps energies. The maximum error in the exonent is only 0%. Taing into account a wide energy range, during which relevant microscoic hysics radically changes, this is excellently recise. Let us tae d as resulted, i.e. d > 0, and assume that the resulting c can be aroximated to unity. Accordingly, the resulting entroy is characterized by stretched exonential, { [ ( ε (,d,r) (x) = ex dr x ) ]} /d. (7) r Mathematically, this function can be interreted as a fractional ower law. A critical exonent is inserted into an ordinary exonential function. Accordingly, in a disordered system, for instance, the article roduction, this tye of entroy characterizes a delayed relaxation. 4 Conclusions and final remars In the resent wor, we have introduced a systematic study for (non)extensive roerties characterizing relativistic heavyion collisions, which are assumed to hadronize and then i.e. forming hadrons or articles. With the nowadays detector technologies, the latter in turn can very recisely be detected. In doing this, we imly GNS on the wellnown statisticalthermal models, such as the HRG model. It is assumed that the degree of (non)extensivity can be determined by the hase sace characterizing the (dis)ordered system of interest. In other words, the roosed aroach is able to determine whether the system of interest has extensive or nonextensive roerties. Furthermore, more details about the nonextensivity can be also determined. The article ratios at energies ranging between 3.8 and 2760 GeV are best reroduced by GNS, where the equivalence classes (c,d) range between 0.9 and. This leads to a crucial conclusion that the statistics describing the article roduction at a wide range of beam energies remarably differ from extensive BG (, ) or nonextensive Tsallis (q, 0). It is certainly interesting to see how the economic aroach (extensive BG), as it deals with smallest free arameters, wors quite well. In fact, this is the main conclusion of this resent wor, either generic statistical aroach with (c,d) = (, ) and the extensive BG aroach reroduces well various article ratios in a wide range of energies. Therefore, when comaring our results with the qentroies, we can roose to refute all revious relevant wors in the literature. Such a simlification will strengthen the discussion. The advantage of the resent wor is now simly to reassure that the article ratios have extensive statistical nature but to illustrate that the roosed generic aroach indeed manifests the degree of (non)extensivity deending on the resulting equivalence classes (c, d). In a revious wor, we have examined this with the lattice QCD thermodynamics, which is er definition assumes additivity (extensivity). Now, we resent another examination with the article ratios roduced in heavyion collisions at various beam energies, which were ambiguously analyzed by thermal models assuming an ac hoc extensive statistics. Furthermore, the resulting equivalence classes imly that the system of interest can best be described by an extended exonential entroy, which basically differs from the wellnow extensive BG and from the wellnow nonextensive Tsallis entroy. Both are very secial cases defined by secific values of the equivalence classes (, ) and (q, 0), resectively. We conclude that the article roduction, in terms of the roduced articleratios, in relativistic heavyion collisions seems to be originated from a dynamically disordered system. Our results roose a lausible interretation why the resulting Tsallis freezeouttemerature, which is called hysical temerature, differs from the resulting BG freezeouttemerature, which is called the temerature. We believe that the insist to imly Tsallistye nonextensivity to the highenergy article roductions seems to be accomanied with a high rice; a radical change in intensive and/or extensive thermodynamic 8

9 quantities such as temerature and internal energy. Our results roose that the resulting freezeout arameters are comatible with the ones obtained when BG statistics is imlied. Accordingly, the wellnown freezeout hasediagram excellently agrees with the recent lattice redictions and with the freezeout hasediagram based on extensive BG statistics. The latter is very obvious, esecially when comaring the equivalence classes obtained in our calculations and the ones characterizing BG. We conclude that the nonextensivity characterizing the highenergy article roductions is solely through two critical exonents defining equivalence classes (c, d). A few final remar on the excellent fits reorted in this aer is now in order. First, as discussed, that the GNS artition function, Eq. (), was imlemented, in which an additional air of free arameters, (c, d), was added, shouldn t be thought as an exlanation based on statistical recision, which liely increases with increasing the number of the free arameters. The other air of free arameters, T ch and µ b, agreeswell with the extensive fits, while the resulting (c, d) can be aroximated to unity, i.e. almost extensive. Thus, the GNS artition function, Eq. (), with its four free arameters is not just a statistical recision. It manifests a generic statistical aroach that seems to characterize various asects for the article roductions. Because of its generic nature, both secial cases, BG and Tsallis, are also included. Its equivalence classes(c, d) define the degree of (non)extensivity. It is obvious that the statistical hysics aroaches to the highenergy henomena have to be understood as an economic descrition of the henomena of interest in terms of a small number of variables including the thermodynamic variables. This is true as long as the statistical nature hasn t be changed. Tsallis statistics is not only introducing an additional arameter, q, but also an secial concet of nonextensivity. In Euclidean field theory or lattice field theory, change of exonential factor imlies change of quantum theory, which does not seem to be ossible as the extensivity is fundamentally imosed. The resent aer resents an interretation to the highenergy henomena, the roduction of various article ratios, through suggesting changes of the foundations of statistical hysics. References [] A Tawfi Int J Mod Phys A (204) [2] A Tawfi J Phys G (203) [3] R Hagedorn Nuovo Cimento Sul 3 47 (965) [4] A Tawfi Prog Theor Phys (20) [5] C Bec Physica A (2000) [6] I Bediaga E M F Curado J M de Miranda Physica A (2000) [7] C Tsallis J Stat Phys (988) [8] D Prato and C Tsallis Phys Rev E (999) [9] A Tawfi Eur Phys J A (206) [0] A Tawfi Proceedings, 38th International Conference on High Energy Physics (ICHEP 206): Chicago, IL, USA, August 30, (206) [] A Tawfi H Yassin E R Abo Elyazeed Chin Phys C (207) [2] J Cleymans and D Woru J Phys G (202) [3] M D Azmi and J Cleymans J Phys G (204) [4] S Abe Physica A (2006) [5] S Abe Physica A (200) [6] S Abe Phys Rev E (200) 9

10 [7] T S Biro P Van Phys Rev E (20) [8] R Hanel S Thurner Eurohys Lett (20) [9] R Hanel S Thurner Eurohys Lett (20) [20] A Deman Phys Rev D (206) [2] A Deman Physica A (202) [22] I Sena A Deman AIP Conf Proc (203) [23] T S Biro G G Barnafoldi P Van Eur Phys J A 49 0 (203) [24] E Megias D P Menezes A Deman Physica A 42 5 (205) [25] HT Ding F Karsch S Muherjee Int J Mod Phys E (205) [26] A Tawfi E Abbas Phys Part Nucl Lett 2 52 (205) [27] A Tawfi et al. Int J Mod Phys E (205) [28] F Karsch K Redlich A Tawfi Eur Phys J C (2003) [29] F Karsch K Redlich A Tawfi Phys Lett B (2003) [30] K Redlich F Karsch A Tawfi J Phys G 30 S27 (2004) [3] A Tawfi J Phys G 3 S05 (2005) [32] A Tawfi Phys Rev D (2005) [33] A Tawfi Fizia B 8 4 (2009) [34] A Tawfi D Toublan Phys Lett B (2005) [35] A Tawfi submitted to PEPAN Lett [36] J Letessier J Rafelsi Hadrons and QuarGluon Plasma (UK:Cambridge University Press) T Ericson P V Landshoff (2004) [37] L Ahle et al. (E866 Collaboration and E97 Collaboration) Phys Lett B (2000) [38] J L Klay et al. (E895 Collaboration) Phys Rev Lett (2002) [39] L Ahle et al. (E866 Collaboration and E97 Collaboration) Phys Lett B 476 (2000) [40] A Andronic P BraunMunzinger J Stachel Nucl Phys A (2006) [4] A Tawfi M Y ElBary D M Habashy M T Mohamed E Abbas Int J Mod Phys E (206) [42] A Tawfi H Yassin E R. Abo Elyazeed Phys Rev D (205) [43] A Tawfi Phys Rev C (203) [44] A Tawfi Nucl Phys A (204) [45] A Tawfi Adv High Energy Phys (203) [46] A Tawfi Nucl Phys A (2006) [47] A Tawfi Eurohys Lett (2006) [48] A Tawfi PEPAN Lett 5 (208) [49] P BraunMunzinger J Stachel J P Wessels N Xu Phys Lett B 365 (996) [50] A Tawfi L I AbouSalem A G Shalaby M Hanafy A Sorin O Rogachevsy W Scheinast Eur Phys J A (206) 0

11 snn (GeV) T ch (MeV) µ b (MeV) c d χ 2 /dof Tab. : The freezeout arameters (T ch and µ b ) and the scaling exonents (c, d) as determined from the roosed GNS fits to various exerimental results. 0 (a) s NN = 2760 GeV c =.009 d =.03 0 (b) s NN = 200 GeV c = d = Fig. : Left anel (a): different article ratios measured at 2760 GeV (symbols) are fitted to HRG with GNS (dashed lines). Right anel (b) shows the same but at 200 GeV. The resulting exonents (c, d) and χ 2 are given in to right corners and Tab.. The freezeout temerature and chemical otential are also listed in Tab..

12 0 (a) s NN = 30 GeV c = d = (b) s NN = 62.4 GeV c = d = (c) s NN = 39 GeV c = 0.99 d = (d) s NN = 27 GeV c = d = (e) s NN = 9.6 GeV c = d = (f) s NN =.5 GeV c = d = (g) s NN = 7.7 GeV c = d = Fig. 2: Panel (a): different articleratios deduced from GNS fits (dashed lines) to the exerimental results at 30 GeV (symbols). Panels (b), (c), (d), (e), (f), and (g) show the same but at 62.4, 39, 27, 9.6,.5, and 7.7 GeV, resectively. The scaling exonents (c, d) and χ 2 are given in to right corners and Tab.. 2

13 0 (a) s NN = 7 GeV c = 0.99 d = (b) s NN = 2 GeV c = 0.99 d = (c) s NN = 9 GeV c = 0.99 d = (d) s NN = 8 GeV c = d = (e) s NN = 6.4 GeV c = 0.97 d = (f) s NN = 5 GeV c = d = (g) s NN = 4.2 GeV c = d = (h) s NN = 3.8 GeV c = d = Fig. 3: The same as in Fig. 2 but at 7 (a), 2 (b), 9 (c), 8 (d), 6.4 GeV (e), 5 (f), 4.2 (g), and 3.8 (h). 3

14 T ch [MeV] s/t 3 = const. Andronic et al. Tawfi et al. Extensive fits for UrQMD This wor µ b [GeV] Fig. 4: The freezeout arameters, temerature T ch vs. baryon chemical otential µ b which haven been deduced from GNS fits of HRG (closed circle) to various articleratios, Figs, 2, and 3. The symbols refer to freezeout arameters deduced from extensive fits: Andronic et al. [40], Tawfi et al. [26, 27], and UrQMD [50]. 4