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.

Finite Difference Time Domain Methods for the Dirac Equation Coupled with the Chern–Simons Gauge Field

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

High-oder Symplectic FDTD Scheme for Solving Time-Dependent Schrdinger Equation

Finite-Difference Time-Domain Simulation of the Maxwell–Schrödinger System

A new conformal technique for FDTD (2, 4) scheme for modeling perfectly conducting composites

Simulating Maxwell–Schrödinger Equations by High-Order Symplectic FDTD Algorithm

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

Discontinuous Galerkin Method Using Laguerre Polynomials for Solving a Time-Domain Electric Field Integral Equation

Cauchy problem for Dirac equations with Chern-Simons-Proca gauge field

COMPUTATIONAL INVESTIGATION OF NANOSCALE SEMICONDUCTOR DEVICES AND OPTOELECTRONIC DEVICES FROM THE ELECTROMAGNETICS AND QUANTUM PERSPECTIVES BY THE FINITE DIFFERENCE TIME DOMAIN METHOD (INVITED REVIEW)

Development of an Novel Approach for the Simulation of Nanodevices Using FDTD

A Study of Space-Time Discretizations for the Dirac Equation

A Hybrid Approach for Solving Coupled Maxwell and Schrödinger Equations Arising in the Simulation of Nano-Devices

An Efficient Finite Difference/hermite-Galerkin Spectral Method for Time-Fractional Coupled Sine-Gordon Equations on Multidimensional Unbounded Domains and Its Application in Numerical Simulations of Vector Solitons

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

An Explicit Adaptive Finite Difference Method for the Cahn–Hilliard Equation

A novel high accuracy finite-difference time-domain method

Effective Mass of Composite Fermions and Fermionic Chern-Simons Theory in Temporal Gauge

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

Fundamental solutions for the Dirac equation in curved spacetime and generalized Euler-Poisson-Darboux equation

Exponential Collocation Methods Based on Continuous Finite Element Approximations for Efficiently Solving the Cubic Schrödinger Equation