be changed to
Ψ(r) = Y cn exρ[-¾<∕z, (2.7)
z-j n
. n
We can see now the unitary time evolution in eq. 2.6 has turned into an exponential decay.
The eigenenergy governs the decay rate, and so the eigenstate with the lowest energy,
i.e. the ground state of the system, decays slowest. Therefore, if we start from some
trial wavefunction (which could be a rough guess of the final solution), and apply suitable
renormalization of the wavefunction (e.g fixing the norm and/or the chemical potential)
during the imaginary time propagation, the trial wavefunction will tend towards the ground
state ψo.
2.2 Application to Dipolar Bose-Einstein Condensates
2.2.1 Dipolar Condensatesttheoretical formalism
In this section, we demonstrate a simple application of the split-step Fourier scheme, in
which we look for the ground state of a non-rotating scalar dipolar condensates. We start
from the 3D GP equation, which can be written down directly as following:
= ^^y~ +v(r) + col≠(r)∣2 + ⅛ f dr'~~Γ-ɜ-ɑɑ^- ⅛(r,)∣2 ψ(r), (2.8)
∂t 2m J ∣r - rψ
where c0 = ⅛πħ2a∣M, cd = μ0μ2∕(4π) with μ = 6μs for chromium atoms, one can easily
verify that cd∣co ≈ 3.6∕α(αβ). The trap is assumed to be harmonic and axially symmetric
V = ^rnω2.(x2 +y2 + Λ2z2). Finally, θ is the polar angle of r - r'. Following the standard