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.
122 Macroscopic Dynamics<br />
9.5 Small amplitude collective oscillations <strong>in</strong> the presence of harmonic<br />
trapp<strong>in</strong>g<br />
In current experiments, the external potential Vext is provided by a harmonic trap<br />
Vext = m<br />
2<br />
<br />
ω 2 xx 2 + ω 2 yy 2 + ω 2 zz 2<br />
(9.34)<br />
to which the optical lattice is superimposed. In the follow<strong>in</strong>g, we will assume the groundstate<br />
to be adequately described by the LDA developed <strong>in</strong> section 5.4, imply<strong>in</strong>g that the size of<br />
the condensate is much larger than the lattice spac<strong>in</strong>g d (see Eq.(5.20)). In accordance with<br />
the discussion above, this is a necessary condition for the use of (9.17,9.18). Moreover, we<br />
presuppose the chemical potential <strong>in</strong> absence of the harmonic conf<strong>in</strong>ement to exhibit the l<strong>in</strong>ear<br />
dependence on average density (see Eq.(5.13))<br />
µopt =˜gn + µgn=0 . (9.35)<br />
The validity of this approximation has been discussed <strong>in</strong> section 5.2. With<strong>in</strong> LDA, the groundstate<br />
macroscopic density profile ¯nM is then def<strong>in</strong>ed by the relation (see discussion <strong>in</strong> section<br />
5.4)<br />
µ =˜g¯nM(r)+µgn=0 + Vext . (9.36)<br />
The dynamic equations (9.17,9.18) do not only require the equation of state µopt(nM)<br />
as an <strong>in</strong>put, but also the effective mass m ∗ (nM). As a simplification, we will consider the<br />
effective mass to be density-<strong>in</strong>dependent <strong>and</strong> thus to be given by the s<strong>in</strong>gle particle effective<br />
mass<br />
m ∗ = m ∗ (gn/ER =0). (9.37)<br />
The applicability of this approximation has been discussed <strong>in</strong> section 6.1.<br />
With Vext given by the harmonic potential (9.34) <strong>and</strong> the approximations (9.35,9.37), the<br />
hydrodynamic equations (9.17,9.18) take the form<br />
∂<br />
∂t nM + ∂x(vMxnM)+∂y(vMynM)+ m<br />
m∗ ∂z (vMznM) =0, (9.38)<br />
m ∂<br />
∂t vM<br />
<br />
+ ∇ Vext +˜gnM + µgn=0 + m<br />
2 v2 Mx + m<br />
2 v2 My + m m<br />
nM<br />
m∗ 2 v2 <br />
Mz =0. (9.39)<br />
These equations are expected to give a correct description of the dynamics occurr<strong>in</strong>g on a<br />
length scale of the order of the system size. In particular, we will explore here the limit of<br />
small amplitude collective oscillations which comply with this condition.<br />
Hydrodynamic equations for small amplitude oscillations<br />
Small amplitude collective oscillations are associated with an oscillat<strong>in</strong>g perturbation <strong>in</strong> the<br />
groundstate density ¯nM <strong>and</strong> a small oscillat<strong>in</strong>g velocity field<br />
∆n(r,t) ∝ ∆n(r) e −iωt , (9.40)<br />
∆v(r,t) ∝ ∆v(r) e −iωt . (9.41)