Recognition: 2 theorem links
· Lean TheoremSelf-Consistent Random Phase Approximation from Projective Truncation Approximation Formalism
Pith reviewed 2026-05-17 03:12 UTC · model grok-4.3
The pith
Projective truncation of Green's function equations produces a self-consistent random phase approximation valid at any temperature.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
Applying the projective truncation approximation to the equation of motion of two-time Green's functions closes the hierarchy and produces the self-consistent random phase approximation. The resulting sc-RPA holds for arbitrary temperatures and recovers Rowe's formalism at zero temperature; the same PTA framework also rationalizes the original formula and opens a path to controlled extensions of sc-RPA.
What carries the argument
The projective truncation approximation (PTA) for the equation of motion of two-time Green's functions, which projects the infinite hierarchy onto a finite subspace to obtain closed equations that become the self-consistent RPA.
If this is right
- The sc-RPA equations can now be used at finite temperature without additional ad-hoc approximations.
- Setting temperature to zero in the derived equations recovers Rowe's original formalism exactly.
- The PTA supplies a general framework for extending sc-RPA to include higher-order correlations or different truncation subspaces.
- The 1D spinless-fermion implementation demonstrates that Luttinger-liquid features and bound states appear naturally in the spectral function once N-representability is enforced.
Where Pith is reading between the lines
- The same PTA route could be applied to derive temperature-dependent extensions of other diagrammatic approximations such as the Bethe-Salpeter equation.
- Difficulties of RPA for symmetric states noted in the paper may be traceable to the choice of truncation subspace and could be tested by enlarging that subspace.
- The static-component problem of PTA might be alleviated by retaining a small number of additional projection operators, providing a concrete numerical test for future implementations.
Load-bearing premise
The projective truncation approximation supplies a sufficiently accurate closure to the infinite Green's-function hierarchy, and the N-representability constraints can be imposed without introducing uncontrolled errors in the chosen model.
What would settle it
Numerical mismatch between the sc-RPA ground-state energy or correlation functions for the one-dimensional spinless-fermion chain and independent benchmark values obtained by exact diagonalization or density-matrix renormalization group.
Figures
read the original abstract
We derive the self-consistent random phase approximations (sc-RPA) from the projective truncation approximation (PTA) for the equation of motion of two-time Green's function. The obtained sc-RPA applies to arbitrary temperature and recovers the Rowe's formalism at zero temperature. The PTA formalism not only rationalize Rowe's formula, but also provides a general framework to extend sc-RPA. We implement the sc-RPA calculation for the one-dimensional spinless fermion model in the parameter regime of disordered ground state, with the N-representability constraints enforced. The obtained ground state energy, correlation function, and density spectral function agree well with existing results. The features of the Luttinger liquid ground state and the continuum/bound state in the spectral function are well captured. We discuss several issues concerning the approximations made in RPAs, difficulties of RPA for symmetric state, and the static component problem of PTA.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript derives the self-consistent random phase approximation (sc-RPA) from the projective truncation approximation (PTA) applied to the equation of motion for two-time Green's functions. It claims that the resulting sc-RPA is valid at arbitrary temperatures and recovers Rowe's zero-temperature formalism. The PTA is presented as both a rationalization of Rowe's approach and a general framework for extending sc-RPA. Numerical results are shown for the one-dimensional spinless fermion model in the disordered ground-state regime with N-representability constraints enforced; ground-state energies, correlation functions, and density spectral functions are reported to agree with existing benchmarks while capturing Luttinger-liquid features and continuum/bound-state spectral structure. The discussion addresses limitations of RPA approximations, difficulties for symmetric states, and the static component problem of PTA.
Significance. If the finite-temperature extension can be established, the work would supply a systematic route to thermal sc-RPA calculations within a Green's-function hierarchy framework, potentially useful for correlated systems beyond the zero-temperature limit. The enforcement of N-representability and the reported numerical agreement in the 1D model constitute concrete strengths. At present the significance remains provisional because only ground-state results are provided despite the arbitrary-temperature claim.
major comments (2)
- [Abstract and Discussion] Abstract and Discussion section: The central claim that the derived sc-RPA applies to arbitrary temperature is not supported by the presented evidence. All numerical implementations and benchmarks are restricted to the ground state (T=0) of the 1D spinless-fermion chain in the disordered regime. The manuscript explicitly flags the 'static component problem of PTA' as an open issue, which directly affects the validity of thermal averages and the self-consistency loop at finite T.
- [Implementation] Implementation section: While N-representability constraints are enforced, no explicit description is given of how the self-consistency loop is closed at finite temperature, how sum rules are preserved under the PTA truncation, or what error controls are applied when the static component is present.
minor comments (3)
- The abstract states agreement with 'existing results' but does not specify the exact reference calculations, system sizes, or parameter values used for the comparisons.
- [Derivation] Clarify the precise form of the projective truncation ansatz employed for the two-particle Green's function and how it reduces to Rowe's equations at T=0.
- [Discussion] The discussion of 'difficulties of RPA for symmetric state' is mentioned but not illustrated with any concrete example or diagnostic from the 1D implementation.
Simulated Author's Rebuttal
We thank the referee for the careful reading and constructive comments on our manuscript. We address the major comments point by point below, indicating where revisions will be made.
read point-by-point responses
-
Referee: [Abstract and Discussion] Abstract and Discussion section: The central claim that the derived sc-RPA applies to arbitrary temperature is not supported by the presented evidence. All numerical implementations and benchmarks are restricted to the ground state (T=0) of the 1D spinless-fermion chain in the disordered regime. The manuscript explicitly flags the 'static component problem of PTA' as an open issue, which directly affects the validity of thermal averages and the self-consistency loop at finite T.
Authors: The derivation of sc-RPA from PTA is performed within the finite-temperature two-time Green's function equation-of-motion framework and is shown to recover Rowe's zero-temperature formalism as a special case. We agree that all numerical results and benchmarks are restricted to the ground state (T=0) and that the static component problem of PTA is explicitly identified in the manuscript as an unresolved issue affecting thermal averages. This does limit the strength of the arbitrary-temperature claim in the absence of finite-T numerical demonstrations. In the revised manuscript we will update the abstract and discussion section to state more clearly that the formalism is derived for arbitrary temperature while noting that the current numerical implementation is at T=0 due to the open static component issue. revision: yes
-
Referee: [Implementation] Implementation section: While N-representability constraints are enforced, no explicit description is given of how the self-consistency loop is closed at finite temperature, how sum rules are preserved under the PTA truncation, or what error controls are applied when the static component is present.
Authors: We will expand the Implementation section to provide a more explicit description of the self-consistency loop as realized for the T=0 calculations, including the precise manner in which N-representability constraints are imposed at each iteration. Sum-rule preservation follows directly from the projective character of the PTA truncation applied to the equation-of-motion hierarchy; we will add a concise explanation of this property. Error controls for the static component are discussed in the context of the ground-state results; we will elaborate on these controls while noting that a complete finite-temperature procedure cannot be specified until the static component problem is resolved. revision: partial
- Complete algorithmic details for closing the self-consistency loop at finite temperature, which cannot be provided until the static component problem of PTA is addressed.
Circularity Check
Derivation of sc-RPA from PTA is self-contained with independent truncation closure; no reduction to fitted inputs or self-citation chains.
full rationale
The manuscript derives sc-RPA by applying the projective truncation approximation directly to the two-time Green's function equation-of-motion hierarchy, yielding a closed set of equations that recover Rowe's zero-temperature form as a special case. PTA supplies the truncation ansatz as an external closure rule rather than being defined in terms of the resulting RPA; the subsequent N-representability enforcement and numerical implementation on the 1D spinless-fermion chain are compared against independent benchmarks (ground-state energy, correlations, spectral functions). No equation is shown to equal its own input by construction, no parameter is fitted to a subset and then relabeled a prediction, and any self-citations serve only as background rather than load-bearing justification for the central mapping. The derivation therefore remains non-circular and externally falsifiable through the reported numerical tests.
Axiom & Free-Parameter Ledger
axioms (1)
- domain assumption The equation of motion for two-time Green's functions admits a projective truncation that closes the hierarchy while preserving key physical properties.
Lean theorems connected to this paper
-
IndisputableMonolith/Cost/FunctionalEquation.leanwashburn_uniqueness_aczel unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
We derive the self-consistent random phase approximations (sc-RPA) from the projective truncation approximation (PTA) for the equation of motion of two-time Green's function... basis {a†_α a_β} (α ≠ β)
-
IndisputableMonolith/Foundation/DimensionForcing.leanalexander_duality_circle_linking unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
The obtained ground state energy, correlation function, and density spectral function agree well with existing results... Luttinger liquid ground state
What do these tags mean?
- matches
- The paper's claim is directly supported by a theorem in the formal canon.
- supports
- The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
- extends
- The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
- uses
- The paper appears to rely on the theorem as machinery.
- contradicts
- The paper's claim conflicts with a theorem or certificate in the canon.
- unclear
- Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.
Reference graph
Works this paper leans on
-
[1]
(81) Here, Q† ν is the excitation operator introduced in Rowe’s work [2]
(80) The first equation of Eq.(35) then becomes O1ν = D/ 2∑ i=1 [ Xiν ( ˜Ai)d − Yiν ( ˜Bi)d ] ≈ D/ 2∑ i=1 [ Xiν ˜Ai − Yiν ˜Bi ] = Q† ν . (81) Here, Q† ν is the excitation operator introduced in Rowe’s work [2]. λ 1ν < 0 corresponds to raising energy when O1ν is applied to an eigenstate, being consistent with the assumption in Rowe’s formula that Q† ν |0⟩ =...
-
[2]
∆ ̸= 0 corresponds to a flux applied to the charge of fermions
0 is set as the energy unit. ∆ ̸= 0 corresponds to a flux applied to the charge of fermions. In this work, a tiny ∆ ̸= 0 is applied lift the degeneracy of single particle energy levels. Due to Jordan-Wigner transformation [52], at ∆ = 0, this model is equivalent to the anisotropic Heisenberg model, also known as XXZ model, whose properties are well studied...
-
[3]
This transformation maps the basis operator Aq k → ˜Aq k = A2π − q 2π − k
spatial inversion symmetry at ∆ = 0 For ∆ = 0, the Hamiltonian is invariant with respect to the spatial inversion transformation, ck → ˜ck = c2π − k. This transformation maps the basis operator Aq k → ˜Aq k = A2π − q 2π − k. Due to the invariance of H(∆ = 0) under this transformation, we have the following constraints at ∆ = 0, X q kk′ = X 2π − q 2π − k 2...
-
[4]
In momentum space, it is written as ck → ˜ck = − c† π − k
particle-hole symmetry at ∆ = 0 and N/L = 1/ 2 At ∆ = 0 and half fermion filling N/L = 1 / 2, Hamil- tonian Eq.(96) is invariant with respect to the particle- hole transformation, ci → ˜ci = (− 1)i+1c† i . In momentum space, it is written as ck → ˜ck = − c† π − k. This transfor- mation maps the basis operator Aq k → ˜Aq k = − Aq π − k− q. Due to the invari...
-
[5]
This transformation maps the basis operator Aq k → ˜Aq k = − A− q k+q+π
combined inversion and particle-hole symmetry at N/L = 1/ 2 and arbitrary ∆ For arbitrary ∆ and at half filling, although H(∆) is not invariant under the above two transformations sepa- rately, it is invariant under the combined transformation ck → ˜ck = − c† k+π . This transformation maps the basis operator Aq k → ˜Aq k = − A− q k+q+π . Due to this symmet...
-
[6]
The lines are for guiding eyes
(red circles). The lines are for guiding eyes. The resul ts of sc-RPA are obtained at parameters L = 192, N = 96, t = 1 . 0, T = 10 − 3, ∆ = 10 − 5, while those of bosonization are at L = 2 N = ∞ , t = 1 . 0, T = 0, and ∆ = 0. The long-distance asymptotic power from Eqs.(117) and (118) ar e marked in the figure. on Tk above or below EF . The fractional occ...
-
[7]
For the sc-RPA on basis B1, operators higher than the quadratic order are omitted
-
[8]
Operators a† α aα (1 ≤ α ≤ L) are omitted from the basis
-
[9]
This is exact for T = 0 and no ground state degeneracy
The static contribution appearing in the spectral theorem is neglected, i.e., ⟨(a† α aβ )† 0(a† α ′aβ ′)0⟩ ≈ 0 (α ̸= β and α ′ ̸= β ′). This is exact for T = 0 and no ground state degeneracy. Otherwise it is not exact
-
[10]
For PTA on the restricted B2 basis, higher order terms of the amplitude Y are neglected in the self- consistent calculation of one- and two-particle den- sities
-
[11]
In most practical RPA calculations, the working orbitals are pre-selected orbitals and are approxi- mately treated as natural orbitals
-
[12]
In the renormalized RPA, Hartree-Fock-like decou- pling of two-particle densities are introduced
-
[13]
The single- and two-particle densities are evaluated on the Hartree-Fock ground state
In the standard RPA, the Fock contribution in L is neglected. The single- and two-particle densities are evaluated on the Hartree-Fock ground state. B. The Case of Singular I Matrix At very high temperatures or for systems with degener- ate orbitals (such as the materials with flat energy band), there may be orbitals with very close occupation num- bers, n...
-
[14]
generalized eigen problem First, we summarize the algorithm of solving the gen- eralized eigen problem Eq.(116), i.e., LqUq = IqUqΛ q. (E1) We seek for the matrices Uq and the diagonal Λ q for a given pair of Hermitian matrices Lq ⪯ 0 and Iq. Since the same algorithm applies to every q, below, we drop the superscript q. Given that L is negative semi-defini...
-
[15]
The set of linear equa- tions to be solved for a fixed Mq are Eqs.(104), (106), and (108)
linear equations for C q and {⟨nk⟩} Second, we summarize the algorithm of solving the lin- ear equations for Cq and {⟨nk⟩}. The set of linear equa- tions to be solved for a fixed Mq are Eqs.(104), (106), and (108). That is, Cq = − Iq ( eβ Mq − 1 ) − 1 , ⟨nk⟩ = 1 L − N ∑ k′(̸=k) Ck− k′ k′k′ , ∑ k ⟨nk⟩ = N. (E9) We denote Wq = ( eβ Mq − 1 ) − 1 . Inserting t...
-
[16]
D. Bohm and D. Pines, Phys. Rev. 82, 625 (1951); ibid. 85, 338 (1952); ibid. 92, 609 (1953); D. Pines, Phys. Rev. 92, 626 (1953)
work page 1951
-
[17]
D. J. Rowe, Rev. Mod. Phys. 40, 153 (1968)
work page 1968
- [18]
- [19]
- [20]
-
[21]
G. Giuliani and G. Vignale, Quantum Theory of the Elec- tron Liquid , Cambridge University Press 2005
work page 2005
- [22]
-
[23]
J. Li, N. D. Drummond, P. Schuck, and V. Olevano, Sci- Post Phys. 6, 040 (2019)
work page 2019
- [24]
- [25]
-
[26]
X. Ren, P. Rinke, C. Joas, and M. Scheffler, J. Mater. Sci. 47 (2012)
work page 2012
-
[27]
G. P. Chen, V. K. Voora, M. M. Agee, S. G. Balasubra- mani, and F. Furche, Annu. Rev. Phys. Chem. 68, 421 (2017)
work page 2017
-
[28]
I. Y. Zhang, R. Shi, and X. Ren, Electron. Struct. 7 043002 (2025)
work page 2025
-
[29]
J. P. Perdew and K. Schmidt, Jacob’s ladder of density functional approximations for the exchange-correlation energy in Density Functional Theory and its Application to Materials, edited by V. Van Doren, C. Van Alsenoy, and P. Geerlings (AIP, Melville, 2001)
work page 2001
- [30]
-
[31]
D. C. Langreth and J. P. Perdew Phys. Rev. B 15, 2884 (1977)
work page 1977
- [32]
-
[33]
J. F. Dobson, in Topics in Condensed Matter Physics, Edited by M. P. Das, Ch. 7 (Nova: New York, 1994)
work page 1994
-
[34]
Y. Gao, W. Zhu, and X. Ren, Phys. Rev. B 101, 035113 (2020)
work page 2020
-
[35]
G. E. Scuseria, T. M. Henderson, and D. C. Sorensen, J. Chem. Phys. 129, 231101 (2008)
work page 2008
-
[36]
R. J. Bartlett and M. Musia/suppress l, Rev. Mod. Phys.79, 291 (2007)
work page 2007
- [37]
-
[38]
Zwanzig, in Lectures in Theoretical Physics (Inter- science, New York, 1961), Vol
R. Zwanzig, in Lectures in Theoretical Physics (Inter- science, New York, 1961), Vol. 3
work page 1961
- [39]
- [40]
-
[41]
Yu. A. Tserkovnikov, Theor. Math. Phys. 49, 993 (1981)
work page 1981
-
[42]
Yu. A. Tserkovnikov, Theor. Math. Phys. 118, 85 (1999)
work page 1999
- [43]
-
[44]
P. Fan, N. H. Tong, and Z. G. Zhu, Phys. Rev. B 106, 155130 (2022)
work page 2022
-
[45]
K. H. Ma, Y. J. Guo, L. Wang, and N. H. Tong, Phys. Rev. E 106, 014110 (2022)
work page 2022
-
[46]
H. W. Jia, W. J. Liu, Y. H. Wu, K. H. Ma, L. Wang, and N. H. Tong, Phys. Rev. B 111, 045153 (2025)
work page 2025
-
[47]
K. H. Ma and N. H. Tong, Phys. Rev. B 104, 155116 (2021)
work page 2021
-
[48]
D. J. Rowe, Phys. Rev. 175, 1283 (1968)
work page 1968
-
[49]
D. S. Delion, P. Schuck, and M. Tohyama, Eur. Phys. J. B 89, 45 (2016)
work page 2016
- [50]
-
[51]
P. Bleiziffer, A. Hesselmann, and A. G¨ orling, J. Chem. Phys. 139 084113 (2013)
work page 2013
-
[52]
M. Hellgren, F. Caruso, D. R. Rohr, X. Ren, A. Ru- bio, M. Scheffler, and P. Rinke, Phys. Rev. B 91, 165110 (2015)
work page 2015
-
[53]
E. Trushin, S. Fauser, A. M¨ olkner, J. Erhard, and A. G¨ orling, Phys. Rev. Lett.134, 016402 (2025)
work page 2025
-
[54]
Y. Jin, D. Zhang, Z. Chen, N. Q. Su, and W. Yang, J. Phys. Chem. Lett. 8 4746 (2017)
work page 2017
-
[55]
V. K. Voora, S. G. Balasubramani, and F. Furche, Phys. Rev. A 99, 012518 (2019)
work page 2019
-
[56]
J. M. Yu, B. D. Nguyen, J. Tsai, D. J. Hernandez, and F. Furche, J. Chem. Phys. 155, 040902 (2021)
work page 2021
-
[57]
A. Storozhenko, P. Schuck, J. Dukelsky, G. R¨ opke , and A. Vdovin, Ann. Phys. 307, 308 (2003)
work page 2003
- [58]
- [59]
- [60]
- [61]
- [62]
- [63]
-
[64]
J. S. Caux, H. Konno, M. Sorrell, and R. Weston, Phys. Rev. Lett 106, 217203 (2011)
work page 2011
-
[65]
J. S. Caux, H. Konno, M. Sorrell, and R. Weston, Stat. Mech, P01007 (2012)
work page 2012
- [66]
- [67]
-
[68]
Giamarchi, Quantum Physics in One Dimension Ox- ford University Press, Oxford, 2004, pp
T. Giamarchi, Quantum Physics in One Dimension Ox- ford University Press, Oxford, 2004, pp. 137–170
work page 2004
-
[69]
H. Barghathi, E. Casiano-Diaz, and A. Del Maestro, Phys. Rev. A 100, 022324(2019)
work page 2019
-
[70]
Sachdev, Quantum Phase Transitions Cambridge Uni- versity Press,Cambridge, 2011, pp.412-422
S. Sachdev, Quantum Phase Transitions Cambridge Uni- versity Press,Cambridge, 2011, pp.412-422
work page 2011
- [71]
-
[72]
A. J. Coleman and V. I. Yukalov, Reduced Density Ma- trices: Coulson ’s Challenge Springer, New York, 2000
work page 2000
-
[73]
D. A. Mazziotti, Phys. Rev. Lett. 130, 153001 (2023)
work page 2023
-
[74]
R. R. Li, M. D. Liebenthal, and A. E. DePrince III, J. Chem. Phys. 155, 174110 (2021)
work page 2021
-
[75]
J. G. Li, N. Michel, W. Zuo, and F. R. Xu, Phys. Rev. C 103, 064324 (2021)
work page 2021
- [76]
-
[77]
J. H. Cha, H. Y. Lee, and H. S. Kim, Scientific Reports, 1, 37490 (2025)
work page 2025
-
[78]
J. D. Baktay, A. E. Feiguin, and J. Rinc´ on, Phys. Rev. B 112, L161116 (2025)
work page 2025
- [79]
-
[80]
A. Menczer, K. Kap´ as, M. A. Werner, and ¨O. Legeza, Phys. Rev. B 109, 195148 (2024)
work page 2024
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.