Efficient spectral computation of the stationary states of rotating Bose-Einstein condensates by the preconditioned nonlinear conjugate gradient method
arXiv:1611.02045 · doi:10.1016/j.jcp.2017.04.040
Abstract
We propose a preconditioned nonlinear conjugate gradient method coupled with a spectral spatial dis-cretization scheme for computing the ground states (GS) of rotating Bose-Einstein condensates (BEC), modeled by the Gross-Pitaevskii Equation (GPE). We first start by reviewing the classical gradient flow (also known as imaginary time (IMT)) method which considers the problem from the PDE standpoint, leading to numerically solve a dissipative equation. Based on this IMT equation, we analyze the forward Euler (FE), Crank-Nicolson (CN) and the classical backward Euler (BE) schemes for linear problems and recognize classical power iterations, allowing us to derive convergence rates. By considering the alternative point of view of minimization problems, we propose the preconditioned gradient (PG) and conjugate gradient (PCG) methods for the GS computation of the GPE. We investigate the choice of the preconditioner, which plays a key role in the acceleration of the convergence process. The performance of the new algorithms is tested in 1D, 2D and 3D. We conclude that the PCG method outperforms all the previous methods, most particularly for 2D and 3D fast rotating BECs, while being simple to implement.
Cited by in corpus (24)
- Droplet crystal ground states of a dipolar Bose gas
- Pattern Formation in Quantum Ferrofluids: from Supersolids to Superglasses
- Computation of Ground States of the Gross-Pitaevskii Functional via Riemannian Optimization
- Supersolidity in Two-Dimensional Trapped Dipolar Droplet Arrays
- Mathematical models and numerical methods for spinor Bose-Einstein condensates
- Exciting the Goldstone Modes of a Supersolid Spin-Orbit-Coupled Bose Gas
- Variational theory for the ground state and collective excitations of an elongated dipolar condensate
- Numerical calculation of dipolar quantum droplet stationary states
- Second-order flows for computing the ground states of rotating Bose-Einstein condensates
- Ferromagnetism in an extended coherently-coupled atomic superfluid
- Vortex patterns in the almost-bosonic anyon gas
- Decoupled sound and amplitude modes in trapped dipolar supersolids
- How to realise a homogeneous dipolar Bose gas in the roton regime
- Interplay of Kelvin-Helmholtz and superradiant instabilities of an array of quantized vortices in a two-dimensional Bose--Einstein condensate
- GPU-accelerated solutions of the nonlinear Schrödinger equation for simulating 2D spinor BECs
- Computing the least action ground state of the nonlinear Schrödinger equation by a normalized gradient flow
- Superradiant phononic emission from the analog spin ergoregion in a two-component Bose-Einstein condensate
- Computing the action ground state for the rotating nonlinear Schrödinger equation
- Crystalline Phases of Laser-Driven Dipolar Bose-Einstein Condensates
- A normalized gradient flow method with attractive-repulsive splitting for computing ground states of Bose-Einstein condensates with higher-order interaction
- Convergence analysis of direct minimization and self-consistent iterations
- Computing ground states of Bose-Einstein Condensates with higher order interaction via a regularized density function formulation
- A Quasi-Optimal Factorization Preconditioner for Periodic Schrödinger Eigenstates in Anisotropically Expanding Domains
- Nehari manifold optimization and its application for finding unstable solutions of semilinear elliptic PDEs