An efficient and spectrally accurate numerical method for computing dynamics of rotating Bose-Einstein condensates

TL;DR

This paper introduces a new spectral method for accurately and efficiently simulating the dynamics of rotating Bose-Einstein condensates using a time-splitting approach that is explicit, stable, and conserves key physical properties.

Contribution

The paper develops a novel time-splitting spectral method for the Gross-Pitaevskii equation with rotation, offering spectral accuracy, stability, and conservation properties for BEC simulations.

Findings

Method is explicit and unconditionally stable.

Achieves spectral accuracy in space and second order in time.

Numerical results confirm stability and conservation properties.

Abstract

In this paper, we propose an efficient and spectrally accurate numerical method for computing the dynamics of rotating Bose-Einstein condensates (BEC) in two dimensions (2D) and 3D based on the Gross-Pitaevskii equation (GPE) with an angular momentum rotation term. By applying a time-splitting technique for decoupling the nonlinearity and properly using the alternating direction implicit (ADI) technique for the coupling in the angular momentum rotation term in the GPE, at every time step, the GPE in rotational frame is decoupled into a nonlinear ordinary differential equation (ODE) and two partial differential equations with constant coefficients. This allows us to develop new time-splitting spectral (TSSP) methods for computing the dynamics of BEC in a rotational frame.The new numerical method is explicit, unconditionally stable, and of spectral accuracy in space and second order accuracy in time. Moreover, it is time reversible and time transverse invariant, and conserves the position density in the discretized level if the GPE does. Extensive numerical results are presented to confirm the above properties of the new numerical method for rotating BEC in 2D and 3D.