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.
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 review the dynamical mean-field theory of strongly correlated electron systems which is based on a mapping of lattice models onto quantum impurity models subject to a self-consistency condition. This mapping is exact for models of correlated electrons in the limit of large lattice coordination (or infinite spatial dimensions). It extends the standard mean-field construction from classical statistical mechanics to quantum problems. We discuss the physical ideas underlying this theory and its mathematical derivation. Various analytic and numerical techniques that have been developed recently in order to analyze and solve the dynamical mean-field equations are reviewed and compared to each other. The method can be used for the determination of phase diagrams (by comparing the stability of various types of long-range order), and the calculation of thermodynamic properties, one-particle Green's functions, and response functions. We review in detail the recent progress in understanding the Hubbard model and the Mott metal-insulator transition within this approach, including some comparison to experiments on three-dimensional transition-metal oxides. We present an overview of the rapidly developing field of applications of this method to other systems. The present limitations of the approach, and possible extensions of the formalism are finally discussed. Computer programs for the numerical implementation of this method are also provided with this article.
Since the first investigation of the Hubbard model in the limit of infinite dimensions by Metzner and Vollhardt, dynamical mean-field theory (DMFT) has become a very powerful tool for the investigation of lattice models of correlated electrons. In DMFT the lattice model is mapped on an effective quantum impurity model in a bath which has to be determined self-consistently. This approach led to significant progress in our understanding of typical correlation problems such as the Mott transition; furthermore, the combination of DMFT with ab-initio methods now allows for a realistic treatment of correlated materials. The focus of these lecture notes is on the relation between quantum impurity physics and the physics of lattice models within DMFT. Issues such as the observability of impurity quantum phase transitions in the corresponding lattice models are discussed in detail.
Quantum mechanical lattice models, such as the famous Hubbard model, were originally intended as simplified descriptions of the behavior of electrons in real solids. Today it is clear that even such simple models can present remarkably rich physics. At the same time, obtaining reliable solutions for correlated quantum systems is difficult and requires the development of both computing resources and new algorithms and approximations. Parallel with the theoretical developments, applications of lattice models have also expanded to the study of ultracold atoms in optical lattices. The hallmark of these experiments is the possibility to design different kinds of lattices and to tune the interactions between the atoms, making it possible to do real experiments on theoretical "toy models". \n \nIn this thesis we have used dynamical mean-field theory (DMFT) and its cluster extensions to study superfluidity and topological phenomena in lattice models. The DMFT combines the exact solution of a small reference system with a self-consistent mean-field treatment of the rest of the lattice. After a short introduction to the physical background of the lattice models, we review the formal definition and practical aspects of the DMFT method. We discuss the Baym-Kadanoff-Luttinger-Ward functional formalism that can be used to derive the DMFT equations, and review the numerical exact diagonalization and continuous-time quantum Monte-Carlo algorithms that were used in the implementation of the method. \n \nPublications I and VI discuss superfluid states induced by correlation effects that cannot easily be treated within a static mean-field theory. In Publication I we study a bilayer system where two different types of superfluids can be found depending on the densities of particles in the layers. Publication VI has a focus in the d-wave superconductivity and non-uniform, striped magnetism of the two-dimensional Hubbard model, which are thought to be connected with the high-Tc superconductivity in cuprate materials. The stripe state along with non-Fermi-liquid properties visible only in a beyond-mean-field treatment were observed also in the study of the Hubbard model on the Lieb lattice in Publication V. \n \nTopological and geometric properties of lattice systems, such as the Berry curvature and Chern numbers, play a central role in Publications II, III and IV. In Publication II we used exact diagonalization, DMFT and mean-field theory to study the topological phase diagram of an interacting Haldane model. The methods consistently predict a phase diagram with an interme-diate topological phase not found in the well-known non-interacting case. In Publications III and IV mean-field theory is used to study the connection between Berry curvature and the superfluid weight in lattice models, and DMFT calculations were used to support the mean-field results. Topological lattice systems have recently been realized in ultracold gas experiments and it is hoped that the calculations presented here inspire this line of research.
In this thesis, we describe methods for solving quantum impurity models out of equilibrium, and consider several applications of these methods. Quantum impurity models describe systems in which a small subsystem (the impurity) is coupled to a large external environment (the bath). These models are useful in describing a wide variety of physical systems, including quantum dots, absorbed adatoms, and molecular junctions. They are also important for the study of lattice models, such as the Hubbard model, where they appear as auxiliary problems within the framework of dynamical mean field theory (DMFT). Many interesting phenomena in these systems, such as transport through a quantum dot, or optical pumping of a Mott insulator, occur away from equilibrium. Theoretical description of these phenomena requires methods that are capable of handling both real time evolution and strong correlation effects presenting a formidable challenge. In this thesis, we describe several methods for the solution of non-equilibrium quantum impurity models which meet this challenge. These methods are based on an expansion in the coupling between the impurity and the bath known as the hybridization expansion. Chapter 1 gives an introduction to the problems addressed in this thesis. Chapter 2 provides a brief overview of Keldysh formalism for nonequilibrium problems. Chapter 3 describes the bare hybridization expansion for the Anderson impurity model as well as its bold extensions. In particular we derive the non-crossing approximation (NCA), the one-crossing approximation (OCA), the inchworm expansion, and a bold expansion for the impurity Green's function. Chapter 4 describes the implementation of a quantum Monte Carlo algorithm based on the inchworm expansion. The rest of the thesis describes applications of these algorithms to problems of physical interest. Chapters 5 and 6 apply the methods previously described to understanding the Kondo effect in quantum dot systems. The first chapter examines the behavior of the Kondo effect under an applied voltage bias, a classic non-equilibrium strongly-correlated problem. The second chapter examines different ways of measuring the cloud of electrons which screen the impurity spin in the Kondo regime. Chapters 7 and 8 apply the same methods to the Hubbard model using the Dynamical mean-field theory (DMFT) mapping. The first chapter studies the Hubbard model in equilibrium using the numerically exact inchworm quantum Monte Carlo method. The second chapter studies the possibility of dynamic switching between the metal and Mott-insulator states of the Hubbard model using the one-crossing approximation (OCA). Finally chapter 9 gives some overall conclusions and suggests future directions for the development of hybridization expansion methods.
Open quantum many body systems describe several quantum simulators, ranging from arrays of superconducting circuits hosting correlated states of light, to ultracold atoms in optical lattices in presence of controlled dissipative processes. Their theoretical understanding is hampered by the exponential scaling of their Hilbert space and by their intrinsic nonequilibrium nature, limiting the applicability of many traditional approaches. In this work we extend Dynamical Mean Field Theory (DMFT), a powerful nonperturbative technique developed for strongly correlated quantum many-body systems, to Markovian open quantum systems. Within Open-DMFT, a master equation describing a lattice of bosonic particles is mapped, at large but finite connectivity, onto an impurity problem describing a single site coupled to Markovian and non-Markovian environments, where the latter accounts for quantum fluctuations beyond Gutzwiller mean-field theory. We introduce a new non-perturbative approach to solve the impurity problem, which dresses the impurity, featuring Markovian dissipative channels, with the non-Markovian bath, in a self-consistent scheme based on a resummation of non-crossing diagrams. Finally, we apply our Open-DMFT to solve a driven-dissipative Bose-Hubbard model, which is relevant to current experiments with dissipative ultracold bosons. We show that this model undergoes a dissipative phase transition towards a superfluid phase, which can be equally interpreted as a quantum many-body synchronization transition of an array of quantum van der Pol oscillators. We show how Open-DMFT captures crucial effects due to finite lattice connectivity, such as hopping-induced dissipation and heating, which are completely missed by mean field theories and which lead to a drastic reduction of the broken symmetry phase that we interpret as a desynchronization driven by quantum fluctuations.
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.