Chinese Physics C Vol. xx,No. x (201x) xxxxxx On thermodynamic self-consistency of generic axiomatic-nonextensive statistics Abdel Nasser Tawfik1,2, Hayam Yassin3, Eman R. Abo Elyazeed3 1EgyptianCenterforTheoretical Physics(ECTP),ModernUniversityforTechnology andInformation(MTI),11571Cairo,Egypt 2WorldLaboratoryforCosmologyAndParticlePhysics(WLCAPP),11571Cairo,Egypt 7 3PhysicsDepartment, FacultyofWomenforArts,ScienceandEducation,AinShamsUniversity,11577Cairo,Egypt 1 0 Abstract: Genericaxiomatic-nonextensivestatistics characterizedbytwoasymptoticproperties,toeachofthema 2 scalingfunctionisassigned,characterizedbytheexponentscanddforfirstandsecondscalingproperty,respectively, n is formulated in a grand-canonical ensemble with finite volume in the thermodynamic limit. The thermodynamic a propertiesofarelativisticidealgasofhadronresonancesarestudied,analytically. Itisfoundthatthisgenericstatistics J satisfiestherequirementsoftheequilibriumthermodynamics. Essentialaspectsofthethermodynamicself-consistency 4 1 areclarified. Analyticalexpressionsareproposedforthestatisticalfitsofvarioustransversemomentumdistributions measuredinmost-centralcollisions atdifferentcollision energiesandcollidingsystems. Estimationsforthefreezeout ] temperature (T ) and the baryon chemical potential (µ ) and the exponents c and d are determined. The earlier h ch b are found compatible with the parameters deduced from Boltzmann-Gibbs (BG) statistics (extensive), while the t - latter refer to generic nonextensivities. The resulting equivalence class (c,d) is associated to stretched exponentials, l c where Lambert function reaches its asymptotic stability. In some measurements, the resulting nonextensiveentropy u is linearly composed on extensive entropies. Apart from power-scaling, the particle ratios and yields are excellent n [ quantities to highlighting whether the particle production takes place, (non)extensively. Various particle ratios and yields measured by the STAR experiment in central collisions at 200, 62.4 and 7.7 GeV are fitted with this novel 1 approach. We found that both c and d<1, i.e. referring to neither BG- nor Tsallis-type statistics, but to (c,d)- v 7 entropy,whereLambertfunctionsexponentially raise. Thefreezeout temperatureandbaryonchemicalpotentialare 9 found comparable with the ones deduced from BG statistics (extensive). We conclude that the particle production 6 at STAR energies is likely a nonextensiveprocess but not necessarily BG or Tsallis type. 4 0 Key words: Nonextensivethermodynamical consistency, Boltzmann and Fermi-Dirac statistics . 1 PACS: 05.70.Ln, 05.70.Fh,05.70.Ce 0 7 1 1 Introduction : v i X Intheoryofextensive(suchasBoltzmann-Gibbs)andnonextensive(suchasTsallis)statistics,thethermodynamic consistency gives a phenomenological description for various phenomena in high-energy experiments [1, 2]. The r a earlier was utilized by Hagedorn [3] in 1960’s to prove that the fireballs or heavy resonances lead to a bootstrap approach, i.e. further fireballs, which - in turn - consist of smaller fireballs and so on. The implementation of the nonextensiveTsallisstatisticswasintroducedinRefs. [5–7]. Assumingthatthedistributionfunctiongetsvariations, due to possible symmetrical change, Tawfik applied nonextensive concepts to high-energy particle production [4]. Recently, variouspapers were quite successful in explaining various aspects of the high-energy particle production using thermodynamically consistent nonextensive statistics of Tsallis-type [8–14]. These are basedon the conjecture that replacing the Boltzmann factor by the q-exponential function of the Tsallis statistics, with q > 1, leads to a good agreementwith the experimental measurements at high energies. Recently, Tawfik explained that this method seems to fail to assure a full incorporation of nonextensivity because fluctuations, correlations, interactions among the produced particles besides the possible modification in the phase space of such an interacting system are not properlytakenintoaccount[15]. Again,theTsallisdistributionwaswidelyappliedtodescribethehadronproduction [8–14]. At high transverse momentum spectra (p ), some authors did not obtained power-law,while at low p , they T T obtained an approximate exponential-distribution. Received30August2016 1)E-mail:a.tawfi[email protected] (cid:13)c2013 Chinese Physical Society and the Institute of High Energy Physics of the Chinese Academy of Sciences and the Institute of ModernPhysicsoftheChineseAcademyofSciences andIOPPublishingLtd 010201-1 Chinese Physics C Vol. xx,No. x (201x) xxxxxx In light of such a wide discrepancy and especially to find a unified statistical description for high-energy col- lisions, we introduce generic axiomatic-nonextensive statistics, in which the phase space determines the degree of (non)extensivity. The latter is not necessarily limited to extensive and/or intensive thermodynamic quantities, such as, temperature and baryon chemical potential. Regarding the success of Tsallis-statistics in describing the transverse momentum spectra, especially at high p , T we first recall that Bialas claimed that such good fit would be incomplete [16]. First, it is believed to ignore the contradiction between the applicability of such a statistical thermal approach at high energy and the perturbative QCD. The crucial question, how comes that the statistical thermal approach which describes well the lattice ther- modynamics (non-perturbative)[1] should be assumed to do such claimed excellent job at ultra-relativistic energies, should be answered. Second, such a statistical fit seems to disbelieve the role of statistical cluster decay, which can be scaled as power laws very similar to that from Tsallis statistics. This simply means that the decay of statistical clustersisconjecturedtobecapabletoexplaintheexcellentreproductionofthemeasuredtransversemomentarather than the Tsallis-type nonextensivity [15]. Also, Bialas [16] presented within the statistical cluster-decay model a numerical analysis for the hadronization processes. It was found that the resulting transverse-momentum distribution can be a Tsallis-like one. Only in a very special case, where the fluctuations of the Lorentz factor and the temperature are given by Beta and Gamma distributions, respectively, the well-known Tsallis distribution can be obtained. The origin of these fluctuations was introduced in Ref. [17]. Bialas explained that the produced hadronic clusters-decay, purely thermally, i.e. following Boltzmann-‘Gibbs (BG) statistics [16]. Furthermore, it is also supposed that the produced hadronic clusters move with a fluctuating Lorentz factor in the transverse direction, i.e. a power law. Thus, the production and the decay of such clusters would be be regarded as examples on superstatistics. The latter can be understood as a kind of a superposition of two different types of statistics corresponding to nonequilibrium systems [18, 19]. In other words, it was concluded [17] that even superstatistics should be based on more general distributions than the Gamma type before being applied to multiparticle productions in high energy collisions. For the sake of completeness, we recall that the various power-law distributions have been already implemented to pp-collisions [20–24]. The crucialquestionremains unanswered,namely,what is the originandthe degreeof nonextensivity? And how to determine the degree of nonextensivity in a strongly correlated system as the relativists heavy-ion collisions? In the present work, we analyse the thermodynamic self-consistency of the generic axiomatic-nonextensive approach, which was introduced in [25, 26] and formulated for further implications in high-energy physics in [15, 27]. In thermal equilibrium, the statistical mechanics is a thermodynamically self-consistent theory. It fulfils the requirements of the equilibrium thermodynamics, because thermodynamic potential - in its thermodynamic limit - can be expressed as a first-order homogeneous function of extensive variables [28]. The entropy belongs to the fundamental checksfor the thermodynamic self-consistency. For a nonextensive quantumgas,the entropy shouldbe fully constructedfrombosonsandfermionscontributions; S =SFD+SBE, whereSFD andSBE areFermi-Diracand q q q q q Bose-Einstein nonextensive entropy, respectively [29]. In microcanonical [30], canonical [31], and grand-canonical [32] ensembles ensembles, the thermodynamic self-consistence of Tsallis statistics has been proved. If its entropic variable has extensive property, it was demonstrated that the homogeneity of the thermodynamic potential leads to the zeroth law of thermodynamics, i.e. additivity principle and Euler theorem. Also, the first and second laws of thermodynamics should be fulfilled. An additional ingredient of special importance in the particle production, namely the fireball self-consistency principle, should be guaranteed,as well. The present paper is organized as follows. The generic axiomatic-nonextensive statistics is reviewed in section 2. The thermodynamic self-consistency will be discussed in section 3. This is divided into nonextensive Boltzmann- Gibbsstatistics(section3.1)andgenericaxiomatic-nonextensivequantumstatistics(section3.2). Sections4and4.2 are devoted to the fitting of the transverse momentum distributions and particle ratios and different beam energies and in different system sizes, respectively. Section 5 elaborates the final conclusions. 2 Short reminder to generic axiomatic-nonextensive statistics Based on Hanel-Thurner entropy [25–27] which is fully expressed in Ref. [15], Ω S [p] = Γ(d+1,1 clogp ) p , (1) c,d A − i −B i Xi=1 where p is the probability of i-th state and Γ(a,b) = ∞dtta−1exp( t) is incomplete gamma-function with i b − A and being arbitrary parameters, the generic axiomatic-Rnonextensive partition function for statistical processes in B 010201-2 Chinese Physics C Vol. xx,No. x (201x) xxxxxx high-energy physics was suggested [27]. For a classical gas, NM|B ∞ d3p lnZ (T,µ) = V g ε (x ), (2) cl iZ (2π)3 c,d,r i Xi 0 where V is the fireball volume and x =β[µ E (p)] with E (p)=(p2+m2)1/2 being the dispersion relation of i-th i i− i i i state (particle). Straightforwardly,the quantum gas partition function reads NM|B ∞ d3p lnZ (T) = V g ln[1 ε (x )], (3) FB ± iZ (2π)3 ± c,d,r i Xi 0 where represent fermions (subscript F) and bosons (subscript B), respectively. The distribution function ε (x ) ± c,d,r i is given as [25, 26] d x d1 ε (x)=exp − W B 1 W [B] , (4) c,d,r (cid:20)1 c(cid:18) k(cid:20) (cid:16) −r(cid:17) (cid:21)− k (cid:19)(cid:21) − where W is Lambert W-function which has real solutions at k=0 with d 0 and at k=1 with d<0, k ≥ (1 c)r (1 c)r B= − exp − , (5) 1 (1 c)r (cid:20)1 (1 c)r(cid:21) − − − − with r=[1 c+cd]−1 and c,d are two constants to be elaborated, shortly. Eq. (4) is valid for both classical and − quantum gas. Theequivalentclasses(c,d)standfortwoexponentsgivingestimationsfortwoscalingfunctionswithtwoasymp- totic properties [25, 26]. Statistical systems in their large size limit violating the fourth Shannon-Khinchin axiom are characterized by a unique pair of scaling exponents (c,d). Such systems have two asymptotic properties of their associatedgeneralizedentropies. Bothpropertiesareassociatedwithonescalingfunctioneach. Eachscalingfunction is characterizedby one exponent; c for first and d for second property. These exponents define equivalence relations of entropic forms, i.e. two entropic forms are equivalent if their exponents are the same. The various thermodynamic observables such as pressure (p), number (n=N/V), energy (ǫ=E/V) and entropy density (s=S/V), respectively [33] can be derived from Eq. (9) or Eq. (10) ∂lnZ ∂p T2∂lnZ µT ∂lnZ ∂p p=T , n= , ǫ= + , s= . (6) ∂V ∂µ V ∂T V ∂µ ∂T 3 Thermodynamic self-consistency Here, we examine the thermodynamic self-consistencyof genericaxiomatic-nonextensivestatistics inBoltzmann- Gibbs and quantum gases. The procedure goes as follows. We start with first and second laws of thermodynamic which controlthe system of interest. Then, we determine the thermodynamic properties of that system and confirm that both laws of thermodynamic are verified. The first law of thermodynamics describes the change in energy (dE) in terms of a change in volume (dV) and entropy (dS): dE(V,S)= PdV +TdS, (7) − where T and P are temperature and pressure coefficients, respectively. In the variationof the free energy, F(V,T)= E T S, the dependence on the entropy as given in Eq. (7) is to be replaced by a temperature-dependence − dF(V,T)=dE TdS SdT PdV SdT. (8) − − ≡− − If the system of interest contained a conserved number (N), it turns to be necessary to introduce chemical potential (µ). Eqs. (7) and (8) should be extended dE(V,S,N) = PdV +TdS+µdN, (9) − dF(V,T,µ) = PdV SdT Ndµ. (10) − − − 010201-3 Chinese Physics C Vol. xx,No. x (201x) xxxxxx The pressure is derived as ∂F p= . (11) − ∂V (cid:12) (cid:12)T,µ (cid:12) From Eqs. (9), (10) and (11), thermodynamics relations can(cid:12)be deduced [33]. Their justification proves the thermo- dynamic self-consistency, ∂p ∂ǫ ∂p ∂ǫ n= , T= , s= , µ= . (12) ∂µ(cid:12) ∂s(cid:12) ∂T(cid:12) ∂n(cid:12) (cid:12)T (cid:12)n (cid:12)µ (cid:12)s (cid:12) (cid:12) (cid:12) (cid:12) (cid:12) (cid:12) (cid:12) (cid:12) To verify the second law of thermodynamics one has to prove that ∂s 0. ≥ In the section that follows, all these thermodynamic quantities shall be derived from the generic axiomatic- nonextensive statistics for Boltzmann-Gibbs, Bose-Einstein and Fermi-Dirac statistical ensembles [32]. 3.1 Nonextensive Boltzmann-Gibbs statistics As discussed in Ref. [34], the universality class (c,d) is conjectured not only to characterize the entropy of the systemof interest, entirely, but also to specify the distribution functions ofthat systemin the thermodynamic limit. Thus,itislikelyabletodetermine(non)extensivityofthesystem. Forinstance,if(c,d)=(1,1),thesystemcanbewell describedbyBG statistics,while if (c,d)=(q,0), the systempossessesa Tsallis-typenonextensovity. Furthermore,if (c,d)=(1,d), stretched exponentials characterize that system. To our knowledge, further details about the physical meaning of (c,d) are being worked out by many colleagues. This would be published in the near future. The last case,for instance, requiresthat d>0 and c 1so that lim ε (x)=exp( dr[1 x/r]1/d 1). In lightof this, for → c→1 c,d,r − − − a single particle, the Boltzmann distribution in nonextensive system, Eq. (9) can be expressed as 1 f(x)= , (13) ε (x) c,d,r from which various thermodynamic quantities can deduced gT ∞ p = p2ln[ε (x)]dp, (14) 2π2Z c,d,r 0 n = g ∞ p2W0hB(1−xr)d1i dp, (15) 2π2Z0 (1−c)[r−x](cid:16)1+W0hB(1−xr)d1i(cid:17) ǫ = −g ∞ p2(xT)W0hB(1−xr)d1i dp+ gµ ∞ p2W0hB(1−xr)d1i dp, (16) 2π2Z0 (1−c)[r−x](cid:16)1+W0hB(1−xr)d1i(cid:17) 2π2Z0 (1−c)[r−x](cid:16)1+W0hB(1−xr)d1i(cid:17) s = g ∞p2ln[ε (x)]dp g ∞ p2(xT)W0hB(1−xr)1di dp. (17) 2π2Z0 c,d,r −2π2T Z0 (1−c)[r−x](cid:16)1+W0hB(1−xr)d1i(cid:17) In an ideal gas approach, such as the hadron resonance gas (HRG) model, the thermodynamic quantities get con- tributions from each hadron resonance. Thus, a summation over different hadron resonancesof which the statistical ensemble is consisting, should be added in front of the right-hand side. To check and satisfy the thermodynamic self-consistency, let us rewrite Eq. (17) as p ǫ µn s= + , (18) T T − T and take its derivative with respect to ǫ, then we get ∂s 1 = . (19) ∂ǫ(cid:12) T (cid:12)n (cid:12) (cid:12) Furthermore, we derive from Eq. (9) [33], p= ǫ+Ts+µn. (20) − 010201-4 Chinese Physics C Vol. xx,No. x (201x) xxxxxx At constant T, the derivative of pressure with respect to µ reads ∂p ∂ǫ ∂s ∂n = +T +n+µ . (21) ∂µ(cid:12) −∂µ ∂µ ∂µ (cid:12)T (cid:12) The differentiations of Eqs. (15), (16), and (17(cid:12)) with respect to µ lead to 2 ∂n = g ∞ p2W0hB(1−xr)1di dp ∂µ 2π2T Z0 (d−dc)r2(1−x/r)2(cid:16)1+W0hB(1−xr)d1i(cid:17)3 g ∞ p2W0hB(1−xr)d1i dp − 2π2T Z0 (d−dc)r2(1−x/r)2(cid:16)1+W0hB(1−xr)d1i(cid:17)2 + g ∞ p2W0hB(1−xr)d1i dp, (22) 2π2T Z0 (1−c)r2(1−x/r)2(cid:16)1+W0hB(1−xr)d1i(cid:17) 2 ∂ǫ = gµ ∞ p2 (cid:16)W0hB(1−xr)d1i(cid:17) dp ∂µ 2π2T Z0 (d−dc)r2(1−x/r)2(cid:16)1+W0hB(1−xr)d1i(cid:17)3 gµ ∞ p2W0hB(1−xr)d1i dp − 2π2T Z0 (d−dc)r2(1−x/r)2(cid:16)1+W0hB(1−xr)d1i(cid:17)2 + gµ ∞ p2W0hB(1−xr)d1i dp 2π2T Z0 (1−c)r2(1−x/r)2(cid:16)1+W0hB(1−xr)d1i(cid:17) 2 g ∞ p2(xT)W0hB(1−xr)1di dp − 2π2T Z0 (d−dc)r2(1−x/r)2(cid:16)1+W0hB(1−xr)d1i(cid:17)3 + g ∞ p2(xT)W0hB(1−xr)d1i dp 2π2T Z0 (d−dc)r2(1−x/r)2(cid:16)1+W0hB(1−xr)d1i(cid:17)2 g ∞ p2(xT)W0hB(1−xr)d1i dp, (23) − 2π2T Z0 (1−c)r2(1−x/r)2(cid:16)1+W0hB(1−xr)d1i(cid:17) 2 ∂s = −g ∞ p2(xT)W0hB(1−xr)d1i dp ∂µ 2π2T2Z0 (d−dc)r2(1−x/r)2(cid:16)1+W0hB(1−xr)d1i(cid:17)3 + g ∞ p2(xT)W0hB(1−xr)d1i dp 2π2T2Z0 (d−dc)r2(1−x/r)2(cid:16)1+W0hB(1−xr)d1i(cid:17)2 g ∞ p2(xT)W0hB(1−xr)d1i dp. (24) − 2π2T2Z0 (1−c)r2(1−x/r)2(cid:16)1+W0hB(1−xr)d1i(cid:17) Then, the substitution from Eqs. (15), (22), (23) and (24) into Eq. (21) gives ∂p =n. (25) ∂µ(cid:12) (cid:12)T (cid:12) (cid:12) 010201-5 Chinese Physics C Vol. xx,No. x (201x) xxxxxx Also, from Eq. (18) we can express the energy density as ǫ=Ts p+µn, (26) − Its differentiation with respect to n is ∂ǫ ∂p ∂µ ∂p∂µ ∂µ = +µ+n = +µ+n . (27) ∂n(cid:12) −∂n ∂n −∂µ∂n ∂n (cid:12)s (cid:12) (cid:12) By substituting from Eq. (25) into this previous equation, we get ∂ǫ =µ, (28) ∂n(cid:12) (cid:12)s (cid:12) (cid:12) which is a proof of the thermodynamic self-consistency. To prove the fourth equation in Eq. (12), let us differentiate Eq. (20) with respect to temperature at constant µ ∂p ∂s ∂ǫ ∂n =s+T +µ , (29) ∂T(cid:12) ∂T −∂T ∂T (cid:12)µ (cid:12) (cid:12) The differentiations of Eqs. (15), (16) and (17) with respect to T read 2 ∂n = −g ∞ p2(xT)W0hB(1−xr)d1i dp ∂T 2π2T2Z0 (d−dc)r2(1−x/r)2(cid:16)1+W0hB(1−xr)d1i(cid:17)3 + g ∞ p2(xT)W0hB(1−xr)d1i dp 2π2T2Z0 (d−dc)r2(1−x/r)2(cid:16)1+W0hB(1−xr)d1i(cid:17)2 g ∞ p2(xT)W0hB(1−xr)d1i dp, (30) − 2π2T2Z0 (1−c)r2(1−x/r)2(cid:16)1+W0hB(1−xr)d1i(cid:17) 010201-6 Chinese Physics C Vol. xx,No. x (201x) xxxxxx 2 ∂ǫ = g ∞ p2(xT)2W0hB(1−xr)d1i dp ∂T 2π2T2Z0 (d−dc)r2(1−x/r)2(cid:16)1+W0hB(1−xr)d1i(cid:17)3 g ∞ p2(xT)2W0hB(1−xr)d1i dp − 2π2T2Z0 (d−dc)r2(1−x/r)2(cid:16)1+W0hB(1−xr)d1i(cid:17)2 + g ∞ p2(xT)2W0hB(1−xr)d1i dp 2π2T2Z0 (1−c)r2(1−x/r)2(cid:16)1+W0hB(1−xr)d1i(cid:17) 2 gµ ∞ p2(xT)(cid:16)W0hB(1−xr)d1i(cid:17) dp − 2π2T2Z0 (d−dc)r2(1−x/r)2(cid:16)1+W0hB(1−xr)d1i(cid:17)3 + gµ ∞ p2(xT)W0hB(1−xr)d1i dp 2π2T2Z0 (d−dc)r2(1−x/r)2(cid:16)1+W0hB(1−xr)d1i(cid:17)2 gµ ∞ p2(xT)W0hB(1−xr)d1i dp, (31) − 2π2T2Z0 (1−c)r2(1−x/r)2(cid:16)1+W0hB(1−xr)d1i(cid:17) 2 ∂s = g ∞ p2(xT)2W0hB(1−xr)d1i dp ∂T 2π2T3Z0 (d−dc)r2(1−x/r)2(cid:16)1+W0hB(1−xr)d1i(cid:17)3 g ∞ p2(xT)2W0hB(1−xr)d1i dp − 2π2T3Z0 (d−dc)r2(1−x/r)2(cid:16)1+W0hB(1−xr)d1i(cid:17)2 + g ∞ p2(xT)2W0hB(1−xr)d1i dp. (32) 2π2T3Z0 (1−c)r2(1−x/r)2(cid:16)1+W0hB(1−xr)d1i(cid:17) Then, the substitution from Eqs. (17), (30), (31) and (32) into Eq. (29) leads to ∂p =s. (33) ∂T(cid:12) (cid:12)µ (cid:12) It is apparent that the given definitions of temperature, n(cid:12)umber density, chemical potential and entropy density are thermodynamically consistent. The first law of thermodynamics describes the consequences of heat transfer, Eq. (18), while the second law sets constrains on it, i.e. δs 0, ≥ 1 ∂s [∂ǫ µ∂n], (34) ≃T − where ∂µ/∂T is conjectured to vanish. Therefore, the second law of thermodynamics is fulfilled, if ∂ǫ µ∂n. (35) ≥ This inequality is fulfilled, when comparing Eq. (30) with Eq. (31). Two cases can be classified. Firstly, ∂ǫ>µ∂n is obvious, at arbitrary µ. Secondly, ∂ǫ=µ∂n is also obtained, at c=1 or µ=ǫ or µ=ǫ+Tr. 3.2 Generic axiomatic-nonextensive quantum statistics For quantum statistics, the distribution function reads 1 f(x)= , (36) 1 ε (x) ± c,d,r 010201-7 Chinese Physics C Vol. xx,No. x (201x) xxxxxx Accordingly, the various thermodynamic quantities can be deduced gT ∞ p = p2ln[1 ε (x)]dp, (37) 2π2Z ± c,d,r 0 n = ±g ∞ p2εc,d,r(x)W0hB(1−xr)1di dp, (38) 2π2Z0 (1−c)[1±εc,d,r(x)](r−x)(cid:16)1+W0hB(1−xr)d1i(cid:17) ǫ = ∓g ∞ p2εc,d,r(x)(xT)W0hB(1−xr)d1i dp 2π2Z0 (1−c)[1±εc,d,r(x)](r−x)(cid:16)1+W0hB(1−xr)d1i(cid:17) gµ ∞ p2εc,d,r(x)W0hB(1−xr)d1i dp, (39) ± 2π2Z0 (1−c)[1±εc,d,r(x)](r−x)(cid:16)1+W0hB(1−xr)d1i(cid:17) s = g ∞p2ln[1 ε (x)]dp g ∞ p2εc,d,r(x)(xT)W0hB(1−xr)d1i dp. (40) 2π2Z0 ± c,d,r ±2π2T Z0 (1−c)[1±εc,d,r(x)](r−x)(cid:16)1+W0hB(1−xr)d1i(cid:17) The differentiations of Eqs. (38), (39) and (40) with respect to µ give 2 ∂n = ±g ∞ p2εc,d,r(x)W0hB(1−xr)d1i dp ∂µ 2π2T Z0 (d−dc)[1±εc,d,r(x)]r2(1−x/r)2(cid:16)1+W0hB(1−xr)d1i(cid:17)3 g ∞ p2εc,d,r(x)W0hB(1−xr)d1i dp ∓ 2π2T Z0 (d−dc)[1±εc,d,r(x)]r2(1−x/r)2(cid:16)1+W0hB(1−xr)d1i(cid:17)2 2 g ∞ p2ζc,d,r(x)W0hB(1−xr)d1i dp − 2π2T Z0 (1−c)2[1±εc,d,r(x)]2r2(1−x/r)2(cid:16)1+W0hB(1−xr)d1i(cid:17)2 2 g ∞ p2εc,d,r(x)W0hB(1−xr)d1i dp ± 2π2T Z0 (1−c)2[1±εc,d,r(x)]r2(1−x/r)2(cid:16)1+W0hB(1−xr)d1i(cid:17)2 g ∞ p2εc,d,r(x)W0hB(1−xr)1di dp, (41) ± 2π2T Z0 (1−c)[1±εc,d,r(x)]r2(1−x/r)2(cid:16)1+W0hB(1−xr)d1i(cid:17) where 2d x d1 ζ (x)=exp − W B 1 W [B] . (42) c,d,r (cid:20)1 c(cid:18) k(cid:20) (cid:16) −r(cid:17) (cid:21)− k (cid:19)(cid:21) − 010201-8 Chinese Physics C Vol. xx,No. x (201x) xxxxxx 2 ∂ǫ = ±gµ ∞ p2εc,d,r(x)(cid:16)W0hB(1−xr)d1i(cid:17) dp ∂µ 2π2T Z0 (d−dc)[1±εc,d,r(x)]r2(1−x/r)2(cid:16)1+W0hB(1−xr)d1i(cid:17)3 gµ ∞ p2εc,d,r(x)W0hB(1−xr)d1i dp ∓ 2π2T Z0 (d−dc)[1±εc,d,r(x)]r2(1−x/r)2(cid:16)1+W0hB(1−xr)d1i(cid:17)2 2 gµ ∞ p2ζc,d,r(x)W0hB(1−xr)d1i dp − 2π2T Z0 (1−c)2[1±εc,d,r(x)]r2(1−x/r)2(cid:16)1+W0hB(1−xr)d1i(cid:17)2 2 gµ ∞ p2εc,d,r(x)W0hB(1−xr)d1i dp ± 2π2T Z0 (1−c)2[1±εc,d,r(x)]r2(1−x/r)2(cid:16)1+W0hB(1−xr)d1i(cid:17)2 gµ ∞ p2εc,d,r(x)W0hB(1−xr)d1i dp ± 2π2T Z0 (1−c)[1±εc,d,r(x)]r2(1−x/r)2(cid:16)1+W0hB(1−xr)d1i(cid:17) 2 g ∞ p2εc,d,r(x)(xT)W0hB(1−xr)d1i dp ∓ 2π2T Z0 (d−dc)[1±εc,d,r(x)]r2(1−x/r)2(cid:16)1+W0hB(1−xr)d1i(cid:17)3 g ∞ p2εc,d,r(x)(xT)W0hB(1−xr)d1i dp ± 2π2T Z0 (d−dc)[1±εc,d,r(x)]r2(1−x/r)2(cid:16)1+W0hB(1−xr)d1i(cid:17)2 2 + g ∞ p2ζc,d,r(x)(xT)W0hB(1−xr)1di dp 2π2T Z0 (1−c)2[1±εc,d,r(x)]r2(1−x/r)2(cid:16)1+W0hB(1−xr)d1i(cid:17)2 2 g ∞ p2εc,d,r(x)(xT)W0hB(1−xr)d1i dp ∓ 2π2T Z0 (1−c)2[1±εc,d,r(x)]r2(1−x/r)2(cid:16)1+W0hB(1−xr)d1i(cid:17)2 g ∞ p2εc,d,r(x)(xT)W0hB(1−xr)1di dp, (43) ∓ 2π2T Z0 (1−c)[1±εc,d,r(x)]r2(1−x/r)2(cid:16)1+W0hB(1−xr)d1i(cid:17) 010201-9 Chinese Physics C Vol. xx,No. x (201x) xxxxxx 2 ∂s = ∓g ∞ p2εc,d,r(x)(xT)W0hB(1−xr)d1i dp ∂µ 2π2T2Z0 (d−dc)[1±εc,d,r(x)]r2(1−x/r)2(cid:16)1+W0hB(1−xr)d1i(cid:17)3 g ∞ p2εc,d,r(x)(xT)W0hB(1−xr)d1i dp ± 2π2T2Z0 (d−dc)[1±εc,d,r(x)]r2(1−x/r)2(cid:16)1+W0hB(1−xr)d1i(cid:17)2 2 + g ∞ p2ζc,d,r(x)(xT)W0hB(1−xr)d1i dp 2π2T2Z0 (1−c)2[1±εc,d,r(x)]2r2(1−x/r)2(cid:16)1+W0hB(1−xr)d1i(cid:17)2 2 g ∞ p2εc,d,r(x)(xT)W0hB(1−xr)d1i dp ∓ 2π2T2Z0 (1−c)2[1±εc,d,r(x)]r2(1−x/r)2(cid:16)1+W0hB(1−xr)1di(cid:17)2 g ∞ p2εc,d,r(x)(xT)W0hB(1−xr)1di dp, (44) ∓ 2π2T2Z0 (1−c)[1±εc,d,r(x)]r2(1−x/r)2(cid:16)1+W0hB(1−xr)d1i(cid:17) From the substitution from Eqs. (41), (43), (44) and (55) into Eq. (21), we get ∂p =n. (45) ∂µ(cid:12) (cid:12)T (cid:12) (cid:12) Also, from Eqs. (39) and (40), we get a direct relation between entropy and energy density, p ǫ µn s= + . (46) T T − T Thus, at constant n and p the differentiation of entropy density with respect to energy density leads to ∂s 1 = , (47) ∂ǫ T which proves the thermodynamic consistency, Eqs. (12). Also, from Eq. (46) we can write another form of energy density ǫ=Ts p+µn. (48) − By differentiating the energy density with respect to n, we get ∂ǫ ∂p ∂µ ∂p∂µ ∂µ = +µ+n = +µ+n . (49) ∂n(cid:12) −∂n ∂n −∂µ∂n ∂n (cid:12)s (cid:12) (cid:12) Then, by substituting from Eq. (45) into the previous equation, ∂ǫ =µ, (50) ∂n(cid:12) (cid:12)s (cid:12) (cid:12) which is a proof of the thermodynamic self-consistency. 010201-10