pith. sign in

arxiv: 2604.15942 · v1 · submitted 2026-04-17 · ❄️ cond-mat.str-el

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

Tommaso Maria Mazzocchi , Enrico Arrigoni This is my paper

Pith reviewed 2026-05-10 07:55 UTC · model grok-4.3

classification ❄️ cond-mat.str-el
keywords iterated perturbation theorynonequilibrium steady statescorrelated impurityquantum transportspectral functionsdifferential conductanceauxiliary master equation
0
0 comments X

The pith

The Kajueter-Kotliar iterated perturbation theory extends accurately to nonequilibrium steady states at arbitrary fillings.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper extends the Kajueter-Kotliar iterated perturbation theory (KK-IPT) to nonequilibrium steady states and to fillings away from half filling. It benchmarks the resulting approach against the auxiliary master equation approach (AMEA), whose accuracy is controlled. In equilibrium, the extension reproduces AMEA results for spectral properties and electron densities across different fillings. Out of equilibrium, for quantum transport through a correlated impurity, the differential conductance and spectral functions show very good agreement with AMEA at moderate temperatures and biases. These findings support nonequilibrium KK-IPT as a usable approximate method for steady states away from half filling, especially in regimes where AMEA becomes less reliable.

Core claim

We extend the Kajueter-Kotliar iterated perturbation theory (KK-IPT) away from half filling to nonequilibrium steady states. Benchmarking against the auxiliary master equation approach shows that the extended method reproduces spectral properties and electron densities with high accuracy in equilibrium for various fillings, and yields very good agreement for differential conductance and spectral functions out of equilibrium in the parameter regime where the benchmark is reliable, particularly at moderate temperatures and biases.

What carries the argument

The nonequilibrium extension of the Kajueter-Kotliar iterated perturbation theory (KK-IPT), which adapts the self-consistent perturbative treatment of the self-energy to steady-state conditions at arbitrary fillings.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The numerical stability in hard-to-benchmark regimes suggests the method could fill gaps left by approaches that lose reliability at low bias and temperature.
  • The extension may apply to other impurity or lattice models where nonequilibrium steady states at arbitrary filling need to be computed efficiently.

Load-bearing premise

The perturbative self-consistency loop of iterated perturbation theory remains accurate and stable when extended to nonequilibrium and to fillings away from half filling.

What would settle it

A controlled calculation or measurement of differential conductance through a correlated impurity in the low-temperature, low-bias regime that deviates significantly from the nonequilibrium KK-IPT prediction.

Figures

Figures reproduced from arXiv: 2604.15942 by Enrico Arrigoni, Tommaso Maria Mazzocchi.

Figure 1
Figure 1. Figure 1: FIG. 1. Spectral functions for equilibrium setup (a) E1 [PITH_FULL_IMAGE:figures/full_fig_p004_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. Imaginary part of the retarded SE profiles for the [PITH_FULL_IMAGE:figures/full_fig_p004_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. Nonequilibrium differential conductance, Eq. ( [PITH_FULL_IMAGE:figures/full_fig_p005_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. Nonequilibrium spectral function [PITH_FULL_IMAGE:figures/full_fig_p006_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5. Particle number, see Eq. ( [PITH_FULL_IMAGE:figures/full_fig_p008_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6. Double occupancy, see Eq. ( [PITH_FULL_IMAGE:figures/full_fig_p009_6.png] view at source ↗
read the original abstract

We extend the Kajueter-Kotliar [Phys. Rev. Lett. 77, 131 (1996)] 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. These results support nonequilibrium KK-IPT as an approximate description of nonequilibrium steady states away from half filling. Although a controlled benchmark is not available in the low-temperature, low-bias regime, where AMEA becomes less reliable, nonequilibrium KK-IPT remains numerically stable in those regions, suggesting that it may provide a useful alternative for nonequilibrium calculations in this regime.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 1 minor

Summary. The paper extends the Kajueter-Kotliar iterated perturbation theory (KK-IPT) from half-filling to arbitrary fillings and from equilibrium to nonequilibrium steady states. It benchmarks the resulting nonequilibrium KK-IPT against the auxiliary master equation approach (AMEA) for a correlated impurity in quantum transport, reporting good agreement on spectral functions, electron densities, and differential conductance at moderate temperatures and biases where AMEA is reliable. The authors note that nonequilibrium KK-IPT remains numerically stable in the low-T, low-bias regime (where AMEA loses reliability) and suggest it may serve as a useful approximate tool there.

Significance. If the accuracy of the extension holds beyond the benchmarked regimes, the method offers a computationally lightweight approximate solver for nonequilibrium steady states away from half-filling, where controlled methods are scarce. The explicit agreement with AMEA in moderate regimes and the absence of free parameters in the formulation are strengths that support its potential utility as a practical alternative for exploring correlated transport.

major comments (2)
  1. [Abstract and concluding discussion] The central suggestion that nonequilibrium KK-IPT 'may provide a useful alternative' in the low-temperature, low-bias regime rests only on numerical stability of the self-consistency loop, without any controlled benchmark against an exact or higher-order method. Because the paper itself states that AMEA becomes unreliable precisely in this regime, the extrapolation from moderate-T/bias agreement to this limit is not yet substantiated and weakens the claim of usefulness there.
  2. [Method section (extension procedure)] Implementation details of the nonequilibrium extension of the filling correction and the precise form of the perturbative self-energy in the steady-state Keldysh formalism are not provided. This absence prevents independent verification of whether the self-consistency loop and the arbitrary-filling adjustment preserve the original KK-IPT structure without introducing uncontrolled approximations.
minor comments (1)
  1. [Abstract] The distinction between demonstrated accuracy (moderate T/bias) and suggested applicability (low T/bias) could be made sharper in the abstract and conclusion to avoid over-interpretation of the stability result.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the thorough review and constructive feedback. We have revised the manuscript to address the concerns, moderating claims in the abstract and conclusion while adding the requested methodological details.

read point-by-point responses

  1. Referee: [Abstract and concluding discussion] The central suggestion that nonequilibrium KK-IPT 'may provide a useful alternative' in the low-temperature, low-bias regime rests only on numerical stability of the self-consistency loop, without any controlled benchmark against an exact or higher-order method. Because the paper itself states that AMEA becomes unreliable precisely in this regime, the extrapolation from moderate-T/bias agreement to this limit is not yet substantiated and weakens the claim of usefulness there.

    Authors: We agree that the original wording in the abstract and conclusion overstates the case by suggesting usefulness in the low-T, low-bias regime solely on the basis of numerical stability, without a controlled benchmark. As noted in the manuscript, no such benchmark is available because AMEA loses reliability there. We have therefore revised both the abstract and the concluding discussion to remove the suggestion of utility in that regime. The updated text now highlights the good agreement with AMEA at moderate temperatures and biases (where AMEA is reliable) and simply states that nonequilibrium KK-IPT remains numerically stable at low T and low bias, without claiming it provides a useful alternative there. This change makes the presentation accurate and avoids unsubstantiated extrapolation. revision: yes

  2. Referee: [Method section (extension procedure)] Implementation details of the nonequilibrium extension of the filling correction and the precise form of the perturbative self-energy in the steady-state Keldysh formalism are not provided. This absence prevents independent verification of whether the self-consistency loop and the arbitrary-filling adjustment preserve the original KK-IPT structure without introducing uncontrolled approximations.

    Authors: We thank the referee for identifying this omission. In the revised manuscript we have added a dedicated subsection to the Methods section that supplies the missing implementation details. We now explicitly give the form of the perturbative self-energy on the Keldysh contour for the nonequilibrium steady state, together with the precise adaptation of the filling correction (including the self-consistent determination of the chemical potential shift) that extends the original KK-IPT procedure away from half filling. These additions confirm that the self-consistency loop and the arbitrary-filling adjustment follow the same structure as the equilibrium KK-IPT without introducing further uncontrolled approximations. The new material should enable independent verification and reproduction of the results. revision: yes

Circularity Check

0 steps flagged

No significant circularity; derivation extends independent prior work and benchmarks against external AMEA method

full rationale

The paper extends the Kajueter-Kotliar iterated perturbation theory (KK-IPT, cited from 1996 independent authors) to nonequilibrium steady states and arbitrary fillings. It explicitly benchmarks spectral properties, densities, differential conductance, and spectral functions against the auxiliary master equation approach (AMEA), whose accuracy is stated to be controlled. No load-bearing self-citations, self-definitional equations, fitted parameters renamed as predictions, or ansatzes smuggled via prior author work appear in the abstract or described chain. The central support for the nonequilibrium KK-IPT extension rests on agreement with this independent benchmark in regimes where AMEA is reliable, rather than reducing to the method's own inputs by construction. The noted lack of controlled benchmark at low T/low bias is a limitation on validation range but does not constitute circularity in the derivation itself.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

Abstract-only review provides no explicit details on free parameters, axioms, or invented entities in the extension; the method relies on prior KK-IPT framework whose internal parameters are not specified here.

pith-pipeline@v0.9.0 · 5505 in / 1087 out tokens · 37912 ms · 2026-05-10T07:55:26.544755+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

47 extracted references · 47 canonical work pages

  1. [1]

    To recover the equilibrium condition it is enough to set µl =µ r =µ; in this case,µwill then be set to zero for simplicity. III. NONEQUILIBRIUM ITERA TED PER TURBA TION THEOR Y In this section, we briefly review KK-IPT for arbitrary fillings [26] and present its extension to nonequilibrium steady states, referred to here as nonequilibrium KK- IPT. The KK-...

    work page doi:10.55776/p33165

  2. [2]

    Georges and G

    A. Georges and G. Kotliar, Hubbard model in infinite dimensions, Phys. Rev. B45, 6479 (1992)

    work page 1992

  3. [3]

    Georges, G

    A. Georges, G. Kotliar, W. Krauth, and M. J. Rozen- berg, Dynamical mean-field theory of strongly correlated fermion systems and the limit of infinite dimensions, Rev. Mod. Phys.68, 13 (1996)

    work page 1996

  4. [4]

    H. Aoki, N. Tsuji, M. Eckstein, M. Kollar, T. Oka, and P. Werner, Nonequilibrium dynamical mean-field theory and its applications, Rev. Mod. Phys.86, 779 (2014)

    work page 2014

  5. [5]

    E. Gull, A. J. Millis, A. I. Lichtenstein, A. N. Rubtsov, M. Troyer, and P. Werner, Continuous-time monte carlo methods for quantum impurity models, Rev. Mod. Phys. 83, 349 (2011)

    work page 2011

  6. [6]

    Werner, A

    P. Werner, A. Comanac, L. de’ Medici, M. Troyer, and A. J. Millis, Continuous-time solver for quantum impu- rity models, Phys. Rev. Lett.97, 076405 (2006)

    work page 2006

  7. [7]

    K. G. Wilson, The renormalization group: Critical phe- nomena and the kondo problem, Rev. Mod. Phys.47, 773 (1975)

    work page 1975

  8. [8]

    Bulla, T

    R. Bulla, T. A. Costi, and T. Pruschke, Numerical renor- malization group method for quantum impurity systems, Rev. Mod. Phys.80, 395 (2008)

    work page 2008

  9. [9]

    Caffarel and W

    M. Caffarel and W. Krauth, Exact diagonalization ap- proach to correlated fermions in infinite dimensions: Mott transition and superconductivity, Phys. Rev. Lett. 72, 1545 (1994)

    work page 1994

  10. [10]

    Keiter and J

    H. Keiter and J. C. Kimball, Perturbation technique for the anderson hamiltonian, Phys. Rev. Lett.25, 672 (1970)

    work page 1970

  11. [11]

    Eckstein and P

    M. Eckstein and P. Werner, Nonequilibrium dynamical mean-field calculations based on the noncrossing approx- imation and its generalizations, Phys. Rev. B82, 115115 (2010)

    work page 2010

  12. [12]

    Dorda, M

    A. Dorda, M. Nuss, W. von der Linden, and E. Arrigoni, Auxiliary master equation approach to non – equilibrium correlated impurities, Phys. Rev. B89, 165105 (2014). 10

    work page 2014

  13. [13]

    Werner, J

    D. Werner, J. Lotze, and E. Arrigoni, Configuration interaction based nonequilibrium steady state impurity solver, Phys. Rev. B107, 075119 (2023)

    work page 2023

  14. [14]

    T. M. Mazzocchi, P. Gazzaneo, J. Lotze, and E. Arrigoni, Correlated mott insulators in strong electric fields: Role of phonons in heat dissipation, Phys. Rev. B106, 125123 (2022)

    work page 2022

  15. [15]

    Gazzaneo, T

    P. Gazzaneo, T. M. Mazzocchi, J. Lotze, and E. Arrigoni, Impact ionization processes in a photodriven mott insu- lator: Influence of phononic dissipation, Phys. Rev. B 106, 195140 (2022)

    work page 2022

  16. [16]

    T. M. Mazzocchi, D. Werner, P. Gazzaneo, and E. Ar- rigoni, Correlated mott insulators in a strong electric field: The effects of phonon renormalization, Phys. Rev. B107, 155103 (2023)

    work page 2023

  17. [17]

    T. M. Mazzocchi, D. Werner, M. Aichhorn, and E. Ar- rigoni, Mixed-configuration approximation for multior- bital systems out of equilibrium, Phys. Rev. B112, 155127 (2025)

    work page 2025

  18. [18]

    Cohen, E

    G. Cohen, E. Gull, D. R. Reichman, and A. J. Millis, Taming the dynamical sign problem in real-time evolu- tion of quantum many-body problems, Phys. Rev. Lett. 115, 266802 (2015)

    work page 2015

  19. [19]

    A. E. Antipov, Q. Dong, J. Kleinhenz, G. Cohen, and E. Gull, Currents and green’s functions of impurities out of equilibrium: Results from inchworm quantum monte carlo, Phys. Rev. B95, 085144 (2017)

    work page 2017

  20. [20]

    Erpenbeck, E

    A. Erpenbeck, E. Gull, and G. Cohen, Quantum monte carlo method in the steady state, Phys. Rev. Lett.130, 186301 (2023)

    work page 2023

  21. [21]

    Erpenbeck, T

    A. Erpenbeck, T. Blommel, L. Zhang, W.-T. Lin, G. Co- hen, and E. Gull, Steady-state properties of multi-orbital systems using quantum monte carlo, The Journal of Chemical Physics161, 094104 (2024)

    work page 2024

  22. [22]

    F. B. Anders, Steady-state currents through nanodevices: A scattering-states numerical renormalization-group ap- proach to open quantum systems, Phys. Rev. Lett.101, 066804 (2008)

    work page 2008

  23. [23]

    J. E. Han, Nonequilibrium statistics of biased kondo res- onance (2025), arXiv:2503.14400

    work page arXiv 2025

  24. [24]

    N´ u˜ nez Fern´ andez, M

    Y. N´ u˜ nez Fern´ andez, M. Jeannin, P. T. Dumitrescu, T. Kloss, J. Kaye, O. Parcollet, and X. Waintal, Learning feynman diagrams with tensor trains, Phys. Rev. X12, 041018 (2022)

    work page 2022

  25. [25]

    A. J. Kim and P. Werner, Strong coupling impurity solver based on quantics tensor cross interpolation, Phys. Rev. B111, 125120 (2025)

    work page 2025

  26. [26]

    Georges and W

    A. Georges and W. Krauth, Physical properties of the half-filled hubbard model in infinite dimensions, Phys. Rev. B48, 7167 (1993)

    work page 1993

  27. [27]

    Kajueter and G

    H. Kajueter and G. Kotliar, New iterative perturbation scheme for lattice models with arbitrary filling, Phys. Rev. Lett.77, 131 (1996)

    work page 1996

  28. [28]

    Lee and K

    W.-R. Lee and K. Park, Dielectric breakdown via emer- gent nonequilibrium steady states of the electric-field- driven mott insulator, Phys. Rev. B89, 205126 (2014)

    work page 2014

  29. [29]

    Martin-Rodero, F

    A. Martin-Rodero, F. Flores, M. Baldo, and R. Pucci, A new solution to the anderson-newns hamiltonian of chemisorption, Solid State Communications44, 911 (1982)

    work page 1982

  30. [30]

    Mart´ ın-Rodero, E

    A. Mart´ ın-Rodero, E. Louis, F. Flores, and C. Tejedor, Interpolative solution for the periodic anderson model of mixed-valence compounds, Phys. Rev. B33, 1814 (1986)

    work page 1986

  31. [31]

    Potthoff, T

    M. Potthoff, T. Wegner, and W. Nolting, Interpolating self-energy of the infinite-dimensional hubbard model: Modifying the iterative perturbation theory, Phys. Rev. B55, 16132 (1997)

    work page 1997

  32. [32]

    L.-F. m. c. Arsenault, P. S´ emon, and A.-M. S. Tremblay, Benchmark of a modified iterated perturbation theory approach on the fcc lattice at strong coupling, Phys. Rev. B86, 085133 (2012)

    work page 2012

  33. [33]

    S., A multi-orbital iterated perturbation theory for model hamiltonians and real material-specific calculations of correlated systems, Eur

    Dasari, Nagamalleswararao, Mondal, Wasim Raja, Zhang, Peng, Moreno, Juana, Jarrell, Mark, and Vid- hyadhiraja, N. S., A multi-orbital iterated perturbation theory for model hamiltonians and real material-specific calculations of correlated systems, Eur. Phys. J. B89, 202 (2016)

    work page 2016

  34. [34]

    E. G. C. P. van Loon, Two-particle correlations and the metal-insulator transition: Iterated perturbation theory revisited, Phys. Rev. B105, 245104 (2022)

    work page 2022

  35. [35]

    Tsuji and P

    N. Tsuji and P. Werner, Nonequilibrium dynamical mean-field theory based on weak-coupling perturbation expansions: Application to dynamical symmetry break- ing in the hubbard model, Phys. Rev. B88, 165115 (2013)

    work page 2013

  36. [36]

    Schwinger, Brownian motion of a quantum oscillator, J

    J. Schwinger, Brownian motion of a quantum oscillator, J. Math. Phys.2, 407 (1961)

    work page 1961

  37. [37]

    L. V. Keldysh, Diagram technique for nonequilibrium processes, Sov. Phys. JETP20, 1018 (1965)

    work page 1965

  38. [38]

    Haug and A.-P

    H. Haug and A.-P. Jauho,Quantum Kinetics in Trans- port and Optics of Semiconductors(Springer, Heidelberg, 1998)

    work page 1998

  39. [39]

    Throughout this work the onsite energy is chosen to be spin independent

    work page

  40. [40]

    Potthoff, T

    M. Potthoff, T. Herrmann, T. Wegner, and W. Nolting, The moment sum rule and its consequences for ferro- magnetism in the hubbard model, physica status solidi (b)210, 199 (1998)

    work page 1998

  41. [41]

    T. M. Mazzocchi and E. Arrigoni, Iterated perturbation theory for mott insulators in a static electric field with optical phonons, physica status solidi (b)n/a, 2300486

    work page

  42. [42]

    Nonetheless, benchmark comparisons with the AMEA impurity solver indicate that the IPT-n 0 approximation performs satis- factorily in all cases considered so far

    We note that, out of equilibrium, the Friedel sum rule discussed in previous work [26, 32] does not apply, and its proper generalization remains unclear. Nonetheless, benchmark comparisons with the AMEA impurity solver indicate that the IPT-n 0 approximation performs satis- factorily in all cases considered so far

    work page

  43. [43]

    For details on the AMEA impurity solver, we refer the reader to our previous work [11, 12]

    work page

  44. [44]

    Bulla, A

    R. Bulla, A. C. Hewson, and T. Pruschke, Numerical renormalization group calculations for the self-energy of the impurity anderson model, Journal of Physics: Con- densed Matter10, 8365 (1998)

    work page 1998

  45. [45]

    Meir and N

    Y. Meir and N. S. Wingreen, Landauer formula for the current through an interacting electron region, Phys. Rev. Lett.68, 2512 (1992)

    work page 1992

  46. [46]

    However, for our choice of leads, neither of these param- eters enters the retarded (advanced) reservoir GF

    We point out that the retarded and Keldysh components of the impurity GF as well as the Keldysh component of the lead GF retain a non-trivial dependence on Φ andT. However, for our choice of leads, neither of these param- eters enters the retarded (advanced) reservoir GF

    work page

  47. [47]

    By construction, the generalized distribution functionF reduces to the Fermi-Dirac distributionf FD in the equi- librium case

    work page