Analysis of a Crank-Nicolson finite difference scheme for (2+1)D perturbed nonlinear Schrödinger equations with saturable nonlinearity

Abstract: We analyze a Crank-Nicolson finite difference discretization for the perturbed (2+1)D nonlinear Schr\"odinger 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\"odinger 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.