“…For the numerical part, many efficient and accurate numerical methods have been proposed and analyzed, such as the finite difference time domain (FDTD) methods and spectral methods for the Klein-Gordon (KG) equation [9,34,42], the Dirac equation [6,7,28,32,33,35,37,48], the Klein-Gordon-Schrödinger equations [10,17,18] and the Maxwell-Dirac systems [33], and so forth. Recently, some numerical methods for the KGD Equation (1) have been proposed in the literature [15,36,38,49,50]. For the case 𝜀 = 1, some authors studied FDTD schemes and proved that they are second order in both time and space [15,49].…”