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

Renormalization group transport theory for open quantum systems: Charge fluctuations in multilevel quantum dots in and out of equilibrium

Carsten J. Lindner, Fabian B. Kugler, V. Meden, Herbert Schoeller

Published May 24, 2019• 21 citations
arXivHTMLPDF

Abstract

We present the real-time renormalization group (RTRG) method as a method to describe the stationary state current through generic multilevel quantum dots in nonequilibrium. The employed approach consists of a very rudimentary approximation for the renormalization group (RG) equations which neglects all vertex corrections while it provides a means to compute the effective dot Liouvillian self-consistently. Being based on a weak-coupling expansion in the tunneling between dot and reservoirs, the RTRG approach turns out to reliably describe charge fluctuations in and out of equilibrium for arbitrary coupling strength, even at zero temperature. We confirm this in the linear response regime with a benchmark against highly accurate numerical renormalization group data in the exemplary case of three-level quantum dots. For small to intermediate bias voltages and weak Coulomb interactions, we find an excellent agreement between RTRG and functional renormalization group data, which can be expected to be accurate in this regime. As a consequence, we advertise the presented RTRG approach as an efficient and versatile tool to describe charge fluctuations in quantum dot systems.

DOI
10.1103/physrevb.99.205142

Related papers

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

Kondo effect out of equilibrium

Thomas Korb, Herbert Schoeller

Correlations are a central ingredient of modern condensed matter physics. With the advent of mesoscopic physics, it became possible to study correlated systems in non-equilibrium situations, which can be reached by e.g. driving a constant current through a quantum dot. A basic question is how correlations on the one hand and a non-equilibrium current on the other hand influence each other. Here this interplay is studied for the model system of a quantum dot in the Kondo regime using the real time renormalization group (RTRG) method. This method is a functional renormalization group (RG) method allowing a non-perturbative description of correlated systems in non-equilibrium. For the Kondo quantum dot the method is generalized to the case of double vertices. The RTRG method not only provides equations for the calculation of the renormalized coupling between the quantum dot spin and the reservoir electrons, but also equations for the calculation of the non-equilibrium decoherence rate and the non-equilibrium current.The RG equations are derived to second order. This approximation is expected to be valid as long as the quantum dot does not enter the strong coupling regime, i.e. the effective coupling strength between the quantum dot spin and the electrons does not become too large. It is found that this is the case if the bias voltage applied to the dot or the temperature is much larger than the Kondo temperature.It is shown that the RG equations lead to coupling functions depending on the energy of the incident and the scattered electron. This dependence is simplified to the dependence on the average energy of the incident and the scattered electron in the weak coupling regime.A closer inspection of the RG equations for the couplings shows that they not only contain terms leading to a logarithmic increase of the coupling strength but also subleading terms which are unimportant in the weak coupling regime. The RG equations are solved analytically for the case where the subleading terms are neglected. It is found that the coupling functions develop resonances at energies associated with the chemical potentials of the reservoirs. The height of these resonances is determined by the decoherence rate. However, for bias voltages much larger than the Kondo temperature these resonances do not affect the conductance for a two terminal setup, so the conductance does not provide a measure for the decoherence rate in this case. It is found that the conductance for the two terminal setup is a universal function of voltage/Kondo temperature and temperature/Kondo temperature in the weak coupling regime. This function is calculated analytically in the limits of vanishing voltage or temperature. Furthermore, it is found that for zero temperature the decoherence rate is proportional to the non-equilibrium current flowing through the system. This supports the physical picture that a non-equilibrium current is associated with a decoherence rate, which cuts the logarithmic scaling of the coupling strength.Finally, the RG equations are solved numerically for zero temperature including the subleading corrections found before. As expected, the subleading corrections have no impact on the decoherence rate and the conductance in the weak coupling regime. On the other hand, it is found that they lead to strong deviations from the results of the analytic leading logarithmic approximation if the bias voltage approaches the Kondo temperature. Generally, it is found that the subleading corrections lead to a saturation of the decoherence rate and the conductance to some finite value. However, for the numerical solution a numerical regularization parameter has to be introduced. It turns out that the saturation value found depends on its value, so no statements on the precise value of the decoherence rate and the current can be made in this regime.

•Jan 1, 2006•0 citations

Renormalization-Group Theory for Quantum Dissipative Systems in Nonequilibrium

Keil, Markus

We apply the recently developed real-time renormalization-group (RTRG) method to different quantum dissipative systems. In contrast to other renormalization-group techniques the RTRG approach is formulated on the Keldysh contour and therefore allows for the description of nonequilibrium situations. The method is based on an exact kinetic equation which determines the time-dependent reduced density matrix. The kinetic equation involves an integral kernel which describes the phenomenon of dissipation. We calculate this quantity using the RTRG method: After applying Wick's theorem the arising diagrams are accounted for successively with respect to their time scale. During that process the theory is left invariant by renormalizing the operators describing the theory. As a starting point we discuss the polaron problem. We calculate the renormalized ground-state energy and the effective mass. We compare our results of the RTRG with those obtained from the flow equation formalism, which is based on infinitesimal transformations. In the central part of this thesis we consider the spin-boson model which is applicable to various physical and chemical systems. It models a particle in a double-well potential which is coupled to a bosonic environment. We determine static properties such as the static susceptibility. We also present results for the time evolution of the reduced density matrix starting from an arbitrary nonequilibrium state. The asymptotic behaviour of the reduced density matrix can be characterized by an oscillation frequency and two decay constants, all of which are determined. Finally, we calculate equilibrium correlation functions. We find consistency of our results with quantum Monte Carlo simulations and with the Shiba-relation. Additionally, our results show the correct scaling behaviour. Comparisons are also drawn to results of the noninteracting blip approximation and to results obtained for the anisotropic Kondo model. Our approach is very accurate for not too large couplings to the environment. Furthermore, the stationary tunnel current through a double quantum dot coupled to a phonon bath is calculated and compared to experimental data. Eventually, the running couplings of the two-lead Kondo model are discussed.

•Jan 1, 2002•0 citations

Renormalization group approach to time-dependent transport through correlated quantum dots

Dante Marvin Kennes, Severin G. Jakobs, Christoph Karrasch +1 more

We introduce a real time version of the functional renormalization group, which allows us to study correlation effects on nonequilibrium transport through quantum dots. Our method is equally capable to address (i) the relaxation out of a nonequilibrium initial state into a (potentially) steady state driven by a bias voltage and (ii) the dynamics governed by an explicitly time-dependent Hamiltonian. All time regimes from transient to asymptotic can be tackled; the only approximation is the consistent truncation of the flow equations at a given order. As an application we investigate the relaxation dynamics of the interacting resonant level model, which describes a fermionic quantum dot dominated by charge fluctuations. Moreover, we study incoherence and relaxation phenomena within the ohmic spin-boson model by mapping the latter to the interacting resonant level model.

arXiv:1111.6982•Physical Review B•Feb 21, 2012•77 citations

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

Frithjof B. Anders, Sebastian Schmitt

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 .

Journal of Physics Conference Series•Apr 1, 2010•6 citations

Renormalization group analysis of the interacting resonant-level model at finite bias: Generic analytic study of static properties and quench dynamics

Sabine Andergassen, Mikhail Pletyukhov, Dirk Schuricht +2 more

Using a real-time renormalization group method we study the minimal model of a quantum dot dominated by charge fluctuations, the two-lead interacting resonant level model, at finite bias voltage. We develop a set of RG equations to treat the case of weak and strong charge fluctuations, together with the determination of power-law exponents up to second order in the Coulomb interaction. We derive analytic expressions for the charge susceptibility, the steady-state current, and the conductance in the situation of arbitrary system parameters, in particular away from the particle-hole symmetric point and for asymmetric Coulomb interactions. In the asymmetric situation we find that power laws can be observed for the current only as a function of the level position (gate voltage) but in general not as a function of the voltage except for extremely large voltages. Furthermore, we study the quench dynamics after sudden switch-on of the level-lead couplings. The time evolution of the dot occupation and current is governed by exponential relaxation accompanied by voltage-dependent oscillations and characteristic algebraic decay.

arXiv:1010.5666•Physical Review B•May 6, 2011•69 citations
We present the real-time renormalization group (RTRG) method as a method to describe the stationary state current through generic multilevel quantum dots in nonequilibrium. The employed approach consists of a very rudimentary approximation for the renormalization group (RG) equations which neglects all vertex corrections while it provides a means to compute the effective dot Liouvillian self-consistently. Being based on a weak-coupling expansion in the tunneling between dot and reservoirs, the RTRG approach turns out to reliably describe charge fluctuations in and out of equilibrium for arbitrary coupling strength, even at zero temperature. We confirm this in the linear response regime with a benchmark against highly accurate numerical renormalization group data in the exemplary case of three-level quantum dots. For small to intermediate bias voltages and weak Coulomb interactions, we find an excellent agreement between RTRG and functional renormalization group data, which can be expected to be accurate in this regime. As a consequence, we advertise the presented RTRG approach as an efficient and versatile tool to describe charge fluctuations in quantum dot systems.
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.
Correlations are a central ingredient of modern condensed matter physics. With the advent of mesoscopic physics, it became possible to study correlated systems in non-equilibrium situations, which can be reached by e.g. driving a constant current through a quantum dot. A basic question is how correlations on the one hand and a non-equilibrium current on the other hand influence each other. Here this interplay is studied for the model system of a quantum dot in the Kondo regime using the real time renormalization group (RTRG) method. This method is a functional renormalization group (RG) method allowing a non-perturbative description of correlated systems in non-equilibrium. For the Kondo quantum dot the method is generalized to the case of double vertices. The RTRG method not only provides equations for the calculation of the renormalized coupling between the quantum dot spin and the reservoir electrons, but also equations for the calculation of the non-equilibrium decoherence rate and the non-equilibrium current.The RG equations are derived to second order. This approximation is expected to be valid as long as the quantum dot does not enter the strong coupling regime, i.e. the effective coupling strength between the quantum dot spin and the electrons does not become too large. It is found that this is the case if the bias voltage applied to the dot or the temperature is much larger than the Kondo temperature.It is shown that the RG equations lead to coupling functions depending on the energy of the incident and the scattered electron. This dependence is simplified to the dependence on the average energy of the incident and the scattered electron in the weak coupling regime.A closer inspection of the RG equations for the couplings shows that they not only contain terms leading to a logarithmic increase of the coupling strength but also subleading terms which are unimportant in the weak coupling regime. The RG equations are solved analytically for the case where the subleading terms are neglected. It is found that the coupling functions develop resonances at energies associated with the chemical potentials of the reservoirs. The height of these resonances is determined by the decoherence rate. However, for bias voltages much larger than the Kondo temperature these resonances do not affect the conductance for a two terminal setup, so the conductance does not provide a measure for the decoherence rate in this case. It is found that the conductance for the two terminal setup is a universal function of voltage/Kondo temperature and temperature/Kondo temperature in the weak coupling regime. This function is calculated analytically in the limits of vanishing voltage or temperature. Furthermore, it is found that for zero temperature the decoherence rate is proportional to the non-equilibrium current flowing through the system. This supports the physical picture that a non-equilibrium current is associated with a decoherence rate, which cuts the logarithmic scaling of the coupling strength.Finally, the RG equations are solved numerically for zero temperature including the subleading corrections found before. As expected, the subleading corrections have no impact on the decoherence rate and the conductance in the weak coupling regime. On the other hand, it is found that they lead to strong deviations from the results of the analytic leading logarithmic approximation if the bias voltage approaches the Kondo temperature. Generally, it is found that the subleading corrections lead to a saturation of the decoherence rate and the conductance to some finite value. However, for the numerical solution a numerical regularization parameter has to be introduced. It turns out that the saturation value found depends on its value, so no statements on the precise value of the decoherence rate and the current can be made in this regime.
We apply the recently developed real-time renormalization-group (RTRG) method to different quantum dissipative systems. In contrast to other renormalization-group techniques the RTRG approach is formulated on the Keldysh contour and therefore allows for the description of nonequilibrium situations. The method is based on an exact kinetic equation which determines the time-dependent reduced density matrix. The kinetic equation involves an integral kernel which describes the phenomenon of dissipation. We calculate this quantity using the RTRG method: After applying Wick's theorem the arising diagrams are accounted for successively with respect to their time scale. During that process the theory is left invariant by renormalizing the operators describing the theory. As a starting point we discuss the polaron problem. We calculate the renormalized ground-state energy and the effective mass. We compare our results of the RTRG with those obtained from the flow equation formalism, which is based on infinitesimal transformations. In the central part of this thesis we consider the spin-boson model which is applicable to various physical and chemical systems. It models a particle in a double-well potential which is coupled to a bosonic environment. We determine static properties such as the static susceptibility. We also present results for the time evolution of the reduced density matrix starting from an arbitrary nonequilibrium state. The asymptotic behaviour of the reduced density matrix can be characterized by an oscillation frequency and two decay constants, all of which are determined. Finally, we calculate equilibrium correlation functions. We find consistency of our results with quantum Monte Carlo simulations and with the Shiba-relation. Additionally, our results show the correct scaling behaviour. Comparisons are also drawn to results of the noninteracting blip approximation and to results obtained for the anisotropic Kondo model. Our approach is very accurate for not too large couplings to the environment. Furthermore, the stationary tunnel current through a double quantum dot coupled to a phonon bath is calculated and compared to experimental data. Eventually, the running couplings of the two-lead Kondo model are discussed.
We introduce a real time version of the functional renormalization group, which allows us to study correlation effects on nonequilibrium transport through quantum dots. Our method is equally capable to address (i) the relaxation out of a nonequilibrium initial state into a (potentially) steady state driven by a bias voltage and (ii) the dynamics governed by an explicitly time-dependent Hamiltonian. All time regimes from transient to asymptotic can be tackled; the only approximation is the consistent truncation of the flow equations at a given order. As an application we investigate the relaxation dynamics of the interacting resonant level model, which describes a fermionic quantum dot dominated by charge fluctuations. Moreover, we study incoherence and relaxation phenomena within the ohmic spin-boson model by mapping the latter to the interacting resonant level model.
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 .
Using a real-time renormalization group method we study the minimal model of a quantum dot dominated by charge fluctuations, the two-lead interacting resonant level model, at finite bias voltage. We develop a set of RG equations to treat the case of weak and strong charge fluctuations, together with the determination of power-law exponents up to second order in the Coulomb interaction. We derive analytic expressions for the charge susceptibility, the steady-state current, and the conductance in the situation of arbitrary system parameters, in particular away from the particle-hole symmetric point and for asymmetric Coulomb interactions. In the asymmetric situation we find that power laws can be observed for the current only as a function of the level position (gate voltage) but in general not as a function of the voltage except for extremely large voltages. Furthermore, we study the quench dynamics after sudden switch-on of the level-lead couplings. The time evolution of the dot occupation and current is governed by exponential relaxation accompanied by voltage-dependent oscillations and characteristic algebraic decay.