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.
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