BendS   1-D mode solver for slab waveguide bends
 
---

— (under construction) —

Online bend mode solver

A mode solver for bent integrated optical dielectric multilayer slab waveguides and curved dielectric interfaces with 1-D cross sections. Following the structure definition in terms of bend radius, refractive indices, layer thicknesses, if applicable, and the vacuum wavelength, the script calculates the complex effective indices of the (leaky) modes supported by the bend, or their phase propagation and attenuation constants, respectively, and allows to inspect the corresponding optical field patterns. It is intended as a basic tool for integrated optics design, in particular for purposes of demonstration.

Input

For a curved slab waveguide with N intermediate layers, the input mask receives the vacuum wavelength, a specification of the polarization, the radius R of the outer rim of the bend, further refractive index values ni (interior region), n1, ... , nN (inner layers 1 to N), ne (exterior region), and thicknesses t1, ... , tN of the intermediate layers. Select N=0 to specify a curved dielectric interface only. All dimensions are meant in micrometers. The figure illustrates the relevant geometry:

Bent multilayer waveguide geometry

Polar coordinates r, θ span the Cartesian x-z-plane; the center of the bend is located at the origin of both coordinate systems. Light propagates along the angular coordinate θ. The refractive index profile is independent of θ and piecewise constant in the radial direction r. All electromagnetic fields and the refractive index distribution are assumed to be constant along the y-axis (perpendicular to the x-z- and r-θ-plane, not shown).

Computational parameters

The solver relies on a heuristic procedure for generating initial guesses for the roots of the transverse resonance condition in the complex plane. A complex secant method then converges the initial guesses numerically to actual roots. Further heuristics are applied for the classification and ordering of these roots.

The solver accepts a few parameters that influence this search for valid mode wavenumbers in the complex plane. If the numbers of trials is not directly specified (Auto-checkbox), the solver applies a heuristic based on modal solutions for an equivalent "tilted" effective refractive index profile for a straight graded-index slab. Alternatively, a number of initial guesses on a regular rectangular grid in the complex plane can be specified explicitly.

The range in the complex plane examined (initial guesses) for effective mode indices can be restricted to real parts larger than the smallest refractive index in the structure. The solver then responds with the modes of lowest radial order only (checkbox). Alternatively, an interval from close-to-zero up to the largest refractive index in the structure is examined. Then also modes of higher radial order can be identified. In all cases, attenuation constants in a range between 10-14 and 10-1 are considered.

While the found solutions should be valid bend modes, one must not rely too much on the solver being able to find all existing bend modes, nor on the correct classification of the modes. In many cases, however, these heuristics have been observed to work adequately. The solver should reproduce the results of Ref. [1] with (almost) no changes to the default computational parameters. You might wish to try this out as a test of the solver, and to become familiar with the interface.

Output

A table shows, for each bend mode, the complex bend mode eigenvalue in the form of either

Relations γ = k Neff and ν = γ R apply, where k = 2π/λ is the vacuum wavenumber associated with the specified vacuum wavelength λ. The mode identifier (order) indicates the number of nodes, in the region r < R, in the real part of the principal electric component Ey of the radial profile of TE modes, and in the principal magnetic component Hy of the radial profile of TM modes (where the phase of the mode profile has been adjusted such that the principal components are real at r = R).

One should be aware that the definition of the bend radius R is, in some respect, arbitrary. For a curved slab with one intermediate layer of thickness t1, a frequent choice is to describe the bend with a radius R' = R - t1/2 at the center of the slab (cf. the figure), rather than at the outer rim, as in case of the present solver. This choice should lead to the same solutions, i.e. to the same angular dependence ∼exp(-i ν θ) of the bend modes. These are then characterized by different values γ' = γ R/R' and Neff' = Neff R/R' of propagation constant and effective mode index, respectively.

Referring to the polar coordinate system as introduced above, this concerns optical electric fields E and magnetic fields H that depend on the spatial coordinates r, θ (these span the x-z-plane of interest), and on time t, with angular frequency ω, as

E(r, θ, t) = Re{E(r) exp(i ω t-i γ R θ)},   H(r, θ, t) = Re{H(r) exp(i ω t-i γ R θ)}.

All fields are constant along the y-direction. The profiles of TE bend modes are of the form E(r) = (0, Ey, 0)(r), and H(r) = (Hr, 0, Hθ)(r), where E and H are the electric and magnetic parts of the mode profile, respectively, depending on the radial coordinate r only. Likewise, the profiles of TM bend modes can be written as E(r) = (Er, 0, Eθ)(r), and H(r) = (0, Hy, 0)(r).

Being solutions of eigenvalue problems, the mode profiles are determined up to some complex constant only. No units are shown for their electric or magnetic fields. Still the given values correspond to a (rough) normalization of the modes to unit angular power flow: The radial integral of the angular component Sθ of the Poynting vector evaluates to 1 W/µm (power per lateral (y) unit length). Correspondingly, all electric fields are given in units of V/µm, magnetic fields are measured in A/µm, the components of the Poynting vector S have units of V·A/µm2 = W/µm2, and the electromagnetic energy density w is measured in W·fs/µm3. In this context the vacuum permittivity and permeability, respectively, are ε0 = 8.85·10-3 A·fs/(V·µm) and µ0 = 1.25·103 V·fs/(A·µm).

Mode profile plots show the real and imaginary parts of the complex field profile, its absolute value, or the absolute squared profile. The background shading indicates the dielectric structure, where darker shading means higher refractive index. After selecting "Plot", the extent of the vertical axis is being adjusted such that it covers the maximum values, determined separately for the electric field strength, magnetic field strength, Poynting vector, and the energy density, over all modes (and all their field components) that have been identified by the solver, on a default radial range. This is to make the plots comparable. Select the button labeled "↕" to adjust the vertical plot range to the functions that are actually displayed.

Single components ψ ∈ {Er, …, Hθ} of the complex-valued electromagnetic mode profile relate to time-varying physical fields Ψ(r, t) = Re ψ(r) exp(i ω t) at θ = 0. The animations show the respective component at recurring points in time, equally distributed over the period of 2π/ω.

Alternatively, propagation plots can be displayed. These show a component of the total physical electromagnetic field E, H on a rectangular domain in the Cartesian x-z-plane, for a unit (initial) amplitude of the bend mode at θ = 0. The global phase of the solution can be adjusted; the animations show the respective component at recurring points in time, equally distributed over the period of 2π/ω. Note that these are valid solutions of Maxwells equations for an angular segment of a bent waveguide only, not for an entire circular cavity; the fields exhibit discontinuities in phase and amplitude (attenuation) at θ = ±π.

Reference

[1] K.R. Hiremath, M. Hammer, R. Stoffer, L. Prkna, J. Ctyroky
Analytical approach to dielectric optical bent slab waveguides
Optical and Quantum Electronics 37 (1-3), 37-61 (2005)  Preprint (ps.gz)  Preprint (pdf)  External online source