REVIEW 3 major objections 4 minor 37 references
Improved treatment of fermion-boson vertices and Bethe-Salpeter equations in non-local extensions of dynamical mean field theory
T0 review · 3 major / 4 minor · reviewed 2026-08-14 · deepseek-v4-flash
Pith's one-line read This paper derives analytic formulas that remove most of the finite-frequency-box error in non-local dynamical mean-field theory, improving the truncation error from O(1/ν_max) to O(1/ν_max^4).
desk verdict A genuinely useful, honestly-scoped methodological advance: analytic out-of-box corrections for non-local DMFT vertices, with the advertised O(1/ν^4) scaling resting on an unproven but plausible decay assumption. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The analytic tail correction $X_q$ and the vertex correction factor $Z_{\nu q}$ are the objects that carry the argument. $X_q$ is the sum over frequencies outside the box of the bare bubble plus one insertion of the asymptotic vertex tail; $Z_{\nu q}$ measures how much the three-leg vertex is renormalized by the tail outside the box. Both are computed from the asymptotic form (3)-(4), which uses only local charge, spin, and particle-particle susceptibilities, so they are cheap to evaluate for arbitrarily large frequencies. Eq. (25) is the load-bearing identity: $\Phi_{\nu\nu'q} = \Phi^{\mathrm{box}}_{\nu\nu'q} + U_q - Z_{\nu q}\tilde U_q Z_{\nu'q}$, with $\tilde U_q = U_q/(1-U_q X_q)$. This identity removes the $1/\nu_{\max}$ error and leaves the smaller $1/\nu_{\max}^4$ error.
What would settle it
Take a model in which the asymptotic form is known to fail, such as a non-ladder vertex with retarded non-local interactions, evaluate the full two-particle irreducible vertex on a very large frequency box, and compare the left and right sides of Eq. (25) as the box size grows. If the difference does not shrink like $1/\nu_{\max}^4$ — or if the correction pushes the vertex away from the exact large-box value — the assumed tail is the culprit.
Extended reading notes
Core claim
The central claim is that the out-of-box frequency contributions, which previously required numerical treatment of the asymptotic tail in a larger box, can be evaluated analytically to leading order. Starting from the asymptotic form (3) of the two-particle irreducible vertex, the paper defines $X_q$ and $Z_{\nu q}$ and derives Eqs. (8), (11), (14), (18), and (25) for the full fermion-boson vertex, the non-local susceptibility, the reduced fermion-boson vertex, the irreducible susceptibility, and the two-particle irreducible vertex, respectively. The relation $\Phi_{\nu\nu'q} = \Phi^{\mathrm{box}}_{\nu\nu'q} + U_q - Z_{\nu q}\tilde U_q Z_{\nu'q}$ is the key identity: it converts a matrix inversion performed only inside the box into the physical vertex. In the companion numerical test on the two-dimensional Hubbard model, the corrected vertices change very little with box size, while the uncorrected ones extrapolate slowly; the corrected error is consistent with $1/\nu_{\max}^4$ scaling.
Load-bearing premise
The whole scheme depends on the high-frequency asymptotic form (3) of the two-particle irreducible vertex, together with the assumption that deviations from this form decay at least as $1/\max(|\nu|,|\nu'|)^3$; if that tail is wrong or decays too slowly, the size of the correction terms $X_q$ and $Z_{\nu q}$ would be misestimated.
Editorial extensions
If this is right
- A single calculation with a small frequency box should suffice for ladder-type non-local DMFT calculations, instead of running several box sizes and extrapolating.
- The method gives a direct relation between the local and non-local two-particle irreducible vertices, making it easier to feed local input into non-local Bethe-Salpeter equations.
- The same formulas apply to both charge and spin channels, and they can be adapted to cases where the vertex depends on more than one momentum transfer.
- The claimed accuracy means the main remaining error is set by the physical decay of the irreducible vertex toward its asymptotic form, not by the truncation.
Reading between the lines
- If the same asymptotic form holds in the particle-particle channel, an analogous set of box corrections could be derived for pairing susceptibilities, where the tail is built from the same local susceptibilities.
- The quartic improvement is conditional on the decay of the deviation from Eq. (3); in models with long-range or retarded interactions the effective exponent would be lower, and the size of $X_q$ and $Z_{\nu q}$ would reveal that.
- The construction suggests a general recipe for any truncated Matsubara summation: add back the analytically known tail before inverting, rather than enlarging the box.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript proposes an analytic treatment of fermionic frequency tails in Bethe-Salpeter equations and fermion-boson vertices for non-local extensions of DMFT. The key idea is to split all frequency summations into a small box, where vertices are known numerically, and the outside region, where the asymptotic form of the irreducible vertex, Eq. (3), is used. The outside contribution is expressed in closed form through the quantities X_q and Z_q, yielding explicit formulas for the full and reduced fermion-boson vertices (Eqs. (8), (14), (21)), the full and irreducible susceptibilities (Eqs. (11), (18)), and the two-particle irreducible vertex (Eq. (25)). The method is tested on spin-sector DMFT results for a two-dimensional Hubbard model, with reported agreement with CT-QMC local susceptibilities and a much weaker dependence on frequency-box size than the uncorrected calculation.
Significance. If the central claims hold, the work is a useful contribution: the correction formulas are explicit, parameter-free, and reduce the need for large frequency boxes and for separate numerical treatment of the large box in earlier schemes. The algebraic nature of the derivation and the concrete numerical demonstration of improved convergence with box size are strengths. However, the headline O(1/ν_max^4) error-scaling claim rests on an unproven decay assumption for the residual irreducible vertex, and the numerical support is limited to one channel and one model parameter set, which weakens the generality asserted in the abstract and conclusion.
major comments (3)
- [Section IV, discussion following Eq. (25)] The central claim that the method achieves O(1/ν_max^4) accuracy is not established. The argument assumes that the deviation of the true irreducible vertex from the asymptotic form (3) decays as 1/max(|ν|,|ν'|)^3; this exponent is not derived, and the text describes it only as 'expected.' Because X_q and Z_q are built entirely from the U_q + \barΦ tail, any component of Φ outside that tail with a slower power law will enter the out-of-box sums at lower order. The numerical fits in Sec. V (Fig. 3) use an a + b/ν_max^4 + c/ν_max^5 form, which presupposes the claimed scaling and therefore cannot by itself confirm the exponent. To support the conclusion, the author should either derive (or cite a derivation of) the 1/ν^3 decay of the residual, or directly measure this residual as a function of frequency and demonstrate the absence of 1/ν^2 or 1/ν^3 contamination.
- [Section IV, transition from Eq. (24) to Eq. (25)] The key relation (25) is obtained by 'algebraic transformations' that are not shown, and Eq. (21) similarly refers to 'algebraic manipulations' of Appendix C of Ref. [30]. Since Eq. (25) is the central new result and is used to convert the box-restricted vertex Φ^{box} into the physical Φ, the derivation should be reproduced in detail, for example in an appendix. Without this, the reader cannot verify the consistency of signs, the treatment of the out-of-box sums, or the order at which terms are neglected. The numerical checks in Sec. V are too limited to substitute for this derivation.
- [Section V] The numerical validation covers only the spin channel of a single two-dimensional Hubbard model with U=10t, t'=0.15t, n=0.96, T=0.08t (and one higher-temperature half-filled run). The paper advertises the method for a broad class of non-local DMFT extensions (ladder DΓA, dual fermion/boson, TRILEX, DMF2RG), but the asymptotic form (3) and the assumed 1/ν^3 decay of its correction are model- and approximation-dependent. In addition, the comparison of the local susceptibility from Eq. (16) with the CT-QMC result is partly a consistency check, because the tail functions X_q and Z_q themselves use local susceptibilities obtained from the same CT-QMC solver. At minimum, the authors should state these limitations explicitly in the conclusions and provide evidence in at least one further channel, parameter set, or non-local context.
minor comments (4)
- [Fig. 3] The axis label 'tX0, Z' in panel (b) is unclear; please spell out the plotted quantities (e.g., tX_0^{(1)}, tX_0^{(2)}, Z_{ν,0}-1) and ensure that all axis labels render correctly.
- [Sec. III, after Eq. (8)] The statement that the second term in X_q and the difference Z_q-1 are 'verified numerically' should be qualified: the verification is performed only for the spin channel of the Hubbard model, not for the full range of claimed applications.
- [Sec. V, Fig. 3] The extrapolations shown by dotted lines and the fits with a+b/ν_max^4+c/ν_max^5 are presented without error bars or fit parameters; reporting these numbers would make the claimed scaling easier to assess.
- [Sec. V] The sentence 'we find the results of extrapolation consistent with those for vertices, obtained without account of finite frequency box effects' is not quantified; please provide the extrapolated vertex values and the corresponding fit parameters.
Circularity Check
No significant circularity: the derivation is algebraic and parameter-free, built on an externally sourced vertex asymptotics; self-citations are minor and non-load-bearing.
full rationale
The central derivation chain starts from the assumed asymptotic form of the two-particle irreducible vertex, Eq. (3), and the corresponding asymptotics of the reducible vertex, Eq. (5), both attributed to prior work by Tagliavini et al. (Ref. [33]), not to the present author. The correction functions X_q and Z_q are then obtained as explicit sums over the outside-box region using that asymptotic tail, with no fitted parameters and no quantity being renamed as a prediction. Equations (8), (11), (14), (18), (21), and (25) are algebraic rearrangements of the Bethe-Salpeter equation under the stated asymptotics; none of them reduces by construction to its own input. The self-citation to Ref. [30] is used only for 'algebraic manipulations, similar to those described in Appendix C' and does not carry the load-bearing argument. The numerical section does contain a consistency-check element: the local susceptibilities entering the asymptotic tail in Eq. (4) are also compared with the CT-QMC result in Sec. V, so that particular comparison is partly a self-consistency test rather than an independent benchmark for those susceptibilities. However, this overlap does not make the derivation circular, because the derived vertex correction formulas do not fit those susceptibilities to the final vertices. The main caveat is instead a correctness risk, not circularity: the claimed O(1/ν_max^4) error improvement relies on the statement that deviations from Eq. (3) 'is expected to scale as 1/max(|ν|, |ν'|)^3', which is asserted as an expectation rather than derived, and the numerical verification covers only one model and one parameter set. That is a limitation of the evidence, but not an equivalence between the output and the input by construction.
Assumptions & free parameters
assumptions (5)
- domain assumption Two-particle irreducible vertex Φ_{νν'q} has the asymptotic form U_q^{c(s)} + Φ̄_{νν'ω} for large ν or ν' (Eq. 3).
- domain assumption The tail Φ̄_{νν'ω} is local and expressed through local charge, spin, particle-particle susceptibilities and local interaction v_c(ω), decaying as (|ν|-|ν'|)^{-2} for large frequency differences.
- standard math The bubble χ0_{νq} decays as 1/ν^2 for large fermionic frequency.
- domain assumption The deviation of the actual vertex from the asymptotic form (3) is O(1/max(|ν|,|ν'|)^3), so neglected terms contribute O(1/ν_max^4).
- domain assumption In the ladder approximation the irreducible vertex Φ is local and identical for the local and non-local problems.
Cite this review
Pith. "Pith review of Improved treatment of fermion-boson vertices and Bethe-Salpeter equations in non-local extensions of dynamical mean field theory." pith.science (2026). https://pith.science/paper/RFYEWLBK
@misc{pith2026190803198,
author = {Pith},
title = {Pith review of: Improved treatment of fermion-boson vertices and Bethe-Salpeter equations in non-local extensions of dynamical mean field theory},
year = {2026},
howpublished = {\url{https://pith.science/paper/RFYEWLBK}},
note = {Machine review of arXiv:1908.03198}
}
read the original abstract
We reconsider the procedure of calculation of fermion-boson vertices and numerical solution of Bethe-Salpeter equations, used in non-local extensions of dynamical mean-field theory. Because of the frequency dependence of vertices, finite frequency box for matrix inversions is typically used, which often requires some treatment of asymptotic behaviour of vertices. Recently [Phys. Rev. B 83, 085102 (2011); 97, 235140 (2018)] it was proposed to split the considered frequency box into smaller and larger one; in the smaller frequency box the numerically exact vertices are used, while beyond this box asymptotics of vertices are applied. Yet, this method requires numerical treatment of vertex asymptotics (including corresponding matrix manipulations) in the larger frequency box and/or knowing fermion-boson vertices, which may be not convenient for numerical calculations. In the present paper we derive the formulae which treat analytically contribution of vertices beyond chosen frequency box, such that only numerical operations with vertices in the chosen small frequency box are required. The method is tested on the Hubbard model and can be used in a broad range of applications of non-local extensions of dynamical mean-field theory.
Figures
Reference graph
Works this paper leans on
-
[33]
A. Tagliavini, S. Hummel, N. Wentzell, S. Andergassen, A. Toschi, and G. Rohringer, Phys. Rev. B 97, 235140 (2018)
work page 2018
-
[30]
A. A. Katanin, Phys. Rev. B 99, 115112 (2019)
work page 2019
-
[1]
Fradkin, S
E. Fradkin, S. A. Kivelson, and J. M. Tranquada, Rev. Mod. Phys. 87, 457 (2015)
2015
-
[2]
Dagotto, Rev
E. Dagotto, Rev. Mod. Phys. 66, 763 (1994)
1994
-
[3]
W. E. Pickett, Rev. Mod. Phys. 61, 433 (1989)
work page 1989
-
[4]
K. Haule and G. Kotliar, New J. Phys. 11, 025021 (2009); A. Georges, L. de’ Medici, and J. Mravlje, Ann. Rev. Cond. Matt. Physics 4, 137 (2013); L. Fanfarillo and E. Bascones, Phys. Rev. B 92, 075136 (2015); A. S. Belozerov, A. A. Katanin, and V. I. Anis imov, Phys. Rev. B 97, 115141 (2018)
work page 2009
-
[5]
Metzner and D
W. Metzner and D. Vollhardt, Phys. Rev. Lett. 62, 324 (1989)
1989
-
[6]
A. Georges, G. Kotliar, W. Krauth, and M. Rozenberg, Rev. Mod. Phys. 68, 13 (1996); G. Kotliar and D. Vollhardt, Physics Today 57, 53 (2004)
work page 1996
Show all 37 references
-
[7]
Chitra and G
R. Chitra and G. Kotliar, Phys. Rev. Lett. 84, 3678 (2000); P. Sun and G. Kotliar, Phys. Rev. B 66, 085120 (2002)
2000
-
[8]
Si and J
Q. Si and J. L. Smith, Phys. Rev. Lett. 77, 3391 (1996); J. L. Smith and Q. Si, Phys. Rev. B 61, 5184 (2000)
1996
-
[9]
Maier, M
T. Maier, M. Jarrell, T. Pruschke, and M. H. Hettler, Rev. Mod. Phys. 77, 1027 (2005)
2005
-
[10]
Rohringer, H
G. Rohringer, H. Hafermann, A. Toschi, A. A. Katanin, A. E. Antipov, M. I. Katsnelson, A. I. Lichtenstein, A. N. Rubtsov, K. Held, Rev. Mod. Phys. 90, 025003 (2018)
2018
-
[11]
Toschi, A
A. Toschi, A. Katanin, and K. Held, Phys. Rev. B 75, 045118 (2007)
2007
-
[12]
K. Held, A. A. Katanin, A. Toschi, Prog. Theor. Phys. Sup pl. 176, 117 (2008)
2008
-
[13]
A. A. Katanin, A. Toschi, K. Held, Phys. Rev. B 80, 075104 (2009)
2009
-
[14]
Rohringer, A
G. Rohringer, A. Toschi, A. A. Katanin, K. Held, Phys. Re v. Lett. 107, 256402 (2011)
2011
-
[15]
Slezak, M
C. Slezak, M. Jarrell, Th. Maier, and J. Deisz, cond-mat /0603421 (unpublished); J. Phys.: Condens. Matter 21, 435604 (2009)
2009
-
[16]
Toschi, G
A. Toschi, G. Rohringer, A. A. Katanin, and K. Held, Ann. der Physik, 523, 698 (2011); A. Galler, P. Thunstr¨ om, P. Gunacker, J. M. Tomczak, and K. Hel d, Phys. Rev. B 95, 115107 (2017)
2011
-
[17]
Valli, T
A. Valli, T. Sch¨ afer, P. Thunstr¨ om, G. Rohringer, S. A ndergassen, G. Sangiovanni, K. Held, and A. Toschi, Phys. Rev. B 91, 115115 (2015); A. Kauch, F. H¨ o rbinger, G. Li, and K. Held, 14 arXiv:1901.09743 (unpublished); P. Pudleiner, P. Thunstr ¨ om, A. Valli, A. Kauch, G...
2015 arXiv
-
[18]
A. N. Rubtsov, M. I. Katsnelson, A. I. Lichtenstein, con d-mat/0612196 (unpublished); Phys. Rev. B 77, 033101 (2008)
2008
-
[19]
Hafermann, S
H. Hafermann, S. Brener, A. N. Rubtsov, M. I. Katsnelson , A. I. Lichtenstein, JETP Lett. 86, 677 (2007)
2007
-
[20]
Brener, H
S. Brener, H. Hafermann, A. N. Rubtsov, M. I. Katsnelson , A. I. Lichtenstein, Phys. Rev. B 77, 195105 (2008)
2008
-
[21]
A. N. Rubtsov, M. I. Katsnelson, A. I. Lichtenstein, A. G eorges, Phys.Rev. B 79 045133 (2009)
2009
-
[22]
Otsuki, H
J. Otsuki, H. Hafermann, and A. I. Lichtenstein, Phys. R ev. B 90, 235132 (2014)
2014
-
[23]
A. N. Rubtsov, M. I. Katsnelson, A. I. Lichtenstein, Ann . of Phys. 327, 1320 (2012)
2012
-
[24]
E. G. C. P. van Loon, A. I. Lichtenstein, M. I. Katsnelson , O. Parcollet, and H. Hafermann, Phys. Rev. B 90, 235135 (2014)
2014
-
[25]
E. A. Stepanov, E. G. C. P. van Loon, A. A. Katanin, A. I. Li chtenstein, M. I. Katsnelson, and A. N. Rubtsov, Phys. Rev. B 93, 045107 (2016)
2016
-
[26]
E. A. Stepanov, A. Huber, E. G. C. P. van Loon, A. I. Lichte nstein, and M. I. Katsnelson, Phys. Rev. B 94, 205110 (2016)
2016
-
[27]
Ayral and O
T. Ayral and O. Parcollet, Phys. Rev. B 92, 115109 (2015); 93, 235124 (2016)
2015
-
[28]
Taranto, S
C. Taranto, S. Andergassen, J. Bauer, K. Held, A. Katani n, W. Metzner, G. Rohringer, A. Toschi, Phys. Rev. Lett. 112, 196402 (2014)
2014
-
[29]
Vilardi, C
D. Vilardi, C. Taranto, and W. Metzner, Phys. Rev. B 99, 104501 (2019)
2019
-
[31]
E. G. C. P. van Loon, F. Krien, H. Hafermann, A. I. Lichten stein, M. I. Katsnelson, Phys. Rev. B 98, 205148 (2018)
2018
-
[32]
Jan Kuneˇ s, Phys. Rev. B 83, 085102 (2011)
2011
-
[34]
J. A. Hertz and D. M. Edwards, J. Phys. F: Met. Phys. 3, 2174 (1973)
1973
-
[35]
Krien, Phys
F. Krien, Phys. Rev. B 99, 235106 (2019)
2019
-
[36]
Wang, Zi Yang Meng, L
Li Huang, Y. Wang, Zi Yang Meng, L. Du, P. Werner, and Xi Da i, Comp. Phys. Comm. 195, 15 140 (2015)
2015
-
[37]
Li Huang, Comp. Phys. Comm. 221, 423 (2017). 16
2017
Reviewed August 14, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.