International audienceOptimization 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...