arxiv:cond-mat/ v1 [cond-mat.str-el] 30 Mar 2004

Transcription

1 What is between Fermi-Dirac and Bose-Einstein Statistics? Krzysztof Byczuk a, Jozef Sa lek b,a, Geoffrey Joyce c, Sarben Sarkar c (a) Institute of Theoretical Physics, Warsaw University, ul. Hoża 69, -68 Warszawa, Poland (b) Institute of Physics, Jagellonian University, ul. Reymonta 4, 3-59 Kraków, Poland (c) Wheatstone Physics Laboratory, King s Collage, Strand, London WC2R 2LS, UK arxiv:cond-mat/43735v [cond-mat.str-el] 3 Mar 24 We overwiev the roerties of a quantum gas of articles with the intermediate statistics defined by Haldane. Although this statistics has no direct connection to the symmetry of the multiarticle wave function, the statistical distribution function interolates continuously between the Fermi- Dirac and the Bose-Einstein limits. We resent an exlicit solution of the transcendental equation for the didtribution function in a general case, as well as determine the thermodynamic roerties in both low- and high-temerature limits. Published in Acta Physica Polonica B 26, 267 (995) PACS numbers: INTRODUCTION Statistical roerties of quantum many-body systems are determined by calculating thermal averages with a roer distribution function, which reresents the mean number of articles occuying each of the single-article states. For the three dimensional systems of noninteracting articles this distribution function is given by either the Fermi-Dirac (FD) or the Bose-Einstein (BE) function, deending on whether the sin of the articles is half-integer or integer, resectively []. In the systems with either strong interactions or of low sace dimensionality (d < 3) the situation changes drastically. It turns out that the oulation of single article quantum level can be given by a function which is of neither FD, nor BE form. The reason for that may be twofold. First of all, in the sace of two dimensions the symmetry of the many body wave function is not necessarily either even or odd [2]. When two identical articles are exchanged the wave function acquires an arbitrary hase factor e iα, where α is so called statistical angle. In other words, the roer symmetry grou in this case is the braid grou for which the irreducible reresentations are χ(α) = e iα. The articles obeying the fractional statistics are called anyons. On the other hand, one can imagine that the interactions between articles are sufficiently strong, e.g. singular, that number of available effective single-article occuancies of each of the levels deends on the number of articles already resent in that state. This roerty can also lead to the dearture from FD and BE statistics in arbitrary satial dimensions [3]. There is no universal theory of statistical roerties of the systems described above. However, there are some exactly solvable models of interacting articles for which the momentum distribution functions are modified [4]. One of the best known cases is the one dimensional Calogero-Sutherland model for which the statistical roerties of quasiarticles change from BE to FD limits as the interaction increases [5]. Having in mind those exactly solvable examles and the fact that, in general, the statistics of elementary excitations (quasiarticles) in correlated systems can be arbitrary, we believe it is valuable to study the thermal roerties of articles with nonstandard distribution functions. In articular, in this aer we address the question how thermodynamic roerties evolve when the statistics changes from BE to FD limits. Our aer is organized as follows. In the first Section we introduce three interolation schemes between BE and FD statistics. Then, in Section 2 and 3 we discuss the high and the low temerature roerties of such systems. In Section 5 we resent general hase diagrams for those three systems of articles. ROUTES OF INTERPOLATING BETWEEN FERMI-DIRAC AND BOSE-EINSTEIN DISTRIBUTION FUNCTIONS For noninteracting bosons the momentum distribution function is n(ǫ i ) = whereas for ideal fermions we have n(ǫ i ) = e β(ǫi µ), () e β(ǫi µ) +. (2) In these functions, and hereinafter, ǫ i denotes the energy of a single article state i, β = /k B T, where T is the temerature, and µ is the chemical otential. Probably, the simlest modification of FD and BE functions can be carried out henomenologically in the following manner (route ) n γ (ǫ i ) = e β(ǫi µ) + γ, (3)

2 2 where γ is an arbitrary real number in the range (, ). The limit corresonds to BE statistics and γ = to FD one. For other allowed values of γ we have an intermediate statistics. In articular, the limit γ = corresonds to the classical, i.e. Boltzman, distribution function. Particles with arbitrary γ statistics are called γ-ons. We show later that this ad hoc roosal reresents a high-temerature aroximation to a general distribution function. As far as we know, it is imossible to derive the distribution function (3) exactly [6, 7]. However, starting from the modified Hilbert-Fock sace with creation and annihilation oerators obeying so called q-mutator algebra, i.e. a i a + j eiα a + j a i = for i j, there is a way to recover the function (3) aroximately where γ = cosα [7]. It is also interesting to observe that the distribution function (3) can be obtained if we suose that the number of ways of utting N i articles into G i states is given by w i = ( ) /γ (G i )!. (4) (N i )! (G i γn i )! Then, using the standard combinatoric aroach to quantum statistics [] we recover eq.(3). The thermodynamic otential Ξ γ (z, V, T) for γ-ons is given by Ξ γ (z, V, T) = i ( + γze βǫi ) /γ, (5) where z = e βµ is the fugacity. In eqs.(4,5) we suosed that γ. The case of γ = must be treated searately. The equation of state for such gas of ideal γ-ons is V k B T = ln Ξ γ(z, V, T) = γ ln( + γze βǫi ), (6) where is the ressure and V is the volume of the system. In order to eliminate z from this formula we need another equation which determines the total number of articles N N = z z ln Ξ γ(z, V, T) = i z + γ = n γ (ǫ i ). eβǫi i (7) In deriving the last equation we could also convince ourselves that the thermodynamic otential (5) is really consistent with the distribution function (3). Eqs. (6) and (7) describe the whole thermodynamics of γ-ons. We are going to exlore them in the next sections. On the other hand, one can interolate between FD and BE limits using the state counting arguments. This rovides us with the second ossibility (route 2). Namely, let the number of single article state D i in the i-th quantum level deends on the number of articles in this grou according to the formula [3, 8] i D i ({N j }) = G i g(n i ). (8) Also, suose that the total number of states of N-body system for a given configuration {N i } is [3] W({N i }) = i (G i g(n i ) + N i )! (N i )!(G i g(n i ) )!, (9) then using the Boltzman equation for the entroy S = k B ln Γ(ǫ), where Γ(ǫ) = {N i } W({N i}) is the total number of states, we can derive the momentum distribution function n g (ǫ i ) [8]. In eqs. (8) and (9) G i is the bare number of single article states (when no articles are resent). g is so called statistical interaction, or statistical arameter, which in fact arameterize the distribution function. We note that for g =, the number of states (9) is the same as for the ideal bosons. The case g = corresonds to free fermions. However, for g, we can again obtain an arbitrary distribution function. Physical interretation of g is follows. Ideal bosons do not obey exclusion rincile, so the number of states in the i-th energy grou would not change if we ut another article into this system and, therefore, D i ({N j }) = G i. Fermions obey Pauli exclusion rincile, which states that at most only one fermion can be resent in a given quantum state. Then, if we ut another fermion into the i-th grou the number of accessible state must decrease by one. Hence, D i ({N j }) = G i N i + and g =. For arbitrary g between zero and one we have a artial exclusion rincile and equation (9) directly leads to fractional exclusion statistics in arbitrary satial dimensions. Exlicitly, the distribution function for any g is given by the following transcendental equation [8] n g (ǫ i )e β(ǫi µ) = ( + ( g) n g (ǫ i )) g ( g n g (ǫ i )) g. () Of course, taking g = or g = in this equation we recover distribution functions () and (2). Particles obeying statistics with arbitrary g are called g-ons. From the last equation we find that the distribution function can be written in the following form n g (ǫ i ) = where w g (ǫ i ) must satisfy the equation w g (ǫ i ) + g, () e β(ǫi µ) = (w g (ǫ i ) + ) g (w g (ǫ i )) g. (2) As in the standard statistical mechanics [] we can find the thermodynamic otential for g-ons Ξ g (z, V, T) and, hence, the equation of state is V k B T = i ( + ( g) ng (ǫ i ) ln g n g (ǫ i ) The equation for the total number of g-ons is N = i ). (3) w g (ǫ i ) + g. (4)

3 3 The function n g (ǫ i ) and w g (ǫ i ) are given by eqs. () and (2) resectively. Again, these functions (3) and (4) with () and (2) determine the whole thermodynamics of g-ons. Finally, there is another way of changing FD and BE statistics (route 3). Namely, suose that the number of articles in each single article quantum state can be equal to, where is a natural number, i.e. N i =,, 2,... (5) The case = corresonds to free fermions and = gives BE statistics. The thermodynamic otential for a given is [9] Ξ (z, V, T) = e β(ǫi µ)ni = e (+)β(ǫi µ). e β(ǫi µ) i n i= i (6) To derive this we used the method of grand canonical ensemble. With the aid of the roerty n (ǫ i ) = z / z ln Ξ we can obtain the distribution function n (ǫ i ) = e β(ǫi µ) + e (+)β(ǫi µ). (7) The equation of state for a gas of so called -ons is V k B T = i and equation for their number is N = ( i ln z+ e βǫi(+) ze β ei, (8) e β(ǫi µ) + e (+)β(ǫi µ) ). (9) Let us summarize this Section. We introduced three schemes of modified quantum statistics and, as a result, we obtained three families of new nonstandard statistical distribution functions which are arameterized by some numbers (γ, g, ). In these cases the statistics interolate between FD and BE limits as the statistical arameters are changed. In the next sections we are going to discuss thermal roerties of such systems as a function of these arameters. HIGH TEMPERATURE PROPERTIES Let us consider first the ideal γ-ons with arabolic disersion relation ǫ k = 2 k 2 2m moving in the d-dimensional sace. Then, we can exand the functions (6) and (7) in owers of γz, and erforming the simle Gaussian integrals we find that k B T = λ d γ f d 2 +(γ, z), (2) n = N V = λ d γ f (γ, z), (2) d 2 where the function f k (γ, z) is defined as f k (γ, z) = i+ (γz)i ( ) i k. (22) i= λ = 2π 2 /mk B T is the thermal wavelength, and n is the article density. Of course, one can simly check that these equations reduce to those for ideal fermions or bosons as γ = or γ = resectively. Now, considering the high temerature limit, i.e. a limit for which nλ d γ, we can solve eq. (3) finding the fugacity z nλ d + γ(nλd ) 2 2 d/2. (23) Then, substituting this into eq.(2) we obtain the hightemerature (low-density) limit of the equation of state k B T = n( + Aγ 2n + ), (24) where A γ 2 is a second virial coefficient λd A γ 2 γ. (25) 2d/2+ The character of the first quantum correction to the classical gas deends on the sign of γ. For γ >, A γ 2 is ositive and, effectively, articles reel each other. In this case the ressure is higher than for the classical gas. For γ < the second virial coefficient is negative which means that the ressure is lowered. In this range of statistical arameter γ articles tend to collase because they effectively attract each other. Roughly seaking, for γ > the articles resemble fermions and for γ < they look like bosons. At the boundary (γ = ) between these two limits there is no quantum correction to the classical behavior. We can also calculate the internal energy U = (d/2)v and, hence, the secific heat of γ-ons at high temeratures. The result is C V = k Bd 2 n( + Aγ 2n( d/2) + ). (26) In the one dimensional case the secific heat is decreasing function of temerature for ositive γ. For d = 3 the behavior is oosite whereas, for d = 2 the secific heat does not deend on temerature. Also, we see that as T the secific heat goes to its classical value (d/2)r, where R is a gas constant (R = nk B ), and this result is indeendent of γ. Next, we examine the high temerature roerties of g-ons. The momentum distribution function for them is not given exlicitly but via the solution of eq. (). However, we can find exactly this solution in term of Taylor ower exansion []. To do this we solve eq.

4 4 () by utilizing the Legendre inversion theorem. The distribution function is in the form of the exansion n g (ǫ i ) = m= (mg + g m) m m! ( ) m, (27) ξ m+ i where ξ i = e β(ǫi µ) and (a) m = a(a+) (a+m ) = Γ(a + m)/γ(a). For details see Aendix A. Then it is straightforward to find the equation of state for ideal g- ons k B T = λ d f d/2+(g, z), (28) and the equation for the total number of articles n = N V = λ d f d/2(g, z). (29) Again, λ is the thermal wavelength and function f k (g, z) is defined as follows f k (g, z) = ( ) m (gm + g m) m z m+ m! (m + ) k. (3) m= Eliminating z from eq. (28) with the aid of eq. (29) we find the equation of state in the high temerature limit in the form k B T = n( + Ag 2n + ), (3) where the second virial coefficient for g-ons is A g 2 = (2g ) λ d 2 d/2+. (32) For g > /2 this coefficient is ositive, so g-ons are like fermions. However, for g < /2 they resemble bosons in that sense that they tend to condense. We also see that identifying statistical arameters γ and g, i.e. γ = 2g, (33) these two sorts of articles obeying fractional statistics are identical at high temeratures. Also, the whole results for the internal energy and the secific heat for γ- ons could be reeated for g-ons if we use the relation (33). Finally, we examine the high temerature roerties of -ons. Again, the distribution function (7) can be exanded in ower series and Gaussian integrals can be erformed easily for the ideal-gas case. We find the equation of state and equation for article number in the forms k B T = λ d g d/2+g(, z), (34) n = N V = λ d g d/2(, z). (35) In the resent case the function g k (, z) is following g k (, z) = m= ( ) mzm(+) m k + If we take > then the fugacity z is m= z m m k. (36) z nλ d (nλd ) 2, (37) 2d/2 where we droed all higher order terms. For = the sign in the second term on the right hand side of eq. (37) would be ositive. Substituting (37) into (34) and taking only terms of the order n 2 we have where k B T = n( + A 2n + ), (38) A 2 = 2 λd d/2+. (39) We see that the second virial coefficient for -ons does not deend on and is always negative for >. Therefore, we conclude that at high temeratures these articles resemble ordinary bosons. To summarize this Section, we notice that for the first two statistical distribution functions we obtained the same behavior in the high temerature limit, whereas in the third case the result was different. In other words, g-ons and γ-ons behave similarly whereas -ons always have the boson tye corrections to classical limit. In the next section we check the low temeratures roerties of those systems. LOW TEMPERATURE PROPERTIES Low temerature roerties of γ-ons must be discussed searately for ositive and negative γ. Let us consider first the γ > case. At T = the momentum distribution function reduces to the ste function n g (ǫ i ) = γ θ(µ ǫ i). (4) All states below the Fermi level are occuied with the mean oulation /γ. As γ this number diverges. For finite γ we see that the averaged number of articles at each quantum level could be fractional and greater than one. The states above the Fermi level are emty. For free articles we can find the Fermi energy ǫ γ F = µ(t = ). Namely, using eq. (4) we calculate that for a given density of articles ( ) 2 ǫ γ nd γ F = d, (4) 2Ω

5 5 where Ω = (2π) d/2 Γ(d/2) ( m ) d/2. The length of the Fermi 2 2mǫ vector k F = γ F deends on the value of γ and so 2 does the volume enclosed by the Fermi surface. In articular we see that as γ the Fermi vector disaears and all articles are on the lowest quantum level. The existence of the Fermi surface imlies that the low temerature roerties of a system are affected only by the low energy article-hole excitations close to it. In other words, only the articles from the vicinity of the Fermi surface give contributions to the thermodynamics in this limit. Therefore, we erform the low temerature (Sommerfeld tye) exansion for any thermal quantity. However, in order to do this we have to make some extra modification in the distribution function (3). Let us rewrite it as follows e β(ǫ µ) + γ = γ e β(ǫ µ ) +, (42) where µ µ + k B T lnγ. Now, we can exand the thermodynamic functions around the oint µ. The general scheme of calculation is resented in Aendix B. Here we quote only the final results. The lowest corrections to the chemical otential are µ = ǫ γ F k BT lnγ π2 6 d 2 2ǫ γ F γ (k BT) 2. (43) For ordinary fermions (γ = ) we had the first correction to the chemical otential roortional to T 2. Now, the linear in temerature correction has aeared. This is a direct consequence of the modified form of the distribution function (3). However, this term disaears in the formula for the internal energy. Namely, as was shown in Aendix B the internal energy is U V = ǫ γ + Ω γ ( ) 2 nd γ d (kb T) 2. (44) 2Ω The general form is similar to that for fermions. The only difference is the factor /γ which enhances the density of states at the Fermi level. Also, the secific heat C V = K γ T (45) is linear function of T. The coefficient K γ (/γ) 2/d diverges as γ. So, we roved that the ground state and the low temeratures roerties of γ-ons for γ > resemble roerties of ordinary fermions. However, it is not the case for γ <. For negative statistical arameter γ the distribution function (3) is singular when ǫ i + k B T lnγ µ =. This may be a signal of the Bose-Einstien condensation (BEC) at least in dimensions greater than two. Let us consider this case. The equation for total number of articles is n = n + λ d γ g /d+( γ, z), (46) where n is the number of articles at the first quantum level (with k = ) divided by V. The function g k is given by g k ( γ, z) = ( γ z) m m k. (47) m= BEC occurs when γ z = and then the state k = starts to be microscoically occuied. This gives that in three dimensional sace we have the following condition [] λ 3 n γ = (48) Hence, we find the critical temerature k B T γ c = 2π 2 m ( ) 2/3 n γ. (49) 2.62 This temerature is smaller for articles with fractional bosonic statistics than for ordinary bosons. Below T c the contribution to the internal energy comes from articles which are outside the condensate and we can easily find that U V ( ) d/2 Γ(d/2 + )ξ(d/2 + ) T = n(k B T) Γ(d/2)ξ(d/2) Tc γ. (5) Hence, the secific heat is C V = Γ(d/2 + )ξ(d/2 + ) Γ(d/2)ξ(d/2) and the equation of state is = d 2 ( d ) ( ) d/2 T k B Tc γ, (5) ( ) d/2 Γ(d/2 + )ξ(d/2 + ) T n(k B T) Γ(d/2)ξ(d/2) Tc γ. (52) The deendence on γ is hidden in Tc γ only. On eq. (49) we see that as γ then Tc γ too, and hence, the thermal functions (5) - (52) diverges. This exresses the fact that the limit γ = is a nonanalytic boundary between articles with γ > and with γ <. In other words, this articular value of γ establishes a true hase boundary between two hases of noninteracting γ-ons. In the next art we are going to discuss the low temerature roerties of g-ons. At T = the momentum distribution function reduces to n g (ǫ i ) = g Θ(µ ǫ i), (53) for all < g. This means that for all ositive g we have a true Fermi surface in the recirocal sace which

6 6 searates occuied and unoccuied states. Note also, that the case g = must be treated searately because then n g (ǫ i ) is infinite. Again, as in the case of γ-ons with γ >, the main contribution to the low energy roerties is given by the article-hole excitations around the Fermi surface. Performing the Sommerfeld-tye exansion we find that (see Aendix B and also []) µ = ǫ g F d 2 2ǫ g gc (g)(k B T) 2. (54) F In the resent case there is no linear temerature corrections to the chemical otential. We also calculate the internal energy U V = ǫ g + Ω(ǫ g F ) 2/d C (g)(k B T) 2, (55) where now the Fermi energy is ǫ g F = ( ndg 2Ω Hence, we obtain the secific heat ) 2 d. (56) C V = K g T, (57) where K g (/g) 2/d C (g). The low temerature roerties of g-ons resemble those of ideal fermions with modified density of states by the factor /g. We also see that the case g = can not be reached analytically from the case g >. This is due to the ossibility of BEC of ideal bosons (g = ). As an additional roblem it was interesting to consider BEC for ideal g-ons with g <. As we roved rigorously in the Aendix C, there is no BEC in this case. So, only true bosons can develo a macroscoic occuation of the lowest (k = ) quantum level at finite temeratures. Finally, we consider -ons at low temeratures. At T = the distribution function (7) is n (ǫ i ) = Θ(µ ǫ i ). (58) We see that, although both functions in eq. (7) diverge at ǫ i = µ the final result is well behaved; those two singularities cancel each other. The last result means that -ons ossess a shar Fermi surface with a wavevector k F = 2m 2 ( ) /d nd (59) 2Ω This roerty allows us to erform again the Sommerfeld (low temerature) exansion for thermodynamic quantities. This was done in Aendix B. We note that for any < the secific heat is linear with temerature, i.e. C V = K T, (6) where now K + 2/d. So, we conclude that at very low temeratures -ons resemble free fermions with modified density of states. CONCLUSIONS In this aer we studied the thermodynamic roerties of articles obeying three different statistics. Each of them interolates the distribution function between the BE and the FD limits. We comared both the high and the low temerature roerties of γ-ons, g-ons and -ons. We showed that γ-ons behave like fermions when γ > and like bosons when γ < at all temeratures. However, their roerties are not the same because their thermal functions are arameterized by γ which changes the number of states. We also showed that in between of the boselike and the fermilike regimes, there is a true boundary condition (γ = ), where the statistics is classical. However, when we aroached this limit from the left or the right sides the thermal functions diverge as /γ. In fig. we draw the schematic hase diagram for γ-ons for different values of γ. The thermodynamic roerties of g-ons turned out to be quite different. At high temeratures we noted two classes of behavior: fermilike for g > /2 and boselike for g < /2. For g = /2 the statistics was classical. However, at low temeratures g-ons resemble ordinary fermions with modified density of states for all < g. This means that there must be a crossover temerature when the tyical behavior changes its character. The schematic hase diagram is rovided in fig. 2. Finally, we considered -ons and we showed that for > they look like bosons at very high temerature. On the other hand, they resemble fermions at low temeratures if <. So, again there must exist a crossover temerature from the bosonic to the fermionic behavior. The hase diagram for -ons is shown in fig. 3. In conclusions, we showed that different routes of interolating between FD and BE limits corresonds to different ways of aroaching the statistical mechanics of articles with a nonstandard statistics. Our methods 2 and 3 are based on the combinatorial aroach, which is not directly related to the symmetry roerties of the many-body wave function with resect to the article transosition. In articular, the knowledge of the second virial coefficient is not enough to determine the full thermodynamics. For examle, although we know the virial coefficient for anyons, it is imossible distinguish which statistics describes them at low temeratures, γ-ons or g-ons. We think that none of them is the roer one. However, anyons in the strong magnetic field seem to be described by g-on statistical mechanics [2]. APPENDIX A In this Aendix we show how to solve eq. () using the Lagrange inversion theorem. Suose, that we have

7 7 a transcendental equation y = a + xφ(y), (6) and we want to find y as an unknown variable. The variables x and a are suosed to be given. The solution can be reresented by infinite ower series y = a + n= x n d n n! dy n [φn (a)]. (62) Noting, that eq. () has the same formal structure as (6) with a = we use eq. (62) to obtain eq. (27). arbitrary statistics. The only limitation is that this thermal function must be well behaved in the vicinity of chemical otential. Let us consider the following integral I = H(ǫ)n(ǫ µ)dǫ (63) APPENDIX B In this Aendix we resent the low temerature exansion for any thermodynamic function calculated with Introducing the new variable y = βǫ βµ and maniulate with the integrals we obtain I = β dyh( y β + µ)n(y) + dyh( y β βµ β + µ)(n(y) a) + β βµ dyh( y + µ)a, (64) β where the constant a must be chosen roerly to the statistics. For examle, for γ-ons we should take a =, g-on statistics needs a = /g whereas for -ons a =. At low temeratures we can aroximate the second integral extending the lower limit of integration to. Then, changing the variable in this integral y y and exanding both functions H(µ±k B Ty) around the oint µ we obtain I = a µ H(ǫ)dǫ + i= C i (a) d i H(µ) i! dǫ i (k B T) i+, (65) where the coefficients C i (a) are defined as follows C i (a) = dyy i (n(y) + ( ) i (n( y) a)). (66) If a = and n(y) is γ-ons momentum distribution function we obtain that the even coefficients are zero C 2n () = and the odd coefficients are given by C 2n () = ( ) n (2π) 2n B 2n 2n ( 2 2n ), (67) where B 2n are the Brilouin numbers. So, γ-ons have almost the same exansion as ideal fermions. Similar results we obtain for -ons. Using the distribution function (7) we find that C 2n () = ( ) n (2π) 2n B 2n 2n ( ( + ) 2n ) (68) and C 2n =. Note that for = we obtain the coefficients for the Sommerfeld exansion of ideal bosons rovided that they do not condense. For g-ons it is imossible to find the analytic formula for C k (g). We only know that C (g) = [] and the rest ones are finite. The formalism develoed above is useful in calculating the chemical otential and the internal energy. Namely, the equation for µ is and for the energy where n = N V = Ω dǫǫ d/2 n(ǫ), (69) Ω = E V = Ω dǫǫ d/2 n(ǫ), (7) ( m ) d/2. (7) (2π) d/2 Γ(d/2) 2 Exanding eq. (69) to first nontrivial order as in eq. (65) we have n Ω(a µ dǫǫ d/2 + C (a)(d/2 )(µ) d/2 2 (k B T) 2 ), (72)

8 8 and then suosing that the correction to the chemical otential is small µ = ǫ a F + δµ we can find that ( µ = ǫ a F (d/2 ) C ( ) ) 2 (a) kb T. (73) a Similarly, we find the low temerature internal energy E V = ǫ + Ω(ǫa F )d/2 C (a)(k B T) 2, (74) where ǫ a = 2Ωa d+2 (ǫa F )d/2+. The last equation allows us to calculate the secific heat which is where K a = 2Ω( nd 2aΩ ) 2/d C (a)k 2 B ǫ a F C V = K a T, (75) APPENDIX C In this Aendix we extend our discussion of g-ons to the case with negative g (i.e. we suose that g ). In rincile, the value of g can be arbitrary and it seems very interesting to investigate how the articles look like with attractive statistical interaction. The most exciting question deals with the BEC in such system. However, as we show below, the condensation does not occur in such gas. The equation for the momentum distribution function is following (cf. ()) n g (ǫ i )ξ = ( + ( + η) n g (ǫ i )) +η ( + η n g (ǫ i )) η, (76) where we define g = η, with η >, and ξ = e β(ǫi µ). Using the method develoed in ref. [] (see also Aendix A) we find the solution of eq. (76) in the form of the infinite ower series, i.e. n g (ξ) = m= (mη + η + m) m m!. (77) ξm+ This series is absolutely converged for ξ > (η+2)η+2 (η+) η+. Next, we find the equation for the total article number N V = V n + λ d g d/2(η, z), (78) where the function g k (η, z) is given by the following series g k (η, z) = m= (ηm + η + m) m m! z m+ (m + ) k. (79) n is the oulation of the lowest energy level and we treat it searately in case of BEC. Eq. (78) is absolutely converged for all z in the range (η + )η+ z z R. (8) (η + 2) η+2 On the other hand, momentum distribution function can be written in the form () where w(ξ) obeys the equation w η ( + w) +η = z eβǫ. (8) Condition for BEC is w = η. Substituting this into eq. (76) and taking the maximal value of ossible z = z R we obtain the allowed ǫ βǫ = 2( + η)ln( + η) (2 + η)ln(2 + η) η lnη. (82) We see that in order to have BEC ǫ must be negative but we have suosed that ǫ, so we obtain contradiction and this roves the absence of BEC for g-ons. [] K.Huang, Statistical mechanics, (John Wiley and sons, 963). [2] J.Leinnas, J.Myrheim, Nuovo Cimento B 37, (977). [3] F.D.M.Haldane, Phys.Rev.Lett. 67, 937 (99). [4] J.Sa lek, W.Wojcik, Phys.Rev. B 37, 532 (988), K.Byczuk, J.Sa lek, Phys.Rev. B 5, 43 (994). [5] D.Bernard, Y.S.Wu, in Proc. of the 6th Nakai Worksho, ed. M.L.Ge et al, (World Scientific, 995). [6] R.Acharya, P. Narayana Swamy, J.Phys. A 27, 7247 (994). [7] K.Byczuk, J.Sa lek, unublished work (994). [8] Y.S.Wu, Phys.Rev.Lett. 73, 922 (994). [9] This statistics was introduced by G.Gentile (94), but reroduced by many authors afterwards. [] G.S.Joyce, S.Sarkar, J.Sa lek, K.Byczuk, Phys.Rev. B 53 (in ress) (996). [] C.Nayak, F.Wilczek, Phys.Rev.Lett. 73, 274 (994). [2] M.D.Johnson, G.S.Canright, Phys.Rev. B 49, 2947 (994). FIGURE CAPTIONS Fig.. Schematic hase diagram for γ-ons. Fig.2. Schematic hase diagram for -ons. Fig.3. Schematic hase diagram for -ons.