---
title: Nonlinear Nonlocal Metasurfaces
url: https://www.emergentmind.com/topics/nonlinear-nonlocal-metasurface
type: topic
---

# Nonlinear Nonlocal Metasurfaces

Nonlinear nonlocal metasurfaces are resonant planar nanostructures in which nonlinear optical processes are mediated by collective, lattice-enabled modes rather than only by isolated local resonances. In the recent literature, this category includes guided-mode-resonant lithium niobate gratings for sum-frequency up-conversion imaging, silicon nonlocal phase-gradient metasurfaces based on quasi–bound states in the continuum for third-harmonic generation with wavefront control, lithium-niobate-on-insulator quasi-BIC arrays for electro-optical modulation in both the linear and nonlinear regime, plasmonic–lithium niobate lattices for continuous tuning of second-harmonic chirality, and topologically asymmetric all-dielectric metasurfaces supporting quasi-trapped modes for harmonic generation and Pancharatnam–Berry phase manipulation [2405.17726], [2602.01895], [2412.03422], [2509.12745], [2507.06180].

## 1. Platform definitions and representative implementations

A recurring feature of these systems is that the metasurface periodicity, symmetry breaking, or lattice anisotropy couples otherwise dark or weakly radiative slab-guided modes to free space. The resulting resonances are described as guided-mode resonances, quasi-BICs, q-BICs, or quasi-trapped modes, depending on the structure and symmetry class. In nonlinear operation, these resonances enhance the local field overlap with $\chi^{(2)}$ or $\chi^{(3)}$ nonlinearities and thereby increase frequency-conversion efficiency or enable nonlinear wavefront engineering [2405.17726], [2602.01895], [2412.03422], [2509.12745], [2507.06180].

| Platform | Resonant structure | Demonstrated function |
|---|---|---|
| Etch-free LiNbO$_3$ metasurface | x-cut LiNbO$_3$ film, $t=307$ nm, with 1D SiO$_2$ ridges of period $a=879$ nm, width $d=500$ nm, height $h=200$ nm | SWIR-to-visible SFG imaging and edge detection |
| Silicon NPGM | 260 nm silicon slab on glass with elliptical holes, period $P=750$ nm, q-BIC at $\lambda\approx1530$ nm | THG, polarization-dependent deflection, dual-beam modulation |
| LiNbO$_3$ quasi-BIC nanowire array | x-cut LNOI, period $p=800$ nm, asymmetric 1D array of nanowires | GHz electro-optic modulation and SHG modulation |
| Plasmonic–LiNbO$_3$ anisotropic lattice | Gold nanodisks on x-cut LN, $A_y=463.5$ nm and $A_z=436.5$ nm | Continuous SH chirality tuning |
| a-Si QTM metasurface | a-Si cylinders, $D=610$ nm, $H=100$ nm, off-center hole, square lattice | THG enhancement and resonant PB-phase control |

The material set is correspondingly diverse. Lithium niobate appears in both purely dielectric and hybrid plasmonic–dielectric realizations because it combines Pockels tunability, a wide transparency window, and strong second-order nonlinearity. Silicon and amorphous silicon appear in third-order implementations because the q-BIC or QTM field enhancement boosts the cubic nonlinear polarization. A plausible implication is that the phrase “nonlinear nonlocal metasurface” now denotes a design class rather than a single geometry.

## 2. Resonant physics, nonlocality, and quality factor

In guided-mode-resonant implementations, the resonance is set by momentum matching between free-space illumination and a slab mode. For the LiNbO$_3$ up-conversion metasurface, the resonance condition is
$$
k_y + G_m = \beta(\omega),
$$
with $G_m=2\pi m/a$ and $\beta(\omega)$ the guided-mode propagation constant. Near normal incidence, the resonance disperses as
$$
\omega_{\mathrm{res}}(k_y)\approx \omega_0 \pm v_g k_y,
$$
so even a small in-plane wavevector shifts the resonance frequency. Because the guided mode propagates in-plane, an oblique incidence of even $\pm 1^\circ$ splits the resonance into a bright and a dark branch, which is the specific nonlocality emphasized in that system [2405.17726].

Temporal coupled-mode theory is used repeatedly to describe linewidth and loading. In the LiNbO$_3$ imaging platform, the transmission is written as
$$
t(\omega,k_y)= t_0 \cdot \frac{\omega-\omega_{\mathrm{res}}(k_y)+i(\gamma_e+\gamma_i)}{\omega-\omega_{\mathrm{res}}(k_y)+i\gamma_i},
$$
with
$$
Q=\frac{\omega_{\mathrm{res}}}{2(\gamma_e+\gamma_i)}.
$$
Measured values were $Q\approx 40$–$200$ depending on fabrication, and the highest-Q samples reached approximately $360$ [2405.17726]. In the lithium-niobate electro-optic modulator, the corresponding quasi-BIC linewidth was below $0.2$ nm, with $Q_{\rm sim}\approx11000$ and $Q_{\rm exp}\gtrsim8000$ around $\lambda_0\approx1550$ nm [2412.03422].

The symmetry mechanism differs across platforms but serves the same role. In the LNOI modulator, a folded guided mode is symmetry-protected at perfect mirror symmetry and becomes a leaky quasi-BIC once the unit cell is made asymmetric. In the a-Si harmonic-generation platform, an off-center cylindrical hole breaks the native $C_\infty$ symmetry to $C_1$, enabling a dark trapped mode to couple weakly to free space and become a high-Q quasi-trapped mode. There, temporal coupled-mode theory gives
$$
\frac{da}{dt}=(j\omega_0-\gamma_r-\gamma_a)a+\kappa s_+,\qquad
Q=\frac{\omega_0}{2(\gamma_r+\gamma_a)},
$$
and experiment reported $Q$ as high as $500$ for $r_h=50$ nm and $A=50$ nm [2507.06180].

The chirality-control platform realizes nonlocality differently, through two orthogonally propagating guided-mode resonances enabled by lattice anisotropy and LN birefringence. Under normal incidence, momentum matching excites a TE mode along $y$ at $\lambda_y\simeq833.5$ nm and a TE mode along $z$ at $\lambda_z\simeq830.0$ nm, with simulated quality factors $Q_y\simeq54$ and $Q_z\simeq244$ [2509.12745]. This two-resonance structure is central to continuous polarization control at the harmonic frequency.

## 3. Nonlinear frequency conversion mechanisms

The nonlinear processes demonstrated in this class span sum-frequency generation, second-harmonic generation, and third-harmonic generation. In the LiNbO$_3$ infrared-imaging platform, the mechanism is SFG with
$$
\omega_3=\omega_1+\omega_2,
$$
where the pump is near $860$ nm, the SWIR signal is resonant at $1528$ nm, and the visible output is near $550$ nm. The nonlinear polarization is written as
$$
P^{(2)}(r,\omega_3)=\epsilon_0\chi^{(2)}(r)\cdot E(r,\omega_1)\cdot E(r,\omega_2),
$$
and, in the undepleted-pump thin-film limit,
$$
E_3(x,y)\propto \epsilon_0\chi^{(2)} E_1(x,y) * E_2(x,y).
$$
Because the film is sub-wavelength in $z$, longitudinal phase mismatch is negligible, while the grating supplies the in-plane momentum required for coupling [2405.17726].

The same work gives a standard approximate scaling for resonant SFG in a triply resonant cavity,
$$
\eta \equiv \frac{P_3}{P_1 P_2}\propto
\frac{\omega_3^2 \epsilon_0^2 |d_{\rm eff}|^2}{n_1^2 n_2^2 n_3 V_{\rm eff}^2}
(Q_1Q_2Q_3),
$$
with $d_{\rm eff}\sim \chi^{(2)}$. Experimentally, the metasurface with $Q\approx40$ reached a normalized efficiency of $1.93\times10^{-5}$ cm$^2$/GW and produced a $458\times$ stronger SFG signal than a bare LiNbO$_3$ film [2405.17726].

In the LNOI quasi-BIC modulator, the second-harmonic source is
$$
P_i^{(2)}(2\omega)=\epsilon_0\sum_{jk}\chi^{(2)}_{ijk}E_j(\omega)E_k(\omega),
$$
with radiated intensity scaling as
$$
I_{2\omega}\propto |\chi_{\rm eff}^{(2)}|^2 |E_\omega|^4.
$$
By electrically shifting the resonance and therefore the local fundamental field, the device produced up to a factor-5 change in SHG on resonance and over one order of magnitude on the resonance slopes, corresponding to $\eta_{\rm SHG}\approx0.12$ V$^{-1}$ for $\Delta V=9$ V [2412.03422].

Third-order realizations use analogous field-enhancement logic. In the silicon NPGM, the resonantly enhanced cubic polarization is
$$
P^{3\omega}(r)=\epsilon_0\chi^{(3)}[E^\omega(r)\cdot E^\omega(r)]E^\omega(r),
$$
and the measured THG conversion efficiency reached $1.45\times10^{-4}$ at $\lambda_{\rm TH}\approx511$ nm under a pump intensity of $1$ GW/cm$^2$; a log–log slope of $3$ confirmed third-order scaling [2602.01895]. In the all-dielectric a-Si QTM metasurface, the third-order polarization is written
$$
P^{(3)}(3\omega;r)=\epsilon_0\chi^{(3)}:E(\omega;r)E(\omega;r)E(\omega;r),
$$
with simulated field enhancement $|E_{\rm loc}|/|E_0|\sim20$ and an experimental THG yield enhancement of approximately $10^3$ relative to an unstructured $100$ nm a-Si film [2507.06180].

## 4. Wavefront engineering, geometric phase, and nonlinear chirality

One of the central technical themes in this area is the attempt to combine the efficiency of nonlocal resonances with the phase control usually associated with local metasurfaces. The silicon NPGM implements this through a nonlocal nonlinear geometric phase. Rotating each elliptical hole by an angle $\theta$ gives, for right-circularly polarized pumping, nonlinear polarizations with phase factors $e^{i4\theta}$ and $e^{i8\theta}$ for the co-polarized and cross-polarized THG components, respectively. In the notation of that work, the nonlinear phase is therefore $\phi_{\rm NL}=4\sigma\theta$ for the co-polarized component and $\phi_{\rm NL}=8\sigma\theta$ for the cross-polarized component. Using 12 supercells with $\theta$ varying in $15^\circ$ steps over a full $360^\circ$ cycle, the TH light is directed into discrete diffraction orders: under RCP pumping, the co-polarized THG appears in the $+2$nd order and the cross-polarized THG in the $+4$th order, while LCP reverses the sign to the negative orders [2602.01895].

The all-dielectric a-Si QTM platform uses a related but distinct PB-phase framework. For an $n$th-order nonlinear process under circularly polarized pump of helicity $\sigma=\pm1$, a rotated meta-atom imparts
$$
\phi_{\rm co}^{(n)}=\sigma(n-1)\theta,\qquad
\phi_{\rm cross}^{(n)}=\sigma(n+1)\theta.
$$
For THG this becomes
$$
\phi_{\rm co}^{(3)}=2\sigma\theta,\qquad
\phi_{\rm cross}^{(3)}=4\sigma\theta.
$$
By rotating $\theta$ from $0$ to $\pi$, the co-polarized TH acquires full $2\pi$ phase coverage and the cross-polarized TH full $4\pi$ coverage. A specific result of that work is that a slight boundary perturbation switches the PB response on only at resonance: off resonance the phase is approximately zero, whereas at $\omega\approx\omega_0$ the local hole region dominates and the nonlinear geometric phase follows the rotation angle [2507.06180].

Continuous control of nonlinear polarization is demonstrated most explicitly in the plasmonic–LiNbO$_3$ chirality platform. There, the orthogonally propagating guided-mode resonances are described by steady-state amplitudes
$$
a(\theta)\propto \frac{\kappa_y E_0\cos\theta}{\Delta\omega_y+i\gamma_y},\qquad
b(\theta)\propto \frac{\kappa_z E_0\sin\theta}{\Delta\omega_z+i\gamma_z},
$$
which feed the tensorial $\chi^{(2)}$ response of x-cut LN. The emitted SH chirality is quantified by
$$
S_3=\frac{I({\rm RCP})-I({\rm LCP})}{I({\rm RCP})+I({\rm LCP})}
=\frac{2\,{\rm Im}\{E_z(2\omega)E_y^*(2\omega)\}}{|E_z(2\omega)|^2+|E_y(2\omega)|^2}.
$$
Under the analytical model,
$$
S_3(\theta)=\frac{2R(\theta)\sin\phi(\theta)}{1+R(\theta)^2}.
$$
At $\lambda_{\rm pump}=831$ nm in simulation, $S_3$ varies from $+0.991$ at $\theta=58^\circ$ to $-0.993$ at $\theta=133^\circ$ while the SH intensity remains within $\pm10\%$ of its mean. Experimentally, at $\lambda_{\rm pump}=828$ nm, $S_3$ was tuned from $+0.920$ at $\theta=60^\circ$ to $-0.815$ at $\theta=120^\circ$, with degree of polarization approximately $0.96$ and intensity variation below $20\%$ [2509.12745].

## 5. Imaging, electro-optical control, and analog processing

Nonlinear nonlocal metasurfaces are not limited to isolated frequency-conversion experiments; they also operate as image transducers and electrically tunable optical elements. In the SWIR up-conversion imaging system, the object is illuminated by the SWIR beam, lens $L_1$ Fourier-images the object onto the metasurface so that the transfer function is uniform in momentum space, and $L_2$ decodes the SFG back to real space on a standard silicon CMOS camera. This architecture was chosen because strong angular dispersion would otherwise suppress or distort spatial frequencies. Experimentally, images with high conversion efficiency and resolution quality were obtained despite strong nonlocality; the spatial resolution was approximately $50\,\mu$m on a Siemens star, and the limit was attributed to the 4-lens imaging chain rather than to the metasurface itself. The signal-to-noise ratio was reported to be on par with direct InGaAs detection and significantly better than a thin-film up-converter [2405.17726].

The same platform also demonstrates simultaneous imaging and edge detection. A spatially variant nonlocal 1D metagrating with a topological phase dislocation produces direct wavefront-preserving up-conversion in the zeroth SFG order,
$$
E_3^{(0)}(x,y)\propto E_2(x,y),
$$
while the first diffraction orders implement a first-derivative filter,
$$
E_3^{(\pm1)}(x,y)\propto \partial_x E_2(x,y),
$$
up to convolution with the pump. Numerically, the direct image is carried by $m=0$ and high-contrast edges by $m=\pm1$ [2405.17726].

Electrical control is realized in the LNOI nonlocal metasurface through the Pockels effect. The first-order change of the inverse-index tensor is
$$
\Delta(n^{-2})_{ij}=r_{ijk}E_k,
$$
or, in scalar form for a principal direction,
$$
\Delta n=-\tfrac12 n^3 r_{ij}E_j.
$$
For x-cut LiNbO$_3$, the largest coefficient is $r_{33}\approx34$ pm/V when the field $E_z$ is aligned to the extraordinary axis. The resonance shift changes the reflectance according to
$$
\Delta\lambda_0=\Bigl(\frac{\partial\lambda_0}{\partial n}\Bigr)\Delta n,\qquad
\Delta R\approx \frac{dR}{d\lambda}\Delta\lambda_0.
$$
The device yielded $\Delta R/R\approx0.1$, with peak-to-peak relative modulation approximately $0.24$ at $20$ V$_{\rm pp}$, a linear modulation efficiency $\eta_{\rm lin}\approx0.015$ V$^{-1}$ for $V_{\rm pp}\le10$ V, a 3 dB bandwidth around $800$ MHz, and detectable modulation up to $1.4$ GHz [2412.03422].

The silicon NPGM extends the control dimension from beam steering to active nonlinear intensity modulation by dual-beam interference. With two fundamental beams,
$$
E_s^\omega = E_1 e^{i\phi_1}+E_2 e^{i\phi_2},
$$
the cubic nonlinearity gives
$$
I_{3\omega}\propto |E_1+E_2 e^{i\Delta\phi}|^6.
$$
By varying relative phase, polarization, and intensity ratio, the THG efficiency was tuned from $3.9\times10^{-9}$ to $5.5\times10^{-3}$, with near-unity modulation depth [2602.01895].

## 6. Trade-offs, misconceptions, and design directions

A persistent design trade-off in the field is the relation between quality factor and usable bandwidth. In the infrared up-conversion metasurface, higher-$Q$ samples near $360$ produced only a modest additional gain of approximately $1.3$–$2$ over lower-$Q$ samples because narrow resonances filtered the femtosecond pulse spectrum. For pulse durations around $200$ fs, that work identifies an optimal range of approximately $Q\approx200$–$400$ [2405.17726]. This directly qualifies the common assumption that larger $Q$ is always better for nonlinear conversion.

A second recurrent issue is the perceived incompatibility between nonlocal efficiency and meta-atom-level phase control. Both the silicon NPGM and the a-Si QTM platform explicitly target this point. One states that nonlocal metasurfaces enable high nonlinear conversion efficiency while local ones offer versatile wavefront control, yet achieving both within a single metasurface remains challenging; the other states that existing designs suffer from a trade off between the high efficiency of nonlocal metasurfaces and the precise wavefront control enabled by local ones [2602.01895], [2507.06180]. Their demonstrations of nonlinear geometric phase, order-selective deflection, and resonant-only PB control show that the two capabilities need not remain decoupled.

Strong nonlocality is also sometimes taken to imply poor imaging fidelity. The LiNbO$_3$ up-conversion system provides a counterexample under a precise condition: Fourier-plane up-conversion is required to avoid loss of spatial frequencies. Within that architecture, direct imaging and derivative filtering coexist on the same device [2405.17726]. A related practical constraint is diffraction at the generated wavelength; in the same work, higher orders at $\lambda_3$ had to be filtered by the collection objective.

The present design directions are material, algorithmic, and system-level. Proposed materials include higher-$\chi^{(2)}$, transparent platforms such as GaP and InGaP for stronger conversion. Proposed methodologies include inverse design and adjoint optimization for tailored dispersion and polarization independence, as well as 2D nonlocal metasurfaces for full-Stokes conversion and multi-channel analog computing [2405.17726]. The chirality-control platform identifies compatibility with etched rib or ridge waveguides and large-area fabrication by nanoimprint or nano-transfer [2509.12745]. The electro-optic LNOI device points toward reconfigurable free-space optical interconnects, LiDAR, spatial light modulators, compact scanners, multiplexers, and photonic neural nets [2412.03422]. The THG phase-gradient platforms identify on-chip frequency converters, switches, modulators, logic gates, nonlinear holography, vortex-beam generation, and possible extension to SPDC or four-wave mixing for quantum nonlinear optics [2602.01895], [2507.06180].

Taken together, these results define nonlinear nonlocal metasurfaces as a research area centered on resonantly enhanced nonlinear response, strong lattice-mediated mode engineering, and increasingly sophisticated control over amplitude, phase, momentum, and polarization at the generated frequency. The current literature suggests that the field is moving from single-function enhancement devices toward multifunctional platforms that combine frequency conversion, analog processing, chirality control, and electro-optic tunability on the same ultrathin photonic interface.

Source: https://www.emergentmind.com/topics/nonlinear-nonlocal-metasurface