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

Extension of the iterated perturbation theory at arbitrary fillings to nonequilibrium steady states

Anonymous, Enrico Arrigoni

Published Aug 31, 2026• 0 citations
arXivHTMLPDF

Abstract

We extend the Kajueter-Kotliar [] iterated perturbation theory (KK-IPT) away from half-filling to nonequilibrium steady states. We benchmark the resulting nonequilibrium KK-IPT approach against the auxiliary master equation approach (AMEA), whose accuracy is controlled in and out of equilibrium. As expected, in equilibrium, KK-IPT reproduces the AMEA results for different fillings with high accuracy at the level of both spectral properties and electron densities. Out of equilibrium, we study quantum transport across a correlated impurity and compute the differential conductance and spectral functions. We find very good agreement between nonequilibrium KK-IPT and AMEA in the parameter regime where the latter is reliable, in particular at moderate temperatures and biases. Although a controlled benchmark is not available in the low-temperature, low-bias regime, where AMEA becomes less reliable, we show that this nonequilibrium KK-IPT impurity solver satisfies the exact spectral sum rules for the first and second moments to high accuracy throughout the entire parameter range studied. These results support nonequilibrium KK-IPT as an approximate description of nonequilibrium steady states away from half-filling. At the same time, comparing against AMEA the double occupancy obtained from the nonequilibrium Galitskii-Migdal expression for the interaction energy shows that the deviation from AMEA remains small near half-filling for moderate and large values of the bias, but grows markedly away from half-filling, delineating the regime in which the method can be trusted quantitatively rather than merely qualitatively.

DOI
10.1103/lwg4-6tq8

Related papers

Configuration interaction based nonequilibrium steady state impurity solver

Daniel Werner, Jan Lotze, Enrico Arrigoni

We present a solver for correlated impurity problems out of equilibrium based on a combination of the so-called auxiliary master equation approach (AMEA) and the configuration interaction expansion. Within AMEA one maps the original impurity model onto an auxiliary open quantum system with a restricted number of bath sites which can be addressed by numerical many-body approaches such as Lanczos/Arnoldi exact diagonalization (ED) or matrix product states (MPS). While the mapping becomes exponentially more accurate with increasing number of bath sites, ED implementations are severely limited due to the fast increase of the Hilbert space dimension for open systems, and the MPS solver typically requires rather long runtimes. Here, we propose to adopt a configuration interaction approach augmented by active space extension to solve numerically the correlated auxiliary open quantum system. This allows access to a larger number of bath sites at lower computational costs than for plain ED. We benchmark the approach with numerical renormalization group results in equilibrium and with MPS out of equilibrium. In particular, we evaluate the current, the conductance, as well as the Kondo peak and its splitting as a function of increasing bias voltage below the Kondo temperature ${T}_{\text{K}}$. We obtain a rather accurate scaling of the conductance as a function of the bias voltage and temperature rescaled by ${T}_{\text{K}}$ for moderate to strong interactions in a wide range of parameters. The approach combines the fast runtime of ED with an accuracy close to the one achieved by MPS making it an attractive solver for nonequilibrium dynamical mean field theory.

arXiv:2210.09623•Physical review. B./Physical review. B•Feb 8, 2023•33 citations

Transport through a correlated interface: Auxiliary master equation approach

Irakli Titvinidze, Antonius Dorda, Wolfgang von der Linden +1 more

We present improvements of a recently introduced numerical method [E. Arrigoni et al., Phys. Rev. Lett. 110, 086403 (2013)] to compute steady-state properties of strongly correlated electronic systems out of equilibrium. The method can be considered as a nonequilibrium generalization of exact diagonalization based dynamical mean-field theory (DMFT). The key modification for the nonequilibrium situation consists in addressing the DMFT impurity problem within an auxiliary system consisting of the correlated impurity, ${N}_{b}$ uncorrelated bath sites, and two Markovian environments (sink and reservoir). Algorithmic improvements in the impurity solver allow to treat efficiently larger values of ${N}_{b}$ than previously in DMFT. This increases the accuracy of the results and is crucial for a correct description of the physical behavior of the system in the relevant parameter range including a semiquantitative description of the Kondo regime. To illustrate the approach, we consider a monoatomic layer of correlated orbitals, described by the single-band Hubbard model, attached to two metallic leads. The nonequilibrium situation is driven by a bias voltage applied to the leads. For this system, we investigate the spectral function and the steady-state current-voltage characteristics in the weakly as well as in the strongly interacting limit. In particular, we investigate the nonequilibrium behavior of quasiparticle excitations within the Mott gap of the correlated layer. We find for low-bias voltage Kondo-type behavior in the vicinity of the insulating phase. In particular, we observe a splitting of the Kondo resonance as a function of the bias voltage.

arXiv:1508.02953•Physical Review B•Dec 21, 2015•24 citations

Auxiliary master equation approach within matrix product states: Spectral properties of the nonequilibrium Anderson impurity model

Antonius Dorda, Martin Ganahl, Hans Gerd Evertz +2 more

Within the recently introduced auxiliary master equation approach it is possible to address steady state properties of strongly correlated impurity models, small molecules or clusters efficiently and with high accuracy. It is particularly suited for dynamical mean field theory in the nonequilibrium as well as in the equilibrium case. The method is based on the solution of an auxiliary open quantum system, which can be made quickly equivalent to the original impurity problem. In its first implementation a Krylov space method was employed. Here, we aim at extending the capabilities of the approach by adopting matrix product states for the solution of the corresponding auxiliary quantum master equation. This allows for a drastic increase in accuracy and permits us to access the Kondo regime for large values of the interaction. In particular, we investigate the nonequilibrium steady state of a single impurity Anderson model and focus on the spectral properties for temperatures T below the Kondo temperature T_K and for small bias voltages phi. For the two cases considered, with T=T_K/4 and T=T_K/10 we find a clear splitting of the Kondo resonance into a two-peak structure for phi close above T_K. In the equilibrium case (phi=0) and for T=T_K/4, the obtained spectral function essentially coincides with the one from numerical renormalization group.

arXiv:1507.02982•Physical Review B•Sep 24, 2015•79 citations

Benchmark of a modified iterated perturbation theory approach on the fcc lattice at strong coupling

Louis-François Arsenault, P. Sémon, A.–M. S. Tremblay

The dynamical mean-field theory approach to the Hubbard model requires a method to solve the problem of a quantum impurity in a bath of noninteracting electrons. Iterated perturbation theory (IPT) has proven its effectiveness as a solver in many cases of interest. Based on general principles and on comparisons with an essentially exact continuous-time quantum Monte Carlo (CTQMC) solver, here we show that the standard implementation of IPT fails away from half-filling when the interaction strength is much larger than the bandwidth. We propose a slight modification to the IPT algorithm that replaces one of the equations by the requirement that double occupancy calculated with IPT gives the correct value. We call this method IPT-$D$. We recover the Fermi liquid ground state away from half-filling. The Fermi liquid parameters, density of states, chemical potential, energy, and specific heat on the fcc lattice are calculated with both IPT-$D$ and CTQMC as benchmark examples. We also calculated the resistivity and the optical conductivity within IPT-$D$. Particle-hole asymmetry persists even at coupling twice the bandwidth. A generalization to the multiorbital case is suggested. Several algorithms that speed up the calculations are described in appendixes.

arXiv:1202.5814•Physical Review B•Aug 22, 2012•44 citations

Auxiliary master equation approach to nonequilibrium correlated impurities

Antonius Dorda, Martin Nuss, Wolfgang von der Linden +1 more

We present a numerical method for the study of correlated quantum impurity problems out of equilibrium, which is particularly suited to address steady-state properties within dynamical mean field theory. The approach, recently introduced by Arrigoni et al. [Phys. Rev. Lett. 110, 086403 (2013)], is based upon a mapping of the original impurity problem onto an auxiliary open quantum system, consisting of the interacting impurity coupled to bath sites as well as to a Markovian environment. The dynamics of the auxiliary system is governed by a Lindblad master equation whose parameters are used to optimize the mapping. The accuracy of the results can be readily estimated and systematically improved by increasing the number of auxiliary bath sites, or by introducing a linear correction. Here, we focus on a detailed discussion of the proposed approach including technical remarks. To solve for the Green's functions of the auxiliary impurity problem, a non-Hermitian Lanczos diagonalization is applied. As a benchmark, results for the steady-state current-voltage characteristics of the single-impurity Anderson model are presented. Furthermore, the bias dependence of the single-particle spectral function and the splitting of the Kondo resonance are discussed. In its present form, the method is fast, efficient, and features a controlled accuracy.

arXiv:1312.4586•Physical Review B•Apr 7, 2014•117 citations

Nonequilibrium Dynamical Mean-Field Theory: An Auxiliary Quantum Master Equation Approach

Enrico Arrigoni, Michael Knap, Wolfgang von der Linden

We introduce a versatile method to compute electronic steady-state properties of strongly correlated extended quantum systems out of equilibrium. The approach is based on dynamical mean-field theory (DMFT), in which the original system is mapped onto an auxiliary nonequilibrium impurity problem imbedded in a Markovian environment. The steady-state Green's function of the auxiliary system is solved by full diagonalization of the corresponding Lindblad equation. The approach can be regarded as the nontrivial extension of the exact-diagonalization-based DMFT to the nonequilibrium case. As a first application, we consider an interacting Hubbard layer attached to two metallic leads and present results for the steady-state current and the nonequilibrium density of states.

arXiv:1210.4167•Physical Review Letters•Feb 20, 2013•141 citations
We extend the Kajueter-Kotliar [] iterated perturbation theory (KK-IPT) away from half-filling to nonequilibrium steady states. We benchmark the resulting nonequilibrium KK-IPT approach against the auxiliary master equation approach (AMEA), whose accuracy is controlled in and out of equilibrium. As expected, in equilibrium, KK-IPT reproduces the AMEA results for different fillings with high accuracy at the level of both spectral properties and electron densities. Out of equilibrium, we study quantum transport across a correlated impurity and compute the differential conductance and spectral functions. We find very good agreement between nonequilibrium KK-IPT and AMEA in the parameter regime where the latter is reliable, in particular at moderate temperatures and biases. Although a controlled benchmark is not available in the low-temperature, low-bias regime, where AMEA becomes less reliable, we show that this nonequilibrium KK-IPT impurity solver satisfies the exact spectral sum rules for the first and second moments to high accuracy throughout the entire parameter range studied. These results support nonequilibrium KK-IPT as an approximate description of nonequilibrium steady states away from half-filling. At the same time, comparing against AMEA the double occupancy obtained from the nonequilibrium Galitskii-Migdal expression for the interaction energy shows that the deviation from AMEA remains small near half-filling for moderate and large values of the bias, but grows markedly away from half-filling, delineating the regime in which the method can be trusted quantitatively rather than merely qualitatively.
We present a solver for correlated impurity problems out of equilibrium based on a combination of the so-called auxiliary master equation approach (AMEA) and the configuration interaction expansion. Within AMEA one maps the original impurity model onto an auxiliary open quantum system with a restricted number of bath sites which can be addressed by numerical many-body approaches such as Lanczos/Arnoldi exact diagonalization (ED) or matrix product states (MPS). While the mapping becomes exponentially more accurate with increasing number of bath sites, ED implementations are severely limited due to the fast increase of the Hilbert space dimension for open systems, and the MPS solver typically requires rather long runtimes. Here, we propose to adopt a configuration interaction approach augmented by active space extension to solve numerically the correlated auxiliary open quantum system. This allows access to a larger number of bath sites at lower computational costs than for plain ED. We benchmark the approach with numerical renormalization group results in equilibrium and with MPS out of equilibrium. In particular, we evaluate the current, the conductance, as well as the Kondo peak and its splitting as a function of increasing bias voltage below the Kondo temperature TK{T}_{\text{K}}TK​. We obtain a rather accurate scaling of the conductance as a function of the bias voltage and temperature rescaled by TK{T}_{\text{K}}TK​ for moderate to strong interactions in a wide range of parameters. The approach combines the fast runtime of ED with an accuracy close to the one achieved by MPS making it an attractive solver for nonequilibrium dynamical mean field theory.
We present improvements of a recently introduced numerical method [E. Arrigoni et al., Phys. Rev. Lett. 110, 086403 (2013)] to compute steady-state properties of strongly correlated electronic systems out of equilibrium. The method can be considered as a nonequilibrium generalization of exact diagonalization based dynamical mean-field theory (DMFT). The key modification for the nonequilibrium situation consists in addressing the DMFT impurity problem within an auxiliary system consisting of the correlated impurity, Nb{N}_{b}Nb​ uncorrelated bath sites, and two Markovian environments (sink and reservoir). Algorithmic improvements in the impurity solver allow to treat efficiently larger values of Nb{N}_{b}Nb​ than previously in DMFT. This increases the accuracy of the results and is crucial for a correct description of the physical behavior of the system in the relevant parameter range including a semiquantitative description of the Kondo regime. To illustrate the approach, we consider a monoatomic layer of correlated orbitals, described by the single-band Hubbard model, attached to two metallic leads. The nonequilibrium situation is driven by a bias voltage applied to the leads. For this system, we investigate the spectral function and the steady-state current-voltage characteristics in the weakly as well as in the strongly interacting limit. In particular, we investigate the nonequilibrium behavior of quasiparticle excitations within the Mott gap of the correlated layer. We find for low-bias voltage Kondo-type behavior in the vicinity of the insulating phase. In particular, we observe a splitting of the Kondo resonance as a function of the bias voltage.
Within the recently introduced auxiliary master equation approach it is possible to address steady state properties of strongly correlated impurity models, small molecules or clusters efficiently and with high accuracy. It is particularly suited for dynamical mean field theory in the nonequilibrium as well as in the equilibrium case. The method is based on the solution of an auxiliary open quantum system, which can be made quickly equivalent to the original impurity problem. In its first implementation a Krylov space method was employed. Here, we aim at extending the capabilities of the approach by adopting matrix product states for the solution of the corresponding auxiliary quantum master equation. This allows for a drastic increase in accuracy and permits us to access the Kondo regime for large values of the interaction. In particular, we investigate the nonequilibrium steady state of a single impurity Anderson model and focus on the spectral properties for temperatures T below the Kondo temperature T_K and for small bias voltages phi. For the two cases considered, with T=T_K/4 and T=T_K/10 we find a clear splitting of the Kondo resonance into a two-peak structure for phi close above T_K. In the equilibrium case (phi=0) and for T=T_K/4, the obtained spectral function essentially coincides with the one from numerical renormalization group.
The dynamical mean-field theory approach to the Hubbard model requires a method to solve the problem of a quantum impurity in a bath of noninteracting electrons. Iterated perturbation theory (IPT) has proven its effectiveness as a solver in many cases of interest. Based on general principles and on comparisons with an essentially exact continuous-time quantum Monte Carlo (CTQMC) solver, here we show that the standard implementation of IPT fails away from half-filling when the interaction strength is much larger than the bandwidth. We propose a slight modification to the IPT algorithm that replaces one of the equations by the requirement that double occupancy calculated with IPT gives the correct value. We call this method IPT-DDD. We recover the Fermi liquid ground state away from half-filling. The Fermi liquid parameters, density of states, chemical potential, energy, and specific heat on the fcc lattice are calculated with both IPT-DDD and CTQMC as benchmark examples. We also calculated the resistivity and the optical conductivity within IPT-DDD. Particle-hole asymmetry persists even at coupling twice the bandwidth. A generalization to the multiorbital case is suggested. Several algorithms that speed up the calculations are described in appendixes.
We present a numerical method for the study of correlated quantum impurity problems out of equilibrium, which is particularly suited to address steady-state properties within dynamical mean field theory. The approach, recently introduced by Arrigoni et al. [Phys. Rev. Lett. 110, 086403 (2013)], is based upon a mapping of the original impurity problem onto an auxiliary open quantum system, consisting of the interacting impurity coupled to bath sites as well as to a Markovian environment. The dynamics of the auxiliary system is governed by a Lindblad master equation whose parameters are used to optimize the mapping. The accuracy of the results can be readily estimated and systematically improved by increasing the number of auxiliary bath sites, or by introducing a linear correction. Here, we focus on a detailed discussion of the proposed approach including technical remarks. To solve for the Green's functions of the auxiliary impurity problem, a non-Hermitian Lanczos diagonalization is applied. As a benchmark, results for the steady-state current-voltage characteristics of the single-impurity Anderson model are presented. Furthermore, the bias dependence of the single-particle spectral function and the splitting of the Kondo resonance are discussed. In its present form, the method is fast, efficient, and features a controlled accuracy.
We introduce a versatile method to compute electronic steady-state properties of strongly correlated extended quantum systems out of equilibrium. The approach is based on dynamical mean-field theory (DMFT), in which the original system is mapped onto an auxiliary nonequilibrium impurity problem imbedded in a Markovian environment. The steady-state Green's function of the auxiliary system is solved by full diagonalization of the corresponding Lindblad equation. The approach can be regarded as the nontrivial extension of the exact-diagonalization-based DMFT to the nonequilibrium case. As a first application, we consider an interacting Hubbard layer attached to two metallic leads and present results for the steady-state current and the nonequilibrium density of states.