New Methods for Computing Transmission and Reflection Spectra of Photonic Crystals and Diffraction Gratings
Yumao Wu
2010-01-01
Abstract:Due to their unusual optical properties and significant potentials in applications, photonic crystals (PhCs) have been extensively studied both theoretically and experimentally. Many PhCs devices, such as waveguide bends, branches, frequency filters and waveguide couplers, are important building blocks of integrated optical circuits. Threedimensional (3D) PhCs with a complete bandgap in the optical wavelength region have potential applications, such as ultrahigh quality-factor cavities and zero-threshold lasers. The woodpile structures composed of alternating layers of rods have attracted much attention due to their relatively simple fabrication process compared with other 3D PhCs. Diffraction gratings have been extensively studied for many years, and they are important in practical applications, such as monochromators, spectrometers, lasers, wavelength division multiplexing devices, optical pulse compressing devices, etc. Numerical methods are essential in the design, analysis and optimization of photonic crystals and diffraction gratings. Some of these methods, such as the finitedifference time-domain (FDTD) method, are general methods that can be used to study various aspects of PhCs and diffraction gratings, but their accuracy and efficiency are often limited. FDTD requires a small mesh size to resolve curved material interfaces and often has difficulties truncating periodic structures that extend to infinity. The Fourier modal method (FMM) is suitable for diffraction gratings with uniform layers, but they are not so efficient and may have convergence problems when a general grating with sloping interfaces or a photonic crystal composed of cylinders must be approximated. The boundary integral equation (BIE) method is somewhat complicated to implement, since the integral operators are related to the quasiperiodic Green’s function which requires sophisticated lattice sums techniques to evaluate. The finite element method (FEM) is very general, but it gives rise to large, complex and indefinite linear systems that are expensive to solve. In this thesis, efficient numerical methods based on the Dirichlet-to-Neumann (DtN) or Neumann-to-Dirichlet (NtD) maps (of unit cells or homogeneous sub-domains) are developed for accurate simulations of two- and three-dimensional photonic crystals and diffraction gratings. For photonic crystals composed of interpenetrating cylinders, i.e. when the radius of the cylinders is larger than √3/4 of the lattice constant, an efficient DtN map method is developed for computing the reflection and transmission spectra. Our method manipulates a pair of operators defined on a set of curves. It is efficient since the wave field in the interiors of the unit cells are never calculated. This is achieved by using the DtN maps which map the wave field on the boundaries of the unit cells to its normal derivative. The DtN map method is also developed for two-dimensional photonic crystals with oblique incident waves. In that case, the DtN operator maps the two…