Simulation of Photonic Metasurfaces using a Surface Integral Equation Solver with Generalized Sheet Transition Conditions

Photonic metasurfaces are efficiently and accurately simulated using a surface integral equation (SIE) solver. This solver models the metasurface as an infinitesimally thin sheet on which generalized sheet transition conditions (GSTCs) are enforced. GSTCs connect the electromagnetic fields on both sides of the metasurface using equivalent susceptibility tensors and avoid modelling of geometrically intricate metasurface unit cells by the SIE solver. Susceptibility tensors are obtained from the transmission and the reflection coefficients associated with a unit cell. A numerical example, involving a metasurface that supports high-quality factor resonances and excited by a Gaussian beam, is presented to demonstrate the accuracy of the proposed approach.