Green's Function-Free Formalism of Projective Truncation Approximation

Kou-Han Ma; Ning-Hua Tong; Yue-Hong Wu

arxiv: 2605.19404 · v1 · pith:CPUYU3IWnew · submitted 2026-05-19 · ❄️ cond-mat.str-el

Green's Function-Free Formalism of Projective Truncation Approximation

Kou-Han Ma , Yue-Hong Wu , Ning-Hua Tong This is my paper

Pith reviewed 2026-05-20 03:13 UTC · model grok-4.3

classification ❄️ cond-mat.str-el

keywords projective truncation approximationreduced density matricesGreen's functionsequation of motionmany-body systemsvariational methods

0 comments

The pith

Projective truncation approximation is reformulated as a self-consistent theory for reduced density matrices without using Green's functions.

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

The paper takes the projective truncation approximation, which was previously used to truncate the equation of motion for Green's functions in many-body systems, and rewrites it entirely in terms of reduced density matrices. This new formalism solves the equations by turning them into an over-constrained optimization problem that links directly to variational theories of density matrices. A reader might care because it removes dependence on Green's functions and their spectral theorem, while clarifying properties of the dynamical matrix and allowing discussion of generalized theorems like virial and Wick's in the density matrix language. It maintains the systematic nature of the truncation while opening connections to other RDM-based methods.

Core claim

In this work PTA is reformulated as a self-consistent theory for the reduced density matrices without reference to Green's functions. The properties of the dynamical matrix M are clarified and the solution of PTA equations is cast into an over-constrained optimization problem that connects to the variational RDM theory. Issues discussed include the scheme of alternative inner product, the generalized virial theorem, the generalized Wick's theorem, and the static component problem of PTA.

What carries the argument

The dynamical matrix M whose properties allow PTA equations to be solved as an over-constrained optimization problem for reduced density matrices.

If this is right

The PTA equations can now be solved self-consistently using only reduced density matrices.
The method connects directly to variational reduced density matrix theory.
Generalized virial theorem and generalized Wick's theorem apply within this RDM formalism.
The static component problem of PTA can be addressed under the new formulation.

Where Pith is reading between the lines

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

This reformulation may facilitate numerical implementations that optimize density matrices directly without computing Green's functions.
It could allow combining PTA with other variational RDM techniques for better approximations in strongly correlated electron systems.

Load-bearing premise

The dynamical matrix M possesses well-defined properties enabling the PTA equations to be treated as an over-constrained optimization without losing the controlled truncation.

What would settle it

Computing the physical quantities for a simple model Hamiltonian using both the original Green's function PTA and the new RDM optimization and finding inconsistent results would falsify the reformulation's validity.

read the original abstract

In previous works, the projected truncation approximation (PTA) was developed as a systematic and controlled method to truncate the equation of motion of Green's functions (GFs) for a given quantum or classical many-body Hamiltonian. The static averages are obtained self-consistently with the GF through the spectral theorem. In this work, PTA is reformulated as a self-consistent theory for the reduced density matrices (RDMs) without reference to GF. We separately discuss the issues of determining the dynamical matrix ${\bf M}$ and solving the physical quantities from it. The properties of ${\bf M}$ is clarified and the solution of PTA equations is cast into an over-constrained optimization problem. This makes connection of the present theory to the variational RDM theory. We discuss various issues of PTA under this formalism, including the scheme of alternative inner product, the generalized virial theorem, the generalized Wick's theorem, and the static component problem of PTA.

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.

Desk Editor's Note private letter to a colleague

This paper reformulates PTA as a Green's function-free RDM optimization but does not demonstrate that it recovers the same results as the original spectral-theorem version.

read the letter

Here's the quick take on arXiv:2605.19404. The paper takes the projected truncation approximation that was previously done with Green's functions and rewrites it as a Green's function-free formalism using reduced density matrices. The key move is turning the PTA equations into an over-constrained optimization problem that connects to variational RDM theory. They do a good job laying out the issues with determining the dynamical matrix M and solving for physical quantities from it. The discussion of properties of M, the scheme of alternative inner product, the generalized virial theorem, the generalized Wick's theorem, and the static component problem adds some clarity to how the approximation works under this new view. Where it falls short is in showing that this reformulation actually gives equivalent results to the original method. The stress-test highlights that equivalence for the dynamical matrix M is asserted but not demonstrated with proof or numerical check. Without that, it's not clear if the over-constrained problem keeps the systematic and controlled character of the truncation or if it introduces new issues. This paper is for people who work on many-body theory and truncation schemes in condensed matter. Someone following the PTA papers or looking at connections between different approximation methods might get something out of it. It shows honest engagement with the structure of the method. I recommend putting it through peer review. It's a solid technical reformulation that could benefit from referee comments on the details and the equivalence question.

Referee Report

2 major / 1 minor

Summary. The manuscript reformulates the projected truncation approximation (PTA) as a Green's function-free self-consistent theory for reduced density matrices (RDMs). It separates the determination of the dynamical matrix M from solving physical quantities, clarifies properties of M, and recasts the PTA equations as an over-constrained optimization problem that connects to variational RDM theory. Additional discussions cover alternative inner products, a generalized virial theorem, a generalized Wick's theorem, and the static component problem.

Significance. If the claimed equivalence holds and the optimization recovers the original controlled truncation, the reformulation could usefully link GF-based PTA to variational RDM methods, potentially enabling new solution strategies. The clarification of M's properties and the over-constrained framing are potentially valuable, but the significance is limited by the absence of explicit verification that static averages and hierarchy truncation are preserved.

major comments (2)

[discussion of dynamical matrix M] Discussion of determining the dynamical matrix M and solving physical quantities from it: the central claim that the over-constrained RDM optimization via M preserves the systematic truncation and static averages of the original spectral-theorem PTA is asserted without an explicit proof or numerical demonstration that solutions of the optimization recover the same hierarchy truncation or averages. This equivalence is load-bearing for the GF-free reformulation.
[over-constrained optimization] Section on the over-constrained optimization problem: it is unclear whether the additional constraints introduced by the optimization alter the controlled and systematic character of the original PTA truncation, as no concrete mapping back to the GF equations of motion or spectral theorem is provided to confirm invariance of the approximation.

minor comments (1)

[abstract] The abstract and introduction would benefit from a brief explicit statement of how the new RDM optimization reduces to the known PTA results in a simple test case (e.g., a two-site model) to anchor the reformulation.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading of our manuscript and the constructive comments. We address each major point below and will revise the manuscript accordingly to provide the requested explicit mappings and verifications.

read point-by-point responses

Referee: [discussion of dynamical matrix M] Discussion of determining the dynamical matrix M and solving physical quantities from it: the central claim that the over-constrained RDM optimization via M preserves the systematic truncation and static averages of the original spectral-theorem PTA is asserted without an explicit proof or numerical demonstration that solutions of the optimization recover the same hierarchy truncation or averages. This equivalence is load-bearing for the GF-free reformulation.

Authors: We agree that the manuscript would benefit from an explicit demonstration. The reformulation begins from the same truncated equations of motion used in the original PTA, with M constructed directly from the projected hierarchy; the optimization then enforces the identical algebraic relations that the spectral theorem would impose on the static averages. To make this transparent, we will add a dedicated subsection (or appendix) that derives the stationarity conditions of the optimization and shows term-by-term equivalence to the original GF equations of motion and the resulting static averages. This will include a short algebraic proof that any solution of the over-constrained problem satisfies the truncated hierarchy and the spectral-theorem relations. revision: yes
Referee: [over-constrained optimization] Section on the over-constrained optimization problem: it is unclear whether the additional constraints introduced by the optimization alter the controlled and systematic character of the original PTA truncation, as no concrete mapping back to the GF equations of motion or spectral theorem is provided to confirm invariance of the approximation.

Authors: We acknowledge that the current text does not supply an explicit back-mapping. The extra constraints in the optimization are not arbitrary; they are precisely the linear relations that arise when the original PTA truncation is substituted into the RDM equations. In the revision we will insert a concrete mapping that starts from the GF equations of motion, applies the same truncation, converts to RDM language, and arrives at the stationarity conditions of the optimization. This will demonstrate that the controlled and systematic character of the truncation is unchanged. revision: yes

Circularity Check

0 steps flagged

Reformulation of PTA to RDM optimization is self-contained without circular reduction

full rationale

The paper reformulates the projected truncation approximation as a GF-free self-consistent theory for reduced density matrices, separately addressing determination of the dynamical matrix M and casting the PTA equations as an over-constrained optimization problem that connects to variational RDM theory. Properties of M are clarified, and issues such as alternative inner products, generalized virial theorem, and generalized Wick's theorem are discussed under the new formalism. No load-bearing step reduces by construction to its inputs, fitted parameters renamed as predictions, or self-citation chains; the derivation develops the optimization independently of the original spectral-theorem approach while preserving the truncation hierarchy. The central claim remains independent of the referenced prior PTA works.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The paper introduces a dynamical matrix M and casts the solution as an over-constrained optimization, but details on any free parameters or new entities are not in the abstract.

axioms (1)

domain assumption Static averages are obtained self-consistently with the GF through the spectral theorem.
Stated in the abstract as part of the previous PTA works.

pith-pipeline@v0.9.0 · 5693 in / 1146 out tokens · 56471 ms · 2026-05-20T03:13:16.693809+00:00 · methodology

discussion (0)

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

IndisputableMonolith/Cost/FunctionalEquation.lean washburn_uniqueness_aczel unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

PTA is reformulated as a self-consistent theory for the reduced density matrices (RDMs) without reference to GF... cast into an over-constrained optimization problem. This makes connection of the present theory to the variational RDM theory.

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

71 extracted references · 71 canonical work pages · 1 internal anchor

[1]

Our present work represents a general and systematic theory in this direction

used the variational density matrix theory to produce the matrices equivalent to K and LB in the RPA equa- tions, providing an example of the application of the non- self-consistent inner product in the RPA context. Our present work represents a general and systematic theory in this direction. In the complete basis limit, the M obtained from Eq.(67) is st...

work page
[2]

The exact dynamical matrix M in this subspace reads M =    ǫ1 − t 0 0 − t ǫ 2 0 0 V 0 ǫ1 + V − t 0 V − t ǫ 2 + V   

Basis B1 = {a1, a 2, n 2a1, n 1a2} Writing down the commutator relations [a1, H ] = ǫ1a1 − ta2 + V n2a1, [a2, H ] = ǫ2a2 − ta1 + V n1a2, [n2a1, H ] = ( ǫ1 + V )n2a1 − tn1a2, [n1a2, H ] = ( ǫ2 + V )n1a2 − tn2a1, (80) we ﬁnd that this basis spans a closed subspace for L. The exact dynamical matrix M in this subspace reads M =    ǫ1 − t 0 0 − t ǫ 2 0 0 V ...

work page
[3]

Here, m ≡ n1− n2, h ≡ a† 1a2 − a† 2a1, p ≡ a† 1a2 + a† 2a1, and n = n1 + n2

Basis B2 = {m, h, p, n } Basis B2 that we consider is composed of the particle- hole excitation operators used in RPA. Here, m ≡ n1− n2, h ≡ a† 1a2 − a† 2a1, p ≡ a† 1a2 + a† 2a1, and n = n1 + n2. On this basis, the dynamical matrix is M =    0 − 2t 0 0 − 2t 0 ǫ2 − ǫ1 0 0 ǫ2 − ǫ1 0 0 0 0 0 0    . (85) The eigen excitation operators and the correspond...

work page
[4]

This equation is applicable to both the canonical and the grand canonical ensembles

O3 and O4 gives trivial equations since they are Her- mitian conserved quantities. This equation is applicable to both the canonical and the grand canonical ensembles. The conserving relations MT ⟨ ⃗A⟩ = 0 and CM = M†C give out Eq.(84). Therefore, for basis B2, PTA produces three independent equations for the ﬁve quantities. We can only solve them in the ...

work page
[5]

We have [a1a2, H ] = ( ǫ1 + ǫ2 + V )a1a2, (88) which gives M = ǫ1 + ǫ2 + V

Basis B3 = {a1a2} This one-dimensional basis forms a closed subspace for L of the spinless dimer model. We have [a1a2, H ] = ( ǫ1 + ǫ2 + V )a1a2, (88) which gives M = ǫ1 + ǫ2 + V . The equation Ceβ M = DT produces ⟨n1n2⟩eβ (ǫ1+ǫ2+V ) = 1 − ⟨ n1⟩ − ⟨ n2⟩ + ⟨n1n2⟩. (89) It applies only to the grand canonical ensemble. Equa- tions MT ⟨ ⃗A⟩ = 0 and CM = M†C a...

work page
[6]

But for even larger interaction, the iterative solution does not converge

0 or so, depending on the detailed iteration scheme. But for even larger interaction, the iterative solution does not converge. Apart from the factors of indeﬁnite com- mutator inner product and the static component prob- lem, we believe that the main reason for this is the over- constraint problem which somehow becomes more severe in the large interactio...

work page
[7]

A. L. Fetter and J. D. Walecka, Quantum theory of many- particle systems (Courier Corporation, 2012)

work page 2012
[8]

Aryasetiawan and O

F. Aryasetiawan and O. Gunnarsson, The GW method, Rep. Prog. Phys. 61, 237 (1998)

work page 1998
[9]

Onida, L

G. Onida, L. Reining, and A. Rubio, Electronic exci- tations: density-functional versus many-body Green’s- function approaches, Rev. Mod. Phys. 74, 601 (2002)

work page 2002
[10]

Mori, Transport, collective motion, and Brownian mo- tion, Prog

H. Mori, Transport, collective motion, and Brownian mo- tion, Prog. Theor. Phys. 33, 423 (1965)

work page 1965
[11]

Mori, A continued-fraction representation of the tim e- correlation functions, Prog

H. Mori, A continued-fraction representation of the tim e- correlation functions, Prog. Theor. Phys. 34, 399 (1965)

work page 1965
[12]

Y. A. Tserkovnikov, A method of solving inﬁnite sys- tems of equations for two-time thermal Green’s functions, Theor. Math. Phys. 49, 993 (1981)

work page 1981
[13]

M. H. Lee, Qrthogonalization process by recurrence rela - tions, Phys. Rev. Lett. 49, 1072 (1982)

work page 1982
[14]

Florencio and M

J. Florencio and M. H. L. Lee, Relaxation functions, memory functions, and random forces in the one- dimensional spin- 1 2 XY and transverse Ising model, Phys. Rev. B 35, 1835 (1987)

work page 1987
[15]

V. S. Viswanath and G. M¨ uller, The recursion method (Springer, 1994)

work page 1994
[16]

Kakehashi and P

Y. Kakehashi and P. Fulde, Coherent potential approxi- mation and projection operators for interacting electrons , Phys. Rev. B 69, 045101 (2004)

work page 2004
[17]

Kakehashi and P

Y. Kakehashi and P. Fulde, Self-consistent projection operator approach to nonlocal excitation spectra, Phys. Rev. B 70, 195102 (2004)

work page 2004
[18]

Mancini and A

F. Mancini and A. Avella, The Hubbard model within the equations of motion approach, Adv. Phys. 53, 537 (2004)

work page 2004
[19]

P. Fan, K. Yang, K. H. Ma, and N. H. Tong, Projective truncation approximation for equations of motion of two- time Green’s functions, Phys. Rev. B 97, 165140 (2018)

work page 2018
[20]

K. H. Ma, Y. J. Guo, L. Wang, and N. H. Tong, Projective-truncation-approximation study of the one- dimensional φ 4 lattice model, Phys. Rev. E 106, 014110 (2022). 21

work page 2022
[21]

D. J. Rowe, Equations-of-motion method and the ex- tended shell model, Rev. Mod. Phys. 40, 153 (1968)

work page 1968
[22]

D. J. Rowe, Methods for calculating ground-state corre la- tions of vibrational nuclei, Phys. Rev. 175, 1283 (1968)

work page 1968
[23]

Ring and P

P. Ring and P. Schuck, The Nuclear Many-Body Problem (Springer, New York, 1980)

work page 1980
[24]

Catara, G

F. Catara, G. Piccitto, and S. M., Towards a self- consistent random-phase approximation for Fermi sys- tems, Phys. Rev. B 54, 17536 (1996)

work page 1996
[25]

Chatterjee and K

K. Chatterjee and K. Perna, Excitation energies from ex - tended random phase approximation employed with ap- proximate one- and two-electron reduced density matri- ces, J. Chem. Phys. 137, 204109 (2012)

work page 2012
[26]

Fan and N

P. Fan and N. H. Tong, Controllable precision of the pro- jective truncation approximation for Green’s functions, Chin. Phys. B 28, 047102 (2019)

work page 2019
[27]

K. H. Ma and N. H. Tong, Interacting spinless fermions on the square lattice: Charge order, phase separation, and superconductivity, Phys. Rev. B 104, 155116 (2021)

work page 2021
[28]

P. Fan, N. H. Tong, and Z. G. Zhu, Spatial spin-spin correlations of the single-impurity Anderson model with a ferromagnetic bath, Phys. Rev. B 106, 155130 (2022)

work page 2022
[29]

H. W. Jia, W. J. Liu, Y. H. Wu, K. H. Ma, L. Wang, and N. H. Tong, Thermal broadening of the phonon spectral function in classical lattice models: Projective truncati on approximation, Phys. Rev. B 111, 045153 (2025)

work page 2025
[30]

Y. H. Wu, X. Ren, and N. H. Tong, Self-consistent ran- dom phase approximation from projective truncation ap- proximation formalism, Phys. Rev. B 113, 165131 (2026)

work page 2026
[31]

S. L. Brunton, M. Mudiˇ si´c, E. Kaiser, and J. N. Kutz, Modern Koopman theory for dyanmical systems, SIAM Review 64, 229 (2022)

work page 2022
[32]

Mauroy, I

A. Mauroy, I. Mezi´c, and Y. Suzuki, eds., The Koopman operator in systems and control: Concepts, methodolo- gies, and applications (Springer, 2020)

work page 2020
[33]

Schuck, D

P. Schuck, D. S. Delion, J. Dukelsky, M. Jemai, E. Litvi- nova, G. R¨ opke, and M. Tohyama, Equation of motion method to strongly correlated fermi systems and ex- tended RPA approaches, Phys. Reports 929, 1 (2021)

work page 2021
[34]

Dagotto, Correlated electrons in high-temperature su- perconductors, Rev

E. Dagotto, Correlated electrons in high-temperature su- perconductors, Rev. Mod. Phys. 66, 763 (1994)

work page 1994
[35]

Weiße and H

A. Weiße and H. Fehske, Exact diagonalization techniques, in Computational Many-Particle Physics (Springer, 2008) pp. 529–544

work page 2008
[36]

S. R. White, Density matrix formulation for quantum renormalization groups, Phys. Rev. Lett. 69, 2863 (1992)

work page 1992
[37]

Schollw¨ ock, The density-matrix renormalization group in the age of matrix product states, Ann

U. Schollw¨ ock, The density-matrix renormalization group in the age of matrix product states, Ann. Phys. 326, 96 (2011)

work page 2011
[38]

Levin and C

M. Levin and C. P. Nave, Tensor renormalization group approach to two-dimensional classical lattice models, Phys. Rev. Lett. 99, 120601 (2007)

work page 2007
[39]

Evenbly and G

G. Evenbly and G. Vidal, Tensor network renormaliza- tion, Phys. Rev. Lett. 115, 180405 (2015)

work page 2015
[40]

D. A. Mazziotti, Realization of quantum chemistry with - out wave functions through ﬁrst-order semideﬁnite pro- bramming, Phys. Rev. Lett. 93, 213001 (2004)

work page 2004
[41]

Mazziotti, Reduced-Density Matrix Mechanics: With Application to Many-Electron Atoms and Molecules, Ad- vances in Chemical Physis (Wiley, 2007)

D. Mazziotti, Reduced-Density Matrix Mechanics: With Application to Many-Electron Atoms and Molecules, Ad- vances in Chemical Physis (Wiley, 2007)

work page 2007
[42]

Eugene DePrince III, Variational determination of t he two-electron reduced density matrix: A tutorial review, WIRES COMPUT MOL SCI 14, e1702 (2024)

A. Eugene DePrince III, Variational determination of t he two-electron reduced density matrix: A tutorial review, WIRES COMPUT MOL SCI 14, e1702 (2024)

work page 2024
[43]

Li, Geometric constraints on two-electron reduced density matrices, Phys

Y. Li, Geometric constraints on two-electron reduced density matrices, Phys. Rev. A 103, 052202 (2021)

work page 2021
[44]

D. A. Mazziotti, Quantum many-body theory from a solution of the N -representability problem, Phys. Rev. Lett. 130, 153001 (2023)

work page 2023
[45]

A. O. Schouten, L. M. Sager Smith, and D. A. Mazz- iotti, Computation of exciton binding energies in exciton condensation, Phys. Rev. B 110, 035110 (2024)

work page 2024
[46]

A. O. Schouten, S. Ewing, and D. A. Mazziotti, Two- electron reduced density matrix method for the elec- tronic structure of correlated materials, Phys. Rev. B 112, 235130 (2025)

work page 2025
[47]

F. G. Scholtz, H. B. Geyer, and F. J. W. Nahne, Quasi- hermitian operators in quantum mechanics and the vari- ational principle, Ann. Phys. 213, 74 (1992)

work page 1992
[48]

R. J. Bartlett and M. Musia/suppress l, Coupled-cluster theory in quantum chemistry, Rev. Mod. Phys. 79, 291 (2007)

work page 2007
[49]

S. J. Wang and W. Cassing, Explicit treatment of N - body correlations within a density-matrix formalism, Ann. Phys. 159, 328 (1985)

work page 1985
[50]

D. N. Zubarev, Double-time Green functions in statisti cal physics, Sov. Phys. Usp. 3, 320 (1960)

work page 1960
[51]

Tohyama and P

M. Tohyama and P. Schuck, Spurious modes in extended RPA theories, Eur. Phys. J. A 19, 203 (2004)

work page 2004
[52]

A. J. Coleman, Structure of fermion density matrices, Rev. Mod. Phys. 35, 668 (1963)

work page 1963
[53]

Hassi and K

S. Hassi and K. Nordstr¨ orn, On projections in a space with an indeﬁnite metric, Linear Algebra and Its Appli- cations 208/209, 401 (1994)

work page 1994
[54]

Indefinite-metric quantum field theory and operator algebra

K. Kawamura, Indeﬁnite-metric quantum ﬁeld theory and operator algebra (2006), arXiv:math/0608076

work page internal anchor Pith review Pith/arXiv arXiv 2006
[55]

Belkasri and L

A. Belkasri and L. Richard, Damping of spin waves for a doped antiferromagnet, Phys. Rev. B 50, 12896 (1994)

work page 1994
[56]

Winterfeldt and D

S. Winterfeldt and D. Ihle, Theory of antiferromag- netic short-range order in the two-dimensional heisenberg model, Phys. Rev. B 56, 5535 (1997)

work page 1997
[57]

L. M. Roth, New method for linearizing many-body equa- tions of motion in statistical mechanics, Phys. Rev. Lett. 20, 1431 (1968)

work page 1968
[58]

Onoda and M

S. Onoda and M. Imada, Operator projection method ap- plied to the single-particle Green’s function in the Hub- bard model, J. Phys. Soc. Jpn. 70, 632 (2001)

work page 2001
[59]

Geipel and W

G. Geipel and W. Nolting, Ferromagnetism in the strongly correlated Hubbard model, Phys. Rev. B 38, 2608 (1988)

work page 1988
[60]

Nolting and W

W. Nolting and W. Borgiel, Band magnetism in the Hub- bard model, Phys. Rev. B 39, 6962 (1989)

work page 1989
[61]

Nolting and W

W. Nolting and W. D. Brewer, Fundamentals of many- body physics (Springer, 2009) Chap. 4

work page 2009
[62]

J. Li, N. D. Drummond, P. Schuck, and V. Olevano, Com- paring many-body approaches against the helium atom exact solution, Sci-Post Phys. 6, 040 (2019)

work page 2019
[63]

S. A. Ponomarenko, M. E. Sherril, D. P. Kilcrease, and G. Csanak, Statistical mean-ﬁeld theory of ﬁnite quan- tum systems: canonical ensemble formulation, J. Phys. A: Math. Gen. 39, L499 (2006)

work page 2006
[64]

Tsutsui and T

K. Tsutsui and T. Kita, Statistical mean-ﬁeld theory of ﬁnite quantum systems: canonical ensemble formulation, J. Phys. Soc. Jpn. 85, 114603 (2016)

work page 2016
[65]

Barghathi, J

H. Barghathi, J. Yu, and A. D. Maestro, Theory of nonin- teracting fermions and bosons in the canonical ensemble, Phys. Rev. Res. 2, 043206 (2020)

work page 2020
[66]

J. G. Hirsch, A. Mariano, J. Dukelsky, and P. Schuck, Fully self-consistent RPA description of the many level pairing model, Annals of Physics 296, 187 (2002). 22

work page 2002
[67]

Schuck, , and M

P. Schuck, , and M. Tohyama, Progress in many-body theory with the equation of motion method: Time- dependent density matrix meets self-consistent RPA and applications to solvable models, Phys. Rev. B 93, 165117 (2016)

work page 2016
[68]

Srednicki, The approach to thermal equilibrium in quantized chaotic systems, J

M. Srednicki, The approach to thermal equilibrium in quantized chaotic systems, J. Phys. A: Math. Gen. 32, 1163 (1999)

work page 1999
[69]

D’Alessio, Y

L. D’Alessio, Y. Kafri, A. Polkovnikov, and M. Rigol, From quantum chaos and eigenstate thermalization to statistical mechanics and thermodynamics, Adv. Phys. 65, 239 (2016)

work page 2016
[70]

Rekkedal, S

J. Rekkedal, S. Coriani, M. F. Iozzi, A. M. Teale, T. Hel- gaker, and T. B. Pedersen, Communication: Analytic gradients in the random-phase approximation, J. Chem. Phys. 139, 081101 (2013)

work page 2013
[71]

R. Song, X. Gong, and H. Z. Ye, Unphysical solutions in coupled-cluster random phase approximation and how to avoid them (2025), arXiv: 2508.11051

work page arXiv 2025

[1] [1]

Our present work represents a general and systematic theory in this direction

used the variational density matrix theory to produce the matrices equivalent to K and LB in the RPA equa- tions, providing an example of the application of the non- self-consistent inner product in the RPA context. Our present work represents a general and systematic theory in this direction. In the complete basis limit, the M obtained from Eq.(67) is st...

work page

[2] [2]

The exact dynamical matrix M in this subspace reads M =    ǫ1 − t 0 0 − t ǫ 2 0 0 V 0 ǫ1 + V − t 0 V − t ǫ 2 + V   

Basis B1 = {a1, a 2, n 2a1, n 1a2} Writing down the commutator relations [a1, H ] = ǫ1a1 − ta2 + V n2a1, [a2, H ] = ǫ2a2 − ta1 + V n1a2, [n2a1, H ] = ( ǫ1 + V )n2a1 − tn1a2, [n1a2, H ] = ( ǫ2 + V )n1a2 − tn2a1, (80) we ﬁnd that this basis spans a closed subspace for L. The exact dynamical matrix M in this subspace reads M =    ǫ1 − t 0 0 − t ǫ 2 0 0 V ...

work page

[3] [3]

Here, m ≡ n1− n2, h ≡ a† 1a2 − a† 2a1, p ≡ a† 1a2 + a† 2a1, and n = n1 + n2

Basis B2 = {m, h, p, n } Basis B2 that we consider is composed of the particle- hole excitation operators used in RPA. Here, m ≡ n1− n2, h ≡ a† 1a2 − a† 2a1, p ≡ a† 1a2 + a† 2a1, and n = n1 + n2. On this basis, the dynamical matrix is M =    0 − 2t 0 0 − 2t 0 ǫ2 − ǫ1 0 0 ǫ2 − ǫ1 0 0 0 0 0 0    . (85) The eigen excitation operators and the correspond...

work page

[4] [4]

This equation is applicable to both the canonical and the grand canonical ensembles

O3 and O4 gives trivial equations since they are Her- mitian conserved quantities. This equation is applicable to both the canonical and the grand canonical ensembles. The conserving relations MT ⟨ ⃗A⟩ = 0 and CM = M†C give out Eq.(84). Therefore, for basis B2, PTA produces three independent equations for the ﬁve quantities. We can only solve them in the ...

work page

[5] [5]

We have [a1a2, H ] = ( ǫ1 + ǫ2 + V )a1a2, (88) which gives M = ǫ1 + ǫ2 + V

Basis B3 = {a1a2} This one-dimensional basis forms a closed subspace for L of the spinless dimer model. We have [a1a2, H ] = ( ǫ1 + ǫ2 + V )a1a2, (88) which gives M = ǫ1 + ǫ2 + V . The equation Ceβ M = DT produces ⟨n1n2⟩eβ (ǫ1+ǫ2+V ) = 1 − ⟨ n1⟩ − ⟨ n2⟩ + ⟨n1n2⟩. (89) It applies only to the grand canonical ensemble. Equa- tions MT ⟨ ⃗A⟩ = 0 and CM = M†C a...

work page

[6] [6]

But for even larger interaction, the iterative solution does not converge

0 or so, depending on the detailed iteration scheme. But for even larger interaction, the iterative solution does not converge. Apart from the factors of indeﬁnite com- mutator inner product and the static component prob- lem, we believe that the main reason for this is the over- constraint problem which somehow becomes more severe in the large interactio...

work page

[7] [7]

A. L. Fetter and J. D. Walecka, Quantum theory of many- particle systems (Courier Corporation, 2012)

work page 2012

[8] [8]

Aryasetiawan and O

F. Aryasetiawan and O. Gunnarsson, The GW method, Rep. Prog. Phys. 61, 237 (1998)

work page 1998

[9] [9]

Onida, L

G. Onida, L. Reining, and A. Rubio, Electronic exci- tations: density-functional versus many-body Green’s- function approaches, Rev. Mod. Phys. 74, 601 (2002)

work page 2002

[10] [10]

Mori, Transport, collective motion, and Brownian mo- tion, Prog

H. Mori, Transport, collective motion, and Brownian mo- tion, Prog. Theor. Phys. 33, 423 (1965)

work page 1965

[11] [11]

Mori, A continued-fraction representation of the tim e- correlation functions, Prog

H. Mori, A continued-fraction representation of the tim e- correlation functions, Prog. Theor. Phys. 34, 399 (1965)

work page 1965

[12] [12]

Y. A. Tserkovnikov, A method of solving inﬁnite sys- tems of equations for two-time thermal Green’s functions, Theor. Math. Phys. 49, 993 (1981)

work page 1981

[13] [13]

M. H. Lee, Qrthogonalization process by recurrence rela - tions, Phys. Rev. Lett. 49, 1072 (1982)

work page 1982

[14] [14]

Florencio and M

J. Florencio and M. H. L. Lee, Relaxation functions, memory functions, and random forces in the one- dimensional spin- 1 2 XY and transverse Ising model, Phys. Rev. B 35, 1835 (1987)

work page 1987

[15] [15]

V. S. Viswanath and G. M¨ uller, The recursion method (Springer, 1994)

work page 1994

[16] [16]

Kakehashi and P

Y. Kakehashi and P. Fulde, Coherent potential approxi- mation and projection operators for interacting electrons , Phys. Rev. B 69, 045101 (2004)

work page 2004

[17] [17]

Kakehashi and P

Y. Kakehashi and P. Fulde, Self-consistent projection operator approach to nonlocal excitation spectra, Phys. Rev. B 70, 195102 (2004)

work page 2004

[18] [18]

Mancini and A

F. Mancini and A. Avella, The Hubbard model within the equations of motion approach, Adv. Phys. 53, 537 (2004)

work page 2004

[19] [19]

P. Fan, K. Yang, K. H. Ma, and N. H. Tong, Projective truncation approximation for equations of motion of two- time Green’s functions, Phys. Rev. B 97, 165140 (2018)

work page 2018

[20] [20]

K. H. Ma, Y. J. Guo, L. Wang, and N. H. Tong, Projective-truncation-approximation study of the one- dimensional φ 4 lattice model, Phys. Rev. E 106, 014110 (2022). 21

work page 2022

[21] [21]

D. J. Rowe, Equations-of-motion method and the ex- tended shell model, Rev. Mod. Phys. 40, 153 (1968)

work page 1968

[22] [22]

D. J. Rowe, Methods for calculating ground-state corre la- tions of vibrational nuclei, Phys. Rev. 175, 1283 (1968)

work page 1968

[23] [23]

Ring and P

P. Ring and P. Schuck, The Nuclear Many-Body Problem (Springer, New York, 1980)

work page 1980

[24] [24]

Catara, G

F. Catara, G. Piccitto, and S. M., Towards a self- consistent random-phase approximation for Fermi sys- tems, Phys. Rev. B 54, 17536 (1996)

work page 1996

[25] [25]

Chatterjee and K

K. Chatterjee and K. Perna, Excitation energies from ex - tended random phase approximation employed with ap- proximate one- and two-electron reduced density matri- ces, J. Chem. Phys. 137, 204109 (2012)

work page 2012

[26] [26]

Fan and N

P. Fan and N. H. Tong, Controllable precision of the pro- jective truncation approximation for Green’s functions, Chin. Phys. B 28, 047102 (2019)

work page 2019

[27] [27]

K. H. Ma and N. H. Tong, Interacting spinless fermions on the square lattice: Charge order, phase separation, and superconductivity, Phys. Rev. B 104, 155116 (2021)

work page 2021

[28] [28]

P. Fan, N. H. Tong, and Z. G. Zhu, Spatial spin-spin correlations of the single-impurity Anderson model with a ferromagnetic bath, Phys. Rev. B 106, 155130 (2022)

work page 2022

[29] [29]

H. W. Jia, W. J. Liu, Y. H. Wu, K. H. Ma, L. Wang, and N. H. Tong, Thermal broadening of the phonon spectral function in classical lattice models: Projective truncati on approximation, Phys. Rev. B 111, 045153 (2025)

work page 2025

[30] [30]

Y. H. Wu, X. Ren, and N. H. Tong, Self-consistent ran- dom phase approximation from projective truncation ap- proximation formalism, Phys. Rev. B 113, 165131 (2026)

work page 2026

[31] [31]

S. L. Brunton, M. Mudiˇ si´c, E. Kaiser, and J. N. Kutz, Modern Koopman theory for dyanmical systems, SIAM Review 64, 229 (2022)

work page 2022

[32] [32]

Mauroy, I

A. Mauroy, I. Mezi´c, and Y. Suzuki, eds., The Koopman operator in systems and control: Concepts, methodolo- gies, and applications (Springer, 2020)

work page 2020

[33] [33]

Schuck, D

P. Schuck, D. S. Delion, J. Dukelsky, M. Jemai, E. Litvi- nova, G. R¨ opke, and M. Tohyama, Equation of motion method to strongly correlated fermi systems and ex- tended RPA approaches, Phys. Reports 929, 1 (2021)

work page 2021

[34] [34]

Dagotto, Correlated electrons in high-temperature su- perconductors, Rev

E. Dagotto, Correlated electrons in high-temperature su- perconductors, Rev. Mod. Phys. 66, 763 (1994)

work page 1994

[35] [35]

Weiße and H

A. Weiße and H. Fehske, Exact diagonalization techniques, in Computational Many-Particle Physics (Springer, 2008) pp. 529–544

work page 2008

[36] [36]

S. R. White, Density matrix formulation for quantum renormalization groups, Phys. Rev. Lett. 69, 2863 (1992)

work page 1992

[37] [37]

Schollw¨ ock, The density-matrix renormalization group in the age of matrix product states, Ann

U. Schollw¨ ock, The density-matrix renormalization group in the age of matrix product states, Ann. Phys. 326, 96 (2011)

work page 2011

[38] [38]

Levin and C

M. Levin and C. P. Nave, Tensor renormalization group approach to two-dimensional classical lattice models, Phys. Rev. Lett. 99, 120601 (2007)

work page 2007

[39] [39]

Evenbly and G

G. Evenbly and G. Vidal, Tensor network renormaliza- tion, Phys. Rev. Lett. 115, 180405 (2015)

work page 2015

[40] [40]

D. A. Mazziotti, Realization of quantum chemistry with - out wave functions through ﬁrst-order semideﬁnite pro- bramming, Phys. Rev. Lett. 93, 213001 (2004)

work page 2004

[41] [41]

Mazziotti, Reduced-Density Matrix Mechanics: With Application to Many-Electron Atoms and Molecules, Ad- vances in Chemical Physis (Wiley, 2007)

D. Mazziotti, Reduced-Density Matrix Mechanics: With Application to Many-Electron Atoms and Molecules, Ad- vances in Chemical Physis (Wiley, 2007)

work page 2007

[42] [42]

Eugene DePrince III, Variational determination of t he two-electron reduced density matrix: A tutorial review, WIRES COMPUT MOL SCI 14, e1702 (2024)

A. Eugene DePrince III, Variational determination of t he two-electron reduced density matrix: A tutorial review, WIRES COMPUT MOL SCI 14, e1702 (2024)

work page 2024

[43] [43]

Li, Geometric constraints on two-electron reduced density matrices, Phys

Y. Li, Geometric constraints on two-electron reduced density matrices, Phys. Rev. A 103, 052202 (2021)

work page 2021

[44] [44]

D. A. Mazziotti, Quantum many-body theory from a solution of the N -representability problem, Phys. Rev. Lett. 130, 153001 (2023)

work page 2023

[45] [45]

A. O. Schouten, L. M. Sager Smith, and D. A. Mazz- iotti, Computation of exciton binding energies in exciton condensation, Phys. Rev. B 110, 035110 (2024)

work page 2024

[46] [46]

A. O. Schouten, S. Ewing, and D. A. Mazziotti, Two- electron reduced density matrix method for the elec- tronic structure of correlated materials, Phys. Rev. B 112, 235130 (2025)

work page 2025

[47] [47]

F. G. Scholtz, H. B. Geyer, and F. J. W. Nahne, Quasi- hermitian operators in quantum mechanics and the vari- ational principle, Ann. Phys. 213, 74 (1992)

work page 1992

[48] [48]

R. J. Bartlett and M. Musia/suppress l, Coupled-cluster theory in quantum chemistry, Rev. Mod. Phys. 79, 291 (2007)

work page 2007

[49] [49]

S. J. Wang and W. Cassing, Explicit treatment of N - body correlations within a density-matrix formalism, Ann. Phys. 159, 328 (1985)

work page 1985

[50] [50]

D. N. Zubarev, Double-time Green functions in statisti cal physics, Sov. Phys. Usp. 3, 320 (1960)

work page 1960

[51] [51]

Tohyama and P

M. Tohyama and P. Schuck, Spurious modes in extended RPA theories, Eur. Phys. J. A 19, 203 (2004)

work page 2004

[52] [52]

A. J. Coleman, Structure of fermion density matrices, Rev. Mod. Phys. 35, 668 (1963)

work page 1963

[53] [53]

Hassi and K

S. Hassi and K. Nordstr¨ orn, On projections in a space with an indeﬁnite metric, Linear Algebra and Its Appli- cations 208/209, 401 (1994)

work page 1994

[54] [54]

Indefinite-metric quantum field theory and operator algebra

K. Kawamura, Indeﬁnite-metric quantum ﬁeld theory and operator algebra (2006), arXiv:math/0608076

work page internal anchor Pith review Pith/arXiv arXiv 2006

[55] [55]

Belkasri and L

A. Belkasri and L. Richard, Damping of spin waves for a doped antiferromagnet, Phys. Rev. B 50, 12896 (1994)

work page 1994

[56] [56]

Winterfeldt and D

S. Winterfeldt and D. Ihle, Theory of antiferromag- netic short-range order in the two-dimensional heisenberg model, Phys. Rev. B 56, 5535 (1997)

work page 1997

[57] [57]

L. M. Roth, New method for linearizing many-body equa- tions of motion in statistical mechanics, Phys. Rev. Lett. 20, 1431 (1968)

work page 1968

[58] [58]

Onoda and M

S. Onoda and M. Imada, Operator projection method ap- plied to the single-particle Green’s function in the Hub- bard model, J. Phys. Soc. Jpn. 70, 632 (2001)

work page 2001

[59] [59]

Geipel and W

G. Geipel and W. Nolting, Ferromagnetism in the strongly correlated Hubbard model, Phys. Rev. B 38, 2608 (1988)

work page 1988

[60] [60]

Nolting and W

W. Nolting and W. Borgiel, Band magnetism in the Hub- bard model, Phys. Rev. B 39, 6962 (1989)

work page 1989

[61] [61]

Nolting and W

W. Nolting and W. D. Brewer, Fundamentals of many- body physics (Springer, 2009) Chap. 4

work page 2009

[62] [62]

J. Li, N. D. Drummond, P. Schuck, and V. Olevano, Com- paring many-body approaches against the helium atom exact solution, Sci-Post Phys. 6, 040 (2019)

work page 2019

[63] [63]

S. A. Ponomarenko, M. E. Sherril, D. P. Kilcrease, and G. Csanak, Statistical mean-ﬁeld theory of ﬁnite quan- tum systems: canonical ensemble formulation, J. Phys. A: Math. Gen. 39, L499 (2006)

work page 2006

[64] [64]

Tsutsui and T

K. Tsutsui and T. Kita, Statistical mean-ﬁeld theory of ﬁnite quantum systems: canonical ensemble formulation, J. Phys. Soc. Jpn. 85, 114603 (2016)

work page 2016

[65] [65]

Barghathi, J

H. Barghathi, J. Yu, and A. D. Maestro, Theory of nonin- teracting fermions and bosons in the canonical ensemble, Phys. Rev. Res. 2, 043206 (2020)

work page 2020

[66] [66]

J. G. Hirsch, A. Mariano, J. Dukelsky, and P. Schuck, Fully self-consistent RPA description of the many level pairing model, Annals of Physics 296, 187 (2002). 22

work page 2002

[67] [67]

Schuck, , and M

P. Schuck, , and M. Tohyama, Progress in many-body theory with the equation of motion method: Time- dependent density matrix meets self-consistent RPA and applications to solvable models, Phys. Rev. B 93, 165117 (2016)

work page 2016

[68] [68]

Srednicki, The approach to thermal equilibrium in quantized chaotic systems, J

M. Srednicki, The approach to thermal equilibrium in quantized chaotic systems, J. Phys. A: Math. Gen. 32, 1163 (1999)

work page 1999

[69] [69]

D’Alessio, Y

L. D’Alessio, Y. Kafri, A. Polkovnikov, and M. Rigol, From quantum chaos and eigenstate thermalization to statistical mechanics and thermodynamics, Adv. Phys. 65, 239 (2016)

work page 2016

[70] [70]

Rekkedal, S

J. Rekkedal, S. Coriani, M. F. Iozzi, A. M. Teale, T. Hel- gaker, and T. B. Pedersen, Communication: Analytic gradients in the random-phase approximation, J. Chem. Phys. 139, 081101 (2013)

work page 2013

[71] [71]

R. Song, X. Gong, and H. Z. Ye, Unphysical solutions in coupled-cluster random phase approximation and how to avoid them (2025), arXiv: 2508.11051

work page arXiv 2025