Plasmonic transmission lines mode solver based on the Method of Moments

In this paper, we apply the integral equation formulation to calculate the mode characteristics of plasmonic transmission lines. These lines can be constructed from any number of arbitrary shaped metallic strips placed inside a stack of planar dielectric or metallic layers. The number of the metallic strips determines the number of the propagating modes. Each mode is characterized by its attenuation constant, phase constant, and field distribution. The adopted formulation is based on solving Maxwell´s equations using the Method of Moments (MoM) technique. This formulation is applied for single and coupled strips of various cross-sections over the frequency range 150-450 THz. Results are compared to Lumerical and CST software for the aim of verification. Good agreement between our solver and these tools has been observed.