We present a formalism based on nonequilibrium dynamical mean-field theory (DMFT) which allows to compute the time-resolved x-ray absorption spectrum (XAS) of photoexcited solids. By applying this formalism to the photodoped half-filled and quarter-filled two-orbital Hubbard models in the Mott insulating regime we clarify how the time-resolved XAS signal reflects the nonequilibrium population of different local states. Apart from the missing broadening associated with continuum excitations, the atomic XAS spectrum computed with the nonthermal state populations provides a good approximation to the full nonequilibrium DMFT result. This suggests a route to combine the accurate DMFT description of nonequilibrium states of solids with cluster calculations of the XAS signal.
We present a formalism based on nonequilibrium dynamical mean-field theory (DMFT) which allows to compute the time-resolved x-ray absorption spectrum (XAS) of photoexcited solids. By applying this formalism to the photodoped half-filled and quarter-filled two-orbital Hubbard models in the Mott insulating regime we clarify how the time-resolved XAS signal reflects the nonequilibrium population of different local states. Apart from the missing broadening associated with continuum excitations, the atomic XAS spectrum computed with the nonthermal state populations provides a good approximation to the full nonequilibrium DMFT result. This suggests a route to combine the accurate DMFT description of nonequilibrium states of solids with cluster calculations of the XAS signal.
Recent seminal experiments have utilized time-resolved x-ray absorption spectroscopy (XAS) to investigate the ultrafast photoinduced renormalization of the electron interaction (“Hubbard U”) in Mott and charge transfer insulators. In this paper, we theoretically analyze the change of interactions due to dynamical screening as it is encoded in the XAS signal, using the nonequilibrium GW+EDMFT formalism. Our study shows that XAS is well-suited for measuring this change, but two aspects should be kept in mind if the screening processes are not substantially faster than the valence electron dynamics: (i) Screening in a photoexcited system can affect both the position and the lineshape of the absorption lines. (ii) In general, the effect cannot be captured by the modification of a single interaction parameter. Specifically, an estimate for Δ� extracted from the shift of the XAS lines does not necessarily describe the related shift of the upper Hubbard band. We clarify these aspects using a minimal cluster model and the three-band Emery model for a charge transfer insulator.
Abstract Resonant inelastic X-ray scattering (RIXS) detects various types of high- and low-energy elementary excitations in correlated solids, and this tool will play an increasingly important role in investigations of time-dependent phenomena in photo-excited systems. While theoretical frameworks for the computation of equilibrium RIXS spectra are well established, the development of appropriate methods for nonequilibrium simulations is an active research field. Here, we apply a recently developed nonequilibrium dynamical mean field theory (DMFT) based approach to compute the RIXS response of photo-excited two-orbital Mott insulators. The results demonstrate the feasibility of multi-orbital nonequilibrium RIXS calculations and the sensitivity of the quasi-elastic fluorescence-like features and excitation peaks on the nonequilibrium population of the Hubbard bands.
Raman spectroscopy uses light scattering to extract information on low-energy excitations of solids. The Raman process is described by diagrams which are fourth order in the light-matter interaction and in particular the resonant contribution, which involves four different space-time arguments, is difficult to evaluate. If one instead simulates explicitly the incoming (classical) light pulse, the Raman signal is given by the outgoing photon flux and can be determined from a two-point correlation function. Such a formalism can be used to compute the time-resolved Raman spectrum of nonequilibrium systems, as well as nonlinear signals which are higher order in the incoming field, such as hyper-Raman scattering. Here we explain how to implement this time-dependent formalism within the dynamical mean field theory framework. The method is illustrated with applications to the Holstein-Hubbard model in the strong electron-phonon coupling regime. We demonstrate hyper-Raman scattering in measurements with strong probe fields and frequency mixing signals in the presence of a pump field and simulate the evolution of Stokes and anti-Stokes features after photoexcitations of metallic and Mott insulating systems.
The study of nonequilibrium phenomena in correlated lattice systems has developed into one of the most active and exciting branches of condensed matter physics. This research field provides rich new insights that could not be obtained from the study of equilibrium situations, and the theoretical understanding of the physics often requires the development of new concepts and methods. On the experimental side, ultrafast pump-probe spectroscopies enable studies of excitation and relaxation phenomena in correlated electron systems, while ultracold atoms in optical lattices provide a new way to control and measure the time evolution of interacting lattice systems with a vastly different characteristic time scale compared to electron systems. A theoretical description of these phenomena is challenging because, first, the quantum-mechanical time evolution of many-body systems out of equilibrium must be computed and second, strong-correlation effects which can be of a nonperturbative nature must be addressed. This review discusses the nonequilibrium extension of the dynamical mean field theory (DMFT), which treats quantum fluctuations in the time domain and works directly in the thermodynamic limit. The method reduces the complexity of the calculation via a mapping to a self-consistent impurity problem, which becomes exact in infinite dimensions. Particular emphasis is placed on a detailed derivation of the formalism, and on a discussion of numerical techniques, which enable solutions of the effective nonequilibrium DMFT impurity problem. Insights gained into the properties of the infinite-dimensional Hubbard model under strong nonequilibrium conditions are summarized. These examples illustrate the current ability of the theoretical framework to reproduce and understand fundamental nonequilibrium phenomena, such as the dielectric breakdown of Mott insulators, photodoping, and collapse-and-revival oscillations in quenched systems. Furthermore, remarkable novel phenomena have been predicted by the nonequilibrium DMFT simulations of correlated lattice systems, including dynamical phase transitions and field-induced repulsion-to-attraction conversions.
We study the relaxation dynamics of photocarriers in the paramagnetic Mott insulating phase of the half-filled two-band Hubbard model. Using nonequilibrium dynamical mean-field theory, we excite charge carriers across the Mott gap by a short hopping modulation, and simulate the evolution of the photodoped population within the Hubbard bands. We observe an ultrafast charge-carrier relaxation driven by the emission of local spin excitations with an inverse relaxation time proportional to the Hund's coupling. The photodoping generates additional side-bands in the spectral function, and for strong Hund's coupling, the photodoped population also splits into several resonances. The dynamics of the local many-body states reveals two effects, thermal blocking and kinetic freezing, which manifest themselves when the Hund's coupling becomes of the order of the temperature or the bandwidth, respectively. These effects, which are absent in the single-band Hubbard model, should be relevant for the interpretation of experiments on correlated materials with multiple active orbitals. In particular, the features revealed in the nonequilibrium energy distribution of the photocarriers are experimentally accessible, and provide information on the role of the Hund's coupling in these materials.
Photodoped states are widely observed in laser-excited Mott insulators, in which charge excitations are quickly created and can exist beyond the duration of the external driving. Despite the fruitful experimental explorations, theoretical studies on the microscopic models face the challenge to simultaneously deal with exponentially separated time scales, especially in multiband systems, where the longtime behaviors are often well beyond the reach of state-of-the-art numerical tools. Here, we address this difficulty by introducing a steady-state description of photodoped Mott insulators using an open-system setup, where the photodoped system is stabilized as a nonequilibrium steady state by a weak external driving. Taking advantage of the stationarity, we implement and discuss the details of an efficient numerical tool using the steady-state dynamical mean-field theory combined with the noncrossing approximation. We demonstrate that these stationary photodoped states exhibit the same properties of their transient counterparts, while being solvable with reasonable computational efforts. Furthermore, they can be parametrized by just a few physical quantities, including the effective temperature and the density of charge excitations, which confirms the universal nature of photodoped states indeed independent of the excitation protocols. As a first application, we consider the stationary photodoped states in a two-band Hubbard model with intertwined spin-and-orbital ordering and find a family of hidden phases unknown from the previous studies, implying an apparently unexplored time regime of the relaxation of the intertwined orders.