Bose-Einstein Condensates in Rotating Traps and Optical ... - BEC
Bose-Einstein Condensates in Rotating Traps and Optical ... - BEC
Bose-Einstein Condensates in Rotating Traps and Optical ... - BEC
Create successful ePaper yourself
Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.
162 Condensate fraction<br />
This expression can be evaluated analytically. One obta<strong>in</strong>s the well known result<br />
∆Ntot<br />
=<br />
Ntot<br />
8 1<br />
3 π1/2 (a3n) 1/2 . (12.20)<br />
In view of the discussion of the case s = 0below, it is important to po<strong>in</strong>t out that the <strong>in</strong>tegr<strong>and</strong><br />
<strong>in</strong> (12.19) contributes most at p ∼ mc where the dispersion (12.17) passes from the phononic<br />
regime ¯hω ∼ p to the s<strong>in</strong>gle particle regime ¯hω ∼ p2 . Yet the convergence is very slow <strong>and</strong><br />
the <strong>in</strong>tegral is saturated by momenta much larger than mc [156], where elementary excitations<br />
have s<strong>in</strong>gle particle character.<br />
Quantum depletion of the 2D uniform gas<br />
Let us consider a uniform system of axial size d <strong>and</strong> transverse size L <strong>and</strong> let us assume that<br />
themotionisfrozen<strong>in</strong>thez-direction. We set Ψ(r) = 1<br />
√ d Ψ2(x, y). Under these conditions<br />
the stationary GPE takes the form<br />
<br />
− ¯h2<br />
The 2D density is given by<br />
The 2D coupl<strong>in</strong>g constant emerges as<br />
2m ∇2⊥ + g<br />
<br />
2<br />
|Ψ2|<br />
d<br />
Ψ2 = µΨ2 . (12.21)<br />
n2 = |Ψ2| 2 = Ntot/L 2 . (12.22)<br />
g2 = g ¯h2 a<br />
=4π . (12.23)<br />
d m d<br />
The groundstate solution is given by<br />
µ = g2n2 . (12.24)<br />
The calculation of the Bogoliubov spectrum is analogous to the 3D case. The result for<br />
the depletion changes only because of the change <strong>in</strong> the dimensionality of the <strong>in</strong>tegral over<br />
momenta. We rewrite the derivation here for clarity. The Bogoliubov equations for the system<br />
<strong>in</strong> the groundstate read<br />
<br />
p2 2m +2g2|Ψ2| 2 <br />
− µ up + g2|Ψ2| 2 vp , =¯hω(p)up<br />
(12.25)<br />
<br />
p2 2m +2g2|Ψ2| 2 <br />
− µ vp + g2|Ψ2| 2 up = −¯hω(p)vp , (12.26)<br />
where p 2 = p 2 x + p 2 y <strong>and</strong> we have assumed that the Bogoliubov amplitudes take the form<br />
1 1<br />
up(x, y) =Up √<br />
d L ei(pxx+pyy)/¯h , (12.27)<br />
1 1<br />
vp(x, y) =Vp √<br />
d L ei(pxx+pyy)/¯h . (12.28)