Multilevel Fast Multipole Algorithm for Electromagnetic Scattering by Large Metasurfaces Using Static Mode Representation

Metasurfaces, consisting of large arrays of interacting subwavelength scatterers, pose significant challenges for general-purpose computational methods due to their large electrical size and multiscale nature. In this communication, leveraging the Poggio–Miller–Chang–Harrington–Wu–Tsai (PMCHWT) formulation, we combine the multilevel fast multipole algorithm (MLFMA) with a static mode representation (SMR) of the unknown equivalent surface current density. These static modes are entire-domain basis functions that depend solely on the object’s shape and are independent of material properties and frequency by Forestiere et al. (2023). By compressing the number of unknowns using SMR and exploiting the efficient O(N log N) scaling of the MLFMA matrix-vector products, our MLFMA-SMR method substantially reduces CPU time and memory requirements compared to traditional MLFMA implementations using Rao–Wilton–Glisson (RWG) basis functions. We assess the method’s accuracy and computational cost (time and memory) through several test cases, including the full-wave simulation of a 100 × 100λ canonical metalens. Overall, the MLFMA-SMR method offers substantial benefits for the analysis and optimization of large-scale metasurfaces and metalenses.