High-Order Accurate Integral Equation Based Mode Solver for Layered Nanophotonic Waveguides

We present a novel high-order boundary integral equation based mode solver for layered nanophotonic waveguides featuring corners and junctions. The proposed solver leverages an indirect single-trace formulation for composite dielectrics which leads to a Muller-like integral equation posed on the cross-section of the material interfaces, and a high-order Nystrom discretization for domains with corners. A variety of numerical examples demonstrate the superior accuracy of the proposed methodology as compared to standard commercial finite element and finite difference based mode solvers.