A new method is proposed to predict the topological properties of 1D periodic structures in wave physics, including quantum mechanics. From Bloch waves, a unique complex valued function is constructed, exhibiting poles and zeros. The sequence of poles and zeros of this function is a topological invariant that can be linked to the Berry–Zak phase. Since the characterization of the topological properties is done in the complex plane, it can easily be extended to the case of non‐Hermitian systems. The sequence of poles and zeros allows to predict topological phase transitions.

Beginning with Fermat’s principle, we provide a detailed derivation of the generalized laws of refraction and reflection for a geometry realizing a metasurface. We first solve the Euler–Lagrange equations for a light ray propagating across the metasurface. The ray-path equation is found analytically, and the results are supported by numerical calculations. We get generalized laws of refraction and reflection that have three main features: (i) They are relevant in gradient-index optics and in geometrical optics; (ii) A collection of rays emerges from the metasurface as a result of multiple reflections inside the metasurface; and (iii) The laws, although derived from Fermat’s principle, differ from previously published results.

Modelling of multi-wavelength metasurfaces relies on adjusting the phase of indi-vidual nanoresonators at several wavelengths.The traditional procedure neglects thenear-field coupling between the nanoresonators, which dramatically reduces the over-all diffraction efficiency, bandwidth, numerical aperture and device diameter.Anotheralternative design strategy is to combine a numerical optimization technique with full-wave simulations to mitigate this problem and optimize the entire metasurface at once.Here, we present a global multiobjective optimization technique that utilizes statisticallearning method to optimize RGB spherical metalenses at the visible wavelengths. Theoptimization procedure, coupled to a high-order full-wave solver, accounts for the nearfield coupling between the resonators. High numerical aperture RGB lenses(NA= 0.47and NA= 0.56) of 8μm and 10μm diameters are optimized with numerical average1focusing efficiencies of 55% and 45%, with an average focusing error smaller than 6%for the RGB colors. The fabricated and experimentally characterized devices present44.16% and 31.5% respective efficiencies. The reported performances represent thehighest focusing efficiencies for highNA >0.5 RGB metalenses obtained so far. Theintegration of multi-wavelength metasurfaces in portable and wearable electronic de-vices requires high performances to offer a variety of applications ranging from classicalimaging to virtual and augmented reality.

We clarify some analytical expressions existing in the literature, • Within the correct formulae we conclude that there is no need for an ad hoc improvement on the opposite to the title paper, • We highlight the symmetry properties of the function to be integrated in order to agree with the usual assumptions made to derive the Kramers-Kronig relations, • The analytical formula we provide may be used to increase the accuracy of the "Poor Man's Kramers-Kronig analysis" method and the "Kramers-Kronig constrained variational analysis" method.