Solution of large-scale plasmonic problems with the multilevel fast multipole algorithm.

A surface integral equation together with the multilevel fast multipole algorithm is successfully applied to fast and accurate resolution of plasmonic problems involving a large number of unknowns. The absorption, scattering, and extinction efficiencies of several plasmonic gold spheres of increasing size are efficiently obtained solving the electric and magnetic current combined-field integral equation. The numerical predictions are compared with reference analytic results to demonstrate the accuracy, suitability, and capabilities of this approach when dealing with large-scale plasmonic problems.