Multiphysics Simulation of Plasmonic Photoconductive Antennas using Discontinuous Galerkin Methods
share
Plasmonic nanostructures significantly improve the performance of photoconductive antennas (PCAs) in generating terahertz radiation. However, they are geometrically intricate and result in complicated electromagnetic (EM) field and carrier interactions under a bias voltage and upon excitation by an optical EM wave. These lead to new challenges in simulations of plasmonic PCAs, which cannot be addressed by existing numerical frameworks. In this work, a multiphysics framework making use of discontinuous Galerkin (DG) methods is developed to address these challenges. The operation of the PCA is analyzed in stationary and transient states, which are described by coupled systems of the Poisson and stationary drift-diffusion (DD) equations and the time-dependent Maxwell and DD equations, respectively. Both systems are discretized using DG schemes. The nonlinearity of the stationary system is accounted for using the Gummel iterative method while the nonlinear coupling between the time-dependent Maxwell and DD equations is tackled during time integration. The DG-based discretization and the explicit time marching help in handling space and time characteristic scales that are associated with different physical processes and differ by several orders of magnitude. The accuracy and applicability of the resulting multiphysics framework are demonstrated via simulations of conventional and plasmonic PCAs.
READ FULL TEXT