Super-resolution of time-splitting methods for the Dirac equation in the nonrelativistic regime

Abstract

We establish error bounds of the Lie-Trotter splitting (S1) and Strang splitting (S2) for the Dirac equation in the nonrelativistic limit regime in the absence of external magnetic potentials, with a small parameter 0<≤ 1 inversely proportional to the speed of light. In this limit regime, the solution propagates waves with O(2) wavelength in time. Surprisingly, we find out that the splitting methods exhibit super-resolution, in the sense of breaking the resolution constraint under the Shannon's sampling theorem, i.e. the methods can capture the solutions accurately even if the time step size τ is much larger than the sampled wavelength at O(2). S1 shows 1/2 order convergence uniformly with respect to , by establishing that there are two independent error bounds τ + and τ + τ/. Moreover, if τ is non-resonant, i.e. τ is away from certain region determined by , S1 would yield an improved uniform first order O(τ) error bound. In addition, we show S2 is uniformly convergent with 1/2 order rate for general time step size τ and uniformly convergent with 3/2 order rate for non-resonant time step size. Finally, numerical examples are reported to validate our findings.

0

Turn this paper into a lesson

ArcXiv compiles a structured reading guide from this paper's metadata: plain-English importance, contributions, prerequisite concepts, which sections to read first, flashcards, and a quiz. Grounded in the abstract, never invented.

Discussion (0)

Sign in to join the discussion.

Loading comments…