Unconditional optimal-order error estimates of linear relaxation compact difference scheme for the coupled nonlinear Schrödinger system

Abstract: This paper presents a linear, decoupled, mass- and energy-conserving numerical scheme for the multi-dimensional coupled nonlinear Schrödinger (CNLS) system. The scheme combines the fourth-order compact difference approximation in space with the relaxation technique in a time-staggered mesh framework, solving the primal unknowns and introduced auxiliary relaxation variables sequentially with high efficiency and high-order accuracy. We show the unique solvability and discrete conservation laws of the developed scheme. In particular, for the first time, leveraging an auxiliary error equation system combined with the cut-off technique, optimal-order error estimates in the discrete H1-norm for the primal variables at the time nodes, and in the discrete L2-norm for the auxiliary relaxation variables at the intermediate time nodes, are rigorously proved without any coupling mesh conditions, which contribute to the primary theoretical contribution of this paper for multi-dimensional CNLS system. Numerical experiments demonstrate convincingly the strong performance of the proposed scheme in long-term simulations, maintaining both physical invariants and high-order accuracy.