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.
48 Groundstate of a <strong>BEC</strong> <strong>in</strong> an optical lattice<br />
the effects of harmonic trapp<strong>in</strong>g have been presented <strong>in</strong> [102]. In this thesis, also results for the<br />
density profile are <strong>in</strong>cluded <strong>and</strong> the dependence on lattice depth of the density at the center<br />
of the harmonic trap averaged over one lattice period is described. We briefly comment on<br />
the generalization to 2D lattices. Based on an ansatz for the order parameter valid <strong>in</strong> deep<br />
lattices, the effective coupl<strong>in</strong>g constant ˜g <strong>and</strong> its effect on the harmonically trapped system<br />
have been previously discussed <strong>in</strong> [70].<br />
5.1 Density profile, energy <strong>and</strong> chemical potential<br />
When look<strong>in</strong>g for stationary solutions, the GP-equation (3.14) takes the form<br />
<br />
− ¯h2 ∂<br />
2m<br />
2<br />
∂z2 + sER s<strong>in</strong> 2<br />
<br />
πz<br />
+ g |Ψ(z)|<br />
d<br />
2<br />
<br />
Ψ(z) =µΨ(z) . (5.1)<br />
The groundstate is given by the solution of this equation with the lowest energy.<br />
Let us rewrite Eq.(5.1) <strong>in</strong> a more convenient form. First, we <strong>in</strong>troduce the rescaled order<br />
parameter<br />
<br />
ϕ(z) =<br />
L2 Ψ(z) ,<br />
N<br />
(5.2)<br />
with L the transverse size of the system <strong>and</strong> N the number of particles per lattice site such<br />
that d/2<br />
dx |ϕ(z)|<br />
−d/2<br />
2 The GPE (5.1) then reads<br />
=1. (5.3)<br />
<br />
− ¯h2 ∂<br />
2m<br />
2<br />
∂z2 + sER s<strong>in</strong> 2<br />
<br />
πz<br />
+ dgn|ϕ(z)|<br />
d<br />
2<br />
<br />
ϕ(z) =µϕ(z) , (5.4)<br />
where we have <strong>in</strong>troduced the average density<br />
n = N<br />
. (5.5)<br />
dL2 Cast<strong>in</strong>g Eq. (5.4) <strong>in</strong> dimensionless form, <strong>in</strong> dimensionless form it is possible to identify the<br />
govern<strong>in</strong>g parameters of the problem. We shall measure length <strong>in</strong> units of d/π, momentum<br />
<strong>in</strong> units of the Bragg-momentum qB =¯hπ/d <strong>and</strong> energy <strong>in</strong> units of the recoil energy ER =<br />
¯h 2 π 2 /2md 2 . In this way, we obta<strong>in</strong> the dimensionless GPE<br />
<br />
− ∂2<br />
∂z 2 + s s<strong>in</strong>2 (z)+π gn<br />
|ϕ(z)|<br />
ER<br />
2<br />
<br />
ϕ(z) = µ<br />
ER<br />
ϕ(z) . (5.6)<br />
In contrast to the case of a s<strong>in</strong>gle particle discussed <strong>in</strong> chapter 4, there are now two govern<strong>in</strong>g<br />
parameters: the lattice depth s <strong>and</strong> the <strong>in</strong>teraction parameter gn/ER. The latter quantity can<br />
be changed by vary<strong>in</strong>g the average density or the scatter<strong>in</strong>g length by means of a Feshbach<br />
resonance.