▎ 摘　　要

In this paper, we present a method to incorporate the intraband and interband terms of the surface conductivity of graphene into the finite-difference time-domain (FDTD) method. The method is based on approximating the surface resistivity of graphene by a series of partial fractions in terms of real or complex conjugate pole-residue pairs. Then, a discrete time-domain surface boundary condition at the graphene sheet is generated, which is then incorporated within the FDTD method using the infinitesimally thin sheet formulation. Numerical examples are presented to validate and demonstrate the capabilities and advantages of the proposed approach.