Uniform error bounds of time-splitting methods for the nonlinear Dirac equation in the nonrelativistic limit regime
share
Super-resolution of the Lie-Trotter splitting (S_1) and Strang splitting (S_2) is rigorously analyzed for the nonlinear Dirac equation without external magnetic potentials in the nonrelativistic limit regime with a small parameter 0<ε≤ 1 inversely proportional to the speed of light. In this limit regime, the solution highly oscillates in time with wavelength at O(ε^2) in time. The splitting methods surprisingly show super-resolution, in the sense of breaking the resolution constraint under the Shannon's sampling theorem, i.e. the methods can capture the solution accurately even if the time step size τ is much larger than the sampled wavelength at O(ε^2). Similar to the linear case, S_1 and S_2 both exhibit 1/2 order convergence uniformly with respect to ε. Moreover, if τ is non-resonant, i.e. τ is away from certain region determined by ε, S_1 would yield an improved uniform first order O(τ) error bound, while S_2 would give improved uniform 3/2 order convergence. Numerical results are reported to confirm these rigorous results. Furthermore, we note that super-resolution is still valid for higher order splitting methods.
READ FULL TEXT