Efficient high-harmonic generation in van der Waals ferroelectric NbOI$_2$ crystals

Layered NbOX$_2$ ($X=\mathrm{Cl,,Br,,I}$), a member of the van der Waals ferroelectric family, exhibits intrinsic ferroelectricity and pronounced nonlinear optical responses, making it a promising candidate for integrated nanophotonics applications. While previous studies have emphasized the material’s strong second-order nonlinear responses, higher-order nonlinear responses are still mostly unexplored. This work systematically investigates NbOI$_2$ using high harmonic generation (HHG) spectroscopy. Driven by an intense mid-infrared laser field centered at $\sim4μ\mathrm{m}$ wavelength, highly anisotropic odd- and even-order harmonics up to the 16th order are generated at a low peak intensity of $0.4\mathrm{TW,cm^{-2}}$, extending beyond the material’s bandgap. Both bulk and flake forms of NbOI$_2$ display pronounced harmonic emission from the near-infrared to the deep-ultraviolet spectral region, with a notably high overall conversion efficiency compared to other known materials. Polarization-resolved measurements reveal that even-order harmonics remain aligned with the crystal polar axis regardless of the driving-field orientation, whereas odd-order harmonics are dynamically affected. First-principles calculations suggest that the flat valence band associated with Peierls dimerization enhances HHG efficiency through electron correlation. These findings provide fresh perspectives on HHG in van der Waals ferroelectric materials and facilitate the development of compact and tunable quantum light sources.

💡 Research Summary

This work presents a comprehensive study of high‑harmonic generation (HHG) in the van Waals ferroelectric crystal NbOI₂, exploring both bulk and exfoliated thin‑film forms. Using intense mid‑infrared femtosecond pulses centered at ~4 µm (photon energy ≈0.31 eV) with a peak intensity of only 0.4 TW cm⁻², the authors generate a series of odd and even harmonics ranging from the 4th to the 16th order (≈1.2 eV to ≈5 eV). The harmonic emission extends well beyond the indirect bandgap of NbOI₂ (2.24 eV) and shows no fluorescence background, confirming a pure nonlinear process.

The crystal structure of NbOI₂ belongs to the monoclinic C2 space group and features in‑plane ferroelectric polarization along the b‑axis together with one‑dimensional Nb‑Nb dimer chains that produce a flat valence band across the Brillouin zone. This flat band, originating from Peierls dimerization, yields a high density of states and strong electron correlation, which the authors identify as a key factor enhancing HHG efficiency.

Intensity‑dependent measurements reveal a transition from perturbative scaling (I^q for low orders at low intensity) to a universal cubic scaling (I³) for all harmonics as the driving field increases, indicating that the process is governed by strong‑field electron dynamics rather than conventional perturbative nonlinear optics. The conversion efficiency of the 4th harmonic reaches 2.2 × 10⁻⁵, while the 4th–6th harmonics achieve efficiencies in the 10⁻⁵–10⁻⁶ range—orders of magnitude higher than those reported for typical solid‑state HHG media such as ZnO or MoS₂ under comparable conditions.

Thickness dependence is systematically investigated using samples from a few hundred nanometers up to tens of micrometers. While bulk crystals (≈30 µm) produce broader harmonic peaks, thin films (especially ~150 nm) exhibit reduced self‑absorption, leading to enhanced efficiency for higher‑order harmonics. Low‑order harmonics display oscillatory intensity variations with thickness, attributed to interference effects and changes in coherence length.

Polarization‑resolved studies show a striking dichotomy: even‑order harmonics are always linearly polarized along the crystal’s polar axis, regardless of the driving‑field orientation, reflecting the broken inversion symmetry along that direction. Odd‑order harmonics, however, follow the driving polarization only loosely; their polarization aligns with whichever crystal axis (polar or non‑polar) is closer to the driving field direction, and it flips when the driving polarization crosses ±45°. This behavior suggests that electron–hole trajectories are constrained by the underlying lattice symmetry, and that odd and even harmonics arise from distinct microscopic mechanisms.

First‑principles density‑functional calculations confirm the presence of a nearly dispersionless valence band associated with Nb‑Nb dimers. Time‑dependent simulations of the driven electron dynamics demonstrate that the flat band prolongs the electron acceleration phase and suppresses rapid recombination, thereby boosting harmonic yield.

Overall, NbOI₂ emerges as an exceptionally efficient HHG platform among van Waals materials, combining strong ferroelectric polarity, a flat correlated band, and low‑threshold driving fields. The ability to generate high‑order harmonics up to the deep‑UV with high conversion efficiency, together with the clear polarization control linked to crystal symmetry, opens pathways for compact on‑chip attosecond light sources, tunable quantum‑light emitters, and novel nonlinear photonic devices.