Simple strategy for the simulation of axially symmetric metasurfaces of large area

Metalenses are composed of nanostructures for focusing light and have been widely explored in many exciting applications. However, their expanding dimensions pose simulation challenges. We propose a method to simulate metalenses in a timely manner using vectorial wave and ray tracing models. We sample the metalens' radial phase gradient and locally approximate it by a linear phase response. Each sampling point is modeled as a binary blazed grating, employing the chosen nanostructure, to build a transfer function set. The metalens transmission or reflection is then obtained by applying the corresponding transfer function to the incoming field on the regions surrounding each sampling point. Fourier optics is used to calculate the scattered fields under arbitrary illumination for the vectorial wave method and a Monte Carlo algorithm is used in the ray tracing formalism. We validated our method against finite difference time domain simulations at 632 nm and we were able to simulate metalenses larger than 3000λ₀ in diameter on a personal computer.