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:1406.0254 · Physical Review B

From thermal equilibrium to nonequilibrium quench dynamics: A conserving approximation for the interacting resonant level

Yuval Vinkler-Aviv, Avraham Schiller, Frithjof B. Anders

Published Oct 8, 2014• 14 citations
arXivHTMLPDF

Abstract

We develop a low-order conserving approximation for the interacting resonant-level model (IRLM), and apply it to (i) thermal equilibrium, (ii) nonequilibrium steady state, and (iii) nonequilibrium quench dynamics. Thermal equilibrium is first used to carefully gauge the quality of the approximation by comparing the results with other well-studied methods, and finding good agreement for small values of the interaction. We analytically show that the power-law exponent of the renormalized level width usually derived using renormalization group approaches can also be correctly obtained in our approach in the weak interaction limit. A closed expression for the nonequilibrium steady-state current is derived and analytically and numerically evaluated. We find a negative differential conductance at large voltages, and the exponent of the power-law suppression of the steady-state current is calculated analytically at zero temperature. The response of the system to quenches is investigated for a single lead as well as for two-lead setup at finite voltage bias at particle-hole symmetry using a self-consistent two-times Keldysh Green function approach, and results are presented for the time-dependent current for different bias and contact interaction strength.

DOI
10.1103/physrevb.90.155110

Related papers

Functional renormalization group study of the interacting resonant level model in and out of equilibrium

Christoph Karrasch, Mikhail Pletyukhov, L. Borda +1 more

We investigate equilibrium and steady-state nonequilibrium transport properties of a spinless resonant level locally coupled to two conduction bands of width $\ensuremath{\sim}\ensuremath{\Gamma}$ via a Coulomb interaction $U$ and a hybridization ${t}^{\ensuremath{'}}$. In order to study the effects of finite bias voltages beyond linear response, a generalization of the functional renormalization group to Keldysh frequency space is employed. Being mostly unexplored in the context of quantum impurity systems out of equilibrium, we benchmark this method against recently published time-dependent density matrix renormalization group data. We thoroughly investigate the scaling limit $\ensuremath{\Gamma}\ensuremath{\rightarrow}\ensuremath{\infty}$ characterized by the appearance of power laws. Most importantly, at the particle-hole symmetric point the steady-state current decays like $J\ensuremath{\sim}{V}^{\ensuremath{-}{\ensuremath{\alpha}}_{J}}$ as a function of the bias voltage $V⪢{t}^{\ensuremath{'}}$, with an exponent ${\ensuremath{\alpha}}_{J}(U)$ that we calculate to leading order in the Coulomb interaction strength. In contrast, we do not observe a pure power-law (but more complex) current-voltage-relation if the energy $ϵ$ of the resonant level is pinned close to either one of the chemical potentials $\ifmmode\pm\else\textpm\fi{}V/2$.

arXiv:0911.5165•Physical Review B•Mar 24, 2010•70 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

Nonequilibrium Green's functions and their relation to the negative differential conductance in the interacting resonant level model

Max Erich Sorantin, Wolfgang von der Linden, Roman Lucrezi +1 more

We evaluate the nonequilibrium single-particle Green's functions in the steady state of the interacting resonant level model (IRLM) under the effect of an applied bias voltage. Employing the so-called auxiliary master equation approach, we present accurate nonperturbative results for the nonequilibrium spectral and effective distribution functions, as well as for the current-voltage characteristics. We find a drastic change of these spectral properties between the regimes of low- and high-bias voltages and discuss the relation of these changes to the negative differential conductance (NDC), a prominent feature in the nonequilibrium IRLM. The anomalous evolution of the effective distribution function next to the impurity shown by our calculations suggests a mechanism whereby the impurity gets effectively decoupled from the leads at voltages where the NDC sets in, in agreement with previous renormalization group approaches. This scenario is qualitatively confirmed by a Hartree-Fock treatment of the model.

arXiv:1812.08618•Physical review. B./Physical review. B•Feb 20, 2019•2 citations

Nonequilibrium Steady-State Transport in Quantum Impurity Models: A Thermofield and Quantum Quench Approach Using Matrix Product States

Frauke Schwarz, Ireneusz Weymann, Jan von Delft +1 more

The numerical renormalization group (NRG) is tailored to describe interacting impurity models in equilibrium, but it faces limitations for steady-state nonequilibrium, arising, e.g., due to an applied bias voltage. We show that these limitations can be overcome by describing the thermal leads using a thermofield approach, integrating out high energy modes using NRG, and then treating the nonequilibrium dynamics at low energies using a quench protocol, implemented using the time-dependent density matrix renormalization group. This yields quantitatively reliable results for the current (with errors ≲3%) down to the exponentially small energy scales characteristic of impurity models. We present results of benchmark quality for the temperature and magnetic field dependence of the zero-bias conductance peak for the single-impurity Anderson model.

arXiv:1708.06315•Physical Review Letters•Sep 26, 2018•81 citations

Nonequilibrium functional renormalization group with frequency-dependent vertex function: A study of the single-impurity Anderson model

Severin G. Jakobs, Mikhail Pletyukhov, Herbert Schoeller

We investigate nonequilibrium properties of the single-impurity Anderson model by means of the functional renormalization group (fRG) within Keldysh formalism. We present how the level broadening $\ensuremath{\Gamma}/2$ can be used as flow parameter for the fRG. This choice preserves important aspects of the Fermi-liquid behavior that the model exhibits in case of particle-hole symmetry. An approximation scheme for the Keldysh fRG is developed which accounts for the frequency dependence of the two-particle vertex in a way similar but not equivalent to a recently published approximation to the equilibrium Matsubara fRG. Our method turns out to be a flexible tool for the study of weak to intermediate on-site interactions $U\ensuremath{\lesssim}3\ensuremath{\Gamma}$. In equilibrium we find excellent agreement with numerical RG results for the linear conductance at finite gate voltage, magnetic field, and temperature. In nonequilibrium, our results for the current agree well with time-dependent density-matrix RG. For the nonlinear conductance as function of the bias voltage, we propose reliable results at finite magnetic field and finite temperature. Furthermore, we demonstrate the exponentially small scale of the Kondo temperature to appear in the second-order derivative of the self-energy. We show that the approximation is, however, not able to reproduce the scaling of the effective mass at large interactions.

arXiv:0911.5502•Physical Review B•May 12, 2010•119 citations

New Method for Studying Steady States in Quantum Impurity Problems: The Interacting Resonant Level Model

Benjamin Doyon

We develop a new perturbative method for studying any steady states of quantum impurities, in or out of equilibrium. We show that steady-state averages are completely fixed by basic properties of the steady-state (Hershfield's) density matrix along with dynamical "impurity conditions." This gives the full perturbative expansion without Feynman diagrams (matrix products instead are used), and "resums" into an equilibrium average that may lend itself to numerical procedures. We calculate the universal current in the interacting resonant level model (IRLM) at finite bias V to first order in Coulomb repulsion U for all V and temperatures. We find that the bias, like the temperature, cuts off low-energy processes. In the IRLM, this implies a power-law decay of the current at large V (also recently observed by Boulat and Saleur at some finite value of U).

arXiv:cond-mat/0703249•Physical Review Letters•Aug 17, 2007•71 citations
We develop a low-order conserving approximation for the interacting resonant-level model (IRLM), and apply it to (i) thermal equilibrium, (ii) nonequilibrium steady state, and (iii) nonequilibrium quench dynamics. Thermal equilibrium is first used to carefully gauge the quality of the approximation by comparing the results with other well-studied methods, and finding good agreement for small values of the interaction. We analytically show that the power-law exponent of the renormalized level width usually derived using renormalization group approaches can also be correctly obtained in our approach in the weak interaction limit. A closed expression for the nonequilibrium steady-state current is derived and analytically and numerically evaluated. We find a negative differential conductance at large voltages, and the exponent of the power-law suppression of the steady-state current is calculated analytically at zero temperature. The response of the system to quenches is investigated for a single lead as well as for two-lead setup at finite voltage bias at particle-hole symmetry using a self-consistent two-times Keldysh Green function approach, and results are presented for the time-dependent current for different bias and contact interaction strength.
We investigate equilibrium and steady-state nonequilibrium transport properties of a spinless resonant level locally coupled to two conduction bands of width $\ensuremath{\sim}\ensuremath{\Gamma}$ via a Coulomb interaction UUU and a hybridization ${t}^{\ensuremath{'}}$. In order to study the effects of finite bias voltages beyond linear response, a generalization of the functional renormalization group to Keldysh frequency space is employed. Being mostly unexplored in the context of quantum impurity systems out of equilibrium, we benchmark this method against recently published time-dependent density matrix renormalization group data. We thoroughly investigate the scaling limit $\ensuremath{\Gamma}\ensuremath{\rightarrow}\ensuremath{\infty}$ characterized by the appearance of power laws. Most importantly, at the particle-hole symmetric point the steady-state current decays like $J\ensuremath{\sim}{V}^{\ensuremath{-}{\ensuremath{\alpha}}_{J}}$ as a function of the bias voltage $V⪢{t}^{\ensuremath{'}}$, with an exponent ${\ensuremath{\alpha}}_{J}(U)$ that we calculate to leading order in the Coulomb interaction strength. In contrast, we do not observe a pure power-law (but more complex) current-voltage-relation if the energy ϵϵϵ of the resonant level is pinned close to either one of the chemical potentials $\ifmmode\pm\else\textpm\fi{}V/2$.
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.
We evaluate the nonequilibrium single-particle Green's functions in the steady state of the interacting resonant level model (IRLM) under the effect of an applied bias voltage. Employing the so-called auxiliary master equation approach, we present accurate nonperturbative results for the nonequilibrium spectral and effective distribution functions, as well as for the current-voltage characteristics. We find a drastic change of these spectral properties between the regimes of low- and high-bias voltages and discuss the relation of these changes to the negative differential conductance (NDC), a prominent feature in the nonequilibrium IRLM. The anomalous evolution of the effective distribution function next to the impurity shown by our calculations suggests a mechanism whereby the impurity gets effectively decoupled from the leads at voltages where the NDC sets in, in agreement with previous renormalization group approaches. This scenario is qualitatively confirmed by a Hartree-Fock treatment of the model.
The numerical renormalization group (NRG) is tailored to describe interacting impurity models in equilibrium, but it faces limitations for steady-state nonequilibrium, arising, e.g., due to an applied bias voltage. We show that these limitations can be overcome by describing the thermal leads using a thermofield approach, integrating out high energy modes using NRG, and then treating the nonequilibrium dynamics at low energies using a quench protocol, implemented using the time-dependent density matrix renormalization group. This yields quantitatively reliable results for the current (with errors ≲3%) down to the exponentially small energy scales characteristic of impurity models. We present results of benchmark quality for the temperature and magnetic field dependence of the zero-bias conductance peak for the single-impurity Anderson model.
We investigate nonequilibrium properties of the single-impurity Anderson model by means of the functional renormalization group (fRG) within Keldysh formalism. We present how the level broadening $\ensuremath{\Gamma}/2$ can be used as flow parameter for the fRG. This choice preserves important aspects of the Fermi-liquid behavior that the model exhibits in case of particle-hole symmetry. An approximation scheme for the Keldysh fRG is developed which accounts for the frequency dependence of the two-particle vertex in a way similar but not equivalent to a recently published approximation to the equilibrium Matsubara fRG. Our method turns out to be a flexible tool for the study of weak to intermediate on-site interactions $U\ensuremath{\lesssim}3\ensuremath{\Gamma}$. In equilibrium we find excellent agreement with numerical RG results for the linear conductance at finite gate voltage, magnetic field, and temperature. In nonequilibrium, our results for the current agree well with time-dependent density-matrix RG. For the nonlinear conductance as function of the bias voltage, we propose reliable results at finite magnetic field and finite temperature. Furthermore, we demonstrate the exponentially small scale of the Kondo temperature to appear in the second-order derivative of the self-energy. We show that the approximation is, however, not able to reproduce the scaling of the effective mass at large interactions.
We develop a new perturbative method for studying any steady states of quantum impurities, in or out of equilibrium. We show that steady-state averages are completely fixed by basic properties of the steady-state (Hershfield's) density matrix along with dynamical "impurity conditions." This gives the full perturbative expansion without Feynman diagrams (matrix products instead are used), and "resums" into an equilibrium average that may lend itself to numerical procedures. We calculate the universal current in the interacting resonant level model (IRLM) at finite bias V to first order in Coulomb repulsion U for all V and temperatures. We find that the bias, like the temperature, cuts off low-energy processes. In the IRLM, this implies a power-law decay of the current at large V (also recently observed by Boulat and Saleur at some finite value of U).