Nextgeneration active integrated optic subsystems
Information society technologies programme
of the European Commission
Project IST200028018, Workpackage 3: Design
Online collections of illustrative example results
Some noncommercial pieces of software that may be of interest for the participating groups
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, UTMESA^{+}
 
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 bentstraight 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, UTMESA^{+}
 

Mode solvers for bent waveguides and ringresonators in 2D and 3D,
related to the implementation of a 3D Film Mode Matching technique:
Ladislav Prkna, IREE Prague

The AAMP Finite Difference Time Domain program provides an environment
to perform twodimensional electromagnetic simulations on nondispersive,
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, UTMESA^{+}
 
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, UTMESA^{+}
 
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 guidedwave
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, UTMESA^{+}
 
A quasianalytic mode solver for dielectric integrated optical
waveguides with rectangular 2D cross sections (Wave Matching Method).
The
htmlmanual
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, UTMESA^{+}
 
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 gradedindex waveguide (propagation constants and mode field distributions). Complex refractive indices are allowed. Method used: transversal impedance method using RungeKutta 3th and 4th 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
 

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, UTMESA^{+}

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 BlochFloquet 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
Manfred Hammer (11.2005, adapted 11.2021)