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

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

Jan Böker, Frithjof B. Anders

Published Jun 21, 2022• 1 citation
arXivHTMLPDF

Abstract

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.

DOI
10.1103/physrevb.105.235127

Related papers

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

Frithjof B. Anders

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.

arXiv:0803.3004•Journal of Physics Condensed Matter•Apr 11, 2008•42 citations

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

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

Restoring the continuum limit in the time-dependent numerical renormalization group approach

Jan Böker, Frithjof B. Anders

The continuous coupling function in quantum impurity problems is exactly partitioned into a part represented by a finite-size Wilson chain and a part represented by a set of additional reservoirs, each coupled to one Wilson chain site. These additional reservoirs represent high-energy modes of the environment neglected by the numerical renormalization group and are required to restore the continuum limit of the original problem. We present a hybrid time-dependent numerical renormalization group approach which combines an accurate numerical renormalization group treatment of the nonequilibrium dynamics on the finite-size Wilson chain with a Bloch-Redfield formalism to include the effect of these additional reservoirs. Our approach overcomes the intrinsic shortcoming of the time-dependent numerical renormalization group approach induced by the bath discretization with a Wilson parameter $\mathrm{\ensuremath{\Lambda}}>1$. We analytically prove that for a system with a single chemical potential, the thermal equilibrium reduced density operator is the steady-state solution of the Bloch-Redfield master equation. For the numerical solution of this master equation, a Lanczos method is employed which couples all energy shells of the numerical renormalization group. The presented hybrid approach is applied to the real-time dynamics in correlated fermionic QISs. An analytical solution of the resonant-level model serves as a benchmark for the accuracy of the method which is then applied to nontrivial models, such as the interacting resonant-level model and the single-impurity Anderson model.

arXiv:2006.07936•Physical review. B./Physical review. B•Aug 31, 2020•2 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

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
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.
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.
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.
The continuous coupling function in quantum impurity problems is exactly partitioned into a part represented by a finite-size Wilson chain and a part represented by a set of additional reservoirs, each coupled to one Wilson chain site. These additional reservoirs represent high-energy modes of the environment neglected by the numerical renormalization group and are required to restore the continuum limit of the original problem. We present a hybrid time-dependent numerical renormalization group approach which combines an accurate numerical renormalization group treatment of the nonequilibrium dynamics on the finite-size Wilson chain with a Bloch-Redfield formalism to include the effect of these additional reservoirs. Our approach overcomes the intrinsic shortcoming of the time-dependent numerical renormalization group approach induced by the bath discretization with a Wilson parameter $\mathrm{\ensuremath{\Lambda}}>1$. We analytically prove that for a system with a single chemical potential, the thermal equilibrium reduced density operator is the steady-state solution of the Bloch-Redfield master equation. For the numerical solution of this master equation, a Lanczos method is employed which couples all energy shells of the numerical renormalization group. The presented hybrid approach is applied to the real-time dynamics in correlated fermionic QISs. An analytical solution of the resonant-level model serves as a benchmark for the accuracy of the method which is then applied to nontrivial models, such as the interacting resonant-level model and the single-impurity Anderson model.
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 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.