REVIEW 4 major objections 4 minor 33 references
Unconventional high-harmonic generation in resonant membrane metasurfaces
T0 review · 4 major / 4 minor · reviewed 2026-08-06 · deepseek-v4-flash
Pith's one-line read A resonant membrane metasurface supporting a quasi-bound state in the continuum produces high harmonics whose intensity scales with non-integer powers of the pump intensity, breaking the conventional integer-power law of high-harmonic…
desk verdict Solid experiment with a real non-integer scaling observation, but the 'break the principle' framing is overstated and the theoretical mechanism is plausible but not nailed down. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The central object is the quasi-bound state in the continuum (qBIC): a symmetry-protected dark mode of a hexagonal lattice of elliptical apertures in a 1 µm silicon membrane, made weakly radiating by breaking the circular symmetry of the apertures. It produces a Fano transmission resonance at 3.96 µm with a simulated Q factor of 580 (measured 480), and at the resonance it localises the pump field strongly inside the silicon. That local field is the mechanism that drives the unconventional scaling: it boosts higher-order nonlinear susceptibilities (up to ninth order) relative to lower-order ones, so the generated harmonic power becomes a mix of several integer power laws, yielding apparent non-integer slopes. The supporting numerical model solves the nonlinear wave equation in the frequency domain, retaining all nonlinear terms driven by the pump field, including pump depletion and cross-phase modulation among harmonics, while neglecting terms that depend only on harmonic fields.
What would settle it
Measure the third-harmonic power slope over at least two decades of pump power while tuning the pump wavelength to, say, 100 nm away from the qBIC resonance in the same metasurface. The paper predicts the slope returns to the integer value 3 off-resonance and to about 2.3 on resonance; if the non-integer slope persists far from the resonance, the qBIC field-enhancement mechanism is not the cause. Alternatively, a comparative scan of the slope across metasurfaces with different Q factors would test the claim that the deviation tracks the resonant field enhancement.
Extended reading notes
Core claim
The central claim is that a qBIC resonance in a membrane metasurface drives high-harmonic generation into a regime where the generated harmonic power no longer scales as the integer power of the pump intensity, as it does in bulk solids and unpatterned membranes. Experimentally, the third, fifth, seventh, and ninth harmonics from the metasurface follow best-fit slopes of 2.3, 3, 4.6, and 4.5, respectively, while the unpatterned membrane shows slopes matching the harmonic orders. The same metasurface pumped off-resonance returns to integer scaling, so the deviation is tied to the qBIC excitation. The authors reproduce the behaviour numerically with a generalised perturbative model that keeps nonlinear susceptibilities up to ninth order and includes pump depletion, self- and cross-phase modulation, and feedback of harmonics onto the pump; the model captures the observed slopes and conversion efficiencies. The interpretation is that the enhanced local fields at the resonance make higher-order susceptibility contributions comparable to lower-order ones, so the effective nonlinearity of the system becomes a mixture of orders.
Load-bearing premise
The numerical story depends on assumed values for the material's higher-order nonlinear susceptibilities and on neglecting nonlinear terms that involve only harmonic fields, so if those assumptions are wrong the simulated slopes could match the data for the wrong reason.
Editorial extensions
If this is right
- Harmonic order no longer fixes the intensity scaling exponent in resonant metasurfaces; the effective exponent becomes a resonance-dependent quantity.
- Resonant field enhancement lets high harmonics up to the ninth be observed at pump fluences around 2 mJ cm⁻², where an unpatterned membrane of the same material produces no detectable ninth harmonic.
- The observed enhancement grows with harmonic order, exceeding three orders of magnitude for the seventh harmonic relative to the unpatterned membrane.
- A full perturbative treatment including higher-order susceptibilities and pump feedback is sufficient to describe the non-integer scaling, so the effect does not require a non-perturbative strong-field mechanism.
- Coupling efficiency into the qBIC matters: metasurfaces with Q factors above roughly 1000 show no enhanced HHG in this experiment, meaning the resonance must be matched to the pump linewidth.
Reading between the lines
- If the effective nonlinear order is set by the local field strength, the power-law slope could become a designable parameter, tunable through resonance Q, detuning, or aperture geometry.
- The same mechanism should appear in other resonant platforms, such as dielectric resonators or epsilon-near-zero films, whenever local fields make high-order susceptibility terms compete with low-order ones; measuring the slope could serve as a generic probe of resonant field enhancement.
- A testable extension is to pump with two beams or polarization states to see whether the non-integer slopes are accompanied by additional inter-order mixing, which the feedback model would predict.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper reports high-harmonic generation (HHG) from a free-standing silicon membrane metasurface supporting a quasi-bound-state-in-the-continuum (qBIC) resonance at 3.96 µm. Experiments show that resonant pumping enhances harmonic signals—by more than three orders of magnitude for the seventh harmonic—and enables observation of the ninth harmonic only on resonance. The authors also report non-integer intensity scaling exponents for the harmonic powers (3ω: 2.3, 5ω: 3, 7ω: 4.6, 9ω: 4.5), in contrast to the integer slopes observed for an unpatterned membrane. A frequency-domain nonlinear wave equation model, including susceptibilities up to ninth order and feedback terms at the pump frequency, is used to reproduce the trends. The central claimed result is that high-Q resonances modify the effective nonlinearity, producing unconventional power scaling.
Significance. If the observation is correct, this is a significant result for nonlinear nanophotonics: it would demonstrate that resonant field enhancement can substantially modify the effective order of a nonlinear process, breaking the conventional integer-power-law scaling of HHG. The experimental platform—free-standing membrane metasurfaces—is clean, and the comparison with unpatterned membranes is a useful control. The reported enhancements and the first observation of ninth-harmonic generation from this platform are valuable. However, the theoretical explanation relies on higher-order susceptibilities whose values are assumed from prior literature or an atomic scaling model, and the finite spectral bandwidth of the pump is not included in the model. These issues leave the mechanism open to alternative explanations, so the paper requires additional validation before the central claim can be accepted.
major comments (4)
- [Section III.A, Figure 4d] The non-integer slopes are the central experimental claim, but Figure 4d shows no error bars or confidence intervals, and the text does not describe the fitting procedure (fit range, weighting, number of points, or how the 'inflection at lower powers' is handled). Without uncertainties, the difference between slopes such as 2.3 and 3 for the third harmonic cannot be assessed as statistically meaningful. Please provide error bars, fit details, and the off-resonance slopes for the patterned membrane that are mentioned in the text but not shown in a figure.
- [Methods V.D and Section III.A] The numerical model expands the field only over harmonic components and therefore treats the pump as effectively monochromatic. This is a serious limitation because the 250 fs pump has a Fourier-limited bandwidth of roughly 90 nm, while the qBIC resonance has a measured Q of 480, corresponding to a linewidth of about 8 nm. The measured harmonic signal is an integral over the pump spectrum weighted by the resonance enhancement, so any power-dependent spectral reshaping (self-phase modulation, thermal effects, free-carrier dispersion) or power-dependent resonance shift could produce apparent non-integer scaling that is unrelated to higher-order susceptibilities. The model cannot capture this effect. Please include a quantitative estimate of the spectral-overlap contribution or provide a control experiment with a spectrally resolved pump and a measurement of the transmitted pump spectrum as a function of input power.
- [Methods V.D, Figure 5b] The theoretical explanation depends on the values of χ(5), χ(7), and χ(9), which are taken from an atomic field scaling model rather than measured in this work. The paper states that reproducing the trends 'necessitated' these higher-order terms and the feedback effects, but it does not report the assumed values, their dispersion, or a sensitivity analysis. With several adjustable higher-order susceptibilities and cross-phase-modulation terms, it is not clear that the non-integer slopes are a robust prediction rather than a fit. Please provide the explicit susceptibility values used, show how the predicted slopes change when they are varied within a plausible range, and justify the truncation at ninth order.
- [Section III.A, paragraphs after Figure 4] The text states that 'we observe the same integer power scaling for the off-resonance excitation of the patterned membrane,' but no such data are presented in the main text or figures. This control is essential for ruling out a trivial explanation in which the resonance merely filters the pump spectrum or introduces a power-dependent coupling. Please show the off-resonance power dependencies for the patterned membrane, ideally alongside the on-resonance data, with the same fitting and error analysis.
minor comments (4)
- [Section III.A, Figure 4 caption] The caption text says the unpatterned membrane is pumped at '4.96 µm,' but the surrounding text and the experimental tuning range indicate this should be 3.95–3.96 µm. Please correct the typo.
- [Introduction, reference [23]] The abbreviation 'CdT' appears to be a typo for CdTe; please correct it.
- [Section III.A, Figure 3a] The description of the THG spectrum says the narrow peak 'aligns precisely with the tripled frequency of the qBIC mode,' but it would be useful to state explicitly how the spectrometer wavelength calibration and the qBIC resonance position were cross-checked, since the pump is tuned over a range much broader than the resonance linewidth.
- [Conclusions] The phrase 'for the first time' in the conclusions should be qualified in light of previous reports of non-integer or non-standard power laws in resonant metasurfaces (e.g., even-harmonic generation in Ref. [21]). Please discuss how the present odd-harmonic observation differs from those prior results.
Circularity Check
No significant circularity: the non-integer power scaling is an experimental observation, and the theoretical model uses independent literature susceptibilities rather than fitting the reported slopes.
full rationale
The paper's central claim—non-integer intensity power dependencies for qBIC metasurface HHG—is established experimentally (Figs. 3 and 4) in direct comparison with an unpatterned membrane that shows integer slopes. The theoretical section (Sec. III.B and Methods V.D) solves the nonlinear wave equation with a perturbative expansion that includes susceptibilities up to ninth order and pump-feedback terms. These susceptibilities are not fitted to the measured non-integer exponents; they are taken from prior silicon dispersion/nonlinear characterizations (refs. 26, 31–33) and an atomic field scaling model. The statement 'reproducing these trends necessitated the inclusion of nonlinear susceptibilities up to the ninth order' is a post-hoc modeling requirement, not a parameter fit to the slopes. No equation in the paper defines the predicted scaling in terms of the measured slopes, and no fitted parameter is renamed as a prediction. Although refs. 26/31–33 share some authors (Vincenti, Scalora, Kivshar), the cited results are external silicon-membrane measurements and ellipsometry, not this qBIC experiment's power slopes, so the self-citation is provenance, not load-bearing. The finite pump bandwidth and spectral-overlap alternative is a scientific validity concern, not a circularity. The derivation chain is therefore self-contained with respect to the reported observation.
Assumptions & free parameters
free parameters (1)
- Effective higher-order susceptibilities χ(5), χ(7), χ(9) =
Not measured; taken from atomic field scaling model (Refs [26,31-33])
assumptions (2)
- domain assumption Nonlinear susceptibilities up to ninth order are isotropic and follow the dispersion from an atomic field scaling approach.
- domain assumption Only nonlinear terms that depend on powers of the fundamental pump field are retained; terms depending solely on higher harmonics are neglected.
Cite this review
Pith. "Pith review of Unconventional high-harmonic generation in resonant membrane metasurfaces." pith.science (2026). https://pith.science/paper/WBPNZYHS
@misc{pith2026250620873,
author = {Pith},
title = {Pith review of: Unconventional high-harmonic generation in resonant membrane metasurfaces},
year = {2026},
howpublished = {\url{https://pith.science/paper/WBPNZYHS}},
note = {Machine review of arXiv:2506.20873}
}
read the original abstract
High-harmonic generation (HHG) in solids has rapidly emerged as a promising platform for creating compact attosecond sources and probing ultrafast electron dynamics. Resonant metasurfaces are essential for enhancement of the otherwise small harmonic generation efficiency through local field enhancement and are essential to circumvent the need of phase matching constraints. Until now, the metasurface-enhanced HHG was believed to follow the conventional integer-power scaling laws that hold for non-resonant bulk HHG. Here, we discover that highly resonant metasurfaces driven by quasi-bound states in the continuum break this principle, manifesting non-integer intensity dependencies of the generated harmonic powers. We show experimentally and theoretically that these unconventional nonlinearities arise from the high-Q resonances that generate local fields strong enough to substantially alter the contribution of higher order susceptibility tensors to the effective nonlinearities of the system. Our findings reveal how harmonic generation rooted in resonant field-driven modification of effective nonlinear susceptibilities can reshape our understanding of light-matter interaction at the nanoscale.
Figures
Reference graph
Works this paper leans on
-
[21]
P. Tonkaev, F. Lai, S. Kruk, Q. Song, M. Scalora, K. Koshelev, and Y. Kivshar, Even-order optical har- monics generated from centrosymmetric-material meta- surfaces, Phys. Rev. Res. 6, 033073 (2024)
work page 2024
-
[1]
Ghimire, A
S. Ghimire, A. D. DiChiara, E. Sistrunk, P. Agostini, L. F. DiMauro, and D. A. Reis, Observation of high-order harmonic generation in a bulk crystal, Nat. Phys. 7, 138 (2011)
2011
-
[2]
Y. Zhao, Z. Chen, C. Wang, Y. Yang, and H.-B. Sun, Efficient second- and higher-order harmonic generation from LiNbO3 metasurfaces, Nanoscale 15, 12926 (2023)
work page 2023
-
[3]
H. Liu, C. Guo, G. Vampa, J. L. Zhang, T. Sarmiento, M. Xiao, P. H. Bucksbaum, J. Vuˇ ckovi´ c, S. Fan, and D. A. Reis, Enhanced high-harmonic generation from an all-dielectric metasurface, Nat. Phys. 14, 1006 (2018)
work page 2018
-
[4]
M. R. Shcherbakov, H. Zhang, M. Tripepi, G. Sartorello, N. Talisa, A. AlShafey, Z. Fan, J. Twardowski, L. A. Krivitsky, A. I. Kuznetsov, E. Chowdhury, and G. Shvets, Generation of even and odd high harmonics in resonant metasurfaces using single and multiple ultra-intense laser pulses, Nat. Commun. 12, 1 (2021)
work page 2021
-
[5]
A. Zalogina, L. Carletti, A. Rudenko, J. V. Moloney, A. Tripathi, H.-C. Lee, I. Shadrivov, H.-G. Park, Y. Kivshar, and S. S. Kruk, High-harmonic generation from a subwavelength dielectric resonator, Sci. Adv. 9, 10.1126/sciadv.adg2655 (2023)
-
[6]
D. Boukhaoui, A. Mikhneva, S. Idlahcen, J. Houard, T. Godin, L. Guiramand, I. Blum, F. Amrani, F. G´ erˆ ome, F. Benabid, D. Gauthier, W. Boutu, H. Merdji, A. Vella, and A. Hideur, High-harmonic gen- eration in solids from a high-energy fiber laser system, APL Photonics 10, 10.1063/5.0244415 (2025)
-
[7]
A. A. Lanin, E. A. Stepanov, A. B. Fedotov, and A. M. Zheltikov, Mapping the electron band structure by intra- band high-harmonic generation in solids, Optica 4, 516 (2017)
2017
Show all 33 references
-
[8]
Tancogne-Dejean, O
N. Tancogne-Dejean, O. D. M¨ ucke, F. X. K¨ artner, and A. Rubio, Impact of the electronic band structure in high- harmonic generation spectra of solids, Physical Review Letters 118, 10.1103/physrevlett.118.087403 (2017)
2017 doi
-
[9]
Ghimire, A
S. Ghimire, A. D. DiChiara, E. Sistrunk, U. B. Szafruga, P. Agostini, L. F. DiMauro, and D. A. Reis, Redshift in the optical absorption of zno single crystals in the pres- ence of an intense midinfrared laser field, Physical Review Letters 107, 10.1103/physrevlett.107.167407 (2011)
2011 doi
-
[10]
Apostolova and B
T. Apostolova and B. Obreshkov, High harmonic gener- ation from bulk diamond driven by intense femtosecond laser pulse, Diamond and Related Materials 82, 165–172 (2018)
2018
-
[11]
Vabishchevich, P
P. Vabishchevich, P. Vabishchevich, P. Vabishchevich, and Y. Kivshar, Nonlinear photonics with metasurfaces, Photonics Res. 11, B50 (2023)
2023
-
[12]
Koshelev, S
K. Koshelev, S. Lepeshov, M. Liu, A. Bogdanov, and Y. Kivshar, Asymmetric Metasurfaces with High- Q Res- onances Governed by Bound States in the Continuum, Phys. Rev. Lett. 121, 193903 (2018)
2018
-
[13]
Carletti, K
L. Carletti, K. Koshelev, C. De Angelis, and Y. Kivshar, Giant Nonlinear Response at the Nanoscale Driven by Bound States in the Continuum, Phys. Rev. Lett. 121, 033903 (2018)
2018
-
[14]
Koshelev, Y
K. Koshelev, Y. Tang, K. Li, D.-Y. Choi, G. Li, and Y. Kivshar, Nonlinear Metasurfaces Governed by Bound States in the Continuum, ACS Photonics 6, 1639 (2019)
2019
-
[15]
Kravtsov, E
V. Kravtsov, E. Khestanova, F. A. Benimetskiy, T. Ivanova, A. K. Samusev, I. S. Sinev, D. Pidgayko, A. M. Mozharov, I. S. Mukhin, M. S. Lozhkin, Y. V. Kapitonov, A. S. Brichkin, V. D. Kulakovskii, I. A. She- lykh, A. I. Tartakovskii, P. M. Walker, M. S. Skolnick, D. N. Krizhan...
2020
-
[16]
Koshelev, S
K. Koshelev, S. Kruk, E. Melik-Gaykazyan, J.-H. Choi, A. Bogdanov, H.-G. Park, and Y. Kivshar, Subwave- length dielectric resonators for nonlinear nanophotonics, Science 367, 288 (2020)
2020
-
[17]
J. Gao, M. A. Vincenti, J. A. Frantz, A. Clabeau, X. Qiao, L. Feng, M. Scalora, and N. Litchinitser, Over- coming losses through phase locking in nonlinear quasi- bound states in the continuum metasurfaces, ACS Ap- plied Nano Materials 7, 21445 (2024)
2024
-
[18]
Carletti, S
L. Carletti, S. S. Kruk, A. A. Bogdanov, C. De Ange- lis, and Y. Kivshar, High-harmonic generation at the nanoscale boosted by bound states in the continuum, Phys. Rev. Res. 1, 023016 (2019). 8
2019
-
[19]
Zograf, K
G. Zograf, K. Koshelev, A. Zalogina, V. Korolev, R. Hollinger, D.-Y. Choi, M. Zuerch, C. Spielmann, B. Luther-Davies, D. Kartashov, S. V. Makarov, S. S. Kruk, and Y. Kivshar, High-Harmonic Generation from Resonant Dielectric Metasurfaces Empowered by Bound States in the Contin...
2022
-
[20]
Tonkaev, K
P. Tonkaev, K. Koshelev, M. A. Masharin, S. V. Makarov, S. S. Kruk, and Y. Kivshar, Observation of En- hanced Generation of a Fifth Harmonic from Halide Per- ovskite Nonlocal Metasurfaces, ACS Photonics 10, 1367 (2023)
2023
-
[22]
Korolev, A
V. Korolev, A. D. Sinelnik, M. V. Rybin, P. Lazarenko, O. M. Kushchenko, V. Glukhenkaya, S. Kozyukhin, M. Zuerch, C. Spielmann, T. Pertsch, I. Staude, and D. Kartashov, Tunable high-order harmonic generation in GeSbTe nano-films, Nanophotonics 13, 3411 (2024)
2024
-
[23]
Z. Long, H. Yang, K. Tian, L. He, R. Qin, Z.-Y. Chen, Q. J. Wang, and H. Liang, High-harmonic generation in CdTe with ultra-low pump intensity and high photon flux, Commun. Phys. 6, 1 (2023)
2023
-
[24]
Sivis, M
M. Sivis, M. Taucer, G. Vampa, K. Johnston, A. Staudte, A. Yu. Naumov, D. M. Villeneuve, C. Ropers, and P. B. Corkum, Tailored semiconductors for high-harmonic op- toelectronics, Science 357, 303 (2017)
2017
-
[25]
Y. Yang, J. Lu, A. Manjavacas, T. S. Luk, H. Liu, K. Kelley, J.-P. Maria, E. L. Runnerstrom, M. B. Sin- clair, S. Ghimire, and I. Brener, High-harmonic genera- tion from an epsilon-near-zero material, Nat. Phys. 15, 1022 (2019)
2019
-
[26]
Hallman, S
K. Hallman, S. Stengel, W. Jaffray, F. Belli, M. Ferrera, M. A. Vincenti, D. de Ceglia, Y. Kivshar, N. Akozbek, S. Mukhopadhyay, J. Trull, C. Cojocaru, and M. Scalora, High-harmonic generation from subwavelength silicon films, Nanophotonics 10.1515/nanoph-2024-0468 (2025)
2025 doi
-
[27]
W. Adi, S. Rosas, A. Beisenova, S. K. Biswas, H. Mei, D. A. Czaplewski, and F. Yesilkoy, Trapping light in air with membrane metasurfaces for vibrational strong cou- pling, Nat. Commun. 15, 1 (2024)
2024
-
[28]
Rosas, W
S. Rosas, W. Adi, A. Beisenova, S. K. Biswas, F. Ku- ruoglu, F. Kuruoglu, H. Mei, M. A. Kats, D. A. Czaplewski, Y. S. Kivshar, and F. Yesilkoy, Enhanced biochemical sensing with high-Q transmission resonances in free-standing membrane metasurfaces, Optica 12, 178 (2025)
2025
-
[29]
A. C. Overvig, S. C. Malek, M. J. Carter, S. Shrestha, and N. Yu, Selection rules for quasibound states in the continuum, Physical Review B 102, 035434 (2020)
2020
-
[30]
Shakirova, A
D. Shakirova, A. C. Valero, D. Riabov, H. Altug, A. Bog- danov, and T. Weiss, Molecular chiral response enhanced by crosstalking quasi-bound states in the continuum, arXiv preprint arXiv:2505.24563 (2025)
2025 arXiv
-
[31]
Rodr ´ ıguez-Sun´ e, J
L. Rodr ´ ıguez-Sun´ e, J. Trull, N. Akozbek, D. De Ceglia, M. Vincenti, M. Scalora, and C. Cojocaru, Retrieving linear and nonlinear optical dispersions of matter: com- bined experiment-numerical ellipsometry in silicon, gold and indium tin oxide, Frontiers in Photonics 2, 74...
2021
-
[32]
Boyd, Nonlinear optics, academic, San Diego, Calif 19922, 39 (2008)
R. Boyd, Nonlinear optics, academic, San Diego, Calif 19922, 39 (2008)
2008
-
[33]
Hallman, L
K. Hallman, L. Rodr ´ ıguez-Sun´ e, J. Trull, C. Cojocaru, M. A. Vincenti, N. Akozbek, R. Vilaseca, and M. Scalora, Harmonic generation from silicon membranes at visi- ble and ultraviolet wavelengths, Optics express 31, 792 (2023)
2023
Reviewed August 6, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.