Higher-order splitting algorithms for solving the nonlinear Schrödinger equation and their instabilities

Abstract

Since the kinetic and the potential energy term of the real time nonlinear Schrödinger equation can each be solved exactly, the entire equation can be solved to any order via splitting algorithms. We verified the fourth-order convergence of some well known algorithms by solving the Gross-Pitaevskii equation numerically. All such splitting algorithms suffer from a latent numerical instability even when the total energy is very well conserved. A detail error analysis reveals that the noise, or elementary excitations of the nonlinear Schrödinger, obeys the Bogoliubov spectrum and the instability is due to the exponential growth of high wave number noises caused by the splitting process. For a continuum wave function, this instability is unavoidable no matter how small the time step. For a discrete wave function, the instability can be avoided only for $\dt k_{max}^2{<\atop\sim}2 π$, where $k_{max}=\pi /\Delta x$.

10 pages, 8 figures, submitted to Phys. Rev. E

References in corpus (3)

• Split-step Fourier methods for the Gross-Pitaevskii equation
• Exact evolution of time-reversible symplectic integrators and their phase error for the harmonic oscillator
• Any-order propagation of the nonlinear Schroedinger equation

Cited by in corpus (10)

• High order Chin actions in path integral Monte Carlo
• Explicit symplectic algorithms based on generating functions for charged particle dynamics
• Numerical method for evolving the Projected Gross-Pitaevskii equation
• Efficient and accurate methods for solving the time-dependent spin-1 Gross-Pitaevskii equation
• Analog cosmological reheating in an ultracold Bose gas
• Numerical method for evolving the dipolar projected Gross-Pitaevskii equation
• Optimal nonlinear coherent mode transitions in Bose-Einstein Condensates utilizing spatio-temporal controls
• Towards a QMC-based density functional including finite-range effects: excitation modes of a ${}^{39}$K quantum droplet
• Solving the spin-2 Gross-Pitaevskii equation using exact nonlinear dynamics and symplectic composition
• Ab-initio Quantum Monte Carlo study of ultracold atomic mixtures