A Resistive Boundary Condition Enhanced DGTD Scheme for the Transient Analysis of Graphene

In this paper, the electromagnetic (EM) features of graphene are characterized by a discontinuous Galerkin time-domain (DGTD) algorithm with a resistive boundary condition (RBC). The atomically thick graphene is equivalently modeled using an RBC by regarding the graphene as an infinitesimally thin conductive sheet. To incorporate RBC into the DGTD analysis, the surface conductivity of the graphene composed of contributions from both intraband and interband terms is first approximated by rational basis functions using the fast-relaxation vector-fitting (FRVF) method in the Laplace domain. Next, through the inverse Laplace transform, the corresponding time-domain matrix equations in integral can be obtained. Finally, these matrix equations are solved by time-domain finite integral technique (FIT). For elements not touching the graphene sheet, however, the well-known Runge-Kutta (RK) method is employed to solve the two first-order time-derivative Maxwell´s equations. The application of the surface boundary condition significantly alleviates the memory consuming and the limitation of time step size required by Courant-Friedrichs-Lewy (CFL) condition. To validate the proposed algorithm, various numerical examples are presented and compared with available references.