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.

80 Stationary states of a <strong>BEC</strong> <strong>in</strong> an optical lattice<br />

that the effective coupl<strong>in</strong>g constant (5.15) <strong>in</strong> the tight b<strong>in</strong>d<strong>in</strong>g regime reads<br />

L/2<br />

˜g = gd<br />

−L/2<br />

f 4 gn=0(z) . (6.45)<br />

A quantitative estimate of the validity of this expression will be given <strong>in</strong> the follow<strong>in</strong>g section<br />

us<strong>in</strong>g a gaussian ansatz for the Wannier function.<br />

Gaussian approximation to the Wannier function of the lowest b<strong>and</strong><br />

In the limit of zero tunnel<strong>in</strong>g (s →∞), Wannier functions become eigenstates. In particular,<br />

a s<strong>in</strong>gle particle <strong>in</strong> the groundstate can be described by the correspond<strong>in</strong>g harmonic oscillator<br />

wavefunction. If s ≫ 1 <strong>and</strong> as long as the s<strong>in</strong>gle-well condensate is far from be<strong>in</strong>g correctly<br />

described by the TF-approximation <strong>in</strong>side each well, it is reasonable to approximate the Wannier<br />

function by a gaussian<br />

ϕ(z) ≡ f(z) =<br />

1<br />

π 1/4√ σ exp(−z2 /2σ 2 ) . (6.46)<br />

The gaussian (6.46) is a useful variational ansatz for the calculation of on-site quantities<br />

as the compressibility, but fails to describe properties which crucially depend on the overlap<br />

between neighbour<strong>in</strong>g Wannier functions, as the effective mass. In fact, when tunnel<strong>in</strong>g is<br />

possible, but its effects only small, we expect the Wannier function of the lowest b<strong>and</strong> to be<br />

still similar to the gaussian (6.46) <strong>in</strong>side a well, but to have oscillat<strong>in</strong>g tails <strong>in</strong> the barrier region<br />

<strong>in</strong> order to ensure orthogonality.<br />

At a given lattice depth s <strong>and</strong> <strong>in</strong>teraction gn/ER, the value of the gaussian width σ is fixed<br />

by requir<strong>in</strong>g the ansatz (6.46) to m<strong>in</strong>imize the energy of the system. S<strong>in</strong>ce contributions due<br />

to tunnel<strong>in</strong>g go beyond the accuracy of this description, it is consistent to consider only the<br />

on-site energy ε0 (see Eq.(6.30)). Moreover, we exp<strong>and</strong> the lattice potential (3.6) around its<br />

m<strong>in</strong>ima<br />

V (z)/ER ≈ s<br />

2 πz<br />

−<br />

d<br />

s<br />

3<br />

4 πz<br />

, (6.47)<br />

d<br />

where the first term corresponds to a harmonic potential of oscillator frequency<br />

˜ω =2 √ sER/¯h, (6.48)<br />

while the second term allows for anharmonicity effects of O(z4 ).Thesolutionσhas to satisfy<br />

the equation<br />

− d3<br />

π3 1<br />

+ sπ σ − sπ3<br />

σ3 d d3 σ3 − 1<br />

<br />

gn π d<br />

2 ER 2<br />

2<br />

π2 1<br />

=0. (6.49)<br />

σ2 S<strong>in</strong>ce <strong>in</strong> we are <strong>in</strong>terested <strong>in</strong> the large-s limit, we neglect the <strong>in</strong>teraction term. Us<strong>in</strong>g s ≫<br />

1, σπ/d

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

Saved successfully!

Ooh no, something went wrong!