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