Dynamical mean-field theory has become a method of choice for calculating spectral properties of strongly correlated systems. In this chapter the concept and the formalism of dynamical mean-field theory are described using a one-band Hubbard model as an illustration. The Anderson impurity model, which forms the auxiliary system, is used as a starting point to develop the formal dynamical mean-field equations. Functional integral formulation, well-suited for quantum Monte Carlo simulations, is also presented and analogy with density functional theory is made via the Legendre transform formulation. The success of the theory is illustrated by some generic examples.
Dynamical mean-field theory has become a method of choice for calculating spectral properties of strongly correlated systems. In this chapter the concept and the formalism of dynamical mean-field theory are described using a one-band Hubbard model as an illustration. The Anderson impurity model, which forms the auxiliary system, is used as a starting point to develop the formal dynamical mean-field equations. Functional integral formulation, well-suited for quantum Monte Carlo simulations, is also presented and analogy with density functional theory is made via the Legendre transform formulation. The success of the theory is illustrated by some generic examples.
Think globally, act locally. Summary Dynamical mean-field theory is designed to treat systems with local effective interactions that are strong compared with the independent-particle terms that lead to delocalized band-like states. Interactions are taken into account by a many-body calculation for an auxiliary system, a site embedded in a dynamical mean field, that is chosen to best represent the coupling to the rest of the crystal. The methods are constructed to be exact in three limits: interacting electrons on isolated sites, a lattice with no interactions, and the limit of infinite dimensions d → ∞ where mean-field theory is exact. This chapter is devoted to the general formulation, the single-site approximation where the calculation of the self-energy is mapped onto a self-consistent quantum impurity problem, and instructive examples for the Hubbard model on a d → ∞ Bethe lattice. Further developments and applications are the topics of Chs. 17–21. One of the most rewarding features of condensed matter theory is the ability to address difficult problems from different points of view. The preceding Chs. 9–15 present an approach based on perturbation expansions in the Coulomb interaction. In particular, the GW approximation for the self-energy has proven to be extremely successful in describing electronic spectra of many materials, as described in Ch. 13. The methods can be applied to the ordered states of materials with d and f states, for example, ferromagnetic Ni and anti-ferromagnetic NiO, as described in Secs. 13.4 and 20.7. However, the GW and related approximations have difficulties in treating cases with degenerate or nearly degenerate states and low-energy excitations. Present-day methods do not describe phenomena like the fluctuations of local moments in a ferromagnetic material above the Curie temperature or the insulating character of NiO in the paramagnetic phase; more generally, they have difficulty describing strong correlation. The topic of this chapter and Chs. 17–21 is dynamical mean-field theory, which is also a Green's function method in which the key quantity is the self-energy. However, it is designed to treat strong interactions for electrons in localized atomic-like states, such as the d and f states in transition metals, lanthanide and actinide elements and compounds.
An essential role for the description of correlated quantum many-particle systems is played by impurity problems. These consist of a small number of localized orbitals, where the electrons are subject to a Coulomb-interaction, coupled to noninteracting bath degrees of freedom. Impurity models are used for the description of magnetic impurities in metals, of quantum points in nanostructures, and in the context of the dynamical mean field theory (DMFT). For their solution a large number of numerical methods exists. We mention here quantum Monte Carlo methods, the numerical renormalization group, or the exact diagonalization method. In this thesis we introduce a new approach to impurity problems, which is based on the functional renormalization group (fRG). In contrast to the methods mentioned above this approach is approximative, but it can be applied with comparatively lower numerical effort. The aim of this thesis is to examine the advantages and limits of this method. An important role in this connection is played by the single impurity Anderson model (SIAM), which is already understood very well with other methods and is seen as a prototype for models with strong local correlations. The spin fluctuations of this model are governed by the so-called Kondo energy scale and a central question of this thesis is whether this scale can be reproduced with our approach. We begin with a short introduction to the Green's function formalism for correlated quantum many particle systems. Then the fRG flow equations for the one-particle irreducible vertex functions are deduced, and several approximation schemes are discussed. In chapter 4 we introduce the fRG scheme used in this thesis. At first the SIAM is mapped to a semi-infinite chain, in which the interacting orbital is given by the first site of this chain. The system is then subdivided into two parts. The first part, which is called 'core', consists of the interacting orbital and the first L bath sites. The remaining bath sites form the second part. At the beginning both parts are decoupled from each other such that the core can be solved exactly. Starting with this exact solution the coupling between the core and the remaining bath sites is switched on slowly in a renormalization group flow. This way it is possible to calculate the one-particle and two-particle correlation functions on the correlated site. We call this flow scheme 'hybridization flow'. The flow equations are formulated in an effective theory on the first bath site outside the core, which turns out to be advantageous compared to other implementations of the hybridization flow. In chapter 5 this hybridization flow method is applied to the SIAM with semi-elliptic bath density of states. We discuss the differences which arise for different core-sizes (L = 0,1,2,3) and for different truncations of the fRG flow equations. In the cases L=0 and L=2 the local Fermi liquid properties of the SIAM and the associated Kondo scale can not be reproduced in both approximation schemes. For the core sizes L=1 and L=3 the ground state of the core is already a spin singlet state in the beginning of the flow and the Fermi liquid properties are reproduced. The Kondo energy scale, which determines the width of the central quasi-particle peak and the size of the local spin-susceptibility, is however not accurately resolved. This is partly due to the finite temperature used in our implementation such that for larger interaction strength the temperature is above the Kondo energy scale. The application of the hybridization flow method to quantum impurity problems in the context of the DMFT is discussed in chapter 6. In this approach lattice models with a local interaction, like the Hubbard model, are mapped to an effective impurity problem for a single lattice site embedded in a dynamical mean field representing the influence of the other electrons. This scheme is exact in infinitely many dimensions and a nonperturbative approximation method for finite-dimensional systems. In particular the Mott metal-insulator transition in the Hubbard model is successfully described by this method. One class of quantities, which is easily accessible from the hybridization flow scheme, is given by the local two-particle correlation functions. These play an essential role in non-local extensions of the DMFT. Therefore we mainly focus on these quantities. At first we apply our flow scheme to the Hubbard model on the Bethe lattice in infinite dimensions and calculate the local one-particle irreducible vertex functions in the insulating and the metallic phase in good agreement with previous calculations, which use the exact diagonalization method as impurity solver. After this we extend our flow scheme to a cluster DMFT method, which includes short-range antiferromagnetic correlations and calculate again the local and non-local vertex functions for the two-dimensional Hubbard model. We compare these to the results of the single-site DMFT.
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.
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 develop a projective quantum Monte Carlo algorithm of the Hirsch-Fye type for obtaining ground state properties of the Anderson impurity model. This method is employed to solve the self-consistency equations of dynamical mean field theory. It is shown that the approach converges rapidly to the ground state so that reliable zero-temperature results are obtained. As a first application, we study the Mott-Hubbard metal-insulator transition of the frustrated one-band Hubbard model, reconfirming the numerical renormalization group results.
To explore correlated electrons in the presence of local and nonlocal disorder, the Blackman-Esterling-Berk method for averaging over off-diagonal disorder is implemented into dynamical mean-field theory using tensor notation. The impurity model combining disorder and correlations is solved using the recently developed fork tensor-product state solver, which allows one to calculate the single particle spectral functions on the real-frequency axis. In the absence of off-diagonal hopping, we establish exact bounds of the spectral function of the noninteracting Bethe lattice with coordination number Z. In the presence of interaction, the Mott insulating paramagnetic phase of the one-band Hubbard model is computed at zero temperature in alloys with site- and off-diagonal disorder. When the Hubbard U parameter is increased, transitions from an alloy band insulator through a correlated metal into a Mott insulating phase are found to take place.