Arxiv:Cond-Mat/9705180V2 [Cond-Mat.Stat-Mech] 18 Jun 1997 Ehd Frarve E 5

Home , Superdiamagnetism

arXiv:cond-mat/9705180v2 [cond-mat.stat-mech] 18 Jun 1997 ehd Frarve e 5. hsis This [5].) see action review effective a the (For using method. found as potential modynamic of features general systems. illustrate dimensions magnetized to spatial in used of transitions be study phase can the 3 Thus than other gas. Bose free the for that dimensions feature spatial general other These for temperatures. high persist very of “dimen- at one-dimensio values a except a undergoes large as gas effectively gas For behaves the spe- and that the reduction” curve. show sional for we field curves field free the magnetic the the reduced approach is heat field s we cific magnetic phase field, the magnetic a the as of de- of how absence values to all the used for Although be persists can i transition temperature. which there critical heat transition, a specific phase the fine a in there of maximum although sense a di- that the still show in three we condensation the free. no and is for otherwise clarified, is is condensation case which Bose-Einstein mensional but of field notion magnetic The constant a with oegsi eea pta dimension spatial general magnetized the (This a at in looked to gas state. [4] Subsequent May Bose ground [2,3].) Schafroth, Ref. of the work in in there the emphasized particles that been imply has of point not build-up infinite does no the this (in is However transition phase the limit). no for is volume situation there the condensation gas, Unlike what Bose Bose-Einstein gas. free emphasize of Bose to magnetized absence important the the is for by It temperature meant gas. condensation is Bose the free to related the scale of a which at effect in Meissner-Ochsenfeld set magne- the exhibited the gas that tized gas. demonstrated free Schafroth the by Furthermore, exhibited famil- the condensation destroyed were Bose-Einstein field the iar magnetic field that constant showed a magnetic He of constant presence [1]. Schafroth a by evaluated in correctly gas Bose relativistic effect. Meissner-Ochsenfeld for ∗ † -aladdress: E-mail -aladdress: E-mail ewl ei ihteepeso o h oa ther- total the for expression the with begin will We h hroyaia rpriso hre non- charged a of properties thermodynamical The esuytecagdnnrltvsi oegsinteracting gas Bose non-relativistic charged the study We D .Fnlyw xmn h antzto n the and magnetization the examine we Finally 5. = D ≥ hr a hs rniinlk htin that like transition phase a was there 5 oeEnti odnaino h antzdielBs ga Bose ideal magnetized the of condensation Bose-Einstein [email protected] [email protected] eateto hsc,Uiest fNwateUo Tyne, Upon Newcastle of University Physics, of Department D ecsl pnTn,Uie igo E 7RU NE1 Kingdom United Tyne, Upon Newcastle n eso results show we and u .Standen B. Guy D n showed and Otbr9 2018) 9, (October how nal s s ∗ n ai .Toms J. David and 1 tr,adi sti eaiu hc a iers oa The to rise ( transition. give condensate phase a can of a which case behaviour signalling this temperature is critical it and ature, uthave must h aueo h ouin for solutions the of nature The then D htwe antcfil speettecniinfor condition ( the gas free present the is field magnetic Ψ know a we when [4,7] work that previous From distinguished. be must h rttr nΩi h cascl atwihaccounts by which described part condensate “classical” possible the a is for Ω in term first The fi sitrrtda hs rniinwt non-zero a with transition phase a as condensation interpreted Bose-Einstein is destroys it field if magnetic a of hsi bandb xrmzn 1 ihrsetto respect with (1) extremizing by obtained is This emi h tnadepeso o h thermodynamic Here the mechanics. for statistical expression from standard potential the is term fteeeg eesfracagdpril naconstant a in particle in charged knowledge Although a a [6]. from for field evaluated magnetic levels be energy can the (1) of in Ω system. the steceia oeta,and potential, chemical the is ses oso that show to easy is ota h rttr nΩmksn otiuint the to for contribution solve no can makes we Ω if in and term charge, first the that so edn opnn,w ilol osdrtepresence here. the component consider inde- single one only a than will of more we have component, general pendent in may field magnetic SeRf 3 o eea ramn. fi spsil for possible is it If treatment.) general a µ for [3] Ref. (See eB/m sdtrie yteslto to solution the by determined is oetenergy. lowest where h oeteeg ievle o h aeo constant a of case strength the of by For field magnetic determined eigenvalue. potential energy chemical lowest the the for value critical a oa charge total ¯ orahteciia value critical the reach to ntepyial neetn aeof case interesting physically the In hte rntteei o-eocondensate non-zero a is there not or Whether hte rnot or Whether ≥ is 0 6= o o-eocondensate. non-zero a for 3 sncsaiynnzr n ehv condensate. a have we and non-zero necessarily is Ψ ¯ nti aeteei onntiilcondensate. non-trivial no is there case this In . f = Ω 0 ( D x µ < µ steegnucincrepnigt the to corresponding eigenfunction the is ) ≥ B Z Q + )wihrqie h pta dimension spatial the requires which 0) = − = .Ti otat ihtecniinfor condition the with contrasts This 5. d † c D β 1 oee fi sntpsil od this do to possible not is it if However . − 2 x 1 m µ X if 0 = Ψ ∂ ∂µ ¯ n Ω D a reach can Ψ 2 ¯ fw e h condensate the set we If . ln 1 m 2 )adn odnae( condensate no and 0) 6= Ψ ¯ h | B D 1 − ehave we − Ψ µ < µ ¯ µ eµ | c E e 2 hnw have we then , − µ 0 = Ψ D n ¯ − µ µ β c ntrsof terms in r h nrylvl of levels energy the are ( c pta iesosthe dimensions spatial eµ E hne ihtemper- with changes sdtrie ythe by determined is where n E | − D Ψ 0 ¯ s . eµ ,adtesecond the and Ψ, ¯ | = 2 h presence the 3 = ) i ω/ eµ β . ( = c where 2 Q = Ψ ¯ hnwe then , kT E ∝ 0) = Ψ 0, = Ψ ¯ ¯ 0 Ψ ) ¯ f .It Ψ. − gives ¯ 0 ω 1 0 6= ( (1) (2) , x = µ ) condensate [1]. The presence of a magnetic field, no mat- the curve starts to resemble that for the free Bose gas. ter how small, results in the charged Bose gas behaving in The maximum in the specific heat becomes sharper as B a different manner from the same system when no mag- is decreased with the temperature of the maximum ap- netic field is present. However on physical grounds one proaching the temperature T0 for the free Bose gas. Sim- might expect that for a very small magnetic field the sys- ilar results hold if the number of particles in the ground tem should behave in almost the same way as it would if state is computed. Thus although for B 6= 0 the sys- B = 0. This is one of the issues addressed in the present tem does not have a phase transition characterized by a paper. In addition the approach of Ref. [3] is clarified critical temperature and a non-zero condensate Ψ¯ and is with a particular example. therefore very different from the same system with B = 0, For the free Bose gas a key feature of Bose-Einstein quantities of direct physical interest, such as the specific condensation is the behaviour of the specific heat at the heat or ground state particle number, are closely approx- critical temperature T0. For D = 3 the specific heat has imated by the free Bose gas result for small enough values a maximum at T = T0, and is continuous there with a of B. We will obtain approximate analytical expressions discontinuous first derivative [8]. For D = 4 the specific for small B later. heat and its first derivative are continuous at T = T0, and

for D ≥ 5 the speciﬁc heat is discontinuous at T = T0 [9].

free

In all cases the non-smooth behaviour of the speciﬁc heat corresponds to a maximum. As already stated the pres-

k N

ence of a magnetic ﬁeld suppresses the phase transition if

D =3, 4. For D ≥ 5 the speciﬁc heat for the magnetized

gas would be expected to show a non-smooth behaviour

at a critical temperature Tc corresponding to the exis-

tence of a phase transition; however for D = 3, 4 the

speciﬁc heat should be perfectly smooth since no phase sp ecic heat in units of

transition occurs. We will concentrate initially on D = 3.

By ﬁrst calculating the internal energy deﬁned by U =

∂Ω/(∂β) with µβ,V,B held ﬁxed, and then the speciﬁc

heat at constant volume by CV = ∂U/(∂T ) with V,B

x in units of x and the total charge Q held ﬁxed, it is possible to show 0 FIG. 1. This shows the speciﬁc heat for the magnetized that Bose gas in units of kN. The curve labelled free is the result

CV 1 3 for the uncharged Bose gas. The x coordinate in the figure is = xψ(1, 1, 2)+ ψ(0, 0, 1) ~ ~ kN ψ(1, 0, 1) 4 x = ω/(kT ) in units of x0 = ω/(kT0) with T0 the critical temperature for the free Bose gas. The numbers labelling the ψ2(2, 1, 2) +x2(ψ(2, 1, 2)+2ψ(2, 2, 3)) − x2 curves are the values for x0. The smaller numbers indicate ψ(2, 0, 1) lower values for the magnetic field. 1 ψ(1, 0, 1) ψ(2, 1, 2) − − x , (3) 4 ψ(2, 0, 1) ψ(2, 0, 1) It is also of interest to observe what happens as the value of the magnetic field is increased. It is clear from where the curves in fig. 1 with larger values of x0 that with increasing B the peak in CV shifts to smaller values of Q m 3/2 x (which corresponds to higher temperatures) and de- N = = V xψ(1, 0, 1) , (4) e 2πβ creases in magnitude. (Although the resolution of the graphs do not show it, there is always a maximum in the is the total number of particles present. Here ψ(a,b,c) specific heat.) As x approaches 0, corresponding to in- denotes a class of sums defined by finitely high temperatures, the specific heat approaches the classical Maxwell-Boltzmann result of 1.5 regardless ∞ a−D/2 −l(ǫ+b)x −lx −c of the magnitude of B, as would be expected. For large B ψ(a,b,c)= l e (1 − e ) . (5) the specific heat actually resembles that for a free Bose

Xl=1 gas in one spatial dimension. This can be understood (We take D = 3 in eqs.(3) and (4).) heuristically by the fact that classically the orbit of a It is possible to solve (4) for ǫ where µ = ω(1/2 − ǫ) charged particle in a constant magnetic field is a spiral around the direction specified by the magnetic field. The and then to calculate CV . The result of this is shown in fig. 1. For comparison we have also indicated the re- size of the orbit is inversely proportional to the magnetic sult for the free Bose gas. As expected the specific heat field; thus as B is increased the motion of the particle be- is perfectly continuous for B 6= 0. As B is decreased comes more and more one-dimensional. This substanti- ates the general point of view used in Ref. [3] in which the (following the curves with smaller values of x0) although leading behaviour of thermodynamic quantities is well CV remains continuous the peak starts to increase and

2 2ζ(1/2) 1 approximated by studying the lowest energy solution. In ǫ ≃ a + − a x1/2 this case as B is increased the gap between the ground π1/2ζ (3/2,a) 2 0 H and first excited state becomes increasingly large and the 1/2 6ζ(3/2) x ground state makes the leading contribution. + 1/2 1 − , (8) 1/2 x0 Although for B 6= 0 with D = 3 there is no phase π ζH (3/2,a)x0 " # transition, Schafroth [1] showed that the magnetized gas 1/2 1/2 exhibited the Meissner-Ochsenfeld effect. He found which is good if 1 − (x/x0) << x0 . Here a is just a 3/2 constant defined by ζH (1/2,a) = 0 with ζH the Hurwitz ρ x0 M ≃− 1 − (6) ζ-function. (Numerically we find a ≃ 0.302721829.) As 2m x x0 is decreased the range over which the expansion holds where ρ is the charge density and x0 = ω/(kT0). The ap- becomes increasingly small, but it is always good at x = proximation assumes that T ≤ T0. Using M = ∂Ω/(∂B), x0 (corresponding to T = T0). The agreement between the exact expression may be computed to be the approximation (8) and a numerical evaluation of ǫ is very good for small x0. ρ 2ψ(0, 0, 1) − 2xψ(1, 1, 2) M = − 1 − . (7)

2m xψ(1, 0, 1)

free

This expression holds for any temperature.

k N

sp ecic heat in units of

ratio

x in units of x

FIG. 3. This shows the speciﬁc heat for small values of the

magnetic ﬁeld. The labels on the curve indicate the results

for x0 = 0.1, 0.5, 1. For comparison the result for no magnetic

x in units of x 0 field is shown in the discontinuous curve. FIG. 2. This shows the ratio of the magnetization computed using the Schafroth approximation in eq. (6) to the exact result. The labels 100,10,1,0.1 on the curves indicate Although it might be thought that there was very little the value of x0, the smaller numbers corresponding to lower difference between our result of ǫ ≃ 0.3 and the result values for the magnetic field. ǫ ≃ 0.5 of Ref. [10], in fact the magnetization and other thermodynamic quantities are extremely sensitive to the The ratio of the Schafroth approximation to the true value of ǫ. Using our expansion (8) we find result is plotted in fig. 2. As B is lowered the Schafroth 1/2 approximation becomes better over a wider range of tem- ρ 6π ζH (−1/2,a) 1/2 M ≃− x0 perature, and as T → T0. As noted by Schafroth [1], the 2m ζ(3/2) n approximation will break down when T becomes too close ζ(1/2) 1 2 −2 ( − a + a )x0 , (9) to T0, and this is borne out by the detailed calculation. ζ(3/2) 6 For T ≃ T0/10 the agreement is within 2%. o 1/2 It is possible to derive the Schafroth approximation in for small x0. This shows the −B behaviour found in (6) from the exact result (7) by approximating for the Ref. [10], but with a slightly different numerical factor class of sums (5). A powerful technique for doing this due to the lower temperature range we are using. makes use of the Mellin-Barnes integral representation In addition to the magnetization we have obtained an for the exponential function [10,3]. In Ref. [10] it was approximate expression for the specific heat using (8) shown that for ǫ ≃ 1/2 (corresponding to µ ≃ 0) there which is

was a superdiamagnetic regime where the magnetization M ∝−B1/2. (A previous paper which studied the mag- CV 15 ζ(5/2) 9ζ(3/2) 1/2 ≃ − 1/2 x0 , (10) netization is Ref. [11], although the superdiamagnetism kN 4 ζ(3/2) 2π ζH (3/2,a) was not found.) The temperature range at which ǫ ≃ 1/2 is close to but slightly above T0. We have obtained re- at T = T0. The first term in (10) which survives in the sults in a temperature range which includes T = T0. By B → 0 limit may be observed [8] to be the specific heat solving (4) for ǫ we find for the free Bose gas at the critical temperature T0. The

3 presence of a non-zero magnetic field is seen to lower the In conclusion, we have presented the results of a study specific heat from the free field result which is consistent of the ideal charged Bose gas in a constant magnetic field. with the results shown in fig. 1. It is possible to system- For certain values of the temperature and magnetic field atically determine higher order terms in the expansions we have obtained analytical results and compared them (8–10) and increase the accuracy of the approximation. with numerical ones. By examining the specific heat it This and other details will be presented elsewhere [12]. was shown how in a strong magnetic field the system As an example of a magnetic system where a phase exhibits a “dimensional reduction” where the system be- transition does occur we will consider the case D = 5. haves as if the spatial dimension was lowered by 2. We When B = 0 the specific heat is discontinuous at the have also examined the way in which the free Bose gas critical temperature T0. For B 6= 0 the specific heat is results are approached for small B. Although there is given by an expression similar to (3), and is perfectly a significant difference between the case B = 0 (which continuous. The results are shown in fig. 3. has a phase transition) and the B 6= 0 gas (which has The critical temperature Tc for the phase transition no phase transition), for small values of B there is very cannot be computed analytically if B 6= 0; however for little difference between the two specific heat curves. Fi- small values of the magnetic field it is possible to show nally for spatial dimensions D ≥ 5 a proper account of that [5] the condensate is essential for obtaining correct results, and this can be done in an efficient manner by using the ζ(3/2) eB effective action method. kTc ≃ kT0 − . (11) 5ζ(5/2) m GBS is grateful to the EPSRC for grant 94004194. DJT would like to thank K. Kirsten, F. Laloe, and S. The critical temperature is lowered from that for the free Ouvry for helpful discussions. We are both grateful to J. Bose gas when there is a non-zero magnetic field. The ap- Daicic and N. Frankel for helping us to clarify comments proximation in (11) is in good agreement with the exact concerning Refs. [10,11]. result for small B but breaks down as B becomes large [12]. The results in fig. 3 show that as B is decreased the curve for the specific heat starts to resemble the discontinuous one for the free gas. The critical temperature approaches the free field result, which is clear from the small B approximation in (11). The height of the specific heat maximum is seen to increase in magnitude and [1] M. R. Schafroth, Phys. Rev. 100, 463 (1955). the slope of the curve for T >Tc becomes steeper as [2] H. P. Rojas, Phys. Lett. B 379, 148 (1996). B is reduced. On the other hand as the magnetic field [3] K. Kirsten and D. J. Toms, Phys. Lett. B 368, 119 is increased in magnitude the specific heat undergoes a (1996); Phys. Rev. D, in press (cond-mat/9703161). marked change. The results for larger values of the mag- [4] R. M. May, J. Math. Phys. 6, 1462 (1965). netic field are shown in fig. 4. [5] D. J. Toms, unpublished, cond-mat/9612003. [6] L. D. Landau and E. M. Lifshitz, Statistical Physics,

(Pergammon, London, 1969). [7] D. J. Toms, Phys. Rev. D 51, 1886 (1995); Phys. Lett. B

343, 259 (1995).

k N [8] F. London, Superﬂuids, Vol. II, (John Wiley, New York,1954). [9] R. M. May, Phys. Rev. A 135, 1515 (1964); R. M. Ziﬀ,

G. E. Uhlenbeck, and M. Kac, Phys. Rep. 32C, 169

(1977).

[10] J. Daicic and N. E. Frankel, Phys. Rev. D 53, 5745

(1996). sp ecic heat in units of [11] T. A. Arias and J. D. Joannopoulos, Phys. Rev. B 39, 4071 (1989).

[12] G. B. Standen and D. J. Toms, paper in preparation.

x in units of x c FIG. 4. This shows the speciﬁc heat for small values of the magnetic ﬁeld. The labels on the curve indicate the results for x0 = 1, 2, 3, 5, 10.

For very large values of B the curves for the speciﬁc heat resemble the result found for the free Bose gas in three spatial dimensions. This provides a further and striking example of the “dimensional reduction” in the presence of a magnetic ﬁeld.