Analysis of a Crank-Nicolson finite difference scheme for (2+1)D perturbed nonlinear Schrödinger equations with saturable nonlinearity
We analyze a Crank-Nicolson finite difference discretization for the perturbed (2+1)D nonlinear Schrödinger equation with saturable nonlinearity and a perturbation of cubic loss. We show the boundedness, the existence and uniqueness of a numerical solution. We establish the error bound to prove the convergence of the numerical solution. Moreover, we find that the convergence rate is at the second order in both time step and spatial mesh size under a mild assumption. The numerical scheme is validated by the extensive simulations of the (2+1)D saturable nonlinear Schrödinger model with cubic loss. The simulations for travelling solitons are implemented by using an accelerated imaginary-time evolution scheme and the Crank-Nicolson finite difference method.
READ FULL TEXT