Symplectic integrator for nonlinear high order Schrödinger equation with a trapped term

Authors:
Linghua Kong;Jialin Hong;Lan Wang;Fangfang Fu
Affiliations:
School of Mathematics and Information Science, Jiangxi Normal University, Nanchang, Jiangxi, 330022, China;State Key Laboratory of Scientific and Engineering Computing, Institute of Computational Mathematics and Scientific/Engineering Computing, AMSS, CAS, P.O. Box 2719, Beijing, 100190, China;School of Mathematics and Information Science, Jiangxi Normal University, Nanchang, Jiangxi, 330022, China;School of Mathematics and Information Science, Jiangxi Normal University, Nanchang, Jiangxi, 330022, China
Venue:
Journal of Computational and Applied Mathematics
Year:
2009

Citing 6
Cited 1

Difference schemes for solving the generalized nonlinear Schrödinger equation

Journal of Computational Physics
Nonlinear solitary waves with Gaussian tails

Physica D
On the preservation of phase space structure under multisymplectic discretization

Journal of Computational Physics
Multi-symplectic Runge-Kutta methods for nonlinear Dirac equations

Journal of Computational Physics
Multi-symplectic Runge-Kutta-Nyström methods for nonlinear Schrödinger equations with variable coefficients

Journal of Computational Physics
Globally conservative properties and error estimation of a multi-symplectic scheme for Schrödinger equations with variable coefficients

Applied Numerical Mathematics

Local path fitting: A new approach to variational integrators

Journal of Computational and Applied Mathematics

Quantified Score

Hi-index	7.29

Visualization

Abstract

In the manuscript, we discuss the symplectic integrator for the numerical solution of a kind of high order Schrodinger equation with trapped terms. The Hamiltonian formulism is discovered for it. We first discretize the Hamiltonian system in space to reduce it to a finite-dimensional one. Then the symplectic midpoint scheme is applied to the temporal discretization. The symplectic scheme we devise is of second order accuracy in time and 2lth order accuracy in space. It is proved that it preserves the charge of the original equation veraciously. The energy is not preserved explicitly for the exception of the linear case. However, after some computation, the energy transit formula in the temporal direction is obtained. In the numerical part, we compare our scheme with some existing schemes, including the leap frog scheme and the energy-preserving scheme. From the numerical evidence, we find that our numerical schemes are efficient and available.The numerical results are consistent with the theoretical analysis.