The multisymplectic numerical method for Gross–Pitaevskii equation

For a Bose–Einstein Condensate placed in a rotating trap and confined in the z-axis, a multisymplectic difference scheme was constructed to investigate the evolution of vortices in this paper. First, we look for a steady state solution of the imaginary time G-P equation. Then, we numerically study the vortices's development in real time, starting with the solution in imaginary time as initial value.