High-Harmonic Spectroscopy (HHS)

• High-Harmonic Spectroscopy (HHS) is a strong-field technique that uses high-harmonic generation to capture ultrafast electron, nuclear, and lattice dynamics in various media.
• It decodes information encoded in harmonic amplitude, phase, polarization, and cutoff behavior to reveal electronic structures, charge migration, and phase transitions.
• Recent developments integrate multidimensional control and advanced retrieval algorithms to isolate specific microscopic processes with femtosecond to attosecond resolution.

High-Harmonic Spectroscopy (HHS) is the use of high-harmonic generation (HHG) as a spectroscopic probe of matter under strong-field driving. In gases and molecules, HHS exploits the coherent sequence of tunnel ionization, continuum acceleration, and recombination; in solids, it reads out coupled intraband and interband dynamics; in liquids, it exposes scattering-limited electron motion. Across these settings, harmonic amplitude, phase, polarization, helicity, alignment dependence, sidebands, and cutoff behavior encode electronic structure, nuclear and lattice dynamics, charge migration, phase transitions, and low-energy electron transport on femtosecond to attosecond time scales (Negro et al., 2014, Klemke et al., 2018, Mondal et al., 2022).

1. Physical basis and defining observables

In the canonical gas-phase picture, HHG follows the three-step model: an electron ionizes from a bound state, is accelerated by the laser field, and recombines with the parent ion, emitting a photon whose energy is set by the instantaneous kinetic energy plus the ionization potential. This semiclassical description is closely connected to the Strong Field Approximation and gives the familiar cutoff relation $E_c = I_p + 3.17 U_p$ in the above-threshold regime. Near threshold, however, Coulomb forces, resonances, and excitation-assisted pathways become essential, so the standard cutoff law becomes less predictive and the harmonic response becomes spectroscopically rich (Soifer et al., 2010, Negro et al., 2014).

The emitted spectrum is typically obtained from a Fourier transform of the dipole acceleration, dipole moment, or current. In molecular HHG, one may write

$$S(\omega)\propto \left|\int dt\, a(t)e^{i\omega t}\right|^2, \qquad S(\omega)\propto \omega^2\left|\int dt\, d(t)e^{i\omega t}\right|^2.$$

In solids, the emitted field is commonly expressed through the nonlinear current, with

$$\mathbf{E}(\omega)\propto i\omega \mathbf{J}(\omega), \qquad \mathbf{J}(t)=\mathbf{J}_{\mathrm{intra}}(t)+\mathbf{J}_{\mathrm{inter}}(t),$$

so that the polarization state and spectral distribution of the harmonics directly reflect the relative phases and amplitudes of interband polarization and intraband current (Hamer et al., 2022, Klemke et al., 2018).

What distinguishes HHS from HHG viewed purely as a source is that the harmonic signal is treated as a structured observable rather than a by-product. In molecules, alignment dependence, ellipticity response, spectral minima, and harmonic phase encode recombination dipoles, orbital symmetry, and ultrafast dynamics. In solids, crystal symmetry, Berry curvature, interband coupling, and band anisotropy appear in harmonic selection rules, polarization states, and cutoff scaling. In liquids, the cutoff can become a characteristic property of the medium because electron trajectories are limited by scattering to a characteristic length connected to the mean free path (Soifer et al., 2010, Mrudul, 2022, Mondal et al., 2022).

2. Molecular and gas-phase spectroscopy

Molecular HHS first developed as a structural probe of aligned molecules. In aligned N$_2$ and O$_2$, near-threshold HHS revealed two distinct channels: a three-step ionization–acceleration–recombination channel associated with long trajectories, and an excitation-assisted multiphoton channel associated with short trajectories. Their separation required combined control of alignment, ellipticity, and trajectory selection, and it established that near-threshold harmonics can probe excited-state symmetries and resonant couplings rather than only ground-state recombination (Soifer et al., 2010).

The extension of HHS to molecular imaging used alignment-resolved HHG, single-molecule emission reconstruction, and tomographic inversion of the recombination dipole to recover highest occupied molecular orbitals in polyatomics. Mid-infrared driving at $1450$ nm, combined with phase retrieval and an iterative Kaczmarz algorithm, enabled orbital tomography in N$_2$O and C$_2$H$_2$. In N$_2$O, a pronounced amplitude minimum near $$S(\omega)\propto \left|\int dt\, a(t)e^{i\omega t}\right|^2, \qquad S(\omega)\propto \omega^2\left|\int dt\, d(t)e^{i\omega t}\right|^2.$$0–$$S(\omega)\propto \left|\int dt\, a(t)e^{i\omega t}\right|^2, \qquad S(\omega)\propto \omega^2\left|\int dt\, d(t)e^{i\omega t}\right|^2.$$1 eV and an associated $$S(\omega)\propto \left|\int dt\, a(t)e^{i\omega t}\right|^2, \qquad S(\omega)\propto \omega^2\left|\int dt\, d(t)e^{i\omega t}\right|^2.$$2 rad phase jump were intensity-invariant across $$S(\omega)\propto \left|\int dt\, a(t)e^{i\omega t}\right|^2, \qquad S(\omega)\propto \omega^2\left|\int dt\, d(t)e^{i\omega t}\right|^2.$$3–$$S(\omega)\propto \left|\int dt\, a(t)e^{i\omega t}\right|^2, \qquad S(\omega)\propto \omega^2\left|\int dt\, d(t)e^{i\omega t}\right|^2.$$4 W/cm$$S(\omega)\propto \left|\int dt\, a(t)e^{i\omega t}\right|^2, \qquad S(\omega)\propto \omega^2\left|\int dt\, d(t)e^{i\omega t}\right|^2.$$5, supporting a HOMO-dominated interpretation; in C$$S(\omega)\propto \left|\int dt\, a(t)e^{i\omega t}\right|^2, \qquad S(\omega)\propto \omega^2\left|\int dt\, d(t)e^{i\omega t}\right|^2.$$6H$$S(\omega)\propto \left|\int dt\, a(t)e^{i\omega t}\right|^2, \qquad S(\omega)\propto \omega^2\left|\int dt\, d(t)e^{i\omega t}\right|^2.$$7, the reconstruction was presented as preliminary because multi-electron influence could not be cleanly excluded (Negro et al., 2014).

Time-resolved molecular HHS broadened the field from structure to dynamics. In NO$$S(\omega)\propto \left|\int dt\, a(t)e^{i\omega t}\right|^2, \qquad S(\omega)\propto \omega^2\left|\int dt\, d(t)e^{i\omega t}\right|^2.$$8, high-harmonic transient grating spectroscopy used a spatial excitation grating to convert pump-induced vibronic motion into background-free diffracted harmonics. The first-order diffracted signal was governed by the coherent difference between emissions from excited and ground-state molecules and was particularly sensitive to bending-induced phase shifts $$S(\omega)\propto \left|\int dt\, a(t)e^{i\omega t}\right|^2, \qquad S(\omega)\propto \omega^2\left|\int dt\, d(t)e^{i\omega t}\right|^2.$$9 on the first excited diabatic surface, thereby tracking wavepacket reshaping through nonadiabatic coupling (Ruf et al., 2012). In 1,3-cyclohexadiene, time-resolved HHS followed a bond-making photochemical reaction: the electronically excited state relaxed almost completely within $$\mathbf{E}(\omega)\propto i\omega \mathbf{J}(\omega), \qquad \mathbf{J}(t)=\mathbf{J}_{\mathrm{intra}}(t)+\mathbf{J}_{\mathrm{inter}}(t),$$0 fs to vibrationally excited ground-state CHD, ring opening occurred around $$\mathbf{E}(\omega)\propto i\omega \mathbf{J}(\omega), \qquad \mathbf{J}(t)=\mathbf{J}_{\mathrm{intra}}(t)+\mathbf{J}_{\mathrm{inter}}(t),$$1 fs, and subsequent conformer changes appeared on the $$\mathbf{E}(\omega)\propto i\omega \mathbf{J}(\omega), \qquad \mathbf{J}(t)=\mathbf{J}_{\mathrm{intra}}(t)+\mathbf{J}_{\mathrm{inter}}(t),$$2–$$\mathbf{E}(\omega)\propto i\omega \mathbf{J}(\omega), \qquad \mathbf{J}(t)=\mathbf{J}_{\mathrm{intra}}(t)+\mathbf{J}_{\mathrm{inter}}(t),$$3 ps timescale (Kaneshima et al., 2017).

Angle-resolved HHS also became a route to reconstructing rotational quantum dynamics. He et al. combined angle-resolved harmonic measurements with a machine-learning-enabled inversion based on simulated annealing to retrieve the full spatiotemporal evolution of a rotational wave packet in N$$\mathbf{E}(\omega)\propto i\omega \mathbf{J}(\omega), \qquad \mathbf{J}(t)=\mathbf{J}_{\mathrm{intra}}(t)+\mathbf{J}_{\mathrm{inter}}(t),$$4. The reconstructed molecular axis distribution, the delay-dependent populations, and the linear phase evolution of individual rotational components were in good agreement with numerical simulations, yielding a genuine molecular rotation movie over the rotational revival (He et al., 2019).

3. Condensed-phase implementations

In solids, HHS probes the coupled dynamics of carriers in periodic bands rather than recollision with an isolated parent ion. Polarization-state-resolved HHS in silicon and quartz showed that the polarization of solid-state harmonics is not fixed by crystal symmetry alone: it can be dynamically controlled by driver ellipticity, intensity, and crystal orientation because interband and intraband currents interfere vectorially within a single optical cycle. Under suitable conditions, elliptically driven fields produced circular harmonics that were substantially brighter than those generated by circular drivers, with measured improvements of $$\mathbf{E}(\omega)\propto i\omega \mathbf{J}(\omega), \qquad \mathbf{J}(t)=\mathbf{J}_{\mathrm{intra}}(t)+\mathbf{J}_{\mathrm{inter}}(t),$$5 for HH5 and HH7 and $$\mathbf{E}(\omega)\propto i\omega \mathbf{J}(\omega), \qquad \mathbf{J}(t)=\mathbf{J}_{\mathrm{intra}}(t)+\mathbf{J}_{\mathrm{inter}}(t),$$6 for HH9 under specific conditions in Si (Klemke et al., 2018).

Time-resolved HHS in correlated materials established a different spectroscopic role. In VO$$\mathbf{E}(\omega)\propto i\omega \mathbf{J}(\omega), \qquad \mathbf{J}(t)=\mathbf{J}_{\mathrm{intra}}(t)+\mathbf{J}_{\mathrm{inter}}(t),$$7, harmonics below the bandgap were dominated by intraband dynamics and followed the photoinduced insulator-to-metal transition. The harmonic yield separated a fast recovery channel of about $$\mathbf{E}(\omega)\propto i\omega \mathbf{J}(\omega), \qquad \mathbf{J}(t)=\mathbf{J}_{\mathrm{intra}}(t)+\mathbf{J}_{\mathrm{inter}}(t),$$8–$$\mathbf{E}(\omega)\propto i\omega \mathbf{J}(\omega), \qquad \mathbf{J}(t)=\mathbf{J}_{\mathrm{intra}}(t)+\mathbf{J}_{\mathrm{inter}}(t),$$9 fs from a slower channel of about $_2$0–$_2$1 ps associated with formation of the metastable monoclinic metal, and it exhibited oscillations at $_2$2 THz corresponding to a coherent phonon of the monoclinic phase (Bionta et al., 2020).

Coherent lattice dynamics became a major solid-state application. In monolayer graphene, coherent excitation of the in-plane $_2$3 phonons produced sidebands around harmonic peaks at $_2$4, with $_2$5 meV. Their polarization followed dynamical-symmetry selection rules: for the iLO mode, even-$_2$6 sidebands were parallel and odd-$_2$7 sidebands perpendicular to the probe, whereas for iTO all sidebands were parallel in the linearly polarized scheme (Rana et al., 2022). With a circularly polarized probe, coherent phonons reduced graphene’s dynamical symmetry from $_2$8 to $_2$9, lifted the $_2$0 selection rule to allow all odd harmonics, and generated sidebands whose ellipticity and phase identified whether iLO or iTO had been excited (Rana et al., 2022).

In $_2$1-quartz, HHS was used as a background-free, energy-domain probe of coherent phonons. By tracking absolute shifts of harmonic line centers, rather than relative intensity modulations, the method extracted optical phonons at $_2$2 cm$_2$3 and $_2$4 cm$_2$5, determined deformation-potential-calibrated oxygen displacements of tens of picometers, and resolved channel-specific four-phonon scattering through time-windowed Gabor analysis (Zhang et al., 2024).

Excitonic and topological responses further extended solid-state HHS. Jensen, Madsen, Rubio, and Tancogne-Dejean showed that a finite exciton population imprints an exciton line at $_2$6 and sidebands at $_2$7 into the harmonic spectrum of a one-dimensional solid, and that the integrated yields of these features scale as $_2$8, allowing exciton-population retrieval (Jensen et al., 2023). In two-dimensional materials, HHS was shown to encode interband coupling, band anisotropy, Berry-curvature-driven anomalous velocity, valley selectivity, and defect-mediated many-body effects (Mrudul, 2022). In topological materials, helicity, circular dichroism, and ellipticity dependence captured topological phases in Chern and $_2$9 insulators, whereas no clear harmonic signatures of topology emerged in the breathing Kagome HOTI bulk response, even though the topological semimetal phase displayed a harmonic enhancement by two to three orders of magnitude (Lorenzo et al., 21 Aug 2025).

Liquids represent a distinct condensed phase. In bulk liquids, HHS revealed a scattering-limited mechanism in which the cutoff becomes wavelength-independent beyond a threshold intensity and is governed instead by a characteristic excursion length associated with the electron mean free path. Water and ethanol exhibited cutoffs of $1450$0 eV and $1450$1 eV, respectively, and the results were validated by TDDFT calculations on liquid supercells and large clusters (Mondal et al., 2022).

4. Specialized modalities and multidimensional extensions

A recurrent theme in HHS is the design of observables that separate otherwise entangled ionization, propagation, and recombination physics. High-harmonic sideband spectroscopy (HHSS) was introduced as a modality for pump-launched charge migration in molecules. In bromobutadiyne, a localized Br-centered hole oscillating at $1450$2 eV modulated the HHG dipole driven by a delayed mid-infrared field polarized perpendicular to the molecular backbone. This geometry suppressed MIR-triggered charge migration and produced sidebands at $1450$3 and $1450$4, with stacked-spectrum offsets $1450$5 and $1450$6. Laser-frequency scans and delay-domain Fourier analysis retrieved both $1450$7 and $1450$8, and the presence of the second component was interpreted as a signature of particle-like charge migration (Hamer et al., 2022).

Bi-circular HHS developed an alternative separation strategy through symmetry. In counter-rotating $1450$9 and $_2$0 fields, the threefold dynamical symmetry of the combined driver suppresses $_2$1 harmonics and permits only $_2$2 orders, with opposite helicities for the two series. The intensity ratio $_2$3 was shown to encode angular-momentum propensity rules, atom-specific recombination matrix elements, and interference of co- and counter-rotating channels, thereby turning helicity-resolved harmonic pairs into a quantitative spectroscopic observable (Baykusheva et al., 2017).

Smirnova and Serbinenko proposed a more general multidimensional HHS framework in which orthogonally polarized multicolor control fields modulate the transverse motion of the electron and thereby disentangle ionization and recombination. In that formulation, optimized odd- and even-harmonic signals give access to ionization delays, recombination times, and sub-cycle ionization rates, while additional control dimensions lift degeneracies that remain hidden in two-dimensional schemes (Serbinenko et al., 2013). A complementary development is multidimensional high-harmonic echo spectroscopy, where sequences of few-cycle strong infrared pulses generate phase-matched EUV echo signals. Because the effective strong-field dipole is broadband and time dependent, the echo can display an ultra-delayed partial rephasing at $_2$4, a behavior not found in standard 2D optical spectroscopy (Jiang et al., 2023).

These modalities share a common logic: HHS becomes more informative when harmonic emission is not treated only as an intensity spectrum, but as a multidimensional, phase-sensitive observable in which sidebands, polarization, helicity, and delay-domain structure are engineered to isolate specific microscopic processes.

5. Instrumentation, calibration, and retrieval

The expansion of HHS into a quantitative spectroscopy has depended on equally substantial advances in apparatus design and inverse methods. A cryogenic solid-state HHG spectrometer based on a vacuum chamber and a closed-cycle helium cryostat demonstrated that the interaction region can remain cryogenic under high-intensity femtosecond driving, enabling temperature-dependent sHHG down to $_2$5 K. In ZnO, in situ optical thermometry based on the Varshni relation yielded $_2$6 eV and verified that the local sample temperature remained cryogenic even during nonperturbative MIR irradiation (Kohrell et al., 2023).

A later precision bulk-solids apparatus emphasized quantitative field and orientation control. It combined dispersion-neutral intensity control for few-cycle pulses, a vacuum HHG module with sub-micrometer and sub-degree sample positioning, an imaging assembly for focal-spot stabilization and spatial filtering, synchronized UV/VUV and EUV spectrometers, and absolute field calibration transferred from gas-phase attosecond streaking. Demonstrations included field-dependent cutoffs in a $_2$7 nm fused-silica membrane and $_2$8-periodic orientation dependence in MgO(100), with statistical errors of only a few percent in angular scans (Mandal et al., 15 Dec 2025).

Equally important are retrieval algorithms. In orbital tomography, phase retrieval combined with an iterative Kaczmarz algorithm reconstructed single-molecule amplitude and phase from revival-resolved macroscopic spectra (Negro et al., 2014). In rotational-wave-packet imaging, simulated annealing inverted angle-resolved harmonic yields into full spatiotemporal molecular axis distributions (He et al., 2019). In phonon spectroscopy, straightforward time-windowed Gabor analysis exposed mode frequencies, lifetimes, and multi-phonon channels directly in the energy domain (Zhang et al., 2024). In HHSS, k-means clustering of sideband peaks and delay-domain Fourier transforms enabled precision retrieval of charge-migration frequencies (Hamer et al., 2022). HHS has therefore evolved not only through new physical observables, but through increasingly explicit inversion protocols that connect those observables back to microscopic dynamics.

6. Interpretation, limitations, and frontier questions

Several interpretive limits recur across the field. Near threshold, Coulomb-free SFA pictures break down because below-threshold emission can arise from Coulomb-enabled returns and excitation-assisted pathways, so channel attribution requires trajectory selection, ellipticity control, and Coulomb-inclusive modeling (Soifer et al., 2010). In molecular imaging, intensity-invariant spectral minima support a single-active-electron interpretation in some systems, but not in all; the C$_2$9H$_2$0 tomography example remained provisional because multi-electron effects could not be excluded (Negro et al., 2014).

In time-resolved modalities, coherence loss is often decisive. In HHSS, sidebands persist for CM lifetimes down to about $_2$1 fs, but unaligned samples average them out, and finite alignment dispersion reduces contrast even when retrieval remains possible up to $_2$2 FWHM (Hamer et al., 2022). In liquids, the interpretation of cutoffs as mean-free-path signatures relies on a scattering-limited trajectory picture that is supported by TDDFT but still abstracts away channel-specific scattering processes (Mondal et al., 2022). In solids, above-bandgap propagation, detector bandwidth, and dephasing complicate polarization-resolved and phonon-resolved measurements, even when thin samples and ab initio simulations indicate that the dominant signatures are microscopic (Klemke et al., 2018, Zhang et al., 2024).

A broader conceptual caution concerns topology. HHS can diagnose Chern phases and some $_2$3 responses through helicity, circular dichroism, and ellipticity dependence, but the same observables do not automatically reveal higher-order topology in bulk-only calculations. The breathing Kagome results therefore argue against treating any single harmonic observable as a universal topological invariant and instead motivate more selective schemes that couple directly to corner or edge states (Lorenzo et al., 21 Aug 2025).

The field nonetheless points toward a consistent frontier. Phase-resolved HHSS can access $_2$4 with sub-cycle precision (Hamer et al., 2022). Cryogenic sHHG opens low-temperature correlated phases to quantitative study (Kohrell et al., 2023). Liquid-phase HHS extends attosecond spectroscopy to native condensed environments (Mondal et al., 2022). Multidimensional echo schemes and structured control fields aim to separate ionization, propagation, recombination, and many-body rearrangement within a single experiment (Serbinenko et al., 2013, Jiang et al., 2023). Taken together, these developments define HHS not as a single technique but as a family of strong-field spectroscopies whose unifying principle is that the emitted high harmonics are a phase-sensitive record of ultrafast matter dynamics.