High-order time splitting methods for the nonlinear Gross-Pitaevskii equation

Fuente: arXiv
Saved in:
Bibliographic Details
Main Authors: Ben, Roberto, Besteiro, Agustín, Rial, Diego
Format: Preprint
Published: 2026
Subjects:
Online Access:
Tags: Add Tag
No Tags, Be the first to tag this record!
_version_ 1866913169719427072
author Ben, Roberto
Besteiro, Agustín
Rial, Diego
author_facet Ben, Roberto
Besteiro, Agustín
Rial, Diego
contents 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.
format Preprint
id arxiv_https___arxiv_org_abs_2605_29362
institution arXiv
publishDate 2026
record_format arxiv
spellingShingle High-order time splitting methods for the nonlinear Gross-Pitaevskii equation
Ben, Roberto
Besteiro, Agustín
Rial, Diego
Numerical Analysis
Analysis of PDEs
35Q55 - 65M70 - 35Q40
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.
title High-order time splitting methods for the nonlinear Gross-Pitaevskii equation
topic Numerical Analysis
Analysis of PDEs
35Q55 - 65M70 - 35Q40
url https://arxiv.org/abs/2605.29362