Comparison of iterative solver performances on multiple surface integral equation formulations for plasmonic scatterers

In computational electromagnetics, five surface integral equation (SIE) formulations are typically used to predict the electromagnetic behavior of plasmonic structures. These SIE formulations are discretized into a matrix form by the method of moments (MoM) approach. The derived linear system of equations needs to be typically resolved by iterative algorithms, since the direct inversion of the MoM matrix is unfeasible for large computational problems. This paper compares the performance of four different iterative solvers (GMRES, TFQMR, CGS, and BICGSTAB) for five different SIE-MoM formulations (PMCHWT, JMCFIE, CTF, CNF, and MNMF) used for the electromagnetic analysis of plasmonic structures.