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:0803.3004 · Journal of Physics Condensed Matter

A numerical renormalization group approach to non-equilibrium Green functions for quantum impurity models

Frithjof B. Anders

Published Apr 11, 2008• 42 citations
arXivHTMLPDF

Abstract

We present a method for the calculation of dynamical correlation functions of quantum impurity systems out of equilibrium using Wilson’s numerical renormalization group (NRG). Our formulation is based on a complete basis set of the Wilson chain and embeds the recently derived algorithm for equilibrium spectral functions. Our method fulfils the spectral weight conserving sum-rule exactly by construction. A local Coulomb repulsion U >0 is switched on at t = 0, and the asymptotic steady-state spectral functions are obtained for various values of U as well as magnetic field strength H and temperature T . These benchmark tests show excellent agreement between the time-evolved and the directly calculated equilibrium NRG spectra for finite U . This method could be used for calculating steady-state non-equilibrium spectral functions at finite bias through interacting nanodevices.

DOI
10.1088/0953-8984/20/19/195216

Related papers

Numerical renormalization group approach to Green’s functions for quantum impurity models

Robert Peters, Thomas Pruschke, Frithjof B. Anders

We present a technique for the calculation of dynamical correlation functions of quantum impurity systems in equilibrium with Wilson's numerical renormalization group. Our formulation is based on a complete basis set of the Wilson chain. In contrast to all previous methods, it does not suffer from overcounting of excitations. By construction, it always fulfills sum rules for spectral functions. Furthermore, it accurately reproduces local thermodynamic expectation values, such as occupancy and magnetization, obtained directly from the numerical renormalization group calculations.

arXiv:cond-mat/0607494•Physical Review B•Dec 18, 2006•254 citations

Open Wilson chain numerical renormalization group approach to steady-state nonequilibrium quantum transport

Anand Manaparambil, Frithjof B. Anders

The numerical renormalization group (NRG) approach was developed to identify and quantify different equilibrium regimes of quantum impurity systems (QISs) with unprecedented accuracy by a tailored finite size representation. Out of equilibrium, the steady-state density operator is not of the Boltzmannian form but one that is determined by the imposed boundary conditions. We extend the NRG to the nonequilibrium setting by augmenting each Wilson site with a reservoir, whose coupling functions are calculated via a continuous fraction expansion in order to recover the continuum limit exactly. The nonequilibrium parameters such as a finite bias as well as a finite temperature gradient enters through the Bloch-Redfield tensor (BRT), whose zero eigenvector gives the steady-state density operator. We used the resulting open chain full density matrix (OC-FDM) approach with an effective single lead description to investigate the charge and spin transport through a quantum dot (QD) under finite bias and temperature gradient. The influence of lead asymmetry and an external magnetic field on the transport properties are also studied. For completeness, we have also investigated local properties of the QD, such as charge fluctuations and find excellent agreement with real-time quantum Monte Carlo (RT-QMC) data. The OC-FDM approach was able to explore Kondo energy scales as low as $T_K/D \approx 10^{-8}$ in the non-equilibrium regime, as well as show convergence with the established equilibrium benchmarks, such as a quantitative agreement with full density matrix numerical renormalization group (FDM-NRG) and Fermi-liquid scaling at small bias and temperature. Owing to the effective single lead description, a single OC-FDM data point takes orders of magnitude less time on a standard laptop, compared to other state-of-the-art numerical methods.

arXiv:2606.14635•Physical review. B./Physical review. B•Oct 23, 2026•0 citations

Variational matrix-product-state approach to quantum impurity models

Andreas Weichselbaum, Frank Verstraete, Ulrich Schollwöck +2 more

We present a unified framework for renormalization group methods, including Wilson's numerical renormalization group (NRG), and White's density-matrix renormalization group (DMRG), within the language of matrix-product-states. This allows improvements over Wilson's NRG for quantum impurity models, as we illustrate for the one-channel Kondo model. Moreover, we use a variational method for evaluating Green's functions. The proposed method is more flexible in its description of spectral properties at finite frequencies, opening the way to time-dependent, out-of-equilibrium impurity problems. It also substantially improves computational efficiency for one-channel impurity problems, suggesting potentially linear scaling of complexity for $n$-channel problems.

arXiv:cond-mat/0504305•Physical Review B•Oct 13, 2009•143 citations

Open Wilson chain numerical renormalization group approach to Green's functions

Jan Böker, Frithjof B. Anders

By combining Wilson's numerical renormalization group (NRG) with a modified Bloch-Redfield approach (BRA), we are able to eliminate the artificial broadening of the Lehmann representation of quantum impurity spectral functions required by the standard NRG algorithm. Our approach is based on the exact reproduction of the continuous coupling function in the original quantum impurity model. It augments each chain site of the Wilson chain by a coupling to an additional reservoir. This open Wilson chain is constructed by a continuous fraction expansion, and the coupling function is treated in second order in the context of the BRA. The eigenvalues of the resulting Bloch-Redfield tensor generate a finite lifetime of the NRG excitations that leads to a natural broadening of the spectral functions. We combine this approach with z-averaging and an analytical exact expression for the correlation part of the self-energy to obtain an accurate representation of the spectral function of the original continuum model in the absence and presence of an external magnetic field.

arXiv:2204.03453•Physical review. B./Physical review. B•Jun 21, 2022•1 citation

Dynamical spectral function from numerical renormalization group: A full excitation approach

Ke Yang, Ning-Hua Tong

For a given quantum impurity model, Wilson's numerical renormalization group (NRG) naturally defines a NRG Hamiltonian whose exact eigenstates and eigenenergies are obtainable. We give exact expressions for the free energy, static, as well as dynamical quantities of the NRG Hamiltonian. The dynamical spectral function from this approach contains full excitations including intra- and intershell excitations. For the spin-boson model, we compare the spectral function obtained from the present method and the full density matrix (FDM) method, showing that while both guarantee rigorous sum rule, the full excitation approach avoids the causality problem of FDM method.

arXiv:2006.16488•Physical review. B./Physical review. B•Aug 13, 2020•3 citations

Manual for the Flexible DM-NRG code

Örs Legeza, Cătălin Paşcu Moca, Anna I. Tóth +2 more

Quantum impurity models describe interactions between some local degrees of freedom and a continuum of non-interacting fermionic or bosonic states. The investigation of quantum impurity models is a starting point towards the understanding of more complex strongly correlated systems, but quantum impurity models also provide the description of various correlated mesoscopic structures, biological and chemical processes, atomic physics and describe phenomena such as dissipation or dephasing. Prototypes of these models are the Anderson impurity model, or the single- and multi-channel Kondo models. The numerical renormalization group method (NRG) proposed by Wilson in mid 70's has been used in its original form for a longtime as one of the most accurate and powerful methods to deal with quatum impurity problems. Recently, a number of new developments took place: First, a spectral sum-conserving density matrix NRG approach (DM-NRG) has been developed, which has also been generalized for non-Abelian symmetries. In this manual we introduce some of the basic concepts of the NRG method and present recently developed Flexible DM-NRG code. This code uses user-defined non-Abelian symmetries dynamically, computes spectral functions, expectation values of local operators for user-defined impurity models. The code can also use a uniform density of states as well as a user-defined density of states. The current version of the code assumes fermionic bath's and it uses any number of U(1), SU(2) charge SU(2) or Z(2) symmetries. The Flexible DM-NRG code can be downloaded from http://www.phy.bme.hu/~dmnrg

arXiv:0809.3143•Sep 18, 2008•23 citations
We present a method for the calculation of dynamical correlation functions of quantum impurity systems out of equilibrium using Wilson’s numerical renormalization group (NRG). Our formulation is based on a complete basis set of the Wilson chain and embeds the recently derived algorithm for equilibrium spectral functions. Our method fulfils the spectral weight conserving sum-rule exactly by construction. A local Coulomb repulsion U >0 is switched on at t = 0, and the asymptotic steady-state spectral functions are obtained for various values of U as well as magnetic field strength H and temperature T . These benchmark tests show excellent agreement between the time-evolved and the directly calculated equilibrium NRG spectra for finite U . This method could be used for calculating steady-state non-equilibrium spectral functions at finite bias through interacting nanodevices.
We present a technique for the calculation of dynamical correlation functions of quantum impurity systems in equilibrium with Wilson's numerical renormalization group. Our formulation is based on a complete basis set of the Wilson chain. In contrast to all previous methods, it does not suffer from overcounting of excitations. By construction, it always fulfills sum rules for spectral functions. Furthermore, it accurately reproduces local thermodynamic expectation values, such as occupancy and magnetization, obtained directly from the numerical renormalization group calculations.
The numerical renormalization group (NRG) approach was developed to identify and quantify different equilibrium regimes of quantum impurity systems (QISs) with unprecedented accuracy by a tailored finite size representation. Out of equilibrium, the steady-state density operator is not of the Boltzmannian form but one that is determined by the imposed boundary conditions. We extend the NRG to the nonequilibrium setting by augmenting each Wilson site with a reservoir, whose coupling functions are calculated via a continuous fraction expansion in order to recover the continuum limit exactly. The nonequilibrium parameters such as a finite bias as well as a finite temperature gradient enters through the Bloch-Redfield tensor (BRT), whose zero eigenvector gives the steady-state density operator. We used the resulting open chain full density matrix (OC-FDM) approach with an effective single lead description to investigate the charge and spin transport through a quantum dot (QD) under finite bias and temperature gradient. The influence of lead asymmetry and an external magnetic field on the transport properties are also studied. For completeness, we have also investigated local properties of the QD, such as charge fluctuations and find excellent agreement with real-time quantum Monte Carlo (RT-QMC) data. The OC-FDM approach was able to explore Kondo energy scales as low as TK/D≈10−8T_K/D \approx 10^{-8}TK​/D≈10−8 in the non-equilibrium regime, as well as show convergence with the established equilibrium benchmarks, such as a quantitative agreement with full density matrix numerical renormalization group (FDM-NRG) and Fermi-liquid scaling at small bias and temperature. Owing to the effective single lead description, a single OC-FDM data point takes orders of magnitude less time on a standard laptop, compared to other state-of-the-art numerical methods.
We present a unified framework for renormalization group methods, including Wilson's numerical renormalization group (NRG), and White's density-matrix renormalization group (DMRG), within the language of matrix-product-states. This allows improvements over Wilson's NRG for quantum impurity models, as we illustrate for the one-channel Kondo model. Moreover, we use a variational method for evaluating Green's functions. The proposed method is more flexible in its description of spectral properties at finite frequencies, opening the way to time-dependent, out-of-equilibrium impurity problems. It also substantially improves computational efficiency for one-channel impurity problems, suggesting potentially linear scaling of complexity for nnn-channel problems.
By combining Wilson's numerical renormalization group (NRG) with a modified Bloch-Redfield approach (BRA), we are able to eliminate the artificial broadening of the Lehmann representation of quantum impurity spectral functions required by the standard NRG algorithm. Our approach is based on the exact reproduction of the continuous coupling function in the original quantum impurity model. It augments each chain site of the Wilson chain by a coupling to an additional reservoir. This open Wilson chain is constructed by a continuous fraction expansion, and the coupling function is treated in second order in the context of the BRA. The eigenvalues of the resulting Bloch-Redfield tensor generate a finite lifetime of the NRG excitations that leads to a natural broadening of the spectral functions. We combine this approach with z-averaging and an analytical exact expression for the correlation part of the self-energy to obtain an accurate representation of the spectral function of the original continuum model in the absence and presence of an external magnetic field.
For a given quantum impurity model, Wilson's numerical renormalization group (NRG) naturally defines a NRG Hamiltonian whose exact eigenstates and eigenenergies are obtainable. We give exact expressions for the free energy, static, as well as dynamical quantities of the NRG Hamiltonian. The dynamical spectral function from this approach contains full excitations including intra- and intershell excitations. For the spin-boson model, we compare the spectral function obtained from the present method and the full density matrix (FDM) method, showing that while both guarantee rigorous sum rule, the full excitation approach avoids the causality problem of FDM method.
Quantum impurity models describe interactions between some local degrees of freedom and a continuum of non-interacting fermionic or bosonic states. The investigation of quantum impurity models is a starting point towards the understanding of more complex strongly correlated systems, but quantum impurity models also provide the description of various correlated mesoscopic structures, biological and chemical processes, atomic physics and describe phenomena such as dissipation or dephasing. Prototypes of these models are the Anderson impurity model, or the single- and multi-channel Kondo models. The numerical renormalization group method (NRG) proposed by Wilson in mid 70's has been used in its original form for a longtime as one of the most accurate and powerful methods to deal with quatum impurity problems. Recently, a number of new developments took place: First, a spectral sum-conserving density matrix NRG approach (DM-NRG) has been developed, which has also been generalized for non-Abelian symmetries. In this manual we introduce some of the basic concepts of the NRG method and present recently developed Flexible DM-NRG code. This code uses user-defined non-Abelian symmetries dynamically, computes spectral functions, expectation values of local operators for user-defined impurity models. The code can also use a uniform density of states as well as a user-defined density of states. The current version of the code assumes fermionic bath's and it uses any number of U(1), SU(2) charge SU(2) or Z(2) symmetries. The Flexible DM-NRG code can be downloaded from http://www.phy.bme.hu/~dmnrg