Second-order flows for computing the ground states of rotating Bose-Einstein condensates

05/02/2022
by   Haifan Chen, et al.
0

Second-order flows in this paper refer to some artificial evolutionary differential equations involving second-order time derivatives distinguished from gradient flows which are considered to be first-order flows. This is a popular topic due to the recent advances of inertial dynamics with damping in convex optimization. Mathematically, the ground state of a rotating Bose-Einstein condensate (BEC) can be modelled as a minimizer of the Gross-Pitaevskii energy functional with angular momentum rotational term under the normalization constraint. We introduce two types of second-order flows as energy minimization strategies for this constrained non-convex optimization problem, in order to approach the ground state. The proposed artificial dynamics are novel second-order nonlinear hyperbolic partial differential equations with dissipation. Several numerical discretization schemes are discussed, including explicit and semi-implicit methods for temporal discretization, combined with a Fourier pseudospectral method for spatial discretization. These provide us a series of efficient and robust algorithms for computing the ground states of rotating BECs. Particularly, the newly developed algorithms turn out to be superior to the state-of-the-art numerical methods based on the gradient flow. This is judged by the fact that larger size of time steps can be applied in explicit methods, and much smaller number of time steps are needed in both explicit and semi-implicit discretizations, while every time step encounters almost the same computational complexity in comparison with the respective gradient flow type approaches. Rich and detailed numerical examples are documented for verification and comparison.

READ FULL TEXT

page 27

page 29

page 30

research
07/04/2019

Second-order semi-implicit projection methods for micromagnetics simulations

Micromagnetics simulations require accurate approximation of the magneti...
research
08/06/2021

A second-order semi-implicit method for the inertial Landau-Lifshitz-Gilbert equation

Recent theoretical and experimental advances show that the inertia of ma...
research
08/26/2019

Stochastic dynamical modeling of turbulent flows

Advanced measurement techniques and high performance computing have made...
research
03/09/2022

Generalized SAV-exponential integrator schemes for Allen-Cahn type gradient flows

The energy dissipation law and the maximum bound principle (MBP) are two...
research
09/01/2019

An efficient method for computing stationary states of phase field crystal models

Computing stationary states is an important topic for phase field crysta...
research
01/11/2022

Path differentiability of ODE flows

We consider flows of ordinary differential equations (ODEs) driven by pa...
research
08/24/2019

Computing ground states of Bose-Einstein Condensates with higher order interaction via a regularized density function formulation

We propose and analyze a new numerical method for computing the ground s...

Please sign up or login with your details

Forgot password? Click here to reset