Dipolar Modeling of Multipolar Metasurfaces

Abstract: Multipolar decomposition is a powerful tool for analyzing and designing metasurfaces, but its practical application is often limited by the mathematical complexity that arises when a large number of multipole moments must be taken in to account. To minimize this modeling complexity without sacrificing accuracy, we present an efficient method that exploits the coordinates origin dependence of spherical multipole moments. We show that the optimal origins for minimizing higher-order contributions, such as quadrupoles and octupoles, depend strictly on the spatial parity of the electromagnetic response. This is achieved by modeling a metasurface response using multipolar generalized sheet transition conditions (GSTCs). By separating the GSTCs into independent even and odd parity components, we can evaluate the electric and magnetic discontinuities at distinct physical positions. This parity-splitting framework allows us to systematically suppress unwanted higher-order terms and reconstruct the complete scattering parameters using only the dipole moments. We validate our analytical approach using two numerical examples: vertically asymmetric dielectric cones on a substrate, and a horizontally symmetry-broken metasurface supporting a double quasi-bound state in the continuum resonance. In both cases, the retrieved scattering parameters show excellent agreement with full-wave simulations. This method provides a simple, physically intuitive framework that simplifies the modeling of geometrically complex and non-local metasurfaces down to a purely dipolar level.