Refine
Language
- English (4)
Keywords
- group velocity (2)
- high-order FEM (2)
- non-linear eigenvalue problem (2)
- photonic crystal wave-guide (2)
- Chebyshev interpolation (1)
- Dirichlet-to-Neumann map (1)
- Newton method (1)
- Robin-to-Robin map (1)
- Taylor expansion (1)
- adaptivity (1)
- band structure calculation (1)
- photonic crystals (1)
- super-cell method (1)
- surface modes (1)
Project
- D26 (4)
Application Area
- D (4)
The computation of guided modes in photonic crystal wave-guides is a key issue in the process of designing devices in photonic communications. Existing methods, such as the super-cell method, provide an efficient computation of well-confined modes. However, if the modes are not well-confined, the modelling error of the super-cell method becomes prohibitive and advanced methods applying transparent boundary conditions for periodic media are needed. In this work we demonstrate the numerical realization of a recently proposed Dirichlet-to-Neumann approach and compare the results with those of the super-cell method. For the resulting non-linear eigenvalue problem we propose an iterative solution based on Newton's method and a direct solution using Chebyshev interpolation of the non-linear operator. Based on the Dirichlet-to-Neumann approach, we present a formula for the group velocity of guided modes that can serve as an objective function in the optimization of photonic crystal wave-guides.
In this paper we present an efficient algorithm for the calculation of photonic crystal band structures and band structures of photonic crystal waveguides. Our method relies on the fact that the dispersion curves of the band structure are smooth functions of the quasi-momentum in the one-dimensional Brillouin zone. We show the derivation and computation of the group velocity, the group velocity dispersion, and any higher derivative of the dispersion curves. These derivatives are then employed in a Taylor expansion of the dispersion curves. We control the error of the Taylor expansion with the help of a residual estimate and introduce an adaptive scheme for the selection of nodes in the one-dimensional Brillouin zone at which we solve the underlying eigenvalue problem and compute the derivatives of the dispersion curves. The proposed algorithm is not only advantageous as it decreases the computational effort to compute the band structure but also because it allows
for the identification of crossings and anti-crossings of dispersion curves, respectively. This identification is not possible with the standard approach of solving the underlying eigenvalue problem at a discrete set of values of the quasi-momentum without taking the mode parity into account.
The efficient and reliable computation of guided modes in photonic crystal wave-guides is of great importance for designing optical devices. Transparent boundary conditions based on Dirichlet-to-Neumann operators allow for an exact computation of well-confined modes and modes close to the band edge in the sense that no modelling error is introduced. The well-known super-cell method, on the other hand, introduces a modelling error which may become prohibitively large for guided modes that are not well-confined. The Dirichlet-to-Neumann transparent boundary conditions are, however, not applicable for all frequencies as they are not uniquely defined and their computation is unstable for a countable set of frequencies that correspond to so called Dirichlet eigenvalues. In this work we describe how to overcome this theoretical difficulty introducing Robin-to-Robin transparent boundary conditions whose construction do not exhibit those forbidden frequencies. They seem, hence, well suited for an exact and reliable computation of guided modes in photonic crystal wave-guides.
In this work we present a complete algorithm for the exact computation of the guided mode band structure in photonic crystal (PhC) wave-guides. In contrast to the supercell method, the used approach does not introduce any modelling error and is hence independent of the confinement of the modes. The approach is based on Dirichlet-to-Neumann (DtN) transparent boundary conditions that yield a nonlinear eigenvalue problem. For the solution of this nonlinear eigenvalue problem we present a direct technique using Chebyshev interpolation that requires a band gap calculation of the PhC in advance. For this band gap calculation –- we introduce as a very efficient tool –- a Taylor expansion of the PhC band structure. We show that our algorithm –- like the supercell method –- converges exponentially, however, its computational costs –- in comparison to the supercell method –- only increase moderately since the size of the matrix to be inverted remains constant.