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:1806.01580 · Physical review. B./Physical review. B

Ground-state properties of the symmetric single-impurity Anderson model on a ring from density-matrix renormalization group, Hartree-Fock, and Gutzwiller theory

Gergely Barcza, Florian Gebhard, Thorben Linneweber, Örs Legeza

Published Apr 19, 2019• 8 citations
arXivHTMLPDF

Abstract

We analyze the ground-state energy, magnetization, magnetic susceptibility, and Kondo screening cloud of the symmetric single-impurity Anderson model (SIAM) that is characterized by the bandwidth $W$, the impurity interaction strength $U$, and the local hybridization $V$. We compare Gutzwiller variational and magnetic Hartree-Fock results in the thermodynamic limit with numerically exact data from the density-matrix renormalization group (DMRG) method on large rings. To improve the DMRG performance, we use a canonical transformation to map the SIAM onto a chain with half the system size and open boundary conditions. We compare to Bethe ansatz results for the ground-state energy, magnetization, and spin susceptibility that become exact in the wide-band limit. Our detailed comparison shows that the field-theoretical description is applicable to the SIAM on a ring for a broad parameter range. Hartree-Fock theory gives an excellent ground-state energy and local moment for intermediate and strong interactions. However, it lacks spin fluctuations and thus cannot screen the impurity spin. The Gutzwiller variational energy bound becomes very poor for large interactions because it does not describe properly the charge fluctuations. Nevertheless, the Gutzwiller approach provides a qualitatively correct description of the zero-field susceptibility and the Kondo screening cloud. The DMRG provides excellent data for the ground-state energy and the magnetization for finite external fields. At strong interactions, finite-size effects make it extremely difficult to recover the exponentially large zero-field susceptibility and the mesoscopically large Kondo screening cloud.

DOI
10.1103/physrevb.99.165130

Related papers

Symmetric single-impurity Kondo model on a tight-binding chain: Comparison of analytical and numerical ground-state approaches

Gergely Barcza, Kevin Bauerbach, Fabian Eickhoff +3 more

The Kondo model describes a magnetic impurity in a metallic host. Here, the authors use variational approaches (Lanczos, Yosida, Gutzwiller) and quasi-exact numerical techniques (numerical renormalization group, density matrix renormalization group) to study the ground-state energy and magnetization, and compare these properties to Bethe ansatz results. The Gutzwiller wave function is an excellent trial state for strong coupling. For weak coupling, the numerical renormalization group gives the best description of the magnetic properties. A combination of Bethe ansatz and perturbation theory provides new benchmark results for the zero-field susceptibility.

arXiv:1911.08279•Physical review. B./Physical review. B•Feb 24, 2020•9 citations

Density matrix renormalization group study in energy space for a single-impurity Anderson model and an impurity quantum phase transition

Tomonori Shirakawa, Seiji Yunoki

The density matrix renormalization group method is introduced in energy space to study Anderson impurity models. The method allows for calculations in the thermodynamic limit and is advantageous for studying not only the dynamical properties, but also the quantum entanglement of the ground state at the vicinity of an impurity quantum phase transition. This method is applied to obtain numerically exactly the ground-state phase diagram of the single-impurity Anderson model on the honeycomb lattice at half-filling. The calculation of local static quantities shows that the phase diagram contains two distinct phases, the local moment (LM) phase and the asymmetric strong coupling (ASC) phase, but no Kondo screening phase. These results are supported by the local spin and charge excitation spectra, which exhibit qualitatively different behavior in these two phases and also reveal the existence of the valence fluctuating point at the phase boundary. For comparison, we also study the low-energy effective pseudogap Anderson model using the method introduced here. Although the high-energy excitations are obviously different, we find that the ground-state phase diagram and the asymptotically low-energy excitations are in good quantitative agreement with those for the single-impurity Anderson model on the honeycomb lattice, thus providing a quantitative justification for the previous studies based on low-energy approximate approaches. Furthermore, we find that the lowest entanglement level is doubly degenerate for the LM phase, whereas it is singlet for the ASC phase and is accidentally threefold degenerate at the valence fluctuating point. This should be contrasted with the degeneracy of the energy spectrum because the ground state is found to be always singlet. Our results therefore clearly demonstrate that the low-lying entanglement spectrum can be used to determine with high accuracy the phase boundary of the impurity quantum phase transition.

arXiv:1604.00721•Physical review. B./Physical review. B•May 13, 2016•10 citations

Charge susceptibility in Kondo systems at half filling: DMRG study

Peter Schmitteckert, Ferdinand Evers

Abstract In this brief communication, the charge susceptibility, χc,imp, is calculated for the single impurity Anderson model in the Kondo regime with the density matrix renormalization group theory (DMRG). Excellent agreement with the Bethe‐Ansatz (BA) results is found in the appropriate limit where the (effective) level broadening Γ☆ and the charging energy, U, are much smaller than the conduction electron bandwidth, D. Deviations occur for large interactions U>D, where we find that χc,imp decays faster, approximately by a factor D/U, as compared to the BA‐prediction. Our work is of methodological relevance, e.g., for studies of the Kondo‐effect in molecular systems, in the sense that it prepares treatments of (spinful) few‐orbital models in the Kondo regime.

Annalen der Physik•Feb 2, 2012•3 citations

Ground state magnetic susceptibility of Kondo impurities

Bedoor Alkurtass

Impurity models describe localized impurities interacting with the spins of electrons. Strong correlations lead to non-trivial behaviours of the model. In this work we present 3 different impurity models: (1) the single-impurity Kondo model; (2) the two-impurity Kondo model; (3) the two-channel Kondo model. For each model we present a spin chain emulation and utilize the density matrix renormalization group (DMRG). We study the magnetization and zero-temperature susceptibility of the impurity when a local magnetic field is applied. We get the expected relations theoretically.

Kuwait Journal of Science•Jun 5, 2024•1 citation

Block Lanczos density-matrix renormalization group method for general Anderson impurity models: Application to magnetic impurity problems in graphene

Tomonori Shirakawa, Seiji Yunoki

We introduce a block Lanczos (BL) recursive technique to construct quasi-one-dimensional models, suitable for density-matrix renormalization group (DMRG) calculations, from single- as well as multiple-impurity Anderson models in any spatial dimensions. This new scheme, named BL-DMRG method, allows us to calculate not only local but also spatially dependent static and dynamical quantities of the ground state for general Anderson impurity models without losing elaborate geometrical information of the lattice. We show that the BL-DMRG method can be easily extended to treat a multiorbital Anderson impurity model where not only inter- and intraorbital Coulomb interactions but also Hund's coupling and pair hopping interactions are included. We also show that the symmetry adapted BL bases can be utilized, when it is appropriate, to reduce the computational cost. As a demonstration, we apply the BL-DMRG method to three different models for graphene with a structural defect and with a single hydrogen or fluorine absorbed, where a single Anderson impurity is coupled to conduction electrons in the honeycomb lattice. These models include (i) a single adatom on the honeycomb lattice, (ii) a substitutional impurity in the honeycomb lattice, and (iii) an effective model for a single carbon vacancy in graphene. Our analysis of the local dynamical magnetic susceptibility and the local density of states at the impurity site reveals that, for the particle-hole symmetric case at half-filling of electron density, the ground state of model (i) behaves as an isolated magnetic impurity with no Kondo screening, while the ground state of the other two models forms a spin-singlet state where the impurity moment is screened by the conduction electrons. We also calculate the real-space dependence of the spin-spin correlation functions between the impurity site and the conduction sites for these three models. Our results clearly show that, reflecting the presence or absence of unscreened magnetic moment at the impurity site, the spin-spin correlation functions decay as $\ensuremath{\propto}$ ${r}^{\ensuremath{-}3}$, differently from the noninteracting limit ($\ensuremath{\propto}$ ${r}^{\ensuremath{-}2}$), for model (i) and as $\ensuremath{\propto}$ ${r}^{\ensuremath{-}4}$, exactly the same as the noninteracting limit, for models (ii) and (iii) in the asymptotic $r$, where $r$ is the distance between the impurity site and the conduction site. Finally, based on our results, we shed light on recent experiments on graphene where the formation of local magnetic moments as well as the Kondo-like behavior have been observed.

arXiv:1409.7995•Physical Review B•Nov 5, 2014•33 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
We analyze the ground-state energy, magnetization, magnetic susceptibility, and Kondo screening cloud of the symmetric single-impurity Anderson model (SIAM) that is characterized by the bandwidth WWW, the impurity interaction strength UUU, and the local hybridization VVV. We compare Gutzwiller variational and magnetic Hartree-Fock results in the thermodynamic limit with numerically exact data from the density-matrix renormalization group (DMRG) method on large rings. To improve the DMRG performance, we use a canonical transformation to map the SIAM onto a chain with half the system size and open boundary conditions. We compare to Bethe ansatz results for the ground-state energy, magnetization, and spin susceptibility that become exact in the wide-band limit. Our detailed comparison shows that the field-theoretical description is applicable to the SIAM on a ring for a broad parameter range. Hartree-Fock theory gives an excellent ground-state energy and local moment for intermediate and strong interactions. However, it lacks spin fluctuations and thus cannot screen the impurity spin. The Gutzwiller variational energy bound becomes very poor for large interactions because it does not describe properly the charge fluctuations. Nevertheless, the Gutzwiller approach provides a qualitatively correct description of the zero-field susceptibility and the Kondo screening cloud. The DMRG provides excellent data for the ground-state energy and the magnetization for finite external fields. At strong interactions, finite-size effects make it extremely difficult to recover the exponentially large zero-field susceptibility and the mesoscopically large Kondo screening cloud.
The Kondo model describes a magnetic impurity in a metallic host. Here, the authors use variational approaches (Lanczos, Yosida, Gutzwiller) and quasi-exact numerical techniques (numerical renormalization group, density matrix renormalization group) to study the ground-state energy and magnetization, and compare these properties to Bethe ansatz results. The Gutzwiller wave function is an excellent trial state for strong coupling. For weak coupling, the numerical renormalization group gives the best description of the magnetic properties. A combination of Bethe ansatz and perturbation theory provides new benchmark results for the zero-field susceptibility.
The density matrix renormalization group method is introduced in energy space to study Anderson impurity models. The method allows for calculations in the thermodynamic limit and is advantageous for studying not only the dynamical properties, but also the quantum entanglement of the ground state at the vicinity of an impurity quantum phase transition. This method is applied to obtain numerically exactly the ground-state phase diagram of the single-impurity Anderson model on the honeycomb lattice at half-filling. The calculation of local static quantities shows that the phase diagram contains two distinct phases, the local moment (LM) phase and the asymmetric strong coupling (ASC) phase, but no Kondo screening phase. These results are supported by the local spin and charge excitation spectra, which exhibit qualitatively different behavior in these two phases and also reveal the existence of the valence fluctuating point at the phase boundary. For comparison, we also study the low-energy effective pseudogap Anderson model using the method introduced here. Although the high-energy excitations are obviously different, we find that the ground-state phase diagram and the asymptotically low-energy excitations are in good quantitative agreement with those for the single-impurity Anderson model on the honeycomb lattice, thus providing a quantitative justification for the previous studies based on low-energy approximate approaches. Furthermore, we find that the lowest entanglement level is doubly degenerate for the LM phase, whereas it is singlet for the ASC phase and is accidentally threefold degenerate at the valence fluctuating point. This should be contrasted with the degeneracy of the energy spectrum because the ground state is found to be always singlet. Our results therefore clearly demonstrate that the low-lying entanglement spectrum can be used to determine with high accuracy the phase boundary of the impurity quantum phase transition.
Abstract In this brief communication, the charge susceptibility, χc,imp, is calculated for the single impurity Anderson model in the Kondo regime with the density matrix renormalization group theory (DMRG). Excellent agreement with the Bethe‐Ansatz (BA) results is found in the appropriate limit where the (effective) level broadening Γ☆ and the charging energy, U, are much smaller than the conduction electron bandwidth, D. Deviations occur for large interactions U>D, where we find that χc,imp decays faster, approximately by a factor D/U, as compared to the BA‐prediction. Our work is of methodological relevance, e.g., for studies of the Kondo‐effect in molecular systems, in the sense that it prepares treatments of (spinful) few‐orbital models in the Kondo regime.
Impurity models describe localized impurities interacting with the spins of electrons. Strong correlations lead to non-trivial behaviours of the model. In this work we present 3 different impurity models: (1) the single-impurity Kondo model; (2) the two-impurity Kondo model; (3) the two-channel Kondo model. For each model we present a spin chain emulation and utilize the density matrix renormalization group (DMRG). We study the magnetization and zero-temperature susceptibility of the impurity when a local magnetic field is applied. We get the expected relations theoretically.
We introduce a block Lanczos (BL) recursive technique to construct quasi-one-dimensional models, suitable for density-matrix renormalization group (DMRG) calculations, from single- as well as multiple-impurity Anderson models in any spatial dimensions. This new scheme, named BL-DMRG method, allows us to calculate not only local but also spatially dependent static and dynamical quantities of the ground state for general Anderson impurity models without losing elaborate geometrical information of the lattice. We show that the BL-DMRG method can be easily extended to treat a multiorbital Anderson impurity model where not only inter- and intraorbital Coulomb interactions but also Hund's coupling and pair hopping interactions are included. We also show that the symmetry adapted BL bases can be utilized, when it is appropriate, to reduce the computational cost. As a demonstration, we apply the BL-DMRG method to three different models for graphene with a structural defect and with a single hydrogen or fluorine absorbed, where a single Anderson impurity is coupled to conduction electrons in the honeycomb lattice. These models include (i) a single adatom on the honeycomb lattice, (ii) a substitutional impurity in the honeycomb lattice, and (iii) an effective model for a single carbon vacancy in graphene. Our analysis of the local dynamical magnetic susceptibility and the local density of states at the impurity site reveals that, for the particle-hole symmetric case at half-filling of electron density, the ground state of model (i) behaves as an isolated magnetic impurity with no Kondo screening, while the ground state of the other two models forms a spin-singlet state where the impurity moment is screened by the conduction electrons. We also calculate the real-space dependence of the spin-spin correlation functions between the impurity site and the conduction sites for these three models. Our results clearly show that, reflecting the presence or absence of unscreened magnetic moment at the impurity site, the spin-spin correlation functions decay as $\ensuremath{\propto}$ ${r}^{\ensuremath{-}3}$, differently from the noninteracting limit ($\ensuremath{\propto}$ ${r}^{\ensuremath{-}2}$), for model (i) and as $\ensuremath{\propto}$ ${r}^{\ensuremath{-}4}$, exactly the same as the noninteracting limit, for models (ii) and (iii) in the asymptotic rrr, where rrr is the distance between the impurity site and the conduction site. Finally, based on our results, we shed light on recent experiments on graphene where the formation of local magnetic moments as well as the Kondo-like behavior have been observed.
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.