Maximization of the quality factor of an optical resonator

Abstract We consider resonance phenomena for the scalar wave equation in an inhomogeneous medium. Resonance is a solution to the wave equation which is spatially localized while its time dependence is harmonic except for decay due to radiation. The decay rate, which is inversely proportional to the qualify factor, depends on the material properties of the medium. In this work, the problem of designing a resonator which has high quality factor (low loss) is considered. The design variable is the index of refraction of the medium. High quality resonators are desirable in a variety of applications, including photonic band gap devices. Finding resonance in a linear wave equation with radiation boundary condition involves solving a nonlinear eigenvalue problem. The magnitude of the ratio between real and imaginary part of the eigenvalue is proportional to the quality factor Q . The optimization we perform is finding a structure which possesses an eigenvalue with largest possible Q . We present a numerical approach for solving this problem. The method consists of first finding a resonance eigenvalue and eigenfunction for a non-optimal structure. The gradient of Q with respect to index of refraction at that state is calculated. Ascent steps are taken in order to increase the quality factor Q . We demonstrate how this approach can be implemented and present numerical examples of high Q structures.

[1]  Alexander Figotin,et al.  Localization of light in lossless inhomogeneous dielectrics , 1998 .

[2]  A. Majda,et al.  Absorbing boundary conditions for the numerical simulation of waves , 1977 .

[3]  J. Gillis,et al.  Classical dynamics of particles and systems , 1965 .

[4]  Anne Greenbaum,et al.  Iterative methods for solving linear systems , 1997, Frontiers in applied mathematics.

[5]  Yu Chen,et al.  High-order corrected trapezoidal quadrature rules for functions with a logarithmic singularity in 2-D , 2002 .

[6]  F. Tisseur Backward error and condition of polynomial eigenvalue problems , 2000 .

[7]  Fadil Santosa,et al.  Optimal Localization of Eigenfunctions in an Inhomogeneous Medium , 2004, SIAM J. Appl. Math..

[8]  D. Larkman,et al.  Photonic crystals , 1999, International Conference on Transparent Optical Networks (Cat. No. 99EX350).

[9]  Yu Chen,et al.  A high-order, fast algorithm for scattering calculation in two dimensions , 2004 .

[10]  J. Zhang,et al.  An Approximate Method for Scattering by Thin Structures , 2005, SIAM J. Appl. Math..

[11]  T. Asano,et al.  High-Q photonic nanocavity in a two-dimensional photonic crystal , 2003, Nature.

[12]  Wolfgang Porod,et al.  An eigenvalue method for open‐boundary quantum transmission problems , 1995 .

[13]  Axel Ruhe ALGORITHMS FOR THE NONLINEAR EIGENVALUE PROBLEM , 1973 .

[14]  J. Keller,et al.  Exact non-reflecting boundary conditions , 1989 .

[15]  Karl Meerbergen,et al.  The Quadratic Eigenvalue Problem , 2001, SIAM Rev..

[16]  Kostas D. Kokkotas,et al.  Quasi-Normal Modes of Stars and Black Holes , 1999, Living reviews in relativity.

[17]  Steven G. Johnson,et al.  Photonic Crystals: The Road from Theory to Practice , 2001 .

[18]  Jean-Pierre Berenger,et al.  A perfectly matched layer for the absorption of electromagnetic waves , 1994 .

[19]  Andrew J. Majda,et al.  Absorbing Boundary Conditions for Numerical Simulation of Waves , 1977 .