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:1206.3201 · Physical Review B

Numerical renormalization group calculation of impurity internal energy and specific heat of quantum impurity models

L. Merker, Theo A. Costi

Published Aug 28, 2012• 12 citations
arXivHTMLPDF

Abstract

We introduce a method to obtain the specific heat of quantum impurity models via a direct calculation of the impurity internal energy requiring only the evaluation of local quantities within a single numerical renormalization group (NRG) calculation for the total system. For the Anderson impurity model we show that the impurity internal energy can be expressed as a sum of purely local static correlation functions and a term that involves also the impurity Green function. The temperature dependence of the latter can be neglected in many cases, thereby allowing the impurity specific heat ${C}_{\mathrm{imp}}$ to be calculated accurately from local static correlation functions; specifically via ${C}_{\mathrm{imp}}=\frac{\ensuremath{\partial}{E}_{\mathrm{ionic}}}{\ensuremath{\partial}T}+\frac{1}{2}\frac{\ensuremath{\partial}{E}_{\mathrm{hyb}}}{\ensuremath{\partial}T}$, where ${E}_{\mathrm{ionic}}$ and ${E}_{\mathrm{hyb}}$ are the energies of the (embedded) impurity and the hybridization energy, respectively. The term involving the Green function can also be evaluated in cases where its temperature dependence is non-negligible, adding an extra term to ${C}_{\mathrm{imp}}$. For the nondegenerate Anderson impurity model, we show by comparison with exact Bethe ansatz calculations that the results recover accurately both the Kondo induced peak in the specific heat at low temperatures as well as the high-temperature peak due to the resonant level. The approach applies to multiorbital and multichannel Anderson impurity models with arbitrary local Coulomb interactions. An application to the Ohmic two-state system and the anisotropic Kondo model is also given, with comparisons to Bethe ansatz calculations. The approach could also be of interest within other impurity solvers, for example, within quantum Monte Carlo techniques.

DOI
10.1103/physrevb.86.075150

Related papers

Full density-matrix numerical renormalization group calculation of impurity susceptibility and specific heat of the Anderson impurity model

L. Merker, Andreas Weichselbaum, Theo A. Costi

Recent developments in the numerical renormalization group (NRG) allow the construction of the full density matrix (FDM) of quantum impurity models [see A. Weichselbaum and J. von Delft, Phys. Rev. Lett. 99, 076402 (2007)] by using the completeness of the eliminated states introduced by F. B. Anders and A. Schiller [F. B. Anders and A. Schiller, Phys. Rev. Lett. 95, 196801 (2005)]. While these developments prove particularly useful in the calculation of transient response and finite-temperature Green's functions of quantum impurity models, they may also be used to calculate thermodynamic properties. In this paper, we assess the FDM approach to thermodynamic properties by applying it to the Anderson impurity model. We compare the results for the susceptibility and specific heat to both the conventional approach within NRG and to exact Bethe ansatz results. We also point out a subtlety in the calculation of the susceptibility (in a uniform field) within the FDM approach. Finally, we show numerically that for the Anderson model, the susceptibilities in response to a local and a uniform magnetic field coincide in the wide-band limit, in accordance with the Clogston-Anderson compensation theorem.

arXiv:1207.2631•Physical Review B•Aug 30, 2012•30 citations

Self-energy method for time-dependent spectral functions of the Anderson impurity model within the time-dependent numerical renormalization group approach

Hoa T. M. Nghiem, Theo A. Costi

The self-energy method for quantum impurity models expresses the correlation part of the self-energy in terms of the ratio of two Green's functions and allows for a more accurate calculation of equilibrium spectral functions than is possible directly from the one-particle Green's function [Bulla et al., J. Phys.: Condens. Matter 10, 8365 (1998)], for example, within the numerical renormalization group method. In addition, the self-energy itself is a central quantity required in the dynamical mean field theory of strongly correlated lattice models. Here, we show how to generalize the self-energy method to the time-dependent situation for the prototype model of strong correlations, the Anderson impurity model. We use the equation-of-motion method to obtain closed expressions for the local Green's function in terms of a time-dependent correlation self-energy, with the latter being given as a ratio of a one-particle time-dependent Green's function and a higher-order correlation function. We benchmark this self-energy approach to time-dependent spectral functions against the direct approach within the time-dependent numerical renormalization group method. The self-energy approach improves the accuracy of time-dependent spectral function calculations, and the closed-form expressions for the Green's function allow for a clear picture of the time evolution of spectral features at the different characteristic time scales. The self-energy approach is of potential interest also for other quantum impurity solvers for real-time evolution, including time-dependent density matrix renormalization group and continuous-time quantum Monte Carlo techniques.

arXiv:2109.01580•Physical review. B./Physical review. B•Nov 12, 2021•3 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

Padé approximant approach for obtaining finite-temperature spectral functions of quantum impurity models using the numerical renormalization group technique

Žiga Osolin, Rok Žitko

We introduce an alternative approach for obtaining smooth finite-temperature spectral functions of quantum impurity models using the numerical renormalization group (NRG) technique. It is based on calculating first the Green's function on the imaginary-frequency axis, followed by an analytic continuation to the real-frequency axis using Pad\'e approximants. The arbitrariness in choosing a suitable kernel in the conventional broadening approach is thereby removed and, furthermore, we find that the Pad\'e method is able to resolve fine details in spectral functions with less artifacts on the scale of $\ensuremath{\omega}\ensuremath{\sim}T$. We discuss the convergence properties with respect to the NRG calculation parameters (discretization $\ensuremath{\Lambda}$, $z$-averaging, truncation cutoff) and the number of Matsubara points taken into account in the analytic continuation. We test the technique on the single-impurity Anderson model and the Hubbard model (within the dynamical mean-field theory). For the Anderson impurity model, we discuss the shape of the Kondo resonance and its temperature dependence. For the Hubbard model, we discuss the inner structure of the Hubbard bands in metallic and insulating solutions at half-filling, as well as in the doped Mott insulator. Based on these test cases, we conclude that the Pad\'e approximant approach provides improved results for spectral functions at low-frequency scales of $\ensuremath{\omega}\ensuremath{\lesssim}T$ and that it is capable of resolving sharp spectral features also at high frequencies.

arXiv:1302.3334•Physical Review B•Jun 28, 2013•35 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

Spectral function of the Anderson impurity model at finite temperatures

Aldo Isidori, David Roosen, Lorenz Bartosch +2 more

Using the functional renormalization group (FRG) and the numerical renormalization group (NRG), we calculate the spectral function of the Anderson impurity model at zero and finite temperatures. In our FRG scheme spin fluctuations are treated nonperturbatively via a suitable Hubbard-Stratonovich field, but vertex corrections are neglected. A comparison with our highly accurate NRG results shows that this FRG scheme provides a good description of the spectral line shape at zero and finite temperatures both in the weak- and strong-coupling regimes, although at zero temperature the FRG is not able to reproduce the known exponential narrowing of the Kondo resonance at strong coupling.

arXiv:1003.3188•Physical Review B•Jun 17, 2010•20 citations
We introduce a method to obtain the specific heat of quantum impurity models via a direct calculation of the impurity internal energy requiring only the evaluation of local quantities within a single numerical renormalization group (NRG) calculation for the total system. For the Anderson impurity model we show that the impurity internal energy can be expressed as a sum of purely local static correlation functions and a term that involves also the impurity Green function. The temperature dependence of the latter can be neglected in many cases, thereby allowing the impurity specific heat Cimp{C}_{\mathrm{imp}}Cimp​ to be calculated accurately from local static correlation functions; specifically via Cimp=\ensuremath∂Eionic\ensuremath∂T+12\ensuremath∂Ehyb\ensuremath∂T{C}_{\mathrm{imp}}=\frac{\ensuremath{\partial}{E}_{\mathrm{ionic}}}{\ensuremath{\partial}T}+\frac{1}{2}\frac{\ensuremath{\partial}{E}_{\mathrm{hyb}}}{\ensuremath{\partial}T}Cimp​=\ensuremath∂T\ensuremath∂Eionic​​+21​\ensuremath∂T\ensuremath∂Ehyb​​, where Eionic{E}_{\mathrm{ionic}}Eionic​ and Ehyb{E}_{\mathrm{hyb}}Ehyb​ are the energies of the (embedded) impurity and the hybridization energy, respectively. The term involving the Green function can also be evaluated in cases where its temperature dependence is non-negligible, adding an extra term to Cimp{C}_{\mathrm{imp}}Cimp​. For the nondegenerate Anderson impurity model, we show by comparison with exact Bethe ansatz calculations that the results recover accurately both the Kondo induced peak in the specific heat at low temperatures as well as the high-temperature peak due to the resonant level. The approach applies to multiorbital and multichannel Anderson impurity models with arbitrary local Coulomb interactions. An application to the Ohmic two-state system and the anisotropic Kondo model is also given, with comparisons to Bethe ansatz calculations. The approach could also be of interest within other impurity solvers, for example, within quantum Monte Carlo techniques.
Recent developments in the numerical renormalization group (NRG) allow the construction of the full density matrix (FDM) of quantum impurity models [see A. Weichselbaum and J. von Delft, Phys. Rev. Lett. 99, 076402 (2007)] by using the completeness of the eliminated states introduced by F. B. Anders and A. Schiller [F. B. Anders and A. Schiller, Phys. Rev. Lett. 95, 196801 (2005)]. While these developments prove particularly useful in the calculation of transient response and finite-temperature Green's functions of quantum impurity models, they may also be used to calculate thermodynamic properties. In this paper, we assess the FDM approach to thermodynamic properties by applying it to the Anderson impurity model. We compare the results for the susceptibility and specific heat to both the conventional approach within NRG and to exact Bethe ansatz results. We also point out a subtlety in the calculation of the susceptibility (in a uniform field) within the FDM approach. Finally, we show numerically that for the Anderson model, the susceptibilities in response to a local and a uniform magnetic field coincide in the wide-band limit, in accordance with the Clogston-Anderson compensation theorem.
The self-energy method for quantum impurity models expresses the correlation part of the self-energy in terms of the ratio of two Green's functions and allows for a more accurate calculation of equilibrium spectral functions than is possible directly from the one-particle Green's function [Bulla et al., J. Phys.: Condens. Matter 10, 8365 (1998)], for example, within the numerical renormalization group method. In addition, the self-energy itself is a central quantity required in the dynamical mean field theory of strongly correlated lattice models. Here, we show how to generalize the self-energy method to the time-dependent situation for the prototype model of strong correlations, the Anderson impurity model. We use the equation-of-motion method to obtain closed expressions for the local Green's function in terms of a time-dependent correlation self-energy, with the latter being given as a ratio of a one-particle time-dependent Green's function and a higher-order correlation function. We benchmark this self-energy approach to time-dependent spectral functions against the direct approach within the time-dependent numerical renormalization group method. The self-energy approach improves the accuracy of time-dependent spectral function calculations, and the closed-form expressions for the Green's function allow for a clear picture of the time evolution of spectral features at the different characteristic time scales. The self-energy approach is of potential interest also for other quantum impurity solvers for real-time evolution, including time-dependent density matrix renormalization group and continuous-time quantum Monte Carlo techniques.
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 introduce an alternative approach for obtaining smooth finite-temperature spectral functions of quantum impurity models using the numerical renormalization group (NRG) technique. It is based on calculating first the Green's function on the imaginary-frequency axis, followed by an analytic continuation to the real-frequency axis using Pad\'e approximants. The arbitrariness in choosing a suitable kernel in the conventional broadening approach is thereby removed and, furthermore, we find that the Pad\'e method is able to resolve fine details in spectral functions with less artifacts on the scale of $\ensuremath{\omega}\ensuremath{\sim}T$. We discuss the convergence properties with respect to the NRG calculation parameters (discretization $\ensuremath{\Lambda}$, zzz-averaging, truncation cutoff) and the number of Matsubara points taken into account in the analytic continuation. We test the technique on the single-impurity Anderson model and the Hubbard model (within the dynamical mean-field theory). For the Anderson impurity model, we discuss the shape of the Kondo resonance and its temperature dependence. For the Hubbard model, we discuss the inner structure of the Hubbard bands in metallic and insulating solutions at half-filling, as well as in the doped Mott insulator. Based on these test cases, we conclude that the Pad\'e approximant approach provides improved results for spectral functions at low-frequency scales of $\ensuremath{\omega}\ensuremath{\lesssim}T$ and that it is capable of resolving sharp spectral features also at high frequencies.
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.
Using the functional renormalization group (FRG) and the numerical renormalization group (NRG), we calculate the spectral function of the Anderson impurity model at zero and finite temperatures. In our FRG scheme spin fluctuations are treated nonperturbatively via a suitable Hubbard-Stratonovich field, but vertex corrections are neglected. A comparison with our highly accurate NRG results shows that this FRG scheme provides a good description of the spectral line shape at zero and finite temperatures both in the weak- and strong-coupling regimes, although at zero temperature the FRG is not able to reproduce the known exponential narrowing of the Kondo resonance at strong coupling.