We establish the error bounds of fourth‐order compact finite difference (4cFD) methods for the Dirac equation in the massless and nonrelativistic regime, which involves a small dimensionless parameter
inversely proportional to the speed of light. In this regime, the solution propagates waves with wavelength
in time and
...