Efficient numerical method for linear stability analysis of solitary waves Original Research Article

In this paper we use a numerical relaxation algorithm to improve and generalize the obtainment of the perturbation eigenstates of nonlinear systems. As a model problem we consider the linear stability analysis of the vortex eigenstates of the cubic–quintic nonlinear Schrödinger equation. It is shown by numerical calculations that the relaxation algorithm permits accurate tracing of complex perturbation eigenvalues.