Quantum transport through a quantum dot: combining the scattering-states numerical renormalization group with nonequilibrium Green functions

Scattering states fulfill the correct boundary conditions of a current carrying open quantum system. Discretizing the energy continuum of these states allows for employing Wilson's numerical renormalization group approach without violating the boundary conditions by using a finite size system. We evolve the analytically known steady-state density operator for a non-interacting quantum-system at finite bias to the full interacting problem by the time-dependent numerical renormalization group after switching on the local charging energy. Using a newly developed algorithm for steady-state nonequilibrium Green functions, we can calculate the current I as function of bias voltage V for arbitrary temperature and magnetic field. A comparison with second-order and GW Kadanoff-Baym-Keldysh results shows excellent agreement for weak interaction strength U.