English

Uniform error bounds of exponential wave integrator methods for the long-time dynamics of the Dirac equation with small potentials

Numerical Analysis 2021-06-29 v1 Numerical Analysis

Abstract

Two exponential wave integrator Fourier pseudospectral (EWI-FP) methods are presented and analyzed for the long-time dynamics of the Dirac equation with small potentials characterized by ε(0,1]\varepsilon \in (0, 1] a dimensionless parameter. Based on the (symmetric) exponential wave integrator for temporal derivatives in phase space followed by applying the Fourier pseudospectral discretization for spatial derivatives, the EWI-FP methods are explicit and of spectral accuracy in space and second-order accuracy in time for any fixed ε=ε0\varepsilon = \varepsilon_0. Uniform error bounds are rigorously carried out at O(hm0+τ2)O(h^{m_0}+\tau^2) up to the time at O(1/ε)O(1/\varepsilon) with the mesh size hh, time step τ\tau and m0m_0 an integer depending on the regularity of the solution. Extensive numerical results are reported to confirm our error bounds and comparisons of two methods are shown. Finally, dynamics of the Dirac equation in 2D are presented to validate the numerical schemes.

Keywords

Cite

@article{arxiv.2106.14107,
  title  = {Uniform error bounds of exponential wave integrator methods for the long-time dynamics of the Dirac equation with small potentials},
  author = {Yue Feng and Jia Yin},
  journal= {arXiv preprint arXiv:2106.14107},
  year   = {2021}
}

Comments

19 pages, 7 figures, 1 table. arXiv admin note: text overlap with arXiv:1504.02881 by other authors