An online solver for the whispering gallery resonances of microrings, microdisks, or more general 2-D circular dielectric multilayer cavities in integrated optics / photonics. Following the structure definition in terms of cavity radius, refractive indices, and layer thicknesses, if applicable, and the specification of a target vacuum wavelength or of an interval of potential resonance wavelengths, the script calculates the complex angular eigenfrequencies of the (leaky) whispering gallery modes (WGMs) supported by the cavity, and their resonance wavelengths and linewidths, resonance frequencies, cavity lifetimes, and quality (Q-) factors, respectively. Facilities for detailed inspection of the corresponding mode profiles and resonant field patterns, and for exporting data and figures, are provided.


For a circular cavity with N intermediate layers, the input mask receives the radius R of the outer rim of the cavity, a specification of the polarization, 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 uniform circular dielectric disk. All dimensions are meant in micrometers. The figure illustrates the relevant geometry:

Circular dielectric multilayer cavity

Polar coordinates r, θ span the Cartesian x-z-plane; the center of the cavity is located at the origin of both coordinate systems. 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 tries to identify whispering gallery modes (WGMs) with resonance wavelengths that are either close to a target vacuum wavelength, or fall within a specific interval of vacuum wavelengths (Interval-checkbox). A somewhat heuristic procedure is applied for generating initial guesses for the roots of the transverse resonance condition in the complex plane. A complex secant method converges these 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 influences the process to determine the integer angular mode orders that relate to the resonances of interest. If the numbers of trials is not directly specified (Auto-checkbox), the solver applies a heuristic based on bend mode solutions for an equivalent "tilted" effective refractive index profile for a straight graded-index slab. Alternatively, a number of initial guesses for regularly spaced angular mode numbers and attenuation levels can be specified explicitly.

The range in the complex plane examined (initial guesses) for angular mode orders can be restricted to values that corresponded to angular light propagation with effective indices 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-2 are considered.

While the found solutions should be valid whispering gallery modes, one must not rely too much on the solver being able to find all existing WGMs, nor on the correct classification of the modes. In many cases, however, these heuristics have been observed to work adequately. The solver should reproduce those results of Ref. [1] that concern resonances of separate, uncoupled rings or disks 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.


A table shows, for each whispering gallery mode, the complex WGM eigenvalue in the form of either

Relations τ = 1/(2 α), Q = ω/(2 α) = λ/(Δλ) , λ = 2 π c/ω, Δω = 2 α = ω/Q = 1/τ, T = 2 π/ω, and f = 1/T = ω/(2 π) apply, where c is the vacuum speed of light.

The WGM identifiers comprise two indices. A first index indicates the radial order of the mode, determined as 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). The second index gives the angular order m of the WGM, the number of full exponential cycles in the principal fields on an angular interval of 2π.

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, as

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

All fields are constant along the y-direction. The profiles of TE polarized WGMs 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 polarized WGMs can be written as E(r) = (Er, 0, Eθ)(r), and H(r) = (0, Hy, 0)(r).

Being solutions of eigenvalue problems, the whispering gallery modes are determined up to some complex constant only. No units are shown for their electric or magnetic fields. The given values correspond to a level of one in the principal field component at the rim r = R of the cavity. Units of V/µm can be assumed for all electric fields, 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 ωc t) at θ = 0. The animations show the respective component at recurring points in time, including the decay effected by the imaginary part of ωc.

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 whispering gallery mode at t = 0. The global phase of the solution can be adjusted; the animations show the respective component at recurring points in time, taking into account the field decay.

Referring to the coordinates as introduced in the figure, the WGMs as calculated by the solver, with positive angular mode order m, represent fields that propagate clockwise in postive θ-direction. For each mode there exists an eigenfunction with the same complex eigenfrequency and negative angular mode order -m, that corresponds to an anticlockwise propagating resonant wave. The solver permits to inspect the respective fields, if the Rev-checkbox is selected. Note that the profiles of the reversed modes differ from the profiles of the mode with positive angular orders in the signs of the radial components. Superpositions of those two degenerate WGMs, with arbitrary complex coefficients, constitute further valid solutions for that particular resonance frequency. Options for inspecting the fields of those superimposed states are available for the propagation plots.

Select points on the planes of the plots to inspect precise local field levels. Clicks outside the actual axes close the coordinate displays. In case of a propagation plot, the "▯"-button toggles a colorbar. Select a field level on the colorbar to superimpose the field plot with a pseudo-contour at that level. Also here the contours are removed by a click in the colorbar area outside the axes.


[1] E.F. Franchimon, K.R. Hiremath, R. Stoffer, M. Hammer
Interaction of whispering gallery modes in integrated optical micro-ring or -disk circuits: Hybrid CMT model
Journal of the Optical Society of America B 30 (4), 1048-1057 (2013)  Preprint (ps.gz)  Preprint (pdf)  External online source