We propose an out-of-equilibrium impurity model for the dynamical mean-field description of the Hubbard model driven by a finite electric field. The out-of-equilibrium impurity environment is represented by a collection of equilibrium reservoirs at different chemical potentials. We discuss the validity of the impurity model and propose a non-perturbative method, based on a quantum Monte Carlo solver, which provides the steady-state solutions of the impurity and original lattice problems. We discuss the relevance of this approach to other non-equilibrium steady-state contexts.
We propose an out-of-equilibrium impurity model for the dynamical mean-field description of the Hubbard model driven by a finite electric field. The out-of-equilibrium impurity environment is represented by a collection of equilibrium reservoirs at different chemical potentials. We discuss the validity of the impurity model and propose a non-perturbative method, based on a quantum Monte Carlo solver, which provides the steady-state solutions of the impurity and original lattice problems. We discuss the relevance of this approach to other non-equilibrium steady-state contexts.
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.
The Dynamical Mean Field Theory (DMFT) is a powerful tool for calculating highly correlated systems (both bosonic and fermionic) in a state of thermodynamic equilibrium. However, in the case of non-equilibrium states, the method has significant limitations that do not allow obtaining correct results. The stumbling block here is the impurity solver: a method for calculating the dynamics of an open system. In this work we present the prototype of a universal impurity solver for the DMFT method which can solve both equilibrium and non-equilibrium problems. We analyse spectral functions of Bose-Hubbard model on the Bethe lattice in different regimes and show that our solver gives correct and physical results.
This chapter, like the previous one, describes quantum Monte Carlo methods for the simulation of a system of interacting electrons at nonzero temperature. Here, we target a particular class of problems, characterized by a small number of correlated sites or orbitals, that we call impurity problems. While we could apply the methods of the previous chapter to such impurity models, the algorithms discussed here are more advantageous not only because they work in continuous imaginary time and hence lack the Trotter approximation error inherent to the previous methods, but also because they manipulate determinants of matrices of smaller size, which can be handled more efficiently. At the same time, the continuous-time approach is applicable to broader classes of impurity problems. Driving the development of these continuous-time impurity solvers was the desire to perform more efficient simulations of correlated lattice models within the framework of dynamical mean-field theory. This formalism maps the lattice problem onto an impurity problem, whose parameters are determined by a self-consistency condition. In some cases, the self-consistent solution gives a good approximation of the properties of the original lattice problem. When we combine this dynamical mean-field approximation with the local-density approximation for electronic structure calculations, we obtain a powerful scheme for electronic structure calculations of strongly correlated materials. In this chapter we discuss the most important classes of impurity models, sketch the basics of the dynamical mean-field approximation, and detail the different variants of continuous-time impurity solvers. Quantum impurity models A quantum impurity model describes an atom or molecule embedded in some host with which it exchanges electrons, spin, and energy. This exchange allows the impurity to make transitions between different quantum states, and in the presence of interactions within or between the impurity orbitals, these transitions lead to nontrivial dynamical properties. Quantum impurity models play a prominent role, for example, in the theoretical description of the magnetic and electric properties of dilute metal alloys and in theoretical studies of quantum dots and molecular conductors. These models also appear as an auxiliary problem whose solution yields the dynamical mean-field description of correlated lattice models.
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.
In this work, we develop a non-equilibrium steady-state non-crossing approximation (NESS-NCA) impurity solver applicable to general impurity problems. The choice of the NCA as the impurity solver enables both a more accurate description of correlation effects with larger Coulomb interaction and scalability to multi-orbital systems. Based on this development, we investigate strongly correlated non-equilibrium states of a dissipative lattice system under constant electric fields. Both the electronic Coulomb interaction and the electric field are treated non-perturbatively using dynamical mean-field theory in its non-equilibrium steady-state form (NESS-DMFT) with the NESS-NCA impurity solver. We validate our implementation using a half-filled single-band Hubbard model attached to a fictitious free Fermion reservoir, which prevents temperature divergence. As a result, we identify metallic and insulating phases as functions of the electric field and the Coulomb interaction along with a phase coexistence region amid the metal-to-insulator transition (MIT). We find that the MIT driven by the electric field is qualitatively similar to the equilibrium MIT as a function of temperature, differing from results in previous studies using the iterative perturbation theory (IPT) impurity solver. Finally, we highlight the importance of the morphology of a correlated system under the influence of an electric field.
We present a numerically exact steady-state inchworm Monte Carlo method for nonequilibrium quantum impurity models. Rather than propagating an initial state to long times, the method is directly formulated in the steady state. This eliminates any need to traverse the transient dynamics and grants access to a much larger range of parameter regimes at vastly reduced computational costs. We benchmark the method on equilibrium Green's functions of quantum dots in the noninteracting limit and in the unitary limit of the Kondo regime. We then consider correlated materials described with dynamical mean field theory and driven away from equilibrium by a bias voltage. We show that the response of a correlated material to a bias voltage differs qualitatively from the splitting of the Kondo resonance observed in bias-driven quantum dots.