Abstract
In this paper, we propose a general boundary integral equation (BIE) approach for solving both exterior (e.g., scattering) and interior (e.g., guided wave propagation) problems involving general plasmonic waveguiding structures with arbitrary cross-sectional geometries and with a continuous translational symmetry in the direction of wave propagation. Contrary to the field-based approach which deals with the electric and magnetic fields, we employ a potential-based formalism instead, involving vector and scalar potentials, and match them at the media interfaces. The proposed approach can handle not only conventional plasmonic waveguides on the order of a few hundred nanometers but also those whose critical sizes are a few nanometers and in which the nonlocal hydrodynamic effects need to be accounted for. The BIEs describing the interaction of light with the plasmonic waveguide are solved by the method of moments (MoM) algorithm. Two illustrative examples, the first one of which deals with the scattering problem, while the second computes the dispersion diagram of a plasmonic cylinder, are considered for both local and nonlocal cases. An excellent agreement is observed between the numerical and theoretical results, with the latter being derived from the generalized Mie theory.