Do dielectric bilayer metasurfaces behave as a stack of decoupled single-layer metasurfaces?

Flat optics or metasurfaces have opened new frontiers in wavefront shaping and its applications. Polarization optics is one prominent area which has greatly benefited from the shape-birefringence of metasurfaces. However, flat optics comprising a single layer of meta-atoms can only perform a subset of polarization transformations, constrained by a symmetric Jones matrix. This limitation can be tackled using metasurfaces composed of bilayer meta-atoms but exhausting all possible combinations of geometries to build a bilayer metasurface library is a very daunting task. Consequently, bilayer metasurfaces have been widely treated as a cascade (product) of two decoupled single-layer metasurfaces. Here, we test the validity of this assumption for dielectric metasurfaces by considering a metasurface made of titanium dioxide on fused silica substrate at a design wavelength of 532 nm. We explore regions in the design space where the coupling between the top and bottom layers can be neglected, i.e., producing a far-field response which approximates that of two decoupled single-layer metasurfaces. We complement this picture with the near-field analysis to explore the underlying physics in regions where both layers are strongly coupled. We also show the generality of our analysis by applying it to silicon metasurfaces at telecom wavelengths. Our unified approach allows the designer to efficiently build a multi-layer dielectric metasurface, either in transmission or reflection, by only running one full-wave simulation for a single-layer metasurface.