We study the high harmonic generation (HHG) in Mott insulators using Floquet dynamical mean-field theory (DMFT). We show that the main origin of the HHG in Mott insulators is the doublon-holon recombination, and that the character of the HHG spectrum differs depending on the field strength. In the weaker-field regime, the HHG spectrum shows a single plateau as in the HHG from gases, and its cut-off energy
εcut scales linearly with the field strength
E0 as
εcut=Δgap+αE0, where
Δgap is the Mott gap. On the other hand, in the stronger-field regime, multiple plateaus emerge and the
m-th cut-off scales as
εcut,m=U+mE0. We show that this difference originates from the different dynamics of the doublons and holons in the weak- and strong-field regimes. We also comment on the similarities and differences between HHG from Mott insulators and from semiconductors. This proceedings paper complements our recent work, Phys. Rev. Lett. 121, 057405 (2018), with additional results and analyses.
Using Floquet dynamical mean-field theory, we study the high-harmonic generation in the time-periodic steady states of wide-gap Mott insulators under ac driving. In the strong-field regime, the harmonic intensity exhibits multiple plateaus, whose cutoff energies ε_{cut}=U+mE_{0} scale with the Coulomb interaction U and the maximum field strength E_{0}. In this regime, the created doublons and holons are localized because of the strong field and the mth plateau originates from the recombination of mth nearest-neighbor doublon-holon pairs. In the weak-field regime, there is only a single plateau in the intensity, which originates from the recombination of itinerant doublons and holons. Here, ε_{cut}=Δ_{gap}+αE_{0}, with Δ_{gap} the band gap and α>1. We demonstrate that the Mott insulator shows a stronger high-harmonic intensity than a semiconductor model with the same dispersion as the Mott insulator, even if the semiconductor bands are broadened by impurity scattering to mimic the incoherent scattering in the Mott insulator.
To study excitonic effects on high-harmonic generation (HHG) in Mott insulators, we investigate pumped nonequilibrium dynamics in the one-dimensional extended Hubbard model. By employing time-dependent calculations based on the exact diagonalization and infinite time-evolving block decimation methods, we find the strong enhancement of the HHG intensity around the exciton energy. The subcycle analysis in the sub-Mott-gap regime shows that the intensity region of the time-resolved spectrum around the exciton energy splits into two levels and oscillates following the driving electric field. This excitonic dynamics is qualitatively different from the dynamics of free doublon and holon but favorably contributes to HHG in the Mott insulator.