15.04.2013 Views

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

SHOW MORE
SHOW LESS

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)

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!