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.

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

One-dimension Periodic Potentials in Schrödinger Equation Solved by the Finite Difference Method

Solutions of the Dirac equation in one fixed and one moving wall well

Delta-Function Spin-1/2 Fermions in A One-Dimensional Potential Well

Analytical solutions of the Dirac equation using the Tridiagonal Representation Approach: General study, limitations, and possible applications

Analytic Continuation in the Coupling Constant Method for the Dirac Equation

Explicit summation of the constituent WKB series and new approximate wave functions

Bound states of the one-dimensional Dirac equation for scalar and vector double square-well potentials

Absorbing Boundary Conditions for Variable Potential Schrödinger Equations via Titchmarsh-Weyl Theory

High-accuracy solution of the Dirac equation for three-dimensional hydrogen atom using a variationally optimized diagonalization method

Bound states of a one-dimensional Dirac equation with multiple delta-potentials

The Energy Eigenvalue for the Singular Wave Function of the Three Dimensional Dirac Delta Schrodinger Potential via Distributionally Generalized Quantum Mechanics

Exact Solution of the Problem of a Two-Dimensional Dirac Particle with Flux 1/2.

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

A deep neural network approach to solve the Dirac equation

One-dimensional quantum scattering from multiple Dirac delta potentials: A Python-based solution

A Study of Space-Time Discretizations for the Dirac Equation

Large time well-posedness for a Dirac–Klein–Gordon system

Employing an operator form of the Rodrigues formula to calculate wavefunctions without differential equations

Gaussian Beam Methods for the Dirac Equation in the Semi-classical Regime