Multilayer waveguides: efficient numerical analysis of general structures

An efficient numerical method for accurately determining the real and/or complex propagation constants of guided modes and leaky waves in general multilayer waveguides is presented. 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 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.< >