Arxiv:1705.09511V1 [Cond-Mat.Quant-Gas]

Home , Bose gas

arXiv:1705.09511v1 [cond-mat.quant-gas] 26 May 2017 etto h aitna fteipretFrigshas gas form Fermi repre- following imperfect number the the occupation of Hamiltonian the the In sentation spinless [5–12]. imperfect, so-called gases the quantum of model the considering interaction. dimension- the of the strength the on on and between system depends the balance of statistics ality possible more this accompanying and how The to interaction is relation question [4]. in effect and detailed crossover investigated resemblant superconductivity BCS-BEC intensely A of the been BCS-theory gas. has the boson recently in repulsive known similar a is behavior of thermodynamic that the to but to attraction lead Fermi the the counterbalance effectively for only that not expect thermodynamics will one statistics of can And, existence fermions. to the constructive attractive investigate a provides to plays and statistics is can role quantum study the repulsion, this extent of effective interparticle the what object the the The for counterbalance and to attraction. already role fermions, this decisive leads For which a in play gas principle, system. and Fermi the Pauli rules ideal equilibrium of the bosons stable collapse between however, the answer of attraction leads the existence sense pure bosons the the the other For out On each simple: attract a different. particles is of the becomes exists. if danger situation description hand, well no other is thermodynamic is the limit there the thermodynamic the other defined simplest fermions provided each the whether and repel collapse In particles, quantum bosons, point therein. or the references literature, if huge and has case, [1–3] has and e.g., interest gases see utmost quantum of interacting been always for thermodynamics of ntefloigw iea nwrt hs usin by questions these to answer an give we following the In existence the of problem subtle and fundamental The H F,imp = pnespitprilswoettlitratceintera interparticle total whose particles point spinless seatytesm steoeo h eusv mefc Bose imperfect repulsive the of one the as same the exactly is a hroyaisof thermodynamics ASnmes 53.k 53.p 67.10.Fj 05.30.Jp, 05.30.Fk, numbers: PACS gas. Bose ideal the for X = ecnie w-iesoa mefc trcieFrian Fermi attractive imperfect two-dimensional consider We k a .INTRODUCTION I. − F ~ a 2 = 2 F m k a ≤ 2 0 n mefc trcieFriadrplieBs gases Bose repulsive and Fermi attractive imperfect o emos and fermions, for 0 n bevstetemdnmceuvlneo h attracti the of equivalence thermodynamic the observes one k hroyai qiaec ftwo-dimensional of equivalence Thermodynamic − a 2 nvriyo asw atua5 203Wra,Poland Warsaw, 02-093 5, Pasteura Warsaw, of University F d N V mefc em a ean eldfie o 0 for defined well remains gas Fermi imperfect 2 = nttt fTertclPyis aut fPhysics of Faculty Physics, Theoretical of Institute 2 ae aiokwk n ao a Piasecki Jaros law Napiórkowski and Marek a , F a ≥ = Dtd etme 9 2018) 19, September (Dated: 0 a B , ≥ o oos eso hti pt fteatato the attraction the of spite in that show We bosons. for 0 (1) as d where d momentum o h aitna fterplieipretBs gas Bose imperfect constant repulsive coupling the the with of Hamiltonian the for atrexp factor xsec fΞ of existence n where esrstesrnt ftema-edatatv po- attractive mean-field the energy of tential strength the measures h rn aoia atto function partition of evaluation canonical The grand the conditions. boundary periodic assume ieipretBs a n hwterthermodynamic their repul- dimensions. show two the and in and gas equivalence discuss gas Bose we Fermi imperfect paper imperfect sive this attractive In the both discussed. [5–12] were the and bosons forces condensation Casimir repulsive Bose-Einstein the to alia, inter applied where, succesfully been has ytemperature by n dmninlielFrigsa xdvlm behaves volume fixed at gas Fermi ideal -dimensional k k ≤ ewr ntegadcnnclesml parametrized ensemble canonical grand the in work We to nryi ersne by represented is energy ction E N steocpto ubro n-atcesaewith state one-particle of number occupation the is 0 = X ∞ h rn aoia atto ucinmyexist. may function partition canonical grand the 2 0 =0 V ∼ β , I H MEFC EM GAS FERMI IMPERFECT THE II. { X 1 n a with gas eoe h oueo h system, the of volume the denotes 1 = . . . , N eusv oegsscnitn of consisting gases Bose repulsive d k } 1+ ~ βa ′ k ∞ /k F,imp exp 2 F d 2 − , V N h oe fipretqatmgases quantum imperfect of model The . B T a n hc it ttepsiiiyta for that possibility the at hints which 2 a T F " k hmclpotential chemical , eipretFrigsand gas Fermi imperfect ve B ( ≤ − N ssyidb h rsneo the of presence the by stymied is , ,V µ V, T, 0 = = β hc al noqeto h very the question into calls which a 2 / F

a 2 0 ≤ X V , − k .Teculn constant coupling The 1. .Tegon tt nryof energy state ground The ). nlgu xrsinholds expression Analogous . a a a 0 n F F = k nparticular, In . ( elcdby replaced aN h ǫ k 2 / 2 − 2 Ξ / πm 2 F,imp µ V ) µ and , − with volume , ( a ,V µ V, T, 2 − F N a N V B 2 = < !# = ) V P and 0 k and , , n (2) a k F , 2

In order to resolve this issue we use the identity and rewrite Ξ(T,V,µ) as

βaF Ξ (T,V,µ)= (4) exp N 2 = F,imp 2V 1 ∞ V 2 V q2 q 1 ∞ 2 2 dq exp − ΞF,id T,V,µ − , V V q 2πβaF 2βaF β dq exp − − Nq (3) −∞Z 2πβaF 2βaF −∞Z where ΞF,id(T,V,µ) denotes the grand canonical partition function of an ideal Fermi gas. For large V Eqs (4) can be represented as

ΞF,imp(T,V,µ) = 1 ∞ V 2 dq exp [−V ϕF (q; T,µ)], (5) 2πβaF −∞Z with

∞ d 1 λd−2 a 1 x 2 ϕ (q; T,µ)= 0 (q − βµ)2 − dx , (6) F λd 2 a d 1+ ex−q  F Γ( 2 + 1)  Z0  

2 where a0 = h /2πm and λ is the thermal de Broglie wave- III. THERMODYNAMICS OF THE 2 length, βa0 = λ . Thus the problem of the existence and TWO-DIMENSIONAL IMPERFECT FERMI GAS evaluation of the series in Eq. (2) has been rephrased as the question of the existence of the integral in Eq. (5). Once the conditions for its existence are settled the bulk thermodynamics can be extracted from it via the method of steepest descent. The thermodynamic limit of the grand canonical The parameter a0 turns out to play a prominent role ΩF,imp(T,V,µ) potential density lim = ωF (T,µ) = in our analysis and for this reason its presence in V →∞ V Eq. (6) is explicitly exposed. The existence of the −pF (T,µ) can be calculated using the method of steepest integral in Eq. (5) depends on the behavior of func- descent. It gives tion ϕF (q; T,µ) for large |q|. When q → −∞ one 2 2 ωF (T,µ)= kB T ϕF (q0(T,µ); T,µ), (7) has ϕF (q; T,µ)/q ∼ a0/2aF λ and this regime poses no problem. On the other hand for q → ∞ one has d where q0(T,µ) minimizes ϕF (q; T,µ) and fulﬁlls the equa- 2 2 −1 d d 2 ϕF (q; T,µ)/q ∼−q /λ Γ( 2 +2) + a0/2aF λ . Thus tion for d> 2 the grand canonical partition function ceases to aF q0 exist (note that this observation refers to fermions with q0 = βµ + ln (1 + e ). (8) a purely attractive interactions). For d = 2, when q → ∞ 0 2 2 one ﬁnds ϕF (q; T,µ)/q ∼ (a0 −aF )/2aF λ which delim- ∂ωF The number density nF (T,µ) = − ∂µ is thus re- its the allowed values of the coupling constant aF to the T lated to q0 via range 0 ≤ aF ≤ a0. For the boundary value aF = a0 βµ one has ϕF (q; T,µ) → − 2 q which means that in this λ q0(T,µ)= β(µ + aF nF (T,µ)) (9) particular case only negative values of chemical potential are allowed. This limitation will be shown to have a while the pressure pF (T,µ) is simple physical interpretation. Note that in two dimensions parameter a has a straightforward interpretation: aF 2 0 pF (T,µ)= − nF (T,µ)+ pF,id(T,µ + aF nF (T,µ)), the ground state energy per particle e = n(a − a )/2, 2 0 0 F (10) where n is the density, becomes negative for aF >a0. where pF,id(T,µ) denotes the pressure of the ideal Fermi gas

∞ k T p (T,µ)= B dx ln 1+ eβµ−x . (11) F,id λ2 Z 0 3

and 40 2 2 ∞ r r kBT nλ (−1) nλ2 30 pF (T,n)= − 2 + 2 e − 1 . λ " 2 r # nF r=1 X (17) n 20 F,id One notes that remarkable identities follow from the above formulas. They relate the imperfect Fermi gas at 10 aF = a0 and the ideal Bose gas [22] (both deﬁned for µ< 0) for which 0 -3 -2 -1 0 1 2 3 1 n (T,µ)= − ln 1 − eβµ (18) βμ B,id λ2 and FIG. 1. The plot of the imperfect Fermi gas density nF nor- ∞ βµ(T,n)r malized by the ideal Fermi gas density nF,id as a function of kBT e pB,id(T,n)= = βµ for a/a0 = 0.7, 0.85, 0.95, 0.99 (larger a/a0 values corre- λ2 r2 spond to larger densities). r=1 X 2 r ∞ −nλ k T 1 − e B . (19) λ2 r2 r=1 X The density nF (T,µ) is obtained by solving Eq. (8) for It follows from Eqs (16) and (18) that for a = a0 q0(T,µ) and inserting the result into Eq. (9) nF (T,µ)= nB,id(T,µ) (20) 1 q0(T,µ) nF (T,µ)= ln 1+ e , (12) λ2 while from Eqs (17) and (19) one ﬁnds see Fig. 1. The equation of state, pF (T,n) takes the form pF (T,n)= pB,id(T,n). (21) aF 2 pF (T,n)= − n + pF,id(T,n), (13) The above equality can be checked using the integral 2 representation of function pB,id(T,µ) in Eq. (19) and where pF,id(T,µ) in Eq. (14). A straightforward calculation

2 r leads to ∞ 1 − enλ kBT ∂ pF,id(T,n)= − 2 2 . (14) [p (T,n) − p (T,n)]=0 (22) λ r ∂n F B,id r=1 T X

The coefficient of isothermal compressibility from which Eq. (21) follows because 1 ∂n [pF (T,n) − pB,id(T,n)] |n=0 = 0. χT (T,n) = of the imperfect Fermi gas n ∂p T has the following form Thus according to Eqs (20) and (21) the two- 2 −1 dimensional imperfect Fermi gas at a = a0 and 2 n a0 χT (T,n)= n (a − aF )+ 2 . (15) the ideal Bose gas, both in the same state charac- 0 enλ − 1 terized by arbitrary T and µ < 0 have the same densities nF (T,µ) = nB,id(T,µ) and pressures For aF < a0 it remains finite. However, for aF = a0 the compressibility becomes infinitely large for n → ∞ pF (T,n)= pB,id(T,n). at fixed T , or equivalently for µ → 0− at fixed T (in this 2 limit nF λ = − ln (−βµ)). In order to put this remarkable equivalence into a broader perspective [14–19] we analyze the two- dimensional imperfect, spinless Bose gas with repulsive interactions characterized by the coupling constant aB > IV. THERMODYNAMIC EQUIVALENCE OF 0 TWO-DIMENSIONAL IMPERFECT FERMI AND 2 2 2 ~ k aB N BOSE GASES HB,imp = nk + . (23) k 2m 2 V X For the special case aF = a0, it follows from Eqs (8-11) As far as the existence of thermodynamics is concerned that there is no upper bound on the coupling constant aB. 1 The repulsive imperfect Bose gas has been intensely dis- n (T,µ)= − ln 1 − eβµ (16) F λ2 cussed in the literature, see e.g., [5–12, 21]. For d > 2 4 it shows Bose-Einstein condensation taking place for and thus the parameter ζ0 fulfills the same equation as −d µ>aBζ(d/2)λ with the critical indices belonging to parameter s0, see Eq. (26), provided the mean-spherical model universality class. The formal- ism used in [11, 12] to evaluate the partition function aB = a0 − aF . (33) is analogous to the one employed here with the identity In other words, if the above relation holds then in Eq. (3) taking the role of the Hubbard-Stratonovich the fermion and boson densities are identical: transformation [13] used in [11, 12]. Thus we only quote nF (T,µ) = nB(T,µ). Similarly, it follows from the relevant formulas. Analogously to the Fermi case, Eqs (10),(13),(19),(29) that the grand canonical free energy density ωB(T,µ) can be n2 obtained as the minimum of p (T,n)= (a − a )+ p (T,n)= F 2 0 F B,id 2 1 a0 2 n ϕB(s; T,µ) = − 2 (s − βµ) + g2 (exp(s)) (a0 − aF − aB)+ pB(T,n) (34) λ 2aB 2 (24) and thus when the relation (33) is fulfilled one has with respect to variable s, where g (z) is the Bose func- 2 p (T,n)= p (T,n). Thus for each value of the coupling tion, for details see [11, 12]. The thermodynamics follows F B constant a ≤ a there exists a repulsive imperfect Bose from the following set of equations F 0 gas characterized by the coupling constant aB = a0 −aF ; ωB(T,µ)= −pB(T,µ)= kBT ϕB(s0(T,µ); T,µ) (25) it has the same density and pressure as the Fermi gas. In particular, when aF = a0 one has aB = 0 and the where the parameter s0(T,µ) minimizing the function previously proven result nF (T,n) = nB,id(T,n) and ϕB(s; T,µ) solves the equation pF (T,n) = pB,id(T,n), see Eqs (20,21) is recovered. In this case the attractive mean field suppresses the aB s0 s0 − βµ = ln (1 − e ). (26) effect of the Fermi statistics and makes the pres- a 0 sure of the Fermi gas approach zero for T → 0 at The density nB(T,µ) is related to s0(T,µ) via any density n. On the other hand, for aF = 0 one 2 has pF,id(T,n) = a0n /2+ pB,id(T,n). These are spe- s0(T,µ)= β(µ − aBnB(T,µ)) (27) cial cases of the general equivalence condition in Eq. (33). and can be rewritten as To summarize, we have shown that thermodynam- 2 s0(T,µ) nB(T,µ)λ = − ln 1 − e (28) ics exists for the two-dimensional attractive imperfect Fermi gas provided the coupling constant aF measuring which gives pB(T,n) in the following form the strength of the attractive energy is not too large, 0 ≤ a ≤ a . No thermodynamics exists for three and a k T 2 F 0 p (T,n)= B n2 + B g 1 − e−nλ . (29) higher dimensional imperfect attractive Fermi systems. B 2 λ2 2 We showed that in two dimensions the density and pressure of the attractive Fermi gas are the same as those of the repulsive imperfect Bose gas provided the coupling The above set of equation Eqs (25-29) can be con- constant aB measuring the strength of repulsion fulfills fronted with the corresponding set Eqs (7-11) for the at- the relation aB = a0 − aF . Needless to say, the above tractive imperfect Fermi gas. In order to find the relation derived equivalence of thermodynamics of imperfect at- between the imperfect attractive Fermi and the imperfect tractive Fermi and repulsive Bose gases is restricted to repulsive Bose gas [4, 14–20] we rewrite Eq. (12) in the two-dimensions where the imperfect Bose gas does not following form suffer the Bose-Einstein condensation. 2 ζ0 nF λ = − ln 1 − e (30) where the parameter ζ0 is defined via

−q0 ζ0 = − ln 1+ e . (31) ACKNOWLEDGMENTS It follows from Eqs (31) and (8) that a M.N. acknowledges the support from National Science ζ − βµ = 1 − F ln 1 − eζ0 (32) Center, Poland via grant 2014/15/B/ST3/02212. 0 a 0

[1] J. L. Lebowitz and E. H. Lieb, Phys. Rev. Lett. 22, 631, [2] W. Thirring ed., The Stability of Matter: From Atoms (1969). to Stars, Selecta of Elliot H. Lieb, Springer, Heidelberg 5

(2005). [13] J. Hubbard, Phys. Rev. Lett. 3, 77, (1959). [3] E. H. Lieb and R. Seiringer, The Stability of Matter in [14] M. Randeria, Ji-Min Duan, and Lih-Yir Shieh, Phys. Quantum Mechanics, Cambridge University Press, Cam- Rev. Lett. 62, 981 (1989); and Phys. Rev. B 41, 327 bridge (2010). (1990). [4] W. Zwerger ed., The BCS-BEC Crossover and the Uni- [15] M. Drechsler and W. Zwerger, Ann. Phys. 1, 15 (1992). tary Fermi Gas, Lecture Notes in Physics 836, Springer, [16] G. Bertaina and S. Giorgini, Phys. Rev. Lett. 106, Heidelberg (2012). 110403 (2011). [5] K. Huang, Statistical Mechanics, Wiley, New York [17] V. Makhalov, K. Martiyanov, and A. Turlapov, Phys. (1963). Rev. Lett. 112, 045301 (2014). [6] E. B. Davies, Commun. Math. Phys. 28, 69 (1972). [18] M. Bauer, M. M. Paarish, and T. Enss, Phys. Rev. Lett. [7] E. Buffet and J. V. Pulè, J. Math. Phys. 24, 1608 (1983). 115, 135302 (2014). [8] V. A. Zagrebnov and J.-B. Bru, Phys. Rep. 350, 291 [19] P. A. Murthy, I. Boettcher, L. Bayha, M. Holzmann, D. (2001). Kedar, M. Neidig, M. G. Ries, A. N. Wenz, G. Zürn, and [9] J. T. Lewis, Statistical Mechanics and Field Theory: S. Jochim, Phys. Rev. Lett. 115, 010401 (2015). Mathematical Aspects, Lecture Notes in Physics 257, [20] E. R. Anderson and J. E. Drut, Phys. Rev. Lett. 115, 234, Springer, New York (1986). 115301, (2015). [10] M. van der Berg, J. T.Lewis, and P. de Smedt, J. Stat. [21] H. W. Diehl and S. B. Rutkevich, arXiv:1704.00507 Phys. 37, 697 (1984). (2017). [11] M. Napiórkowski and J. Piasecki, Phys. Rev. E 84, [22] R. M. Ziff, G. E. Uhlenbeck, and M. Kac, Phys. Rep. 32, 061105 (2011). 169 (1977). [12] M. Napiórkowski, P. Jakubczyk and K. Nowak, J. Stat. Mech. P06015 (2013).