Fast and accurate analysis of optical metamaterials using surface integral equations and the parallel multilevel fast multipole algorithm

We present fast and accurate simulations of optical metamaterials using surface integral equations and the multilevel fast multipole algorithm (MLFMA). Problems are formulated with the electric and magnetic current combined-field integral equation and solved iteratively with MLFMA, which is parallelized using the hierarchical strategy on distributed-memory architectures. Realistic metamaterials involving dielectric, perfectly conducting, and plasmonic regions of finite extents are solved rigorously with the developed implementation without any periodicity assumptions.