On the Numerical Simulation of Metasurfaces With Impedance Boundary Condition Integral Equations

Metasurfaces are thin metamaterial layers characterized by unusual dispersion properties of surface/guided wave and/or reflection properties of otherwise incident plane waves. At the scales intervening in their design, metasurfaces can be described through a surface impedance boundary condition. The impedance, possibly tensorial, is often “modulated,” i.e., it can vary from place to place on the surface (by design). We investigate on different integral equation formulations of the problem, with special attention to the stability properties of the resulting system matrix.