A finite difference scheme for (2+1)D cubic-quintic nonlinear Schrödinger equations with nonlinear damping.

Solitons of the purely cubic nonlinear Schrödinger equation in a space dimension of n ≥ 2 suffer critical and supercritical collapses. These solitons can be stabilized in a cubic-quintic nonlinear medium. In this paper, we analyze the Crank-Nicolson finite difference scheme for the (2+1)D cubic-quintic nonlinear Schrödinger equation with cubic damping. We show that both the discrete solution, in the discrete L 2 -norm, and discrete energy are bounded. By using appropriate settings and estimations, the existence and the uniqueness of the numerical solution are proved. In addition, the error estimations are established in terms of second order for both space and time in discrete L 2 -norm and H 1 -norm. Numerical simulations for the (2+1)D cubic-quintic nonlinear Schrödinger equation with cubic damping are conducted to validate the convergence. • The (2+1)D cubic-quintic NLS equation can against 2D-soliton collapses. • A Crank-Nicolson finite difference scheme for the damped (2+1)D cubic-quintic NLS equation is proposed. • The existence and uniqueness of solutions are demonstrated using the boundedness of discrete energy and mass. • The scheme converges second order in time and space. • The numerical experiments are demonstrated to verify the theoretical results. [ABSTRACT FROM AUTHOR]