Abstract:We analyze rigorously error estimates and compare numerically spatial/temporal resolution of various numerical methods for the discretization of the Dirac equation in the nonrelativistic limit regime, involving a small dimensionless parameter 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} which is inversely proportional to the speed of light. In this limit regime, the solution is highly oscillatory in time, i.e. there are propagating waves with wavelength O(ε2)\documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$O(\varepsilon ^2)$$\end{document} and O(1) in time and space, respectively. We begin with several frequently used finite difference time domain (FDTD) methods and obtain rigorously their error estimates in the nonrelativistic limit regime by paying particular attention to how error bounds depend explicitly on mesh size h and time step τ\documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\tau $$\end{document} as well as the small parameter ε\documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\varepsilon $$\end{document}. Based on the error bounds, in order to obtain ‘correct’ numerical solutions in the nonrelativistic limit regime, i.e. 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}, the FDTD methods share the same ε\documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\varepsilon $$\end{document}-scalability on time step and mesh size as: τ=O(ε3)\documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\tau =O(\varepsilon ^3)$$\end{document} and h=O(ε)\documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$h=O(\sqrt{\varepsilon })$$\end{document}. Then we propose and analyze two numerical methods for the discretization of the Dirac equation by using the Fourier spectral discretization for spatial derivatives combined with the symmetric exponential wave integrator and time-splitting technique for temporal derivatives, respectively. Rigorous error bounds for the two numerical methods show that their ε\documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\varepsilon $$\end{document}-scalability is improved to τ=O(ε2)\documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\tau =O(\varepsilon ^2)$$\end{document} and h=O(1)\documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$h=O(1)$$\end{document} when 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}. Extensive numerical results are reported to support our error estimates.

First Attempt to Overcome the Disaster of Dirac Sea in Imaginary Time Step Method

Avoid the Tsunami of the Dirac Sea in the Imaginary Time Step Method

Solving the Dirac equation with nonlocal potential by Imaginary Time Step method

Invariant-Related Unitary Transformation Method and Exact Solutions for the Quantum Dirac Field in A Time-Dependent Spatially Homogeneous Electric Field

An Improved MOD Discontinuous Galerkin Method for Solving Time-Domain Electric Field Integral Equation

Optimization of the Imaginary Time Step Evolution for the Dirac Equation

Convergence for Imaginary Time Step Evolution in the Fermi and Dirac Seas

Solving Dirac Equations on A 3d Lattice with Inverse Hamiltonian and Spectral Methods

Solving the relativistic Hartree-Bogoliubov equation with the finite-difference method

The Extended Wronskian Determinant Approach and the Iterative Solutions of One-Dimensional Dirac Equation

Real Stabilization Method for Single Particle Resonances

Improved uniform error estimates for the two-dimensional nonlinear space fractional Dirac equation with small potentials over long-time dynamics

Real Stabilization Method and Its Application in Calculating Resonance States

Optimized Dirac Woods-Saxon basis for covariant density functional theory

Numerical Methods and Comparison for the Dirac Equation in the Nonrelativistic Limit Regime

A Hybrid Explicit-Implicit Scheme for the Time-Dependent Wigner Equation

Simulation of Zitterbewegung by Modelling the Dirac Equation in Metamaterials

Bohm'S Formulation in Imaginary Time: Estimation of Energy Eigenvalues

Symplectic Interior Penalty Discontinuous Galerkin Method for Solving the Seismic Scalar Wave Equation

Analytic Continuation in the Coupling Constant Method for the Dirac Equation

Riemann-Hilbert problem for the defocusing Lakshmanan-Porsezian-Daniel equation with fully asymmetric non-zero boundary conditions