945 research outputs found
Computing the ground state solution of Bose-Einstein condensates by a normalized gradient flow
In this paper, we prove the energy diminishing of a normalized gradient flow
which provides a mathematical justification of the imaginary time method used
in physical literatures to compute the ground state solution of Bose-Einstein
condensates (BEC). We also investigate the energy diminishing property for the
discretization of the normalized gradient flow. Two numerical methods are
proposed for such discretizations: one is the backward Euler centered finite
difference (BEFD), the other one is an explicit time-splitting sine-spectral
(TSSP) method. Energy diminishing for BEFD and TSSP for linear case, and
monotonicity for BEFD for both linear and nonlinear cases are proven.
Comparison between the two methods and existing methods, e.g. Crank-Nicolson
finite difference (CNFD) or forward Euler finite difference (FEFD), shows that
BEFD and TSSP are much better in terms of preserving energy diminishing
property of the normalized gradient flow. Numerical results in 1d, 2d and 3d
with magnetic trap confinement potential, as well as a potential of a stirrer
corresponding to a far-blue detuned Gaussian laser beam are reported to
demonstrate the effectiveness of BEFD and TSSP methods. Furthermore we observe
that the normalized gradient flow can also be applied directly to compute the
first excited state solution in BEC when the initial data is chosen as an odd
function.Comment: 28 pages, 6 figure
A finite element method with mesh adaptivity for computing vortex states in fast-rotating Bose-Einstein condensates
Numerical computations of stationary states of fast-rotating Bose-Einstein
condensates require high spatial resolution due to the presence of a large
number of quantized vortices. In this paper we propose a low-order finite
element method with mesh adaptivity by metric control, as an alternative
approach to the commonly used high order (finite difference or spectral)
approximation methods. The mesh adaptivity is used with two different numerical
algorithms to compute stationary vortex states: an imaginary time propagation
method and a Sobolev gradient descent method. We first address the basic issue
of the choice of the variable used to compute new metrics for the mesh
adaptivity and show that simultaneously refinement using the real and imaginary
part of the solution is successful. Mesh refinement using only the modulus of
the solution as adaptivity variable fails for complicated test cases. Then we
suggest an optimized algorithm for adapting the mesh during the evolution of
the solution towards the equilibrium state. Considerable computational time
saving is obtained compared to uniform mesh computations. The new method is
applied to compute difficult cases relevant for physical experiments (large
nonlinear interaction constant and high rotation rates).Comment: to appear in J. Computational Physic
Computation of Ground States of the Gross-Pitaevskii Functional via Riemannian Optimization
In this paper we combine concepts from Riemannian Optimization and the theory
of Sobolev gradients to derive a new conjugate gradient method for direct
minimization of the Gross-Pitaevskii energy functional with rotation. The
conservation of the number of particles constrains the minimizers to lie on a
manifold corresponding to the unit norm. The idea developed here is to
transform the original constrained optimization problem to an unconstrained
problem on this (spherical) Riemannian manifold, so that fast minimization
algorithms can be applied as alternatives to more standard constrained
formulations. First, we obtain Sobolev gradients using an equivalent definition
of an inner product which takes into account rotation. Then, the
Riemannian gradient (RG) steepest descent method is derived based on projected
gradients and retraction of an intermediate solution back to the constraint
manifold. Finally, we use the concept of the Riemannian vector transport to
propose a Riemannian conjugate gradient (RCG) method for this problem. It is
derived at the continuous level based on the "optimize-then-discretize"
paradigm instead of the usual "discretize-then-optimize" approach, as this
ensures robustness of the method when adaptive mesh refinement is performed in
computations. We evaluate various design choices inherent in the formulation of
the method and conclude with recommendations concerning selection of the best
options. Numerical tests demonstrate that the proposed RCG method outperforms
the simple gradient descent (RG) method in terms of rate of convergence. While
on simple problems a Newton-type method implemented in the {\tt Ipopt} library
exhibits a faster convergence than the (RCG) approach, the two methods perform
similarly on more complex problems requiring the use of mesh adaptation. At the
same time the (RCG) approach has far fewer tunable parameters.Comment: 28 pages, 13 figure
A Regularized Newton Method for Computing Ground States of Bose-Einstein condensates
In this paper, we propose a regularized Newton method for computing ground
states of Bose-Einstein condensates (BECs), which can be formulated as an
energy minimization problem with a spherical constraint. The energy functional
and constraint are discretized by either the finite difference, or sine or
Fourier pseudospectral discretization schemes and thus the original infinite
dimensional nonconvex minimization problem is approximated by a finite
dimensional constrained nonconvex minimization problem. Then an initial
solution is first constructed by using a feasible gradient type method, which
is an explicit scheme and maintains the spherical constraint automatically. To
accelerate the convergence of the gradient type method, we approximate the
energy functional by its second-order Taylor expansion with a regularized term
at each Newton iteration and adopt a cascadic multigrid technique for selecting
initial data. It leads to a standard trust-region subproblem and we solve it
again by the feasible gradient type method. The convergence of the regularized
Newton method is established by adjusting the regularization parameter as the
standard trust-region strategy. Extensive numerical experiments on challenging
examples, including a BEC in three dimensions with an optical lattice potential
and rotating BECs in two dimensions with rapid rotation and strongly repulsive
interaction, show that our method is efficient, accurate and robust.Comment: 25 pages, 6 figure
A new Sobolev gradient method for direct minimization of the Gross-Pitaevskii energy with rotation
In this paper we improve traditional steepest descent methods for the direct
minimization of the Gross-Pitaevskii (GP) energy with rotation at two levels.
We first define a new inner product to equip the Sobolev space and derive
the corresponding gradient. Secondly, for the treatment of the mass
conservation constraint, we use a projection method that avoids more
complicated approaches based on modified energy functionals or traditional
normalization methods. The descent method with these two new ingredients is
studied theoretically in a Hilbert space setting and we give a proof of the
global existence and convergence in the asymptotic limit to a minimizer of the
GP energy. The new method is implemented in both finite difference and finite
element two-dimensional settings and used to compute various complex
configurations with vortices of rotating Bose-Einstein condensates. The new
Sobolev gradient method shows better numerical performances compared to
classical or gradient methods, especially when high rotation rates
are considered.Comment: to appear in SIAM J Sci Computin
Accurate and efficient numerical methods for computing ground states and dynamics of dipolar Bose-Einstein condensates via the nonuniform FFT
In this paper, we propose efficient and accurate numerical methods for
computing the ground state and dynamics of the dipolar Bose-Einstein
condensates utilising a newly developed dipole-dipole interaction (DDI) solver
that is implemented with the non-uniform fast Fourier transform (NUFFT)
algorithm. We begin with the three-dimensional (3D) Gross-Pitaevskii equation
(GPE) with a DDI term and present the corresponding two-dimensional (2D) model
under a strongly anisotropic confining potential. Different from existing
methods, the NUFFT based DDI solver removes the singularity by adopting the
spherical/polar coordinates in Fourier space in 3D/2D, respectively, thus it
can achieve spectral accuracy in space and simultaneously maintain high
efficiency by making full use of FFT and NUFFT whenever it is necessary and/or
needed. Then, we incorporate this solver into existing successful methods for
computing the ground state and dynamics of GPE with a DDI for dipolar BEC.
Extensive numerical comparisons with existing methods are carried out for
computing the DDI, ground states and dynamics of the dipolar BEC. Numerical
results show that our new methods outperform existing methods in terms of both
accuracy and efficiency.Comment: 26 pages, 5 figure
Recommended from our members
Anisotropic collapse in three-dimensional dipolar Bose-Einstein condensates
- …