Continuous-discontinuous Galerkin time domain (CDGTD) method with generalized dispersive material (GDM) model for computational photonics

The discontinuous Galerkin time domain (DGTD) method and its recent flavor, the continuous-discontinuous Galerkin time domain (CDGTD) method, have been extensively applied to simulations in the radio frequency (RF) and microwave (MW) regimes due to their inherent ability to efficiently model multiscale problems. We propose to extend the CDGTD method to nanophotonics while exploiting its advantages which have already been established in the RF and MW regimes, such as domain decomposition, non-conformal meshing, high-order elements, and hp-refinement. However, at optical frequencies many materials are highly dispersive, so the modeling of nanophotonic devices requires accurate handling of different dielectric functions, including those of plasmonic elements, dielectrics, and tunable materials. In this paper, we propose a CDGTD method that incorporates a generalized dispersive material (GDM) model which is an efficient way to implement a wide range of optical dispersion models with a universal analytic function. Physics-based dispersion models, such as the Drude, Debye, Lorentz, and critical points as well as more complicated behavior founded on ab-initio principles can all be obtained as special cases of the universal GDM approach. The accuracy and convergence of this GDM-incorporated CDGTD are verified by numerical examples. The CDGTD method, equipped with the GDM model, paves the way to the efficient design and optimization of large scale photonic devices with a diverse range of optical dispersive materials.