{"id":106336,"date":"2025-07-31T01:53:15","date_gmt":"2025-07-31T01:53:15","guid":{"rendered":"https:\/\/www.europesays.com\/us\/106336\/"},"modified":"2025-07-31T01:53:15","modified_gmt":"2025-07-31T01:53:15","slug":"efficient-ghz-electro-optical-modulation-with-a-nonlocal-lithium-niobate-metasurface-in-the-linear-and-nonlinear-regime","status":"publish","type":"post","link":"https:\/\/www.europesays.com\/us\/106336\/","title":{"rendered":"Efficient GHz electro-optical modulation with a nonlocal lithium niobate metasurface in the linear and nonlinear regime"},"content":{"rendered":"<p>Design and linear characterization<\/p>\n<p>Figure\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-62072-7#Fig1\" rel=\"nofollow noopener\" target=\"_blank\">1<\/a> illustrates the design and working principle of our device. We investigate a nonlocal lithium niobate on insulator (LNOI) metasurface made by asymmetric, periodically arranged nanostripes obtained from an x-cut LiNbO3 thin film on a finite SiO2 layer on top of a Si substrate (Fig.\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-62072-7#Fig1\" rel=\"nofollow noopener\" target=\"_blank\">1a<\/a>). We aim at observing efficient and fast electric modulation of the metasurface optical properties induced by a control voltage, VEO applied to the metasurface. To this aim, in-plane, gold contacts generate an electric field (<b>E<\/b>EO) aligned along the z-axis to exploit the largest component of the electro-optic tensor r33 of LiNbO3. To obtain a significant intensity modulation, the structure must sustain resonant modes, which are spectrally shifted when the refractive index is modulated by an applied voltage (Fig.\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-62072-7#Fig1\" rel=\"nofollow noopener\" target=\"_blank\">1b<\/a>). Furthermore, at a given wavelength \u03bb0, a device with large quality factor Q\u2009=\u2009\u03bb0\/\u03b4\u03bb0 features a resonance with a narrow spectral width \u03b4\u03bb0 and, consequently, a large sensitivity to small variations of the refractive index induced by the applied voltage.<\/p>\n<p><b id=\"Fig1\" class=\"c-article-section__figure-caption\" data-test=\"figure-caption-text\">Fig. 1: Linear and nonlinear electro-optic modulation with a lithium niobate metasurface.<\/b><a class=\"c-article-section__figure-link\" data-test=\"img-link\" data-track=\"click\" data-track-label=\"image\" data-track-action=\"view figure\" href=\"https:\/\/www.nature.com\/articles\/s41467-025-62072-7\/figures\/1\" rel=\"nofollow noopener\" target=\"_blank\"><img decoding=\"async\" aria-describedby=\"Fig1\" src=\"https:\/\/www.europesays.com\/us\/wp-content\/uploads\/2025\/07\/41467_2025_62072_Fig1_HTML.png\" alt=\"figure 1\" loading=\"lazy\" width=\"685\" height=\"513\"\/><\/a><\/p>\n<p><b>a<\/b> Sketch of modulated reflected power of a beam at fundamental frequency \\({P}_{\\omega }\\) and of its second harmonic \\({P}_{2\\omega }\\) from a x-cut LiNbO3 on a SiO2-Si substrate. The modulation is induced by an applied voltage (VEO), which produces a driving electric field <b>E<\/b>EO parallel to the optical field <b>E<\/b>\u03c9. <b>b<\/b> Illustration of a resonance shift induced by electro-optic effect.<\/p>\n<p>To this aim, we exploit high Q resonances hosted by a dielectric waveguide slab \u2212 namely, a high refractive index layer (LiNbO3), comprised between two low-index claddings (air and SiO2)\u2014patterned with periodically arranged nanowires (i.e. a diffraction grating)<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 34\" title=\"Huang, L. et al. Ultrahigh-Q guided mode resonances in an all-dielectric metasurface. Nat. Commun. 14, &#010;                  https:\/\/doi.org\/10.1038\/s41467-023-39227-5&#010;                  &#010;                 (2023).\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-62072-7#ref-CR34\" id=\"ref-link-section-d50345937e1093\" rel=\"nofollow noopener\" target=\"_blank\">34<\/a>. Figure\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-62072-7#Fig2\" rel=\"nofollow noopener\" target=\"_blank\">2a<\/a> shows a sketch of the elementary cell of the metasurface with the relevant geometrical parameters. A planar dielectric waveguide is known to host guided modes whose dispersion lies below the light line, fully confining the radiation in the slab volume with no access from external radiation sources. The reduction of translational symmetry introduced by the grating causes a folding of the guided modes dispersion in the first Brillouin zone, ending up above the light line, and consequently evolving in optically accessible guided mode resonances (GMRs)<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 34\" title=\"Huang, L. et al. Ultrahigh-Q guided mode resonances in an all-dielectric metasurface. Nat. Commun. 14, &#010;                  https:\/\/doi.org\/10.1038\/s41467-023-39227-5&#010;                  &#010;                 (2023).\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-62072-7#ref-CR34\" id=\"ref-link-section-d50345937e1101\" rel=\"nofollow noopener\" target=\"_blank\">34<\/a>,<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 35\" title=\"Sun, K. et al. Infinite-Q guided modes radiate in the continuum. Phys. Rev. B. 107, 115415 (2023).\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-62072-7#ref-CR35\" id=\"ref-link-section-d50345937e1104\" rel=\"nofollow noopener\" target=\"_blank\">35<\/a>. A mirror symmetry of the grating with respect to the xz plane yields a dark fundamental mode at the \u0393 point (i.e., at normal incidence) due to the odd electric field distribution (see Fig.\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-62072-7#Fig2\" rel=\"nofollow noopener\" target=\"_blank\">2b<\/a>), namely becoming a symmetry-protected BIC. By introducing a small asymmetry in the nanowires\u2019 widths, we further reduce the translational symmetry of the system, turning the BIC into a leaky, quasi-BIC with finite Q (see Fig.\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-62072-7#Fig2\" rel=\"nofollow noopener\" target=\"_blank\">2c<\/a>), which can now be excited at the \u0393 point.<\/p>\n<p><b id=\"Fig2\" class=\"c-article-section__figure-caption\" data-test=\"figure-caption-text\">Fig. 2: Sample layout and linear optical characterization.<\/b><a class=\"c-article-section__figure-link\" data-test=\"img-link\" data-track=\"click\" data-track-label=\"image\" data-track-action=\"view figure\" href=\"https:\/\/www.nature.com\/articles\/s41467-025-62072-7\/figures\/2\" rel=\"nofollow noopener\" target=\"_blank\"><img decoding=\"async\" aria-describedby=\"Fig2\" src=\"https:\/\/www.europesays.com\/us\/wp-content\/uploads\/2025\/07\/41467_2025_62072_Fig2_HTML.png\" alt=\"figure 2\" loading=\"lazy\" width=\"685\" height=\"750\"\/><\/a><\/p>\n<p><b>a<\/b> Schematic cross-section of a single unit cell of the metasurface. All the samples have the same periodicity p\u2009=\u2009800\u2009nm and LiNbO3 height before (hLN\u2009=\u2009600\u2009nm) and after patterning (h\u2019LN\u2009=\u2009450\u2009nm) with a substrate thickness, hsub\u2009=\u20092\u2009\u03bcm. The angle \u03b8\u2009=\u200965\u00b0 accounts for the sidewall slopes introduced by the physical etching. The variable parameters from sample to sample are the filling factor \\(\\left.{FF}=1-({D}_{1}\\,+{D}_{2})\/p\\right)\\) and the asymmetry parameter \\(\\alpha=({D}_{2}-{D}_{1})\/({D}_{2}+{D}_{1})\\). Red and blue arrows indicate the polarization of the incident wave that excites the TM (parallel to the LN ordinary axis) and TE (parallel to the LN extraordinary axis) mode, respectively. <b>b<\/b> COMSOL simulation of the optical field z-component of the optical field <b>E<\/b>opt, excited at the quasi-BIC resonance for FF\u2009=\u20090.3 and \u03b1\u2009=\u20090. Incident optical field <b>E<\/b>opt is TE-polarized and propagates with <b>k<\/b> || (\u2013x). <b>c<\/b> COMSOL plane-wave simulation of the metasurface reflectance depending on the asymmetry parameter \\(\\alpha\\), with fixed FF\u2009=\u20090.3. <b>d<\/b> Dark field and bright-field images in real colours of the LiNbO3 chip. The LiNbO3 metasurfaces are located in the areas with strong scattering (bright areas in dark field). The white scale bar on the bottom right is equal to 100 \u03bcm. The orange lines in the bright field image identify the edges of the interdigitated electrodes. <b>e<\/b> Scanning-electron micrograph showing the metasurface obtained by patterning the LiNbO3 thin film. The scale bar on the bottom right is equal to 500\u2009nm. The white, dashed line indicates the LiNbO3 extraordinary axis. <b>f<\/b> Normalized reflection spectra of the sample with FF\u2009=\u20090.3 and \u03b1\u2009=\u20090.19. The solid lines show experimental data obtained by sweeping the laser wavelength in steps of 10 p.m., with electric field exciting the TE (blue solid line) and TM (red solid line) mode. The dashed lines are the corresponding numerical simulations. <b>g<\/b>, <b>h<\/b> Zoom on spectra shown in panel f (thin lines) fitted with Fano profiles (thick lines) for TM and TE polarization, respectively.<\/p>\n<p>The design optimization has been performed initially by means of rigorous coupled-wave analysis, typically applied to solve scattering from periodic dielectric structures<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 36\" title=\"Moharam, M. G. &amp; Gaylord, T. K. Rigorous coupled-wave analysis of planar-grating diffraction. J. Opt. Soc. Am. 71, 811 (1981).\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-62072-7#ref-CR36\" id=\"ref-link-section-d50345937e1412\" rel=\"nofollow noopener\" target=\"_blank\">36<\/a>. The EO effect and the SHG are then modelled with a finite-element commercial software (COMSOL Multiphysics\u00ae, see Methods and Supplementary Information\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"supplementary material anchor\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-62072-7#MOESM1\" rel=\"nofollow noopener\" target=\"_blank\">S2<\/a>). By setting the unit cell periodicity as p\u2009=\u2009800\u2009nm (see Fig.\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-62072-7#Fig2\" rel=\"nofollow noopener\" target=\"_blank\">2a<\/a>), the metasurface features the fundamental TE21 resonant mode at the \u0393 point in the communication C band (i.e, 1530 to 1565\u2009nm). This mode shows up in the reflectance spectrum when the structure asymmetry \\(\\alpha=\\,({D}_{2}-\\,{D}_{1})\/\\,({D}_{2}\\,+\\,{D}_{1})\\) is nonzero, with a linewidth broadening with increasing \u03b1 (see Fig.\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-62072-7#Fig2\" rel=\"nofollow noopener\" target=\"_blank\">2c<\/a>)<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 37\" title=\"Koshelev, K. et al. Nonlinear metasurfaces governed by bound states in the continuum. ACS Photonics 6, 1639&#x2013;1644 (2019).\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-62072-7#ref-CR37\" id=\"ref-link-section-d50345937e1542\" rel=\"nofollow noopener\" target=\"_blank\">37<\/a>. We set \\({D}_{1}\\) and \\({D}_{2}\\) to have a fill factor \\({FF}=1-({D}_{1}\\,+{D}_{2})\/p=0.3\\), corresponding to a fundamental TE mode around 1550\u2009nm, and \u03b1\u2009=\u20090.2, to realize a theoretical Q\u2009=\u200911000. Such a value for the Q factor is selected since it allows more tolerance in terms of fabrication accuracy and better resonance coupling in the experiment. Nevertheless, Q is sufficiently high to achieve efficient EO modulation, since the associated mode linewidth (\u03b4\u03bb0 \\(\\simeq\\) 0.15\u2009nm) is comparable with the expected EO shift induced by CMOS-compatible voltages, assuming a tuning sensitivity \u0394\u03bbEO\/\u0394VEO \\(\\approx\\) 0.01\u2009nm\/V<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 28\" title=\"Weiss, A. et al. Tunable metasurface using thin-film lithium niobate in the telecom regime. ACS Photonics 9, 605&#x2013;612 (2022).\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-62072-7#ref-CR28\" id=\"ref-link-section-d50345937e1744\" rel=\"nofollow noopener\" target=\"_blank\">28<\/a>,<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 30\" title=\"Gao, B., Ren, M., Wu, W., Cai, W. &amp; Xu, J. Electro-optic lithium niobate metasurfaces. Sci. Chin. Phys. Mech. Astron. 64, 240362 (2021).\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-62072-7#ref-CR30\" id=\"ref-link-section-d50345937e1747\" rel=\"nofollow noopener\" target=\"_blank\">30<\/a>.<\/p>\n<p>The metasurfaces were realized on a commercial LNOI stack with 600\u2009nm-thick LiNbO3 film on top of a SiO2 layer with thickness hsub\u2009=\u20092\u2009\u03bcm (see Fig.\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-62072-7#Fig2\" rel=\"nofollow noopener\" target=\"_blank\">2d, e<\/a>) following a top-down approach (see Methods and ref. <a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 38\" title=\"Kaufmann, F., Finco, G., Maeder, A. &amp; Grange, R. Redeposition-free inductively-coupled plasma etching of lithium niobate for integrated photonics. Nanophotonics 12, 1601&#x2013;1611 (2023).\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-62072-7#ref-CR38\" id=\"ref-link-section-d50345937e1765\" rel=\"nofollow noopener\" target=\"_blank\">38<\/a>). The choice of planar Au electrodes ensures a uniform distribution of the driving electric field in the LiNbO3 film (see Fig.\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"supplementary material anchor\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-62072-7#MOESM1\" rel=\"nofollow noopener\" target=\"_blank\">S3<\/a>, Supplementary Information), which increases the overlap with the field of the optical mode, also located mostly in the LiNbO3 volume (see Fig.\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-62072-7#Fig2\" rel=\"nofollow noopener\" target=\"_blank\">2b<\/a>). The uniform electric field distribution is advantageous for electro-optic (EO) modulation, as effectively leveraging the Pockels effect depends on the overlap between the low-frequency electric field and the optical field within the LiNbO3 material. The elongated footprint of the metasurface ensures a large active area to fully exploit the high Q of the nonlocal mode, which extends in the direction where the translational symmetry breaking occurs. On the other hand, the metasurface width is reduced to 10\u2009\u03bcm to increase the driving field amplitude. Concurrently, the large electrode width minimizes the contact resistance, hence reducing the impact on the device bandwidth. We realized a set of replicas with slight variation of D1 and D2 (see Methods for details and Supplementary Information\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"supplementary material anchor\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-62072-7#MOESM1\" rel=\"nofollow noopener\" target=\"_blank\">S3<\/a>) to target the sample geometry closer to the desired quality factor and resonance position, finally identifying a fabricated replica with \u03b1\u2009=\u20090.19 and FF\u2009=\u20090.3 as the one matching the original design specifications. Additionally, to demonstrate the generality of this approach, we also investigated a sample with FF\u2009=\u20090.5, which is characterized by a fundamental mode around 1560\u2009nm (see Supplementary Information, sections S<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"supplementary material anchor\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-62072-7#MOESM1\" rel=\"nofollow noopener\" target=\"_blank\">3<\/a> and S<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"supplementary material anchor\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-62072-7#MOESM1\" rel=\"nofollow noopener\" target=\"_blank\">6<\/a>).<\/p>\n<p>The linear properties and the EO response of the sample are characterized in an epi-reflection configuration by impinging with a linearly polarized, continuous-wave (CW) diode laser tuneable in the optical communication C band (1460\u2009nm\u20131570\u2009nm) and detecting the reflected signal with a photodiode (details in Methods and Supplementary Information\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"supplementary material anchor\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-62072-7#MOESM1\" rel=\"nofollow noopener\" target=\"_blank\">S1<\/a>). The quasi-BIC resonance wavelength shifts when interrogated away from the \u0393 point of the photonic band structure, typically by about 10\u2009nm\/deg. Therefore, to mitigate the resonance linewidth broadening resulting from the angular spread of the light from the objective, we use a cylindrical lens to focus the excitation beam onto the objective back-aperture. This approach effectively reduces the numerical aperture, enabling nearly collimated illumination in the direction of the nonlocal mode extension. At the same time, due to the one-dimensional character of the metasurface, we tightly focus the beam on the direction parallel to the nanowires, hence preserving a sizeable beam fluence (see Supplementary Information, Fig.\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"supplementary material anchor\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-62072-7#MOESM1\" rel=\"nofollow noopener\" target=\"_blank\">S1a<\/a>). We determine the resonance position by sweeping the laser wavelength, obtaining the reflectance spectra depicted in Fig.\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-62072-7#Fig2\" rel=\"nofollow noopener\" target=\"_blank\">2f<\/a>. The sharp quasi-BIC peak is superimposed to a broad Fabry\u2013P\u00e9rot fringe caused by the 2\u2009\u03bcm-thick SiO2 layer. A polarization parallel or perpendicular with respect to the LiNbO3 extraordinary index ne excites either a transverse-electric (TE) or transverse-magnetic (TM) mode, respectively, resulting in a reflection peak at 1553\u2009nm or 1492\u2009nm, respectively. We also report an experimental redshift of about 5\u2009nm with respect to the COMSOL finite-element simulations of the same structures (see Fig.\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-62072-7#Fig2\" rel=\"nofollow noopener\" target=\"_blank\">2f<\/a>). This slight discrepancy could be attributed to deviations of the fabricated nanostructure from its nominal design. Due to the asymmetry in the resonance lineshape arising from the interaction between the sharp quasi-BIC and the broad Fabry-Perot fringe, we fitted the experimental data using a Fano profile (see Fig.\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-62072-7#Fig2\" rel=\"nofollow noopener\" target=\"_blank\">2g, h<\/a>)<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 39\" title=\"Limonov, M. F. Fano resonance for applications. Adv. Opt. Photonics 13, 703 (2021).\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-62072-7#ref-CR39\" id=\"ref-link-section-d50345937e1835\" rel=\"nofollow noopener\" target=\"_blank\">39<\/a>. We evaluate the resonance quality factor of the TE mode Qexp to be larger than 7500, where \\(\\delta {\\lambda }_{0}^{\\exp }\\,\\)\u22430.20\u2009nm is the resonance linewidth retrieved from the fitting. This result is in good agreement with the value calculated from the simulations (Qsim \\(\\simeq\\) 11000 with \\(\\delta {\\lambda }_{0}^{{{\\rm{sim}}}}\\)\u2009=\u20090.15\u2009nm). The lower value of Qexp with respect to Qsim might be related to fabrication tolerances and to the difference in excitation between the simulated platform (perfectly planar wavefront) and the real system (partially collimated beam).<\/p>\n<p>We investigate the static EO response of the device by applying a bias VEO\u2009=\u2009VDC by means of a waveform generator. The shift of the quasi-BIC resonance can be clearly identified in the reflectance spectra of Fig.\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-62072-7#Fig3\" rel=\"nofollow noopener\" target=\"_blank\">3a<\/a>, superimposed to an unmodulated instrumental artifact. From the central wavelength \u03bb0, retrieved by the Fano fit, we estimate an EO induced shift of about \u0394\u03bbEO \\(\\simeq\\) 0.05\u2009nm when \u0394VDC\u2009=\u20099\u2009V, close to the value retrieved from the simulated spectra in Fig.\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-62072-7#Fig3\" rel=\"nofollow noopener\" target=\"_blank\">3d<\/a> (\\(\\Delta {\\lambda }_{{{\\rm{EO}}}}^{{{\\rm{sim}}}}\\simeq\\) 0.06\u2009nm). We, therefore, estimate a device tuning sensitivity \\(\\Delta {\\lambda }_{{{\\rm{EO}}}}\/\\Delta {V}_{{{\\rm{EO}}}}\\)\u2009=\u20095.6\u2009pm\/V. To observe the dynamic modulation of the reflectance, we apply a sinusoidal driving potential \\({V}_{{{\\rm{EO}}}}\\left(t\\right)\\,=\\,\\frac{{V}_{{{\\rm{pp}}}}}{2}\\sin \\left(2\\pi {f}_{{\\mathrm{mod}}}t\\right)\\), characterized by a variable peak-to-peak amplitude Vpp and modulation frequency fmod. We use a lock-in amplifier, referenced by the same waveform generator, to demodulate at f\u2009=\u2009fmod, obtaining the signal oscillation amplitude \\({P}_{\\omega }\\left(f={f}_{{\\mathrm{mod}}}\\right).\\) The DC component of the signal \\({P}_{\\omega }\\left(f=0\\right)\\) (namely, the reflectance) is synchronously detected in order to track modulation efficiency variations while changing excitation wavelength. We modulate the signal at fmod\u2009=\u2009100\u2009kHz, well below the cutoff of our detection system (600\u2009MHz). The relative modulation is calculated by evaluating the ratio \\({P}_{\\omega }(f={f}_{{\\mathrm{mod}}})\/{P}_{\\omega }\\left(f=0\\right)\\), which corresponds to \\(\\Delta {P}_{\\omega }\/\\,{P}_{\\omega }=({P}_{\\omega }(V={V}_{{{\\rm{pp}}}}\/\\,2)-{P}_{\\omega }(V=0))\/{P}_{\\omega }(V=0)\\). The dispersive nature of the reflection spectrum will cause a dependence on \u03bb of the modulation amplitude, which is maximized by exciting at the resonance derivative extrema. The maximum modulation amplitude that is obtained from the data in Fig.\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-62072-7#Fig3\" rel=\"nofollow noopener\" target=\"_blank\">3b<\/a> is \\(\\varDelta {P}_{{\\it{\\omega }}}\/{P}_{{\\it{\\omega }}}\\cong 0.12\\) (i.e., 0.24 of peak-to-peak relative modulation), obtained with an input voltage as low as Vpp\u2009=\u200920\u2009V. This is approximately 5 times smaller than the value obtained from simulations (see Fig.\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-62072-7#Fig3\" rel=\"nofollow noopener\" target=\"_blank\">3e<\/a>). Such a discrepancy is mostly ascribed to two experimental limitations. On the one hand, the experimental reflectance background is larger than the simulated one (see Fig.\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-62072-7#Fig2\" rel=\"nofollow noopener\" target=\"_blank\">2f<\/a>), hence reducing the resonance contrast and, consequently, the relative modulation. On the other hand, the presence of fabrication tolerances (e.g., periodicity, asymmetry) and defects in the realized sample, as well as deviations from ideal plane-wave illumination, both contribute to the reduction of Qexp. Figure\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-62072-7#Fig3\" rel=\"nofollow noopener\" target=\"_blank\">3c<\/a> illustrates how electro-optically induced modulation varies with Vpp. At low modulating amplitudes \u0394\u03bbEO \u00ab \u03b4\u03bb and \u03b7 remains linear with respect to Vpp. In this regime, the modulation is proportional to the local derivative of the reflectance spectrum, and the linearity of \u0394n with respect to <b>E<\/b>EO implies the linearity of \u03b7 versus Vpp<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 27\" title=\"Weigand, H. et al. Enhanced electro-optic modulation in resonant metasurfaces of lithium niobate. ACS Photonics 8, 3004&#x2013;3009 (2021).\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-62072-7#ref-CR27\" id=\"ref-link-section-d50345937e2667\" rel=\"nofollow noopener\" target=\"_blank\">27<\/a>. In this regime, we reach a modulation efficiency\u2014defined as relative modulation per applied Volt \u2013 \\(\\eta=\\left(\\frac{\\Delta {P}_{\\omega }}{{P}_{\\omega }}\\right)\/\\left(\\frac{{V}_{{{\\rm{pp}}}}}{2}\\right)\\) of about 0.015\u2009V\u22121 for Vpp\u2009\u03bbEO is no longer negligible with respect to the resonance width, resulting in a deviation from linearity. We point out that the investigation shows a modulation amplitude of the order of 0.07 for a bias-voltage Vpp\u2009=\u200910\u2009V. Similar efficiency was reported in the literature only from devices sacrificing the modulation bandwidth to attain larger <b>E<\/b>EO<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 28\" title=\"Weiss, A. et al. Tunable metasurface using thin-film lithium niobate in the telecom regime. ACS Photonics 9, 605&#x2013;612 (2022).\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-62072-7#ref-CR28\" id=\"ref-link-section-d50345937e2791\" rel=\"nofollow noopener\" target=\"_blank\">28<\/a>. We note that the modulation efficiency of the metasurface can be further enhanced via two key strategies: (i) by increasing the quality factor of the quasi-BIC resonance by further reducing the asymmetry factor; and (ii) by suppressing the Fabry-P\u00e9rot resonance, which now limits the reflectance contrast of the TM mode, either by removing the silicon substrate or by adjusting the thickness of the SiO\u2082 layer.<\/p>\n<p><b id=\"Fig3\" class=\"c-article-section__figure-caption\" data-test=\"figure-caption-text\">Fig. 3: Characterization of the reflected signal modulation by applying DC and AC electrical stimuli.<\/b><a class=\"c-article-section__figure-link\" data-test=\"img-link\" data-track=\"click\" data-track-label=\"image\" data-track-action=\"view figure\" href=\"https:\/\/www.nature.com\/articles\/s41467-025-62072-7\/figures\/3\" rel=\"nofollow noopener\" target=\"_blank\"><img decoding=\"async\" aria-describedby=\"Fig3\" src=\"https:\/\/www.europesays.com\/us\/wp-content\/uploads\/2025\/07\/41467_2025_62072_Fig3_HTML.png\" alt=\"figure 3\" loading=\"lazy\" width=\"685\" height=\"385\"\/><\/a><\/p>\n<p><b>a<\/b> Resonance shift upon driving the metasurface with a static bias VDC. The dashed lines correspond to the reflected intensity for three different values of VDC, the solid lines are obtained by fitting the data with a Fano profile. <b>b<\/b> Relative modulation (in orange), defined as the ratio \\({P}_{\\omega }(f={f}_{{\\mathrm{mod}}})\/{P}_{\\omega }(f=0)\\), where \\({P}_{\\omega }(f={f}_{{\\mathrm{mod}}})\\) is the signal demodulated at the modulation frequency fmod and \\({P}_{\\omega }(f=0)\\) is the DC component of the signal, corresponding to the reflectance. The sinusoidal driving field is characterized by a peak-to-peak amplitude Vpp\u2009=\u200920\u2009V and a frequency fmod\u2009=\u2009100\u2009kHz. The bias VDC is set to 0. In blue, the reflectance spectrum (dashed) and its Fano profile fit (solid). The wavelength is swept in steps of 10\u2009pm. <b>c<\/b> Experimental (orange) and simulated (in blue) variation of the maximum modulation with respect to the peak-to-peak modulating amplitude Vpp. Experimental points are obtained by modulating with fmod\u2009=\u2009100\u2009kHz. <b>d<\/b> Reflectance for different applied VDC values and <b>e<\/b> reflectance spectrum and relative modulation (orange dots) simulated using COMSOL Multiphysics. The relative modulation is calculated as \\([{P}_{\\omega }({V=V}_{{pp}}\/2)-{P}_{\\omega }\\left(V=0\\right)]\/{P}_{\\omega }(V=0)\\), with Vpp\u2009=\u200920\u2009V. <b>f<\/b> Device frequency response obtained by calculating the parameter \\({S}_{21}=10{\\log }_{10}[{P}_{\\omega }(f)\/{P}_{\\omega }^{\\max }(f)]\\) with a vector-network analyzer (VNA). The VNA oscillator drives the sample with Vpp\u2009=\u20091.5\u2009V. All the panels refer to theoretical and experimental investigations on the sample with FF\u2009=\u20090.3 and \u03b1\u2009=\u20090.19. The optical field is polarized along the z axis (see Fig.\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-62072-7#Fig2\" rel=\"nofollow noopener\" target=\"_blank\">2<\/a>).<\/p>\n<p>We performed the analysis of the modulation efficiency on several fabricated designs, always obtaining values between 0.05 and 0.1 for Vpp\u2009=\u200910\u2009V (see Supplementary Information\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"supplementary material anchor\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-62072-7#MOESM1\" rel=\"nofollow noopener\" target=\"_blank\">S3<\/a>). To complete the investigation, we studied the efficiency of the TM mode, obtaining a reduced modulation depth (see Fig.\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"supplementary material anchor\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-62072-7#MOESM1\" rel=\"nofollow noopener\" target=\"_blank\">S5<\/a>, Supplementary Information) caused both by the smaller electro-optic coefficient r13\u2009=\u200910.3\u2009pm\/V<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 16\" title=\"Holmes, R. J., Kim, Y. S., Brandle, C. D. &amp; Smyth, D. M. Evaluation of crystals of LiNbo3 doped with MgO or TiO2 for electrooptic devices. Ferroelectrics 51, 41&#x2013;45 (1983).\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-62072-7#ref-CR16\" id=\"ref-link-section-d50345937e3333\" rel=\"nofollow noopener\" target=\"_blank\">16<\/a> and by the slightly larger resonance width of 0.5\u2009nm (see Fig.\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-62072-7#Fig2\" rel=\"nofollow noopener\" target=\"_blank\">2h<\/a>).<\/p>\n<p>Finally, we assessed the modulation speed of the device. To overcome the limitation set by the lock-in amplifier bandwidth (600\u2009MHz), we employed a vector-network analyzer (VNA) operating at frequencies up to 20\u2009GHz and delivering a bias peak-to-peak amplitude Vpp\u2009=\u20091.5\u2009V (see Fig.\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"supplementary material anchor\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-62072-7#MOESM1\" rel=\"nofollow noopener\" target=\"_blank\">S1b<\/a>, Supplementary Information). Figure\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-62072-7#Fig3\" rel=\"nofollow noopener\" target=\"_blank\">3f<\/a> displays the resulting VNA transmission parameter, defined as \\({S}_{21}=10{\\log }_{10}[{P}_{\\omega }(f)\/{P}_{{{\\rm{\\omega }}}}^{\\max }(f)]\\), where\\(\\,{P}_{\\omega }(f)\\) is the AC component of the modulated optical signal and \\({P}_{{{\\rm{s}}}}^{\\max }(f)\\) is its maximum value in the considered frequency range. We show a \u20133\u2009dB (\u22126\u2009dB) attenuation for a driving frequency f\u22123dB \\(\\simeq\\) 800\u2009MHz (f\u22126dB \\(\\simeq\\) 1.4\u2009GHz). Possibly, a larger bandwidth could be obtained by reducing the length of the interdigitated electrodes, which is currently 1.5\u2009mm, since their capacitance scales linearly with the length. The limit for the increased frequency will be imposed by a trade-off between the electrode length and the metasurface area, which needs a minimal size to sustain high Q optical resonances due to the nonlocal character of the mode.<\/p>\n<p>SHG modulation<\/p>\n<p>The observed modulation shown in the previous section is inherently limited by the imbalance between the high reflected power \\({P}_{\\omega }(V=0)\\) and the effective modulated power \\({\\Delta P}_{\\omega }(\\Delta V={V}_{{{\\rm{pp}}}}\/2)\\). This imposes a natural upper limit on the modulation amplitude, governed by the intrinsic resonance depth, which corresponds to the ratio between the maximum and the minimum values of the Fano profile over the background. For example, considering the measured reflectance spectrum shown in Fig.\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-62072-7#Fig3\" rel=\"nofollow noopener\" target=\"_blank\">3c<\/a>, the maximum modulation achievable by a complete shift from the resonance peak \\({P}_{\\omega }^{\\max }\\) to off-resonance \\({P}_{\\omega }^{{{\\rm{off}}}}\\) would yield \\(({P}_{\\omega }^{\\max }{-P}_{\\omega }^{{{\\rm{off}}}})\\,\/{P}_{\\omega }^{\\max }\\,\\simeq\\) 0.3.<\/p>\n<p>An advantage of modulating resonantly enhanced upconverted signals, such as SHG, is the increased sensitivity due to the quadratic dependence on the fundamental electric field, resulting in narrower linewidths. However, a significant drawback is the relatively low absolute conversion efficiency arising from the perturbative nature of the nonlinear process and the small interaction volume inherent to subwavelength devices. To generate substantial SH signals typically requires high peak intensities and, therefore, pulsed lasers. Yet, the broad bandwidths associated with pulsed sources often exceed the resonance wavelength shift \u0394\u03bbEO achievable in ultrathin devices, making efficient EO modulation of upconverted signals a challenging task. As a result, very large bias voltages are typically required, which are incompatible with CMOS integration<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 40\" title=\"He, Z. et al. Electro-optically modulated nonlinear metasurfaces. Nano Lett. 24, 14215&#x2013;14221 (2024).\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-62072-7#ref-CR40\" id=\"ref-link-section-d50345937e3895\" rel=\"nofollow noopener\" target=\"_blank\">40<\/a>. This limitation can only be circumvented by exploiting laser sources with linewidths significantly narrower than \u0394\u03bbEO, combined with high-Q resonances. Recently, Anthur et al. demonstrated that the field enhancements associated with the narrow resonances of quasi-BIC modes in a III\u2013V semiconductor metasurface allow to obtain sizeable SHG using continuous-wave (CW) pump sources, effectively compensating the reduced peak intensity<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 41\" title=\"Anthur, A. P. et al. Continuous wave second harmonic generation enabled by quasi-bound-states in the continuum on gallium phosphide metasurfaces. Nano Lett. 20, 8745&#x2013;8751 (2020).\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-62072-7#ref-CR41\" id=\"ref-link-section-d50345937e3907\" rel=\"nofollow noopener\" target=\"_blank\">41<\/a>.<\/p>\n<p>For this reason, we employed the same tuneable CW source used for the linear characterization also as pump for the SHG. We set the laser polarization along the largest component of the \u03c7(2) tensor d33, corresponding to TE illumination. The generated nonlinear signal is then spectrally filtered and detected by means of a dispersive spectrometer equipped with a charge-coupled device (CCD) (see Fig.\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"supplementary material anchor\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-62072-7#MOESM1\" rel=\"nofollow noopener\" target=\"_blank\">S1a<\/a>, Supplementary Information). The spectrum of the collected signal clearly shows the SH peak (see Fig.\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"supplementary material anchor\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-62072-7#MOESM1\" rel=\"nofollow noopener\" target=\"_blank\">S6<\/a>, Supplementary Information) without any incoherent background. We characterize the SH signal as a function of the wavelength of the pump, with a pumping intensity of 12\u2009kW\/cm2 (see Fig.\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-62072-7#Fig4\" rel=\"nofollow noopener\" target=\"_blank\">4b<\/a>).<\/p>\n<p><b id=\"Fig4\" class=\"c-article-section__figure-caption\" data-test=\"figure-caption-text\">Fig. 4: Second harmonic modulation by CD electrical stimuli.<\/b><a class=\"c-article-section__figure-link\" data-test=\"img-link\" data-track=\"click\" data-track-label=\"image\" data-track-action=\"view figure\" href=\"https:\/\/www.nature.com\/articles\/s41467-025-62072-7\/figures\/4\" rel=\"nofollow noopener\" target=\"_blank\"><img decoding=\"async\" aria-describedby=\"Fig4\" src=\"https:\/\/www.europesays.com\/us\/wp-content\/uploads\/2025\/07\/41467_2025_62072_Fig4_HTML.png\" alt=\"figure 4\" loading=\"lazy\" width=\"685\" height=\"298\"\/><\/a><\/p>\n<p><b>a<\/b> Spectrum of the reflected signal from the sample with FF\u2009=\u20090.3 and \u03b1\u2009=\u20090.19 (see also the Supplementary Information, section S<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"supplementary material anchor\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-62072-7#MOESM1\" rel=\"nofollow noopener\" target=\"_blank\">6<\/a> for the sample with FF\u2009=\u20090.5). <b>b<\/b> Normalized SHG excitation spectrum, obtained by sweeping the wavelength of the fundamental beam with a constant optical intensity on the sample of 12\u2009kW\/cm2. <b>c\u2013f<\/b> EO modulation of SHG, excited with four different wavelengths (identified by vertical dashed lines in <b>a<\/b>, <b>b<\/b> depending on the applied DC voltage VDC.<\/p>\n<p>The data reveal a full width at half maximum of the SHG excitation spectrum of \\(\\delta {\\lambda }_{{{\\rm{SHG}}}}\\,=\\,0.14\\,{{\\rm{nm}}}\\simeq \\frac{\\,\\delta {\\lambda }_{0}^{\\exp }}{\\sqrt{2}}\\). The SHG peak aligns with the maximum of the electromagnetic energy stored in LiNbO3, which notably does not correspond to the maximum of the reflectance spectrum shown in Fig.\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-62072-7#Fig4\" rel=\"nofollow noopener\" target=\"_blank\">4a<\/a> due to the asymmetry in the Fano line shape (see Fig.\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"supplementary material anchor\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-62072-7#MOESM1\" rel=\"nofollow noopener\" target=\"_blank\">S7<\/a>, Supplementary Information). To tune the SHG emission via the Pockels effect, we selected various excitation wavelengths (represented by vertical, dashed lines in Fig.\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-62072-7#Fig4\" rel=\"nofollow noopener\" target=\"_blank\">4b<\/a>, corresponding to the SHG peak, its inflection points, and to off-resonant excitation), while varying the applied static voltage VDC. In case of on-resonance pumping (as in Fig.\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-62072-7#Fig4\" rel=\"nofollow noopener\" target=\"_blank\">4e<\/a>), we report a reduction up to a factor of 5 upon application of \u0394VDC\u2009=\u2009\u00b1\u20099\u2009V. The change in intensity becomes monotonic when the excitation is tuned on the slopes of the SHG peak (Fig.\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-62072-7#Fig4\" rel=\"nofollow noopener\" target=\"_blank\">4d, f<\/a>), with \\([{P}_{{SHG}}\\left(9{V}\\right)-{P}_{{SHG}}(0{V})]\/{P}_{{SHG}}(0{V}) &gt; 1.1\\) of relative SHG intensity variation by applying \u0394VDC\u2009=\u20099\u2009V, which is in good agreement with the nonlinear simulation under same excitation conditions (a factor 2 higher, see Fig.\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"supplementary material anchor\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-62072-7#MOESM1\" rel=\"nofollow noopener\" target=\"_blank\">S8<\/a>, Supplementary Information). Notably, we reach a modulation efficiency \\(\\eta=\\,\\frac{\\Delta {P}_{{{\\rm{SHG}}}}}{{P}_{{{\\rm{SHG}}}}}\/{\\Delta V}_{{DC}} &gt; 0.12\\,{V}^{-1}\\) The opposing trends shown in Fig.\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-62072-7#Fig4\" rel=\"nofollow noopener\" target=\"_blank\">4d, f<\/a> with respect to the sign of VDC confirm the genuine electro-optic origin of the change in the SHG power (namely, VDC\u2009&gt;\u20090\u2009V leads to a blue-shift of the resonance, enhancing\/suppressing the SHG when excited on the rising\/falling slope of the SHG excitation spectrum, respectively). Off-resonant excitation, on the other hand, shows low SHG signal intensity and no electro-optic modulation thereof (Fig.\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-62072-7#Fig4\" rel=\"nofollow noopener\" target=\"_blank\">4c<\/a>). Importantly, the low SHG signal outside the resonance leads to a high modulation efficiency, close to the simulated value. This is in striking contrast with the linear case, where the modulation is limited by a large reflectance baseline. We note that the maximum modulation that could be obtained with a complete shift of the resonance is only limited by a non-zero off-resonant SHG contribution.<\/p>\n","protected":false},"excerpt":{"rendered":"Design and linear characterization Figure\u00a01 illustrates the design and working principle of our device. We investigate a nonlocal&hellip;\n","protected":false},"author":3,"featured_media":106337,"comment_status":"","ping_status":"","sticky":false,"template":"","format":"standard","meta":{"footnotes":""},"categories":[25],"tags":[10046,16250,10047,16251,16163,492,159,67,132,68],"class_list":{"0":"post-106336","1":"post","2":"type-post","3":"status-publish","4":"format-standard","5":"has-post-thumbnail","7":"category-physics","8":"tag-humanities-and-social-sciences","9":"tag-metamaterials","10":"tag-multidisciplinary","11":"tag-nanophotonics-and-plasmonics","12":"tag-nonlinear-optics","13":"tag-physics","14":"tag-science","15":"tag-united-states","16":"tag-unitedstates","17":"tag-us"},"share_on_mastodon":{"url":"https:\/\/pubeurope.com\/@us\/114945353623852613","error":""},"_links":{"self":[{"href":"https:\/\/www.europesays.com\/us\/wp-json\/wp\/v2\/posts\/106336","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/www.europesays.com\/us\/wp-json\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/www.europesays.com\/us\/wp-json\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/www.europesays.com\/us\/wp-json\/wp\/v2\/users\/3"}],"replies":[{"embeddable":true,"href":"https:\/\/www.europesays.com\/us\/wp-json\/wp\/v2\/comments?post=106336"}],"version-history":[{"count":0,"href":"https:\/\/www.europesays.com\/us\/wp-json\/wp\/v2\/posts\/106336\/revisions"}],"wp:featuredmedia":[{"embeddable":true,"href":"https:\/\/www.europesays.com\/us\/wp-json\/wp\/v2\/media\/106337"}],"wp:attachment":[{"href":"https:\/\/www.europesays.com\/us\/wp-json\/wp\/v2\/media?parent=106336"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/www.europesays.com\/us\/wp-json\/wp\/v2\/categories?post=106336"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/www.europesays.com\/us\/wp-json\/wp\/v2\/tags?post=106336"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}