Ultrafast Third-Harmonic Spectral Modulation and Self-Action in Resonant Nonlocal Metasurfaces

Abstract: Quasi-bound states in the continuum enable exceptional field confinement, strongly reducing the pump intensity threshold for nonlinear light-matter interaction in dielectric metasurfaces. As a result, nonlinear self-action effects, often elusive in bulk nonlinear media, emerge at moderate excitation intensities. Here, we show how pulse duration and resonant coupling govern the nonlinear self-action mechanism in resonantly enhanced third-harmonic (TH) generation from a nonlocal metasurface. We identify two excitation regimes that interact differently with the resonant mode, revealing complementary intensity-dependent responses. Under spectrally narrow picosecond excitation, resonance-enhanced TH generation shows pronounced deviations from cubic scaling at high intensities. In contrast, broadband femtosecond excitation transiently drives the resonance, encoding the nonlinear response in the spectral reshaping and broadening of the TH signal. Simulations reproduce both regimes: continuous-wave modeling captures picosecond power scaling and higher-harmonic interactions, while time-domain simulations resolve femtosecond dynamics. These results clarify nonlinear self-action in metasurfaces featuring quasi-bound states in the continuum, linking strong field confinement to conversion efficiency, scaling behavior, and distinct spectral dynamics under different excitation conditions. This work sheds light on the interplay between pulse duration, bandwidth, and resonant coupling in high-Q nonlocal dielectric metasurfaces, advancing their use in ultrafast and nonlinear nanophotonics.

• The paper demonstrates that picosecond excitation leads to sub-cubic THG scaling due to Kerr-induced refractive index changes and self-action effects.
• The femtosecond regime shows dynamic spectral modulation and broadening, reflecting transient coupling between the pump spectrum and the qBIC resonance.
• A unified modeling framework combining CW and time-domain hydrodynamic simulations validates the nonlinear dynamics in high-Q metasurfaces.

Ultrafast Third-Harmonic Spectral Modulation and Nonlinear Self-Action in Resonant Nonlocal Metasurfaces

Introduction and Motivation

This paper presents an in-depth analysis of ultrafast third-harmonic generation (THG) and nonlinear self-action phenomena in high-Q dielectric metasurfaces supporting quasi-bound states in the continuum (qBICs). The motivation is grounded in the necessity to harness strong nonlinear optical processes at the nanoscale with moderate excitation intensities, circumventing the shortcomings of bulk materials—primarily their limited nonlinear response, which traditionally demands either high pump powers or extended propagation lengths.

Dielectric metasurfaces with high-Q qBIC resonances offer a promising modality, leveraging field localization to enhance nonlinear responses dramatically. However, as field confinement intensifies, nonlinear self-action effects become prominent, significantly affecting the spectral properties and efficiency of nonlinear frequency conversion processes. This work systematically investigates how the interplay between excitation pulse duration, bandwidth, and resonant coupling dynamics crucially shapes the nonlinear response and its manifestations across different excitation regimes.

Experimental Platform: Metasurface Design and Linear Characterization

The experimental platform consists of a one-dimensional periodic array of amorphous silicon nanobars with engineered in-plane asymmetry. The symmetry breaking within each unit cell enables efficient free-space coupling to qBICs under TE polarization. The qBIC supports a spatially structured near-field profile conducive to strong field enhancement while maintaining low radiative losses, as confirmed by both COMSOL simulations and transmission spectroscopy.

Spectral characterization reveals a pronounced Fano resonance aligned with the qBIC mode and a sharply polarization-dependent response, with the TH intensity exhibiting a pronounced resonant peak for TE polarization. The metasurface's design parameters—nanobar widths, heights, gap, and periodicity—are meticulously controlled to facilitate robust and reproducible excitation of the targeted qBIC mode in the relevant spectral window.

Nonlinear Regimes: Picosecond vs. Femtosecond Excitation

Picosecond Excitation: Resonance-Resolved Nonlinear Dynamics

Under picosecond excitation, the pump bandwidth is markedly narrower than the qBIC linewidth. This regime allows quasi-stationary, resonant driving over the full mode lifetime and is amenable to continuous-wave (CW) modeling. Experimental measurements and simulations reveal that resonantly enhanced THG exhibits pronounced deviation from the expected cubic scaling (i.e., $I_{THG} \propto I_{pump}^3$) at higher pump intensities. In particular, as the excitation approaches the resonance peak, power-dependent redshift and spectral broadening of the TH response are observed, with the scaling exponent decreasing towards nearly linear behavior at the highest intensities.

This sub-cubic scaling is attributed to intensity-dependent refractive index changes (Kerr effect), self- and cross-phase modulation between the fundamental and harmonic fields, and the participation of higher-order nonlinear susceptibilities. The modeling incorporates all relevant nonlinear polarization contributions. The data unambiguously demonstrate that self-action arises from field confinement and not from pump depletion or material damage.

Femtosecond Excitation: Ultrafast Transient Coupling and Spectral Modulation

In the femtosecond regime, the pump bandwidth exceeds the qBIC resonance width, resulting in a transient, spectrally diverse interaction. The nonlinear response is now encoded in the spectral structure of the TH emission itself. Spectrometer-resolved measurements, supported by full time-domain Maxwell-Lorentz hydrodynamic simulations, reveal that only a subset of the pump’s spectrum couples efficiently to the resonance. As the pump is tuned through the resonance and intensity increases, the TH spectrum exhibits pronounced modulation, broadening, and reshaping. This transient nonlinear spectral modulation is a direct manifestation of power-driven refractive index shifts dynamically altering the qBIC conditions during pulse interaction.

Notably, the onset of observable self-action effects in this regime requires significantly higher intensities than in the picosecond case, reflecting the reduced spectral overlap with the resonance. The impact on power-law scaling is less direct than in the quasi-stationary case, with the main observable being the time-dependent spectral evolution of the outgoing harmonics.

Theoretical Modeling and Simulation Framework

The analysis employs two complementary numerical approaches:

• CW, frequency-domain modeling for the picosecond regime, allowing efficient calculation of resonance-driven nonlinearities and spectral transformations under quasi-continuous excitation.
• Hydrodynamic, Maxwell-Lorentz time-domain modeling for the femtosecond regime, accounting for the full ultrafast dynamics, nonlinear refractive index changes, and interactions between all field components.

In both cases, the macroscopic polarization dynamics are solved self-consistently with Maxwell’s equations, encompassing both the linear and nonlinear (primarily third-order) response of amorphous silicon. The modeling framework robustly captures the spectral, temporal, and scaling properties observed experimentally, validating the underlying physical mechanisms.

Strong Numerical Results and Bold Claims

• Deviation from cubic scaling: Near qBIC resonance, the THG process transitions from $I_{THG} \sim I_{pump}^3$ to nearly linear scaling ($I_{THG} \sim I_{pump}^1$) at high input powers under picosecond excitation, confirming the dominance of nonlinear self-action over the conventional third-order response.
• Spectral reshaping in femtosecond regime: The observed intensity-dependent spectral broadening and modulation of TH emission occur at intensities two to three orders of magnitude higher than in the picosecond case and only for pump wavelengths spectrally aligned with the qBIC resonance.
• Unified framework: The nonlinear self-action arises solely from intensity-induced refractive index changes and field-driven resonance modulation, not from inherent property limitations or extrinsic loss mechanisms.

Implications and Prospects for Nonlinear Nanophotonics

These results provide a comprehensive foundation for understanding nonlinear self-action in high-Q nanostructures, emphasizing the sensitivity of these platforms to excitation temporal structure. Practically, the findings furnish concrete guidelines for the design of ultrafast nanophotonic devices: the ability to toggle between power-scaling modulation and spectral filtering/reshaping regimes by adjusting the pump parameters enables new classes of active devices for frequency conversion, ultrafast modulation, and signal processing.

Theoretically, the work situates qBIC-based metasurfaces within the broader context of time-dependent photonic systems. The intrinsic time-varying refractive properties under strong excitation imply that many functionalities conventionally attributed to externally modulated systems (e.g., nonreciprocity, parametric amplification) can emerge natively from the nonlinear dynamics of high-Q metasurfaces, obviating the need for complex material engineering.

Looking forward, the demonstrated sensitivity to pump temporal and spectral structure suggests further opportunities:

• Advanced pulse shaping (chirping, multiple color schemes, pump-probe techniques) could be used to selectively control resonance dynamics and harmonic emission.
• Strong-coupling and spatiotemporal modulation effects could be explored in cascaded or multiplexed metasurface arrays.
• The design paradigm could extend to other nonlinear processes, potentially scaling to frequency comb generation or ultrafast switching.

Conclusion

This study establishes that field-confinement-driven nonlinear self-action in qBIC-based dielectric metasurfaces fundamentally shapes both efficiency and the scaling of third-harmonic signals, with distinct observable phenomena arising in narrowband and broadband excitation regimes. The findings support a unified view where time-dependent, intensity-driven resonance modulation is intrinsic to the operation of strongly nonlinear, high-Q nanoscale devices. These insights will inform the engineering of next-generation ultrafast and nonlinear functional photonic systems.