Modeling large metasurfaces comprised of nonuniform plasmonic nanorods arrays

In this work, we present a numerically efficient technique, called the Characteristics Basis Function Method (CBFM), to model large non-uniform plasmonic arrays of nanorods. The use of the CBFM dramatically reduces the number of unknowns without sacrificing the computational accuracy and enables us to handle large, truncated and non-uniform arrays. The method is rigorous and includes the mutual coupling effects to obtain accurate results. Furthermore, the method is readily parallelizable and is capable of handling arrays comprised of nanoantennas with complex shapes, without any difficulty.