Multi-symplectic Runge-Kutta methods for nonlinear Dirac equations

In this paper, we consider the multi-symplectic Runge-Kutta (MSRK) methods applied to the nonlinear Dirac equation in relativistic quantum physics, based on a discovery of the multi-symplecticity of the equation. In particular, the conservation of energy, momentum and charge under MSRK discretizations is investigated by means of numerical experiments and numerical comparisons with non-MSRK methods. Numerical experiments presented reveal that MSRK methods applied to the nonlinear Dirac equation preserve exactly conservation laws of charge and momentum, and conserve the energy conservation in the corresponding numerical accuracy to the method utilized. It is verified numerically that MSRK methods are stable and convergent with respect to the conservation laws of energy, momentum and charge, and MSRK methods preserve not only the inner geometric structure of the equation, but also some crucial conservative properties in quantum physics. A remarkable advantage of MSRK methods applied to the nonlinear Dirac equation is the precise preservation of charge conservation law.

[1]  S. Reich,et al.  Multi-symplectic integrators: numerical schemes for Hamiltonian PDEs that conserve symplecticity , 2001 .

[2]  Alberto Alvarez,et al.  Linearized Crank-Nicholson scheme for nonlinear Dirac equations , 1992 .

[3]  G. Akrivis.  A First Course In The Numerical Analysis Of Differential Equations [Book News & Reviews] , 1998, IEEE Computational Science and Engineering.

[4]  E MooreBrian,et al.  Multi-symplectic integration methods for Hamiltonian PDEs , 2003 .

[5]  J. Marsden,et al.  Introduction to mechanics and symmetry , 1994 .

[6]  Ying Liu,et al.  A novel numerical approach to simulating nonlinear Schro"dinger equations with varying coefficients , 2003, Appl. Math. Lett..

[7]  S. Reich.  Multi-Symplectic Runge—Kutta Collocation Methods for Hamiltonian Wave Equations , 2000 .

[8]  Geng Sun,et al.  The multi-symplecticity of partitioned Runge-Kutta methods for Hamiltonian PDEs , 2005, Math. Comput..

[9]  Brian E. Moore,et al.  Multi-symplectic integration methods for Hamiltonian PDEs , 2003, Future Gener. Comput. Syst..

[10]  Kuo Pen-Yu,et al.  The numerical study of a nonlinear one-dimensional Dirac equation , 1983 .

[11]  C. Schober,et al.  Geometric integrators for the nonlinear Schrödinger equation , 2001 .

[12]  Matthew West,et al.  Approximate momentum conservation for spatial semidiscretizations of semilinear wave equations , 2004, Numerische Mathematik.

[13]  J. Marsden,et al.  Multisymplectic Geometry, Variational Integrators, and Nonlinear PDEs , 1998, math/9807080.

[14]  B. Carreras,et al.  Interaction dynamics for the solitary waves of a nonlinear Dirac model , 1981 .

[15]  Brian E. Moore,et al.  Backward error analysis for multi-symplectic integration methods , 2003, Numerische Mathematik.

[16]  J. M. Sanz-Serna,et al.  Numerical Hamiltonian Problems , 1994 .

[17]  J. M. Sanz-Serna,et al.  Split-step spectral schemes for nonlinear Dirac systems , 1989 .

[18]  C. Schober,et al.  On the preservation of phase space structure under multisymplectic discretization , 2004 .

[19]  Ying Liu,et al.  Multisymplecticity of the centred box scheme for a class of hamiltonian PDEs and an application to quasi-periodically solitary waves , 2004 .

[20]  Jerrold E. Marsden,et al.  Variational Methods, Multisymplectic Geometry and Continuum Mechanics , 2001 .

[21]  Ying Liu,et al.  Globally conservative properties and error estimation of a multi-symplectic scheme for Schrödinger equations with variable coefficients , 2006 .