Alternating direction implicit method for solving two-dimensional cubic nonlinear Schrödinger equation

In this paper, four alternating direction implicit (ADI) schemes are presented for solving two-dimensional cubic nonlinear Schrödinger equations. Firstly, we give a Crank–Nicolson ADI scheme and a linearized ADI scheme both with accuracy O(Δt2+h2)O(Δt2+h2), with the same method, use fourth-order Padé compact difference approximation for the spatial discretization; two HOC-ADI schemes with accuracy O(Δt2+h4)O(Δt2+h4) are given. The two linearized ADI schemes apply extrapolation technique to the real coefficient of the nonlinear term to avoid iterating to solve. Unconditionally stable character is verified by linear Fourier analysis. The solution procedure consists of a number of tridiagonal matrix equations which make the computation cost effective. Numerical experiments are conducted to demonstrate the efficiency and accuracy, and linearized ADI schemes show less computational cost. All schemes given in this paper also can be used for two-dimensional linear Schrödinger equations.