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
You also want an ePaper? Increase the reach of your titles
YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.
10.1 Current-phase dynamics <strong>in</strong> the tight b<strong>in</strong>d<strong>in</strong>g regime 135<br />
the values (10.2,10.3) <strong>and</strong> allow them to undergo a general time-dependent evolution<br />
Sl = kld − µ(k)<br />
¯h t → Sl(t) , (10.6)<br />
nl = n → nl(t) . (10.7)<br />
Furthermore, we require the shape of the wavefunction at site l to be well approximated by the<br />
Wannier function fl obta<strong>in</strong>ed for a Bloch state (10.1) with average density n = nl(t). This<br />
presupposes an adiabatic adaptation of the shape of the wavefunction to the <strong>in</strong>stantaneous<br />
value of nl(t). Given these assumptions the time-dependent state Ψ(z,t) is written <strong>in</strong> the<br />
form<br />
Ψ(z,t) = <br />
l<br />
<br />
fl(z; nl(t)) nl(t)de iSl(t)<br />
. (10.8)<br />
This ansatz must ensure that the quantities nl <strong>in</strong>deed have the mean<strong>in</strong>g of the average density<br />
at site l. Thus, (10.8) must satisfy the equation<br />
ld+d/2<br />
dz |Ψ(z,t)| 2 = nld. (10.9)<br />
ld−d/2<br />
This is achieved by requr<strong>in</strong>g fl(z; nl) to be very well localized at the site l such that fl±1<br />
have a negligible contribution to the population of site l ( ld+d/2<br />
ld−d/2 dzf 2 l±1 =0). In addition, we<br />
require the orthogonalization condition<br />
<br />
dzfl(z; nl(t))fl ′(z; nl ′(t)) = δl,l ′ , (10.10)<br />
to be approximately satisfied. Exact orthogonality is guaranteed of course only for a stationary<br />
state where nl = nl ′.<br />
Dynamical equations phases <strong>and</strong> site occupations<br />
The wavefunction Ψ(z,t) evolves accord<strong>in</strong>g to the time-dependent GP-equation (3.14). To<br />
obta<strong>in</strong> dynamical equations for the phases Sl <strong>and</strong> average densities nl, we <strong>in</strong>sert (10.8) <strong>in</strong>to<br />
the GPE. Upon multiplication of (3.14) by Ψ ∗ <strong>and</strong> <strong>in</strong>tegration over the total volume, us<strong>in</strong>g<br />
(10.10) we obta<strong>in</strong><br />
where<br />
˙nl =<br />
˙Sl = − µl<br />
¯h<br />
<br />
l ′ =l+1,l−1<br />
<br />
µl = fl<br />
+ <br />
δl,l′ √<br />
nlnl<br />
¯h<br />
′ s<strong>in</strong>(Sl − Sl ′) , (10.11)<br />
l ′ =l+1,l−1<br />
− ¯h2 ∂ 2 z<br />
2m<br />
δl,l′ <br />
µ nl<br />
2¯h<br />
′<br />
nl<br />
+ V (z)+gnld|fl| 2<br />
cos(Sl − Sl ′) , (10.12)<br />
<br />
fldz , (10.13)<br />
<strong>and</strong> δ l,l′<br />
µ are directly related to the overlap<br />
while the time-dependent tunnel<strong>in</strong>g parameters δl,l′ between two neighbour<strong>in</strong>g Wannier functions<br />
δ l,l′<br />
<br />
= −2 dz fl − ¯h2 ∂2 <br />
z<br />
+ V fl<br />
2m ′ + gnldfl|fl| 2 fl ′ + gnl ′dfl|fl ′|2fl ′<br />
<br />
,<br />
δ l,l′<br />
µ = δ l,l′<br />
<br />
− 4gnld fl|fl| 2 fl ′dz. (10.14)