On numerical methods for the semi-nonrelativistic limit system of the nonlinear Dirac equation

Solving the nonlinear Dirac equation in the nonrelativistic limit regime numerically is difficult, because the solution oscillates in time with frequency of Oε-2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\mathscr {O}} \! \left( \varepsilon ^{-2}\right) $$\end{document}, where 0<ε≪1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$0<\varepsilon \ll 1$$\end{document} is inversely proportional to the speed of light. Yongyong Cai and Yan Wang have shown, however, that such solutions can be approximated up to an error of Oε2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\mathscr {O}} \! \left( \varepsilon ^2\right) $$\end{document} by solving the semi-nonrelativistic limit system, which is a non-oscillatory problem. For this system, we construct a two-step method, called the explicit exponential midpoint rule, and prove second-order convergence of the semi-discretization in time. Furthermore, we construct a benchmark method based on standard techniques and compare the efficiency of both methods. Numerical experiments show that the new integrator reduces the computational costs per time step to 40% and within a given runtime improves the accuracy significantly.