We investigate the potential of a quantum Boltzmann equation without momentum conservation for description of strongly correlated electron systems out of equilibrium. In a spirit similar to dynamical mean field theory (DMFT), the momentum conservation of the electron-electron scattering is neglected, which yields a time-dependent occupation function for the equilibrium spectral function, even in cases where well-defined quasiparticles do not exist. The main assumption of this method is that the spectral function remains sufficiently rigid under the nonequilibrium evolution. We compare the result of the quantum Boltzmann equation to nonequilibrium DMFT simulations for the case of photocarrier relaxation in Mott insulators, where processes on very different timescales emerge, i.e., impact ionization, intra-Hubbard-band thermalization, and full thermalization. Since quantum Boltzmann simulations without momentum conservation are computationally cheaper than nonequilibrium DMFT, this method allows the simulation of more complicated systems or devices, and to access much longer times.
We investigate the potential of a quantum Boltzmann equation without momentum conservation for description of strongly correlated electron systems out of equilibrium. In a spirit similar to dynamical mean field theory (DMFT), the momentum conservation of the electron-electron scattering is neglected, which yields a time-dependent occupation function for the equilibrium spectral function, even in cases where well-defined quasiparticles do not exist. The main assumption of this method is that the spectral function remains sufficiently rigid under the nonequilibrium evolution. We compare the result of the quantum Boltzmann equation to nonequilibrium DMFT simulations for the case of photocarrier relaxation in Mott insulators, where processes on very different timescales emerge, i.e., impact ionization, intra-Hubbard-band thermalization, and full thermalization. Since quantum Boltzmann simulations without momentum conservation are computationally cheaper than nonequilibrium DMFT, this method allows the simulation of more complicated systems or devices, and to access much longer times.
We investigate the potential of a quantum Boltzmann equation without momentum conservation, i.e., a kinetic equation for the local distribution function, for describing correlated electron systems out of equilibrium. The equation can be derived either for the dynamics of well-defined quasi- particles, neglecting the momentum conservation of the scattering in the spirit of dynamical mean field theory (DMFT), or, equivalently, as the quantum Boltzmann equation for the local occupied density of states in non-equilibrium DMFT. The approach may thus be applied even in cases where well-defined quasiparticles do not exist, with the main assumption that the spectral function remains rigid under the non-equilibrium evolution. We compare the result of the Boltzmann equation to non-equilibrium DMFT for the case of photo-carrier relaxation in Mott insulators, where processes on very different timescales emerge, i.e., impact ionization, intra-Hubbard-band thermalization, and full thermalization. Since the Boltzmann without momentum conservation is much simpler, it may allow us to simulate more complicated systems, whole devices, and much longer times.
Collective orders and photoinduced phase transitions in quantum matter can evolve on timescales which are orders of magnitude slower than the femtosecond processes related to electronic motion in the solid. Quantum Boltzmann equations can potentially resolve this separation of timescales, but are often constructed by assuming the existence of quasiparticles. Here we derive a quantum Boltzmann equation which only assumes a separation of timescales (taken into account through the gradient approximation for convolutions in time), but is based on a nonperturbative scattering integral, and makes no assumption on the spectral function such as the quasiparticle approximation. In particular, a scattering integral corresponding to nonequilibrium dynamical mean-field theory is evaluated in terms of an Anderson impurity model in a nonequilibrium steady state with prescribed distribution functions. This opens the possibility to investigate dynamical processes in correlated solids with quantum impurity solvers designed for the study of nonequilibrium steady states.
We compute scattering rates for electrons in the two-dimensional Hubbard model for a one-orbital metal and Mott insulator and a two-orbital band insulator by means of dynamical mean-field theory (DMFT) and the Boltzmann scattering equation (BSE), which for the considered equilibrium scattering rates is equivalent to second order perturbation theory. As an intermediate method between both, we also consider the BSE without momentum conservation. In the weak interaction regime and for the band insulator, the last two agree to very good accuracy. The BSE with momentum conservation, on the other hand, shows slightly larger scattering rates, and a momentum differentiation of these on the Fermi surface. For the Mott insulator at strong interaction, the DMFT electron scattering rates are much larger and defy a BSE description. Noteworthy, the scattering rates for the band insulator are exceedingly small because---in contrast to the Mott insulator---impact ionization is not mediated by strong (density-density) Coulomb matrix elements.
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.
This work is motivated by the fact that the investigation of non-equilibrium phenomena in strongly correlated electron systems has developed into one of the most active and exciting branches of condensed matter physics as it provides rich new insights that could not be obtained from the study of equilibrium situations. However, a theoretical description of those phenomena is missing. Therefore, in this thesis, we develop a numerical method that can be used to study two minimal models -- the Hubbard model and the Anderson impurity model with general parameter range and time dependence. We begin by introducing the theoretical framework and the general features of the Hubbard model. We then describe the dynamical mean field theory (DMFT), which was first invented by Georges in 1992. It provides a feasible way to approach strongly correlated electron systems and reduces the complexity of the calculations via a mapping of lattice models onto quantum impurity models subject to a self-consistency condition. We employ the non-equilibrium extension of DMFT and map the Hubbard model to the single impurity Anderson model (SIAM). Since the fundamental component of the DMFT method is a solver of the single impurity Anderson model, we continue with a description of the formalism to study the real-time dynamics of the impurity model staring at its thermal equilibrium state. We utilize the non-equilibrium strong-coupling perturbation theory and derive semi-analytical approximation methods such as the non-crossing approximation (NCA) and the one-crossing approximation (OCA). We then use the Quantum Monte-Carlo method (QMC) as a numerically exact method and present proper measurements of local observables, current and Green's functions. We perform simulations of the current after a quantum quench from equilibrium by rapidly applying a bias voltage in a wide range of initial temperatures. The current exhibits short equilibrium times and saturates upon the decrease of temperature at all times, indicating Kondo behavior both in the transient regime and in the steady state. However, this bare QMC solver suffers from a dynamical sign problem for long time propagations. To overcome the limitations of this bare treatment, we introduce the ``Inchworm algorithm'', based on iteratively reusing the information obtained in previous steps to extend the propagation to longer times and stabilize the calculations. We show that this algorithm greatly reduces the required order for each simulation and re-scales the exponential challenge to quadratic in time. We introduce a method to compute Green's functions, spectral functions, and currents for inchworm Monte Carlo and show how systematic error assessments in real time can be obtained. We illustrate the capabilities of the algorithm with a study of the behavior of quantum impurities after an instantaneous voltage quench from a thermal equilibrium state. We conclude with the applications of the unbiased inchworm impurity solver to DMFT calculations. We employ the methods for a study of the one-band paramagnetic Hubbard model on the Bethe lattice in equilibrium, where the DMFT approximation becomes exact. We begin with a brief introduction of the Mott metal insulator phase diagram. We present the results of both real time Green's functions and spectral functions from our nonequilibrium calculations. We observe the metal-insulator crossover as the on-site interaction is increased and the formation of a quasi-particle peak as the temperature is lowered. We also illustrate the convergence of our algorithms in different aspects.
We present a review of our recent work in extending the successful dynamical mean-field theory from the equilibrium case to nonequilibrium cases. In particular, we focus on the problem of turning on a spatially uniform, but possibly time varying, electric field (neglecting all magnetic field effects). We show how to work with a manifestly gauge-invariant formalism, and compare numerical calculations from a transient-response formalism to different types of approximate treatments, including the semiclassical Boltzmann equation and perturbation theory in the interaction. In this review, we solve the nonequilibrium problem for the Falicov-Kimball model, which is the simplest many-body model and the easiest problem to illustrate the nonequilibrium behavior in both diffusive metals and Mott insulators. Due to space restrictions, we assume the reader already has some familiarity both with the Kadanoff-Baym-Keldysh nonequilibrium formalism and with equilibrium dynamical mean-field theory; we provide a guide to the literature where additional details can be found.