Design and analysis of chiral and achiral metasurfaces with the finite element method.

The rise of metasurfaces to manipulate the polarization states of light motivates the development of versatile numerical methods able to model and analyze their polarimetric properties. Here we make use of a scattered-field formulation well suited to the Finite Element Method (FEM) to compute the Stokes-Mueller matrix of metasurfaces. The major advantage of the FEM lies in its versatility and its ability to compute the optical properties of structures with arbitrary and realistic shapes, and rounded edges and corners.