Multilayer waveguides: efficient numerical analysis of general structures

Abstract
An efficient numerical method is presented for accurately determining the real and/or complex propagation constants of guided modes and leaky waves in general multilayer waveguides. The method is applicable to any lossless and/or lossy (dielectric, semiconductor, metallic) waveguide structure. The method is based on the argument principle theorem and is capable of extracting all of the zeros of any analytic function in the complex plane. It is applied in the present work to solving the multilayer waveguide dispersion equation derived from the well known thin-film transfer matrix theory. Excellent agreement is found with seven previously published results and with results from two limiting cases where the propagating constants can be obtained analytically.