Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:1608.02434 · Physical review. B./Physical review. B

Kondo physics of the Anderson impurity model by distributional exact diagonalization

S. Motahari, Ryan Requist, David Jacob

Published Dec 14, 2016• 11 citations
arXivHTMLPDF

Abstract

The distributional exact diagonalization (DED) scheme is applied to the description of Kondo physics in the Anderson impurity model. DED maps Anderson's problem of an interacting impurity level coupled to an infinite bath onto an ensemble of finite Anderson models, each of which can be solved by exact diagonalization. An approximation to the self-energy of the original infinite model is then obtained from the ensemble-averaged self-energy. Using Friedel's sum rule, we show that the particle number constraint, a central ingredient of the DED scheme, ultimately imposes Fermi liquid behavior on the ensemble-averaged self-energy, and thus is essential for the description of Kondo physics within DED. Using the numerical renormalization group (NRG) method as a benchmark, we show that DED yields excellent spectra, both inside and outside the Kondo regime for a moderate number of bath sites. Only for very strong correlations ($U/\mathrm{\ensuremath{\Gamma}}\ensuremath{\gg}10$) does the number of bath sites needed to achieve good quantitative agreement become too large to be computationally feasible.

DOI
10.1103/physrevb.94.235133

Related papers

Distributional exact diagonalization formalism for quantum impurity models

Mats Granath, Hugo U. R. Strand

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.

arXiv:1201.6160•Physical Review B•Sep 10, 2012•36 citations

Single impurity Anderson model and dynamical mean field theory : a functional renormalization group study

Michael Kinza, Carsten Honerkamp A

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.

•Jan 1, 2014•0 citations

Superperturbation solver for quantum impurity models

Hartmut Hafermann, Jung, C., S. Brener +3 more

We present a very efficient solver for the general Anderson impurity problem. It is based on the perturbation around a solution obtained from exact diagonalization, using a small number of bath sites. We formulate a perturbation theory which is valid for both weak and strong coupling and interpolates between these limits. Good agreement with numerically exact quantum Monte Carlo results is found for a single bath site over a wide range of parameters. In particular, the Kondo resonance in the intermediate-coupling regime is well reproduced for a single bath site and the lowest-order correction. The method is particularly suited for low temperatures, alleviates analytical continuation of imaginary time data due to the absence of statistical noise and can be generalized to obtain dynamical quantities directly on the real axis.

•Jan 1, 2009•36 citations

Superperturbation solver for quantum impurity models

Hartmut Hafermann, C. Jung, S. Brener +3 more

We present a very efficient solver for the general Anderson impurity problem. It is based on the perturbation around a solution obtained from exact diagonalization, using a small number of bath sites. We formulate a perturbation theory which is valid for both weak and strong coupling and interpolates between these limits. Good agreement with numerically exact quantum Monte Carlo results is found for a single bath site over a wide range of parameters. In particular, the Kondo resonance in the intermediate-coupling regime is well reproduced for a single bath site and the lowest-order correction. The method is particularly suited for low temperatures, alleviates analytical continuation of imaginary time data due to the absence of statistical noise and can be generalized to obtain dynamical quantities directly on the real axis.

arXiv:0809.1051•Europhysics Letters (EPL)•Jan 1, 2009•51 citations

Fermi-liquid theory for the single-impurity Anderson model

Christophe Mora, Cătălin Paşcu Moca, Jan von Delft +1 more

We generalize Nozi\`eres' Fermi-liquid theory for the low-energy behavior of the Kondo model to that of the single-impurity Anderson model. In addition to the electrons' phase shift at the Fermi energy, the low-energy Fermi-liquid theory is characterized by four Fermi-liquid parameters: the two given by Nozi\`eres that enter to first order in the excitation energy, and two additional ones that enter to second order and are needed away from particle-hole symmetry. We express all four parameters in terms of zero-temperature physical observables, namely the local charge and spin susceptibilities and their derivatives with respect to the local level position. We determine these in terms of the bare parameters of the Anderson model using Bethe ansatz and numerical renormalization group (NRG) calculations. Our low-energy Fermi-liquid theory applies throughout the crossover from the strong-coupling Kondo regime via the mixed-valence regime to the empty-orbital regime. From the Fermi-liquid theory, we determine the conductance through a quantum dot symmetrically coupled to two leads in the regime of small magnetic field, low temperature, and small bias voltage, and compute the coefficients of the $\ensuremath{\sim}{B}^{2}, \ensuremath{\sim}{T}^{2}$, and $\ensuremath{\sim}{V}^{2}$ terms exactly in terms of the Fermi-liquid parameters. The coefficients of ${T}^{2}, {V}^{2}$, and ${B}^{2}$ are found to change sign during the Kondo to empty-orbital crossover. The crossover becomes universal in the limit that the local interaction is much larger than the level width. For completeness, we also compute the shot noise and discuss the resulting Fano factor.

arXiv:1409.3451•Physical Review B•Aug 10, 2015•54 citations

Single Particle Dynamics of Anderson-like Impurity Models: A Functional Renormalization Group Study

R Hedden

The Anderson Impurity Model is a generic model to describe systems in which one or more energy levels are coupled to an electronic host. On the energy levels the electrons are subject to a Coulomb interaction. Possible applications of this model are magnetic atoms in a metal or quantum dot systems. Already for one impurity it displays very interesting physics including the well known Kondo effect. For two and more impurities the properties of the system in the local moment regime are determined by an interplay of the Kondo effect and the coupling between the impurity spins. Due to the interplay of multiple energy scales, this model poses a serious task for any computational tool. For a single impurity the model can be very efficiently treated using Wilson"s numerical renormalization group (NRG). Since the numerical effort of NRG calculations rises exponentially with the number of impurities, generically it cannot be applied to systems with more than two impurities. The functional renormalization group (fRG) has been very successfully applied to quantum impurity problems in the past. In this work it will be used to calculate the single-particle dynamics for Anderson-like quantum impurity systems in equilibrium. A comparison to NRG in the case of a single impurity will serve as a benchmark before systems with two impurities and quantum dot systems with two and three dots will be studied. The interplay between the Kondo effect and magnetic and orbital ordering will be investigated in these systems. Although the truncation in the presented fRG scheme is motivated by second order perturbation theory, it will be shown that fRG works also in those cases where perturbation theory fails. Generally it will turn out that fRG is capable of treating a broad variety of Anderson impurity systems within the same framework, including asymmetric impurity levels and interaction between multiple impurities.

•Jan 1, 2007•1 citation
The distributional exact diagonalization (DED) scheme is applied to the description of Kondo physics in the Anderson impurity model. DED maps Anderson's problem of an interacting impurity level coupled to an infinite bath onto an ensemble of finite Anderson models, each of which can be solved by exact diagonalization. An approximation to the self-energy of the original infinite model is then obtained from the ensemble-averaged self-energy. Using Friedel's sum rule, we show that the particle number constraint, a central ingredient of the DED scheme, ultimately imposes Fermi liquid behavior on the ensemble-averaged self-energy, and thus is essential for the description of Kondo physics within DED. Using the numerical renormalization group (NRG) method as a benchmark, we show that DED yields excellent spectra, both inside and outside the Kondo regime for a moderate number of bath sites. Only for very strong correlations (U/\ensuremathΓ\ensuremath≫10U/\mathrm{\ensuremath{\Gamma}}\ensuremath{\gg}10U/\ensuremathΓ\ensuremath≫10) does the number of bath sites needed to achieve good quantitative agreement become too large to be computationally feasible.
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.
We present a very efficient solver for the general Anderson impurity problem. It is based on the perturbation around a solution obtained from exact diagonalization, using a small number of bath sites. We formulate a perturbation theory which is valid for both weak and strong coupling and interpolates between these limits. Good agreement with numerically exact quantum Monte Carlo results is found for a single bath site over a wide range of parameters. In particular, the Kondo resonance in the intermediate-coupling regime is well reproduced for a single bath site and the lowest-order correction. The method is particularly suited for low temperatures, alleviates analytical continuation of imaginary time data due to the absence of statistical noise and can be generalized to obtain dynamical quantities directly on the real axis.
We present a very efficient solver for the general Anderson impurity problem. It is based on the perturbation around a solution obtained from exact diagonalization, using a small number of bath sites. We formulate a perturbation theory which is valid for both weak and strong coupling and interpolates between these limits. Good agreement with numerically exact quantum Monte Carlo results is found for a single bath site over a wide range of parameters. In particular, the Kondo resonance in the intermediate-coupling regime is well reproduced for a single bath site and the lowest-order correction. The method is particularly suited for low temperatures, alleviates analytical continuation of imaginary time data due to the absence of statistical noise and can be generalized to obtain dynamical quantities directly on the real axis.
We generalize Nozi\`eres' Fermi-liquid theory for the low-energy behavior of the Kondo model to that of the single-impurity Anderson model. In addition to the electrons' phase shift at the Fermi energy, the low-energy Fermi-liquid theory is characterized by four Fermi-liquid parameters: the two given by Nozi\`eres that enter to first order in the excitation energy, and two additional ones that enter to second order and are needed away from particle-hole symmetry. We express all four parameters in terms of zero-temperature physical observables, namely the local charge and spin susceptibilities and their derivatives with respect to the local level position. We determine these in terms of the bare parameters of the Anderson model using Bethe ansatz and numerical renormalization group (NRG) calculations. Our low-energy Fermi-liquid theory applies throughout the crossover from the strong-coupling Kondo regime via the mixed-valence regime to the empty-orbital regime. From the Fermi-liquid theory, we determine the conductance through a quantum dot symmetrically coupled to two leads in the regime of small magnetic field, low temperature, and small bias voltage, and compute the coefficients of the $\ensuremath{\sim}{B}^{2}, \ensuremath{\sim}{T}^{2}$, and $\ensuremath{\sim}{V}^{2}$ terms exactly in terms of the Fermi-liquid parameters. The coefficients of T2,V2{T}^{2}, {V}^{2}T2,V2, and B2{B}^{2}B2 are found to change sign during the Kondo to empty-orbital crossover. The crossover becomes universal in the limit that the local interaction is much larger than the level width. For completeness, we also compute the shot noise and discuss the resulting Fano factor.
The Anderson Impurity Model is a generic model to describe systems in which one or more energy levels are coupled to an electronic host. On the energy levels the electrons are subject to a Coulomb interaction. Possible applications of this model are magnetic atoms in a metal or quantum dot systems. Already for one impurity it displays very interesting physics including the well known Kondo effect. For two and more impurities the properties of the system in the local moment regime are determined by an interplay of the Kondo effect and the coupling between the impurity spins. Due to the interplay of multiple energy scales, this model poses a serious task for any computational tool. For a single impurity the model can be very efficiently treated using Wilson"s numerical renormalization group (NRG). Since the numerical effort of NRG calculations rises exponentially with the number of impurities, generically it cannot be applied to systems with more than two impurities. The functional renormalization group (fRG) has been very successfully applied to quantum impurity problems in the past. In this work it will be used to calculate the single-particle dynamics for Anderson-like quantum impurity systems in equilibrium. A comparison to NRG in the case of a single impurity will serve as a benchmark before systems with two impurities and quantum dot systems with two and three dots will be studied. The interplay between the Kondo effect and magnetic and orbital ordering will be investigated in these systems. Although the truncation in the presented fRG scheme is motivated by second order perturbation theory, it will be shown that fRG works also in those cases where perturbation theory fails. Generally it will turn out that fRG is capable of treating a broad variety of Anderson impurity systems within the same framework, including asymmetric impurity levels and interaction between multiple impurities.