Efficient and accurate methods for solving the time-dependent spin-1 Gross-Pitaevskii equation

TL;DR

This paper introduces new symplectic numerical methods for accurately solving the time-dependent spin-1 Gross-Pitaevskii equation, enabling better simulations of spinor Bose-Einstein condensates.

Contribution

It presents the first fully symplectic integrators for spin-1 condensate evolution, using a novel two-way splitting approach for exact flow solutions.

Findings

Second-order and fourth-order symplectic methods developed

Methods outperform existing approaches in numerical tests

Enhanced accuracy and stability in simulating spinor condensates

Abstract

We develop a numerical method for solving the spin-1 Gross-Pitaevskii equation. The basis of our work is a two-way splitting of the spin-1 evolution equation that leads to two exactly solvable flows. We use this to implement a second-order and a fourth-order symplectic integration method. These are the first fully symplectic methods for evolving spin-1 condensates. We develop two non-trivial numerical tests to compare our methods against two other approaches.