Kinetic Theory of a Spin-12 Bose-Condensed Gas

a rXiv:h ep-th/9213v11O ct1992Bose-Einstein condensation of scalar fields on hyperbolic manifolds Guido Cognola and Luciano Vanzo Dipartimento di Fisica -Universit`a di Trento ∗,Italia and Istituto Nazionale di Fisica Nucleare,Gruppo Collegato di Trento september 1992PACS numbers:03.70Theory of quantized fields05.90Other topics in statistical physics and thermodynamics1IntroductionBose-Einstein condensation for a non relativistic ideal gas has a long history[1].The physical phenomenon is well described in many text books (see for example ref.[2])and a rigorous mathematical discussion of it was given by many authors[3,4].The generalization to a relativistic idel Bose gas is non trivial and only recently has been discussed in a series of papers[5,6,7].It is well known that in the thermodynamic limit(infinite volume andfixed density)there is a phase transition of thefirst kind in correspondence of the critical temperature at which the condensation manifests itself.At that temperature,the first derivative of some continuous thermodynamic quantities has a jump.If the volume is keepedfinite there is no phase transition,nevertheless the phenomenon of condensation still occurs,but the critical temperature in this case is not well defined.For manifolds with compact hyperbolic spatial part of the kind H N/Γ,Γbe-ing a discrete group of isometries for the N-dimensional Lobachevsky space H N, zero temperature effects as well asfinite temperature effects induced by non-trivial topology,have been recently studied in some detail[8,9,10,11,12,?,14,15,16]. To our knowledge,a similar analysis has not yet been carried out for non compact hyperbolic manifolds.Hyperbolic spaces have remarkable properties.For example,the continuous spec-trum of the Laplace-Beltrami operator has a gap determined by the curvature ra-dius of H N,implying that masslessfields have correlation functions exponentially decreasing at infinity(such a gap is not present for the Dirac operator).For that reason,H4was recently proposed as an excellent infrared regulator for massless quantumfield theory and QCD[17].Critical behaviour is even more striking.In two flat dimension vortex configurations of a complex scalarfield,the XY model for He4films have energy logarithmically divergent with distance,while on H2it isfinite. This implies that the XY model is disordered at anyfinite temperature on H2.Even quantum mechanics on H2has been the subject of extensive investigations[18].The manifold H4is also of interest as it is the Euclidean section of anti-de Sitter space which emerges as the ground state of extended supergravity theories.The stress tensor on this manifold has been recently computed for both boson and fermion fields using zeta-function methods[19].In the present paper we shall discussfinite temperature effects and in particular the Bose-Einstein condensation for a relativistic ideal gas in a3+1dimensional ultrastatic space-time M=R×H3.We focus our attention just on H3,because such a manifold could be really relevant for cosmological and astrophysical applications.To this aim we shall derive the thermodynamic potential for a charged scalarfield of mass m on M,using zeta function,which on H3is exactly known.We shall see that the thermodynamic potential has two branch points when the chemical potentialµriches±ωo,ω2o=κ+m2being the lower bound of the spectrum of the operator L m=−△+m2and−κthe negative constant curvature of H3.The values±ωo will be riched byµ=µ(T)of course for T=0,but also for T=T c>0.This is the critical temperature at which the Bose gas condensates.The paper is organized as follows.In section2we study the elementary properties of the Laplace-Beltrami operator on H3;in particular we derive its spectrum and build up from it the related zeta-function.In section3we briefly recall how zeta-function can be used in order to regularize the partition function and we derive the regularized expression for the thermodynamic potential.In section4we discuss the Bose-Einstein condensation and derive the critical temperatures in both the cases of low and high temperatures.In section5we consider in detail the low and high temperature limits and derive the jump of thefirst derivative of the specific heat.The paper end with some considerations on the results obtained and some suggestions for further developments.2The spectrum and the zeta function of Laplace-Beltrami operator on H3For the aims of the present paper,the3-dimensional Lobachevsky space H3can be seen as a Riemannian manifold of constant negative curvature−κ,with hyperbolic metric dl2=d̺2+sinh2̺(dϑ2+sin2ϑdϑ2)and measure dΩ=sinh2̺d̺dΣ,dΣbeing the measure on S2.For convenience,here we normalize the curvature−κto−1.In these coordinates,the Laplace-Beltrami operator△reads△=∂2∂̺+1so,in order to derive it,it is sufficient to study radial wave functions of−△,that is solution of equationd2ud̺+λu=0(2) which reduces tod2vνsinh̺(4) Now,the L2(dΩ)scalar product for uν(̺)is(uν,uν′)=4πν2δ(ν−ν′)(5)from which the density of states̺(ν)=Vν2/2π2directly follows.As usual,we have introduced the large,finite volume V to avoid divergences.When possible,the limit V→∞shall be understood.At this point the computation of zeta function is straightforward.As we shall see in the following,what we are really interested in,is the zeta function related to the operators Q±=L1/2m±µ.The eigenvalues of L m areω2(ν)=ν2+a2=ν2+κ+m2,then we getζ(s;Q±)=V(4π)3/2Γ(s−1/2)(2a)s−3F(s+1,s−3;s−12a)(6)where F(α,β;γ;z)is the hypergeometric function.For its properties and its integral representations see for example ref.[20].It has to be noted that eq.(6)is the very same one has on aflat space for a massivefield with mass equal to a.Here in fact, the curvature plays the role of an effective mass.As we see from eq.(6),the zeta function related to the pseudo-differential oper-ators Q±has simple poles at the points s n=3,2,1,−1,−2,−3,...with residues b n(±µ)=Res(ζ(s;Q±),s n)given byb3(±µ)=Vπ2;b1(±µ)=V2a);(8)c−n=(−1)n nV(2a)n+3dz F(α,β;γ;z)=αβgdV(11)where Lµ=−(∂τ−µ)2+L m and the Wick rotationτ=ix0has to be understood. In eq.(11)the integration has to be taken over allfieldsφ(τ,x a)withβ-periodicity with respect toτ.The eigenvalues of the whole operator Lµ,sayµn,νreadµn,ν= 2πnlog det(ℓ−2Lµ)=−1β[logℓ2ζ(0;Lµ)+ζ′(0;Lµ)](13)βℓbeing an arbitrary normalization parameter coming from the scalar path-integral measure.Note thatℓ,which has the dimensions of a mass,is necessary in order to keep the zeta-function dimensionless for all s.Thefinite temperature andµdependent part of the thermodynamic potential does not suffer of the presence of such an arbitrary parameter.On the contrary,ℓenters in the regularized expression of vacuum energy and this creates an ambiguity[12],which is proportional to the heat kernel expansion coefficient K N(L m)related to L m(in general,K N(L m)=0). When the theory has a natural scale parameter,like the mass of the particle or the constant curvature of the manifold,the ambiguity can be removed by an”ad hoc”choice ofℓ[23].Here we would like to study the behaviour of thermodynamic quantities,then we are only interested in theµand T dependent part of the thermodynamic potential; that is a well defined quantity,which does not need regularization.To compute it, it is not necessary to use all the analytic properties of zeta function(for a careful derivation of vacuum energy see for example refs.[12,14]).Then we can proceed in a formal way and directly compute log det Lµdisregarding the vacuum energy divergent term.First of all we observe that∞ n=−∞log(ω2+(2πn/β+iµ)2)=∞ n=−∞ dω24ω cothβ2(ω−µ) dω2(14)=−log 1−e−β(ω+µ) −log 1−e−β(ω−µ) −βωUsing eq.(13),recalling thatω2=ν2+a2and by integrating overνwith the state density that we have derived in the previous section,we get the standard result1Ω(β,µ)=−2π2β log 1−e−β(ω(ν)+µ) +log 1−e−β(ω(ν)−µ) ν2dν(15)V+∞ n=1cosh nβµK2(anβ)π2∞ n=0µ2n2;L m)β−s ds(17)πi1E(β,µ)=−where K2is the modified Bessel function,c is a sufficiently large real number and ζR(s)is the usual Riemann zeta-function.The integral representations(17)and (18),which are valid for|µ|<a,are useful for high temperature expansion.On the contrary,the representation(16)in terms of modified Bessel functions is more useful for the low temperature expansion,since the asymptotics of Kνis well known.4Bose-Einstein condensationIn order to discuss Bose-Einstein condensation we have to analyze the behaviour of the charge density∂Ω(V,β,z)ρ=zV(expβωj−z)(20) and the activity z=expβµhas been introduced.Theωj in the sum are meant to be the Dirichlet eigenvalues for any normal domain V⊂H3.That is,V is a smooth connected submanifold of H3with non empty piecewise C∞boundary.By the infinite volume limit we shall mean that a nested sequence of normal domains V k has been choosen together with Dirichlet boundary conditions and such that V k≡H3.The reason for this choice is the following theorem due to Mac Kean (see for example[24]):—ifωok denotes the smallest Dirichlet eigenvalue for any sequence of normal domains V kfilling all of H3thenωok≥a and lim k→∞ωok=a.(Although the above inequality is also true for Neumann boundary conditions,the existence of the limit in not assured to the authors knowledge).Now we can show the convergence of thefinite volume activity z k to a limit point ¯z as k→∞.Tofix ideas,let us supposeρ≥0:then z k∈(1,expβωok).Since ρ(V,β,z)is an increasing function of z such thatρ(V,β,1)=0andρ(V,β,∞)=∞, for eachfixed V k there is a unique z k(¯ρ,β)∈(1,βexpωok)such that¯ρ=ρ(V k,β,z k).By compactness,the sequence z k must have at least one fixed point ¯z and as ωok →a 2as k goes to infinity,by Mc Kean theorem,¯z ∈[1,exp βa ].From this point on,the mathematical analysis of the infinite volume limit exactly parallels the one in flat space for non relativistic systems,as it is done in various references [25,26,4].Inparticular,there is a critical temperature T c over which there are no particles in the ground state.T c is the unique solution of the equation̺=sinh βa cosh βaV ∂µ= ∞0 1e β(ω(ν)−µ)−1 ν22π2 ∞0ν2dνκ+m 2.Thevery difference between flat and hyperbolic spaces occurs for massless particles.We shall return on this important point in a moment.Solutions of eq.(21)can be easily obtained in the two cases βa ≫1and βa ≪1(in the case of massive bosons these correspond to non relativistic and ultrarela-tivistic limits respectively).We have in fact̺≃T 3e x 2/2a −1= aT 2π2 ∞0x 2dx3;βa ≪1(24)from which we get the corresponding critical temperaturesT c =2πζR (3/2)2/3;βa ≫1(25)T c = 3̺∂µ̺(T,µ)(29)and since ∂µ̺diverges for µ=a we obtain µ′(T +c )=0.This is not the case of µ′′.In fact we shall see that µ′′(T +c )is different from zero and therefore µ′′(T )isadiscontinuous function of temperature.This implay that thefirst derivative of the specific heat C V has a jump for T=T c given bydC VdT T−c=µ′′(T+c)∂U(T,µ)∂TT=T+c(30)U(T,µ)being the internal energy,which can be derived by means of equation U(β,µ)=−µ̺V+∂πΓ(k+1)Γ(−k+5/2)(2s)−k(32) Then,for small T we haveE(β,µ)≃−a4Vanβ5/2∞ k=0Γ(k+5/2) 2π 3/2∞ n=1e−nβ(a−|µ|)T A2;T∂̺2̺(35)where A=2.363and C=−2.612are two coefficients of the expansion ∞ n=1e−nxNow,using eq.(30),we have the standard resultdC VdT T−c=3̺C2T c(37) The high temperature expansion could be obtained by using eq.(17),like in ref.[15].Here we shall use eq.(18),because for the aim of the present paper it is more ing the properties ofζ(s;Q±),which we have discussed in section2,we see that the integrand function in eq.(18)ζR(s+1)Γ(s)[ζ(s;Q+)+ζ(s;Q−)]β−(s+1)(38) has simple poles at s=3,1,0,−3,−5,−7,...and a double pole at s=−1.In-tegrating this function on a closed path containing all the poles,we get the high temperature expansion,valid for T>T c(hereγis the Euler-Mascheroni constant)E(β,µ)≃−V π212β2(a2−2µ2)+(a2−µ2)3/224π2(3a2−µ2)+a44π+γ−3(−2π)n(2n+1)where we have used the formulaζ′R(−2n)=Γ(2n+1)ζR(2n+1)3+µT(a2−µ2)1/212π2(41)−2(−2π)nForµ=a,the leading term of this expression gives again the result(24).From eq.(41),by a strightforward computation and taking only the leading terms into account,one getsµ′′(T+c)≃−12π2andfinally,from eqs.(30)and(31) dC VdT T−c≃−32̺π2pands adiabatically,can represent a manifold of the form we have considered.The problem we have studied then canfind physical applications in the standard model of the universe.References[1]A.Einstein.Berl.Ber.,22,261,(1924).[2]K.Huang.Statistical Mechanics.J.Wiley and Sons,Inc.,New York,(1963).[3]H.Araki and E.J.Woods.J.Math.Phys.,4,637,(1963).[4]ndau and m.Math.Phys.,70,43,(1979).[5]H.E.Haber and H.A.Weldom.Phys.Rev.Lett.,46,1497,(1981).[6]H.E.Haber and H.A.Weldom.J.Math.Phys.,23,1852,(1982).[7]H.E.Haber and H.A.Weldom.Phys.Rev.D,25,502,(1982).[8]A.A.Bytsenko and Yu.P.Goncharov.Mod.Phys.Lett.A,6,669,(1991).[9]Yu.P.Goncharov and A.A.Bytsenko.Class.Quantum Grav.,8,L211,(1991).[10]A.A.Bytsenko and Y.P.Goncharov.Class.Quantum Grav.,8,2269,(1991).[11]A.A.Bytsenko and S.Zerbini.Class.Quant.Gravity,9,1365,(1992).[12]G.Cognola,L.Vanzo and S.Zerbini.J.Math.Phys.,33,222,(1992).[13]A.A.Bytsenko,L.Vanzo and S.Zerbini.Mod.Phys.Lett.A,7,397,(1992).[14]A.A.Bytsenko,G.Cognola and L.Vanzo.Vacuum energy for3+1dimesionalspace-time with compact hyperbolic spatial part.Technical Report,Universit´a di Trento,UTF255,(1992).to appear in JMP.[15]G.Cognola and L.Vanzo.Thermodynamic potential for scalarfields in space-time with hyperbolic spatial part.Technical Report,Universit´a di Trento,UTF 258,(199).to be published.[16]A,A.Bytsenko,L.Vanzo and S.Zerbini.Zeta-function regularization approachtofinite temperature effects in Kaluza-Klein space-time.Technical Report,Uni-versit´a di Trento,UTF259,(1992).to be published.[17]C.G.Callan and F.Wilczek.Nucl.Phys.B,340,366,(1990).[18]N.Balasz and C.Voros.Phys.Rep.,143,109,(1986).[19]R.Camporesi and A.Higuchi.Phys.Rev.D,45,3591,(1992).[20]I.S.Gradshteyn and I.M.Ryzhik.Table of integrals,series and products.Aca-demic press,Inc.,New York,(1980).[21]A.Actor.Phys.Lett.B,157,53,(1985).[22]m.Math.Phys.,55,133,(1977).[23]J.S.Dowker and J.P.Schofield.Nucl.Phys.B,327,267,(1989).[24]I.Chavel.Eingenvalues in Riemannian Geometry.Accademic Press,(1984).[25]R.Ziff,G.E.Uhlenbeck and M.Kac.Phys.Rep.,32,169,(1977).[26]J.T.Lewis and J.V.Pul´m.Math.Phys.,36,1,(1974).[27]J.M.Blatt and S.T.Butler.Phys.Rev.Lett.,100,476,(1955).。

a r X i v :n u c l -t h /9308015v 1 19 A u g 1993Bose-condensation through resonance decayU.Ornik 1∗,M.Pl¨u mer 2†and D.Strottman 3‡1Theory Group,Gesellschaft f¨u r Schwerionenforschung (GSI),Darmstadt,F.R.Germany 2Physics Department,University of Marburg,FRG 3Theoretical Division,Los Alamos National Laboratory,Los Alamos,USA 87545February 9,2008Abstract We show that a system described by an equation of state which contains a high number of degrees of freedom (resonances)can create a considerable amount of super-fluid (condensed)pions through the decay of short-lived resonances,if baryon number and entropy are large and the dense matter decouples from chemical equilibrium ear-lier than from thermal equilibrium.The system cools down faster in the presence of a condensate,an effect that may partially compensate the enhancement of the lifetime expected in the case of quark-gluon-plasma formation.M@VAX.HRZ.UNI-MARBURG.DE‡E.Mail:DDS@1The investigation of hot hadronic matter through heavy ion collisions has shown that an understanding of medium effects andfinal state interactions is of great importance for the interpretation of the particle spectra in order to distinguish new physical phenomena like the formation of quark-gluon plasma(QGP)from background phenomena caused by “conventional”physical effects.We refer here,e.g.,to the strangeness enhancement,the J/Ψsuppression and the soft-pion puzzle,which can in principle be explained either by the formation of a QGP or from purely hadronic origin(for a recent review,see ref.[1]).An interesting explanation for the soft pion puzzle was proposed in[2].The authors argued that the baryonic and mesonic resonances may decouple from chemical equilibrium but still remain in thermal equilibrium.Since the resonance production rates fall drastically after their decoupling from chemical equilibrium,the remaining short-lived resonances would then rapidly decay into pions.As a consequence,the pions may acquire a non-zero chemical potentialµπwhich leads to a softer p⊥-spectrum.A different possible explanation for the absence of local chemical equilibrium of pions in thefinal(decoupling)stage,which should not be confused with the scenario discussed in ref.[2],was proposed in refs.[3,4]where it is assumed that the pions are initially created out of local chemical(and even thermal!)equilibrium and it turns out that the system never reaches complete local chemical equilibrium[3,4].In the present letter,we investigate the scenario of ref.[2],where it is assumed that the dense matter reaches local thermodynamic equilibrium in an early stage of the expansion. The authors of[2]considered only the case that the chemical potential stays below the pion mass.Moreover,they did not attempt to determine the value of the pion chemical potential in terms of the temperature and baryon number density at the point of decou-pling from local chemical equilibrium.Below,we address the problem if and under what conditions the decay of short-lived resonances can lead to the formation of a pionic Bose condensate,and what happens to the condensate in the subsequent expansion until the system decouples from local thermal equilibrium and thefinal state particles are emitted. To do so,we explicitly consider the transition from a resonance gas to a gas of stable and long-lived particles for a hadronic equation of state which contains the known resonances up to masses of2GeV,imposing constraints from the conservation of energy-momentum,2baryon number and strangeness.In particular,we shall show that for a system rich in entropy and baryon numberµπmay reach the pion mass and,consequently,the pions may form a Bose condensate in the central region in relativistic heavy ion collisions.For a hydrodynamically expanding system,the formation of such a condensate implies the presence of a superfluid component.The idea that the hadronic matter might be superfluid was put forward a long time ago in a somewhat different context[5].This could explain in a natural way the presence of a coherent component in multiparticle production.The implications for pion interferometry measurements will be discussed in a separate paper[6].We note that although the creation of a remarkable amount of condensed pions is interesting by itself,these effects could also be used to answer the question of whether pions and resonances are in chemical equilibrium before they decouple from thermal equilibrium.Because thefinal stage in the hadronic phase can have a significant influence on the particle spectrum,the question of whether the freeze-out occurs in chemical and thermal equilibrium,as assumed in[7,8]or only in thermal equilibrium[2]is of great importance.In order to investigate the conditions necessary to form a condensate,let usfirst consider a system in local thermodynamic(i.e.,thermal and chemical)equilibrium.The medium is then completely described in terms of its equation of state(EOS)which can be expressed in the formεfl=f(T fl,µk fl)(1)whereεfl is the energy-density of thefluid,T fl the temperature,and the quantitiesµk fl are the chemical potentials related to conserved charges Q k(such as baryon number B and strangeness S).The index fl refers to quantities which describe thefluid.During a heavy ion collision we expect under certain conditions the formation of a system of hot hadronic matter or a QGP described by eq.(1)which is in chemical and thermal equilibrium.The system then undergoes a phase of hydrodynamic expansion until it has become so dilute that the interactions are no longer strong enough to maintain local equi-3librium.Thisfinal stage,when the collective behaviour terminates and particles decouple from the dense matter,is referred to as freeze-out.The scenario of ref.[2]differs from most of the other hydrodynamic descriptions of heavy ion collisions in so far as in[2]the dense matter is assumed to go out of chemical equilibrium but remain in thermal equilibrium for some time.Or,to put it differently,there will be a“chemical freeze-out”independent of and prior to the“thermal freeze-out”.As was mentioned above,the transition out of chemical equilibrium leads to the decay of short-lived resonances and the appearance of a non-zero chemical potential of the pions.Below,we shall assume that the chemical freeze-out,just like the thermal freeze-out,occurs “instantaneously”(i.e.,on a three-dimensional hypersurface characterized by a condition such as T(x,t)=const.).This seems a reasonable assumption since the relevant time scales are the decay times of the short-lived resonances which are on the order of∼1 fm/c.To illustrate the rapidness of such a transition,we consider the simplified case of a system of pions and rho-mesons in local thermal equilibrium that undergoes a longitudinal scaling[9]expansion.In addition to the decayρ→ππwe also take into account the inverse processππ→ρ+X.The equations which describe the time evolution of the system aredετ(2)dnπτ+2nρdτ=−nρτρ+ σv n2π(4)whereε=ερ+επand P=Pρ+Pπare the energy density and pressure,n i(i=ρ,π) are the number densities andτ≡√transition is even faster.Thefigure also shows that the decay of rho-mesons alone is not sufficient to drive the pions into the condensate(for that,one needs to take into account the decay contributions of all the other short-lived resonances as well).We now return to the description of the expanding hot and dense matter at the point of chemical freeze-out.It will be assumed that thefluid freezes out into“stable”π’s,K’s, nucleons,Λ’s and the long living resonancesω’s andη’s.The term“stable”here means that the decay of the particle or resonance occurs only after the complete(thermal and chemical)freeze-out of the system.In particular,the contribution ofω’s andη’s to the pion chemical potential is zero.Nevertheless,their later decay leads to a non-thermal component of the pion spectrum which is determined by the decay kinematics[7].The system of chemically frozen out particles can be described by a chemical decoupling temperature T ch.f.,chemical potentials for baryons and strange particles,µB andµS, and a pion chemical potential,µπ,which describes the overpopulation of pions due to the decoupling of resonances.The chemical potentials are determined by the requirement that the energy densityεfl,the baryonic density b fl and the strangeness density s fl=0of the fluid are equal to those of the system after chemical freeze-out.We obtain the following system of equations:εfl=εthermπ(µπ,T ch.f.)+mπn conπΘ(mπ−µπ)+i=K,N,Λ,ω,η,...εthermi(µS,µB,T ch.f.)(5)b fl= k B k n therm k(µS,µB,µπ,T ch.f.)(6)s fl= k S k n therm k(µS,µB,µπ,T ch.f.)=0(7) whereεtherm i (µS,µB,µπ,T)=g iexp E i−˜µi2π2p2dp1T ±1(9)5are the thermal parts of the energy density and the number density of particle species i, and where we have used the notation˜µi= k Q k iµk(10) for the chemical potentials.Note that on the r.h.s.of eq.(10)the charges Q k include the pion number Nπwhich is a conserved quantity after chemical freeze-out.In eq.(5)we have introduced a term n conπthat describes the condensed component of pions which will appear if the chemical equilibrium for pions reaches the pion mass.In this case an overpopulation of pion states occurs and the remaining pions will be forced into the ground state,thereby forming a Bose condensate.Thefluid then consists of a thermal component and a superfluid component.The latter moves with the rapidity of thefluid, whereas the thermal pions have on top of thefluid rapidity a thermal distribution.The resultant single inclusive distribution of pions emitted from a hydrodynamically expanding source1dydp⊥=gπexp p iµuµ−˜µπresonances and a treatment of compression effects,thus retaining the essential features of nuclear matter near the ground state[11].In addition to the work presented in[11]we take into account the effects offinite baryonic and mesonic masses,use Fermi statistics for the baryons and include mesons up to masses of2GeV.Fig.2shows the thermal and the condensate component of pions at chemical freeze-out, and the chemical potentialsµπ,µb andµs,as functions of the temperature T ch.f.that characterizes the decoupling from chemical equilibrium,for three different values of the baryon number density b.It can be seen how the baryonic and strange chemical potentials,µb andµs,decrease with increasing temperature T ch.f..Clearly,the strange chemical potential induced by thefinite baryon number never execeeds150MeV.As the amount of heavier resonances increases with temperature,the pion chemical potentialµπgrows with T ch.f..It becomes equal to the pion mass in the temperature region170-200MeV;the value of T ch.f.whereµπ=mπdepends on the baryonic density.This behaviour is also reflected in the dependence of the thermal and condensate densities of pions on the temperature and baryon number density at chemical freeze-out.It can be seen that at sufficiently high baryonic densities and temperatures the condensate contribution to the chemically decoupled pions becomes significantly large.It is noteworthy that the condensation effect is enhanced if the density of the strange particles has not yet reached its equilibrium value. In this case the remaining energy is distributed among the non strange resonances,and this leads to an enhanced production of particles in the condensate.Finally,we need to discuss the question whether or not the condensate survives the time period between the chemical and the thermal freeze-out.This is of particular importance if the dense matter decouples from chemical equilibrium at temperatures T ch.f.∼170−200 MeV which are considerably higher than those usually associated with the thermal freeze-out,T th.f.∼mπ.After decoupling from chemical equilibrium,the system of stable and long-lived particles continues to expand until it has cooled down to the temperature T th.f..For a one-dimensional scaling expansion,the evolution of the system is determined by the equations7s τ=const .(12)b τ=const .(13)s e τ=const .(14)n πτ=const .(15)n ατ=const .(16)which describe the conservation of entropy S e ,baryon number B ,strangeness S and pion number N π(s e ,b ,s and n πbeing the corresponding densities,respectively).The index αlabels those of the stable particle species (K,N,Λ,ω,η,...)which also have decoupled from chemical equilibrium.For the gas of stable particles and long-lived resonances,the thermodynamic quantities are given by the expressions on the r.h.s.of eqs.(5)–(10),with s e =i (εi +P i −˜µi n i )/T .Energy density and number density of the pions consistof a thermal and a condensate component,επ=εtherm π+m πn con πand n π=n therm π+n con π,respectively.Note that the condensate does not contribute to the entropy.Consequently,the cooling of the system –i.e.,the function T (τ)–does not depend on the superfluid component.Fig.3shows the cooling curves T (τ)of a pion gas,for different values of the pion chemical potential.For comparison we have also included the results for an ideal gas of massless pions.It can be seen that the cooling rate increases with increasing pion chemical potential and becomes maximal in the presence of a condensate.This implies that Bose condensates may reduce the lifetime of the fireball.In particular,we are interested in the time dependence of the fraction of pions in thecondensate,f con ≡n con π/(n therm π+n con π).Eqs.(12,15)imply that for the fraction of thermal pions,f therm =1−f con ,f th.f.therm =f ch.f.therm n therm π/s e ) th.f.the pion gas,it turns out that in the temperature range between200MeV and150MeV, f therm varies by less than10%.To summarize,we have shown that the excited hadronic matter created in ultrarelativistic heavy ion collisions could form a pion condensate if(i)the system is rich in baryon number and entropy density,(ii)the hadrons decouple from local chemical equilibrium earlier than from thermal equilibrium,as suggested in[2],and(iii)the loss of chemical equilibrium occurs at temperatures of∼180MeV or higher.The effect depends on the chemical decoupling temperature,on the baryonic density and,rather sensitively,on the resonance contribution to the EOS.1The effect is increased if the chemical equilibrium for strange particles is not complete,i.e.if the amount of strangeness which is initially zero and starts to increase during the collision and expansion process has not yet reached its equilibrium value.It is not the purpose of the present study to prove that such a decoupling process occurs. To do so would require a detailed knowledge of the density and temperature dependence of the inelastic and elastic cross sections of hadrons in the dense matter.Rather,we would like to point out here that,if this decoupling happens in a certain temperature region(like it was quoted in[2]),then a pion condensate should appear.If a pion condensate is created in a heavy ion collision,it may be distinguishable from the case of a resonance gas in thermal and chemical equilibrium essentially by the following effects.•It should lead to a coherent component which then would appear in the two-particle correlation functions of identical pions through a reduction of the intercept of the correlation function and the appearance of a two exponent behaviour of the correla-tion function[13].•The condensate component is moving with thefluid velocity and has no additional thermal component which smears out the distribution.This might lead to character-istic bumps and shoulders in the rapidity and transverse momentum distributions of pions2.One also can expect that the coherent component disappears if one considers particles of velocities that exceed the maximumfluid velocities.•The lifetime of thefireball may be reduced.This could at least partially compensate the enhancement of the lifetime expected[15]in the case that quark-gluon-plasma is formed in the intial stage.Instructive and helpful discussions with R.M.Weiner,F.Navarra and H.Leutwyler are gratefully acknowledged.This work was supported in part by the Deutsche Forschungsge-meinschaft(DFG),the Federal Minister for Research and Technology(BMFT)under the contract no.06MR731,the Gesellschaft f¨u r Schwerionenforschung(GSI)and the Los Alamos National Laboratory.References[1]Proceedings of Quark Matter91,eds.T.C.Awes,F.E.Obenshain,F.Plasil,M.R.Strayer and C.Y.Wong,Nucl.Phys.A544(1990)513.[3]S.Gavin and P.V.Ruuskanen,Phys.Lett.B262(1992)1403.[5]A.Mann and R.M.Weiner,Nouvo Cim.10A(1972)383;S.Eliezer and R.M.Weiner,Phys.Rev D13(1991)503.[8]J.Sollfrank,P.Koch and U.Heinz,Phys.Lett.B252,(1982)140.[10]F.Cooper and G.Frye,Phys.Rev.D10(1989)3735.[12]U.Ornik,F.W.Pottag and R.M.Weiner,Phys.Rev.Lett.63(1989)278;B242(1990)181.[15]G.Bertsch,M.Gong,L.McLerran,V.Ruuskanen and E.Sarkkinen,Phys.Rev.D37(1988)1896;G.Bertsch,Nucl.Phys.A498FIGURE CAPTIONSFigure1:Left column:Dependence of the chemical potentials for baryons(dashed), strangeness(dotted)and pions(solid)on the chemical freeze-out temperature T ch.f.,for three different values of the baryon number density b.Right column:The thermal(solid) and condensate(dashed)components of the pion number densities at chemical freeze-out.Figure3:。

2022年自考专业(英语)英语科技文选考试真题及答案一、阅读理解题Directions: Read through the following passages. Choose the best answer and put the letter in the bracket. (20%)1、 (A) With the recent award of the Nobel Prize in physics, the spectacular work on Bose-Einstein condensation in a dilute gas of atoms has been honored. In such a Bose-Einstein condensate, close to temperatures of absolute zero, the atoms lose their individuality and a wave-like state of matter is created that can be compared in many ways to laser light. Based on such a Bose-Einstein condensate researchers in Munich together with a colleague from the ETH Zurich have now been able to reach a new state of matter in atomic physics. In order to reach this new phase for ultracold atoms, the scientists store a Bose-Einstein condensate in a three-dimensional lattice of microscopic light traps. By increasing the strength of the lattice, the researchers are able to dramatically alter the properties of the gas of atoms and can induce a quantum phase transition from the superfluid phase of a Bose-Einsteincondensate to a Mott insulator phase. In this new state of matter it should now be possible to investigate fundamental problems of solid-state physics, quantum optics and atomic physics. For a weak optical lattice the atoms form a superfluid phase of a Bose-Einstein condensate. In this phase, each atom is spread out over the entire lattice in a wave-like manner as predicted by quantum mechanics. The gas of atoms may then move freely through the lattice. For a strong optical lattice the researchers observe a transition to an insulating phase, with an exact number of atoms at each lattice site. Now the movement of the atoms through the lattice is blocked due to therepulsive interactions between them. Some physicists have been able to show that it is possible to reversibly cross the phase transition between these two states of matter. The transition is called a quantum phase transition because it is driven by quantum fluctuations and can take place even at temperatures of absolute zero. These quantum fluctuations are a direct consequence of Heisenberg’s uncertainty relation. Normally phase transitions are driven by thermal fluctuations, which are absent at zero temperature. With their experiment, the researchers in Munich have been able to enter a new phase in the physics of ultracold atoms. In the Mott insulator state theatoms can no longer be described by the highly successful theories for Bose-Einstein condensates. Now theories are required that take into account the dominating interactions between the atoms and which are far less understood. Here the Mott insulator state may help in solving fundamental questions of strongly correlated systems, which are the basis for our understanding of superconductivity. Furthermore, the Mott insulator state opens many exciting perspectives for precision matter-wave interferometry and quantum computing.What does the passage mainly discuss?A.Bose-Einstein condensation.B.Quantum phase transitions.C.The Mott insulator state.D.Optical lattices.2、What will the scientists possibly do by reaching the new state of matter in atomic physics?A.Store a Bose-Einstein condensate in three-dimensional lattice of microscopic light traps.B.Increase the strength of the lattice.C.Alter the properties of the gas of atoms.D.Examine fundamental problems of atomic physics.3、Which of the following is NOT mentioned in relation to aweak optical lattice?A.The atoms form a superfluid phase of a Bose-Einstein condensate.B.Each atom is spread out over the entire lattice.C.The gas of atoms may move freely through the lattice.D.The superfluid phase changes into an insulating phase.4、What can be said about the quantum phase transition?A.It can take place at temperatures of absolute zero.B.It cannot take place above the temperatures of absolute zero.C.It is driven by thermal fluctuations.D.It is driven by the repulsive interactions between atoms.5、The author implies all the following about the Mott insulator state EXCEPT that______.A.the theory of Bose-Einstein condensation can’t possibly account for the atoms in the Mott insulator stateB.not much is known about the dominating interactions between the atoms in the Mott insulator stateC.it offers new approaches to exact quantum computingD.it forms a superfluid phase of a Bose-Einstein condensate6、 (B) Gene therapy and gene-based drugs are two ways we would benefit from our growing mastery of genetic science. But therewill be others as well. Here is one of the remarkable therapies on the cutting edge of genetic research that could make their way into mainstream medicine in the c oming years. While it’s true that just about every cell in the body has the instructions to make a complete human, most of those instructions are inactivated, and with good reason: the last thing you want for your brain cells is to start churning out stomach acid or your nose to turn into a kidney. The only time cells truly have the potential to turn into any and all body parts is very early in a pregnancy, when so-called stem cells haven’t begun to specialize. Most diseases involve the death of healthy cells—brain cells in Alzheimer’s, cardiac cells in heart disease, pancreatic cells in diabetes, to name a few; if doctors could isolate stem cells, then direct their growth, they might be able to furnish patients with healthy replacement tissue. It was incredibly difficult, but last fall scientists at the University of Wisconsin managed to isolate stem cells and get them to grow into neural, gut, muscle and bone cells. The process still can’t be controlled, and may have unforeseen limitations; but if efforts to understand and master stem-cell development prove successful, doctors will have a therapeutic tool of incredible power. The same applies to cloning, whichis really just the other side of the coin; true cloning, as first shown, with the sheep Dolly two years ago, involves taking a developed cell and reactivating the genome within, resenting its developmental instructions to a pristine state. Once that happens, the rejuvenated cell can develop into a full-fledged animal, genetically identical to its parent. For agriculture, in which purely physical characteristics like milk production in a cow or low fat in a hog have real market value, biological carbon copies could become routine within a few years. This past year scientists have done for mice and cows what Ian Wilmut did for Dolly, and other creatures are bound to join the cloned menagerie in the coming year. Human cloning, on the other hand, may be technically feasible but legally and emotionally more difficult. Still, one day it will happen. The ability to reset body cells to a pristine, undeveloped state could give doctors exactly the same advantages they would get from stem cells: the potential to make healthy body tissues of all sorts. And thus to cure disease.That could prove to be a true “miracle cu re”.What is the passage mainly about?A.Tomorrow’s tissue factory.B.A terrific boon to medicine.C.Human cloning.D.Genetic research.7、 According to the passage, it can be inferred that which of the following reflects the author’s opinion?A.There will inevitably be human cloning in the coming year.B.The potential to make healthy body tissues is undoubtedly a boon to human beings.C.It is illegal to clone any kind of creatures in the world.D.It is legal to clone any kind of creatures in the world except human.8、Which of the following is NOT true according to the passage?A.Nearly every cell in the human brain has the instructions to make a complete human.B.It is impossible for a cell in your nose to turn into a kidney.C.It is possible to turn out healthy replacement tissues with isolated stem cells.D.There will certainly appear some new kind of cloned animal in the near future.9、All of the following are steps involved in true cloning EXCEPT_______.A.selecting a stem cellB.taking a developed cellC.reactivating the genome within the developed cellD.resetting the developmental instructions in the cell to its original state10、The word “rejuvenated” in para. 5 is closest in meaning to_______.A.rescuedB.reactivatedC.recalledD.regulated参考答案：【一、阅读理解题】1~5CDDAD6~10DBBA。

a r X i v :c o n d -m a t /9810197v 1 [c o n d -m a t .s t a t -m e c h ] 16 O c t 1998Accepted to PHYSICAL REVIEW A for publicationBose-Einstein condensation in a one-dimensional interacting system due to power-lawtrapping potentialsM.Bayindir,B.Tanatar,and Z.GedikDepartment of Physics,Bilkent University,Bilkent,06533Ankara,TurkeyWe examine the possibility of Bose-Einstein condensation in one-dimensional interacting Bose gas subjected to confining potentials of the form V ext (x )=V 0(|x |/a )γ,in which γ<2,by solving the Gross-Pitaevskii equation within the semi-classical two-fluid model.The condensate fraction,chemical potential,ground state energy,and specific heat of the system are calculated for various values of interaction strengths.Our results show that a significant fraction of the particles is in the lowest energy state for finite number of particles at low temperature indicating a phase transition for weakly interacting systems.PACS numbers:03.75.Fi,05.30.Jp,67.40.Kh,64.60.-i,32.80.PjI.INTRODUCTIONThe recent observations of Bose-Einstein condensation (BEC)in trapped atomic gases [1–5]have renewed inter-est in bosonic systems [6,7].BEC is characterized by a macroscopic occupation of the ground state for T <T 0,where T 0depends on the system parameters.The success of experimental manipulation of externally applied trap potentials bring about the possibility of examining two or even one-dimensional Bose-Einstein condensates.Since the transition temperature T 0increases with decreasing system dimension,it was suggested that BEC may be achieved more favorably in low-dimensional systems [8].The possibility of BEC in one -(1D)and two-dimensional (2D)homogeneous Bose gases is ruled out by the Hohen-berg theorem [9].However,due to spatially varying po-tentials which break the translational invariance,BEC can occur in low-dimensional inhomogeneous systems.The existence of BEC is shown in a 1D noninteracting Bose gas in the presence of a gravitational field [10],an attractive-δimpurity [11],and power-law trapping po-tentials [12].Recently,many authors have discussed the possibility of BEC in 1D trapped Bose gases relevant to the magnetically trapped ultracold alkali-metal atoms [13–18].Pearson and his co-workers [19]studied the in-teracting Bose gas in 1D power-law potentials employing the path-integral Monte Carlo (PIMC)method.They have found that a macroscopically large number of atoms occupy the lowest single-particle state in a finite system of hard-core bosons at some critical temperature.It is important to note that the recent BEC experiments are carried out with finite number of atoms (ranging from several thousands to several millions),therefore the ther-modynamic limit argument in some theoretical studies [15]does not apply here [8].The aim of this paper is to study the two-body interac-tion effects on the BEC in 1D systems under power-law trap potentials.For ideal bosons in harmonic oscillator traps transition to a condensed state is prohibited.It is anticipated that the external potentials more confin-ing than the harmonic oscillator type would be possible experimentally.It was also argued [15]that in the ther-modynamic limit there can be no BEC phase transition for nonideal bosons in 1D.Since the realistic systems are weakly interacting and contain finite number of particles,we employ the mean-field theory [20,21]as applied to a two-fluid model.Such an approach has been shown to capture the essential physics in 3D systems [21].The 2D version [22]is also in qualitative agreement with the results of PIMC simulations on hard-core bosons [23].In the remaining sections we outline the two-fluid model and present our results for an interacting 1D Bose gas in power-law potentials.II.THEORYIn this paper we shall investigate the Bose-Einstein condensation phenomenon for 1D interacting Bose gas confined in a power-law potential:V ext (x )=V 0|x |κF (γ)G (γ)2γ/(2+γ),(2)andN 0/N =1−TF (γ)=1x 1/γ−1dx1−x,(4)and G (γ)=∞x 1/γ−1/2dxNk B T 0=Γ(1/γ+3/2)ζ(1/γ+3/2)T 01/γ+3/2.(6)Figure 1shows the variation of the critical temperature T 0as a function of the exponent γin the trapping po-tential.It should be noted that T 0vanishes for harmonic potential due to the divergence of the function G (γ=2).It appears that the maximum T 0is attained for γ≈0.5,and for a constant trap potential (i.e.V ext (x )=V 0)the BEC disappears consistent with the Hohenberg theorem.0.00.5 1.0 1.5 2.0γ0.00.20.40.6k B T 0 (A r . U n .)FIG.1.The variation of the critical temperature T 0withthe external potential exponent γ.We are interested in how the short-range interactioneffects modify the picture presented above.To this end,we employ the mean-field formalism and describe the col-lective dynamics of a Bose condensate by its macroscopictime-dependent wave function Υ(x,t )=Ψ(x )exp (−iµt ),where µis the chemical potential.The condensate wavefunction Ψ(x )satisfies the Gross-Pitaevskii (GP)equa-tion [24,25]−¯h 2dx 2+V ext (x )+2gn 1(x )+g Ψ2(x )Ψ(x )=µΨ(x ),(7)where g is the repulsive,short-range interaction strength,and n 1(x )is the average noncondensed particle distribu-tion function.We treat the interaction strength g as a phenomenological parameter without going into the de-tails of actually relating it to any microscopic descrip-tion [26].In the semi-classical two-fluid model [27,28]the noncondensed particles can be treated as bosons in an effective potential [21,29]V eff(x )=V ext (x )+2gn 1(x )+2g Ψ2(x ).(8)The density distribution function is given byn 1(x )=dpexp {[p 2/2m +V eff(x )−µ]/k B T }−1,(9)and the total number of particles N fixes the chemical potential through the relationN =N 0+ρ(E )dE2mgθ[µ−V ext (x )−2gn 1(x )],(12)where θ[x ]is the unit step function.More precisely,the Thomas-Fermi approximation [7,20,30]would be valid when the interaction energy ∼gN 0/Λ,far exceeds the kinetic energy ¯h 2/2m Λ2,where Λis the spatial extent of the condensate cloud.For a linear trap potential (i.e.γ=1),a variational estimate for Λis given by Λ= ¯h 2/2m (π/2)1/22a/V 0 1/3.We note that the Thomas-Fermi approximation would breakdown for tem-peratures close to T 0where N 0is expected to become very small.The above set of equations [Eqs.(9)-(12)]need to be solved self-consistently to obtain the various physical quantities such as the chemical potential µ(N,T ),the condensate fraction N 0/N ,and the effective potential V eff.In a 3D system,Minguzzi et al .[21]solved a simi-lar system of equations numerically and also introduced an approximate semi-analytical solution by treating the interaction effects perturbatively.Motivated by the suc-cess [21,22]of the perturbative approach we consider aweakly interacting system in1D.To zero-order in gn1(r), the effective potential becomesV eff(x)= V ext(x)ifµ<V ext(x)2µ−V ext(x)ifµ>V ext(x).(13) Figure2displays the typical form of the effective po-tential within our semi-analytic approximation scheme. The most noteworthy aspect is that the effective poten-tial as seen by the bosons acquire a double-well shape because of the interactions.We can explain this result by a simple argument.Let the number of particles in the left and right wells be N L and N R,respectively,so that N=N L+N R.The nonlinear or interaction term in the GP equation may be approximately regarded as V=N2L+N2R.Therefore,the problem reduces to the minimization of the interaction potential V,which is achieved for N L=N R.FIG.2.Effective potential V eff(x)in the presence of in-teraction(x0=(µ/V0)1/γa).Thick dotted line represents external potential V ext(x).The number of condensed atoms is calculated to beN0=2γa√ze x−1+ 2µ/k B Tµ/k B TH(γ,µ,xk B T)(2µ/k B T−x)1/γ−1/2dxexp[(E−µ)/k B T]−1=κ(k B T)1/γ+1/2J(γ,µ,T),(18) whereJ(γ,µ,T)= ∞2µ/k B T x1/γ+1/2dxze x−1.and Ecis the energy of the particles in the condensateE c=g(1+γ)(2γ+1)gV1/γ.(19)The kinetic energy of the condensed particles is neglected within our Thomas-Fermi approximation to the GP equa-tion.III.RESULTS AND DISCUSSIONUp to now we have based our formulation for arbitrary γ,but in the rest of this work we shall present our re-sults forγ=1.Our calculations show that the results for other values ofγare qualitatively similar.In Figs. 3and4we calculate the condensate fraction as a func-tion of temperature for various values of the interaction strengthη=g/V0a(at constant N=105)and different number of particles(at constantη=0.001),respectively. We observe that as the interaction strengthηis increased, the depletion of the condensate becomes more apprecia-ble(Fig.3).As shown in the correspondingfigures,a significant fraction of the particles occupies the ground state of the system for T<T0.The temperature depen-dence of the chemical potential is plotted in Figs.5and 6for various interaction strengths(constant N=105) and different number of particles(constantη=0.001) respectively.0.00.20.40.60.8 1.0T/T 00.00.20.40.60.81.0N 0/NN 0/N=1−(T/T 0)3/2η=10−5η=10−3η=10−1η=10FIG.3.The condensate fraction N 0/N versus temperature T /T 0for N =105and for various interaction strengths η.Effects of interactions on µ(N,T )are seen as large de-viations from the noninteracting behavior for T <T 0.In Fig.7we show the ground state energy of an interacting 1D system of bosons as a function of temperature for dif-ferent interaction strengths.For small η,and T <T 0, E is similar to that in a noninteracting system.As ηincreases,some differences start to become noticeable,and for η≈1we observe a small bump developing in E .This may indicate the breakdown of our approxi-mate scheme for large enough interaction strengths,as we can find no fundamental reason for such behavior.It is also possible that the Thomas-Fermi approximation em-ployed is violated as the transition to a condensed state is approached.0.00.20.40.60.8 1.0T/T 00.00.20.40.60.81.0N 0/NN 0/N=1−(T/T 0)3/2N=108N=105N=103N=101FIG.4.The condensed fraction N 0/N versus temperature T /T 0for η=0.001and for different number of particles N .0.00.20.40.60.8 1.0 1.2T/T 0−100100200300400µ/V 0η=1η=0.1η=0.001η=0.00001FIG.5.The temperature dependence of the chemical potential µ(N,T )for various interaction strength and for N =105particles.Although it is conceivable to imagine the full solution of the mean-field equations [Eq.(9)-(12)]may remedy the situation for larger values of η,the PIMC simulations [19]also seem to indicate that the condensation is inhibited for strongly interacting systems.The results for the spe-cific heat calculated from the total energy curves,i.e.C V =d E /dT ,are depicted in Fig.8.The sharp peak at T =T 0tends to be smoothed out with increasing in-teraction strength.It is known that the effects of finite number of particles are also responsible for such a be-havior [20].In our treatment these two effects are not disentangled.It was pointed out by Ingold and Lam-brecht [14]that the identification of the BEC should also be based on the behavior of C V around T ≈T 0.0.00.20.40.60.8 1.0 1.2T/T 0−5050100µ/V 0N=107N=105N=103N=101FIG.6.The temperature dependence of the chemical po-tential µ(N,T )for different number of particles N and for η=0.001.0.00.20.40.60.8 1.0 1.2T/T 00.00.20.40.60.8<E >/N k B T 0η=0η=0.001η=0.1η=1Maxwell−BoltzmannFIG.7.The temperature dependence of the total energy of 1D Bose gas for various interaction strengths ηand N =105particles.Our calculations indicate that the peak structure of C V remains even in the presence of weak interactions,thus we are led to conclude that a true transition to a Bose-Einstein condensed state is predicted within the present approach.0.00.20.40.60.81.01.2T/T 00.00.20.40.60.81.0C V /N k Bη=0η=0.001η=0.1Maxwell−BoltzmannFIG.8.The temperature dependence of the specific heat C V for various interaction strengths ηand N =105particles.IV.CONCLUDING REMARKSIn this work we have applied the mean-field,semi-classical two-fluid model to interacting bosons in 1D power-law trap potentials.We have found that for a range of interaction strengths the behavior of the thermo-dynamic quantities resembles to that of non-interactingbosons.Thus,BEC in the sense of macroscopic occu-pation of the ground state,occurs when the short-range interparticle interactions are not too strong.Our results are in qualitative agreement with the recent PIMC sim-ulations [19]of similar systems.Both 2D and 1D sim-ulation results [19,23]indicate a phase transition for a finite number system,in contrast to the situation in the thermodynamic limit.Since systems of much larger size can be studied within the present approach,our work complements the PIMC calculations.The possibility of studying the tunneling phenomenon of condensed bosons in spatially different regions sepa-rated by a barrier has recently attracted some attention [31–34].In particular,Dalfovo et al .[32]have shown that a Josephson-type tunneling current may exist for bosons under the influence of a double-well trap potential.Za-pata et al .[34]have estimated the Josephson coupling energy in terms of the condensate density.It is inter-esting to speculate on such a possibility in the present case,since the effective potential in our description is of the form of a double-well potential (cf.Fig.2).In our treatment,the interaction effects modify the single-well trap potential into one which exhibits two minima.Thus if we think of this effective potential as the one seen by the condensed bosons and according to the general ar-guments [31–34]based on two weakly connected systems we should have an oscillating flux of particles when the chemical potential in the two wells is different.Any con-figuration with N L =N R which is always the case for odd number of bosons will result in an oscillatory mo-tion.It would be interesting to explore these ideas in future work.ACKNOWLEDGMENTSThis work was supported by the Scientific and Techni-cal Research Council of Turkey (TUBITAK)under Grant No.TBAG-1736and TBAG-1662.We gratefully ac-knowledge useful discussions with Prof.C.Yalabık and E.Demirel.[5]D.J.Han,R.H.Wynar,Ph.Courteille,and D.J.Heinzen,Phys.Rev.A57,R4114(1998).[6]I.F.Silvera,in Bose-Einstein Condensation,Ed.by A.Griffin,D.W.Snoke,and S.Stringari(Cambridge Uni-versity Press,Cambridge,1995).[7]F.Dalfovo,S.Giorgini,L.P.Pitaevskii,and S.Stringari,preprint,cond-mat/9806038(to be published in Reviews of Modern Physics);A.S.Parkins and D.F.Walls,Phys.Rep.303,1(1998).[8]W.Ketterle and N.J.van Druten,Phys.Rev.A54,656(1996).[9]P.C.Hohenberg,Phys.Rev.158,383(1967).[10]A.Widom,Phys.Rev.176,254(1968).[11]L.C.Ioriatti,Jr.,S.G.Rosa,Jr.,and O.Hipolito,Am.J.Phys.44,744(1976).[12]V.Bagnato and D.Kleppner,Phys.Rev.A44,7439(1991).[13]T.Haugset and H.Haugerud,Phys.Rev.A57,3809(1998).[14]G.-L.Ingold and mbrecht,Eur.Phys.J.D1,29(1998).[15]W.J.Mullin,J.Low Temp.Phys.110,167(1998)[16]W.Deng and P.M.Hui,Solid State Commun.104,729(1997).[17]H.Monien,M.Linn,and N.Elstner,preprint,cond-mat/9711178.[18]F.Brosens,J.T.Devreese,and L.F.Lemmens,SolidState Commun.100,123(1996).[19]S.Pearson,T.Pang,and C.Chen,Phys.Rev.A58,1485(1998).[20]S.Giorgini,L.P.Pitaevskii,and S.Stringari,Phys.Rev.A,54,4633(1996);J.Low Temp.Phys.109,309(1997).[21]A.Minguzzi,S.Conti,and M.P.Tosi,J.Phys.Cond.Matter9,L33(1997).[22]M.Bayindir and B.Tanatar,Phys.Rev.A58,3134(1998).[23]S.Heinrichs and W.J.Mullin,preprint,cond-mat/9807331.[24]E.P.Gross,Nuovo Cimento20,454(1961).[25]L.P.Pitaevskii,Zh.Eksp.Teor.Fiz.40,646(1961)[Sov.Phys.JETP13,451(1961)].[26]M.Olshanii,Phys.Rev.Lett.81,938(1998).[27]R.J.Dodd,K.Burnett,M.Edwards,and C.W.Clark,Acta Phys.Pol.A93,45(1998).[28]M.Naraschewski and D.M.Stamper-Kurn,Phys.Rev.A58,2423(1998).[29]V.Bagnato,D.E.Pritchard,and D.Kleppner,Phys.Rev.A35,4354(1987).[30]G.Baym and C.J.Pethick,Phys.Rev.Lett.76,6(1996).[31]S.Grossmann and M.Holthaus,Z.Naturforsch.50a,323(1995).[32]F.Dalfovo,L.P.Pitaevskii,and S.Stringari,Phys.Rev.A54,4213(1996).[33]A.Smerzi,S.Fantoni,S.Giovanazzi,and S.R.Shenoy,Phys.Rev.Lett.79,4950(1997).[34]I.Zapata,F.Sols,and A.J.Leggett,Phys.Rev.A57,28(1998).。

In the July 14, 1995 issue of Science magazine, researchers from JILA reported achieving a temperature far lower than had ever been produced before and creating an entirely new state of matter predicted decades ago by Albert Einstein and Indian physicist Satyendra Nath Bose. Cooling rubidium atoms to less than 170 billionths of a degree above absolute zero caused the individual atoms to condense into a "superatom" behaving as a single ent ity. The graphic shows three-dimensional successive snap shots in time in which the atoms condensed from less dense red, yellow and green areas into very dense blue to white areas. JILA is jointly operated by NIST and the University of Colorado at Boulder.A Bose–Einstein condensate (BEC) is a state of matter of a dilute gas of weakly interacting bosons confined in an external potential and cooled to temperatures very near to absolute zero (0 K, −273.15 °C, or −459.67 °F). Under such conditions, a large fraction of the bosons occupy the lowest quantum state of the external potential, and all wave functions overlap each other, at which point quantum effects become apparent on a macroscopic scale.This state of matter was first predicted by Satyendra Nath Bose and Albert Einstein in 1924–25. Bose first sent a paper to Einstein on the quantum statistics of light quanta (now called photons). Einstein was impressed, translated the paper himself from English to German and submitted it for Bose to the Zeitschrift für Physik which published it. Einstein then extended Bose's ideas to material particles (or matter) in two other papers.[1]Seventy years later, the first gaseous condensate was produced by Eric Cornell and Carl Wieman in 1995 at the University of Colorado at Boulder NIST-JILA lab, using a gas of rubidium atoms cooled to 170 nanokelvin (nK) [2](1.7×10−7 K). For their acheivments Cornell, Wieman, and Wolfgang Ketterle at MIT received the 2001 Nobel Prize in Physics. [3]TheoryThe slowing of atoms by use of cooling apparatus produces a singular quantum state known as a Bose condensate or Bose–Einstein condensate. This phenomenon was predicted in 1925 by generalizing Satyendra Nath Bose's work on the statistical mechanics of (massless) photons to (massive) atoms. (The Einstein manuscript, once believed to be lost, was found in a library at Leiden University in 2005.[4]) The result of the efforts of Bose and Einstein is the concept of a Bose gas, governed by Bose–Einstein statistics, which describes the statistical distribution of identical particles with integer spin, now known as bosons. Bosonic particles, which include the photon as well as atoms such as helium-4, are allowed to share quantum states with each other. Einstein demonstrated that cooling bosonic atoms to a very low temperature would cause them to fall (or "condense") into the lowest accessible quantum state, resulting in a new form of matter.This transition occurs below a critical temperature, which for a uniform three-dimensional gas consisting of non-interacting particles with no apparent internal degrees of freedom is given by:where:is t he critical temperature,is t he particle density,is t he mass per boson,is t he reduced Planck constant,is t he Boltzmann constant, andthe Riemann zeta function; (sequenceA078434 in OEIS)Einstein's argumentConsider a collection of N noninteracting particles which can each be in one of two quantum states, and . If the two states are equal in energy, each different configuration is equally likely.If we can tell which particle is which, there are 2N differentconfigurations, since each particle can be in or independently. In almost all the configurations, about half the particles are in and the other half in . The balance is a statistical effect, the number of configurations is largest when the particles are divided equally.If the particles are indistinguishable, however, there are only N+1 different configurations. If there are K particles in state , there are N-K particles in state . Whether any particular particle is in state or in state cannot be determined, so each value of K determines a unique quantum state for the whole system. If all these states are equally likely, there is no statistical spreading out; it is just as likely for all the particles to sit in as for the particles to be split half and half.Suppose now that the energy of state is slightly greater than the energy of state by an amount E. At temperature T, a particle will have a lesser probability to be in state by exp(-E/T). In the distinguishable case, the particle distribution will be biased slightly towards state and the distribution will be slightly different from half and half. But in the indistinguishable case, since there is no statistical pressure toward equal numbers, the most likely outcome is that most of the particles will collapse into state .In the distinguishable case, for large N, the fraction in state can be computed. It is the same as coin flipping with a coin which has probability p = exp(-E/T) to land tails. The fraction of heads is 1/(1+p), which is a smooth function of p, of the energy.In the indistinguishable case, each value of K is a single state, which has its own separate Boltzmann probability. So the probability distribution is exponential:For large N, the normalization constant C is (1-p). The expected total number of particles which are not in the lowest energy state, in the limitthat , is equal to . It doesn't grow when N islarge, it just approaches a constant. This will be a negligible fraction of the total number of particles. So a collection of enough bose particles in thermal equilibrium will mostly be in the ground state, with only a few in any excited state, no matter how small the energy difference.Consider now a gas of particles, which can be in different momentum states labelled . If the number of particles is less than the number of thermally accessible states, for high temperatures and low densities, the particles will all be in different states. In this limit the gas is classical. As the density increases or the temperature decreases, the number of accessible states per particle becomes smaller, and at some point more particles will be forced into a single state than the maximum allowed for that state by statistical weighting. From this point on, any extra particle added will go into the ground state.To calculate the transition temperature at any density, integrate over all momentum states the expression for maximum number of excited particles p/(1-p):When the integral is evaluated with the factors of k B and ℏ restored by dimensional analysis, it gives the critical temperature formula of the preceding section. Therefore, this integral defines the critical temperature and particle number corresponding to the conditions of zero chemical potential (μ = 0 in the Bose–Einstein statistics distribution).Gross–Pitaevskii equationMain article: Gross–Pitaevskii equationThe state of the BEC can be described by the wavefunction of the condensate . For a system of this nature, is interpreted as the particledensity, so the total number of atoms isProvided essentially all atoms are in the condensate (that is, have condensed to the ground state), and treating the bosons using mean fieldtheory, the energy (E) associated with the state is:Minimising this energy with respect to infinitesimal variations in ,and holding the number of atoms constant, yields the Gross-Pitaevski equation (GPE) (also a non-linear Schrödinger equation):where:is the mass of the bosons,is the external potential,is representative of the inter-particle interactions.The GPE provides a good description of the behavior of BEC's and is thus often applied for theoretical analysis.DiscoveryIn 1938, Pyotr Kapitsa, John Allen and Don Misener discovered that helium-4 became a new kind of fluid, now known as a superfluid, at temperatures less than 2.17 K (the lambda point). Superfluid helium has many unusual properties, including zero viscosity (the ability to flow without dissipating energy) and the existence of quantized vortices. It was quickly realized that the superfluidity was due to partialBose–Einstein condensation of the liquid. In fact, many of the properties of superfluid helium also appear in the gaseous Bose–Einstein condensates created by Cornell, Wieman and Ketterle (see below). Superfluid helium-4 is a liquid rather than a gas, which means that the interactions between the atoms are relatively strong; the original theory of Bose–Einstein condensation must be heavily modified in order to describe it. Bose–Einstein condensation remains, however, fundamental to the superfluid properties of helium-4. Note that helium-3, consisting of fermions instead of bosons, also enters a superfluid phase at lowtemperature, which can be explained by the formation of bosonic Cooper pairs of two atoms each (see also fermionic condensate).The first "pure" Bose–Einstein condensate was created by Eric Cornell, Carl Wieman, and co-workers at JILA on June 5, 1995. They did this by cooling a dilute vapor consisting of approximately two thousand rubidium-87 atoms to below 170 nK using a combination of laser cooling (a technique that won its inventors Steven Chu, Claude Cohen-Tannoudji, and William D. Phillips the 1997 Nobel Prize in Physics) and magnetic evaporative cooling. About four months later, an independent effort led by Wolfgang Ketterle at MIT created a condensate made of sodium-23. Ketterle's condensate had about a hundred times more atoms, allowing him to obtain several important results such as the observation of quantum mechanical interference between two different condensates. Cornell, Wieman and Ketterle won the 2001 Nobel Prize in Physics for their achievement.[5]The Bose–Einstein condensation also applies to quasiparticles in solids.A magnon in an antiferromagnet carries spin 1 and thus obeysBose–Einstein statistics. The density of magnons is controlled by an external magnetic field, which plays the role of the magnon chemical potential. This technique provides access to a wide range of boson densities from the limit of a dilute Bose gas to that of a strongly interacting Bose liquid. A magnetic ordering observed at the point of condensation is the analog of superfluidity. In 1999 Bose condensation of magnons was demonstrated in the antiferromagnet TlCuCl3.[6] The condensation was observed at temperatures as large as 14 K. Such a high transition temperature (relative to that of atomic gases) is due to the greater density achievable with magnons and the smaller mass (roughly equal to the mass of an electron). In 2006, condensation of magnons in ferromagnets was even shown at room temperature,[7]where the authors used pumping techniques.Velocity-distribution data graph（速度分布数据图）Velocity-distribution data of a gas of rubidium atoms, confirming the discovery of a new phase of matter, the Bose–Einstein condensate. Left: just before the appearance of a Bose–Einstein condensate. Center: just after the appearance of the condensate. Right: after further evaporation, leaving a sample of nearly pure condensate.In the image accompanying this article, the velocity-distribution data indicates the formation of a Bose–Einstein condensate out of a gas of rubidium atoms. The false colors indicate the number of atoms at each velocity, with red being the fewest and white being the most. The areas appearing white and light blue are at the lowest velocities. The peak is not infinitely narrow because of the Heisenberg uncertainty principle: since the atoms are trapped in a particular region of space, their velocity distribution necessarily possesses a certain minimum width. This width is given by the curvature of the magnetic trapping potential in the given direction. More tightly confined directions have bigger widths in the ballistic velocity distribution. This anisotropy of the peak on the right is a purely quantum-mechanical effect and does not exist in the thermal distribution on the left. This famous graph served as the cover-design for 1999 textbook Thermal Physics by Ralph Baierlein.[8]Vortices（涡旋）As in many other systems, vortices can exist in BECs. These can be created, for example, by 'stirring' the condensate with lasers, or rotating the confining trap. The vortex created will be a quantum vortex. Thesephenomena are allowed for by the non-linear term in the GPE. As the vortices must have quantised angular momentum, the wavefunction will be of the form where ρ,z and θ are as in thecylindrical coordinate system, and is the angular number. To determine φ(ρ,z), the energy of must be minimised, according to the constraint . This is usually done computationally, however in a uniform medium the analytic formwhere:is d ensity far from the vortex,is h ealing length of the condensate.demonstrates the correct behavior, and is a good approximation.A singly-charged vortex () is in the ground state, with its energy εv given bywhere:is t he farthest distance from the vortex considered.(To obtain an energy which is well defined it is necessary to include this boundary b.)For multiply-charged vortices () the energy is approximated bywhich is greater than that of singly-charged vortices, indicating that these multiply-charged vortices are unstable to decay. Research has,however, indicated they are metastable states, so may have relatively long lifetimes.Closely related to the creation of vortices in BECs is the generation of so-called dark solitons in one-dimensional BECs. These topological objects feature a phase gradient across their nodal plane, which stabilizes their shape even in propagation and interaction. Although solitons carry no charge and are thus prone to decay, relativelylong-lived dark solitons have been produced and studied extensively.[9] Unusual characteristicsFurther experimentation by the JILA team in 2000 uncovered a hitherto unknown property of Bose–Einstein condensates. Cornell, Wieman, and their coworkers originally used rubidium-87, an isotope whose atoms naturally repel each other, making a more stable condensate. The JILA team instrumentation now had better control over the condensate so experimentation was made on naturally attracting atoms of another rubidium isotope, rubidium-85 (having negative atom-atom scattering length). Through a process called Feshbach resonance involving a sweep of the magnetic field causing spin flip collisions, the JILA researchers lowered the characteristic, discrete energies at which the rubidium atoms bond into molecules, making their Rb-85 atoms repulsive and creating a stable condensate. The reversible flip from attraction to repulsion stems from quantum interference among condensate atoms which behave as waves.When the scientists raised the magnetic field strength still further, the condensate suddenly reverted back to attraction, imploded and shrank beyond detection, and then exploded, blowing off about two-thirds of its 10,000 or so atoms. About half of the atoms in the condensate seemed to have disappeared from the experiment altogether, not being seen either in the cold remnant or the expanding gas cloud.[10] Carl Wieman explained that under current atomic theory this characteristic of Bose–Einstein condensate could not be explained because the energy state of an atom near absolute zero should not be enough to cause an implosion; however, subsequent mean field theories have been proposed to explain it.Because supernova explosions are also preceded by an implosion, the explosion of a collapsing Bose–Einstein condensate was named "bosenova", a pun on the musical style bossa nova.The atoms that seem to have disappeared almost certainly still exist in some form, just not in a form that could be accounted for in that experiment. Most likely they formed molecules consisting of two bonded rubidium atoms.The energy gained by making this transition imparts a velocity sufficient for them to leave the trap without being detected.Current researchCompared to more commonly encountered states of matter, Bose–Einstein condensates are extremely fragile. The slightest interaction with the outside world can be enough to warm them past the condensation threshold, eliminating their interesting properties and forming a normal gas. It is likely to be some time before any practical applications are developed.Nevertheless, they have proven useful in exploring a wide range of questions in fundamental physics, and the years since the initial discoveries by the JILA and MIT groups have seen an explosion in experimental and theoretical activity. Examples include experiments that have demonstrated interference between condensates due to wave-particle duality,[11]the study of superfluidity and quantized vortices,[12]and the slowing of light pulses to very low speeds using electromagnetically induced transparency.[13]Vortices in Bose–Einstein condensates are also currently the subject of analogue gravity research, studying the possibility of modeling black holes and their related phenomena in such environments in the lab. Experimentalists have also realized "optical lattices", where the interference pattern from overlapping lasers provides a periodic potential for the condensate. These have been used to explore the transition between a superfluid and a Mott insulator,[14] and may be useful in studying Bose–Einstein condensation in fewer than three dimensions, for example the Tonks-Girardeau gas.Bose–Einstein condensates composed of a wide range of isotopes have been produced.[15]Related experiments in cooling fermions rather than bosons to extremely low temperatures have created degenerate gases, where the atoms do not congregate in a single state due to the Pauli exclusion principle. To exhibit Bose–Einstein condensation, the fermions must "pair up" to form compound particles (e.g. molecules or Cooper pairs) that are bosons. The first molecular Bose–Einstein condensates were created in November 2003 by the groups of Rudolf Grimm at the University of Innsbruck, Deborah S. Jin at the University of Colorado at Boulder and Wolfgang Ketterle at MIT. Jin quickly went on to create the first fermionic condensate composed of Cooper pairs.[16]In 1999, Danish physicist Lene Vestergaard Hau led a team from Harvard University which succeeded in slowing a beam of light to about 17 metresper second and, in 2001, was able to momentarily stop a beam. She was able to achieve this by using a superfluid. Hau and her associates at Harvard University have since successfully transformed light into matter and back into light using Bose–Einstein condensates: details of the experiment are discussed in an article in the journal Nature, 8 February 2007.[17]Subtleties(微妙之处)Up to 2004, using the above-mentioned "ultralow temperatures",Bose–Einstein condensates had been obtained for a multitude of isotopes involving mainly alkaline and alkaline earth atoms (7Li, 23Na, 41K, 52Cr,84St, 85Rb, 87Rb, 133Cs and 174Yb). Not astonishingly, condensation research was finally successful even with hydrogen, although with the aid of special methods. In contrast, the superfluid state of the bosonic 4He at temperatures below the temperature of 2.17 K is not a good example of Bose–Einstein condensation, because the interaction between the 4He bosons is simply too strong, so that at zero temperature, contrary to Bose–Einstein theory, only 8% rather than 100% of the atoms are in the ground state. Even the fact that the above-mentioned alkaline gases show bosonic, rather than fermionic behaviour, as solid state physicists or chemists would expect, is based on a subtle interplay of electronic and nuclear spins: at ultralow temperatures and corresponding excitation energies, the half-integer (in units of ) total spin of the electronic shell and the also half-integer total spin of the nucleus of the atom are coupled by the (very weak) hyperfine interaction to the integer (!) total spin of the atom. Only the fact that this last-mentioned total spin is integral causes the ultralow-temperature behaviour of the atom to be bosonic, whereas the "chemistry" of the systems at room temperature is determined by the electronic properties, i.e. is essentially fermionic, since at room temperature thermal excitations have typical energies which are much higher than the hyperfine values. (Here one should remember the spin-statistics theorem of Wolfgang Pauli, which states that half-integer spins lead to fermionic behaviour, e.g., the Pauli exclusion principle forbidding that more than two electrons possess the same energy, whereas integer spins lead to bosonic behaviour, e.g., condensation of identical bosonic particles in a common ground state.)The ultralow temperature requirement of Bose–Einstein condensates of alkali metals does not generalize to all types of Bose–Einstein condensates. In 2006, physicists under S. Demokritov in Münster, Germany,[18]found Bose–Einstein condensation of magnons (i.e. quantized spinwaves) at room temperature, admittedly by the application of pump processes。

Chapter19Bose-Einstein CondensationAbstract Bose-Einstein condensation(BEC)refers to a prediction of quantum sta-tistical mechanics(Bose[1],Einstein[2])where an ideal gas of identical bosons undergoes a phase transition when the thermal de Broglie wavelength exceeds the mean spacing between the particles.Under these conditions,bosons are stimulated by the presence of other bosons in the lowest energy state to occupy that state as well,resulting in a macroscopic occupation of a single quantum state.The con-densate that forms constitutes a macroscopic quantum-mechanical object.BEC was first observed in1995,seventy years after the initial predictions,and resulted in the award of2001Nobel Prize in Physics to Cornell,Ketterle and Weiman.The exper-imental observation of BEC was achieved in a dilute gas of alkali atoms in a mag-netic trap.Thefirst experiments used87Rb atoms[3],23Na[4],7Li[5],and H[6] more recently metastable He has been condensed[7].The list of BEC atoms now includes molecular systems such as Rb2[8],Li2[9]and Cs2[10].In order to cool the atoms to the required temperature(∼200nK)and densities(1013–1014cm−3) for the observation of BEC a combination of optical cooling and evaporative cooling were employed.Early experiments used magnetic traps but now optical dipole traps are also common.Condensates containing up to5×109atoms have been achieved for atoms with a positive scattering length(repulsive interaction),but small con-densates have also been achieved with only a few hundred atoms.In recent years Fermi degenerate gases have been produced[11],but we will not discuss these in this chapter.BECs are now routinely produced in dozens of laboratories around the world. They have provided a wonderful test bed for condensed matter physics with stunning experimental demonstrations of,among other things,interference between conden-sates,superfluidity and vortices.More recently they have been used to create opti-cally nonlinear media to demonstrate electromagnetically induced transparency and neutral atom arrays in an optical lattice via a Mott insulator transition.Many experiments on BECs are well described by a semiclassical theory dis-cussed below.Typically these involve condensates with a large number of atoms, and in some ways are analogous to describing a laser in terms of a semiclassi-cal meanfield.More recent experiments however have begun to probe quantum39739819Bose-Einstein Condensation properties of the condensate,and are related to the fundamental discreteness of the field and nonlinear quantum dynamics.In this chapter,we discuss some of these quantum properties of the condensate.We shall make use of“few mode”approxi-mations which treat only essential condensate modes and ignore all noncondensate modes.This enables us to use techniques developed for treating quantum optical systems described in earlier chapters of this book.19.1Hamiltonian:Binary Collision ModelThe effects of interparticle interactions are of fundamental importance in the study of dilute–gas Bose–Einstein condensates.Although the actual interaction potential between atoms is typically very complex,the regime of operation of current exper-iments is such that interactions can in fact be treated very accurately with a much–simplified model.In particular,at very low temperature the de Broglie wavelengths of the atoms are very large compared to the range of the interatomic potential.This, together with the fact that the density and energy of the atoms are so low that they rarely approach each other very closely,means that atom–atom interactions are ef-fectively weak and dominated by(elastic)s–wave scattering.It follows also that to a good approximation one need only consider binary collisions(i.e.,three–body processes can be neglected)in the theoretical model.The s–wave scattering is characterised by the s–wave scattering length,a,the sign of which depends sensitively on the precise details of the interatomic potential [a>0(a<0)for repulsive(attractive)interactions].Given the conditions described above,the interaction potential can be approximated byU(r−r )=U0δ(r−r ),(19.1) (i.e.,a hard sphere potential)with U0the interaction“strength,”given byU0=4π¯h2am,(19.2)and the Hamiltonian for the system of weakly interacting bosons in an external potential,V trap(r),can be written in the second quantised form asˆH=d3rˆΨ†(r)−¯h22m∇2+V trap(r)ˆΨ(r)+12d3rd3r ˆΨ†(r)ˆΨ†(r )U(r−r )ˆΨ(r )ˆΨ(r)(19.3)whereˆΨ(r)andˆΨ†(r)are the bosonfield operators that annihilate or create a par-ticle at the position r,respectively.19.2Mean–Field Theory —Gross-Pitaevskii Equation 399To put a quantitative estimate on the applicability of the model,if ρis the density of bosons,then a necessary condition is that a 3ρ 1(for a >0).This condition is indeed satisfied in the alkali gas BEC experiments [3,4],where achieved densities of the order of 1012−1013cm −3correspond to a 3ρ 10−5−10−6.19.2Mean–Field Theory —Gross-Pitaevskii EquationThe Heisenberg equation of motion for ˆΨ(r )is derived as i¯h ∂ˆΨ(r ,t )∂t = −¯h 22m ∇2+V trap (r ) ˆΨ(r ,t )+U 0ˆΨ†(r ,t )ˆΨ(r ,t )ˆΨ(r ,t ),(19.4)which cannot in general be solved.In the mean–field approach,however,the expec-tation value of (19.4)is taken and the field operator decomposed asˆΨ(r ,t )=Ψ(r ,t )+˜Ψ(r ,t ),(19.5)where Ψ(r ,t )= ˆΨ(r ,t ) is the “condensate wave function”and ˜Ψ(r )describes quantum and thermal fluctuations around this mean value.The quantity Ψ(r ,t )is in fact a classical field possessing a well–defined phase,reflecting a broken gauge sym-metry associated with the condensation process.The expectation value of ˜Ψ(r ,t )is zero and,in the mean–field theory,its effects are assumed to be small,amounting to the assumption of the thermodynamic limit,where the number of particles tends to infinity while the density is held fixed.For the effects of ˜Ψ(r )to be negligibly small in the equation for Ψ(r )also amounts to an assumption of zero temperature (i.e.,pure condensate).Given that this is so,and using the normalisationd 3r |Ψ(r ,t )|2=1,(19.6)one is lead to the nonlinear Schr¨o dinger equation,or “Gross–Pitaevskii equation”(GP equation),for the condensate wave function Ψ(r ,t )[13],i¯h ∂Ψ(r ,t )∂t = −¯h 22m ∇2+V trap (r )+NU 0|Ψ(r ,t )|2 Ψ(r ,t ),(19.7)where N is the mean number of particles in the condensate.The nonlinear interaction term (or mean–field pseudo–potential)is proportional to the number of atoms in the condensate and to the s –wave scattering length through the parameter U 0.A stationary solution forthe condensate wavefunction may be found by substi-tuting ψ(r ,t )=exp −i μt ¯h ψ(r )into (19.7)(where μis the chemical potential of the condensate).This yields the time independent equation,40019Bose-Einstein Condensation−¯h2 2m ∇2+V trap(r)+NU0|ψ(r)|2ψ(r)=μψ(r).(19.8)The GP equation has proved most successful in describing many of the meanfield properties of the condensate.The reader is referred to the review articles listed in further reading for a comprehensive list of references.In this chapter we shall focus on the quantum properties of the condensate and to facilitate our investigations we shall go to a single mode model.19.3Single Mode ApproximationThe study of the quantum statistical properties of the condensate(at T=0)can be reduced to a relatively simple model by using a mode expansion and subsequent truncation to just a single mode(the“condensate mode”).In particular,one writes the Heisenberg atomicfield annihilation operator as a mode expansion over single–particle states,ˆΨ(r,t)=∑αaα(t)ψα(r)exp−iμαt/¯h=a0(t)ψ0(r)exp−iμ0t/¯h+˜Ψ(r,t),(19.9) where[aα(t),a†β(t)]=δαβand{ψα(r)}are a complete orthonormal basis set and {μα}the corresponding eigenvalues.Thefirst term in the second line of(19.9)acts only on the condensate state vector,withψ0(r)chosen as a solution of the station-ary GP equation(19.8)(with chemical potentialμ0and mean number of condensate atoms N).The second term,˜Ψ(r,t),accounts for non–condensate atoms.Substitut-ing this mode expansion into the HamiltonianˆH=d3rˆΨ†(r)−¯h22m∇2+V trap(r)ˆΨ(r)+(U0/2)d3rˆΨ†(r)ˆΨ†(r)ˆΨ(r)ˆΨ(r),(19.10)and retaining only condensate terms,one arrives at the single–mode effective Hamil-tonianˆH=¯h˜ω0a †a0+¯hκa†0a†0a0a0,(19.11)where¯h˜ω0=d3rψ∗0(r)−¯h22m∇2+V trap(r)ψ0(r),(19.12)and¯hκ=U02d3r|ψ0(r)|4.(19.13)19.5Quantum Phase Diffusion:Collapses and Revivals of the Condensate Phase401 We have assumed that the state is prepared slowly,with damping and pumping rates vanishingly small compared to the trap frequencies and collision rates.This means that the condensate remains in thermodynamic equilibrium throughout its prepara-tion.Finally,the atom number distribution is assumed to be sufficiently narrow that the parameters˜ω0andκ,which of course depend on the atom number,can be re-garded as constants(evaluated at the mean atom number).In practice,this proves to be a very good approximation.19.4Quantum State of the CondensateA Bose-Einstein condensate(BEC)is often viewed as a coherent state of the atomic field with a definite phase.The Hamiltonian for the atomicfield is independent of the condensate phase(see Exercise19.1)so it is often convenient to invoke a symmetry breaking Bogoliubovfield to select a particular phase.In addition,a coherent state implies a superposition of number states,whereas in a single trap experiment there is afixed number of atoms in the trap(even if we are ignorant of that number)and the state of a simple trapped condensate must be a number state(or,more precisely, a mixture of number states as we do not know the number in the trap from one preparation to the next).These problems may be bypassed by considering a system of two condensates for which the total number of atoms N isfixed.Then,a general state of the system is a superposition of number difference states of the form,|ψ =N∑k=0c k|k,N−k (19.14)As we have a well defined superposition state,we can legitimately consider the relative phase of the two condensates which is a Hermitian observable.We describe in Sect.19.6how a particular relative phase is established due to the measurement process.The identification of the condensate state as a coherent state must be modified in the presence of collisions except in the case of very strong damping.19.5Quantum Phase Diffusion:Collapsesand Revivals of the Condensate PhaseThe macroscopic wavefunction for the condensate for a relatively strong number of atoms will exhibit collapses and revivals arising from the quantum evolution of an initial state with a spread in atom number[21].The initial collapse has been described as quantum phase diffusion[20].The origins of the collapses and revivals may be seen straightforwardly from the single–mode model.From the Hamiltonian40219Bose-Einstein CondensationˆH =¯h ˜ω0a †0a 0+¯h κa †0a †0a 0a 0,(19.15)the Heisenberg equation of motion for the condensate mode operator follows as˙a 0(t )=−i ¯h [a 0,H ]=−i ˜ω0a 0+2κa †0a 0a 0 ,(19.16)for which a solution can be written in the form a 0(t )=exp −i ˜ω0+2κa †0a 0 t a 0(0).(19.17)Writing the initial state of the condensate,|i ,as a superposition of number states,|i =∑n c n |n ,(19.18)the expectation value i |a 0(t )|i is given byi |a 0(t )|i =∑n c ∗n −1c n √n exp {−i [˜ω0+2κ(n −1)]t }=∑nc ∗n −1c n √n exp −i μt ¯h exp {−2i κ(n −N )t },(19.19)where the relationship μ=¯h ˜ω0+2¯h κ(N −1),(19.20)has been used [this expression for μuses the approximation n 2 =N 2+(Δn )2≈N 2].The factor exp (−i μt /¯h )describes the deterministic motion of the condensate mode in phase space and can be removed by transforming to a rotating frame of reference,allowing one to writei |a 0(t )|i =∑nc ∗n −1c n √n {cos [2κ(n −N )t ]−isin [2κ(n −N )t ]}.(19.21)This expression consists of a weighted sum of trigonometric functions with different frequencies.With time,these functions alternately “dephase”and “rephase,”giving rise to collapses and revivals,respectively,in analogy with the behaviour of the Jaynes–Cummings Model of the interaction of a two–level atom with a single elec-tromagnetic field mode described in Sect.10.2.The period of the revivals follows di-rectly from (19.21)as T =π/κ.The collapse time can be derived by considering the spread of frequencies for particle numbers between n =N +(Δn )and n =N −(Δn ),which yields (ΔΩ)=2κ(Δn );from this one estimates t coll 2π/(ΔΩ)=T /(Δn ),as before.From the expression t coll T /(Δn ),it follows that the time taken for collapse depends on the statistics of the condensate;in particular,on the “width”of the initial distribution.This dependence is illustrated in Fig.19.1,where the real part of a 0(t )19.5Quantum Phase Diffusion:Collapses and Revivals of the Condensate Phase403Fig.19.1The real part ofthe condensate amplitudeversus time,Re { a 0(t ) }for an amplitude–squeezed state,(a )and a coherent state (b )with the same mean numberof atoms,N =250.20.40.60.81-11234560b a is plotted as a function of time for two different initial states:(a)an amplitude–squeezed state,(b)a coherent state.The mean number of atoms is chosen in each case to be N =25.The timescales of the collapses show clear differences;the more strongly number–squeezed the state is,the longer its collapse time.The revival times,how-ever,are independent of the degree of number squeezing and depend only on the interaction parameter,κ.For example,a condensate of Rb 2,000atoms with the ω/2π=60Hz,has revival time of approximately 8s,which lies within the typical lifetime of the experimental condensate (10–20s).One can examine this phenomenon in the context of the interference between a pair of condensates and indeed one finds that the visibility of the interference pat-tern also exhibits collapses and revivals,offering an alternative means of detecting this effect.To see this,consider,as above,that atoms are released from two conden-sates with momenta k 1and k 2respectively.Collisions within each condensate are described by the Hamiltonian (neglecting cross–collisions)ˆH =¯h κ a †1a 1 2+ a †2a 22 ,(19.22)from which the intensity at the detector follows asI (x ,t )=I 0 [a †1(t )exp i k 1x +a †2(t )expi k 2x ][a 1(t )exp −i k 1x +a 2(t )exp −i k 2x ] =I 0 a †1a 1 + a †2a 2+ a †1exp 2i a †1a 1−a †2a 2 κt a 2 exp −i φ(x )+h .c . ,(19.23)where φ(x )=(k 2−k 1)x .If one assumes that each condensate is initially in a coherent state of amplitude |α|,with a relative phase φbetween the two condensates,i.e.,assuming that|ϕ(t =0) =|α |αe −i φ ,(19.24)40419Bose-Einstein Condensation then one obtains for the intensityI(x,t)=I0|α|221+exp2|α|2(cos(2κt)−1)cos[φ(x)−φ].(19.25)From this expression,it is clear that the visibility of the interference pattern under-goes collapses and revivals with a period equal toπ/κ.For short times t 1/2κ, this can be written asI(x,t)=I0|α|221+exp−|α|2κ2t2,(19.26)from which the collapse time can be identified as t coll=1/κ|α|.An experimental demonstration of the collapse and revival of a condensate was done by the group of Bloch in2002[12].In the experiment coherent states of87Rb atoms were prepared in a three dimensional optical lattice where the tunneling is larger than the on-site repulsion.The condensates in each well were phase coherent with constant relative phases between the sites,and the number distribution in each well is close to Poisonnian.As the optical dipole potential is increased the depth of the potential wells increases and the inter-well tunneling decreases producing a sub-Poisson number distribution in each well due to the repulsive interaction between the atoms.After preparing the states in each well,the well depth is rapidly increased to create isolated potential wells.The nonlinear interaction of(19.15)then determines the dynamics in each well.After some time interval,the hold time,the condensate is released from the trap and the resulting interference pattern is imaged.As the meanfield amplitude in each well undergoes a collapse the resulting interference pattern visibility decreases.However as the meanfield revives,the visibility of the interference pattern also revives.The experimental results are shown in Fig.19.2.Fig.19.2The interference pattern imaged from the released condensate after different hold times. In(d)the interference fringes have entirely vanished indicating a complete collapse of the am-plitude of the condensate.In(g),the wait time is now close to the complete revival time for the coherent amplitude and the fringe pattern is restored.From Fig.2of[12]19.6Interference of Two Bose–Einstein Condensates and Measurement–Induced Phase405 19.6Interference of Two Bose–Einstein Condensatesand Measurement–Induced PhaseThe standard approach to a Bose–Einstein condensate assumes that it exhibits a well–defined amplitude,which unavoidably introduces the condensate phase.Is this phase just a formal construct,not relevant to any real measurement,or can one ac-tually observe something in an experiment?Since one needs a phase reference to observe a phase,two options are available for investigation of the above question. One could compare the condensate phase to itself at a different time,thereby ex-amining the condensate phase dynamics,or one could compare the phases of two distinct condensates.This second option has been studied by a number of groups, pioneered by the work of Javanainen and Yoo[23]who consider a pair of statisti-cally independent,physically–separated condensates allowed to drop and,by virtue of their horizontal motion,overlap as they reach the surface of an atomic detec-tor.The essential result of the analysis is that,even though no phase information is initially present(the initial condensates may,for example,be in number states),an interference pattern may be formed and a relative phase established as a result of the measurement.This result may be regarded as a constructive example of sponta-neous symmetry breaking.Every particular measurement produces a certain relative phase between the condensates;however,this phase is random,so that the symme-try of the system,being broken in a single measurement,is restored if an ensemble of measurements is considered.The physical configuration we have just described and the predicted interference between two overlapping condensates was realised in a beautiful experiment per-formed by Andrews et al.[18]at MIT.The observed fringe pattern is shown in Fig.19.8.19.6.1Interference of Two Condensates Initially in Number States To outline this effect,we follow the working of Javanainen and Yoo[23]and consider two condensates made to overlap at the surface of an atom detector.The condensates each contain N/2(noninteracting)atoms of momenta k1and k2,respec-tively,and in the detection region the appropriatefield operator isˆψ(x)=1√2a1+a2exp iφ(x),(19.27)whereφ(x)=(k2−k1)x and a1and a2are the atom annihilation operators for the first and second condensate,respectively.For simplicity,the momenta are set to±π, so thatφ(x)=2πx.The initial state vector is represented simply by|ϕ(0) =|N/2,N/2 .(19.28)40619Bose-Einstein Condensation Assuming destructive measurement of atomic position,whereby none of the atoms interacts with the detector twice,a direct analogy can be drawn with the theory of absorptive photodetection and the joint counting rate R m for m atomic detections at positions {x 1,···,x m }and times {t 1,···,t m }can be defined as the normally–ordered averageR m (x 1,t 1,...,x m ,t m )=K m ˆψ†(x 1,t 1)···ˆψ†(x m ,t m )ˆψ(x m ,t m )···ˆψ(x 1,t 1) .(19.29)Here,K m is a constant that incorporates the sensitivity of the detectors,and R m =0if m >N ,i.e.,no more than N detections can occur.Further assuming that all atoms are in fact detected,the joint probability density for detecting m atoms at positions {x 1,···,x m }follows asp m (x 1,···,x m )=(N −m )!N ! ˆψ†(x 1)···ˆψ†(x m )ˆψ(x m )···ˆψ(x 1) (19.30)The conditional probability density ,which gives the probability of detecting an atom at the position x m given m −1previous detections at positions {x 1,···,x m −1},is defined as p (x m |x 1,···,x m −1)=p m (x 1,···,x m )p m −1(x 1,···,x m −1),(19.31)and offers a straightforward means of directly simulating a sequence of atom detections [23,24].This follows from the fact that,by virtue of the form for p m (x 1,···,x m ),the conditional probabilities can all be expressed in the simple formp (x m |x 1,···,x m −1)=1+βcos (2πx m +ϕ),(19.32)where βand ϕare parameters that depend on {x 1,···,x m −1}.The origin of this form can be seen from the action of each measurement on the previous result,ϕm |ˆψ†(x )ˆψ(x )|ϕm =(N −m )+2A cos [θ−φ(x )],(19.33)with A exp −i θ= ϕm |a †1a 2|ϕm .So,to simulate an experiment,one begins with the distribution p 1(x )=1,i.e.,one chooses the first random number (the position of the first atom detection),x 1,from a uniform distribution in the interval [0,1](obviously,before any measurements are made,there is no information about the phase or visibility of the interference).After this “measurement,”the state of the system is|ϕ1 =ˆψ(x 1)|ϕ0 = N /2 |(N /2)−1,N /2 +|N /2,(N /2)−1 expi φ(x 1) .(19.34)That is,one now has an entangled state containing phase information due to the fact that one does not know from which condensate the detected atom came.The corre-sponding conditional probability density for the second detection can be derived as19.6Interference of Two Bose–Einstein Condensates and Measurement–Induced Phase 407n u m b e r o f a t o m s n u m b e r o f a t o m s 8position Fig.19.3(a )Numerical simulation of 5,000atomic detections for N =10,000(circles).The solid curve is a least-squares fit using the function 1+βcos (2πx +ϕ).The free parameters are the visibility βand the phase ϕ.The detection positions are sorted into 50equally spaced bins.(b )Collisions included (κ=2γgiving a visibility of about one-half of the no collision case.From Wong et al.[24]40819Bose-Einstein Condensationp (x |x 1)=p 2(x 1,x )p 1(x 1)=1N −1 ˆψ†(x 1)ˆψ†(x )ˆψ(x )ˆψ(x 1) ˆψ†(x 1)ˆψ(x 1) (19.35)=12 1+N 2(N −1)cos [φ(x )−φ(x 1)] .(19.36)Hence,after just one measurement the visibility (for large N )is already close to 1/2,with the phase of the interference pattern dependent on the first measurement x 1.The second position,x 2,is chosen from the distribution (19.36).The conditional proba-bility p (x |x 1)has,of course,the form (19.32),with βand ϕtaking simple analytic forms.However,expressions for βand ϕbecome more complicated with increasing m ,and in practice the approach one takes is to simply calculate p (x |x 1,···,x m −1)numerically for two values of x [using the form (19.30)for p m (x 1,...,x m −1,x ),and noting that p m −1(x 1,...,x m −1)is simply a number already determined by the simu-lation]and then,using these values,solve for βand ϕ.This then defines exactly the distribution from which to choose x m .The results of simulations making use of the above procedure are shown in Figs 19.3–19.4.Figure 19.3shows a histogram of 5,000atom detections from condensates initially containing N /2=5,000atoms each with and without colli-sions.From a fit of the data to a function of the form 1+βcos (2πx +ϕ),the visibil-ity of the interference pattern,β,is calculated to be 1.The conditional probability distributions calculated before each detection contain what one can define as a con-000.10.20.30.40.50.60.70.80.91102030405060number of atoms decided 708090100x=0x=1x=2x=4x=6Fig.19.4Averaged conditional visibility as a function of the number of detected atoms.From Wong et al.[13]19.7Quantum Tunneling of a Two Component Condensate40900.51 1.520.500.5Θz ο00.51 1.520.500.5Θx ο(b)1,234elliptic saddle Fig.19.5Fixed point bifurcation diagram of the two mode semiclassical BEC dynamics.(a )z ∗,(b )x ∗.Solid line is stable while dashed line is unstable.ditional visibility .Following the value of this conditional visibility gives a quantita-tive measure of the buildup of the interference pattern as a function of the number of detections.The conditional visibility,averaged over many simulations,is shown as a function of the number of detections in Fig.19.4for N =200.One clearly sees the sudden increase to a value of approximately 0.5after the first detection,followed by a steady rise towards the value 1.0(in the absence of collisions)as each further detection provides more information about the phase of the interference pattern.One can also follow the evolution of the conditional phase contained within the conditional probability distribution.The final phase produced by each individual simulation is,of course,random but the trajectories are seen to stabilise about a particular value after approximately 50detections (for N =200).19.7Quantum Tunneling of a Two Component CondensateA two component condensate in a double well potential is a non trivial nonlinear dynamical model.Suppose the trapping potential in (19.3)is given byV (r )=b (x 2−q 20)2+12m ω2t (y 2+z 2)(19.37)where ωt is the trap frequency in the y –z plane.The potential has elliptic fixed points at r 1=+q 0x ,r 2=−q 0x near which the linearised motion is harmonic withfrequency ω0=q o (8b /m )1/2.For simplicity we set ωt =ω0and scale the length in units of r 0= ¯h /2m ω0,which is the position uncertainty in the harmonic oscillatorground state.The barrier height is B =(¯h ω/8)(q 0/r 0)2.We can justify a two mode expansion of the condensate field by assuming the potential parameters are chosen so that the two lowest single particle energy eigenstates are below the barrier,with41019Bose-Einstein Condensation the next highest energy eigenstate separated from the ground state doublet by a large gap.We will further assume that the interaction term is sufficiently weak that, near zero temperature,the condensate wave functions are well approximated by the single particle wave functions.The potential may be expanded around the two stablefixed points to quadratic orderV(r)=˜V(2)(r−r j)+...(19.38) where j=1,2and˜V(2)(r)=4bq2|r|2(19.39) We can now use as the local mode functions the single particle wave functions for harmonic oscillators ground states,with energy E0,localised in each well,u j(r)=−(−1)j(2πr20)3/4exp−14((x−q0)2+y2+z2)/r20(19.40)These states are almost orthogonal,with the deviation from orthogonality given by the overlap under the barrier,d3r u∗j(r)u k(r)=δj,k+(1−δj,k)ε(19.41) withε=e−12q20/r20.The localised states in(19.40)may be used to approximate the single particle energy(and parity)eigenstates asu±≈1√2[u1(r)±u2(r)](19.42)corresponding to the energy eigenvalues E±=E0±R withR=d3r u∗1(r)[V(r)−˜V(r−r1)]u2(r)(19.43)A localised state is thus an even or odd superposition of the two lowest energy eigenstates.Under time evolution the relative phase of the superposition can change sign after a time T=2π/Ω,the tunneling time,where the tunneling frequency is given byΩ=2R¯h=38ω0q20r2e−q20/2r20(19.44)We now make the two-mode approximation by expanding thefield operator asˆψ(r,t)=c1(t)u1(r)+c2(t)u2(r)(19.45) where。

Bose-Einstein condensationShihao LiBJTU ID#:13276013;UW ID#:20548261School of Science,Beijing Jiaotong University,Beijing,100044,ChinaJune1,20151What is BEC?To answer this question,it has to begin with the fermions and bosons.As is known,matters consist of atoms,atoms are made of protons,neutrons and electrons, and protons and neutrons are made of quarks.Also,there are photons and gluons that works for transferring interaction.All of these particles are microscopic particles and can be classified to two families,the fermion and the boson.A fermion is any particle characterized by Fermi–Dirac statistics.Particles with half-integer spin are fermions,including all quarks,leptons and electrons,as well as any composite particle made of an odd number of these,such as all baryons and many atoms and nuclei.As a consequence of the Pauli exclusion principle,two or more identical fermions cannot occupy the same quantum state at any given time.Differing from fermions,bosons obey Bose-Einstein statistics.Particles with integer spin are bosons,such as photons,gluons,W and Z bosons,the Higgs boson, and the still-theoretical graviton of quantum gravity.It also includes the composite particle made of even number of fermions,such as the nuclei with even number ofnucleons.An important characteristic of bosons is that their statistics do not restrict the number of them that occupy the same quantum state.For a single particle,when the temperature is at the absolute zero,0K,the particle is in the state of lowest energy,the ground state.Supposing that there are many particle,if they are fermions,there will be exactly one of them in the ground state;if they are bosons,most of them will be in the ground state,where these bosons share the same quantum states,and this state is called Bose-Einstein condensate (BEC).Bose–Einstein condensation(BEC)—the macroscopic groundstate accumulation of particles of a dilute gas with integer spin(bosons)at high density and low temperature very close to absolute zero.According to the knowledge of quantum mechanics,all microscopic particles have the wave-particle duality.For an atom in space,it can be expressed as well as a wave function.As is shown in the figure1.1,it tells the distribution but never exact position of atoms.Each distribution corresponds to the de Broglie wavelength of each atom.Lower the temperature is,lower the kinetic energy is,and longer the de Broglie wavelength is.p=mv=h/λ(Eq.1.1)When the distance of atoms is less than the de Broglie wavelength,the distribution of atoms are overlapped,like figure1.2.For the atoms of the same category,the overlapped distribution leads to a integral quantum state.If those atoms are bosons,each member will tend to a particular quantum state,and the whole atomsystem will become the BEC.In BEC,the physical property of all atoms is totally identical,and they are indistinguishable and like one independent atom.Figure1.1Wave functionsFigure1.2Overlapped wave functionWhat should be stressed is that the Bose–Einstein condensate is based on bosons, and BEC is a macroscopic quantum state.The first time people obtained BEC of gaseous rubidium atoms at170nK in lab was1995.Up to now,physicists have found the BEC of eight elements,most of which are alkali metals,calcium,and helium-4 atom.Always,the BEC of atom has some amazing properties which plays a vital role in the application of chip technology,precision measurement,and nano technology. What’s more,as a macroscopic quantum state,Bose–Einstein condensate gives a brand new research approach and field.2Bose and Einstein's papers were published in1924.Why does it take so long before it can be observed experimentally in atoms in1995?The condition of obtaining the BEC is daunting in1924.On the one hand,the temperature has to approach the absolute zero indefinitely;on the other hand,the aimed sample atoms should have relatively high density with few interactions but still keep in gaseous state.However,most categories of atom will easily tend to combine with others and form gaseous molecules or liquid.At first,people focused on the BEC of hydrogen atom,but failed to in the end. Fortunately,after the research,the alkali metal atoms with one electron in the outer shell and odd number of nuclei spin,which can be seen as bosons,were found suitable to obtain BEC in1980s.This is the first reason why it takes so long before BEC can be observed.Then,here’s a problem of cooling atom.Cooling atom make the kinetic energy of atom less.The breakthrough appeared in1960s when the laser was invented.In1975, the idea of laser cooling was advanced by Hänsch and Shallow.Here’s a chart of the development of laser cooling:Year Technique Limit Temperature Contributors 1980～Laser cooling of the atomic beam～mK Phillips,etc. 19853-D Laser cooling～240μK S.Chu,etc. 1989Sisyphus cooling～0.1~1μK Dalibard,etc. 1995Evaporative cooling～100nK S.Chu,etc. 1995The first realization of BEC～20nK JILA group Until1995,people didn’t have the cooling technique which was not perfect enough,so that’s the other answer.By the way,the Nobel Prize in Physics1997wasawarded to Stephen Chu,Claude Cohen－Tannoudji,and William D.Phillips for the contribution on laser cooling and trapping of atoms.3Anything you can add to the BEC phenomena(recent developments,etc.)from your own Reading.Bose–Einstein condensation of photons in an optical microcavity BEC is the state of bosons at extremely low temperature.According to the traditional view,photon does not have static mass,which means lower the temperature is,less the number of photons will be.It's very difficult for scientists to get Bose Einstein condensation of photons.Several German scientists said they obtained the BEC of photon successfully in the journal Nature published on November24th,2011.Their experiment confines photons in a curved-mirror optical microresonator filled with a dye solution,in which photons are repeatedly absorbed and re-emitted by the dye molecules.Those photons could‘heat’the dye molecules and be gradually cooled.The small distance of3.5 optical wavelengths between the mirrors causes a large frequency spacing between adjacent longitudinal modes.By pumping the dye with an external laser we add to a reservoir of electronic excitations that exchanges particles with the photon gas,in the sense of a grand-canonical ensemble.The pumping is maintained throughout the measurement to compensate for losses due to coupling into unconfined optical modes, finite quantum efficiency and mirror losses until they reach a steady state and become a super photons.(Klaers,J.,Schmitt,J.,Vewinger, F.,&Weitz,M.(2010).Bose-einstein condensation of photons in an optical microcavity.Nature,468(7323), 545-548.)With the BEC of photons,a brand new light source is created,which gives a possible to generate laser with extremely short wavelength,such as UV laser and X-ray laser.What’s more,it shows the future of powerful computer chip.Figure3.1Scheme of the experimental setup.4ConclusionA Bose-Einstein condensation(BEC)is a state of matter of a dilute gas of bosons cooled to temperatures very close to absolute zero.Under such conditions,a large fraction of bosons occupy the lowest quantum state,at which point macroscopic quantum phenomena become apparent.This state was first predicted,generally,in1924-25by Satyendra Nath Bose and Albert Einstein.And after70years,the Nobel Prize in Physics2001was awarded jointly to Eric A.Cornell,Wolfgang Ketterle and Carl E.Wieman"for theachievement of Bose-Einstein condensation in dilute gases of alkali atoms,and for early fundamental studies of the properties of the condensates".This achievement is not only related to the BEC theory but also the revolution of atom-cooling technique.5References[1]Pethick,C.,&Smith,H.(2001).Bose-einstein condensation in dilute gases.Bose-Einstein Condensation in Dilute Gases,56(6),414.[2]Klaers J,Schmitt J,Vewinger F,et al.Bose-Einstein condensation of photons in anoptical microcavity[J].Nature,2010,468(7323):545-548.[3]陈徐宗,&陈帅.(2002).物质的新状态——玻色-爱因斯坦凝聚——2001年诺贝尔物理奖介绍.物理,31(3),141-145.[4]Boson(n.d.)In Wikipedia.Retrieved from:</wiki/Boson>[5]Fermion(n.d.)In Wikipedia.Retrieved from:</wiki/Fermion>[6]Bose-einstein condensate(n.d.)In Wikipedia.Retrieved from:</wiki/Bose%E2%80%93Einstein_condensate>[7]玻色-爱因斯坦凝聚态(n.d.)In Baidubaike.Retrieved from:</link?url=5NzWN5riyBWC-qgPhvZ1QBcD2rdd4Tenkcw EyoEcOBhjh7-ofFra6uydj2ChtL-JvkPK78twjkfIC2gG2m_ZdK>。