The many-body formalism for dynamical mean-field theory is extended to treat nonequilibrium problems. We illustrate how the formalism works by examining the transient decay of the oscillating current that is driven by a large electric field turned on at time t=0. We show how the Bloch oscillations are quenched by the electron-electron interactions, and how their character changes dramatically for a Mott insulator.
The many-body formalism for dynamical mean-field theory is extended to treat nonequilibrium problems. We illustrate how the formalism works by examining the transient decay of the oscillating current that is driven by a large electric field turned on at time t=0. We show how the Bloch oscillations are quenched by the electron-electron interactions, and how their character changes dramatically for a Mott insulator.
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.
Dynamical mean-field theory is generalized to solve the nonequilibrium Keldysh boundary problem: a system is started in equilibrium at a temperature T=0.1, a uniform electric field is turned on at t=0, and the system is monitored as it approaches the steady state. The focus here is on the Bloch oscillations of the current and how they decay after their initial appearance near t=0. The system is evolved out to the longest time allowed by our computational resources---in most cases we are unable to reach the steady state. The strongly correlated material is modeled by the spinless Falicov-Kimball model at half-filling on a hypercubic lattice in $d=\ensuremath{\infty}$ dimensions, which has a Mott-like metal-insulator transition at U=2. The computational algorithm employed is highly efficient, parallelizes well, and scales to thousands of processors. For strong fields, we find beats develop with a period of $2\ensuremath{\pi}∕U$, while for strong interactions, the Bloch oscillations are sharply damped and become quite irregular in time.
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 effect of spatially nonlocal correlations on the nonequilibrium dynamics of interacting fermions by constructing the nonequilibrium dynamical cluster theory, a cluster generalization of the nonequilibrium dynamical mean-field theory (DMFT). The formalism is applied to interaction quenches in the Hubbard model in one and two dimensions, and the results are compared with data from single-site DMFT, the time-dependent density matrix renormalization group, and lattice perturbation theory. Both in one and two dimensions the double occupancy quickly thermalizes, while the momentum distribution relaxes only on much longer time scales. For the two-dimensional square lattice we find a strongly momentum-dependent evolution of the momentum distribution around the Fermi energy, with a much faster relaxation near the momenta $(0,\ensuremath{\pi})$ and $(\ensuremath{\pi},0)$ than near $(\ensuremath{\pi}/2,\ensuremath{\pi}/2)$. This result is interpreted as reflecting the momentum-anisotropic quasiparticle lifetime of the marginal Fermi liquid. The method is further applied to the two-dimensional Hubbard model driven by a dc electric field, where the damping of the Bloch oscillation of the current is found to be less effective than predicted by DMFT and lattice perturbation theory.
Understanding strongly correlated systems driven out of equilibrium is a challenging task necessitating the simultaneous treatment of quantum mechanics, dynamical constraints, and strong interactions. A Mott insulator subjected to a uniform and static electric field is prototypical, raising key questions such as the fate of Bloch oscillations with increasing correlation strength, the approach to a steady-state DC transport regime and the role of dissipation in it, and electric-field-driven phase transitions. Despite tremendous efforts over the past decade employing various numerical and analytical approaches, the manner in which a nonequilibrium steady state gets established has remained an unresolved problem. We develop here an effective large-N Keldysh field theory for studying nonequilibrium transport in a regular one-dimensional dissipative Mott insulator system subjected to a uniform electric field. Upon abruptly turning on the electric field (a quench), a transient oscillatory current response reminiscent of Bloch oscillations is found. In the regime of small tunneling conductance the amplitude of these oscillations, over a large time window, decreases as an inverse square power law in time, ultimately going over to an exponential decay beyond a large characteristic time ${\ensuremath{\tau}}_{d}$ that increases with N. Such a relaxation to a steady-state DC response is absent in the dissipation-free Hubbard chain at half filling. The steady-state current at small fields is governed by large-distance cotunneling, a process absent in the equilibrium counterpart. The low-field DC current has a Landau-Zener-Schwinger form but qualitatively differs from the expression for pair-production probability for the dissipation-free counterpart. The breakdown of perturbation theory in the Mott phase possibly signals a nonequilibrium phase transition to a metallic phase. Our study sheds light on the approach of a driven, dissipative strongly correlated system to a nonequilibrium steady state and also provides a general analytic microscopic framework for understanding other nonequilibrium phenomena in these systems.
Using nonequilibrium dynamical mean-field theory, we compute the time evolution of the current in a Mott insulator after a strong electric field is turned on. We observe the formation of a quasistationary state in which the current is almost time independent although the system is constantly excited. At moderately strong fields this state is stable for quite long times. The stationary current exhibits a threshold behavior as a function of the field, in which the threshold increases with the Coulomb interaction and vanishes as the metal-insulator transition is approached.