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

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<\epsilon \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({\epsilon }^2)$ and $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 $h$ and time step $\tau$ as well as the small parameter $0<\epsilon \le 1$. Based on the error bound, in order to obtain `correct' numerical solutions in the nonrelativistic limit regime, i.e. $0<\epsilon \ll 1$, the CNFD method requests the $\epsilon$-scalability: $\tau =O({\epsilon }^3)$ and $h=O(\sqrt{\epsilon })$. 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 $\epsilon$-scalability is improved to $\tau =O({\epsilon }^2)$ and $h=O(1)$ when $0<\epsilon \ll 1$ compared with the CNFD method. Extensive numerical results are reported to confirm our error estimates.

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

References in corpus (13)

• The electronic properties of graphene
• Mathematical theory and numerical methods for Bose-Einstein condensation
• Hidden pseudospin and spin symmetries and their origins in atomic nuclei
• Giant Nonlocality near the Dirac Point in Graphene
• On the eigenvalues of operators with gaps. Application to Dirac operators
• Wave packets in Honeycomb Structures and Two-Dimensional Dirac Equations
• Numerical methods for nonlinear Dirac equation
• Resonantly enhanced pair production in a simple diatomic model
• A dispersion and norm preserving finite difference scheme with transparent boundary conditions for the Dirac equation in (1+1)D
• Solitary wave in the Nonlinear Dirac Equation with arbitrary nonlinearity
• Stability of solitary waves in the nonlinear Dirac equation with arbitrary nonlinearity
• A uniformly accurate (UA) multiscale time integrator pseudospectral method for the Dirac equation in the nonrelativistic limit regime
• A Split-Step Numerical Method for the Time-Dependent Dirac Equation in 3-D Axisymmetric Geometry

Cited by in corpus (5)

• Comparison of numerical methods for the nonlinear Klein-Gordon equation in the nonrelativistic limit regime
• A fourth-order compact time-splitting Fourier pseudospectral method for the Dirac equation
• Uniform error bounds of time-splitting methods for the nonlinear Dirac equation in the nonrelativistic regime without magnetic potential
• A fourth-order compact time-splitting method for the Dirac equation with time-dependent potentials
• Study of instability of the Fourier split-step method for the massive Gross--Neveu model