HHG in Quantum Dots: Confinement & Nonlinearity

• High-order harmonic generation in quantum dots is the process where intense laser fields interact with nanoscale confined carriers, producing coherent multiples of the laser frequency.
• Experimental setups use mid-infrared pulses and size-tuned quantum dots to optimize harmonic yield, cutoff energy, and spectral structure amidst strong-field effects.
• Theoretical frameworks like tight-binding and density-matrix models reveal how quantum confinement, dephasing, and Coulomb interactions govern the nonlinear optical responses in quantum dots.

High-order harmonic generation (HHG) in quantum dots (QDs) represents an intersection of strong-field nanophysics, ultrafast optics, and mesoscale quantum materials science. HHG is the process by which intense laser fields interacting with matter drive the emission of coherent radiation at frequencies that are integer multiples (harmonics) of the fundamental frequency. In QDs, the interplay between spatial confinement, discrete quantum structure, and many-body interactions produces nonlinear optical responses distinct from those observed in both atoms/molecules and extended solids. Recent experimental and theoretical work has elucidated the critical dependencies of the HHG yield, cutoff energy, and spectral structure on QD size, shape, dimensionality, and material composition, as well as the characteristics of the driving field.

1. Physical Principles of HHG in Quantum Dots

HHG in quantum dots is governed by the interaction of a strong optical field with carriers confined in a nanometer-scale potential well. The essential mechanisms involve:

• Nonlinear polarization response: The total nonlinear polarization can be decomposed into interband (transition between quantized valence and conduction states) and intraband (coherent motion within a band) components,

$$P_\mathrm{NL}(t) = P_\mathrm{inter}(t) + P_\mathrm{intra}(t)$$

where $$P_\mathrm{inter}(t) = \sum_{v,c} d_{vc} \rho_{cv}(t)$$ and $P_\mathrm{intra}(t) = \sum_n d_{nn} \rho_{nn}(t)$. The relative strength of these channels is modulated by the QD size and material parameters (Nakagawa et al., 2022).

• Modification of classical strong-field concepts: The standard gas-phase recollision cutoff law ($E_\mathrm{cutoff} = I_p + 3.17 U_p$, $U_p$ the ponderomotive energy) is fundamentally altered in QDs. When the QD diameter, $d$, falls below the classical electron quiver (oscillatory) radius $r_\mathrm{osc} = eE_0/m_e\omega^2$, boundary scattering truncates electron trajectories, reduces attainable kinetic energy, and suppresses high-energy (above-gap) harmonics (Gopalakrishna et al., 2022, Thümmler et al., 1 Dec 2025).
• Quantum confinement and level spacing: For spherical semiconductor QDs, quantized energy level spacing and subband gaps scale as $\sim 1/R^2$ where $R = d/2$ is the dot radius. As the dot size decreases, both the energy gap and the suppression of intraband dipole elements (scaling $\sim R$) play central roles in determining HHG efficiency (Nakagawa et al., 2022).
• Role of dephasing: Wall-induced scattering not only limits carrier acceleration but introduces phase randomization and chaotization, resulting in dephasing and suppressed interband coherence necessary for HHG at high orders (Gopalakrishna et al., 2022).

2. Experimental Methodologies and Laser-QD Interaction Regimes

The paper and exploitation of HHG in QDs require precise control of excitation conditions and QD properties:

• Sample systems: Experimental efforts predominantly employ colloidal II-VI semiconductor QDs (e.g., CdSe, CdS) synthesized with diameters $d=1.8$–$8.2$ nm, which are optically active well below the bulk exciton Bohr radius ($d<11$ nm for CdSe) (Gopalakrishna et al., 2022, Nakagawa et al., 2022).
• Laser systems: Mid-infrared (mid-IR) pulses with tunable wavelengths ($\lambda=3.0$–$5.0$ μm), durations in the 70–100 fs range, and intensities up to $\sim 1$ TW/cm$^2$ are used, entering the strong-field, non-perturbative regime (Gopalakrishna et al., 2022, Nakagawa et al., 2022, Thümmler et al., 1 Dec 2025).
• Harmonic detection: Transmitted harmonics (typically 200–1100 nm) are resolved via spectrometers and cooled CCDs. For orientation-averaged measurements, films are probed with random polarization states (Thümmler et al., 1 Dec 2025).
• HHG in laser-induced plasmas: QDs embedded in two-dimensional matrices (e.g., MoS$_2$ decorated with CdSe or CdSe/V$_2$O$_5$) can also be interrogated in plasma plumes created by femtosecond or picosecond heating pulses; this configuration can substantially raise both HHG yield and cutoff (Konda et al., 23 Apr 2024).

3. Size and Wavelength Scaling of HHG: Suppression, Enhancement, and Design Rules

The dependence of HHG efficiency and cutoff on QD size and laser parameters is nontrivial and sharply demarcated:

• Suppression regime: For $d<r_\mathrm{osc}$, especially when $d \lesssim 2$–3 nm and $\lambda \gtrsim 3.5\,\mu$m, above-gap harmonics ($E> E_g$) are suppressed by factors of $10^2$ or more relative to the bulk, with almost complete quenching for the smallest QDs (Gopalakrishna et al., 2022, Thümmler et al., 1 Dec 2025). The ratio $S(d) = Y_\mathrm{above}(d)/Y_\mathrm{below}(d)$ collapses from $\sim 1$ (bulk) to $\lesssim 0.01$ for $d=2.2$ nm at $\lambda=3.72\,\mu$m.
• Scaling law and optimization: The $N$-th harmonic yield scales as $$I_N(R) \propto |d_\mathrm{intra}(R) E_0|^N f[E_g(R)]$$, with $d_\mathrm{intra} \propto R$ and $E_g(R)\propto 1/R^2$. This means larger QDs favor increased intraband acceleration and thus higher yields, with $3$–$5$ nm delivering maximal intensity for CdSe and CdS QDs (Nakagawa et al., 2022).
• Wavelength dependence: Longer $\lambda$ increases $r_\mathrm{osc}$ (quiver radius), further constricting the effective phase space for high-order emission in small dots. Above-gap suppression becomes even more severe at $\lambda = 4.74\,\mu$m compared to shorter driving wavelengths (Gopalakrishna et al., 2022, Thümmler et al., 1 Dec 2025).
• Ellipticity dependence: The harmonic yield as a function of driving ellipticity is approximately Gaussian. The Full Width at Half Maxima (FWHM) narrows for higher harmonics and longer wavelengths, a signature shared with extended solids (Thümmler et al., 1 Dec 2025).
• Wall scattering and dephasing: Classical and quantum (TDSE, TDDFT) simulations confirm that electrons in small QDs experience strong boundary-induced dephasing, leading to smeared kinetic energy oscillations and quenching of the coherent recombination essential for HHG (Gopalakrishna et al., 2022).

4. Theoretical Models and Computational Approaches

Multiple theoretical formalisms underpin the interpretation and quantitative modeling of HHG in QDs:

• Real-space tight-binding models: Three-dimensional real-space approaches employing parameters derived from density functional theory and maximally localized Wannierization now permit simulations of QDs up to $\sim$3.3 nm (hundreds to thousands of atoms). These models reproduce experimental trends—specifically, the onset of above-gap HHG only for $d \gtrsim 2$ nm and the saturation to bulk yield at $d \gtrsim 3$ nm (Thümmler et al., 1 Dec 2025).
• Density-matrix dynamics: The time evolution of the single-particle density matrix—propagated under a tight-binding Hamiltonian that includes the laser interaction in the length gauge—captures both interband transitions and boundary-induced dephasing, yielding directly the emitted HHG spectrum (Gopalakrishna et al., 2022, Nakagawa et al., 2022).
• Classical trajectory and Monte Carlo methods: These are used to visualize and quantify the effect of wall scattering, with phase randomization at each boundary reflection dampening the periodicity of kinetic energy and hence suppressing above-gap emission (Gopalakrishna et al., 2022).
• Analytical scaling and cutoff behavior: Unlike the quadratic $E_0^2$ dependence ($E_\mathrm{cutoff}^{\rm atom}=I_p+3.17U_p$) in atoms, HHG cutoff energy in graphene quantum dots shows a characteristic $\sqrt{E_0}$ scaling owing to discrete-level structure and long-range correlations (Avetissian et al., 2023), while in discrete-level semiconductor QDs, scaling with QD size and field strength reflects confinement-imposed limits (Ghazaryan et al., 2022, Avchyan et al., 2021).
• Role of Coulomb interactions: Inclusion of long-range Coulomb terms (Hartree–Fock or extended Hubbard) in graphene and semiconductor QDs enhances both the intensity and cutoff of HHG, with pronounced fine structure in the spectra and coherence-dependent multi-plateau formation (Avetissian et al., 2023, Ghazaryan et al., 2022).

5. Material Systems: From Semiconductors to Graphene QDs

Comparative studies across material and geometric classes reveal nuanced effects:

• Semiconductor QDs (e.g., CdSe, CdS): Exhibit sharply tunable quantum confinement; the balance of interband and intraband channels is modulated by size-dependent subband spacing and effective mass (Nakagawa et al., 2022). Surface roughness, ligand shell disorder, and matrix embedding further control dephasing and yield (Gopalakrishna et al., 2022).
• Graphene quantum dots (GQDs): Both rectangular and triangular GQDs display unique nonlinear responses determined by edge geometry (zigzag vs armchair), lateral size, and symmetry. Zigzag-edged GQDs exhibit simultaneous odd and even harmonics due to absence of inversion symmetry, while armchair GQDs emit only odd orders. In GQDs, strongly-correlated dynamics and virtual bound-bound transitions dominate over recollision, with cutoff energy scaling as $\sqrt{E_0}$ (Ghazaryan et al., 2022, Avetissian et al., 2023, Avchyan et al., 2021).
• QD-decorated 2D systems and nanoplasmas: Incorporation of bare or core/shell QDs (e.g., CdSe/V$_2$O$_5$) into 2D materials (e.g., MoS$_2$) modifies the effective plasma properties upon laser ablation. This configuration yields enhanced harmonic intensity (by $1.5\times$–$2.5\times$) and extended cutoff (up to eight harmonic orders higher), attributable to both increased local field strengths and reduced composite ionization potential (Konda et al., 23 Apr 2024).

6. Outlook and Design Guidelines for HHG Optimization in QDs

Synthesis control and theoretical advances provide generalizable design principles:

• QD size selection: To maximize HHG yield and cutoff, the QD diameter should exceed the classical quiver radius for chosen laser parameters ($d \gtrsim r_\mathrm{osc}$). For typical mid-IR fields, this implies $d \gtrsim 3$–$5$ nm for II-VI QDs (Gopalakrishna et al., 2022, Thümmler et al., 1 Dec 2025).
• Tuning intraband resonance: For efficient coupling of intraband and interband pathways, design QDs where subband spacing $\Delta E_\mathrm{sub}$ approaches the driving photon energy $\hbar\omega_0$—generally reached at $R\gtrsim1.5$ nm for CdSe-type materials (Nakagawa et al., 2022).
• Mitigating dephasing: Minimize boundary roughness and surface ligand disorder to reduce extrinsic dephasing. Use shorter driving wavelengths or engineered local field enhancements (e.g., with plasmonic or dielectric antennas) at fixed QD size to decrease $r_\mathrm{osc}$ and favor bulk-like HHG behavior (Gopalakrishna et al., 2022).
• Material and shape engineering: In GQDs, lateral size and edge orientation (zigzag vs armchair) allow for further control of cutoff energy, angular emission profile, and harmonic content. Edge-state engineering permits the tailoring of selection rules (e.g., for even/odd orders) (Ghazaryan et al., 2022, Avchyan et al., 2021).
• Theoretical frameworks: Continued development of efficient real-space tight-binding and density-matrix models, potentially including explicit many-body Coulomb and phonon coupling, is essential for predictive design across the QD/mesoscopic solid-state regime (Thümmler et al., 1 Dec 2025).

7. Comparative Table: Size Dependence of Above-Gap HHG Yield

The following summarizes key size-scaling behaviors for CdSe QDs (at $\lambda=3.72$ μm), concentrating on the normalized ratio $S(d)=Y_{\rm above}(d)/Y_{\rm below}(d)$ (Gopalakrishna et al., 2022):

QD diameter $d$ (nm)	$S(d)$ at 3.72 μm	Qualitative Regime
Bulk (~140 nm)	$\sim 1$	Bulk-like, full HHG plateau
4.2	$\sim 0.8$	Weak spatial confinement
3.0	$\sim 0.2$	Onset significant suppression
2.2	$\lesssim 0.01$	Nearly complete above-gap quench

For $\lambda=4.74$ μm, the suppression becomes even more pronounced, with negligible above-gap harmonic emission for $d<3$ nm. This table highlights the critical size threshold demarcating efficient and quenched HHG regimes due to spatial confinement.

---

The emerging understanding of HHG in quantum dots positions these systems as tunable platforms for nonlinear nanophotonics, extreme ultraviolet (XUV) light sources, and probes of strong-field and many-body quantum dynamics at the mesoscale. The interplay of quantum confinement, boundary effects, and carrier coherence in QDs reveals a new regime—a “nanosolid” intermediate between atomic/molecular and bulk solid-state HHG response (Gopalakrishna et al., 2022, Thümmler et al., 1 Dec 2025).