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.

2.1 S<strong>in</strong>gle vortex l<strong>in</strong>e configuration 11<br />

explicit d-dependence of the heal<strong>in</strong>g length <strong>in</strong>side the logarithm, <strong>in</strong> order to account for the<br />

density dependence of the size of the vortex core.<br />

Eq.(2.6) shows that no excitation energy is carried by the system when d = R⊥ . Of course<br />

this estimate, be<strong>in</strong>g derived us<strong>in</strong>g the TF-approximation, is not accurate if we go too close<br />

to the border. If the surface of the condensate is described by a more realistic density profile<br />

the energy is nevertheless expected to vanish when the vortex l<strong>in</strong>e is sufficiently far outside<br />

the bulk region. With<strong>in</strong> the simplify<strong>in</strong>g assumptions made above we can conclude that the<br />

configuration with d = R⊥ corresponds to the absence of a vortex, while the transition from<br />

d = R⊥ to d =0describes the nucleation path of the vortex.<br />

Angular momentum of the vortex configuration<br />

The <strong>in</strong>clusion of vorticity is not only accompanied by an energy cost but also by the appearance<br />

of angular momentum. The simplest way to calculate the angular momentum associated with<br />

a displaced vortex is to assume axi-symmetric trapp<strong>in</strong>g (ε =0) <strong>and</strong> to work with the TFapproximation.<br />

In this limit the size of the vortex core is small compared to the radius of the<br />

condensate so that one can use the vortex-free expression n(r) =µ[1−(r⊥/R⊥) 2 −(z/Z) 2 ]/g<br />

for the density profile of the condensate, where r2 ⊥ = x2 +y2 <strong>and</strong> g =4πa¯h 2 /m is the coupl<strong>in</strong>g<br />

constant, fixed by the positive scatter<strong>in</strong>g length a. Then, one can write the angular momentum<br />

<strong>in</strong> the form<br />

<br />

<br />

Lz = m dz dr⊥r⊥n(r⊥,z) vvortex · dl , (2.7)<br />

where the l<strong>in</strong>e <strong>in</strong>tegral is taken along a circle of radius r⊥. Use of Stokes’ theorem gives the<br />

result [49]<br />

<br />

d<br />

Lz(d/R⊥) =N¯h 1 −<br />

<br />

2 5/2<br />

, (2.8)<br />

where d is the distance of the vortex l<strong>in</strong>e from the symmetry axis. Eq. (2.8) shows that the<br />

angular momentum per particle is reduced from the value ¯h as soon as the vortex is displaced<br />

from the center <strong>and</strong> becomes zero for d = R⊥.<br />

Vortex energy <strong>in</strong> the rotat<strong>in</strong>g frame<br />

Eq.(2.6) makes evident that a macroscopic energy is required to achieve a transition to a<br />

vortex-state. In a rotat<strong>in</strong>g trap of the form (1.3) the system may nevertheless like to acquire<br />

the vortex configuration. This happens if there is a total energy ga<strong>in</strong> <strong>in</strong> the rotat<strong>in</strong>g frame<br />

where the system has energy<br />

E(Ω) = E − ΩLz . (2.9)<br />

Here E is the energy <strong>in</strong> the laboratory frame, Lz is the angular momentum, <strong>and</strong> Ω is the<br />

angular velocity of the trap around the z-axis.<br />

In an axi-symmetric trap (ε =0) Eq.(2.9) yields the energy [49, 27]<br />

<br />

2 3/2 <br />

2 5/2<br />

d<br />

d<br />

Ev(d/R⊥, Ω,µ)=Ev(d =0,ε=0,µ) 1 −<br />

− ΩN¯h 1 −<br />

. (2.10)<br />

R⊥<br />

R⊥<br />

R⊥

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

Saved successfully!

Ooh no, something went wrong!