A new second-order accurate time discretization method for the nonlinear cubic Schrödinger equation

Abstract

A new second-order accurate semi-analytical time discretization method is introduced for the numerical solution of the one-dimensional nonlinear cubic Schrödinger equation. This method is based on the combination of the method of lines, Crank–Nicolson method, Newton method and Lanczos’ Tau method. It is a self-starting averaged two-time-level scheme that has proved to be stable, accurate and energy conservative for long time integration periods. At each time level, approximate solutions are sought on a segmented spatial interval as finite expansions in terms of a given orthogonal polynomial basis mapped appropriately onto each spatial subsegment. We have carried out numerical simulation concerning several cases for the propagation, collision and the bound states of solitons. Accurate results have been obtained using Chebyshev and Legendre polynomials. These results are well comparable with other published results obtained by the use of various standard numerical methods.