Stability and asymptotic properties of a linearized hydrodynamic medium model for dispersive media in nanophotonics

We analyze the stability of a linearized hydrodynamical model describing the response of nanometric dispersive metallic materials illuminated by optical light waves that is the situation occurring in nanoplasmonics. This model corresponds to the coupling between the Maxwell system and a PDE describing the evolution of the polarization current of the electrons in the metal. We show the well posedness of the system, polynomial stability and optimal energy decay rate. We also investigate the numerical stability for a discontinuous Galerkin type approximation and several explicit time integration schemes.

Discontinuous Galerkin Finite Element Methods for Maxwell\u27s Equations in Dispersive and Metamaterials Media

Discontinuous Galerkin Finite Element Method (DG-FEM) has been further developed in this dissertation. We give a complete proof of stability and error estimate for the DG-FEM combined with Runge Kutta which is commonly used in different fields. The proved error estimate matches those numerical results seen in technical papers. Numerical simulations of metamaterials play a very important role in the design of invisibility cloak, and sub-wavelength imaging. We propose a leap-frog discontinuous Galerkin Finite Element Method to solve the time-dependent Maxwell\u27s equations in metamaterials. The stability and error estimate are proved for this scheme. The proposed algorithm is implemented and numerical results supporting the analysis are provided. The wave propagation simulation in the double negative index metamaterials supplemented with perfectly matched layer (PML) boundary is given with one discontinuous Galerkin time difference method (DGTD), of which the stability and error estimate are proved as well in this dissertation. To illustrate the effectiveness of this DGTD, we present some numerical result tables which show the consistent convergence rate and the simulation of PML in metamaterials is tested in this dissertation as well. Also the wave propagation simulation in metamaterals by this DGTD scheme is consistent with those seen in other papers. Several techniques have appeared for solving the time-dependent Maxwell\u27s equations with periodically varying coefficients. For the first time, I apply the discontinuous Galerkin (DG) method to this homogenization problem in dispersive media. For simplicity, my focus is on obtaining a solution in two-dimensions (2D) using 2D corrector equations. my numerical results show the DG method to be both convergent and efficient. Furthermore, the solution is consistent with previous treatments and theoretical expectations

Stability and asymptotic properties of a linearized hydrodynamic medium model for dispersive media in nanophotonics

International audienceWe analyze the stability of a linearized hydrodynamical model describing the response of nanometric dispersive metallic materials illuminated by optical light waves that is the situation occurring in nanoplasmonics. This model corresponds to the coupling between the Maxwell system and a PDE describing the evolution of the polarization current of the electrons in the metal. We show the well posedness of the system, polynomial stability and optimal energy decay rate. We also investigate the numerical stability for a discontinuous Galerkin type approximation and several explicit time integration schemes.

Spectral Element Method Simulation of Linear and Nonlinear Electromagnetic Field in Semiconductor Nanostructures

<p>In this dissertation, the spectral element method is developed to simulate electromagnetic field in nano-structure consisting of dielectric, metal or semiconductor. The spectral element method is a special kind of high order finite element method, which has spectral accuracy. When the order of the basis function increases, the accuracy increases exponentially. The goal of this dissertation is to implement the spectral element method to calculate the electromagnetic properties of various semiconductor nano-structures, including photonic crystal, photonic crystal slab, finite size photonic crystal block, nano dielectric sphere. The linear electromagnetic characteristics, such as band structure and scattering properties, can be calculated by this method with high accuracy. In addition, I have explored the application of the spectral element method in nonlinear and quantum optics. The effort will focus on second harmonic generation and quantum dot nonlinear dynamics. </p><p>The electromagnetic field can be simulated in both frequency domain and time domain. Each method has different application for research and engineering. In this dissertation, the first half of the dissertation discusses the frequency domain solver, and the second half of the dissertation discusses the time domain solver.</p><p>For frequency domain simulation, the basic equation is the second order vector Helmholtz equation of the electric field. This method is implemented to calculate the band structure of photonic crystals consisting of dielectric material as well as metallic materials. Because the photonic crystal is periodic, only one unit cell need to be simulated in the computational domain, and a periodic boundary condition is applied. The spectral accuracy is inspected. Adding the radiation boundary condition at top and bottom of the computational region, the scattering properties of photonic crystal slab can be calculated. For multiple layers photonic crystal slab, the block-Thomas algorithm is used to increase the efficiency of the calculation. When the simulated photonic crystals are finite size, unlike an infinitely periodic system, the periodic boundary condition does not apply. In order to increase the efficiency of the simulation, the domain decomposition method is implemented. </p><p>The second harmonic generation, which is a kind of nonlinear optical effect, is simulated by the spectral element method. The vector Helmholtz equations of multiple frequencies are solved in parallel and the consistence solution with nonlinear effect is obtained by iterative solver. The sensitivity of the second harmonic generation to the thickness of each layer can be calculated by taking the analytical differential of the equation to the thickness of each element. </p><p>The quantum dot dynamics in semiconductor are described by the Maxwell-Bloch equations. The frequency domain Maxwell-Bloch equations are deduced. The spectral element method is used to solve these equations to inspect the steady state quantum dot dynamic behaviors under the continuous wave electromagnetic excitation.</p><p>For time domain simulation, the first order curl equations in Maxwell equations are the basic equations. A spectral element method based on brick element is implemented to simulate a nano-structure consisting of woodpile photonic crystal. The resonance of a micro-cavity consisting of a point defect in the woodpile photonic crystal block is simulated. In addition, the time domain Maxwell-Bloch equations are implemented in the solver. The spontaneous emission process of quantum dot in the micro-cavity is inspected. </p><p>Another effort is to implement the Maxwell-Bloch equations in a previously implemented domain decomposition spectral element/finite element time domain solver. The solver can handle unstructured mesh, which can simulate complicated structure. The time dependent dynamics of a quantum dot in the middle of a nano-sphere are investigated by this implementation. The population inversion under continuous and pulse excitation is investigated. </p><p>In conclusion, the spectral element method is implemented for frequency domain and time domain solvers. High efficient and accurate solutions for multiple layers nano-structures are obtained. The solvers can be applied to design nano-structures, such as photonic crystal slab resonators, and nano-scale semiconductor lasers.</p>Dissertatio