Instability of the finite-difference split-step method on the background of a soliton of the nonlinear Schrodinger equation

Abstract: We consider the implementation of the split-step method where the linear part of the nonlinear Schr\"odinger equation is solved using a finite-difference discretization of the spatial derivative. The von Neumann analysis predicts that this method is unconditionally stable on the background of a constant-amplitude plane wave. However, simulations show that the method can become unstable on the background of a soliton. We present an analysis explaining this instability. Both this analysis and the instability itself are substantially different from those of the Fourier split-step method, which computes the spatial derivative by spectral discretization. We also found that the modes responsible for the numerical instability are supported by the sides of the soliton, in contrast to unstable modes of linearized nonlinear wave equations, which (the modes) are supported by the soliton's core.