2012/08/02 by Taras I. Lakoba, Lakoba, T. I.
Physics and Astronomy · #Advanced Fiber Laser Technologies #FOS: Mathematics #FOS: Physical sciences #Nonlinear Photonic Systems #Nonlinear Waves and Solitons #Numerical Analysis (math.NA) #Pattern Formation and Solitons (nlin.PS)
paper · pdf · doi:10.48550/arxiv.1208.0578
openalex publication_date 2012/08/02 · openalex created_date 2016/09/16 · openalex updated_date 2026/07/28
We consider the implementation of the split-step method where the linear part of the nonlinear Schrödinger 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.