Symplectic simulation of dark solitons motion for nonlinear Schrodinger equation

摘要

In this paper, we study symplectic simulation of dark solitons motion of nonlinear Schrodinger equation (NLSE). The Ablowitz-Ladik model (A-L model) of NLSE can be expressed as a non-canonical Hamiltonian system. By using splitting technique, we construct explicit splitting K-symplectic methods for the A-L model. On the other hand, the A-L model can be transformed into a canonical system and standard symplectic methods can be employed to perform numerical simulation. A second order K-symplectic method and a second order symplectic method are employed to simulate one dark soliton and two dark solitons motion for the A-L model and its canonicalized system respectively. By comparing with a third-order non-symplectic Runge-Kutta method, we show the superiorities of the two symplectic methods in long-term tracking the motion of dark solitons and preserving the invariants. We also compare the CPU times of K-symplectic methods and standard symplectic methods and show that the former ones are more efficient. The energy-preserving scheme is also applied for non-canonical Hamiltonian systems. The numerical results demonstrate that the K-symplectic methods can nearly preserve the energy, the discrete invariants of A-L model and conserved quantities of NLSE, but the energy-preserving scheme can only exactly preserve the energy.