“…where α and β are constants. The nonlinear Schrödinger equation has been extensively studied by various numerical methods, such as finite element methods, [27] finite difference methods, [28] spectral methods, [29,30] splitting methods, [6,7,18,19,[31][32][33] etc. Among these numerical methods of different categories, the multisymplectic method has attracted special attention for its better numerical stability for long-time computations and perfect performance in preserving the intrinsic properties and conservation laws of nonlinear Schrödinger equations.…”