NAIS project

Next-generation active integrated optic subsystems

Information society technologies programme of the European Commission
Project IST-2000-28018, Workpackage 3: Design

Online collections of illustrative example results

Some noncommercial pieces of software that may be of interest for the participating groups

3D CMT model of circular microresonators A simulator for 3D circular microresonators based on numerically computed profiles for vectorial modes of bent and straight optical channels. A variant of coupled mode theory is used to evaluate the evolution of an arbitrary number of modes that couple as they propagate along a set of waveguides. The C++ implementation includes a reasonably comfortable user interface, in which the structure can be defined, calculations are done, and fields may be visualized. Mode field calculations are not fully integrated so far: Basis mode profiles are computed externally (FMM or OlympIOs) and imported from files in a discretized format. Two coupler configurations are combined with the bends of the ring or disk cavity to compute the wavelength response of the full microresonator. After a license agreement with the University of Twente, the company C2V will take care of the further distribution and support of the programs.
Please contact Remco Stoffer, UT-MESA+
2D CMT model of circular microresonators Circurs: A semianalytical 2D microresonator model. Analytical solutions for modes supported by bend segments and straight channels are combined in a framework of frequency domain coupled mode theory. The resonator is functionally decomposed into two bent-straight waveguide couplers, which are connected to each other by segments of bent waveguides. The C++ implementation provides objects for straight and bent waveguides, coupler configurations, and full microresonator structures. The latter includes methods for the evaluation of the power transmission, the free spectral range, the width of the resonances, and, for given vacuum wavelength, of the full (2D) optical fields. Scans over the wavelength parameter permit to compute the spectral response of the microresonator device. A preliminary version of an online manual is available, including the C++ code and a few example driver files.
Please contact Kiran Hiremath, UT-MESA+
2D multilayer slab mode solver
 
Mode solver for 2D bends and rings/disks
 
FMM mode solver for 3D step-index channels
 
FMM mode solver for 3D step-index bent channel waveguides
 
FMM mode solvers for 3D bends and rings/disks
Mode solvers for bent waveguides and ringresonators in 2D and 3D, related to the implementation of a 3D Film Mode Matching technique:
  • 2D mode solver for slab waveguides (multilayers). Accepts both lossless and lossy layers. Uses either transfer matrix, impedance or reflectance method and bracketing of effective indices for lossless and contour integral evaluation or mode tracking for lossy structures. The waveguide is enclosed with artificial magnetic or electric walls (optionally PMLs) - this discretization of the continuous spectrum allows to locate not only guided but radiation modes, too.
  • 2D mode solver for planar microresonators. Both microring and microdisk structures are supported. The propagation constant (here a natural number to meet the resonant condition) is estimated from the mode of a 2D bent waveguide with the same parameters. A contour line integral method is then used to search the complex plane of frequency for zeros of the dispersion equation, which holds for microresonators as well as for bent waveguides.
  • 3D mode solver for straight waveguides, semivectorial Film Mode Matching approach for lossless waveguides with step-index cross-section. Each slice (horizontally uniform segment - multilayer) is characterized by a set of eigenmodes. These are used to establish the dispersion relation horizontally. The structure is vertically enclosed with suitable boundary conditions to discretize the continuous spectrum of the multilayers.
  • 3D mode solver for bent waveguides, semivectorial Film Mode Matching approach for waveguides with step-index cross-section, provides complex root tracking of bend modes. Taking a straight waveguide mode as initial guess, a complex zero of the dispersion equation of a bent waveguide with long radius is polished. The stepwise decrease of the bend radius yields the mode evolution of a structure being bent.
  • 3D mode solver for microresonators, semivectorial Film Mode Matching approach for microrings and microdisks with step-index cross-section, intended as a different approach to microresonator modeling. If circular symmetry is considered, a microresonator is characterized as a lossy resonant structure with complex eigenfrequencies. Somewhat more time consuming as the slices' eigenmodes must be recomputed for each trial frequency.
Ladislav Prkna, IREE Prague
2D-FDTD The AAMP Finite Difference Time Domain program provides an environment to perform two-dimensional electromagnetic simulations on non-dispersive, linear, nonabsorbing and nonamplifying materials with arbitrary shapes. Perfectly Matched Layers (PML's) are employed to absorb radiation leaving the calculation domain. Input fields can either be generated from behind the PML's, or inside the window by means of the Total Field / Scattered Field approach. Analysis can be performed on e.g. waveguide mode powers in a wide spectral range; fields can be saved to files at arbitrary times. A preliminary version of the FDTD program for MS Windows and a brief user manual are available online.
Please contact Remco Stoffer, UT-MESA+
2D analytic bend mode solver A semianalytical 2D bend mode solver, implemented in C++, based on field matching at radial interfaces of different materials. In the present state, the mode solver is limited to structures with at most one guiding layer. For given real frequency / vacuum wavelength, the program computes complex valued mode profiles and the corresponding complex propagation constants, or effective mode indices and attenuation constants, respectively. Mode profiles are represented in terms of Bessel and Hankel functions with complex order, which are evaluated using their "uniform asymptotic expansions". A secant method is used to find complex roots that correspond to valid propagation constants. Further details are available in a research report. The programs have been included in the Circurs subroutine collection.
Please contact Kiran Hiremath, UT-MESA+
Metric - Mode expasion tool for rectangular integrated optical circuits Metric: Tools for semianalytic modeling in integrated optics / photonics, simulations of lossless 2D configurations with rectangular refractive index distributions. Includes a mode solver engine for dielectric multilayer slab waveguides, and solvers for guided-wave Helmholtz (-scattering) problems based on spectral discretizations along one (BEP: bidirectional eigenmode propagation) or two coordinate axes (QUEP: quadridirectional eigenmode propagation). A series of application examples, embedded in the html manual, complements the downloadable, commented C++ sources.
Manfred Hammer, UT-MESA+
WMM mode solver A quasianalytic mode solver for dielectric integrated optical waveguides with rectangular 2D cross sections (Wave Matching Method). The html-manual includes the downloadable, commented C++ sources, accompanied by several application examples. Although the programs are designed for straight waveguides only, the included features of perturbation theory may prove to be useful when dealing with lossy or anisotropic materials. A collection of additional application files extends the WMM coupled mode theory implementation to adiabatic, threedimensional directional couplers. Accepting the - admittedly crude - approximation of the bend mode of a ring cavity by a profile that corresponds to a straight waveguide with analogous cross section, these procedures can provide a computationally cheap means to tackle the problem of coupling the cavity loop in a cylindrical microresonator to a straight port waveguide. Several documents reporting on the evaluation of this model for the configurations considered in NAIS can be found in the workplace section of the NAIS website.
Manfred Hammer, UT-MESA+
1D mode analysis Modal analysis of multilayer and graded index waveguides with 1D cross sections. Specifically:
(1)  Modal analysis of a general 1D multilayer waveguide (propagation constants and mode field distributions). Complex refractive indices are allowed. In open waveguides, both guided and leaky modes can be calculated. Method used: transfer matrix method.
(2)  Modal analysis of a general 1D graded-index waveguide (propagation constants and mode field distributions). Complex refractive indices are allowed. Method used: transversal impedance method using Runge-Kutta 3-th and 4-th order integration.
(3)  Modal analysis of a general 1D multilayer waveguide (propagation constants and mode field distributions) of waveguides containing metal layers with negative real part of permittivity that support surface plasmons. Method used: transfer matrix method.
Please contact Jiri Ctyroky, IREE Prague
OMS - 1D multilayer slab waveguide mode solver - Online !
 
EIMS - 2D effective index mode solver - Online !
Online mode solvers (Java applets), directly executable in your web browser. A 1D multilayer slab waveguide mode solver and a 2D effective index multilayer waveguide mode solver are available. The programs are meant as a kind of pocket calculator for rough and quick assessment of problems involving slab waveguides, and for purposes of demonstration. The limited functionality is hardly sufficient for more extensive design tasks, where one should resort to other computational means. Analogous C++ programs (Metric) are available as well, without a graphical user interface, but applicable to a substantially larger variety of problems.
Manfred Hammer, UT-MESA+
Mode expansion simulations Uni- and bidirectional eigenmode propagation simulations for segmented planar waveguide structures. More specifically:
(1)  Modally resolved transmission and reflection in a 1D (planar) waveguide consisting of a concatenation of an arbitrary number of sections of longitudinally uniform multilayer waveguides. Complex refractive indices are allowed. Method used: Bidirectional mode expansion and propagation method with PMLs.
(2)  Modally resolved transmission in a 1D (planar) waveguide consisting of a concatenation of a several sections of longitudinally uniform multilayer waveguides. One section can contain metal layer supporting propagation of surface plasmons. Complex refractive indices are allowed. Method used: Unidirectional option of a bidirectional mode expansion and propagation method with PMLs.
(3)  Modally resolved transmission and reflection in a 1D (planar) waveguide consisting of a concatenation of periodically repeated longitudinally uniform sections of multilayer waveguides. Complex refractive indices are generally allowed. Method used: Bidirectional mode expansion and propagation method combined with the Bloch-Floquet theorem and PMLs.
Please contact Jiri Ctyroky, IREE Prague

For information on the C2V software products (OlympIOs), please consult the C2V website directly.

Selected research reports / drafts

Journal articles and contributions to conference proceedings related to WP3 of NAIS

... and two theses

email contact Manfred Hammer   (11.2005, adapted 11.2021)