English

Spatial resolution of different discretizations over long-time for the Dirac equation with small potentials

Numerical Analysis 2021-05-24 v1 Numerical Analysis

Abstract

We compare the long-time error bounds and spatial resolution of finite difference methods with different spatial discretizations for the Dirac equation with small electromagnetic potentials characterized by ε(0,1]\varepsilon \in (0, 1] a dimensionless parameter. We begin with the simple and widely used finite difference time domain (FDTD) methods, and establish rigorous error bounds of them, which are valid up to the time at O(1/ε)O(1/\varepsilon). In the error estimates, we pay particular attention to how the errors depend explicitly on the mesh size hh and time step τ\tau as well as the small parameter ε\varepsilon. Based on the results, in order to obtain "correct" numerical solutions up to the time at O(1/ε)O(1/\varepsilon), the ε\varepsilon-scalability (or meshing strategy requirement) of the FDTD methods should be taken as h=O(ε1/2)h = O(\varepsilon^{1/2}) and τ=O(ε1/2)\tau = O(\varepsilon^{1/2}). To improve the spatial resolution capacity, we apply the Fourier spectral method to discretize the Dirac equation in space. Error bounds of the resulting finite difference Fourier pseudospectral (FDFP) methods show that they exhibit uniform spatial errors in the long-time regime, which are optimal in space as suggested by the Shannon's sampling theorem. Extensive numerical results are reported to confirm the error bounds and demonstrate that they are sharp.

Keywords

Cite

@article{arxiv.2105.10468,
  title  = {Spatial resolution of different discretizations over long-time for the Dirac equation with small potentials},
  author = {Yue Feng and Jia Yin},
  journal= {arXiv preprint arXiv:2105.10468},
  year   = {2021}
}

Comments

23 pages, 3 figures. arXiv admin note: text overlap with arXiv:1504.02881 by other authors

R2 v1 2026-06-24T02:21:05.190Z