Harmonic Spectral Resonance

• Harmonic Spectral Resonance is a phenomenon in nonlinear systems where discrete resonances amplify driving fields and shape harmonic output.
• It governs both amplitude and phase by exploiting interference, Fano asymmetries, and cavity enhancements to control frequency conversion.
• Applications span nanophotonics, spintronics, and plasma physics, enabling tailored optical filtering and efficient multicolor harmonic generation.

Searching arXiv for papers related to harmonic spectral resonance and resonance-enhanced harmonic generation. Harmonic spectral resonance denotes the class of phenomena in which a harmonic response is not determined solely by the order of the nonlinearity, but is strongly shaped by the spectral structure of the underlying system: discrete resonances can amplify the driving field, filter the generated harmonic, interfere with a continuum, or impose phase and timing constraints on the emitted radiation. In nonlinear nanophotonics this includes quasi-bound states in the continuum that enhance third-harmonic generation through $|E_\omega|^6$ scaling and intensity-dependent self-action; in exciton-polariton microcavities it appears as Fano-shaped second-harmonic spectra; in high-harmonic generation it can modify both amplitude and spectral phase through resonant bound states; and analogous harmonic-resonant structures also occur in spintronics, plasma emission, internal-wave dynamics, superconducting resonators, and forced oscillators with unusual spectral content (Nardi et al., 2 Jul 2026, Wang et al., 2019, Camp et al., 2015, David et al., 13 Mar 2025, Deng et al., 28 Jun 2025, Beydeda et al., 2023, Liang et al., 2016, Gendelman et al., 6 Jul 2026).

1. Core definitions and mathematical structure

In nonlinear optics, harmonic generation is frequency conversion in which an input field at frequency $\omega$ drives a nonlinear polarization that radiates at integer multiples $n\omega$. For third-harmonic generation in a centrosymmetric dielectric, the relevant bulk response is third order, with

$$P^{(3)}(3\omega)\propto \chi^{(3)}(3\omega;\omega,\omega,\omega)\,E_\omega^3,$$

so that

$$I_{3\omega}\propto |P^{(3)}(3\omega)|^2\propto |\chi^{(3)}|^2|E_\omega|^6.$$

If the fundamental field is resonantly enhanced by a mode with local-field factor $F(\omega)$, then

$$E_\omega^{\mathrm{loc}}=F(\omega)E_\omega^{\mathrm{inc}},\qquad I_{3\omega}(\omega)\propto |\chi^{(3)}|^2|F(\omega)|^6|E_\omega^{\mathrm{inc}}|^6.$$

This is one canonical form of harmonic spectral resonance: a narrow resonance strongly amplifies the fundamental and sharply shapes the harmonic yield in frequency space (Nardi et al., 2 Jul 2026).

A second canonical form is interference between a discrete resonant state and a broad harmonic continuum. In angle-resolved second-harmonic generation from ZnO microwires, the harmonic intensity near a polariton mode is described by the Fano form

$$I(\omega)=A\frac{(\epsilon+q)^2}{\epsilon^2+1}+I_b, \qquad \epsilon=\frac{2(\omega-\omega_r)}{\Gamma},$$

where $q$ is the asymmetry parameter. Here harmonic spectral resonance is not simply a peak in conversion efficiency; it is an asymmetric line shape produced by constructive and destructive interference between direct SHG and SHG coupled into an exciton-polariton mode (Wang et al., 2019).

A third form is resonance in strong-field harmonic generation. In helium, resonant enhancement of a harmonic occurs when a Stark-shifted bound state satisfies

$|E_{np}-E_0|+U_p=q\hbar\omega,$

so that a specific harmonic order is enhanced by a multiphoton resonance with a dressed excited state (Camp et al., 2015). In Brunel radiation, the harmonic spectrum takes yet another form: the radiated field is proportional to the convolution of the electron-density spectrum and the laser spectrum,

$\omega$0

so resonant frequencies emerge from the spectral structure of both the ionization bursts and the driving field (David et al., 13 Mar 2025).

These formulations indicate that harmonic spectral resonance is a family of mechanisms rather than a single line-shape archetype. Depending on the platform, the controlling object may be a cavity enhancement factor, a Fano asymmetry parameter, a dressed-state resonance condition, or a convolution between a generated density spectrum and an input field spectrum.

2. Resonantly enhanced harmonic generation in metasurfaces and dielectric nanophotonics

High-$\omega$1 dielectric metasurfaces provide a particularly direct realization of harmonic spectral resonance. In a nonlocal amorphous-Si metasurface built from asymmetric dimers, a quasi-bound state in the continuum produces a narrow Fano resonance in transmission and strong confinement of the delocalized resonant mode. Under spectrally narrow picosecond excitation, the resonance is fully resolved and third-harmonic generation exhibits a sharp wavelength-dependent peak aligned with the qBIC resonance; at high intensities that peak redshifts, broadens, and departs from cubic scaling, with the exponent in $\omega$2 dropping from $\omega$3 to $\omega$4 near resonance. Under broadband femtosecond excitation, the same resonance acts as a spectral filter on the harmonic, producing narrow features, distortions, and intensity-dependent broadening of the TH spectrum. Continuous-wave modeling reproduces the picosecond power scaling and higher-harmonic interactions, while time-domain Maxwell-Lorentz simulations reproduce the femtosecond spectral reshaping (Nardi et al., 2 Jul 2026).

A related but distinct implementation is the resonance-gradient metasurface. Here a family of germanium “bone-like” resonators supports a high-$\omega$5 Mie-like mode in the mid-IR, with optimal cutout depth around $\omega$6 nm giving $\omega$7. A spatial scaling factor $\omega$8 shifts the resonance wavelength approximately linearly, so a single $\omega$9 mm strip spans a continuous range of resonant wavelengths. By tuning the pump between $n\omega$0 and $n\omega$1 and shifting the focal spot along the strip, resonantly enhanced third- and fifth-harmonic generation is obtained across the full tuning band. For an unstructured Ge film of the same thickness under identical conditions, no detectable harmonics are observed. At lower powers the measured scaling laws are sub-perturbative, with THG power $n\omega$2 and FHG power $n\omega$3, which the paper attributes to saturation and onset of nonperturbative effects (Jangid et al., 2023).

These two metasurface architectures illustrate complementary strategies. The qBIC system concentrates on extreme field confinement and the resulting self-action of a single high-$n\omega$4 resonance. The gradient metasurface instead replicates a narrow resonant response across space, thereby circumventing the usual high-$n\omega$5/bandwidth trade-off without broadening the resonance itself. This suggests that, in metasurfaces, harmonic spectral resonance can be engineered either as a strongly nonlinear local phenomenon or as a spatially programmable spectral resource (Nardi et al., 2 Jul 2026, Jangid et al., 2023).

3. Interference, phase control, and ultrafast spectral reshaping

In ZnO microwires, second-harmonic spectral resonance takes the form of exciton-polariton Fano interference. The broad SHG band acts as a continuum, while narrow exciton-polariton whispering-gallery modes provide discrete resonances. When the SH frequency is tuned across a polariton mode, the SH spectrum changes from symmetric to strongly asymmetric depending on detection angle, and the asymmetry parameter is directly linked to the phase of the discrete polariton channel by

$n\omega$6

Angle-resolved measurements therefore probe not only spectral enhancement and suppression, but the phase structure of the mode itself (Wang et al., 2019).

In resonance-enhanced high-harmonic generation from a tin plasma, the salient effect is phase rather than amplitude alone. A strong Sn$n\omega$7 autoionizing transition near $n\omega$8 eV enhances harmonic $n\omega$9, but RABITT measurements show that the resonance also changes the relative phase of neighboring harmonics. In the detuned case, sideband $$P^{(3)}(3\omega)\propto \chi^{(3)}(3\omega;\omega,\omega,\omega)\,E_\omega^3,$$0 exhibits a group-delay shift of about $$P^{(3)}(3\omega)\propto \chi^{(3)}(3\omega;\omega,\omega,\omega)\,E_\omega^3,$$1 as experimentally and about $$P^{(3)}(3\omega)\propto \chi^{(3)}(3\omega;\omega,\omega,\omega)\,E_\omega^3,$$2 as in TDSE simulations. In the more resonant case, sidebands $$P^{(3)}(3\omega)\propto \chi^{(3)}(3\omega;\omega,\omega,\omega)\,E_\omega^3,$$3 and $$P^{(3)}(3\omega)\propto \chi^{(3)}(3\omega;\omega,\omega,\omega)\,E_\omega^3,$$4 lose clear $$P^{(3)}(3\omega)\propto \chi^{(3)}(3\omega;\omega,\omega,\omega)\,E_\omega^3,$$5 oscillations, indicating strong intra-line phase variation and loss of phase locking. The analysis attributes this behavior to the phase of the recombination dipole moment rather than to amplitude enhancement alone (Haessler et al., 2012).

A closely related conclusion appears in helium. Harmonics $$P^{(3)}(3\omega)\propto \chi^{(3)}(3\omega;\omega,\omega,\omega)\,E_\omega^3,$$6, $$P^{(3)}(3\omega)\propto \chi^{(3)}(3\omega;\omega,\omega,\omega)\,E_\omega^3,$$7, and $$P^{(3)}(3\omega)\propto \chi^{(3)}(3\omega;\omega,\omega,\omega)\,E_\omega^3,$$8 are enhanced when Stark-shifted $$P^{(3)}(3\omega)\propto \chi^{(3)}(3\omega;\omega,\omega,\omega)\,E_\omega^3,$$9 states are brought into multiphoton resonance, but time-frequency analysis shows that both short and long quantum paths are enhanced and that the long trajectory of harmonic $$I_{3\omega}\propto |P^{(3)}(3\omega)|^2\propto |\chi^{(3)}|^2|E_\omega|^6.$$0 acquires an additional phase shift of approximately $$I_{3\omega}\propto |P^{(3)}(3\omega)|^2\propto |\chi^{(3)}|^2|E_\omega|^6.$$1 on resonance. Macroscopic propagation slightly distorts the time profiles, yet the phase shift remains recognizable. This establishes that harmonic spectral resonance in strong-field HHG couples bound-state structure to quantum-path dynamics rather than merely increasing spectral intensity (Camp et al., 2015).

Two additional cases emphasize how resonance sculpts harmonic spectra through driven-state structure. In porphyrin thin films, the fifth harmonic at about $$I_{3\omega}\propto |P^{(3)}(3\omega)|^2\propto |\chi^{(3)}|^2|E_\omega|^6.$$2 eV lies close to the intense B band associated with the $$I_{3\omega}\propto |P^{(3)}(3\omega)|^2\propto |\chi^{(3)}|^2|E_\omega|^6.$$3 $$I_{3\omega}\propto |P^{(3)}(3\omega)|^2\propto |\chi^{(3)}|^2|E_\omega|^6.$$4-$$I_{3\omega}\propto |P^{(3)}(3\omega)|^2\propto |\chi^{(3)}|^2|E_\omega|^6.$$5 transition. The resonant multiphoton excitation leads to an early onset of non-perturbative behavior for the fifth harmonic, in contrast to higher off-resonant orders (Li et al., 2022). In Brunel radiation driven by two-color fields, the only knowledge of the optical-field extrema in the time domain is sufficient to reproduce the computed spectrum and justify the emergence of various resonance frequencies; in that setting the spectrum is governed by the convolution of the Fourier transforms of the electron density and laser field, and the resonant structure is set by the sequence of ionization bursts at field extrema (David et al., 13 Mar 2025).

Taken together, these results show that harmonic spectral resonance often reorganizes the phase and temporal structure of harmonic emission. This suggests that it is more accurate to regard resonance as a joint amplitude-phase condition, especially in ultrafast and strong-field regimes.

4. Manifestations beyond nonlinear optics

In spintronics, a direct analogue appears in Néel-torque antiferromagnetic resonance. For a collinear two-sublattice antiferromagnet driven by a staggered Néel spin-orbit torque, the dynamical susceptibilities at the first, second, and third harmonics are dramatically enhanced compared with conventional spin-orbit torque. For each sublattice, the NSOT case amplifies $$I_{3\omega}\propto |P^{(3)}(3\omega)|^2\propto |\chi^{(3)}|^2|E_\omega|^6.$$6 by $$I_{3\omega}\propto |P^{(3)}(3\omega)|^2\propto |\chi^{(3)}|^2|E_\omega|^6.$$7, $$I_{3\omega}\propto |P^{(3)}(3\omega)|^2\propto |\chi^{(3)}|^2|E_\omega|^6.$$8 by $$I_{3\omega}\propto |P^{(3)}(3\omega)|^2\propto |\chi^{(3)}|^2|E_\omega|^6.$$9, and $F(\omega)$0 by $F(\omega)$1. For the total susceptibilities, the reported enhancements are $F(\omega)$2 at the first harmonic, $F(\omega)$3 at the second, and $F(\omega)$4 at the third. The spectral resonance of the AFMR therefore drives a ladder of harmonic spectral lines at $F(\omega)$5, $F(\omega)$6, and $F(\omega)$7, enabling frequency conversion in the sub-terahertz regime (Deng et al., 28 Jun 2025).

In auroral plasma physics, the Electron Cyclotron Maser Instability provides a different harmonic-resonant structure. The resonance condition

$F(\omega)$8

defines which electrons couple to X-mode waves at the fundamental and higher cyclotron harmonics. The paper revisiting AKR argues that fundamental $F(\omega)$9 resonance on the confined lower X-mode branch can become hyperbolic rather than elliptic for sufficiently large parallel wavenumber, fitting the loss-cone boundary in the upward current region. Higher harmonics $$E_\omega^{\mathrm{loc}}=F(\omega)E_\omega^{\mathrm{inc}},\qquad I_{3\omega}(\omega)\propto |\chi^{(3)}|^2|F(\omega)|^6|E_\omega^{\mathrm{inc}}|^6.$$0 lie mainly on the free-space upper X-mode branch, with most harmonic radiation below $$E_\omega^{\mathrm{loc}}=F(\omega)E_\omega^{\mathrm{inc}},\qquad I_{3\omega}(\omega)\propto |\chi^{(3)}|^2|F(\omega)|^6|E_\omega^{\mathrm{inc}}|^6.$$1 and only a narrow band above. The second harmonic band below $$E_\omega^{\mathrm{loc}}=F(\omega)E_\omega^{\mathrm{inc}},\qquad I_{3\omega}(\omega)\propto |\chi^{(3)}|^2|F(\omega)|^6|E_\omega^{\mathrm{inc}}|^6.$$2 is argued to arise efficiently through nonlinear wave-wave interaction of confined fundamental X-modes (Baumjohann et al., 2023).

In geophysical fluid dynamics, harmonic-generation resonance produces a one-way spectral energy transfer. Certain internal gravity waves satisfy both $$E_\omega^{\mathrm{loc}}=F(\omega)E_\omega^{\mathrm{inc}},\qquad I_{3\omega}(\omega)\propto |\chi^{(3)}|^2|F(\omega)|^6|E_\omega^{\mathrm{inc}}|^6.$$3 and $$E_\omega^{\mathrm{loc}}=F(\omega)E_\omega^{\mathrm{inc}},\qquad I_{3\omega}(\omega)\propto |\chi^{(3)}|^2|F(\omega)|^6|E_\omega^{\mathrm{inc}}|^6.$$4, so nonlinear free-surface terms resonantly excite the second harmonic. The resulting amplitude equations show irreversible transfer of energy from the fundamental to the second harmonic, rendering such waves inherently unstable. The paper states that there are countably infinite such unstable waves, and that the mechanism does not obtain for a linearly-stratified fluid if a rigid lid or linearized free-surface boundary condition is used (Liang et al., 2016).

Linear oscillators with weak non-standard frequency modulation furnish yet another form. For

$$E_\omega^{\mathrm{loc}}=F(\omega)E_\omega^{\mathrm{inc}},\qquad I_{3\omega}(\omega)\propto |\chi^{(3)}|^2|F(\omega)|^6|E_\omega^{\mathrm{inc}}|^6.$$5

Jacobi-Anger expansion yields an infinite harmonic ladder

$$E_\omega^{\mathrm{loc}}=F(\omega)E_\omega^{\mathrm{inc}},\qquad I_{3\omega}(\omega)\propto |\chi^{(3)}|^2|F(\omega)|^6|E_\omega^{\mathrm{inc}}|^6.$$6

with ever expanding broadband frequency spectra as time progresses. A weakly damped oscillator then exhibits always two transient resonance captures involving two distinct harmonics, while the overall response decays as $$E_\omega^{\mathrm{loc}}=F(\omega)E_\omega^{\mathrm{inc}},\qquad I_{3\omega}(\omega)\propto |\chi^{(3)}|^2|F(\omega)|^6|E_\omega^{\mathrm{inc}}|^6.$$7. In the undamped case, the paper distinguishes simple and non-simple resonances; the latter can grow as $$E_\omega^{\mathrm{loc}}=F(\omega)E_\omega^{\mathrm{inc}},\qquad I_{3\omega}(\omega)\propto |\chi^{(3)}|^2|F(\omega)|^6|E_\omega^{\mathrm{inc}}|^6.$$8 under constructive phase conditions, in contrast to the classical $$E_\omega^{\mathrm{loc}}=F(\omega)E_\omega^{\mathrm{inc}},\qquad I_{3\omega}(\omega)\propto |\chi^{(3)}|^2|F(\omega)|^6|E_\omega^{\mathrm{inc}}|^6.$$9 growth under unmodulated harmonic forcing (Gendelman et al., 6 Jul 2026).

An analogous, explicitly nonphysical use appears in harmonic networks. There, fixed DCT filters decompose local image patches into spatial harmonics and a $$I(\omega)=A\frac{(\epsilon+q)^2}{\epsilon^2+1}+I_b, \qquad \epsilon=\frac{2(\omega-\omega_r)}{\Gamma},$$0 convolution learns how to combine them. This suggests a transfer of the idea from physical resonance to frequency-selective representation learning: the system learns how strongly to emphasize each harmonic mode in a local spectral basis (Ulicny et al., 2018).

5. Identification, simulation, and inverse analysis

Because harmonic spectral resonance often mixes amplitude, phase, and modal structure, its identification relies on specialized diagnostics. In the qBIC metasurface, continuous-wave frequency-domain simulations reproduce the quasi-stationary picosecond regime, while split-step time-domain Maxwell-Lorentz simulations track femtosecond buildup and decay of the quasi-BIC mode, instantaneous Kerr index changes, nonlinear phase accumulation, and broadband harmonic generation. The two methods are therefore complementary rather than interchangeable (Nardi et al., 2 Jul 2026).

In optically poled Si$$I(\omega)=A\frac{(\epsilon+q)^2}{\epsilon^2+1}+I_b, \qquad \epsilon=\frac{2(\omega-\omega_r)}{\Gamma},$$1N$$I(\omega)=A\frac{(\epsilon+q)^2}{\epsilon^2+1}+I_b, \qquad \epsilon=\frac{2(\omega-\omega_r)}{\Gamma},$$2 microresonators, pump-probe spectral mapping resolves both the pump and second-harmonic branches over probe-pump offsets exceeding $$I(\omega)=A\frac{(\epsilon+q)^2}{\epsilon^2+1}+I_b, \qquad \epsilon=\frac{2(\omega-\omega_r)}{\Gamma},$$3 GHz. The coupled-mode model uses intracavity amplitudes $$I(\omega)=A\frac{(\epsilon+q)^2}{\epsilon^2+1}+I_b, \qquad \epsilon=\frac{2(\omega-\omega_r)}{\Gamma},$$4 and $$I(\omega)=A\frac{(\epsilon+q)^2}{\epsilon^2+1}+I_b, \qquad \epsilon=\frac{2(\omega-\omega_r)}{\Gamma},$$5 with effective detunings $$I(\omega)=A\frac{(\epsilon+q)^2}{\epsilon^2+1}+I_b, \qquad \epsilon=\frac{2(\omega-\omega_r)}{\Gamma},$$6 and $$I(\omega)=A\frac{(\epsilon+q)^2}{\epsilon^2+1}+I_b, \qquad \epsilon=\frac{2(\omega-\omega_r)}{\Gamma},$$7, and the measured lock-in response maps the spectral evolution of the pump resonance and SH resonance during all-optical poling. This method provides direct access to a regime where the SH detuning reaches about $$I(\omega)=A\frac{(\epsilon+q)^2}{\epsilon^2+1}+I_b, \qquad \epsilon=\frac{2(\omega-\omega_r)}{\Gamma},$$8 GHz while the SH linewidth is about $$I(\omega)=A\frac{(\epsilon+q)^2}{\epsilon^2+1}+I_b, \qquad \epsilon=\frac{2(\omega-\omega_r)}{\Gamma},$$9 GHz, thereby making it possible to study broadband conversion well outside the narrow doubly resonant limit (Zhou et al., 7 Oct 2025).

Phase-sensitive HHG requires still different tools. In the tin-plasma study, low-resolution XUV-IR cross-correlation measures the femtosecond envelope, while RABITT extracts the relative spectral phases of neighboring harmonics. TDSE and recombination-dipole calculations are then used to attribute the observed group-delay shifts and loss of phase locking to the phase of the resonant recombination dipole moment (Haessler et al., 2012).

A mathematically distinct inverse approach is harmonic inversion. For ordinary resonances, the spectrum is modeled as

$q$0

while at exceptional points the correct representation becomes

$q$1

The corresponding time signal is no longer a sum of pure exponentials, but an exponential multiplied by a polynomial of degree $q$2. Extending nonlinear harmonic inversion to this setting allows one to extract the amplitudes of higher-order poles in the energy domain and the coefficients of the polynomial prefactor in the time domain, thereby identifying exceptional points from non-Lorentzian resonance spectra (Fuchs et al., 2014).

In superconducting coplanar resonators, identification means separating designed harmonic modes from parasitic resonances. The reported strategies include tracking resonance evolution as a function of temperature, magnetic field, and microwave power, and applying minute amounts of dielectric or ESR-active material to probe local electric- and magnetic-field maxima. Designed $q$3 CPW harmonics follow the odd-mode ladder $q$4 with $q$5, while spurious modes arise from standing waves in the substrate or housing and show different temperature, field, and nonlinear responses (Beydeda et al., 2023).

More generally, the forced-oscillator analysis of periodic solutions provides a minimal spectral criterion: for

$q$6

resonance occurs when the forcing period is an integer multiple of the natural period and the forcing has nonzero $q$7-projection on $q$8 or $q$9. Period matching alone is therefore insufficient (Benson et al., 2024).

6. Applications, misconceptions, and open directions

Across platforms, harmonic spectral resonance is used to control conversion efficiency, spectral selectivity, and dynamical response. High-$|E_{np}-E_0|+U_p=q\hbar\omega,$0 dielectric metasurfaces provide compact and efficient third-harmonic sources with controllable spectra for integrated nonlinear optics, multi-color sources, and chip-scale spectroscopy, while the femtosecond regime enables all-optical waveform shaping, temporal gating via resonant filtering, and intensity-dependent spectral encoding (Nardi et al., 2 Jul 2026). Resonance-gradient metasurfaces extend this logic to spectrally tunable third- and fifth-harmonic generation over a broad mid-IR band, and hyperbolic metamaterials with double second-harmonic resonance cones are proposed as platforms for subwavelength second- and higher-harmonic imaging microscopy (Jangid et al., 2023, Ceglia et al., 2013). In ZnO microwires, phase-controlled Fano line shapes point to on-chip nonlinear optical filters and switches based on harmonic resonances (Wang et al., 2019). In antiferromagnetic nano-devices, gigantic harmonic generation in AFMR motivates multilayer frequency amplifiers and converters in the sub-terahertz range (Deng et al., 28 Jun 2025).

Several recurrent misconceptions are explicitly contradicted by the literature. First, harmonic spectral resonance is not only an amplitude effect: in tin-plasma HHG and helium, resonance modifies spectral phase, group delay, quantum-path timing, and even phase locking (Haessler et al., 2012, Camp et al., 2015). Second, high $|E_{np}-E_0|+U_p=q\hbar\omega,$1 does not by itself solve spectral coverage; in uniform resonant metasurfaces it narrows the operating band, and broad coverage is obtained instead by spatial gradients or by dynamically reconfigurable photo-induced $|E_{np}-E_0|+U_p=q\hbar\omega,$2 (Jangid et al., 2023, Zhou et al., 7 Oct 2025). Third, exact doubly resonant alignment is not always required: in Si$|E_{np}-E_0|+U_p=q\hbar\omega,$3N$|E_{np}-E_0|+U_p=q\hbar\omega,$4 microresonators, photo-induced SHG can persist from the preferred doubly resonant condition into a highly detuned state, although stable all-optical poling requires $|E_{np}-E_0|+U_p=q\hbar\omega,$5 and $|E_{np}-E_0|+U_p=q\hbar\omega,$6 (Zhou et al., 7 Oct 2025). Fourth, period matching alone does not guarantee resonance for general forcing; the decisive quantity is the forcing’s projection on the natural mode (Benson et al., 2024). Fifth, in linearly stratified fluids the harmonic-generation instability disappears if the nonlinear free-surface boundary condition is replaced by a rigid lid or linear form (Liang et al., 2016).

Open problems are correspondingly platform-specific but conceptually aligned. In nonlinear metasurfaces, balancing very high $|E_{np}-E_0|+U_p=q\hbar\omega,$7 against fabrication sensitivity and nonlinear instabilities, modeling fully nonperturbative harmonic cascades, and controlling spatio-temporal dynamics remain open design problems (Nardi et al., 2 Jul 2026). In optically poled Si$|E_{np}-E_0|+U_p=q\hbar\omega,$8N$|E_{np}-E_0|+U_p=q\hbar\omega,$9, long-term stability and precise control of the photo-induced $\omega$00 grating remain active issues (Zhou et al., 7 Oct 2025). In multilevel and confined systems, cascade resonance shows that cutoff laws can be redefined by discrete level structure rather than by continuum excursion energy, suggesting broader relevance to nanoribbon and quantum-dot geometries where only a few subbands satisfy frequency-matching conditions (Zhang et al., 2021). In strong-field spin systems, the detuning dependence of multiple harmonic generation in electric dipole spin resonance indicates a need for fuller theories incorporating realistic noise and detailed Landau-Zener dynamics at anticrossings (Stehlik et al., 2013).

A plausible implication is that harmonic spectral resonance is best treated as a unifying spectral-structure principle rather than as a narrowly optical effect. In every realization represented here, harmonic output is determined by how a driven system reorganizes energy among discrete or structured spectral channels under constraints set by resonance, interference, and nonlinearity.