arxiv: v3 [cond-mat.stat-mech] 19 Mar 2015

Transcription

1 Canonical Ensemble in Non-extensive Statistical Mechanics Julius Ruseckas Institute of Theoretical Physics and Astronomy, Vilnius University, A. Goštauto 2, LT-008 Vilnius, Lithuania arxiv: v3 [cond-mat.stat-mech] 9 Mar 205 The framework of non-extensive statistical mechanics, proposed by Tsallis, has been used to describe a variety of systems. The non-extensive statistical mechanics is usually introduced in a formal way, using the maximization of entropy. In this article we investigate the canonical ensemble in the non-extensive statistical mechanics using a more traditional way, by considering a small system interacting with a large reservoir via short-range forces. The reservoir is characterized by generalized entropy instead of the Boltzmann-Gibbs entropy. Assuming eual probabilities for all available microstates we derive the euations of the non-extensive statistical mechanics. Such a procedure can provide deeper insight into applicability of the non-extensive statistics. I. INTRODUCTION Complexity in natural or artificial systems may be caused by long-range interactions, long-range memory, nonergodicity or multifractality. Such systems have exotic thermodynamical properties and are unusual from the point of view of traditional Boltzmann-Gibbs statistical mechanics. Statistical description of complex systems can be provided using the non-extensive statistical mechanics that generalizes the Boltzmann-Gibbs statistics [ 3]. The non-extensive statistical mechanics has been used to describe phenomena in various in high-energy physics [4], spin-glasses [5], cold atoms in optical lattices [6], trapped ions[7], anomalous diffusion [8, 9], dusty plasmas[0], low-dimensional dissipative and conservative maps in the dynamical systems [ 3], turbulent flows [4], Langevin dynamics with fluctuating temperature [5, 6]. Concepts related to the non-extensive statistical mechanics have found applications not only in physics but in chemistry, biology, mathematics, economics, and informatics as well [7 9]. The basis of he non-extensive statistical mechanics is the generalized entropy [] S = µ p(µ), () where p(µ) is the probability of finding the system in the state characterized by the parameters µ; the parameter describes the non-extensiveness of the system. More generalized entropies and distribution functions are introduced in Refs. [20, 2]. The generalized entropy () can be written in a form similar to the Bolzmann-Gibbs entropy S BG = p(µ) ln p(µ) (2) as an average of -logarithm []: where the -logarithm is defined as S = µ µ p(µ)ln p(µ), (3) ln x = x. (4) In the limit the -logarithm becomes an ordinary logarithm, thus the Boltzmann-Gibbs entropy can be obtained from E. () in the limit [, 2]. The inverse function of the -logarithm is the -exponential function exp (x) [+( )x] +, (5) with [x] + = x if x > 0, and [x] + = 0 otherwise. The -exponential and -logarithm appear in many euations of non-extensive statistical mechanics[]. Some properties of -exponential and -logarithm are presented in Appendix B. Electronic address: julius.ruseckas@tfai.vu.lt; URL:

2 The euilibrium of an isolated system consisting of N particles is described by the microcanonical ensemble. In the statistical physics it is assumed that the euilibrium in the microcanonical ensemble corresponds to eually probable microstates [22, 23], therefore in the microcanonical ensemble p(µ) = /W, where W is the number of microstates. Non-extensive statistical mechanics can describe non-ergodic systems where not all available microstates can be reached. In this case W is the effective number of microstates, that is the number of microstates whose probability is not zero. When probabilities are eual, E. () for the generalized entropy takes the simpler form S = ln W. (6) In the systems with long-range interactions and long-range correlations the effective number of microstates W can grow not exponentially with the number of particles in the system N but slower, as a power-law of N. For such a systems the standard Boltzmann-Gibbs entropy (2) is not proportional to the number of particles in the system and thus is not extensive. The extensive uantity is the generalized entropy () for some value of. In general, if the entropy S is proportional to the number of particles N then the number of microstates W grows as W exp N. There are two different cases: (i) < and W N /( ). The number of microstates grows as a power-law. (ii) > and W behaves as ( ( )AN) /( ). In this case there is a maximum value of the number of particles N crit where the number of microstates becomes infinite and thus the macroscopic limit N cannot be taken. Because the of this complication occurring when > in this paper we consider only the case of < ; the value of in all the euations below should be assumed to be less than. The case of > warrants a separate investigation and is outside of the scope of the present paper. The canonical ensemble in the non-extensive statistical mechanics is usually introduced in a formal way, starting from the maximization of the generalized entropy[]. The physical assumptions appear in the maximization procedure in the form of constraints. However, the -averages used for constraints are unusual from the point of view of ordinary, Boltzmann-Gibbs statistics. The physical justification of -averages and escort distributions is not completely clear. Thus a more physically transparent method would be useful for understanding the non-extensive statistical mechanics. The goal of this paper is to investigate the canonical ensemble in the non-extensive statistical mechanics using a more traditional way, by considering a small system interacting with a large reservoir via short-range forces. Consistent investigation of such a situation has not been performed yet. We assume that the generalized entropy () for some value of < instead of the Boltzmann-Gibbs entropy is extensive for the reservoir. In addition, as in the standard statistical mechanics we assume eual probabilities for all available microstates of the combined system consisting of the small system and the reservoir. By doing so we can avoid the critiue of the generalized entropy presented in Refs. [24, 25]. The paper is organized as follows: To make the comparison of the non-extensive statistical mechanics with the standard Boltzmann-Gibbs statistical mechanics easier, in Section II we briefly present the usual construction of the canonical ensemble in the standard statistical mechanics. In Section III we consider the canonical ensemble in the non-extensive statistical mechanics and in Section IV we explore the resulting Legendre transformation structure. Section V summarizes our findings. 2 II. CANONICAL ENSEMBLE IN BOLTZMANN GIBBS STATISTICAL MECHANICS To highlight the differences from the non-extensive statistical mechanics, let us at first briefly review the canonical ensemble in the extensive Boltzmann-Gibbs statistical mechanics. The standard approach [22, 23] is to consider a composite system consisting of a system under investigation S interacting with a large reservoir R. The system S has energy E, the energy of the reservoir R is E R and the energy of the composite system is E tot. Due to the interaction the system S and the reservoir R can exchange energy. The interaction is assumed to be short-range, therefore in the macroscopic limit the energy of the interaction is negligible and the total energy of the composite system is E tot = E + E R. For simplicity we assume that there is no exchange of the particles between the system S and the reservoir R. The number of microstates in the system S having the energy E is W(E) and the number of microstates in the reservoir is W R (E R ). Here it is assumed that the numbers of microstates depend only on the energy. The short range interactions of the system S with the reservoir R do not significantly change the numbers of microstates and thus the total number of microstates in the combined system when the system S has energy E is W(E)W R (E tot E). The full number of microstates W tot (E tot ) of the combined system is obtained summing over all available energies of the system S: W tot (E tot ) = E W(E)W R (E tot E). (7)

3 Introducing the entropy of the system S(E) = lnw(e) and the entropy of the reservoir S R (E R ) = lnw R (E R ) we can write 3 W tot (E tot ) = E e S(E)+ S R(E tot E). (8) In the extensive Boltzmann-Gibbs statistics the entropy of the reservoir S R is proportional to the number of particles N R in the reservoir and is macroscopically large. The sum of large exponentials can be approximated by the largest term, as is described in the Appendix A. In the statistical mechanics it is postulated that in the euilibrium the probability of each microstate is the same and eual /W tot. Thus the most probable state of the composite system corresponds to the largest term in the sum (8). The most probable energy U of the system S corresponding to this largest term can be found from the condition U S(U) E tot S R (E tot U) = 0. (9) This condition allows to introduce the temperature T characterizing the euilibrium: T = U S(U) = S R (E tot U). (0) E tot We can also consider the situation when the Boltzmann-Gibbs entropy of the system S is not necessarily extensive and proportional to the number of particles N in the system. If we introduce the generalized entropy as S (E) = ln W(E) then the sum (8) becomes W tot (E tot ) = E e k S B (E) e k S B R(E tot E) () and the largest term is determined from the condition ( e k S B (U) ) U S (U) E tot S R (E tot U) = 0. (2) Here we have used E. (B4). We can conclude, that the temperature T is related to the generalized entropy S via the euation T = Introducing the auxiliary -temperature by the euation U S (U) + S (U). (3) T = U S (U) (4) we get the relation ( T = T + ) S (U). (5) The auxiliary temperature T in the formulation of the non-extensive statistical mechanics based on maximization of entropy can appear as the inverse of the Lagrange multiplier associated with the energy constraint. Althoug T is not the physical temperature, it can have another physical meaning. For example, such effective temperature is directly related to the density of vortices in type II superconductors [26]. The relation (5) between the physical temperature T and the auxiliary temperature T has been proposed by various authors in Refs. [27 35]. Definitions of the temperature associated with different formulations of the non-extensive statistical mechanics have been analyzed in Ref. [36]. The general reuirement that composition rules of entropy and energy should satisfy to be compatible with zeroth law of thermodynamics has been investigated in Ref. [37]. It has been shown that formal logarithms of the original uantities should be additive. Note, that the statistics of the system S is determined by the reservoir, as we see from E. (8). Therefore, it is more convenient to describe even such a system using the Boltzmann-Gibbs entropy. The same conclusion has been made in Ref. [38]: it has been shown that physical temperature and pressure within the formalism for non-extensive

4 thermostatistics leads to expressions which coincide with those obtained by using the standard formalism of statistical mechanics. According to the central postulate of the statistical mechanics, the probability of the microstate µ µ R where the system S is in the microstate µ and the reservoir is in the microstate µ R is p(µ µ R ) = The probability of the microstate µ of the system S then is W(E tot ). (6) p(µ) = µ R p(µ µ R ). (7) 4 If the energy of the microstate µ is E µ then the energy of the reservoir is E tot E µ and the number of acceptable microstates of the reservoir is W R (E tot E µ ). We obtain that the probability of the microstate is eual to Approximating the number of microstates as p(µ) = W R(E tot E µ ) W(E tot ). (8) W R (E tot E) = e S R(E tot E) e S R(E tot) E E tot S R(E tot) (9) we obtain that the probability of the microstate of the system S is proportional to the Boltzmann factor ( P(E) = exp ) T E. (20) Here we used the definition of temperature (0): E tot S R (E tot ) T. Note, that for the justification of the exponential form of E. (20) it is essential that the Bolzmann-Gibbs entropy of the large reservoir had very small second derivative, 2 EtotS 2 R (E tot ) 0. This reuirement means that the heat capacity of the reservoir C R = T 2 2 S Etot 2 R (E tot ). (2) should be very large, that is the reservoir should be a thermostat. Other possible forms of entropy (for example, the generalized entropy S (R) = ln W R with ) do not have small second derivative and thus do not lead to a good approximation for the probability of microstate. From the Boltzmann factor (20) follows that the normalized probability of the microstate can be written as where p(µ) = Z e T Eµ, (22) Z = µ e T Eµ (23) is the partition function. The distribution of the energy of the system E is obtained multiplying the probability p(µ) by the number of microstates having energy E µ = E. This number is eual to W(E) = e k S(E) B, therefore the distribution of the energy is given by p(e) = Z e S(E) T E. (24) The probability p(e) should be normalized, thus the partition function can be also written as Z = E e S(E) T E. (25)

5 5 In the macroscopic limit the sum of large exponentials can be approximated by the largest term, therefore lnz S(U) T U, (26) where the energy U corresponding to the largest term is obtained from the euation The average energy of the system U S(U) = T. (27) Ū = µ E µ p(µ) (28) can be determined form the partition function: The free energy F is introduced according to the euation The euality defines the average entropy S. Combining Es. (29) (3) we get Ū = T 2 lnz. (29) T F = T lnz. (30) F Ū T S (3) S Ū = T. Due to approximation (26) in the macroscopic limit the average energy Ū coincides with the most probable energy U and the average entropy S coincides with S(U). III. CANONICAL ENSEMBLE IN NON-EXTENSIVE STATISTICAL MECHANICS Now let us consider the composite system where the large reservoir R is not described by the Boltzmann-Gibbs statistics. The Boltzmann-Gibbs entropy is not proportional to the number of particles N R in the reservoir and is not extensive. The extensive uantity is the generalized entropy S (R) for some value of : S (R) N R. Here k we consider only the situation when <. The number of microstates in such a reservoir is W R = e B N R. Therefore, the Boltzmann-Gibbs entropy = lnw R depends on the number of particles in the reservoir as S R lnn R. This expression for the Boltzmann-Gibbs entropy is similar to the entropy of the system consisting of d = uasi-particles, whereas the number of particles N R plays the role of the volume. A simple model of such a system has been presented in Ref. [39]: the model consist of a spin chain containing N R spins; spins next to each other have almost always the same direction, except there are d cases when the next spin has an opposite direction. In the Boltzmann-Gibbs statistic large reservoir has large heat capacity C R. Similarly, here we reuire that the -heat capacity of the reservoir, defined by E. (38), should be large. The approach presented in this Section is similar to the approach in Ref. [3]. However, in Ref. [3] the reservoir is incorrectly interpreted as a heath bath and having large heat capacity. As we have seen in the previous Section, such a reservoir leads to the exponential Boltzmann factor and approximation of the expansion of the number of states as a -exponential is not justified. Interaction of the system S weakly coupled to a finite reservoir having a finite energy has been considered in Ref. [40]. Under the assumption that the number of microstates of the reservoir having energy less than E R grows as a power-law of E R, the -exponential distribution of the energy of the system has been obtained. However, in Ref. [40] the extensivity of the generalized entropy has been not used and the parameter tends to when the number of particles of the reservoir increases. As in the previous Section, the system under consideration S is interacting with the reservoir via short-range interactions, thus the total energy of the composite system in the macroscopic limit is E tot = E + E R and the

6 total number of microstates in the composite system when the system S has energy E is W(E)W R (E tot E). Introducing the generalized entropy of the system S (E) = ln W(E) and the generalized entropy of the reservoir (E R ) = ln W R (E R ) we can write 6 W tot (E tot ) = E e k S B (E) k (E B tot E) e. (32) When < this sum cannot be approximated by the largest term. Approximation of a sum of large -exponentials is investigated in Appendix A. We assume that the postulate of eual probabilities of microstates in the euilibrium is valid also in the non-extensive statistical mechanics. When the postulate of eual probabilities of microstates is assumed, the statistics of the system S is determined by the reservoir according to E. (8). Therefore, even an ordinary system interacting with the reservoir having large -heat capacity is more conveniently described by the -entropy. The largest term in the sum (32) corresponds to the most probable state of the composite system and is found from the euation U S (U) + S (U) E tot S (R) (E tot U) + (E tot U) = 0. Here U is the most-probable energy of the system S. Thus in order to satisfy the zeroth law of thermodynamics we need to define the temperature T as T = U S (U) + S (U) = E tot S (R) (E tot U) + (E tot U). (33) This definition of the temperature is the same as E. (3). If one introduces the entropy of the combined system as S (tot) (E tot ) = ln W(E tot ) then as a conseuence of the impossibility to approximate the sum (32) by the largest term the entropy of the combined system is not a simple combination of the entropies of the system S and the reservoir R: S (tot) (E tot ) S (U) + S (R) (E tot U) + S (U)S (R) (E tot U). Due to this the conclusions of Ref. [4] that the zeroth law of thermodynamics holds only if the energy is also nonadditive does not apply for the situation considered in this paper. On the other hand, if one assumes that the interaction between the system S and R are long range and the energy is not additive then the pseudo-additivity of entropies can be valid [35]. The probability of a microstate of the system S is given by E. (8). Similarly as in the previous Section, assuming that the second derivative of -entropy of the reservoir is very small, 2 S (R) (E tot ) 0, we can approximate W R (E tot E) = e k B (E tot E) E 2 tot k B (E tot U) k (E U) B E (E tot tot U) e. (34) Note, that now the Boltzmann-Gibbs entropy of the reservoir does not have a small second derivative. The condition 2 Etot 2 (E tot ) 0 means that the -heat capacity of the reservoir, defined in E. (38), is very large. The ordinary heat capacity of such a reservoir can be determined as follows: if we increase the energy of the reservoir by a small amount E, the increase of the temperature T, according to E. (33), is T = E T This means that the heat capacity of the reservoir C R = E T is where T (R) C R = 2 ER 2 (E R ) E. (35) E R S (R) (E R ) + T T (R) C (R) is the auxiliary -temperature of the reservoir defined via the euation T (R), (36) = E R (E R ) (37)

7 7 and C (R) = (T (R) ) 2 2 ER 2 (E R ) (38) is the -heat capacity of the reservoir, defined similarly to the physical heat capacity, E. (2). If the second derivative of the generalized entropy is small, 2 ER 2 (E R ) 0, then the heat capacity is C = kb. If we increase the energy of the reservoir with very large -heat capacity by E, the new temperature of the reservoir becomes The expression for the heat capacity C = kb T = T + E. (39) is the same as the heat capacity of a gas consisting of d = uasi-particles. Using the property (B2) of the -exponential function we obtain that the probability of the microstate of the system S is proportional to the factor ) P(E) = exp ( (E U), (40) T(U) where T is the temperature according to E. (33). However, the temperature T depends also on the properties of the system, not only on the reservoir. It is more convenient to introduce the temperature that the reservoir not interacting with the system could have: T(0) = Taking into account that 2 ER 2 (E R ) 0 we obtain T(U) E tot S (R) (E tot ) + (E tot ). (4) + (E tot ) U E tot (E tot ) E tot S (R) (E tot ) = T(0) U. (42) This euation shows that the interaction with the system lowers the temperature of the reservoir. On the other hand, the -temperature of the reservoir, defined by E. (37) does not change. Using E. (42) we get that the probability of the microstate of the system S is proportional to ( P(E) = exp ) T(0) E. (43) An expression similar to E. (43) has been obtained in Ref. [3]. However, in Ref. [3] the temperature that enters P(E) has been interpreted as a physical temperature T, because the reservoir has been assumed to be a thermostat. The correct observation that the energy of the reservoir interacting with the system should decrease has been presented in Ref. [28]. One common objection to E. (43) is that this expression is not invariant to the change of zero of energies []. However, this reflects the physical situation of the system interacting with the reservoir having very large -heat capacity and, conseuently, small physical heat capacity. The zero of the energy of the system S is fixed by the reuirement that the energy of the reservoir should be E tot when E = 0. If we shift the energy zero by E, this is euivalent to the decrease of the energy of the reservoir by E. This decrease of the energy of the reservoir decreases the temperature. The probability of the microstate should remain the same, thus the new factor should be proportional to the old: ) ( P (E) = exp ( T (0) E P(E + E) = exp ) (E + E). (44) T(0) It follows that T (0) = T(0) E. (45)

8 This euation is consistent with E. (39). Similar argument has been presented in Ref. [40] by considering a system S interacting with a finite reservoir. Interesting feature of E. (43) is the presence of the cut-off energy: it follows from the definition of the -exponential function that P(E) becomes zero when E E max where E max = T(0). (46) This property of P(E) ensures that the physical temperature T is always positive. Discussion of possible cut-off prescriptions associated with Tsallis distributions is presented in Ref. [42]. Using the factor (43) we can write the normalized probability of the microstate as where p(µ) = Z e T(0) Eµ, (47) 8 Z = µ e T(0) Eµ (48) is the generalized partition function. The distribution of the energy of the system E is obtained multiplying the k S B (E) probability p(µ) by the number W(E) = e of microstates having energy E µ = E: p(e) = Z e k S B (E) e T(0) E T(E) T(0) S(E) T(0) E = e, (49) Z where T(E) = T(0) E. (50) The probability p(e) should be normalized, thus the partition function can be also written as Z = E T(E) T(0) S(E) T(0) E e. (5) The energy U corresponding to the largest term in the sum (5) is determined by the euation U S (U) + S (U) = T(U). (52) According to E. (42) the temperature T(U) coincides with the physical temperature. IV. GENERALIZED FREE ENERGY The probability proportional to the factor (43) admits several different possibilities to generalize the free energy. First of all, there are three possibilities corresponding to three temperatures: initial temperature of the reservoir T(0), auxiliary -temperature T and the physical temperature T(U). From those three choices only the temperature T(0) depends only on the reservoir and does not depend on the properties of the system. On the other hand, the temperature T(U) has a direct thermodynamical interpretation. In addition, the average energy of the system is connected to the generalized entropy with the parameter 2. The derivative of this generalized entropy with respect to average energy yields another auxiliary temperature T 2 and the corresponding generalized free energy. A. Initial temperature of the reservoir and unnormalized -averages Let us consider first the generalized free energy F corresponding to the temperature T(0). This choice is closely related to the approximation of the sum of large -exponentials and to unnormalized -averages. When < the sum

9 in E. (5) cannot be approximated by the largest term even in the macroscopic limit. The approximate expression for the sum of large -exponentials is obtained in Appendix A. According to the results of Appendix A and E. (5) the -logarithm of Z can be approximated as E ln Z E ( S (E) + S (E µ ) T(0) E ) p(e). (53) For any function of the energy f(e) the following euality holds: f(e)p(e) = ( ) k S f(e µ )e B (E µ) k S e B (E µ) p(µ) = f(e µ ) + E µ S (E µ ) p(µ). (54) Therefore, we can approximate the -logarithm of Z as ln Z ( ) S (E µ ) + µ S (E µ ) T(0) E µ p(µ). (55) This euation suggest to introduce the unnormalized -average energy of the system Ū = E µ p(µ) = ( E + ) S (E) p(e). (56) k µ B E This unnormalized -average of the energy can be determined from the generalized partition function Z using the euation In analogy to E. (30) we introduce the generalized free energy µ Ū = T(0) 2 T(0) ln Z. (57) F = T(0)ln Z. (58) 9 The euation F Ū T(0) S (59) defines the entropy S which is related to the unnormalized -average of the entropy S. Using Es. (57) (59) we obtain S Ū = T(0). (60) Entropy S can be calculated using the probabilities p(µ) according to E. (). Indeed, we have S = T(0) (Ū F ) = µ Expressing the energy from the probability p(µ) we get T(0) E µp(µ) + ln Z. (6) E µ = T(0)ln [p(µ)z ] = T(0) ( ln p(µ)+p(µ) ln Z ). (62) Inserting this expression for the energy E µ into E. (6) and taking into account that µ p(µ) = we obtain µ S = k p(µ) B. (63) This expression is consistent with the approximation (53). According to the approximation (53) the entropy S is S E S (E)p(E) = µ µ S (E µ ) + S (E µ ) p(µ). (64) In the macroscopic limit the entropy S (E) is large and we can approximate p(µ) S k µ p(µ) B. (65) This expression is the same as (63).

10 0 B. -temperature and normalized -averages It can be more convenient to deal with normalized -averages. The normalized -average of the energy is µ U = E µp(µ). (66) µ p(µ) Using E. (54) the normalized -average of the energy can be written as ) E (+ E S (E) p(e) U = ( (67) E + S (E) )p(e) The sums of the form E f(e)p(e) can be approximated as f(u) E p(e) +f (U) E (E U)p(E). The sum E (E U)p(E) is small, since close to the maximum U the probability p(e) is an even function of E U. We obtain that in the macroscopic limit the normalized -average of energy U should be close to the most probable energy U. Since the sum of large -exponentials cannot be approximated by the largest term when <, it is not possible to determine the most probable energy U or the average energy µ E µp(µ) knowing only the generalized partition function Z. However, it is possible to calculate U, which is close to U. On the other hand, the entropy S cannot be approximated by S (U). Using E. (63) the normalized -average U can be related to the unnormalized Ū via the euation U = Ū + S. (68) From E. (68) follows that the introduction of the normalized -average energy U allows to factorize the generalized partition function Z : Z = e The entropy of the combined system S (tot) S (tot) k S B T(0)Ū = e k S B e T(0) U using the expansion (34) can be written as = ln e k B (E tot U) Using E. (69) and assuming that U U we get that the entropy S (tot) combination of entropies from the non-extensive statistical mechanics: S (tot). (69) T(U) U e Z (70) can be expressed as the usual pseudo-additive (E tot U )+ S + (E tot U ) S. (7) However, in this euation the entropy S is not directly connected to the number of microstates of the system S. Let us introduce an auxiliary -temperature T of the system S via the euation T = S U. (72) Using Es. (60) and (68) we get ( T(U ) = T + ) S, (73) where T(U ) = T(0) U is the temperature of the reservoir corresponding to the energy of the system eual to U. Since S > 0, the -temperature is always smaller than the physical temperature T(U ). Note, that only physical temperatures of the system and the reservoir are eual. The -temperature of the system T is not eual to the -temperature of the reservoir T (R). We introduce the -analog of the free energy corresponding to the temperature T : F U T S. (74)

11 Then, using Es. (72) and (74) we get S = F T. (75) We define the -heat capacity of the system as C = U S 2 F = T = T T T T 2. (76) The physical heat capacity C can be determined as the derivative of U with respect to the physical temperature T(U ): C = U T(U ). (77) Using Es. (72) and (73) we get the euation that relates the physical heat capacity with the auxiliary -heat capacity: C =. (78) T(U ) T C + Since T < T(U ), from E. (78) follows that the physical heat capacity C is always smaller than the -heat capacity C. Generalized partition function Z related to the generalized free energy F is Z e F T. (79) Note that Z Z. Using Es. (72), (74) and (79) we get the expression for the energy U : U = T 2 T ln Z. (80) The generalized partition function Z cannot be directly expressed as a sum. However, Z can be connected to a sum of -exponentials as follows: we write the probability of the microstate in the form p(µ) = Z e T(U) (Eµ U), (8) where Z = µ e T(U) (Eµ U) (82) is related to Z via the euation k Z = Z e B T(U) U. (83) Using E. (69) we get k S Z = e B. (84) Therefore, ln Z = ln Z T U. (85)

12 2 C. Physical temperature and Rényi entropy The third possibility is to in introduce the free energy corresponding to the physical temperature T(U ). In order to do this let us consider another entropy, given by the euation k S S = lne B = k ( B ln + ) S. (86) The entropy S is more directly connected to the physical temperature T(U ). Indeed, using Es. (72) and (73) we get that the derivative of the entropy S gives the physical temperature: S = U T(U ). (87) From E. (63) it follows that ( ) S = ln p(µ). (88) µ Thus the entropy S is the Rényi entropy [43, 44]. Using the Rényi entropy we introduce the free energy corresponding to the physical temperature T(U ): F U T(U ) S. (89) Also in this case we retain the Legendre transformation structure. For example, using E. (87) we get F T(U ) = S. (90) D. Average energy of the system It is impossible to exactly determine the average energy of the system Ū = µ E µ p(µ) (9) knowing only the sum of -exponents Z. However, the knowledge of another sum µ ) 2 (e T(0) Eµ allows us to do so. Indeed, using the property of the -exponential function (B4) and the expression for the probability p(µ) E. (47) we get Ū = T(0) 2 (2 )Z T(0) µ ) 2 (e T(0) Eµ. (92) Instead of this sum we can use the generalized entropy () with the parameter = 2 : S 2 = µ p(µ)2. (93) Using the generalized entropy S 2 the expression for the average energy Ū becomes Ū = T(0) 2 ( (2 )Z T(0) Z2 ) S2. (94)

13 3 We can obtain another expression for the average energy Ū by inserting E µ from E. (62) into E. (9): Combining Es. (94) and (95) we obtain Ū = T(0)Z S 2 T(0)ln Z. (95) Therefore, yet another auxiliary temperature T 2, introduced by the euation is eual to Ū T(0) T(0)Z S 2 = 0. (96) 2 T(0) T 2 = S 2 Ū, (97) T 2 = T(0) Z 2. (98) This relation between temperatures is exactly the same as obtained by maximizing the entropy(93) with the constraint (9) []. The temperature T 2, similarly as the temperature T, depends not only on the reservoir but also on the properties of the system. The generalized free energy corresponding to the averageinternal energy of the system Ū and the temperature T 2 is This expression for the generalized free energy F 2 is similar to the expression F 2 = Ū T 2 S 2. (99) F = Ū (2 )T 2 S 2 (00) for the generalized free energy F that follows from E. (95). We see that in general F F 2. V. DISCUSSION In summary, we have demonstrated that a small system interacting with a large reservoir having large -heat capacity can be described by the non-extensive statistical mechanics. From the point of view of the ordinary statistics such a reservoiris similar to a gas of d = /( ) uasi-particles. The probability of the microstate of the system interacting with the reservoir via short-range forces is given by the -exponential function (43), instead of the exponential Boltzmann factor (20). Large -heat capacity of the reservoir leads to a small physical heat capacity, therefore the temperature in the euilibrium T depends both on the properties of the reservoir and the properties of the system. In order to avoid this inconvenience one can consider the temperature T(0) of the reservoir that is not interacting with the system or introduce an auxiliary -temperature T (R) that remains constant due to large -heat capacity of the reservoir. Small heat capacity of the reservoir does not allow to consider it as a thermostat, thus the description using the standard canonical ensemble of the statistical mechanics is not applicable. The treatment of the canonical ensemble presented in this paper allows us to obtain relations between the physical temperature T and the auxiliary -temperature T (73) as well as between the -heat capacity C and the physical heat capacity C (78). Sums of large exponentials often appear in the Boltzmann-Gibbs statistical mechanics. Such sums can be approximated by keeping only the largest term. Similarly, in the non-extensive statistical mechanics appear sums of large -exponentials. However, for such sums taking only the largest term is a very poor approximation. This is because the -exponential function with < does not decrease as fast as the exponential function. As a conseuence, the deviations from the most probable state in the non-extensive statistical mechanics are much larger than the deviations in the standard statistical mechanics. As it is shown in Appendix B, sums of large -exponentials are well approximated using -averages. This fact is one the reasons why -averages play such an important role in the non-extensive statistics. In this paper we considered the reservoir for which the generalized entropy with < is extensive. The case of > is more complicated, because one cannot take the macroscopic limit N. The investigation of the small system interacting with the reservoir characterized by > remains a task for the future.

14 4 One can easily show that the sum of large exponentials Appendix A: Sum of large -exponentials Z = W i= e Nφ(i) can be approximated by the largest term. Indeed, if φ max is the maximum of φ(i) then and e Nφmax Z We Nφmax 0 lnz N φ max lnw N. If W grows slower than exponentially with increasing N, then in the limit of large N the ratio lnw/n vanishes and we have Now let us consider the sum of large -exponentials lnz lim N N = φ max. Z = W i= (A) (A2) (A3) (A4) e Nφ(i), (A5) where <. In contrast to the sum of large exponentials, approximation of the sum of large -exponentials with the largest term is a poor one. We can construct a better approximation as follows: let us introduce the weights and the unnormalized -average By noticing that one can write the difference ln Z N φ as φ p(i) = enφ(i) Z W φ(i)p(i). i= Nφ(i) = ln [p(i)z ] = ln p(i)+p(i) ln Z (A6) (A7) (A8) W i= ln Z N φ = p(i). (A9) The sum i p(i) can have the largest possible value when all weights p(i) are eual. In such a case p(i) = /W and i p(i) = W. Thus If W grows with increasing N as N p and ln Z N φ ln W N. p < then in the limit of large N the ratio ln W/N vanishes and we have ln Z lim N N = φ. This gives the reuired approximation of the sum of -exponentials. Note, that ln W/N as the upper limit of the difference ln Z N φ is the worst case, when all terms in the sum are eual. For sufficiently fast decreasing terms the sum i p(i) can be bounded even for larger W. (A0) (A) (A2)

15 5 Appendix B: Some properties of -exponential function In this paper we have used the following properties of -exponential and -logarithm: multiplication of two - exponentials -logarithm of a product e x ey = e[+( )y]x+y e x+y = e x e y +( )x = e x+[+( )x]y, (B) x = e y +( )y e, (B2) ln xy = [+( )ln y]ln x+ln y = ln x+[+( )ln x]ln y, (B3) the derivatives of -exponential and -logarithm: d dx ex = (e x ), d dx ln x = x. The euations (B) (B5) can be easily derived using the definitions (5) and are presented in the Appendix A of Ref. []. (B4) (B5) [] C. Tsallis, Introduction to Nonextensive Statistical Mechanics Approaching a Complex World (Springer, New York, 2009). [2] C. Tsallis, Braz. J. Phys. 39, 337 (2009). [3] L. Telesca, Tectonophysics 494, 55 (200). [4] A. Adare et al., Phys. Rev. D 83, (20). [5] R. M. Pickup, R. Cywinski, C. Pappas, B. Farago, and P. Fouuet, Phys. Rev. Lett. 02, (2009). [6] E. Lutz and F. Renzoni, Nature Physics 9, 65 (203). [7] R. G. DeVoe, Phys. Rev. Lett. 02, (2009). [8] Z. Huang, G. Su, A. El Kaabouchi, Q. A. Wang, and J. Chen, J. Stat. Mech. 200, L0500 (200). [9] J. Prehl, C. Essex, and K. H. Hoffman, Entropy 4, 70 (202). [0] B. Liu and J. Goree, Phys. Rev. Lett. 00, (2008). [] O. Afsar and U. Tirnakli, EPL 0, (203). [2] U. Tirnakli, C. Tsallis, and C. Beck, Phys. Rev. E 79, (2009). [3] G. Ruiz, T. Bountis, and C. Tsallis, Int. J. Bifurcation Chaos 22, (202). [4] C. Beck and S. Miah, Phys. Rev. E 87, (203). [5] A. A. Budini, Phys. Rev. E 86, 009 (202). [6] J.-L. Du, J. Stat. Mech. 202, P02006 (202). [7] C. M. Gell-Mann and C. Tsallis, Nonextensive Entropy Interdisciplinary Applications (Oxford University Press, New York, 2004). [8] S. Abe, Astrophys. Space Sci. 305, 24 (2006). [9] S. Picoli, R. S. Mendes, L. C. Malacarne, and R. P. B. Santos, Braz. J. Phys. 39, 468 (2009). [20] R. Hanel and S. Thurner, EPL 93, (20). [2] R. Hanel, S. Thurner, and M. Gell-Mann, PNAS 08, 6390 (20). [22] L. D. Landau and E. M. Lifshitz, Statistical Physics, Part (Butterworth-Heinemann, 980). [23] P. T. Landsberg, Thermodynamics and Statistical Mechanics (Dover, New York, 204). [24] Q. A. Wang, Eur. Phys. J. B 26, 357 (2002). [25] M. Neuenberg, Phys. Rev. E 67, 0364 (2003). [26] F. D. Nobre, E. M. F. Curado, A. M. C. Souza, and R. F. S. Andrade, Phys. Rev. E 9, (205). [27] S. Abe, Physica A 269, 403 (999). [28] S. K. Rama, Phys. Lett. A 276, 03 (2000). [29] S. Abe, S. Martínez, F. Pennini, and A. Plastino, Phys. Lett. A 28, 26 (200). [30] S. Abe, Physica A 300, 47 (200). [3] S. Abe and A. K. Rajagopal, Europhys. Lett. 55, 6 (200). [32] S. Martínez, F. Pennini, and A. Plastino, Physica A 295, 246 (200). [33] S. Martínez, F. Pennini, and A. Plastino, Physica A 295, 46 (200). [34] M. Casas, S. Martínez, F. Pennini, and A. Plastino, Physica A 305, 4 (2002). [35] A. M. Scarfone, Phys. Lett. A 374, 270 (200). [36] Q. A. Wang, L. Nivanen, A. Le Méhauté, and M. Pezeril, Europhys. Lett. 65, 606 (2004).

16 [37] T. S. Biró and P. Ván, Phys. Rev. E 83, 0647 (20). [38] R. Toral, Physica A 37, 209 (2003). [39] J. Ruseckas, Phys. Lett. A 379, 654 (205). [40] A. R. Plastino and A. Plastino, Phys. Lett. A 93, 40 (994). [4] C. Ou and J. Chen, Physica A 370, 525 (2006). [42] A. M. Teweldeberhan, A. R. Plastino, and H. G. Miller, Phys. Lett. A 343, 7 (2005). [43] A. Rényi, Probability theory (North-Holland, Amsterdam, 970). [44] C. Tsallis, J. Stat. Phys. 52, 479 (988). 6