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ö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.

arXiv comment28 pages, 1 figures, 2 tables

References in corpus (8)

• Solitons in nonlinear lattices
• Numerical Solution of the Gross-Pitaevskii Equation for Bose-Einstein Condensation
• Computational methods for the dynamics of the nonlinear Schrodinger/Gross-Pitaevskii equations
• Three Dimensional Simulation of Jet Formation in Collapsing Condensates
• Criteria for the experimental observation of multi-dimensional optical solitons in saturable media
• Stabilization of two-dimensional solitons in cubic-saturable nonlinear lattices
• Radiation dynamics in fast soliton collisions in the presence of cubic loss
• Collision-induced amplitude dynamics of fast 2D solitons in saturable nonlinear media with weak nonlinear loss