A normalized gradient flow method with attractive-repulsive splitting for computing ground states of Bose-Einstein condensates with higher-order interaction
arXiv:1711.02506 · doi:10.1016/j.jcp.2018.04.038
Abstract
In this paper, we generalize the normalized gradient flow method to compute the ground states of Bose-Einstein condensates (BEC) with higher order interactions (HOI), which is modelled via the modified Gross-Pitaevskii equation (MGPE). Schemes constructed in naive ways suffer from severe stability problems due to the high restrictions on time steps. To build an efficient and stable scheme, we split the HOI term into two parts with each part treated separately. The part corresponding to a repulsive/positive energy is treated semi-implicitly while the one corresponding to an attractive/negative energy is treated fully explicitly. Based on the splitting, we construct the BEFD-splitting and BESP-splitting schemes. A variety of numerical experiments shows that the splitting will improve the stability of the schemes significantly. Besides, we will show that the methods can be applied to multidimensional problems and to the computation of the first excited state as well.
References in corpus (7)
- Efficient spectral computation of the stationary states of rotating Bose-Einstein condensates by the preconditioned nonlinear conjugate gradient method
- Energy-dependent effective interactions for dilute many-body systems
- Computation of Ground States of the Gross-Pitaevskii Functional via Riemannian Optimization
- Stability of a BEC with Higher-order Interactions near a Feshbach Resonance
- A simple model for interactions and corrections to the Gross-Pitaevskii Equation
- Thomas-Fermi Approximation for a Condensate with Higher-order Interactions
- Ground states of Bose-Einstein condensates with higher order interaction