▎ 摘　　要

A two-step locally one-dimensional finite-difference time-domain method (LOD-FDTD) is developed for modeling a graphene sheet biased by a magnetostatic field. The graphene sheet is considered as a current source characterized by two coupled auxiliary equations. The unconditional stability of the proposed algorithm is analyzed theoretically and verified numerically. Numerical experiments show that the proposed method can be used to efficiently capture the transmission properties of graphene sheets with a good accuracy. Further, the proposed method is applied to investigate the properties of the Faraday and Kerr rotations.