Error estimates of finite difference methods for the Dirac equation in the massless and nonrelativistic regime

Abstract

We present four frequently used finite difference methods and establish the error bounds for the discretization of the Dirac equation in the massless and nonrelativistic regime, involving a small dimensionless parameter $0<\epsilon \ll 1$ inversely proportional to the speed of light. In the massless and nonrelativistic regime, the solution exhibits rapid motion in space and is highly oscillatory in time. Specifically, the wavelength of the propagating waves in time is at $O(\epsilon )$, while in space it is at $O(1)$ with the wave speed at $O({\epsilon }^{-1}).$ We adopt one leap-frog, two semi-implicit, and one conservative Crank-Nicolson finite difference methods to numerically discretize the Dirac equation in one dimension and establish rigorously the error estimates which depend explicitly on the time step $\tau$, mesh size $h$, as well as the small parameter $\epsilon$. The error bounds indicate that, to obtain the `correct' numerical solution in the massless and nonrelativistic regime, i.e. $0<\epsilon \ll 1$, all these finite difference methods share the same $\epsilon$-scalability as time step $\tau =O({\epsilon }^{3/2})$ and mesh size $h=O({\epsilon }^{1/2})$. A large number of numerical results are reported to verify the error estimates.

arXiv comment22 pages, 1 figure and 6 tables

References in corpus (4)

• A semiclassical approach to the Dirac equation
• A Split-Step Numerical Method for the Time-Dependent Dirac Equation in 3-D Axisymmetric Geometry
• A fourth-order compact time-splitting Fourier pseudospectral method for the Dirac equation
• Super-resolution of time-splitting methods for the Dirac equation in the nonrelativistic regime