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.

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

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

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

Optimization of the Imaginary Time Step Evolution for the Dirac Equation

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

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

High-order Symplectic FDTD Scheme for Solving a Time-Dependent Schrödinger Equation

SYMPLECTIC PSEUDOSPECTRAL TIME-DOMAIN SCHEME FOR SOLVING TIME-DEPENDENT SCHRODINGER EQUATION

A Discrete Time-Dependent Method for Metastable Atoms and Molecules in Intense Fields

Perturbative Theory and Quantum Rate Equation of Time-Dependent Double-Well Transport

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

The Dirac Vacuum in Discrete Spacetime

Time-Dependent Quantum Transport: An Efficient Method Based On Liouville-Von-Neumann Equation For Single-Electron Density Matrix

Parallel transport dynamics for mixed quantum states with applications to time-dependent density functional theory

A combined R-matrix eigenstate basis set and finite-differences propagation method for the time-dependent Schr\"{od}dinger equation: the one-electron case

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

Solving coupled Non-linear Schrödinger Equations via Quantum Imaginary Time Evolution

A Dirac sea pilot-wave model for quantum field theory

Perturbative solution approach for computing the two-photon Kapitza-Dirac effect in a Gaussian beam standing light wave

Application of high-order SF-SFDTD scheme for solving a time-dependent Schrödinger equation

A Mathematical Analysis of Dirac Equation Physics