Numerical study of the effect of defect layer on unmagnetized plasma photonic crystals

Abstract

In this paper an explicit finite-difference time domain scheme developed in staggered grids is used to solve the Maxwell's equations in Drude medium. Besides the preservation of discrete zero-divergence condition in electric and magnetic fields, we also aim to conserve the inherent conservation laws in simple medium all the time using the temporally second-order accurate explicit symplectic partitioned Runge-Kutta scheme. Within the framework of a semidiscretized method, the first-order spatial derivative terms in Faraday's and Ampere's equations are approximated to get an accurate numerical dispersion relation equation. The derived numerical angular frequency is accurately related to the wavenumber of Maxwell's equations for the space centered scheme of fourth-order accuracy. The resulting symplectic finite difference scheme developed in the time domain minimizes the difference between the exact and numerical group velocities. This newly proposed scheme is applied to model EM waves in the unmagnetized plasma crystal which contains a defect layer in photonic crystal. Our purpose is to numerically study the effects of defect layers on the propagation insight.