WMM based coupled mode theory: Adiabatic directional couplers
MESA+ Research Institute, University of Twente
(2001-2005, adapted 11.2021)
This page lists a collection of WMM application files that address the problem of guided light propagation through a straight (port) waveguide that is evanescently (vertically or horizontally) coupled to a ring (cavity) waveguide segment. See Ref.  for more details about the coupler geometry, about the implemented algorithms, and for results of the computations for a specific set of parameters. Ref.  summarizes how the properties of the directional couplers are to be embedded in a ringresonator model.
Top view of the directional coupler configuration: A segment of a cylindrical waveguide core (WG2) with radius R, evanescently coupled to a straight waveguide (WG1). Letters A, B, a, b denote the ports of the device, p(z) is the local y-displacement of the ring rim at position z.
Coupler cross section: The port waveguide with core width s, thickness v, and refractive index ng is embedded in a medium with refractive index nb, buried at a distance b below the surface x=0 of that medium. The cavity strip of width w and thickness t with refractive index nr is placed on top of the surface, covered by a material with refractive index na. The horizontal position of the port waveguide is defined by the y-coordinate g (positive or negative values are to be considered) of the core center. p is the horizontal, z-dependent displacement of the ring rim.
Details about the WMM simulation framework can be found in the html manual, including instructions for download, compilation, and execution of the C++ sources provided. The application files are only meant as initial protoypes. The paragraphs on waveguide definition give guidelines how to adapt the programs to other cross section geometries.
Port waveguide, port.c:
Mode properties of an embedded channel waveguide with rectangular core. Program output: Propagation constants / effective mode indices (
Ring waveguide, ring.c:
Mode properties of a rectangular raised strip waveguide. Program output: Propagation constants / effective mode indices (
Ring waveguide, width dependence, ringdisp.c:
Propagation constants versus the width of a rectangular raised strip waveguide. Program output: Propagation constants (
Coupler simulation, cpl.c:
Coupled mode theory analysis of the guided wave propagation along the composite coupler structure. Program output: Properties of the basis modes: Propagation constants and effective indices (
Coupler simulation, parameter dependences, cplsc.c:
CMT analysis of the coupler structure, scattering matrices versus the relative core position. Program output: Basis mode properties: Propagation constants and effective indices (
Approximations and assumptions
Some of the approximations and assumptions the coupled mode approach is based on are mentioned in Ref. . Additionally, the recipes that were originally derived and implemented for composite, longitudinally invariant structures (see Ref.  for details) are here applied to a device with z-dependent refractive index profile. Reasonable results can thus be expected only for an intermediate range of core distances (more precise: observations in ) and for relatively large ring radii (do not forget to adjust Z accordingly). The specific form of the couplers allows in a certain respect to check the consistency of the simulations:
paB+pab.xyfas produced by
cpl.cshould contain only unity values in the second column The same should hold for the sum of the scattering coefficients in the file
Tin the output streams of
cplsc.cought to be symmetric.
Necessarily the simulations satisfy these conditions only approximately. One should not trust the results to a degree that is higher than the accuracy observed in the above relations.
Strictly speaking, the "power coupling" matrix should appear in the expressions that evaluate the modal power. This becomes relevant for local configurations with "strong coupling", with pronounced mode overlaps. Hence, in particular for intermediate positions (z around 0), the interpretation of the squares of the individual amplitudes of the basis modes |Cj(z)|2 (and of their sum) is not straightforward. Nevertheless, assuming "decoupled" waveguides at z=-Z and at z=Z, the interpretation of the squared amplitudes as guided output powers is adequate.
For specific configurations, the relations may be pronouncedly violated, as a consequence of the approximations inherent in the present CMT version. However, in these cases we have observed the deviations bewtween the two offdiagonal elements of the coupler scattering matrix to be of the same order of magnitude as the violation of power conservation, where the average of the two coupling coefficients fits relatively well to the - almost coinciding - transmission coefficients, the diagonal elements of the scatter matrix. The programs cpl.c and cplsc.c therefore declare the average of the offdiagonal elements as the approximation of the coupling coefficient. Particular caution is advisable when using these results.
stderroutput stream. In contrast to
stderraccess is unbuffered, hence the data is more likely to survive an unintended end of the program execution. The sections about filenames and the log information comment about the types of output produced by the WMM programs.
?s replaced by characters that identify the simulation type) that records the evaluations of the least-squares error function the WMM approach relies upon. In particular for the dispersion calculation many of these files appear; delete these after a successful completion of the programs.
CalcBasisModes 0, it should be made sure that the relevant mode definition files are present in the directory where the program is executed.
Relying completely on the properties of straight waveguides, these programs have to be complemented by results from a suitable bend mode solver when aiming at a 3D model for the microring-resonators that are investigated in the NAIS project. This concerns the mode attenuation in the curved cavity waveguides in the first place. In a second place, available bend mode properties can prospectively help to improve a coupler design found by means of the WMM CMT programs: Given the mode profile and propagation constant of the straight basis field for the ring core, one could try to modify the core geometry and position, such that the bend mode calculated for the new core with the appropriate radius coincides as far as possible with the original straight field. This concerns the position of the field maximum, adjustable by a slight shift of the relative core positions, and the effective mode index, which can be matched e.g. by a change in the cavity core width.
Coupled mode model for 3D directional couplers,
NAIS internal note, September 2002.
 Standard ringresonator model, NAIS internal note, April 2003.
 M. Lohmeyer, N. Bahlmann, O. Zhuromskyy, and P. Hertel. Radiatively coupled waveguide polarization splitter simulated by wave-matching based coupled mode theory. Optical and Quantum Electronics 31, 877-891 (1999).