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.
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.
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 ground-state energy, wave-function, hybridization gap and magnetic phasediagram of the Single Site Impurity Anderson Model (SIAM) and the Periodic AndersonModel (PAM) were calculated using Exact-Diagonalization (ED) technique. The EDtechnique was applied to the SIAM and PAM in one dimension before speculating on Ce 4fbeing a test bed of the results, since it has 4f localized electrons in the Heavy fermioncompound CeCu2Si2. Each of the results obtained describe the evolution of the electronicsystem in terms of a simple set of parameters.The parameter regimes studied suggest asmooth phase transition from an antiferromagnetic (AFM) phase to a ferromagnetic (FM)phase at quarter fillings in one dimension when t increases and a magnetic instability wasobserved as U → ∞. The role of hybridization between the conduction electrons and thelocalized f orbital in the Anderson impurity in a semiconductor was also investigated as Vincreases in the SIAM and the PAM.
We develop a method for calculating the self-energy of a quantum impurity coupled to a continuous bath by stochastically generating a distribution of finite Anderson models that are solved by exact diagonalization, using the noninteracting local spectral function as a probability distribution for the sampling. The method enables calculation of the full analytic self-energy and single-particle Green's function in the complex frequency plane, without analytic continuation, and can be used for both finite and zero temperature at arbitrary fillings. Results are in good agreement with imaginary frequency data from continuous-time quantum Monte Carlo calculations for the single-impurity Anderson model and the two-orbital Hubbard model within dynamical mean-field theory (DMFT) as well as real frequency data for self-energy of the single-band Hubbard model within DMFT using the numerical renormalization group. The method should be applicable to a wide range of quantum impurity models and particularly useful when high-precision real frequency results are sought.
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.
Various methods and approximation schemes are used to study many-electron interacting systems. Two important many-particle models, the Anderson model and the Hubbard model, and their electromagnetic properties have been investigated in many parameter regimes, and applied to physical systems. An Anderson single-impurity model Hamiltonian based calculation of the magnetic susceptibility is performed for YbN in the presence of crystal fields using an alteration of the Non-Crossing Approximation proposed by Zwicknagl et.al., incorporating parameters obtained from ab initio band structure calculations. It yields good agreement with experimental data. For the Anderson lattice model, a variational scheme which uses specific many-electron wavefunctions as basis is applied to both one- and two-dimensional systems represented by symmetric Anderson lattice Hamiltonians. Without much computational effort, the ground state energy is well approximated, especially in strong-coupling limit. Some electronic properties are examined using the variational ground state wavefunction. The one-dimensional Hubbard model has been solved exactly for small-size clusters by diagonalizing the Hamiltonian in the basis of many-electron Bloch states. The results for the energy spectrum and eigenfunctions of the ground state and low-lying excited states are presented. Also, mean field calculations of the two-dimensional single-band Hubbard model and Cu-O lattice model (three-band Hubbard model) are carried out for various physical quantities including the energy, occupation probability, staggered magnetization, momentum distribution Fermi surface and density of states, by using a projection operator formalism. To develop a systematic approach to solving many-electron problems, the many-particle partition function for the free electron gas system is explored using a cumulant expansion scheme. Starting from the ground state, the partition function can be approximated to any order in terms of excitation energy. Its application to interacting systems such as the Anderson model and the Hubbard model is briefly discussed.
We investigate an extended version of the periodic Anderson model (the so-called periodic Anderson-Hubbard model) with the aim to understand the role of interaction between conduction electrons in the formation of the heavy-fermion and mixed-valence states. Two methods are used: (i) variational calculation with the Gutzwiller wave function optimizing numerically the ground-state energy and (ii) exact diagonalization of the Hamiltonian for short chains. The f-level occupancy and the renormalization factor of the quasiparticles are calculated as a function of the energy of the f orbital for a wide range of the interaction parameters. The results obtained by the two methods are in reasonably good agreement for the periodic Anderson model. The agreement is maintained even when the interaction between band electrons, Ud, is taken into account, except for the half-filled case. This discrepancy can be explained by the difference between the physics of the one- and higher-dimensional models. We find that this interaction shifts and widens the energy range of the bare f level, where heavy-fermion behavior can be observed. For large-enough Ud this range may lie even above the bare conduction band. The Gutzwiller method indicates a robust transition from Kondo insulator to Mott insulator in the half-filled model, while Ud enhances the quasiparticle mass when the filling is close to half filling.