Bose 2.2 User Manual Page 109

  • Download
  • Add to my manuals
  • Print
  • Page
    / 155
  • Table of contents
  • BOOKMARKS
  • Rated. / 5. Based on customer reviews
Page view 108
6.2.3 Numerical Simulations of the Collapse Dynamics
In a collaboration with the theory group in Hannover, we have examined more closely
the collapse dynamics of the dipolar condensate in the lattice by means of numerical
simulations. The simulations are based on the time-dependent Gross-Pitaevskii equation
(given by Eq. (6.1)),
i~
t
ψ(r, t) =
~
2
2m
2
+ V
ext
(r) i~
L
3
2
N
2
|ψ(r, t)|
4
+ N
4π~
2
a
m
|ψ(r, t)|
2
+ N
Z
d
3
r
0
V
0
dd
(r r
0
) |ψ(r
0
, t)|
2
ψ(r, t),
(6.1)
where the non-unitary term proportional to
L
3
describes three-body atom losses due
to inelastic collisions.
N
is the initial number of atoms in the condensate and
V
ext
(
r
)
is the combined trapping potential of ODT and lattice with
V
ext
(
r
) =
U sin
2
(
πz/d
lat
) +
m
P
i=x,y,z
ω
2
i
r
2
i
/
2. The contact interactions are included via the
s
-wave scattering length
a and the dipole-dipole interaction potential V
0
dd
(r r
0
) is given by Eq. (2.6).
As a first step, the ground-state of the condensate is calculated for a scattering length
a
= 2
a
0
and a lattice depth
U
=
U
init
= 12
.
6
E
R
. This task is performed by integrating
the GPE in imaginary time, with the loss parameter
L
3
set to zero. Then, the simulations
are performed by a real-time evolution of Eq.
(6.1)
, following our experimental sequence
and choosing
100
L
3
= 2
·
10
40
m
6
/s
. In analogy to the experimental imaging procedure,
the time-evolution of the system, shown in Fig. 6.6(a), is illustrated by integrating the
atomic density along the
x
-direction. After sufficiently long expansion of the dBEC in the
TOF, the dilute cloud undergoes a ballistic expansion, since interactions do not play a
role anymore. Calculating the momentum distribution of the system in this regime (see
Fig. 6.6(b)) yields the spatial density distribution of the atomic cloud in the far field of
the TOF, corresponding to the absorption images taken in the experiment. Finally, the
time evolution of the number
N
of coherent atoms is extracted from the simulations, with
the results presented in Fig. 6.6(c). We note that the time origin in the simulations,
t
= 0,
is set at the end of the lattice ramp. Furthermore, a fixed holding time
t
hold
= 0
.
6
ms
is
chosen, before starting the TOF.
Collapse dynamics
We split the discussion of the collapse dynamics, shown in Fig. 6.6, into three parts: First,
we discuss the deconfinement-induced collapse of a single dBEC (
U
= 0
E
R
). Then, we
consider the stable in-trap configurations (
U
= 8
.
2
E
R
and
U
= 12
.
6
E
R
) and finally we
address the unstable situations at the lattice depths U = 3.2 E
R
and U = 6.3 E
R
.
(i)
By choosing the final lattice depth
U
= 0
E
R
, we obtain an unstable in-trap
configuration of the system. At the end of the lattice ramp, the periodic density
100
The loss rate
L
3
is chosen equal to the value that is used in the simulations of the
d
-wave collapse of a
single dBEC [36] (see section 6.1).
109
Page view 108
1 2 ... 104 105 106 107 108 109 110 111 112 113 114 ... 154 155

Comments to this Manuals

No comments