Materials classified as strongly-correlated systems often exhibit complex and fascinating properties due to interactions among electrons as well as between electrons and other constituents of the material. A common model to describe electronic system with strong on-site Coulomb interaction is Hubbard model. One very powerful approximation method for solving Hubbard model having been widely used over the last few decades is dynamical mean-field theory (DMFT). This method maps the original lattice problem into an impurity problem embedded in a self-consistent bath. Apart from the many variants of implementation of DMFT, it relies on using an impurity solver as part of its algorithm. In this work, rather than solving a Hubbard model, we propose to explore the impurity solver itself for solving a problem of metallic host doped with correlated elements, commonly referred to as Anderson impurity model (AIM). We solve the model using the distributional exact diagonalisation method. Our particular aim is to show how the metal-insulator transition (MIT) occurs in the system and how the phenomenon reflects in its optical conductivity for various model parameters.
Materials classified as strongly-correlated systems often exhibit complex and fascinating properties due to interactions among electrons as well as between electrons and other constituents of the material. A common model to describe electronic system with strong on-site Coulomb interaction is Hubbard model. One very powerful approximation method for solving Hubbard model having been widely used over the last few decades is dynamical mean-field theory (DMFT). This method maps the original lattice problem into an impurity problem embedded in a self-consistent bath. Apart from the many variants of implementation of DMFT, it relies on using an impurity solver as part of its algorithm. In this work, rather than solving a Hubbard model, we propose to explore the impurity solver itself for solving a problem of metallic host doped with correlated elements, commonly referred to as Anderson impurity model (AIM). We solve the model using the distributional exact diagonalisation method. Our particular aim is to show how the metal-insulator transition (MIT) occurs in the system and how the phenomenon reflects in its optical conductivity for various model parameters.
The DOS, the dynamic (optical) conductivity, and the phase diagram of a strongly correlated and strongly disordered paramagnetic Anderson-Hubbard model are analyzed within the generalized dynamical mean field theory (DMFT + Σ approximation). Strong correlations are taken into account by the DMFT, and disorder is taken into account via an appropriate generalization of the self-consistent theory of localization. The DMFT effective single-impurity problem is solved by a numerical renormalization group (NRG); we consider the three-dimensional system with a semielliptic DOS. The correlated metal, Mott insulator, and correlated Anderson insulator phases are identified via the evolution of the DOS and dynamic conductivity, demonstrating both the Mott-Hubbard and Anderson metal-insulator transition and allowing the construction of the complete zero-temperature phase diagram of the Anderson-Hubbard model. Rather unusual is the possibility of a disorder-induced Mott insulator-to-metal transition.
The density of states, the dynamic (optical) conductivity, and the phase diagram of the paramagnetic two-dimensional Anderson-Hubbard model with strong correlations and disorder are analyzed within the generalized dynamical mean field theory (DMFT + Σ approximation). Strong correlations are accounted by the DMFT, while disorder is taken into account via the appropriate generalization of the self-consistent theory of localization. We consider the two-dimensional system with the rectangular “bare” density of states (DOS). The DMFT effective single-impurity problem is solved by numerical renormalization group (NRG). The “correlated metal,” Mott insulator, and correlated Anderson insulator phases are identified from the evolution of the density of states, optical conductivity, and localization length, demonstrating both Mott-Hubbard and Anderson metal-insulator transitions in two-dimensional systems of finite size, allowing us to construct the complete zero-temperature phase diagram of the paramagnetic Anderson-Hubbard model. The localization length in our approximation is practically independent of the strength of Hubbard correlations. But the divergence of the localization length in a finite-size two-dimensional system at small disorder signifies the existence of an effective Anderson transition.
Single impurity Anderson model describes a system consisting of non-interacting conduction electrons coupled with a localized orbital having strongly interacting electrons at a particular site. This model has been proven successful to explain the phenomenon of metal-insulator transition through Anderson localization. Despite the well-understood behaviors of the model, little has been explored theoretically on how the model properties gradually evolve as functions of hybridization parameter, interaction energy, impurity concentration, and temperature. Here, we propose to do a theoretical study on those aspects of a single impurity Anderson model using the distributional exact diagonalization method. We solve the model Hamiltonian by randomly generating sampling distribution of some conducting electron energy levels with various number of occupying electrons. The resulting eigenvalues and eigenstates are then used to define the local single-particle Green function for each sampled electron energy distribution using Lehmann representation. Later, we extract the corresponding self-energy of each distribution, then average over all the distributions and construct the local Green function of the system to calculate the density of states. We repeat this procedure for various values of those controllable parameters, and discuss our results in connection with the criteria of the occurrence of metal-insulator transition in this system.
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.
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.
The optical conductivity of the weakly doped two-dimensional repulsive Hubbard model on the square lattice with the nearest and next-nearest hoppings is calculated within the generalized dynamical mean-field $(\mathrm{DMFT}+{\ensuremath{\sum}}_{\mathbf{p}})$ approach, which includes correlation length scale $\ensuremath{\xi}$ into the standard DMFT equations via the momentum dependent self-energy ${\ensuremath{\sum}}_{\mathbf{p}}$, with a full account of appropriate vertex corrections. This approach takes into consideration the nonlocal dynamical correlations induced, e.g., by short-ranged collective spin-density-wavelike antiferromagnetic spin fluctuations, which (at high enough temperatures) can be viewed as a quenched Gaussian random field with finite correlation length $\ensuremath{\xi}$. The DMFT effective single-impurity problem is solved by numerical renormalization group. We consider both the case of correlated metal with the bandwidth $W\ensuremath{\lesssim}U$ and that of doped Mott insulator with U⪢W (U---the value of local Hubbard interaction). The optical conductivity calculated within $\mathrm{DMFT}+{\ensuremath{\sum}}_{\mathbf{p}}$ demonstrates typical pseudogap behavior within the quasiparticle band, in qualitative agreement with experiments in copper oxide superconductors. For large values of U, pseudogap anomalies are effectively suppressed.