Fast Hybrid Method Applying To Metasurfaces Without And With Substrate

A computer-implemented method for full-wave simulations of metasurfaces, comprising:

• receiving, by a computer processor, a metasurface comprised of a plurality of scattering structures residing on a surface of a substrate, where the plurality of scattering structures have same shape;
• determining a T-matrix for the surface of the substrate in vector plane waves including evanescent waves
• determining, by a compute processor, a baseline T-matrix for a given scattering structure of the metasurface;
• generating, by the computer processor, a model for the metasurface by assigning a T-matrix to each scattering structure comprising the metasurface, where each assigned T-matrix is a variant of the baseline T-matrix;
• for each scattering structure in the model, formulating a Foldy-Lax equation for a particular scattering structure, thereby forming a set of Foldy-Lax equations, where each Foldy-Lax equation accounts for scattering from each of remaining scattering structures in the model and accounts for scattering of the surface of the substrate;
• and
• determining a scattered wave for each of the plurality of scattering structures by applying a Fast Fourier transform to the set of Foldy-Lax equations.