No Arabic abstract
The high harmonic spectrum of the Mott insulating Hubbard model has recently been shown to exhibit plateau structures with cutoff energies determined by $n$th nearest neighbor doublon-holon recombination processes. The spectrum thus allows to extract the on-site repulsion $U$. Here, we consider generalizations of the single-band Hubbard model and discuss the signatures of bosonic excitations in high harmonic spectra. Specifically, we study an electron-plasmon model which captures the essential aspects of the dynamically screened Coulomb interaction in solids and a multi-orbital Hubbard model with Hund coupling which allows to analyze the effect of local spin excitations. For the electron-plasmon model, we show that the high harmonic spectrum can reveal information about the screened and bare onsite interaction, the boson frequency, as well as the relation between boson coupling strength and boson frequency. In the multi-orbital case, string states formed by local spin excitations result in an increase of the radiation intensity and cutoff energy associated with higher order recombination processes.
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 $epsilon_{rm 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 $m$-th plateau originates from the recombination of $m$-th 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, $epsilon_{rm cut} = Delta_{rm gap} + alpha E_0$, with $Delta_{rm gap}$ the band gap and $alpha>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.
With a combination of numerical methods, including quantum Monte Carlo, exact diagonalization, and a simplified dynamical mean-field model, we consider the attosecond charge dynamics of electrons induced by strong-field laser pulses in two-dimensional Mott insulators. The necessity to go beyond single-particle approaches in these strongly correlated systems has made the simulation of two-dimensional extended materials challenging, and we contrast their resulting high-harmonic emission with more widely studied one-dimensional analogues. As well as considering the photo-induced breakdown of the Mott insulating state and magnetic order, we also resolve the time and ultra-high frequency domains of emission, which are used to characterize both the photo-transition, and the sub-cycle structure of the electron dynamics. This extends simulation capabilities and understanding of the photo-melting of these Mott insulators in two-dimensions, at the frontier of attosecond non-equilibrium science of correlated materials.
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 $epsilon_{rm cut}$ scales linearly with the field strength $E_0$ as $epsilon_{rm cut}=Delta_{rm gap} + alpha E_0$, where $Delta_{rm 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 $epsilon_{rm cut,m}=U + m E_0$. 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.
The density-matrix renormalization group (DMRG) is employed to calculate optical properties of the half-filled Hubbard model with nearest-neighbor interactions. In order to model the optical excitations of oligoenes, a Peierls dimerization is included whose strength for the single bonds may fluctuate. Systems with up to 100 electrons are investigated, their wave functions are analyzed, and relevant length-scales for the low-lying optical excitations are identified. The presented approach provides a concise picture for the size dependence of the optical absorption in oligoenes.
Angle-resolved photoemission spectroscopy (ARPES) has revealed peculiar properties of mobile dopants in correlated anti-ferromagnets (AFMs). But describing them theoretically, even in simplified toy models, remains a challenge. Here we study ARPES spectra of a single mobile hole in the $t-J$ model. Recent progress in the microscopic description of mobile dopants allows us to use a geometric decoupling of spin and charge fluctuations at strong couplings, from which we conjecture a one-to-one relation of the one-dopant spectral function and the spectrum of a constituting spinon in the emph{undoped} parent AFM. We thoroughly test this hypothesis for a single hole doped into a 2D Heisenberg AFM by comparing our semi-analytical predictions to previous quantum Monte Carlo results and our large-scale time-dependent matrix product state (td-MPS) calculations of the spectral function. Our conclusion is supported by a microscopic trial wavefuntion describing spinon-chargon bound states, which captures the momentum and $t/J$ dependence of the quasiparticle residue. Our conjecture suggests that ARPES measurements in the pseudogap phase of cuprates can directly reveal the Dirac-fermion nature of the constituting spinons. Specifically, we demonstrate that our trial wavefunction provides a microscopic explanation for the sudden drop of spectral weight around the nodal point associated with the formation of Fermi arcs, assuming that additional frustration suppresses long-range AFM ordering. We benchmark our results by studying the cross-over from two to one dimension, where spinons and chargons are confined and deconfined respectively.