An IBC enhanced DGTD scheme for transient analysis of EM interactions with graphene

A discontinuous Galerkin time-domain (DGTD) method is proposed for analyzing electromagnetic field interactions on graphene from microwave to terahertz frequencies. An impedance boundary condition (IBC) is utilized to model the graphene within the DGTD framework. The numerical flux is reformulated to take into account the IBC. Highly dispersive surface conductivity of graphene present in the resulting flux expression is approximated in terms of rational functions using the fast-relaxation vector-fitting technique. Via inverse Laplace transform, this facilitates the time domain matrix equations into an integral form for time variable t, finite integral technique (FIT) with recursive convolution method is employed to discrete and solve the matrix equations. The accuracy and applicability of the proposed IBC-DGTD is verified by numerical experiments.