Optimization of the performance of flat optical components, also dubbed metasurfaces, is a crucial step towards their implementation in realistic optical systems. Yet, most of the design techniques, which rely on large parameter search to calculate the optical scattering response of elementary building blocks, do not account for near-field interactions that strongly influence the device performance. In this work, we exploit two advanced optimization techniques based on statistical learning and evolutionary strategies together with a fullwave high order Discontinuous Galerkin Time-Domain (DGTD) solver to optimize phase gradient metasurfaces. We first review the main features of these optimization techniques and then show that they can outperform most of the available designs proposed in the literature. Statistical learning is particularly interesting for optimizing complex problems containing several global minima/maxima. We then demonstrate optimal designs for GaN semiconductor phase gradient metasurfaces operating at visible wavelengths. Our numerical results reveal that rectangular and cylindrical nanopillar arrays can achieve more than respectively 88% and 85% of diffraction efficiency for TM polarization and both TM and TE polarization respectively, using only 150 fullwave simulations.
Comparison between the classical approach to phase gradient metasurface design and our inverse design approach based on the CMA-ES (Covariance Matrix Adaptation Evolution Strategy) algorithm. Top figures: geometry obtained using the classical approach in which each nanopillar is optimized manually by changing the diameter and finally placed on the substrate in order to obtain the desired phase shift needed to maximize the light deflection for the first order mode at λ = 600 nm. Bottom figures: results obtained using the CMA-ES optimization algorithm.
M. Elsawy, S. Lanteri, R. Duvigneau, G. Briere, M.S. Mohamed and P. Genevet
Global optimization of metasurface designs using statistical learning methods