We compare the long-time error bounds and spatial resolution of finite difference methods with different discretizations for Dirac equation small electromagnetic potentials characterized by ??(0,1] a dimensionless parameter. begin simple widely used time domain (FDTD) methods, establish rigorous them, which are valid up to at O(1/?). In estimates, we pay particular attention how errors depend e...