Fourth-order uniformly accurate integrators with long time near conservations for the nonlinear Dirac equation in the nonrelativistic regime

Abstract: In this paper, we propose two novel fourth-order integrators that exhibit uniformly high accuracy and long-term near conservations for solving the nonlinear Dirac equation (NLDE) in the nonrelativistic regime. In this regime, the solution of the NLDE exhibits highly oscillatory behavior in time, characterized by a wavelength of O($\varepsilon^{2}$) with a small parameter $\varepsilon>0$. To ensure uniform temporal accuracy, we employ a two-scale approach in conjunction with exponential integrators, utilizing operator decomposition techniques for the NLDE. The proposed methods are rigorously proved to achieve fourth-order uniform accuracy in time for all $\varepsilon\in (0,1]$. Furthermore, we successfully incorporate symmetry into the integrator, and the long-term near conservation properties are analyzed through the modulated Fourier expansion. The proposed schemes are readily extendable to linear Dirac equations incorporating magnetic potentials, the dynamics of traveling wave solutions and the two/three-dimensional Dirac equations. The validity of all theoretical ndings and extensions is numerically substantiated through a series of numerical experiments.