Efficient HIE-FDTD method for designing a dual-band anisotropic terahertz absorption structure.

The finite-difference time-domain (FDTD) method is considered to be one of the most accurate and common methods for the simulation of optical devices. However, the conventional FDTD method is subject to the Courant-Friedrich-Levy condition, resulting in extremely low efficiency for calculating two-dimensional materials (2DMs). Recent researches on the hybrid implicit-explicit FDTD (HIE-FDTD) method show that the method can efficiently simulate homogeneous and isotropic 2DMs such as graphene sheet; however, it is inapplicable to the anisotropic medium. In this paper, we propose an in-plane anisotropic HIE-FDTD method to simulate optical devices containing graphene and black phosphorus (BP) sheets. Numerical analysis shows that the proposed method is accurate and efficient. With this method, we present a novel multi-layer graphene-BP-based dual-band anisotropic terahertz absorption structure (GBP-DATAS) and analyze its optical characteristics. Combining the advantages of graphene and BP localized surface plasmons, the GBP-DATAS demonstrates strong anisotropic plasmonic resonance and high absorption rate in the terahertz band.