![Arxiv:1612.02496V1 [Cond-Mat.Quant-Gas]](http://data.docslib.org/img/af815560ce2b3c72917562b665f46cbe-1.webp)
arXiv:1612.02496v1 [cond-mat.quant-gas] 8 Dec 2016 oPuiecuinpicpeado h elgbeeffects negligible the of and of principle exclusion Pauli to 4.Ide,tenal da iuto a enexperi- suppression been the [16– has of of advantage well situation taking as ideal by ones nearly realized sys- ideal mentally the fermion fifteen trapped Indeed, interacting of last of 34]. also the only but not over [12–15] study tems interest, the renewed in [7– years, a gases Fermi triggered atomic 11] trapped in degeneracy quantum behavior. ther- singular the dimensions possible study low of to in model modynamics solvable IFG interesting the an thus becomes any [4–6], towards interaction unstable dimensions, inherently finite IFG. are low an systems in as Fermi degree dramatically since some changes in situations described short The be a can theory through forces liquid repulsively ideal range Fermi interacting gas the Landau’s Fermi fermions of on scope where, ideal relies explicative model the the gas of Fermi of terms success The in (IFG). ones— scales astrophysical sub-nuclear from density, to particle the range of wide values a under- of phenomena—in qualitative physical the several allowed of [3], independently standing Dirac made and [2] [1], Fermi principle by exclusion Pauli’s rate n,a nmn ntne,acmlt nesadn of understanding complete a importance instances, great many of of in shape resulted as based the has and, Techniques changing potential by trapping further entropy [36–39]. the of cool fermions giving-away to the atomic devised experimen- on of trend, being cloud this are a On techniques new achieve interest. to tally of needed to phenomena are models, and the temperatures matter [35], lower condensed variety much blocking probe a though Pauli to test as experiments directly such design to effects set- possibility quantum experimental the of the open on has achieved tings control the Further, ∗ fjsevilla@fisica.unam.mx nteohrhn,teeprmna elzto of realization experimental the hand, other the On incorpo- that statistics quantum the of discovery The p s wv cteigfrtetmeauerne involved. ranges temperature the for scattering -wave wv cteigi pnplrzdfringssdue gases fermion spin-polarized in scattering -wave atcedniya srvae ntecrepnigthermo characterized corresponding the are in gas revealed Fermi is as trapped nonmono density appearance a the particle the of of interpret indicator states to an equilibrium as allows potenti temperature, argument on chemical an potential the thes Such which of at one presented. temperature of vanishing relevance non the the favors physically thermo that the compressibi argument in (isothermal dimensionality, susceptibilities low at thermodynamic emerge, that peratures oe a oetasaeaaye.Priua teto is attention Particular analyzed. are potentials law power h ffcso o iesoaiyo h hroyaiso F a of thermodynamics the on dimensionality low of effects The .INTRODUCTION I. hroyaiso o-iesoa rpe em gases Fermi trapped low-dimensional of Thermodynamics nttt eFıia nvria ainlAtooad M´ Aut´onoma de Nacional F´ısica, Universidad de Instituto po otl2-6,000MeioDF,MEXICO M´exico D.F., 01000 20-364, Postal Apdo. rnic .Sevilla J. Francisco eta eut ob oooi eraigfnto of function decreasing the monotonic from po- diminishing a chemical temperature, the be the to discussed— potential, results box widely tential impenetrable the a con- IFG by For three-dimensional fined lev- the namely different [41–52]. case, at textbook— contexts importance and and/or meaning els its of sion tzr eprtr,t h auso h da classi- ideal the of than values larger much the temperatures for to gas cal temperature, zero at λ where osatdvddb 2 by divided constant ncmaio,teceia oeta fteIGvan- criti- BEC, IBG the of called the temperature temperature, cal [41]. of characteristic potential a Ref. below chemical in ishes Dickerson the and comparison, be- Cook In this by of presented argument is physical havior qualitative, clear, A ature. where Broglie, h perneo h odne hs hntegsis gas the when phase condensed the cooled. of appearance the rcs tacaatrsi temperature characteristic condensation-like where a a suffer 40], at can [28, IFG process the Refs. that in suggested out is pointed it been has interesting as thermo- and precludes behavior density not smooth does particle however, the a this of temperature, shows function IFG as behavior the dynamic (BEC) dimensions, condensation three Bose-Einstein in so-called the suffers of result value. would great atoms fermionic non-interacting trapped t ftesystem the of ity gas. ideal classical the the to of asymptotically values converging temperatures larger for o cu.A hsqik ieso,teceia po- decreases chemical gases, ideal the Bose and dimension, Fermi quirky the both, this does of At tential transition BEC occur. the not therefore finite and any at [53] order temperature off-diagonal-long-range no shows IBG hc enstecaatrsi temperature potential, characteristic chemical the the of defines sign author which of the this change where the of [40], that support suggest Ref. in in presented approach are thermodynamic phenomenon a on based T rl,tesgicneof significance the Truly, ndsicinwt h da oegs(B) which (IBG), gas Bose ideal the with distinction In hspcuecagsdaaial stedimensional- the as dramatically changes picture This iyadseicha) An heat). specific and lity = yaissusceptibilities. dynamics k p ∗ B fatemdnmcrgm,weethe where regime, thermodynamic a of hrceitctmeaue,namely, temperatures, characteristic e yai ucin fsaeadi the in and state of functions dynamic ie odffrn hrceitctem- characteristic different to given 2 ylre utain neeg and energy in fluctuations larger by is value, energy Fermi the reaches al steBlzansconstant, Boltzmann’s the is π ~ oi eedneo h chemical the of dependence tonic 2 /mk T rigstapdb isotropic by trapped gas ermi eoe h ytmsaslt temper- absolute system’s the denotes d B slwrd ntodmnin the dimensions two In lowered. is , T π stetemlwvlnt fde of wavelength thermal the is exico, , T m c , h aso h atceand particle the of mass the µ n erae monotonically decreases and energy-entropy a oiae h discus- the motivated has em energy Fermi ~ T stePlanck’s the is k 0 B − Arguments . 1 ( T ~ 0 2 , /mλ marks , E T 2 F ), , 2 monotonically with temperature essentially in the same particle density of states (DOS) of the form [28, 59] functional way [54], being different only by an additive d/s−1 constant, expressly, the Fermi energy. This results in the g(ε)= Gd,s ε , (1) same temperature dependence of their respective specific where ε denotes the energy, Gd,s and s are positive con- heats at constant volume CV [54–56]. In general, this last outstanding feature occurs whenever the number of stants, the former depends on d and on the specific energy energy levels per energy interval is uniform as in the case spectrum of the system, while the later is determined by of a one dimensional gas in an harmonic trap [57, 58], the particular system dynamics. or the case s = d where s is the exponent of the single- Two instances lead to the power-law dependance in s expression (1): the first one is based on the general- particle energy spectrum of the form ε p , p being the s particle momentum [59]. ∝ ized energy-momentum relation [59, 66] εk = sk , k C k In one dimension, the chemical potential of the IBG being the magnitude of the particle wave-vector and decreases monotonically with temperature, and as in the s > 0 is a constant whose particular form depends onC s. The physical cases s = 2, 1 correspond, respec- two dimensional case, this behavior is related to the im- ~2 possibility of BEC as shown by Hohenberg [53], at fi- tively, to the nonrelativistic IFG with C2 = 2m and ~ nite, non-zero, temperature. In contrast, the chemical to the ultrarelativistic IFG for which C1 = c , c be- potential of the IFG exhibits a nonmonotonic behavior: ing the speed of light. In this case Gd,s takes the form d−1 d/2 d/s starts rising quadratically with T above the Fermi energy /[2 π Γ(d/2)sCs ], with Γ(σ) the gamma function instead of decreasing from it, and returns to its usual Vand = Ld the volume of the system. The second in- monotonic-decreasing behavior at temperatures that can stanceV is based on the d-dimensional IFG trapped by α be as large as twice the Fermi temperature (see Fig. 1 an isotropic potential of the form U(r) = U0 (r/r0) , below, see also Fig. 1 in Ref. [60]). This unexpected, and where U0, r0, are two constants that characterize the en- not well understood behavior, can be exhibited math- ergy and length scales of the trap. This trapping po- ematically by the Sommerfeld expansion [60, 61] or by tential leads, in the semi-classical approximation [67], 2 d/2 other methods [62–64], though no intuitive physical ex- d(1/2−1/s) to G = (2/s−1)Γ[d(1/s−1/2)] mr0 U with planation of it, that predicts its appearance in the more d,s Γ(d/2)Γ[d/s]~d 2 0 general case, seems to have been given before [65]. This s−1 =1/2+ α−1. Notice that in the later case, one can forms the basis for the motivation of the present paper. immediately establish the thermodynamic equivalence between the IBG and the IFG, namely, α = 2d/(2 d), After this excursus, one may conceive dimension two − as a crossover value for which the thermodynamic prop- implying that no such equivalence is possible in dimen- erties of ideal quantum gases are conspicuously distinct sions d> 2 for positive α. The equivalence does occur in for d > 2 than those for d < 2. This can be seen in two dimensions if α , which corresponds to the infi- → ∞ the specific heat, which in the case of the IFG exhibits a nite well potential and in one dimension if α = 2, which no-bump bump transition as dimension is varied from corresponds to the harmonic potential. 3 to 1 [60]→ analogous to the well known cusp no-cusp The thermodynamical properties of the ideal quan- transition of the IBG specific heat. In the later→ case, the tum gases are easily computed from the grand potential Ω(T, ,µ) U TS µN [68, 69], where U, and S denote cusp marks the BEC phase transition while no physical V ≡ − − meaning is yet given for the bump in the former case. the internal energy, and entropy respectively. For the trapped gas, denotes the appropriate thermodynamic In this paper we provide an analysis that attempts to V explain the various features that are observed in the low- variable that generalizes the volume of a fluid in a rigid- dimensional, trapped IFG, focusing in the nonmonotonic walls container (see Ref. [70] for the case of the three- dependence on T of the chemical potential. In section II dimensional harmonic trap), which in this paper is taken 2 d/2 the system under consideration is described, thermody- as = mr0 U d(1/2−1/s) which reduces to = ω−d namics quantities are calculated and characteristics tem- V 2 0 V r ~ 2 peratures are introduced. In section III a heuristic ex- for the isotropic harmonic trap U( ) = ω (r/r0) with 1/2 planation of the nonmonotonic dependence of the IFG r0 = (2~/mω) . For a gas of noninteracting fermions, chemical potential on temperature is given. In sections Ω(T, ,µ) can be written in the thermodynamic limit, IV and V the physical meaning of two relevant charac- N V , with N/ = constant, [71] as teristic temperatures is given. Finally, conclusions and → ∞ V → ∞ V final remarks conform section VI. ∞ Ω(T, ,µ)= k T dε g(ε) V − B × Z0 ln [exp β(ε µ) fF D(ε,T )] , (2) II. GENERAL RELATIONS, CALCULATION { − } OF THE CHEMICAL POTENTIAL AND THE −1 where fF D(ε,T ) = exp[β(ε µ)]+1 is the Fermi- THERMODYNAMICAL SUSCEPTIBILITIES Dirac distribution function{ that− gives} the average occu- pation of the single-particle energy state ε at absolute We consider an IFG of N, conserved, spinless fermions temperature T . As usual, β denotes the inverse of the in arbitrary dimension d > 0. We assume a single- product of T and kB the Boltzmann’s constant . 3
F by d/s = 1/4 a) d/s = 1/2 /E d/s = 3/4
µ 2 d/s = 1 d d/s = 3/2 S/ = k G Γ(d/s) (k T )d/s +1 d/s = 2 V − B d,s B s × d/s = 3 1 µ Li ( eβµ) Li ( eβµ) . (5) d/s+1 − − k T d/s − B 0 In the top panel of Fig. 1 the temperature dependence of the ratio µ/EF is shown for different values of the ra- -1 tio d/s and for N/ fixed, where T denotes the Fermi V F temperature defined through the relation EF = kBTF , s/d d s/d -2 where EF is explicitly given by (N/ ) in 0 1 2 3 4 5 sGd,s V d dimensions. For d/s < 1 the nonmonotonic depen- T/TF dence on temperature is clearly shown (the dashed line 5 F 10 corresponds to the case d/s = 1/2, while d/s = 1/4 is
/T b) presented with the only purpose of making the effects 4 * 10 of the system dimensionality more conspicuous). In the T 3 limit of high temperatures, T TF , the classical result 10 d/s ≫ µ kBT ln (T/TF ) Γ(d/s + 1) is recovered. 2 → 10 As occursh for the 2D ideal gas ini a box potential (s = 1 d = 2), the DOS is a constant whenever s = d, and the 10 0 chemical potential has the well known analytical depen- Tµ/ T T / TF 0 F dence on the temperature µ = E +k T ln 1 e−TF /T . 10 F B − For T TF , the chemical potential lies below the Fermi -1 10 ~ ~ energy≪ by a negligible, exponentially small correction. The low temperature behavior of µ for d = s can be ob- 0 0.2 0.4 0.6 0.8 1 tained approximately as a direct application6 of the Som- d/s merfeld expansion for T TF (see Ref. [73] pp. 45-46), namely ≪ FIG. 1. (Color online) a) Dimensionless chemical potential 2 µ/EF as function of the dimensionless temperature T/TF for 2 π d T 4 different values of d/s. The crossings with the horizontal lines µ EF 1 1 + [T/TF ] . (6) ∗ 0 ≃ − 6 s − TF O µ/EF = 1 and µ/EF = 0 marks the temperatures T and T " # ∗ respectively. b) The temperature T as function of s/d (solid 0 line), additionally the temperatures T (dashed line) and Tµ The power-law dependence on ε in expression (1) is man- (dash-dotted line) are included for comparison. ifested itself in the last expression, where the ratio d/s appears explicitly. Clearly, for d/s < 1, the chemical po- tential rises from the Fermi energy quadratically with T , and the non-monotonousness is a result of the fact that The average number of fermions N(T, ,µ) in the sys- for large enough temperatures, µ(T ) falls down with tem- tem is given by (∂Ω/∂µ) [68] whichV gives − T,V perature to negative values close to those of the classical gas. As a consequence of this “turning around”, µ(T ) d/s βµ develops a maximum at temperature Tµ and equals EF N/ = Gd,s Γ(d/s) (kB T ) Lid/s( e ), (3) V − − at two distinct temperatures, at T ∗ and 0, if d/s < 1, and only at T = 0 otherwise. Thus, the solution to ∞ l σ where Liσ(z)= l=1 z /l is the polylogarithm function the equation µ(T ) = EF as function of the parameter of order σ [72] and Gd,s = Gd,s/ . Expression (3) relates d/s, bifurcates at the critical value d = s as is shown N and µ, and forP fixed N, theV chemical potential is a in the bottom panel of Fig. 1. Note that for s = 2 ∗ function of the system temperature and volume. The and d = 1, T is as large as 1.896 TF and diverges as internal energy U(T, ,µ) per volume is given by d/s 0. This can be shown straightforwardly from Eq. V → ∗ (3) by putting µ = EF , since then, T must satisfies the ∗ −1 −TF /T d/s+1 βµ equation 1 = 1+ e in that limit. U/ = Gd,s Γ(d/s +1)(kB T ) Lid/s+1( e ), V − − In addition, the temperatures T 0 and T , that mark (4) µ the change of sign of µ and its maximum, respectively, are also shown in the bottom panel of Fig. 1 (dashed while the entropy S(T, ,µ)= (∂Ω/∂T ) per volume line and dashed-dotted line). T 0 is determined from the V − V,µ 4
0 1.5 equation µ(T )=0, which explicitly gives a) −s/d 1.25 T 0 = Γ(d/s + 1) ζ(d/s)(1 21−d/s) T , (7) − F /s) 1 h i B this expression gives the approximated values 3.48 TF , d/s = 1/4 1.44 TF and 0.989 TF for d/s = 1/2, 1, 3/2 respectively. 0.75 d/s = 1/2 0 The temperature T diverge as exp (s/d) ln 2 as d/s d/s = 3/4 { s/d } → d/s = 1 d/s = 3/2 0 and goes to zero as [e/(d/s)]/ 2πd/s as d/s 1, / (dNk 0.5 where e is the Euler-Napier number. ≫ V d/s = 2 d/s = 3
p C It is clear from expression (6) that d 2. -1 0 1 2 10 10 10 10 The nonmonototicity of the chemical potential, just re- ferred during the previous paragraphs, is revealed in the T / TF thermodynamic susceptibilities. In this work we focus ∂U FIG. 2. (Color online) Normalized thermodynamic suscep- on the specific heat at constant volume CV = and ∂T V tibilities as function of the dimensionless temperature T/T 1 ∂n F the isothermal compressibility κT = 2 , given by for different values of d/s, to say, 1/4, 1/2, 3/4, 1, 3/2, 2 n ∂µ T and 3. a) Specific heat per particle at constant generalized volume, the circles mark the corresponding values of CV at βµ C d d Li ( e ) TκT that is at the temperature at which κT has a maximum V = +1 d/s+1 − Nk s s Li ( eβµ) values, analogously, the triangles do the same for at Tµ where B d/s − µ has a maximum. b) Isothermal compressibility scaled with 2 d Li ( eβµ) d 2 π Γ(d/2) d/s −(d/s+1) d d/s κ0 = 2 sCs E , rhombus mark the corre- − s 2πd/ F βµ , (8a) − s Lid/s−1( e ) sponding values of κT at TCV that corresponds to the tem- − perature at which κT has a maximum values the triangles do the same as in a). Notice the nonmonotonic dependence on βµ T for d/s < 1. Lid/s−1( e ) κ = V − , (8b) T Nk T Li ( eβµ) B d/s − respectively. 0 ∗ TABLE I. The temperatures T , T , Tµ, TC and Tκ for In Fig. 2 the dimensionless CV (T ) s/dNkB (top panel) ∗ V T which µ(T 0) = 0, µ(T )= E , µ(T ) is maximum, C (T ) and κT /κ0 (bottom panel) are shown as function of the F µ V CV is maximum and κT (Tκ ) is maximum, for three character- dimensionless temperature T/TF for different values of T d/s, clearly, for d/s < 1, both quantities exhibit a non- istic values of d/s, namely 1/4, 1/2 and 3/4 at which a non- monotic behavior is observed. monotonous dependence on T . The specific heat clearly 0 ∗ exhibit the universal linear dependence on T in the low d/s T T Tµ TCV TκT temperature regime and rises with temperature evidenc- 0.25 15.6729 13.2260 5.0286 0.6532 0.3751 ing the effects of dimensionality. In the high tempera- 0.5 3.4797 1.8960 0.9365 0.8632 0.2906 ture regime all the curves converge to the classical re- 0.75 1.9830 0.6666 0.4086 1.3893 0.2080 sult dNkB/s. Analogously, the isothermal compressibil- ity exhibits the universal behavior in the low-temperature 5 regime namely a finite value due to the degeneracy pres- the energy change involved in the process nor the change sure. As temperature rises the effects of dimensionality in temperature, making the nature of the argument just are uncovered but are hidden again in the high temper- qualitative. In fact, the difficulty in quantifying those ature regime, where the classical dependence on temper- quantities arises from the use of the thermodynamic re- ature appears. lation The nonmonotonic dependence with temperature of ∂U both thermodynamic susceptibilities is manifested as a µ(S, ,N)= (10) V ∂N S,V global maximum at the temperatures TCV and TκT , re- spectively (see solid lines in both panels of Fig. 2). One which requires the knowledge of U(S, ,N), rarely con- would be tempted to propose that either of these temper- V atures would distinguish between two distinct behaviors sidered for analysis in the variables S, ,N. From (2) the functions N = N(µ, T, ), Eq. (3), andV S = S(µ, T, ), of the IFG: one where the corresponding susceptibility V V behaves anomalously and other where it behaves stan- Eq. (5), are obtained and solved in order to obtain dardly. Notice nevertheless, that such temperatures do U(S, ,N). In Fig. 3 the internal energy at constant en- tropyV is plotted as function of the particle density N/ not match between them nor with any of the tempera- V ∗ for S/kB =0.1 and for different values of the ratio d/s. ture Tµ or T , as can be quantitatively appreciated in V Table I and in Fig. 2, where solid triangles in both pan- The slope of the curves give the value of the chemical els identify the values of the corresponding susceptibility potential as given by expression (10). Also in the same Fig. 3, but in the bottom panel, the temperature of the evaluated at Tµ, solid circles in panel a) indicate the val- system, scaled with the Fermi temperature, as function ues of CV at TκT and analogously, solid squares in panel of the particle density is shown for S/kB =0.1. Clearly, b) mark the value of κT at TCV for d/s =1/4, 1/2, and V 3/4). Such discrepancy among all these temperatures the systems cools regardless of the ratio d/s, when adding makes difficult to consider them as points that mark the particles to the system in an isentropic way. It is possible to obtain an expression for µ(S, ,N) separation of two distinct thermodynamic behaviors. V For the signal value d/s = 1, expressions in terms of from (10) by the use of the asymptotic behavior of the σ Polylogarithm functions Liσ( z) ln(z) /Γ(σ +1)+ elementary functions are possible (black-dashed lines in 2 σ−2 − − ≃ Fig. 2), namely (π /6) ln(z) /Γ(σ 1)+..., after some algebra we have expressly that in the− degenerate regime βµ CV Li2(e ) −βµ βµ =2 βµ (1 + e ) ln(1 + e ) (9a) 2 2 NkB ln(1 + e ) − 3 S s s µ(S, ,N) EF 1+ 2 V 1 , eβµ V ≃ 2π kB N d d − " V # κT = V βµ βµ (9b) NkBT (1 + e ) ln(1 + e ) (11) where the nonmonotonic dependence on T is evident where we have used that the polylogarithm functions when d/s < 1. On the other hand, for the sake of com- of order 0, 1, 2 correspond to the elementary funtions pleteness we compute the system temperature as function Li0(z) = ln(1 + x), Li1(z) = z/1 + z and Li2(z) = of the particle density, in the degenerate limit, which is z xdx/1+ x respectively. For d/s > 1, the variation 0 given by with temperature of the thermodynamic susceptibilities Ris standard. 3 s S T (S, ,N) T V . (12) V ≃ π2 d F k N BV III. HEURISTIC EXPLANATION OF THE and, as is shown in the lower panel of Fig. 3, decrease as NONMONOTONIC DEPENDENCE OF µ ON T (N/ )−1 . FOR d/s < 1 HowV can we understand the rising of the chemical po- tential when d/s < 1? Consider the number of parti- The monotonic decreasing behavior of the chemical po- cles that can be excited by the energy kB T EF from tential with temperature for d/s 1 is understood from the d-dimensional Fermi sphere. This number≪ is approx- ≥ the argument based on the fact that the internal energy imately given by NkB T/EF while the number of avail- U, diminishes from its zero temperature value EF after able states above the Fermi energy can be approximated adiabatically adding a fermion at the small temperatures by g(EF )kBT . The quotient between both quantities is T TF . Quoting Cook and Dickerson [41], the system exactly s/d. This simple and heuristic argument shows cools≪ by redistributing the particles into the available en- that there are more single-particle excited states than ergy states in such a way that the particle added goes into excitable particles for d/s > 1, which is evident because “. . . a low lying, vacant single particle state, which will be of the monotonic increasing behavior of the DOS. In a little below EF ”. This is a consequence, as we will show principle all the excited particles can be accommodated below, that in the three-dimensional case the change of into the available states without violating Pauli’s princi- the Helmholtz free energy is dominated by the change of ple. The accommodation, however, is not arbitrary. The entropy in the low temperature limit, however, the argu- probability of occupation of the available states in ther- ment provided in [41], does not give neither the amount of mal equilibrium must follow the Fermi-Dirac distribution 6
F 0.8
E ν 4 S/kB = 0.1 d/s = 1/4
ν 10 d/s = 1/2 / d/s = 3/4
U d/s = 1 0.6 d/s = 3/2 2 d/s = 2 10 d/s = 3 ) T
( 0.4 R 0 d/s = 1/4 10 d/s = 1/2 0.2 d/s = 3/4 d/s = 1 d/s = 3/2 d/s = 2 -2 d/s = 3 10 -2 -1 0 10 10 10 0 -2 -1 0 1 8 10 10 10 10
F 10 d/s = 1/4 T S/k ν = 0.1 / B d/s = 1/2 T / T 6 d/s = 3/4 F T 10 d/s = 1 d/s = 3/2 FIG. 4. (Color online) The ratio of the number of particles in d/s = 2 4 10 d/s = 3 the energy interval [EF ,EF + ∆] to the number of available states in the same energy interval, R(T ), see Eq. (13), as function of temperature for different values of d/s. 2 10
0 10 same energy interval,
EF +∆ -2 10 dε g(ε)fF D(ε,T ) -2 -1 0 E 10 10 10 R(T )= F , (13) R EF +∆ N/ν dε g(ε) EF FIG. 3. (Color online) The dependence of the scaled inter- R nal energy U/VEF (top panel) and scaled temperature T/TF This quantity is shown in Fig. 4 as function of tem- (bottom panel) as function of the dimensionless particle den- perature with ratio ∆/EF = 0.001, for different values V V sity N/ , where denotes the systems volume scaled with of d/s. The choice ∆ kBT EF guarantees that a V ∼ ≪ an arbitrary volume 0. negligible number of particles occupy states out of the interval [EF , EF + ∆]. Under this condition we can ap- proximate R(T ) by fF D(EF ,T ) and for temperatures 1 π2 d T and therefore just a fraction of the excitable fermions are 0
≥ / n isotherm µ0(n) is a concave function of the particle den- sity, in other words, isotherms µ(T,n) of higher temper- d/s = 1/4 ature are situated below µ (n) (see Ref. [40] where the core d/s = 1/2
0 n 0.5 d/s = 3/4 case d/s = 3/2 is discussed). In contrast, for d/s < 1, d/s = 1 the zero temperature isotherm is a convex function of n d/s = 3/2 as shown in Fig. 5, and two possibilities may happen: i) d/s = 2 ′ ∗ ′ d/s = 3 if T < T the intersection occurs at ncore < n (dark- 0 broken lines in Fig. 5 for the case d/s =1/2); ii) on the 0 0.2 0.4 0.6 0.8 1 ′ ∗ ′ contrary, if T >T the intersection occurs at ncore >n 0 as shown explicitly in the same figure. The dependence T / T FIG. 6. (Color online) Fraction of particles in the Fermi- 0 F d/s = 1/2 sphere-like condensate as function of temperature scaled with 0 E 2 T for different values of d/s. µ/ isotherms T' > T* µ(T, ,N) from the thermodynamic relation 1 V ncore = n' ∂F µ(T, ,N)= (15) isotherms T' < T* V ∂N T,V ncore < n' n' 0 where F = F (T, ,N) stands for the Helmholtz free en- 0.5 1 1.5 ergy given by FV= k T ln Z (β) = U TS with 0 n/n B N,V T/TF=0.1 0 − − ZN,V(β)= E V exp βEN,V the canonical partition 1 3 3.5 N, {− } -1 function. The sum is made over the energies EN,V of all possible configurationsP with exactly N fermions in the FIG. 5. Isotherms in the plane µ-n for the case d/s = 1/2. volume . An advantage of expression (15) over the use 0 V Chemical potential and particle density are scaled with EF of the relation (10), is that at constant temperature and and n0, respectively, which correspond to the Fermi energy volume, the chemical potential measures the balance be- of an IFG with an arbitrary particle density n0. The thick tween the change of the internal energy and the heat line corresponds to the zero temperature isotherm, while thin- exchanged when the number of particles in the system is dashed lines label the isotherms with scaled temperatures 0 varied from N to N + 1, making it suitable for the use of T/TF = 0.1, 1, 3, and 3.5. For isotherms lying in the re- gion µ > 0 but below the zero temperature isotherm it is an energy-entropy argument [75]. Thus expression (15) provides a suitable operational definition, in the discrete possible to find ncore. case, of the chemical potential when only one particle is added isothermally to the system, namely [47, 73, 76] of the fraction ncore/n on temperature is shown in Fig. 6 for different values of d/s, and is explicitly given by the µ(T, ,N) =∆F F (N +1,T, ) F (N,T, ), (16a) expression V ≡ V − V ZN,V(β) d/s =kBT ln . (16b) n µ(T ) ZN+1,V (β) core = for T T 0. (14) n E ≤ F The rhs of expression (16a) can be explicitly written as ∆U T ∆S, where ∆U = U(T, ,N +1) U(T, ,N) and Notice that the nonmonotonic dependence of µ(T ) for ∆S −= S(T, ,N +1) S(T, ,NV), are the− internalV energy d/s < 1 makes n /n to reach the value 1 at the tem- core change of theV system− and theV heat produced T ∆S when perature T ∗ (see circles in Fig. 6). adding, isothermally, exactly one more fermion [77]. At zero temperature, the chemical potential is given by V. THE ARGUMENT ENERGY-ENTROPY the change in internal energy only, whose value coincides AND THE MEANING OF T ∗ with the Fermi energy of N + 1 fermions, i.e.
µ( ,N)= EF,N+1, (17) We now attempt to give a physical meaning to µ in V the region where is larger than EF , i.e. in the inter- where EF,N denotes the explicit dependence of the Fermi val of temperatures [0,T ∗]. For this purpose we compute energy on the particle number. If this value is subtracted 8 from (16a) we have that TABLE II. The firts 4 expressions for the canonical partition ∆µ = ∆U ′ T ∆S (18) function are shown − where ∆µ = µ(T, ,N) E and ∆U ′ = ∆U N ZN (β) V − F,N+1 − EF,N+1. In this way, if for a given temperature we have 1 Z1(β) 1 2 − 1 that ∆µ 0, i.e., the chemical potential lies below the 2 2 Z1 (β) 2 Z1(2β) Fermi energy,≤ then the relative change in the internal 1 3 − 1 1 3 6 Z1 (β) 2 Z1(β)Z1(2β)+ 3 Z1(3β) energy is smaller than the respective heat exchange by 1 4 − 1 2 1 4 24 Z1 (β) 4 Z1 (β)Z1(2β)+ 3 Z1(β)Z1(3β)+ adding the particle. In other words, the effects of the 1 2 1 Z1 (2β) − Z1(4β) addition of a particle to the system, in an isothermal 8 4 way, are such that the entropic effects dominate over the energetic ones at that T . This argument accounts for the monotonic-decreasing behavior of µ with T , and is and only the algorithm for computing the unrestricted equivalent with argument given in Ref. [41]. Further, partitions of the integer N is needed. Such algorithm if ∆µ > 0 for a given T , then the energetic changes are forms part of the MATHEMATICA software package dis- the ones that dominate over the entropic ones, which give tribution. The computation time and memory require- origin to the rise of µ above the Fermi energy as has been ments grow with number of partitions p(N) of N, which grows asymptotically as exp √n thus limiting compu- shown in the previous section. The temperature that sep- { } ∗ tation to N 10. Suprisingly, the calculation exhibits a arate both regimes coincides with T , which is different ∼ from zero when d/s < 1. This suggest the possibility fast convergence to the well-known result obtained from of interpreting T ∗ as a critical temperature at which a (3) for 64 particles (see Fig. 7 for the box-like trap). phase transition occurs. For the box potential in one dimension, the single- In order to show the validity of these ideas we first particle partition function is given in terms of the Jacobi ∞ n2 use expression (16b) to compute µ in two distinct theta function ϑ3(u, q)=1+2 n=1 q cos(2n u) as 1 −βε0 ∗ one-dimensional systems each consisting of N spinless Z1(β) = 2 ϑ3(0,e ) 1 where the energy scale ε0 − P 2 2 fermions. One corresponds to the IFG trapped by a box- ~ π in the argument of the exponential is 2mL2 . For few par- like potential (d/s = 1/2), the other to the experimen- 2 ticles, the chemical potential does not rise as (T/TF ) tally feasible system of and IFG trapped by a harmonic as is expected from the grand-canonical-ensemble result, trap (d/s = 1). We show that for the former case, µ rises but it grows much more slower as is shown in the inset above EF,N+1 and eventually return to its decreasing be- of Fig. 7. This is consequence of the low-temperature havior as the system temperature is increased from zero. behavior of the partition function, which satisfies that For the later, we show that the µ < EF,N+1 for all T > 0. βEF,N+1 −β(2N+1)ε0 ZN /ZN+1 e 1+ e for T/TF 1, For exactly N non-interacting fermions, the partition ∼ ≪ leading thus to ∆µ k Te−β(2N+1)ε0 . function satisfy the recursive relation [78, 79] B For the harmonic∼ potential in dimension one, an exact N analytical expression for ZN (β) is known [24, 81], namely 1 n+1 ZN (β)= ( 1) Z1(nβ)ZN−n(β), (19) N − N n=1 ~ 2 ω −1 X ZN (β) = exp N β [1 exp( β~ω j)] . where Z (β)= exp βε is the single-particle parti- − 2 − − 1 ε j=1 tion function with ε the{− single-particle} energy spectrum Y P (21) and Z0(β) 1. A direct application of Eq. (16b) leads to Expression≡ (19) can be reduced to the calculation of
Z1(mβ), with m a positive integer, by noting that ZN (β) µ = EF,N+1+kBT ln [1 exp( βEF,N+1) exp( β~ω/2)] . can be written as a sum of the product over the distinct − − − (22) parts of all the partitions (λ1, λ2,...,λr) of N (a par- Clearly (22) is a monotonically decreasing function of T tition is defined as a nonincreasing{ sequence} of positive r agreeing with the grand canonical ensemble result, ex- −T /T integers λ1, λ2,...,λr such that i=1 ηiλi = N, where pression µ = E + k T ln 1 e F is recovered in F B − ηi denotes the multiplicity of the part λi in a given par- the limit N . tition (for instance the partitionP 2 + 2 + 2 + 1 of 7, the → ∞ part 2 has multiplicity 3, see Ref. [80] pp.1), thus r VI. CONCLUSIONS AND FINAL REMARKS ηm N ( 1) ηm ZN (β) = ( 1) −ηm [Z1(λmβ)] , − λm ηm! {(λ1,λ2,...,λr )} m=1 In this paper we have presented a discussion on the X Y (20) meaning of the nonmonotonic dependence on temper- The first four terms can be checked straightforwardly and ature of the thermodynamics properties of low dimen- are shown in Table II. sional, trapped, IFGs, with focus on the chemical po- Computationally, evaluation of expression (20) is faster tential (a similar behavior has been predicted for weakly than evaluating expression (19) since recursion is avoided repulsively interacting bose gases [82] in that a hard core 9
[28, 40], and introduced a new characteristic temperature T ∗ T 0, which corresponds to the nonzero value of the ≤ ∗ temperature at which µ(T , ,N)= EF . 1 We found that T ∗ marksV the temperature at which the particle density of a Fermi-like core ncore, that starts forming T 0, saturates at the value of the total particle F,N+1 ∗ 1.03 density of the system n. This suggest that T can be
E Sommerfeld Approx. N = 2 considered as the relevant temperature of the isotropi- / 0.5 1.02 N = 4
µ N = 8 cally trapped IFG, as is supported by the energy-entropy- 1.01 N = 16 N = 32 like argument presented in Sect. V. The region in the ∗ 1 from Eq. (3) µ-T plane, for which µ > E for T T , represents F ≤ 0 0.1 0.2 the set of thermodynamic states for which the change in 0 0 1 2 3 the Helmholtz free energy, when increasing the particle density of the system, is dominated by the changes of T / T the internal energy and would correspond to an ordered F phase. In the complementary region for which T T ∗, the thermodynamic states are characterized by changes≥ FIG. 7. (Color online) The chemical potential, Eq. (16b), in F (T, ,N) dominated by heat exchange by changing scaled with E as function of the dimensionless tempera- F,N+1 entropy,V and can be considered as a “disordered phase”. ture T/TF , with TF = EF,N+1/kB , for the number of particles N = 2, 4, 8, 16, and 32. Note that the results with 32 parti- Though, heuristic energy-entropy arguments have been cles (double-dashed-dotted line) are close to the grand canoni- used to uncover the possibility of a phase transition [75], cal ensemble result (black line in long dashes) computed from we want to emphasize that we are not claiming the ex- (3). A comparison with the Sommerfeld approximation for istence of a phase transition in the IFG, on the basis low temperatures is shown in the inset. that thermodynamic quantities do not show a singular behavior of the thermodynamic susceptibilities at T ∗. Though the chemical potential is not directly measured in current experiments, development on imaging tech- Bose gas behaves, at least qualitatively, as an ideal Fermi niques of ultracold gases [83–86] has open the possibility gas). The parameter used to characterize the trapping to experimentalists to measure the local particle-density and dimensionality d of the system is merely d/s, that in situ and from the data to extract µ and T . explicitly appears in the single-particle density of states (1). Thus low dimensional trapped systems are charac- terized by values of d/s < 1. In this range of values, ACKNOWLEDGMENTS the chemical potential, the specific heat at constant vol- ume and the isothermal compressibility, exhibit a non- I want to express my gratitude to Mauricio Fortes and monotonic dependence on temperature which have been Miguel Angel Sol´ısfor valuable discussions, and to Omar characterized by the temperatures Tµ, TCV , TκT [60] re- Pi˜na P´erez for helping in generating figure 2. Support is spectively. We also have computed T 0 as function of d/s acknowledge to UNAM-DGAPA PAPIIT-IN113114.
[1] W. Pauli, Zeitschrift f¨ur Physik 31, 765-83 (1925). Gas,” Phys. Rev. Lett. 88, 120405 (2002). [2] E. Fermi, Rend. Acc. Lincei 3, 145 (1926). [10] Z. Hadzibabic, S. Gupta, C. A. Stan, C. H. Schunck, M. [3] P. Dirac, Proc. R. Soc. A 112, 661 (1926). W. Zwierlein, K. Dieckmann, and W. Ketterle, “Fifty- [4] J. S´oyom, “The Fermi gas model of one-dimensional con- fold Improvement in the Number of Quantum Degenerate ductors,” Adv. Phys. 28 (2), 201-303 (1979). Fermionic Atoms,” Phys. Rev. Lett. 91, 160401 (2003). [5] J. Voit, “One-dimensional Fermi liquids,” Rep.Prog. [11] T. Fukuhara, Y. Takasu, M. Kumakura, and Y. Taka- Phys. 58 (9), 977 (1995). hashi, “Degenerate Fermi Gases of Ytterbium,” Phys. [6] X.-W. Guan, M. T. Batchelor, and C. Lee, “Fermi gases Rev. Lett. 98, 030401 (2007). in one dimension: From Bethe ansatz to experiments,” [12] C. Chin, M. Bartenstein, A. Altmeyer, S. Riedl, S. Rev. Mod. Phys. 85 (4), 1633–1691 (2013). Jochim, J. H. Denschlag, and R. Grimm, “Observation [7] B. DeMarco, and D. Jin, “Onset of Fermi Degeneracy in of the Pairing Gap in a Strongly Interacting Fermi Gas,” a Trapped Atomic Gas,” Science 285, 1703-1706 (1999). Science 305, 1128-1130 (2004). [8] A.G. Truscott, K.E. Strecker, W.I. McAlexander, G.B. [13] G. K. Campbell, M. M. Boyd, J. W. Thomsen, M. J. Partridge, R.G. Hule, “Observation of Fermi Pressure Martin, S. Blatt, M. D. Swallows, T. L. Nicholson, T. in a Gas of Trapped Atoms,” Science 291, 2570-2572 Fortier, C. W. Oates, S. A. Diddams, et al., “Probing (2001). Interactions Between Ultracold Fermions,” Science 324, [9] S. R. Granade, M. E. Gehm, K. M. O’Hara, and J. E. 360-363 (2009). Thomas, “All-Optical Production of a Degenerate Fermi [14] M.W. Zwierlein, C.H. Schunck, A. Schirotzek, W. Ket- 10
terle, “Direct Observation of the Superfluid Phase Tran- [34] Mir Mehedi Faruka, and G.M. Bhuiyanb, “Thermody- sition in Ultracold Fermi Gases,” Nature 442, 54-58 namics of Ideal Fermi Gas Under Generic Power Law (2006). Potential in d-dimensions,” A. Phys. Pol. B 46, No. 12, [15] Y. Shin, A. Schirotzek, C.H. Schunck, and W. Ketterle, 2419, (2015). “Realization of a strongly interacting Bose-Fermi mixture [35] B. DeMarco, S.B. Papp, and D.S. Jin, “Pauli Blocking of from a two-component Fermi gas,” Phys. Rev. Lett. 101, Collisions in a Quantum Degenerate Atomic Fermi Gas,” 070404 (2008). Phys. Rev. Lett. 86, 5409-5412 (2001). [16] D.A. Butts, and D.S. Rokhsar, “Trapped Fermi Gases,” [36] D.M. Stamper-Kurn, “Shifting entropy elsewhere,” Phys. Rev. A 55, 4346 (1997). Physics 2, 80 (2009). [17] J. Schneider and H. Wallis, “Mesoscopic Fermi gas in a [37] J. Catani, G. Barontini, G. Lamporesi, F. Rabatti, G. harmonic trap,” Phys. Rev. A 57, 1253-1259 (1998). Thalhammer, F. Minardi, S. Stringari and M. Inguscio, [18] M. Li, Z. Yan, J. Chen, L. Chen, and C. Chen, “Thermo- “Entropy Exchange in a Mixture of Ultracold Atoms,” dynamic properties of an ideal Fermi gas in an external Phys. Rev. Lett. 103, 140401 (2009). potential with U = brt in any dimensional space,” Phys. [38] J.-S. Bernier, C. Kollath, A. Georges, L. De Leo, and Rev. A 58, 1445-1449 (1998). F. Gerbier, “Cooling fermionic atoms in optical lattices [19] A. Thislagam, “Dimensionality dependence of Pauli by shaping the confinement,” Phys. Rev. A 79, 061601 blocking effects in semiconductor quantum wells,”J. (2009). Phys. Chem. Solids 60, 497-502 (1999). [39] T. Paiva, R. Scalettar, M. Randeria and N. Trivedi, [20] P. Vignolo, A. Minguzzi, and M. P. Tosi, “Exact Par- “Fermions in 2D Optical Lattices: Temperature and En- ticle and Kinetic-Energy Densities for One-Dimensional tropy Scales for Observing Antiferromagnetism and Su- Confined Gases of Noninteracting Fermions,” Phys. Rev. perfluidity,” Phys. Rev. Lett. 104, 066406 (2010). Lett. 85, 2850 (2000). [40] V. Romero-Rochin, “Phase diagram of quantum fluids. [21] M. Brack, and B.P. van Zyl, “Simple Analytical Particle The role of the chemical potential and the phenomenon and Kinetic Energy Densities for a Dilute Fermionic Gas of condensation.”arXiv:1109.1815[cond-mat.quant-gas]. in a d-Dimensional Harmonic Trap, ”Phys. Rev. Lett. 86, [41] G. Cook & R.H. Dickerson, “Understanding the Chemi- 1574-1577 (2000). cal Potential,” Am. J. Phys. 63 (8), 737-742 (1995). [22] G.M. Bruun and C.W. Clark, “Ideal gases in time- [42] R. Baierlein, “The elusive chemical potential,” Am. J. dependent traps,” Phys. Rev. A 61, 061601 (2000). Phys. 69 (4), 423-434 (2001). [23] F. Gleisber, W. Wonneberger, U. Schl¨oder and C. Zim- [43] G. Job and F. Herrmann, “Chemical potential a quan- mermann, “Noninteracting fermions in a one-dimensional tity in search of recognition,” Eur. J. Phys. 27, 353-371 harmonic atom trap: Exact one-particle properties at (2006). zero temperature,” Phys. Rev. A 62, 063602 (2000). [44] Carl E. Mungan, ”Chemical potential of one-dimensional [24] M.N. Tran, M.V. Murthy, and R.K. Bhaduri, “Ground- simple harmonic oscillators,” Eur. J. Phys. 30, 1131-1136 state fluctuations in finite Fermi systems,” Phys. Rev. E (2009). 63, 031105 (2001). [45] M.R.A. Shegelski, “New result for the chemical potential [25] Z. Akdeniz, P. Vignolo, A. Minguzzi, and M.P. Tosi, of intrinsic semiconductors: Low-temperature breakdown “Temperature dependence of density profiles for a cloud of the Fermi-Dirac distribution function,” Solid State of noninteracting fermions moving inside a harmonic trap Commun. 58, 351-354 (1986); in one dimension,” Phys. Rev. A 66, 055601 (2002) [46] P.T. Landsberg and D.C. Browne, “The Chemical Po- [26] P. Vignolo and A. Minguzzi, “Shell structure in the den- tential of an Intrinsic Semiconductor near T=0,” Solid sity profiles for noninteracting fermions in anisotropic State Commun. 62, 207-208 (1987); harmonic confinement,” Phys. Rev. A 67, 053601 (2003). [47] M.R.A. Shegelski, “The chemical potential of an ideal in- 1332 (1979). trinsic semiconductor,” Am. J. Phys. 72, 676-678 (2004); [27] M. Grether, M. Fortes, M. de Llano, J.L. del R´ıo, F.J. [48] T.A. Kaplan, “The Chemical Potential,” J. Stat. Phys. Sevilla, M.A. Sol´ıs and A.A. Valladares, “Harmonically 122, 1237-1260 (2006). trapped ideal quantum gases,” Eur. Phys. J. D 23, 117- [49] S. Panda and B. K. Panda, “Chemical potential and 124 (2003). internal energy of the noninteracting Fermi gas in [28] D-V Anghel, “Condensation in ideal Fermi gases,” J. fractional-dimensional space,” Pramana — Journal of Phys. A: Math. Gen. 36, L577-L783 (2003). Physics, 75, 393-402 (2010). [29] M.N. Tran, “Exact ground-state number fluctuations of [50] F.J. Sevilla and L. Olivares-Quiroz, “Chemical potential trapped ideal and interacting fermions,” J. Phys. A: for the interacting classical gas and the ideal quantum Math. Gen. 36, 961 (2003). gas obeying a generalized exclusion principle’,’ Eur. J. [30] B.P. van Zyl, R.K. Bhaduri, A. Suzuki and M. Brack, Phys., 33 (3), 709 (2012). “Some exact results for a trapped quantum gas at finite [51] F. J. Sevilla and O. Pi˜na, arXiv:1407.7187 [cond- temperature,” Phys. Rev. A 67, 023609 (2003). mat.stat-mech]. [31] E.J. Mueller, “Density profile of a Harmonically Trapped [52] P. Salas and M. A. Sols, ”Trapping Effect of Periodic Ideal Fermi Gas in Arbitrary Dimension,” Phys. Rev. Structures on the Thermodynamic Properties of a Fermi Lett. 93, 190404 (2004) Gas,“ Journal of Low Temperature Physics 175, 427-434 [32] D. Anghel, O. Fefelov, and Y.M. Galperin, “Fluctuations (2014). of the Fermi condensate in ideal gases,” J. Phys. A: Math. [53] P. C. Hohenberg, “Existence of long-range order in Gen 38, 9405 (2005). one and two dimensions,” Phys. Rev. 158 (2), 383-390 [33] Dae-Yup Song, “Exact coherent states of a noninteract- (1967). ing Fermi gas in a harmonic trap,” Phys. Rev. A 74, [54] M.H. Lee, “Equivalence of ideal gases in two dimensions 051602(R) (2006). and Landens relations,” Phys. Rev. E 55, 1518 (1997). 11
[55] R. May, “Quantum statistics of ideal gases in two dimen- [74] T.P. Meyrath, F. Schreck, J.L. Hanssen, C.-S. Chuu, and sions,” Phys. Rev. 135, A1515 (1964). M.G. Raizen, “Bose-Einstein condensate in a box,” Phys. [56] D-V Anghel, “Gases in two dimensions: universal ther- Rev. A 71, 041604(R) (2005). modynamics and its consequences,” J. Phys. A: Math. [75] B. Simon and A. D. Sokal, “Rigorous Entropy-Energy Gen. 35, 7255-7267 (2002). Arguments,” J. Stat. Phys. 25 (4), 679-694 (1981). [57] K. Sch¨onhammer and V. Meden, “Fermion-Boson trans- [76] J Tobochnik, H Gould and J Matcha, “Understanding mutation and comparison of statistical ensembles in one temperature and chemical potential using computer sim- dimension,” Am. J. Phys. 64, 1168-1176 (1996). ulations,” Am. J. Phys. 73, 708-716 (2005). [58] M. Crescimanno and A.S. Landsberg, “Spectral equiv- [77] Expression (16a) arises from the forward difference dis- alence of bosons and fermions in one-dimensional har- cretization of the conventional definition (15). The phys- monic potentials,” Phys. Rev A 63, 035601 (2001). ical interpretation is simple and the same for the back- [59] R.K. Pathria, “Similarities and differences between Bose ward difference: they simply give the change in the free and Fermi gases,” Phys. Rev. E 57, 2697-2702 (1998). energy when adding an extra particle to the system. Ex- [60] M. Grether, M. de Llano, and M.A. Sol´ıs, “Anomalous perimental realization of such situation are for instance, behavior of ideal Fermi gas below two dimensions,” Eur. Bose-Einstein condensation at constant temperature by Phys. J. D 25, 287-291 (2003). increasing the particle number of the system [? ] and [61] E. Cetina, F. Maga˜na and A.A. Valladares, “The free- injection of electrons into low-dimensional systems. The electron gas in n dimensions,” Am. J. Phys. 45 (10), negative of such difference gives the change in the free en- 960-963 (1977). ergy when removing a particle from the system. Another [62] M. Howard Lee, “Polylogarithmic analysis of chemical definition of µ is given by the average of the forward and potential and fluctuations in a D-dimensional free Fermi backward differences which is suitable when both pro- gas at low temperatures,” J. Math. Phys. 36 (3), 1217- cesses, adding and removing a particle, are present at 1229 (1995). finite temperatures as occurs in the grand canonical en- [63] M. Grether, M. de Llano, and M. Howard Lee, “Anoma- semble. lous behavior of ideal Fermi gas below 2D: The ideal [78] P. Borrmann and G. Franke, “Recursion Formulas for quantum dot and the Pauli exclusion principle,” Int. J. quantum statistical partition functions,” J. Chem. Phys. Mod. Phys. B 23, 4121-4128 (2009). 98, 2484-2485 (1993). [64] I. Ch´avez, M. Grether and M. de Llano,“Low- [79] S. Pratt, “Canonical and Microcanonical Calculations for dimensional Fermi and Bose gases”, Physica E 44, 394- Fermi Systems,” Phys. Rev. Lett. 84, 4255-4258 (2000). 399 (2011). [80] George E. Andrews, Encyclopedia of Mathematics and its [65] Indeed, the precise argument presented in Ref. 41 is only Applications, G.-C. and P. Tur´an Eds., Vol. 2 The theory valid for the free IFG in three dimensions. of partitions (Addison-Wesley 1976). [66] V.C. aguilera-Navarro, M. de Llano and M.A. Sol´ıs, [81] M.N. Tran, M.V.N. Murthy, and R.K. Bhaduri, “On the “BoseEinstein condensation for general dispersion rela- quantum density of states and partitioning an integer,” tions,” Eur. J. Phys. 20, 177-182 (1999). Ann. Phys. (N.Y.) 311, 204-219 (2004) and references [67] V. Bagnato, D.E. Pritchard, and D. Kleppner, “Bose- therein. Einstein condensation in an external potential,” Phys. [82] Z. Zhang, H. Fu, J. Chen, and M. Li, “Chemical potential Rev. A 35, 4354-4358 (1987). of a trapped interacting Bose gas,” Phys. Rev. A 79, [68] R.K. Pathria, Statistical Mechanics (2nd ed. Oxford 055603 (2009). 1996). [83] T.-L. Ho and Q. Zhou, “Obtaining the phase diagram [69] K. Huang, Statistical Mechanics (2nd ed. John Wiley & and thermodynamic quantities of bulk systems from den- Sons 1987). sities of trapped gases,” Nature Physics 6, FEB 131-134 [70] N. Sandoval-Figueroa and V. Romero-Rochin, Phys. Rev. (2011). E 78, 061129 (2008) and references therein. [84] T. M¨uller, B. Zimmermann, J. Meineke, J.-P. Bruntut, [71] In difference with a gas in a 3D-box potential where ra- T. Esslinger, and H. Moritz, “Local observation of anti- tio N/V is well defined as V being the box volume is bunching in a trapped Fermi gas,” Phys. Rev. Lett. 105 well defined too, in the gas trapped by harmonic po- (4), 040401 (2010). tentials the thermodynamic limit is defined through the [85] S. Nascimb`ene, N. Navon, K.J. Jiang, F. Chevy, and C. limit of ω → 0 and N → ∞ such that the product ω3N Salomon, “Exploring the thermodynamics of a universal is constant as has been discussed in the literature see Fermi gas,” Nature 463, FEB 1057-1060 (2010). for instance S.R. de Groot, G.J. Hooyman and C.A. Ten [86] C. Sanner, E. J. Su, A. Keshet, R. Gommers, Y.-I. Seldam, Proc. R. Soc. London Ser. A 203 266 (1950). Shin, W. Huang, and W. Ketterle, “Suppression of Den- [72] L. Lewin, Dilogarithms and Associated Functions, (Mc- sity Fluctuations in a Quantum Degenerate Fermi Gas,” Donald, London, 1958). Phys. Rev. Lett. 105 (4), 040402 (2010). [73] N.W. Ashcroft and N.D. Mermin, Solid State Physics, (Saunders College 1976).