High-order time splitting methods for the nonlinear Gross-Pitaevskii equation
This provides a more efficient numerical tool for simulating Bose-Einstein condensates, a domain-specific problem in computational physics.
The paper develops high-order time-splitting methods for the Gross-Pitaevskii equation that avoid negative time steps, reducing computational cost growth from exponential to quadratic with order. Numerical tests show accurate ground state and long-time conservation for orders up to 14.
We propose a high-order numerical methodology for computing the ground state and time evolution of the two-dimensional Gross-Pitaevskii equation with harmonic trapping potential. The ground state is obtained by combining normalized gradient flow with time-splitting schemes featuring strictly positive coefficients, which makes them suitable for both dissipative and Hamiltonian regimes and avoids the negative time steps required by classical symplectic methods. The computational cost of these methods grows quadratically with the order, in contrast to the exponential growth of symplectic alternatives. Numerical benchmarks assess ground state convergence under different initializations, long-time preservation of mass and Hamiltonian energy for varying mass constraints, and the cost-accuracy trade-off for orders q = 2, 4,..., 14.