Lobatto IIIA–IIIB discretization of the strongly coupled nonlinear Schrödinger equation

In this paper, we construct a second order semi-explicit multi-symplectic integrator for the strongly coupled nonlinear Schrödinger equation based on the two-stage Lobatto IIIA–IIIB partitioned Runge–Kutta method. Numerical results for different solitary wave solutions including elastic and inelastic collisions, fusion of two solitons and with periodic solutions confirm the excellent long time behavior of the multi-symplectic integrator by preserving global energy, momentum and mass.