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

TL;DR

This paper investigates the unexpected instability of a finite-difference split-step method for the nonlinear Schrödinger equation when applied to soliton backgrounds, revealing modes supported by soliton sides.

Contribution

It provides a detailed analysis explaining the instability mechanism of the finite-difference split-step method on soliton backgrounds, contrasting it with spectral methods.

Findings

Finite-difference split-step method can become unstable on soliton backgrounds.

Unstable modes are supported by the sides of the soliton.

The instability mechanism differs from that of spectral split-step methods.

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.