Robust Location of Optical Fiber Modes via the Argument Principle Method

Abstract

We implement a robust, globally convergent root search method for transcendental equations guaranteed to locate all complex roots within a specified search domain, based on Cauchy's residue theorem. Although several implementations of the argument principle already exist, ours has several advantages: it allows singularities within the search domain and branch points are not fatal to the method. Furthermore, our implementation is simple and is written in MATLAB, fulfilling the need for an easily integrated implementation which can be readily modified to accommodate the many variations of the argument principle method, each of which is suited to a different application. We apply the method to the step index fiber dispersion relation, which has become topical due to the recent proliferation of high index contrast fibers. We also find modes with permittivity as the eigenvalue, catering to recent numerical methods that expand the radiation of sources by eigenmodes.

24 pages, 4 figures

References in corpus (3)

• Strong coupling of single emitters to surface plasmons
• Self-consistent multi-mode lasing theory for complex or random lasing media
• Modal Analysis and Coupling in Metal-Insulator-Metal Waveguides

Cited by in corpus (7)

• Generalizing normal mode expansion of electromagnetic Green's tensor to lossy resonators in open systems
• Global Complex Roots and Poles Finding Algorithm Based on Phase Analysis for Propagation and Radiation Problems
• Revisiting the boundary conditions for second-harmonic generation at metal-dielectric interfaces
• Overcoming the bottleneck for quantum computations of complex nanophotonic structures: Purcell and FRET calculations using a rigorous mode hybridization method
• Surface second-harmonic from metallic nanoparticle configurations - a transformation optics approach
• An efficient solver for the generalized normal modes of non-uniform open optical resonators
• Resolving the Gibbs phenomenon via a discontinuous basis in a mode solver for open optical systems