<figure><img src="image_1.png" /></figure>

HAL
open science

# Advancing wavefront shaping with resonant nonlocal metasurfaces: beyond the limitations of lookup tables

Enzo Isnard, Sébastien Héron, Stéphane Lanteri, Mahmoud Elsawy

► To cite this version:

Enzo Isnard, Sébastien Héron, Stéphane Lanteri, Mahmoud Elsawy. Advancing wavefront shaping with resonant nonlocal metasurfaces: beyond the limitations of lookup tables. Scientific Reports, 2024, 14, pp.1555. 10.1038/s41598-024-51898-8 . inserm-04210828v2

HAL Id: inserm-04210828

https://inserm.hal.science/inserm-04210828v2

Submitted on 25 Jan 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

Public Domain
<figure><img src="image_1.png" /></figure>

HAL
open science

# Advancing wavefront shaping with resonant nonlocal metasurfaces: beyond the limitations of lookup tables

Enzo Isnard, Sébastien Héron, Stéphane Lanteri, Mahmoud Elsawy

► To cite this version:

Enzo Isnard, Sébastien Héron, Stéphane Lanteri, Mahmoud Elsawy. Advancing wavefront shaping with resonant nonlocal metasurfaces: beyond the limitations of lookup tables. Scientific Reports, 2024, 14, pp.1555. 10.1038/s41598-024-51898-8 . inserm-04210828v2

HAL Id: inserm-04210828

https://inserm.hal.science/inserm-04210828v2

Submitted on 25 Jan 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

Public Domain
<figure><img src="image_2.png" /></figure>

OPEN

# Advancing wavefront shaping with resonant nonlocal metasurfaces: beyond the limitations of lookup tables

Enzo Isnard<fnref n="1" />,², Sébastien Héron², Stéphane Lanteri<fnref n="1" /> & Mahmoud Elsawy ID <fnref n="1" />✉

Resonant metasurfaces are of paramount importance in addressing the growing demand for reduced thickness and complexity, while ensuring high optical efficiency. This becomes particularly crucial in overcoming fabrication challenges associated with high aspect ratio structures, thereby enabling seamless integration of metasurfaces with electronic components at an advanced level. However, traditional design approaches relying on lookup tables and local field approximations often fail to achieve optimal performance, especially for nonlocal resonant metasurfaces. In this study, we investigate the use of statistical learning optimization techniques for nonlocal resonant metasurfaces, with a specific emphasis on the role of near-field coupling in wavefront shaping beyond single unit cell simulations. Our study achieves significant advancements in the design theoretical conception of resonant metasurfaces. For transmission-based metasurfaces, a beam steering design outperforms the classical design by achieving an impressive efficiency of 80% compared to the previous 23%. Additionally, our optimized extended depth-of-focus (EDOF) metalens yields a remarkable five-fold increase in focal depth, a four-fold enhancement in focusing power compared to conventional designs and an optical resolution superior to 600 cycle/mm across the focus region. Moreover, our study demonstrates remarkable performance with a wavelength-selected beam steering metagrating in reflection, achieving exceptional efficiency surpassing 85%. This far outperforms classical gradient phase distribution approaches, emphasizing the immense potential for groundbreaking applications in the field of resonant metasurfaces.

Metasurfaces are flat optical components that have the ability to alter electromagnetic wave properties, including phase, amplitude, and polarization<fnref n="1" />⁻⁴. They are constructed using subwavelength-sized elements known as meta-atoms that interact with incoming waves to achieve a specific outcome. In this context, several applications have already been proposed such as ultra-compact metasurface-based deflectors, beam splitters, lenses, waveplates, holograms⁵⁻¹⁰.

In metasurface design, local and nonlocal effects refer to the way the properties of individual resonators affect the behavior of the whole structure and especially of their neighbors. A local metasurface is a type of metasurface in which the individual meta-atom only interacts with external electromagnetic waves that are incident directly on them. In other words, each resonator acts independently, and the electromagnetic response of the metasurface can be determined by a simple lookup table that relates the properties of the resonators to the desired wavefront. Generally speaking this approach could be applied for metasurfaces based on the modification of the effective refractive index of low-quality factor truncated nano-waveguides³,⁶,¹¹,¹². It is also applicable to metasurfaces which rely on the Pancharatnam Berry (PB) phase induced by birefringent meta-atoms.¹³⁻¹⁶.

On the other hand, nonlocal metasurface is a type of metasurface in which the meta-atoms cannot be considered isolated and the response of each meta-atom depends on the state of its neighbors¹⁷,¹⁸. In other words, in a nonlocal metasurface, the meta-atoms are not only interacting with the external waves that are incident directly on them, but also with the waves that are incident on their neighboring meta-atoms. This metasurface interacts

<footnote n="1">¹Université Côte d'Azur, Inria, CNRS, LJAD, 06902 Sophia Antipolis Cedex, France. ²THALES Research and Technology, 1 Avenue Augustin Fresnel, 91120 Palaiseau, France. ✉email: mahmoud.elsawy@inria.fr</footnote>
<figure><img src="image_2.png" /></figure>

OPEN

# Advancing wavefront shaping with resonant nonlocal metasurfaces: beyond the limitations of lookup tables

Enzo Isnard<fnref n="1" />,², Sébastien Héron², Stéphane Lanteri<fnref n="1" /> & Mahmoud Elsawy ID <fnref n="1" />✉

Resonant metasurfaces are of paramount importance in addressing the growing demand for reduced thickness and complexity, while ensuring high optical efficiency. This becomes particularly crucial in overcoming fabrication challenges associated with high aspect ratio structures, thereby enabling seamless integration of metasurfaces with electronic components at an advanced level. However, traditional design approaches relying on lookup tables and local field approximations often fail to achieve optimal performance, especially for nonlocal resonant metasurfaces. In this study, we investigate the use of statistical learning optimization techniques for nonlocal resonant metasurfaces, with a specific emphasis on the role of near-field coupling in wavefront shaping beyond single unit cell simulations. Our study achieves significant advancements in the design theoretical conception of resonant metasurfaces. For transmission-based metasurfaces, a beam steering design outperforms the classical design by achieving an impressive efficiency of 80% compared to the previous 23%. Additionally, our optimized extended depth-of-focus (EDOF) metalens yields a remarkable five-fold increase in focal depth, a four-fold enhancement in focusing power compared to conventional designs and an optical resolution superior to 600 cycle/mm across the focus region. Moreover, our study demonstrates remarkable performance with a wavelength-selected beam steering metagrating in reflection, achieving exceptional efficiency surpassing 85%. This far outperforms classical gradient phase distribution approaches, emphasizing the immense potential for groundbreaking applications in the field of resonant metasurfaces.

Metasurfaces are flat optical components that have the ability to alter electromagnetic wave properties, including phase, amplitude, and polarization<fnref n="1" />⁻⁴. They are constructed using subwavelength-sized elements known as meta-atoms that interact with incoming waves to achieve a specific outcome. In this context, several applications have already been proposed such as ultra-compact metasurface-based deflectors, beam splitters, lenses, waveplates, holograms⁵⁻¹⁰.

In metasurface design, local and nonlocal effects refer to the way the properties of individual resonators affect the behavior of the whole structure and especially of their neighbors. A local metasurface is a type of metasurface in which the individual meta-atom only interacts with external electromagnetic waves that are incident directly on them. In other words, each resonator acts independently, and the electromagnetic response of the metasurface can be determined by a simple lookup table that relates the properties of the resonators to the desired wavefront. Generally speaking this approach could be applied for metasurfaces based on the modification of the effective refractive index of low-quality factor truncated nano-waveguides³,⁶,¹¹,¹². It is also applicable to metasurfaces which rely on the Pancharatnam Berry (PB) phase induced by birefringent meta-atoms.¹³⁻¹⁶.

On the other hand, nonlocal metasurface is a type of metasurface in which the meta-atoms cannot be considered isolated and the response of each meta-atom depends on the state of its neighbors¹⁷,¹⁸. In other words, in a nonlocal metasurface, the meta-atoms are not only interacting with the external waves that are incident directly on them, but also with the waves that are incident on their neighboring meta-atoms. This metasurface interacts

<footnote n="1">¹Université Côte d'Azur, Inria, CNRS, LJAD, 06902 Sophia Antipolis Cedex, France. ²THALES Research and Technology, 1 Avenue Augustin Fresnel, 91120 Palaiseau, France. ✉email: mahmoud.elsawy@inria.fr</footnote>
with light in a way that is not local, meaning that the effect of meta-atoms on the light is not limited to their immediate vicinity. Instead, the behavior of the entire metasurface is determined by the collective behavior of all the meta-atoms, making it difficult to predict the response using a lookup table from the unit cell simulations.

The nonlocal interaction is inherently present for metasurface considering phase addressing mechanism relies on resonant scattering of light, and the associated light scattering modulation properties¹⁹⁻²². Recent works based on lookup table considerations revealed an intriguing connection between the full phase modulation and the topological effects governing the resonant light interaction of nanoparticles²³,²⁴.

Furthermore, the concept of nonlocal metasurfaces has been extended to structures whose governing principles are deeply rooted in the symmetry features of quasi-bound states in the continuum (q-BIC)¹⁷,²⁵,²⁶ for selective wavefront control. Yet, the overall efficiency (ratio of the power in the desired direction compared to the incident power) is modest since most of these studies predict the overall performance response using a simple lookup table from the single-unit cell simulation which neglect the cross talk between adjacent resonant elements.

The importance of resonant nonlocal metasurfaces becomes evident in their ability to facilitate a remarkable reduction in metasurface thickness and meta-atom complexity. This reduction in structural requirements is crucial for several reasons. First, it allows for the integration of metasurfaces within existing electronic components, paving the way for advanced functionalities and enhanced performance in integrated photonic systems. Second, the decreased thickness simplifies the fabrication process, enabling large-scale manufacturing with higher precision and reduced costs.

In this study, we demonstrate that the local approach is not suitable to predict the performance of highly resonant metasurfaces. Moreover, we propose an optimization method for nonlocal metasurfaces based on a global statistical learning optimization method referred as EGO (Efficient Global Optimization) which has been successfully applied to local metasurfaces⁵,²⁷,²⁸.

EGO is based on building a surrogate model of the objective function using a Gaussian process regression (GPR) approach. This model is used to approximate the true objective function and to guide the search for the optimal solution. EGO iteratively updates the surrogate model by adding new observations to it, and then uses an acquisition function to select the next point to evaluate. The acquisition function balances exploration and exploitation to identify the next point that is most likely to improve the objective function's value. EGO has several advantages over other optimization algorithms. First, the surrogate model that is built by GPR is capable of capturing complex and nonlinear relationships between input and output variables. Second, EGO balances exploration and exploitation, enabling it to efficiently search the design space for the global optimum while avoiding getting trapped in local optima. For more details about the optimization methodology, we refer to our recent works⁵,²⁷⁻²⁹.

One of the key strengths of EGO is its ability to handle expensive-to-evaluate black-box functions. In such cases, evaluating the objective function is computationally intensive, and traditional optimization methods may be too slow to converge to a satisfactory solution. EGO mitigates this by using a surrogate model to approximate the objective function, reducing the number of function evaluations required to find the global optimum. By using a Gaussian process regression surrogate model and a carefully designed acquisition function, EGO efficiently explores the design space to identify the global optimum while minimizing the number of function evaluations required. Its ability to balance exploration and exploitation and handle expensive black-box functions make it a popular choice for a wide range of optimization problems. Another key advantage of EGO is that the creation of a surrogate model during the optimization process makes it possible to perform uncertain quantification analysis and to adapt the algorithm to take into account the fabrication imperfections as illustrated in Ref. ²⁸.

We demonstrate the effectiveness of our method by applying it to the design of a nonlocal metasurface for wavefront control with various mechanisms ranging from classical Huygens's metasurface¹⁹ to wavelength selective wavefront metagratings. The results show that our method is able to achieve superior performance compared to traditional synthesizing approach based on locality assumption, with a significant reduction in the number of function evaluations compared to conventional global optimization methods¹⁸. Overall, our study provides a promising approach for the optimization of nonlocal metasurfaces, and highlights the potential of EGO method for solving complex electromagnetic design problems.

## Huygens's metasurface

Huygen's metasurface is a promising approach to reduce the aspect ratio of constitutive elements of metasurfaces, which is crucial to reduce their fabrication complexity. It Huygen's metasurface is a type of resonant metasurface that is designed to manipulate the scattering of light in a highly directional manner. It consists of an array of subwavelength resonators that are carefully engineered to exhibit both magnetic and electric dipole resonances at the same frequency. By controlling the relative phase and amplitude of the magnetic and electric dipoles in each resonator, the metasurface can achieve the required full phase modulation associated with near-unity transmission amplitude¹⁹,³⁰,³¹.

The Huygens's metasurface configuration has diverse potential applications, including optical communication, imaging, and sensing. However, the design of such metasurfaces is challenging due to the intricate interplay between the magnetic and electric dipole resonances, which is determined by the position of the poles and zeros as elaborated in Ref. ²³. Moreover, the nanoresonators have extended evanescent electromagnetic fields that couple their resonances beyond adjacent neighbors, making the conventional design of phase-gradient metasurfaces complicated. In addition, various factors, such as the geometrical parameters of the resonators, the polarization and wavelength of the incident light, and the surrounding environment, can affect the complex interplay between the resonances. Therefore, designing a high-performance Huygen's metasurface necessitates meticulous consideration of these factors and advanced optimization techniques to achieve the desired properties.
with light in a way that is not local, meaning that the effect of meta-atoms on the light is not limited to their immediate vicinity. Instead, the behavior of the entire metasurface is determined by the collective behavior of all the meta-atoms, making it difficult to predict the response using a lookup table from the unit cell simulations.

The nonlocal interaction is inherently present for metasurface considering phase addressing mechanism relies on resonant scattering of light, and the associated light scattering modulation properties¹⁹⁻²². Recent works based on lookup table considerations revealed an intriguing connection between the full phase modulation and the topological effects governing the resonant light interaction of nanoparticles²³,²⁴.

Furthermore, the concept of nonlocal metasurfaces has been extended to structures whose governing principles are deeply rooted in the symmetry features of quasi-bound states in the continuum (q-BIC)¹⁷,²⁵,²⁶ for selective wavefront control. Yet, the overall efficiency (ratio of the power in the desired direction compared to the incident power) is modest since most of these studies predict the overall performance response using a simple lookup table from the single-unit cell simulation which neglect the cross talk between adjacent resonant elements.

The importance of resonant nonlocal metasurfaces becomes evident in their ability to facilitate a remarkable reduction in metasurface thickness and meta-atom complexity. This reduction in structural requirements is crucial for several reasons. First, it allows for the integration of metasurfaces within existing electronic components, paving the way for advanced functionalities and enhanced performance in integrated photonic systems. Second, the decreased thickness simplifies the fabrication process, enabling large-scale manufacturing with higher precision and reduced costs.

In this study, we demonstrate that the local approach is not suitable to predict the performance of highly resonant metasurfaces. Moreover, we propose an optimization method for nonlocal metasurfaces based on a global statistical learning optimization method referred as EGO (Efficient Global Optimization) which has been successfully applied to local metasurfaces⁵,²⁷,²⁸.

EGO is based on building a surrogate model of the objective function using a Gaussian process regression (GPR) approach. This model is used to approximate the true objective function and to guide the search for the optimal solution. EGO iteratively updates the surrogate model by adding new observations to it, and then uses an acquisition function to select the next point to evaluate. The acquisition function balances exploration and exploitation to identify the next point that is most likely to improve the objective function's value. EGO has several advantages over other optimization algorithms. First, the surrogate model that is built by GPR is capable of capturing complex and nonlinear relationships between input and output variables. Second, EGO balances exploration and exploitation, enabling it to efficiently search the design space for the global optimum while avoiding getting trapped in local optima. For more details about the optimization methodology, we refer to our recent works⁵,²⁷⁻²⁹.

One of the key strengths of EGO is its ability to handle expensive-to-evaluate black-box functions. In such cases, evaluating the objective function is computationally intensive, and traditional optimization methods may be too slow to converge to a satisfactory solution. EGO mitigates this by using a surrogate model to approximate the objective function, reducing the number of function evaluations required to find the global optimum. By using a Gaussian process regression surrogate model and a carefully designed acquisition function, EGO efficiently explores the design space to identify the global optimum while minimizing the number of function evaluations required. Its ability to balance exploration and exploitation and handle expensive black-box functions make it a popular choice for a wide range of optimization problems. Another key advantage of EGO is that the creation of a surrogate model during the optimization process makes it possible to perform uncertain quantification analysis and to adapt the algorithm to take into account the fabrication imperfections as illustrated in Ref. ²⁸.

We demonstrate the effectiveness of our method by applying it to the design of a nonlocal metasurface for wavefront control with various mechanisms ranging from classical Huygens's metasurface¹⁹ to wavelength selective wavefront metagratings. The results show that our method is able to achieve superior performance compared to traditional synthesizing approach based on locality assumption, with a significant reduction in the number of function evaluations compared to conventional global optimization methods¹⁸. Overall, our study provides a promising approach for the optimization of nonlocal metasurfaces, and highlights the potential of EGO method for solving complex electromagnetic design problems.

## Huygens's metasurface

Huygen's metasurface is a promising approach to reduce the aspect ratio of constitutive elements of metasurfaces, which is crucial to reduce their fabrication complexity. It Huygen's metasurface is a type of resonant metasurface that is designed to manipulate the scattering of light in a highly directional manner. It consists of an array of subwavelength resonators that are carefully engineered to exhibit both magnetic and electric dipole resonances at the same frequency. By controlling the relative phase and amplitude of the magnetic and electric dipoles in each resonator, the metasurface can achieve the required full phase modulation associated with near-unity transmission amplitude¹⁹,³⁰,³¹.

The Huygens's metasurface configuration has diverse potential applications, including optical communication, imaging, and sensing. However, the design of such metasurfaces is challenging due to the intricate interplay between the magnetic and electric dipole resonances, which is determined by the position of the poles and zeros as elaborated in Ref. ²³. Moreover, the nanoresonators have extended evanescent electromagnetic fields that couple their resonances beyond adjacent neighbors, making the conventional design of phase-gradient metasurfaces complicated. In addition, various factors, such as the geometrical parameters of the resonators, the polarization and wavelength of the incident light, and the surrounding environment, can affect the complex interplay between the resonances. Therefore, designing a high-performance Huygen's metasurface necessitates meticulous consideration of these factors and advanced optimization techniques to achieve the desired properties.
## Huygen's metasurface for beam steering

To demonstrate the effectiveness of our statistical learning optimization method in designing efficient nonlocal metasurfaces, we first consider the resonant structures presented in³² to design beam steering device. The structure is composed of silicon nanocylinders with a refractive index of *n* = 3.56 embedded in a host silica with a refractive index of *n* = 1.45, as shown in the inset in Fig. 1a, with a fixed height of *h* = 170 nm. We initially follow the local approach design methodology by constructing a look-up table starting from a single unit-cell with periodic boundary conditions. The subwavelength period is fixed at *P* = 620 nm. In Fig. 1a, we vary the radius of the nano-resonator to achieve the desired response, i.e., full phase modulation (blue curve) associated with high-transmission response (green color) for a fixed wavelength of *λ* = 1064 nm. To complete the first step, we designed a beam steering design by considering 8 unit cells (to sample 0-2π phase) with a phase gradient distribution. We focus on the Transverse Electric (TE) polarized light case with the electric field oscillating in the y-direction³³ based on this classical synthesizing phase approach, the maximum light deflection efficiency in the first-order mode does not exceed 35% at *λ* = 1064 nm (see blue curve in Fig. 1a). This is due to the strong near-field coupling related to the nonlocal coupling between the resonators, which is not taken into account in the classical design approach (see Fig. 1c).

Unlike the classical design, in order to optimize the near-field coupling, the optimization scenario requires optimizing the radius of each resonator along with the distances between them. For the 8-pixels configuration, we optimize 15 parameters and the objective is to maximize light in the first order mode. This is a challenging optimization scenario, therefore, we rely on our advanced numerical methodology by considering high-level statistical learning optimization technique that requires less iterations compared to classical global methods²⁷ together with a high-fidelity 3D electromagnetic solver from the DIOGENES suite³⁴ to capture effectively the strong-near field coupling. The optimization results are depicted in Fig. 1d. The blue points refer to the Design Of Experiment (DOE) iterations, while the black points indicate the optimization phase. The green curve represents the best values found during the optimization process. For more details about the optimization methodology, we refer to our recent works⁵,²⁷⁻²⁹. The best design can reach a deflection efficiency of 80% as it is indicated in Fig. 1e. The huge improvement in the performance compared to the classical design is due to the near-field coupling, which is now taken into account in the optimization scenario as illustrated in the field profile depicted in Fig. 1f.

Moreover, the optimized design reveals the presence of high-order modes in contrast to the classical design based on local periodic approximation, as depicted in the last element of the period illustrated in Fig. 1f. This occurrence signifies the influence of strong non-local effects, an important factor not considered in the classical

<figure><img src="image_3.png" /></figure>

**Figure 1.** Optimization results for resonant Huygens's metasurface. (a) Refers to phase and transmission response for single unit-cell simulations as a function of the radius of the silicon pillar at $\lambda$ = 1064 nm. The inset in (a) indicates the schematic representation of the structure under consideration. We consider a silicon pillar encapsulated in SiO₂ material. The height of the pillars and the period are fixed as $h$ = 170 nm and $P$ = 620 nm. The radius of the pillars can vary from 130 nm to 240 nm. (b,c) beam steering configuration based on classical synthesizing approach where the radii of the cylinders are varying in a gradient manner based on the look-up table given in (a), for a fixed period $P$ = 620 nm. The field map represents $\Re(E_y)$ at $\lambda$ = 1064 nm. (d) optimization results for 8 pillars by varying the radii of each pillar and the distances in between them (15 parameters). The blue points refer to the design of experiment phase (DOE), the black points indicate the EGO iterations. The green curve is for the best values during the optimization process. (e) Power deflected in the first order mode for the optimal design obtained from the EGO, the corresponding field profile of $\Re(E_y)$ is given in (f).
## Huygen's metasurface for beam steering

To demonstrate the effectiveness of our statistical learning optimization method in designing efficient nonlocal metasurfaces, we first consider the resonant structures presented in³² to design beam steering device. The structure is composed of silicon nanocylinders with a refractive index of *n* = 3.56 embedded in a host silica with a refractive index of *n* = 1.45, as shown in the inset in Fig. 1a, with a fixed height of *h* = 170 nm. We initially follow the local approach design methodology by constructing a look-up table starting from a single unit-cell with periodic boundary conditions. The subwavelength period is fixed at *P* = 620 nm. In Fig. 1a, we vary the radius of the nano-resonator to achieve the desired response, i.e., full phase modulation (blue curve) associated with high-transmission response (green color) for a fixed wavelength of *λ* = 1064 nm. To complete the first step, we designed a beam steering design by considering 8 unit cells (to sample 0-2π phase) with a phase gradient distribution. We focus on the Transverse Electric (TE) polarized light case with the electric field oscillating in the y-direction³³ based on this classical synthesizing phase approach, the maximum light deflection efficiency in the first-order mode does not exceed 35% at *λ* = 1064 nm (see blue curve in Fig. 1a). This is due to the strong near-field coupling related to the nonlocal coupling between the resonators, which is not taken into account in the classical design approach (see Fig. 1c).

Unlike the classical design, in order to optimize the near-field coupling, the optimization scenario requires optimizing the radius of each resonator along with the distances between them. For the 8-pixels configuration, we optimize 15 parameters and the objective is to maximize light in the first order mode. This is a challenging optimization scenario, therefore, we rely on our advanced numerical methodology by considering high-level statistical learning optimization technique that requires less iterations compared to classical global methods²⁷ together with a high-fidelity 3D electromagnetic solver from the DIOGENES suite³⁴ to capture effectively the strong-near field coupling. The optimization results are depicted in Fig. 1d. The blue points refer to the Design Of Experiment (DOE) iterations, while the black points indicate the optimization phase. The green curve represents the best values found during the optimization process. For more details about the optimization methodology, we refer to our recent works⁵,²⁷⁻²⁹. The best design can reach a deflection efficiency of 80% as it is indicated in Fig. 1e. The huge improvement in the performance compared to the classical design is due to the near-field coupling, which is now taken into account in the optimization scenario as illustrated in the field profile depicted in Fig. 1f.

Moreover, the optimized design reveals the presence of high-order modes in contrast to the classical design based on local periodic approximation, as depicted in the last element of the period illustrated in Fig. 1f. This occurrence signifies the influence of strong non-local effects, an important factor not considered in the classical

<figure><img src="image_3.png" /></figure>

**Figure 1.** Optimization results for resonant Huygens's metasurface. (a) Refers to phase and transmission response for single unit-cell simulations as a function of the radius of the silicon pillar at $\lambda$ = 1064 nm. The inset in (a) indicates the schematic representation of the structure under consideration. We consider a silicon pillar encapsulated in SiO₂ material. The height of the pillars and the period are fixed as $h$ = 170 nm and $P$ = 620 nm. The radius of the pillars can vary from 130 nm to 240 nm. (b,c) beam steering configuration based on classical synthesizing approach where the radii of the cylinders are varying in a gradient manner based on the look-up table given in (a), for a fixed period $P$ = 620 nm. The field map represents $\Re(E_y)$ at $\lambda$ = 1064 nm. (d) optimization results for 8 pillars by varying the radii of each pillar and the distances in between them (15 parameters). The blue points refer to the design of experiment phase (DOE), the black points indicate the EGO iterations. The green curve is for the best values during the optimization process. (e) Power deflected in the first order mode for the optimal design obtained from the EGO, the corresponding field profile of $\Re(E_y)$ is given in (f).
design. In order to clarify this concept, in Fig. 2, we present the response of the last element surrounded with periodic boundary condition compared to its response inside the optimized grating. Our results clearly illustrate that the higher order mode is absent in the periodic case. However, due to the coupling between adjacent elements, the higher order mode can be excited, significantly contributing to the enhancement of overall efficiency.

## Extended depth of focus metalens (EDOF)

Metalenses have emerged as a highly promising technology in the field of optics, capable of effectively assuming the role of extended depth of focus (EDOF) lenses. Diverging from traditional lenses, EDOF lenses exhibit a consistent point spread function (PSF) throughout an extended range along the optical axis. By leveraging EDOF lenses, not only can novel functionalities be realized, such as the ability to bring objects at different distances from the lens into sharp focus, but they also alleviate the stringent constraints associated with precisely aligning lenses on top of a sensor. Extending the DOF also increases the spectral operation range because it reduces the impact of the focal shift due to chromatic aberrations on the quality of the images obtained³⁵.

Conventional techniques, including cubic, shifted axicon, log-asphere, and SQUBIC lenses, have been proposed for synthesizing phase masks to achieve EDOF properties³⁶. However, these traditional synthesis approaches encounter difficulties in striking a balance between maximizing focal efficiency and extending the focal depth. The limitations of the lookup table approach become apparent in addressing the complex requirements of achieving the desired EDOF properties. Notably, inverse design methods based on traditional optimization mechanisms have been utilized to maximize EDOF metalens performance, particularly in relying on classical optimization algorithm³⁷,³⁸ with numerical aperture not exceeding 0.44 accompanied with huge number of simulations.

The proposed design for this application involves a cylindrical metalens composed of unit cells described in Fig. 1a. The lens comprises 18 unit cells aligned along the x-axis with periodic boundary conditions along the y-axis to mimic 1D metalens. We start from a design with a NA of 0.6 and we extend the focal depth in a rectangle from the $f$ - DOF to $f + 5$DOF with a width of 2 FWHM where $f$ is the focal distance and DOF is given by the formula in³⁹.

In order to extend the depth of focus and achieve a uniform optical field distribution within the target depth, we define the objective function as the geometric mean of the mean power and the minimum power obtained from 20 equally spaced surfaces above and below the focal point. Using a minimum in the objective is primordial in order to obtain a design with a homogeneous intensity in the target region. This objective function allows us to simultaneously increase the focal depth along the z-axis, decrease the FWHM along the x-axis, and enhance the focusing power efficiency.

The employed geometric parametrization in this study aligns with the methodology presented in Section 2.1, which encompasses the diameter and distances between the resonators. However, to enhance the degrees of freedom, we also explored the variation of the resonator height, ranging from 150 nm to 250 nm, resulting in the optimization of a total of 19 parameters. The optimization process, as depicted in Fig. 3a, initiates with a Design of Experiments (DOE) consisting of 90 designs, indicated by the blue points. The process concludes after 450 simulations, represented by the black points in Fig. 3a. The optimal design achieved a performance of 30% according to our objective function. Furthermore, the electric field intensity ($|\mathbf{E}|^2$) demonstrates an uniform distribution across the entire focal zone, as shown in the first sub-figure of Fig. 3b.

<figure><img src="image_4.png" /><figcaption>**Figure 2.** Fields maps of $\Re(E_y)$ around the last meta-element in the supercell of the optimized grating shown in Fig. 1f when surrounded with periodic boundary conditions or inside the grating. The black rectangle delimits the meta-atoms boundaries. When there is a translational symmetry only the electric and magnetic dipoles are exited while when this symmetry is broken due to the strong coupling between elements, higher order modes can be excited.</figcaption></figure>
design. In order to clarify this concept, in Fig. 2, we present the response of the last element surrounded with periodic boundary condition compared to its response inside the optimized grating. Our results clearly illustrate that the higher order mode is absent in the periodic case. However, due to the coupling between adjacent elements, the higher order mode can be excited, significantly contributing to the enhancement of overall efficiency.

## Extended depth of focus metalens (EDOF)

Metalenses have emerged as a highly promising technology in the field of optics, capable of effectively assuming the role of extended depth of focus (EDOF) lenses. Diverging from traditional lenses, EDOF lenses exhibit a consistent point spread function (PSF) throughout an extended range along the optical axis. By leveraging EDOF lenses, not only can novel functionalities be realized, such as the ability to bring objects at different distances from the lens into sharp focus, but they also alleviate the stringent constraints associated with precisely aligning lenses on top of a sensor. Extending the DOF also increases the spectral operation range because it reduces the impact of the focal shift due to chromatic aberrations on the quality of the images obtained³⁵.

Conventional techniques, including cubic, shifted axicon, log-asphere, and SQUBIC lenses, have been proposed for synthesizing phase masks to achieve EDOF properties³⁶. However, these traditional synthesis approaches encounter difficulties in striking a balance between maximizing focal efficiency and extending the focal depth. The limitations of the lookup table approach become apparent in addressing the complex requirements of achieving the desired EDOF properties. Notably, inverse design methods based on traditional optimization mechanisms have been utilized to maximize EDOF metalens performance, particularly in relying on classical optimization algorithm³⁷,³⁸ with numerical aperture not exceeding 0.44 accompanied with huge number of simulations.

The proposed design for this application involves a cylindrical metalens composed of unit cells described in Fig. 1a. The lens comprises 18 unit cells aligned along the x-axis with periodic boundary conditions along the y-axis to mimic 1D metalens. We start from a design with a NA of 0.6 and we extend the focal depth in a rectangle from the $f$ - DOF to $f + 5$DOF with a width of 2 FWHM where $f$ is the focal distance and DOF is given by the formula in³⁹.

In order to extend the depth of focus and achieve a uniform optical field distribution within the target depth, we define the objective function as the geometric mean of the mean power and the minimum power obtained from 20 equally spaced surfaces above and below the focal point. Using a minimum in the objective is primordial in order to obtain a design with a homogeneous intensity in the target region. This objective function allows us to simultaneously increase the focal depth along the z-axis, decrease the FWHM along the x-axis, and enhance the focusing power efficiency.

The employed geometric parametrization in this study aligns with the methodology presented in Section 2.1, which encompasses the diameter and distances between the resonators. However, to enhance the degrees of freedom, we also explored the variation of the resonator height, ranging from 150 nm to 250 nm, resulting in the optimization of a total of 19 parameters. The optimization process, as depicted in Fig. 3a, initiates with a Design of Experiments (DOE) consisting of 90 designs, indicated by the blue points. The process concludes after 450 simulations, represented by the black points in Fig. 3a. The optimal design achieved a performance of 30% according to our objective function. Furthermore, the electric field intensity ($|\mathbf{E}|^2$) demonstrates an uniform distribution across the entire focal zone, as shown in the first sub-figure of Fig. 3b.

<figure><img src="image_4.png" /><figcaption>**Figure 2.** Fields maps of $\Re(E_y)$ around the last meta-element in the supercell of the optimized grating shown in Fig. 1f when surrounded with periodic boundary conditions or inside the grating. The black rectangle delimits the meta-atoms boundaries. When there is a translational symmetry only the electric and magnetic dipoles are exited while when this symmetry is broken due to the strong coupling between elements, higher order modes can be excited.</figcaption></figure>
<figure><img src="image_5.png" /></figure>

<figure><img src="image_6.png" /></figure>

<figure><img src="image_7.png" /></figure>

<figure><img src="image_8.png" /><figcaption>**Figure 3.** Optimization results for the EDOF metalens. (a) Objective function values as a function of the number of simulations. The insight represents the geometry of half of the metalens in the *x* direction. (b) Field maps of ||*E*||² at 1064 nm of the optimized design with fixed and varying *h*, the parabolic designs and of the log-asphere design. The green rectangle represents the region in which we try to increase the DOF. (c) Plot of the minimum and maximum focus efficiency over 20 equally spaced rectangles in the target region. The focus efficiency is defined as the part the incident power that is transmitted in a focus plane. The focus efficiency for the optimized design with a varying *h* varies from 26 to 35% while for the parabolic designs it varies between 0.05 and 27%. (d) Optical Transfer Functions (OTFs) of each designs at *f*, *f* + 2.5 DOF and *f* + 5 DOF. Optimized designs preserve spatial frequencies up to 600 cycle/mm at the three locations whereas the classical ones both fail at *f* + 5 DOF.</figcaption></figure>

The optimization process converged to a specific height value *h* = 189.5 nm, which differs from the height utilized in the lookup table presented in Fig. 1a. To further investigate this aspect, we performed another optimization run while keeping the height fixed at *h* = 170 nm, as defined in Fig. 1a.

The field profile distribution in the case with a constant height exhibits inferior performance compared to the case with varying height, as depicted in the first two sub-figures of Fig. 3b. Moreover, the traditional parabolic phase configuration yields a narrow focal spot that fails to cover the full depth of focus. In order to provide a comprehensive comparison, the log-sphere synthesizing approach is also considered. The central anular zone and the outer-most one of the log-asphere phase profile are chosen to be respectively *f* − DOF and *f* + 5 DOF. However, this approach only achieves an EDOF of 4 µm accompanied by a significant decrease in power efficiency, as indicated in Fig. 3b. In conclusion, the optimized design with varying height demonstrates an EDOF of approximately 20 µm, which is five times larger than that of the classical design. Additionally, Fig. 3c illustrates the enhanced focusing efficiency of the optimal design, surpassing that of the classical designs by a factor of three in average. This gain is greatly due to a large reduction of the back-scattering of the structure. 50% of the losses of all designs are due to the small aperture of the lenses, which leads to unavoidable diffraction (see Supplementary Informations 1.4 for more details).

Furthermore, in Fig. 3d, we present the Optical Transfert Function (OTF) at three distinct *z* values: the focal point, the final positions within the focal depth zone and a position between the two. Notably, the optimized designs OTFs keep a similar shape at the three locations while the classical ones can not capture spacial frequencies higher than 100 cycle/mm at the second location. This implies that our design demonstrates significantly improved tolerance to changes in the focal plane position while maintaining nearly the same power efficiency and optical resolution. These results confirm that the local lookup table approach is unsuitable for complex metalens objectives due to the strong near-field coupling, which is a crucial factor for achieving the desired objectives. Importantly, our optimization process only required 450 iterations to converge to the globally optimal design, considering 19 parameters. This is made possible by the use of the EGO algorithm, which allows us to achieve
<figure><img src="image_5.png" /></figure>

<figure><img src="image_6.png" /></figure>

<figure><img src="image_7.png" /></figure>

<figure><img src="image_8.png" /><figcaption>**Figure 3.** Optimization results for the EDOF metalens. (a) Objective function values as a function of the number of simulations. The insight represents the geometry of half of the metalens in the *x* direction. (b) Field maps of ||*E*||² at 1064 nm of the optimized design with fixed and varying *h*, the parabolic designs and of the log-asphere design. The green rectangle represents the region in which we try to increase the DOF. (c) Plot of the minimum and maximum focus efficiency over 20 equally spaced rectangles in the target region. The focus efficiency is defined as the part the incident power that is transmitted in a focus plane. The focus efficiency for the optimized design with a varying *h* varies from 26 to 35% while for the parabolic designs it varies between 0.05 and 27%. (d) Optical Transfer Functions (OTFs) of each designs at *f*, *f* + 2.5 DOF and *f* + 5 DOF. Optimized designs preserve spatial frequencies up to 600 cycle/mm at the three locations whereas the classical ones both fail at *f* + 5 DOF.</figcaption></figure>

The optimization process converged to a specific height value *h* = 189.5 nm, which differs from the height utilized in the lookup table presented in Fig. 1a. To further investigate this aspect, we performed another optimization run while keeping the height fixed at *h* = 170 nm, as defined in Fig. 1a.

The field profile distribution in the case with a constant height exhibits inferior performance compared to the case with varying height, as depicted in the first two sub-figures of Fig. 3b. Moreover, the traditional parabolic phase configuration yields a narrow focal spot that fails to cover the full depth of focus. In order to provide a comprehensive comparison, the log-sphere synthesizing approach is also considered. The central anular zone and the outer-most one of the log-asphere phase profile are chosen to be respectively *f* − DOF and *f* + 5 DOF. However, this approach only achieves an EDOF of 4 µm accompanied by a significant decrease in power efficiency, as indicated in Fig. 3b. In conclusion, the optimized design with varying height demonstrates an EDOF of approximately 20 µm, which is five times larger than that of the classical design. Additionally, Fig. 3c illustrates the enhanced focusing efficiency of the optimal design, surpassing that of the classical designs by a factor of three in average. This gain is greatly due to a large reduction of the back-scattering of the structure. 50% of the losses of all designs are due to the small aperture of the lenses, which leads to unavoidable diffraction (see Supplementary Informations 1.4 for more details).

Furthermore, in Fig. 3d, we present the Optical Transfert Function (OTF) at three distinct *z* values: the focal point, the final positions within the focal depth zone and a position between the two. Notably, the optimized designs OTFs keep a similar shape at the three locations while the classical ones can not capture spacial frequencies higher than 100 cycle/mm at the second location. This implies that our design demonstrates significantly improved tolerance to changes in the focal plane position while maintaining nearly the same power efficiency and optical resolution. These results confirm that the local lookup table approach is unsuitable for complex metalens objectives due to the strong near-field coupling, which is a crucial factor for achieving the desired objectives. Importantly, our optimization process only required 450 iterations to converge to the globally optimal design, considering 19 parameters. This is made possible by the use of the EGO algorithm, which allows us to achieve
convergence with a significantly lower number of iterations compared to standard state-of-the-art optimization algorithms³⁷,³⁸.

## Nonlocal resonant metasurface in reflection

In this section, we aim to provide another description of a nonlocal metasurface example that cannot be designed efficiently using the lookup table synthesizing local approach. Specifically, we consider the classical 2D metagrating configuration to design a beam steering device in reflection. To achieve high performance, we rely on the asymmetric-cavity resonator, also known as the Gires-Tournois cavity, to design a reflective metasurface that supports a nonlinear phase response associated with a high reflection amplitude response. The single-unit cell, depicted in the inset of Fig. 4a, comprises a 2D rectangular TiO₂ ridge that is separated from a Silver (Ag) layer by a distance denoted as *D*. In this configuration, the distance *D* is the crucial parameter required to achieve the ideal performance Gires-Tournois condition by manipulating the reflection zeros to obtain full phase modulation with unity reflectivity. Further details on this approach are available in²⁴. Notably, the height of the ridge *H* is chosen as 400 nm, while the distance from the Ag layer is kept constant at *D* = 140 nm in Fig. 4a. The design is configured to operate at a wavelength *λ* = 600 nm. To achieve a step nonlinear phase response with more

<figure><img src="image_9.png" /></figure>

**Figure 4.** The performance of the nonlocal metasurface device with a nonlinear phase response and high reflection amplitude response. (a) Phase and reflection response for single unit cell simulations as a function of the width of the TiO₂ nano ridges. The blue curve and left label represent the constant reflectivity of the device, while the orange curve and right label represent the nonlinear phase response of the device. The inset in (a) indicates the schematic representation of the structure under consideration. We consider a TiO₂ nano ridge (invariant along *y* direction) encapsulated in SiO₂ material. The refractive index of the Ag layer is obtained from Ref.⁴⁰, while the refractive index of TiO₂, and SiO₂ are fixed as 2.5 and 1.45, respectively. The height of the ridge and the period are fixed as *H* = 400 nm and *P* = 300 nm, respectively. The width of the pillar can vary from 130 nm to 220 nm. (b) Beam steering configuration composed of 3 different ridges with varying distances in between. (c) Beam steering response (the power reflected off in the first order mode) based on different optimization scenarios. The blue curve refers to the classical local synthesizing approach where the widths of the ridges are varying in a gradient manner based on the look-up table given in (a), for a fixed period *P* = 300 nm. The orange curve indicates the nonlocal optimization by varying the widths and distances in between, while keeping *H* = 400 nm and *D* = 140 nm as in (a). The green curve is for the full nonlocal optimization, where *H* and *D* are varying along with the widths and distances between ridges. (d) Field profiles for the optimized scenarios depicted in (c) at the desired wavelength *λ* = 600 nm. The optimized parameters in each case are depicted in Table 3 in the the supplementary information section.
convergence with a significantly lower number of iterations compared to standard state-of-the-art optimization algorithms³⁷,³⁸.

## Nonlocal resonant metasurface in reflection

In this section, we aim to provide another description of a nonlocal metasurface example that cannot be designed efficiently using the lookup table synthesizing local approach. Specifically, we consider the classical 2D metagrating configuration to design a beam steering device in reflection. To achieve high performance, we rely on the asymmetric-cavity resonator, also known as the Gires-Tournois cavity, to design a reflective metasurface that supports a nonlinear phase response associated with a high reflection amplitude response. The single-unit cell, depicted in the inset of Fig. 4a, comprises a 2D rectangular TiO₂ ridge that is separated from a Silver (Ag) layer by a distance denoted as *D*. In this configuration, the distance *D* is the crucial parameter required to achieve the ideal performance Gires-Tournois condition by manipulating the reflection zeros to obtain full phase modulation with unity reflectivity. Further details on this approach are available in²⁴. Notably, the height of the ridge *H* is chosen as 400 nm, while the distance from the Ag layer is kept constant at *D* = 140 nm in Fig. 4a. The design is configured to operate at a wavelength *λ* = 600 nm. To achieve a step nonlinear phase response with more

<figure><img src="image_9.png" /></figure>

**Figure 4.** The performance of the nonlocal metasurface device with a nonlinear phase response and high reflection amplitude response. (a) Phase and reflection response for single unit cell simulations as a function of the width of the TiO₂ nano ridges. The blue curve and left label represent the constant reflectivity of the device, while the orange curve and right label represent the nonlinear phase response of the device. The inset in (a) indicates the schematic representation of the structure under consideration. We consider a TiO₂ nano ridge (invariant along *y* direction) encapsulated in SiO₂ material. The refractive index of the Ag layer is obtained from Ref.⁴⁰, while the refractive index of TiO₂, and SiO₂ are fixed as 2.5 and 1.45, respectively. The height of the ridge and the period are fixed as *H* = 400 nm and *P* = 300 nm, respectively. The width of the pillar can vary from 130 nm to 220 nm. (b) Beam steering configuration composed of 3 different ridges with varying distances in between. (c) Beam steering response (the power reflected off in the first order mode) based on different optimization scenarios. The blue curve refers to the classical local synthesizing approach where the widths of the ridges are varying in a gradient manner based on the look-up table given in (a), for a fixed period *P* = 300 nm. The orange curve indicates the nonlocal optimization by varying the widths and distances in between, while keeping *H* = 400 nm and *D* = 140 nm as in (a). The green curve is for the full nonlocal optimization, where *H* and *D* are varying along with the widths and distances between ridges. (d) Field profiles for the optimized scenarios depicted in (c) at the desired wavelength *λ* = 600 nm. The optimized parameters in each case are depicted in Table 3 in the the supplementary information section.
than 2π phase modulation, we vary the width of the TiO₂ ridges from 130 nm to 220 nm while maintaining a constant reflectivity (see blue curve and left label in Fig. 4a). Consequently, the single-unit cell design (with periodic boundary conditions on the slides) indicates that the resulting configuration displays a highly nonlinear phase response with a high reflection amplitude response, indicating that the device can be used effectively for beam steering purposes.

We began our investigation by examining a supercell that consists of three ridges with different widths, arranged in a specific manner to steer reflected light in a SiO₂ layer at a particular angle θ = 27.37°. The supercell configuration is depicted in Fig. 4b. Initially, we designed a classical local configuration, where the *H*, *D*, and *P* were fixed at 400 nm, 140 nm, and 300 nm, respectively, as shown in the lookup table provided in Fig. 4a. The widths of the three ridges were selected to provide a full 2π phase variation in a gradient manner with a step of 120 deg. The corresponding results are shown by the blue curve in Fig. 4a, where the total power reflected in the first order mode did not exceed 25%.

However, when we considered the distances between the ridges as optimization parameters in addition to the widths, we observed a significant improvement in the performance of the metasurface. This resulted in a total reflected power in the first order mode of 75% at the desired wavelength, as indicated by the orange curve in Fig. 4c. This huge improvement in the performance compared to the local approximation was due to the resonance behavior and the strong nonlocal interactions that were accurately accounted for by varying the distances between the ridges, as shown in the field profile given in the first two columns in Fig. 4d. Furthermore, we found that we could enhance the efficiency of the metasurface design by including both *H* and *D* in the nonlocal optimization. In this case, the efficiency reaches 85% (see green curve and field map in Fig. 4c and d, respectively). Interestingly, we also obtain a novel set of parameters where the *D* value became 100 nm instead of 140 nm, while the optimal height was *H* = 400 nm. We obtain the same conclusion when we increased the number of ridges, as depicted in the supplementary information Fig. 3, where the results based on the classical local lookup table provides lower performance compared to the nonlocal optimization scenario.

## Conclusion

In conclusion, our study highlights the importance of considering nonlocal effects and near-field coupling in the design of resonant metasurfaces. The limitations of using a local single unit-cell approximation with periodic boundary conditions have been demonstrated, emphasizing the need for more accurate predictions of metasurface behavior. The nonlocal distribution of fields plays a crucial role in determining key parameters such as resonance wavelength, quality factor, and farfield response of the metasurface. Therefore, relying solely on local lookup tables for metasurface design is insufficient.

To push this limit, we employed a global advanced statistical learning optimization technique, which ensure convergence to the global design using a moderate number of solver calls compared to traditional optimization algorithms. Indeed, state-of-the art optimization algorithms as evolutionary algorithm and local gradient can be utilized, but to the extent of requiring prohibitive number of simulations, getting stuck in a local design or requiring a careful starting point. Three different resonant metasurfaces were considered in this study. The first one operates in transmission and utilizes the properties of resonant Huygens metasurfaces. Numerical results indicate that a beam steering metasurface achieved an efficiency of 80%, surpassing the classical design with only 23% efficiency. Additionally, an extended depth-of-focus (EDOF) metalens was optimized, resulting in a five fold increase in focal depth and a four fold enhancement in focusing power compared to classical designs.

Furthermore, the performance of wavelength-selected metagratings in reflection was investigated based on asymmetric cavity considerations. Notably, non-intuitive designs were obtained, surpassing 85% efficiency, which far exceeds the efficiency achieved with classical gradient phase distribution.

The presence of nonlocal effects arising from the interaction between neighboring unit cells significantly alters the behavior of the fields. A comprehensive understanding of this phenomenon is crucial for designing high-performance resonant metasurfaces. This knowledge empowers researchers and engineers to create metasurfaces with enhanced performance, opening avenues for novel applications in various domains.

## Data availability

The data that support the findings of this study are available from the authors upon reasonable request.

Received: 3 October 2023; Accepted: 10 January 2024

Published online: 18 January 2024

## References

1. Yu, N. et al. Light propagation with phase discontinuities: Generalized laws of reflection and refraction. *Science* **334**, 333–337 (2011).

2. Genevet, P., Capasso, F., Aieta, F., Khorasaninejad, M. & Devlin, R. Recent advances in planar optics: From plasmonic to dielectric metasurfaces. *Optica* **4**, 139–152 (2017).

3. Astilean, S., Lalanne, P., Chavel, P., Cambril, E. & Launois, H. High-efficiency subwavelength diffractive element patterned in a high-refractive-index material for 633 nm. *Opt. Lett.* **23**, 552–554 (1998).

4. Lin, D., Fan, P., Hasman, E. & Brongersma, M. L. Dielectric gradient metasurface optical elements. *Science* **345**, 298–302 (2014).

5. Elsawy, M. M. et al. Multiobjective statistical learning optimization of rgb metalens. *ACS Photon. *8*, 2498–2508 (2021).

6. Wang, S. et al. A broadband achromatic metalens in the visible. *Nat. Nanotechnol.* **13**, 227–232 (2018).

7. Chen, W. T. et al. A broadband achromatic Metalens for focusing and imaging in the visible. *Nat. Nanotechnol.* **13**, 220–226 (2018).

8. Pan, M. et al. Dielectric metalens for miniaturized imaging systems: progress and challenges. *Light: Sci. Appl.* **11**, 1–32 (2022).

9. Song, Q. et al. Ptychography retrieval of fully polarized holograms from geometric-phase metasurfaces. *Nat. Commun.* **11**, 1–8 (2020).
than 2π phase modulation, we vary the width of the TiO₂ ridges from 130 nm to 220 nm while maintaining a constant reflectivity (see blue curve and left label in Fig. 4a). Consequently, the single-unit cell design (with periodic boundary conditions on the slides) indicates that the resulting configuration displays a highly nonlinear phase response with a high reflection amplitude response, indicating that the device can be used effectively for beam steering purposes.

We began our investigation by examining a supercell that consists of three ridges with different widths, arranged in a specific manner to steer reflected light in a SiO₂ layer at a particular angle θ = 27.37°. The supercell configuration is depicted in Fig. 4b. Initially, we designed a classical local configuration, where the *H*, *D*, and *P* were fixed at 400 nm, 140 nm, and 300 nm, respectively, as shown in the lookup table provided in Fig. 4a. The widths of the three ridges were selected to provide a full 2π phase variation in a gradient manner with a step of 120 deg. The corresponding results are shown by the blue curve in Fig. 4a, where the total power reflected in the first order mode did not exceed 25%.

However, when we considered the distances between the ridges as optimization parameters in addition to the widths, we observed a significant improvement in the performance of the metasurface. This resulted in a total reflected power in the first order mode of 75% at the desired wavelength, as indicated by the orange curve in Fig. 4c. This huge improvement in the performance compared to the local approximation was due to the resonance behavior and the strong nonlocal interactions that were accurately accounted for by varying the distances between the ridges, as shown in the field profile given in the first two columns in Fig. 4d. Furthermore, we found that we could enhance the efficiency of the metasurface design by including both *H* and *D* in the nonlocal optimization. In this case, the efficiency reaches 85% (see green curve and field map in Fig. 4c and d, respectively). Interestingly, we also obtain a novel set of parameters where the *D* value became 100 nm instead of 140 nm, while the optimal height was *H* = 400 nm. We obtain the same conclusion when we increased the number of ridges, as depicted in the supplementary information Fig. 3, where the results based on the classical local lookup table provides lower performance compared to the nonlocal optimization scenario.

## Conclusion

In conclusion, our study highlights the importance of considering nonlocal effects and near-field coupling in the design of resonant metasurfaces. The limitations of using a local single unit-cell approximation with periodic boundary conditions have been demonstrated, emphasizing the need for more accurate predictions of metasurface behavior. The nonlocal distribution of fields plays a crucial role in determining key parameters such as resonance wavelength, quality factor, and farfield response of the metasurface. Therefore, relying solely on local lookup tables for metasurface design is insufficient.

To push this limit, we employed a global advanced statistical learning optimization technique, which ensure convergence to the global design using a moderate number of solver calls compared to traditional optimization algorithms. Indeed, state-of-the art optimization algorithms as evolutionary algorithm and local gradient can be utilized, but to the extent of requiring prohibitive number of simulations, getting stuck in a local design or requiring a careful starting point. Three different resonant metasurfaces were considered in this study. The first one operates in transmission and utilizes the properties of resonant Huygens metasurfaces. Numerical results indicate that a beam steering metasurface achieved an efficiency of 80%, surpassing the classical design with only 23% efficiency. Additionally, an extended depth-of-focus (EDOF) metalens was optimized, resulting in a five fold increase in focal depth and a four fold enhancement in focusing power compared to classical designs.

Furthermore, the performance of wavelength-selected metagratings in reflection was investigated based on asymmetric cavity considerations. Notably, non-intuitive designs were obtained, surpassing 85% efficiency, which far exceeds the efficiency achieved with classical gradient phase distribution.

The presence of nonlocal effects arising from the interaction between neighboring unit cells significantly alters the behavior of the fields. A comprehensive understanding of this phenomenon is crucial for designing high-performance resonant metasurfaces. This knowledge empowers researchers and engineers to create metasurfaces with enhanced performance, opening avenues for novel applications in various domains.

## Data availability

The data that support the findings of this study are available from the authors upon reasonable request.

Received: 3 October 2023; Accepted: 10 January 2024

Published online: 18 January 2024

## References

1. Yu, N. et al. Light propagation with phase discontinuities: Generalized laws of reflection and refraction. *Science* **334**, 333–337 (2011).

2. Genevet, P., Capasso, F., Aieta, F., Khorasaninejad, M. & Devlin, R. Recent advances in planar optics: From plasmonic to dielectric metasurfaces. *Optica* **4**, 139–152 (2017).

3. Astilean, S., Lalanne, P., Chavel, P., Cambril, E. & Launois, H. High-efficiency subwavelength diffractive element patterned in a high-refractive-index material for 633 nm. *Opt. Lett.* **23**, 552–554 (1998).

4. Lin, D., Fan, P., Hasman, E. & Brongersma, M. L. Dielectric gradient metasurface optical elements. *Science* **345**, 298–302 (2014).

5. Elsawy, M. M. et al. Multiobjective statistical learning optimization of rgb metalens. *ACS Photon. *8*, 2498–2508 (2021).

6. Wang, S. et al. A broadband achromatic metalens in the visible. *Nat. Nanotechnol.* **13**, 227–232 (2018).

7. Chen, W. T. et al. A broadband achromatic Metalens for focusing and imaging in the visible. *Nat. Nanotechnol.* **13**, 220–226 (2018).

8. Pan, M. et al. Dielectric metalens for miniaturized imaging systems: progress and challenges. *Light: Sci. Appl.* **11**, 1–32 (2022).

9. Song, Q. et al. Ptychography retrieval of fully polarized holograms from geometric-phase metasurfaces. *Nat. Commun.* **11**, 1–8 (2020).
10. Xie, Y.-Y. et al. Metasurface-integrated vertical cavity surface-emitting lasers for programmable directional lasing emissions. *Nat. Nanotechnol.* **15**, 125–130 (2020).

11. Lalanne, P., Astilean, S., Chavel, P., Cambril, E. & Launois, H. Design and fabrication of blazed binary diffractive elements with sampling periods smaller than the structural cutoff. *JOSA A* **16**, 1143–1156 (1999).

12. Khorasaninejad, M. & Capasso, F. Broadband multifunctional efficient meta-gratings based on dielectric waveguide phase shifters. *Nano Lett.* **15**, 6709–6715 (2015).

13. Bomzon, Z., Biener, G., Kleiner, V. & Hasman, E. Space-variant Pancharatnam-berry phase optical elements with computer-generated subwavelength gratings. *Opt. Lett.* **27**, 1141–1143 (2002).

14. Khorasaninejad, M. & Crozier, K. B. Silicon nanofin grating as a miniature chirality-distinguishing beam-splitter. *Nat. Commun.* **5**, 1–6 (2014).

15. Khorasaninejad, M. et al. Metalenses at visible wavelengths: Diffraction-limited focusing and subwavelength resolution imaging. *Science* **352**, 1190–1194 (2016).

16. Tian, Y. et al. Manipulation of the arbitrary scattering angle based on all-dielectric transmissive pancharatnam berry phase coding metasurfaces in the visible range. *Opt. Express* **28**, 32107–32123 (2020).

17. Overvig, A. C., Malek, S. C. & Yu, N. Multifunctional nonlocal metasurfaces. *Phys. Rev. Lett.* **125**, 017402 (2020).

18. Cai, H. et al. Inverse design of metasurfaces with non-local interactions. *npj Comput. Mater.* **6**, 1–8 (2020).

19. Liu, W. & Kivshar, Y. S. Generalized Kerkер effects in nanophotonics and meta-optics. *Opt. Express* **26**, 13085–13105 (2018).

20. Yu, Y. F. et al. High-transmission dielectric metasurface with 2π phase control at visible wavelengths. *Laser Photon. Rev.* **9**, 412–418 (2015).

21. Huang, L. et al. Three-dimensional optical holography using a plasmonic metasurface. *Nat. Commun.* **4**, 1–8 (2013).

22. Bin-Alam, M. S. et al. Ultra-high-q resonances in plasmonic metasurfaces. *Nat. Commun.* **12**, 1–8 (2021).

23. Colom, R. et al. Crossing of the branch cut: The topological origin of a universal 2π-phase retardation in non-Bermitian metasurfaces. *Laser Photon. Rev.* **2200976** (2023).

24. Elsawy, M. et al. Universal active metasurfaces for ultimate wavefront molding by manipulating the reflection singularities. *Laser Photon. Rev.* **2200880** (2022).

25. Overvig, A. C., Malek, S. C., Carter, M. J., Shrestha, S. & Yu, N. Selection rules for quasibound states in the continuum. *Phys. Rev. B* **102**, 035434 (2020).

26. Malek, S. C., Overvig, A. C., Shrestha, S. & Yu, N. Active nonlocal metasurfaces. *Nanophotonics* **10**, 655–665 (2021).

27. Elsawy, M. M. et al. Global optimization of metasurface designs using statistical learning methods. *Sci. Rep.* **9**, 1–15 (2019).

28. Elsawy, M. M., Binois, M., Duvigneau, R., Lanteri, S. & Genevet, P. Optimization of metasurfaces under geometrical uncertainty using statistical learning. *Opt. Express* **29**, 29887–29898 (2021).

29. Elsawy, M. M., Lanteri, S., Duvigneau, R., Fan, J. A. & Genevet, P. Numerical optimization methods for metasurfaces. *Laser Photon. Rev.* **14**, 1900445 (2020).

30. Koshelev, K. & Kivshar, Y. Dielectric resonant metaphotonics. *Acs Photon.* **8**, 102–112 (2020).

31. Rybin, M. V. et al. High-q supercavity modes in subwavelength dielectric resonators. *Phys. Rev. Lett.* **119**, 243901 (2017).

32. Özdemir, A., Hayran, Z., Takashima, Y. & Kurt, H. Polarization independent high transmission large numerical aperture laser beam focusing and deflection by dielectric huygens' metasurfaces. *Opt. Commun.* **401**, 46–53 (2017).

33. Gigli, C. et al. Fundamental limitations of Huygens' metasurfaces for optical beam shaping. *Laser Photon. Rev.* **15**, 2000448 (2021).

34. DIOGENeS: A discontinuous-galerkin based software suite for nano-optics. https://diogenes.inria.fr.

35. Baranikov, V. V. et al. Large field-of-view and multi-color imaging with gap quadratic metalenses. arXiv:2305.16676 (2023).

36. Huang, L., Whitehead, J., Colburn, S. & Majumdar, A. Design and analysis of extended depth of focus metalenses for achromatic computational imaging. *Photon. Res.* **8**, 1613–1623 (2020).

37. Zheng, Y. et al. Designing high-efficiency extended depth-of-focus metalens via topology-shape optimization. *Nanophotonics* **11**, 2967–2975 (2022).

38. Xiao, S. et al. Inverse design of a near-infrared metalens with an extended depth of focus based on double-process genetic algorithm optimization. *Opt. Express* **31**, 8668–8681 (2023).

39. Bayati, E. et al. Inverse designed metalenses with extended depth of focus. *ACS Photon.* **7**, 873–878 (2020).

40. Johnson, P. B. & Christy, R.-W. Optical constants of the noble metals. *Phys. Rev. B* **6**, 4370 (1972).

## Acknowledgements

This work was supported by the French government, through the UCAJEDI Investments in the Future project managed by the National Research Agency (ANR) under reference number ANR-15-IDEX-01. The authors are grateful to the OPAL infrastructure and the Université Côte d'Azur's Center for High-Performance Computing for providing resources and support. Some of the experiments presented in this paper were also carried out using the PlaFRIM experimental testbed, supported by Inria, CNRS (LABRI and IMB), Université de Bordeaux, Bordeaux INP and Conseil Régional d'Aquitaine. This project was also provided with computer and storage resources by GENCI at IDRIS thanks to the grant 2023-A0130610263 on the supercomputer Jean Zay's the CSL partition. We gratefully acknowledge the French National Association for Research and Technology (ANRT, CIFRE grant number 2022/0739) for their financial support.

## Author contributions

M.E. proposed the structures. E.I. and M.E. performed the DGTD simulations and the optimizations of the structures. All authors analyzed the simulations and optimization results. E.I. and M.E. wrote the manuscript. All authors reviewed the manuscript.

## Competing interests

The authors declare no competing interests.

## Additional information

**Supplementary Information** The online version contains supplementary material available at https://doi.org/10.1038/s41598-024-51898-8.

**Correspondence** and requests for materials should be addressed to M.E.

**Reprints and permissions information** is available at www.nature.com/reprints.
10. Xie, Y.-Y. et al. Metasurface-integrated vertical cavity surface-emitting lasers for programmable directional lasing emissions. *Nat. Nanotechnol.* **15**, 125–130 (2020).

11. Lalanne, P., Astilean, S., Chavel, P., Cambril, E. & Launois, H. Design and fabrication of blazed binary diffractive elements with sampling periods smaller than the structural cutoff. *JOSA A* **16**, 1143–1156 (1999).

12. Khorasaninejad, M. & Capasso, F. Broadband multifunctional efficient meta-gratings based on dielectric waveguide phase shifters. *Nano Lett.* **15**, 6709–6715 (2015).

13. Bomzon, Z., Biener, G., Kleiner, V. & Hasman, E. Space-variant Pancharatnam-berry phase optical elements with computer-generated subwavelength gratings. *Opt. Lett.* **27**, 1141–1143 (2002).

14. Khorasaninejad, M. & Crozier, K. B. Silicon nanofin grating as a miniature chirality-distinguishing beam-splitter. *Nat. Commun.* **5**, 1–6 (2014).

15. Khorasaninejad, M. et al. Metalenses at visible wavelengths: Diffraction-limited focusing and subwavelength resolution imaging. *Science* **352**, 1190–1194 (2016).

16. Tian, Y. et al. Manipulation of the arbitrary scattering angle based on all-dielectric transmissive pancharatnam berry phase coding metasurfaces in the visible range. *Opt. Express* **28**, 32107–32123 (2020).

17. Overvig, A. C., Malek, S. C. & Yu, N. Multifunctional nonlocal metasurfaces. *Phys. Rev. Lett.* **125**, 017402 (2020).

18. Cai, H. et al. Inverse design of metasurfaces with non-local interactions. *npj Comput. Mater.* **6**, 1–8 (2020).

19. Liu, W. & Kivshar, Y. S. Generalized Kerkер effects in nanophotonics and meta-optics. *Opt. Express* **26**, 13085–13105 (2018).

20. Yu, Y. F. et al. High-transmission dielectric metasurface with 2π phase control at visible wavelengths. *Laser Photon. Rev.* **9**, 412–418 (2015).

21. Huang, L. et al. Three-dimensional optical holography using a plasmonic metasurface. *Nat. Commun.* **4**, 1–8 (2013).

22. Bin-Alam, M. S. et al. Ultra-high-q resonances in plasmonic metasurfaces. *Nat. Commun.* **12**, 1–8 (2021).

23. Colom, R. et al. Crossing of the branch cut: The topological origin of a universal 2π-phase retardation in non-Bermitian metasurfaces. *Laser Photon. Rev.* **2200976** (2023).

24. Elsawy, M. et al. Universal active metasurfaces for ultimate wavefront molding by manipulating the reflection singularities. *Laser Photon. Rev.* **2200880** (2022).

25. Overvig, A. C., Malek, S. C., Carter, M. J., Shrestha, S. & Yu, N. Selection rules for quasibound states in the continuum. *Phys. Rev. B* **102**, 035434 (2020).

26. Malek, S. C., Overvig, A. C., Shrestha, S. & Yu, N. Active nonlocal metasurfaces. *Nanophotonics* **10**, 655–665 (2021).

27. Elsawy, M. M. et al. Global optimization of metasurface designs using statistical learning methods. *Sci. Rep.* **9**, 1–15 (2019).

28. Elsawy, M. M., Binois, M., Duvigneau, R., Lanteri, S. & Genevet, P. Optimization of metasurfaces under geometrical uncertainty using statistical learning. *Opt. Express* **29**, 29887–29898 (2021).

29. Elsawy, M. M., Lanteri, S., Duvigneau, R., Fan, J. A. & Genevet, P. Numerical optimization methods for metasurfaces. *Laser Photon. Rev.* **14**, 1900445 (2020).

30. Koshelev, K. & Kivshar, Y. Dielectric resonant metaphotonics. *Acs Photon.* **8**, 102–112 (2020).

31. Rybin, M. V. et al. High-q supercavity modes in subwavelength dielectric resonators. *Phys. Rev. Lett.* **119**, 243901 (2017).

32. Özdemir, A., Hayran, Z., Takashima, Y. & Kurt, H. Polarization independent high transmission large numerical aperture laser beam focusing and deflection by dielectric huygens' metasurfaces. *Opt. Commun.* **401**, 46–53 (2017).

33. Gigli, C. et al. Fundamental limitations of Huygens' metasurfaces for optical beam shaping. *Laser Photon. Rev.* **15**, 2000448 (2021).

34. DIOGENeS: A discontinuous-galerkin based software suite for nano-optics. https://diogenes.inria.fr.

35. Baranikov, V. V. et al. Large field-of-view and multi-color imaging with gap quadratic metalenses. arXiv:2305.16676 (2023).

36. Huang, L., Whitehead, J., Colburn, S. & Majumdar, A. Design and analysis of extended depth of focus metalenses for achromatic computational imaging. *Photon. Res.* **8**, 1613–1623 (2020).

37. Zheng, Y. et al. Designing high-efficiency extended depth-of-focus metalens via topology-shape optimization. *Nanophotonics* **11**, 2967–2975 (2022).

38. Xiao, S. et al. Inverse design of a near-infrared metalens with an extended depth of focus based on double-process genetic algorithm optimization. *Opt. Express* **31**, 8668–8681 (2023).

39. Bayati, E. et al. Inverse designed metalenses with extended depth of focus. *ACS Photon.* **7**, 873–878 (2020).

40. Johnson, P. B. & Christy, R.-W. Optical constants of the noble metals. *Phys. Rev. B* **6**, 4370 (1972).

## Acknowledgements

This work was supported by the French government, through the UCAJEDI Investments in the Future project managed by the National Research Agency (ANR) under reference number ANR-15-IDEX-01. The authors are grateful to the OPAL infrastructure and the Université Côte d'Azur's Center for High-Performance Computing for providing resources and support. Some of the experiments presented in this paper were also carried out using the PlaFRIM experimental testbed, supported by Inria, CNRS (LABRI and IMB), Université de Bordeaux, Bordeaux INP and Conseil Régional d'Aquitaine. This project was also provided with computer and storage resources by GENCI at IDRIS thanks to the grant 2023-A0130610263 on the supercomputer Jean Zay's the CSL partition. We gratefully acknowledge the French National Association for Research and Technology (ANRT, CIFRE grant number 2022/0739) for their financial support.

## Author contributions

M.E. proposed the structures. E.I. and M.E. performed the DGTD simulations and the optimizations of the structures. All authors analyzed the simulations and optimization results. E.I. and M.E. wrote the manuscript. All authors reviewed the manuscript.

## Competing interests

The authors declare no competing interests.

## Additional information

**Supplementary Information** The online version contains supplementary material available at https://doi.org/10.1038/s41598-024-51898-8.

**Correspondence** and requests for materials should be addressed to M.E.

**Reprints and permissions information** is available at www.nature.com/reprints.
**Publisher's note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_10.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International

License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2024
**Publisher's note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_10.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International

License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2024