Modelling of multi-wavelength metasurfaces relies on adjusting the phase indi-vidual nanoresonators at several wavelengths.The traditional procedure neglects thenear-field coupling between nanoresonators, which dramatically reduces over-all diffraction efficiency, bandwidth, numerical aperture and device diameter.Anotheralternative design strategy is to combine a optimization technique with ful...