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

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)

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

Saved successfully!

Ooh no, something went wrong!