A normalized gradient flow method with attractive-repulsive splitting for computing ground states of Bose-Einstein condensates with higher-order interaction

TL;DR

This paper introduces a stable and efficient numerical method for computing ground states of Bose-Einstein condensates with higher-order interactions by splitting the interaction term and treating parts explicitly or semi-implicitly.

Contribution

The paper develops a novel splitting scheme for the modified Gross-Pitaevskii equation, improving stability and applicability to multidimensional problems and excited states.

Findings

Splitting improves scheme stability significantly.

Methods are effective for multidimensional problems.

Applicable to computing first excited states.

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.