Modeling open nanophotonic systems using the Fourier modal method: generalization to 3D Cartesian coordinates

Abstract

Recently, an open geometry Fourier modal method based on a new combination of an open boundary condition and a non-uniform k-space discretization was introduced for rotationally symmetric structures providing a more efficient approach for modeling nanowires and micropillar cavities [J. Opt. Soc. Am. A 33, 1298]. Here, we generalize the approach to three-dimensional (3D) Cartesian coordinates allowing for the modeling of rectangular geometries in open space. The open boundary condition is a consequence of having an infinite computational domain described using basis functions that expand the whole space. The strength of the method lies in discretizing the Fourier integrals using a non-uniform circular "dartboard" sampling of the Fourier k space. We show that our sampling technique leads to a more accurate description of the continuum of the radiation modes that leak out from the structure. We also compare our approach to conventional discretization with direct and inverse factorization rules commonly used in established Fourier modal methods. We apply our method to a variety of optical waveguide structures and demonstrate that the method leads to a significantly improved convergence enabling more accurate and efficient modeling of open 3D nanophotonic structures.

Benchmarking five numerical simulation techniques for computing resonance wavelengths and quality factors in photonic crystal membrane line defect cavities

“…The FMM with PMLs included is referred to as aperiodic FMM (aFMM). As a perspective for future FMM work, we note that an infinite-size computational domain version of FMM has recently been developed [55], in which absorbing boundary conditions are avoided.…”

High quality factor GaAs microcavity with buried bullseye defects

“…We use a recently introduced a numerical technique based on an open geometry formulation of the Fourier modal method [28,29]. Here, the geometry is divided into material sections uniform along the propagation z axis.…”

Open-geometry modal method based on transverse electric and transverse magnetic mode expansion for orthogonal curvilinear coordinates

“…In this work, a true open geometry boundary condition based on an innite domain combined with an equidistant discretization of the continuous k-space was implemented. More recently, non-uniform k-space discretization schemes for rotationally symmetric [25] and Cartesian coordinates [26] enabling signicantly improved convergence were proposed. However, the modeling of state-of-the-art singlephoton sources requires a modal method formalism allowing for orthogonal curvilinear coordinate systems, e.g.…”