Open Wilson chain numerical renormalization group approach to steady-state non-equilibrium quantum transport
Abstract
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−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.
Related
- Wilson coefficients from a non-renormalization theorem in 2D SYM
- Numerical Hints for Dyon Condensation at θ=2π via Wilson-'t Hooft Loops in SU(2) Yang-Mills Theory
- Numerical hints for dyon condensation at θ = 2π via Wilson-’t Hooft loops in SU(2) Yang-Mills theory
- Analysis, thermodynamics, and a numerical solver for a pressure-temperature equilibrium closure of the four-equation model
- Optimal Intermediate Hamiltonians for Non-Equilibrium Free Energy Calculations: A Numerical Study of Markov Models
BibTeX
@misc{manaparambil2026open,
title = {{Open Wilson chain numerical renormalization group approach to steady-state non-equilibrium quantum transport}},
author = {Anand Manaparambil and Frithjof B. Anders},
howpublished = {arXiv},
year = {2026},
month = oct,
doi = {10.48550/arxiv.2606.14635},
eprint = {2606.14635},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2606.14635},
note = {CC BY 4.0},
}
RIS
TY - UNPB
TI - Open Wilson chain numerical renormalization group approach to steady-state non-equilibrium quantum transport
AU - Anand Manaparambil
AU - Frithjof B. Anders
PY - 2026
DA - 2026/10/26
PB - arXiv
DO - 10.48550/arxiv.2606.14635
UR - https://arxiv.org/abs/2606.14635
AB - 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 ≈ 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.
ER -