English

Error estimates of numerical methods for the nonlinear Dirac equation in the nonrelativistic limit regime

Numerical Analysis 2017-01-10 v3

Abstract

We present several numerical methods and establish their error estimates for the discretization of the nonlinear Dirac equation in the nonrelativistic limit regime, involving a small dimensionless parameter 0<ε10<\varepsilon\ll 1 which is inversely proportional to the speed of light. In this limit regime, the solution is highly oscillatory in time, i.e. there are propagating waves with wavelength O(ε2)O(\varepsilon^2) and O(1)O(1) in time and space, respectively. We begin with the conservative Crank-Nicolson finite difference (CNFD) method and establish rigorously its error estimate which depends explicitly on the mesh size hh and time step τ\tau as well as the small parameter 0<ε10<\varepsilon\le 1. Based on the error bound, in order to obtain `correct' numerical solutions in the nonrelativistic limit regime, i.e. 0<ε10<\varepsilon\ll 1, the CNFD method requests the ε\varepsilon-scalability: τ=O(ε3)\tau=O(\varepsilon^3) and h=O(ε)h=O(\sqrt{\varepsilon}). Then we propose and analyze two numerical methods for the discretization of the nonlinear Dirac equation by using the Fourier spectral discretization for spatial derivatives combined with the exponential wave integrator and time-splitting technique for temporal derivatives, respectively. Rigorous error bounds for the two numerical methods show that their ε\varepsilon-scalability is improved to τ=O(ε2)\tau=O(\varepsilon^2) and h=O(1)h=O(1) when 0<ε10<\varepsilon\ll 1 compared with the CNFD method. Extensive numerical results are reported to confirm our error estimates.

Keywords

Cite

@article{arxiv.1511.01192,
  title  = {Error estimates of numerical methods for the nonlinear Dirac equation in the nonrelativistic limit regime},
  author = {Weizhu Bao and Yongyong Cai and Xiaowei Jia and Jia Yin},
  journal= {arXiv preprint arXiv:1511.01192},
  year   = {2017}
}

Comments

35 pages. 1 figure. arXiv admin note: substantial text overlap with arXiv:1504.02881