The exact multi-resolvent hierarchy of dense many-body spectra encodes higher-order spectral-weight correlations, but closing it predictively for interacting Hamiltonians remains a challenge. We introduce the multi-resolvent self-consistency (MRSC) equations, which separate an exact algebraic hierarchy from an approximate finite-parameter closure. The exact Feshbach geometric identity organizes r-resolvent cumulants into a ladder over irreducible vertices Gamma^(r), and the projected hierarchy expresses those vertices as functionals of the cumulants. The closure works in finite parameters: a connectivity basis organizes a general Hamiltonian's seed into class amplitudes, and a causal r-fold Stieltjes representation with a moment-field identity reduce moment matching to a finite-dimensional problem closed at a predictive fixed point. On the free-fermion benchmark the tested exact identities are reproduced and the Faddeeva closure is recovered at r=1. In a blind test on the interacting Ising chain, with the exact two-resolvent data withheld, the Hamiltonian-derived closure predicts the tested probes to 0.35%-0.97%. At the observable level the exact two-resolvent sector supplies a sign-correct, energy-resolved correction to the single-resolvent baseline for the off-diagonal eigenstate-thermalization-hypothesis (ETH) variance, of order one relative to it. Closure errors organize hierarchically: a level-by-level deficit decomposition, anchor measurements of the isolatable terms, and an explicit propagation bound. The same moments further yield a diagnostic chain for the window conditions of local rotational invariance -- an exact drift-kernel identity, a proven single-state identifiability boundary, and a transferability census over couplings and realizations -- a finite-order framework, not a proof of universal emergence.
The exact multi-resolvent hierarchy of dense many-body spectra encodes higher-order spectral-weight correlations, but closing it predictively for interacting Hamiltonians remains a challenge. We introduce the multi-resolvent self-consistency (MRSC) equations, which separate an exact algebraic hierarchy from an approximate finite-parameter closure. The exact Feshbach geometric identity organizes r-resolvent cumulants into a ladder over irreducible vertices Gamma^(r), and the projected hierarchy expresses those vertices as functionals of the cumulants. The closure works in finite parameters: a connectivity basis organizes a general Hamiltonian's seed into class amplitudes, and a causal r-fold Stieltjes representation with a moment-field identity reduce moment matching to a finite-dimensional problem closed at a predictive fixed point. On the free-fermion benchmark the tested exact identities are reproduced and the Faddeeva closure is recovered at r=1. In a blind test on the interacting Ising chain, with the exact two-resolvent data withheld, the Hamiltonian-derived closure predicts the tested probes to 0.35%-0.97%. At the observable level the exact two-resolvent sector supplies a sign-correct, energy-resolved correction to the single-resolvent baseline for the off-diagonal eigenstate-thermalization-hypothesis (ETH) variance, of order one relative to it. Closure errors organize hierarchically: a level-by-level deficit decomposition, anchor measurements of the isolatable terms, and an explicit propagation bound. The same moments further yield a diagnostic chain for the window conditions of local rotational invariance -- an exact drift-kernel identity, a proven single-state identifiability boundary, and a transferability census over couplings and realizations -- a finite-order framework, not a proof of universal emergence.
We develop a nonperturbative framework for generic nonintegrable many-body systems that reorganizes the expansion of diagonal Green's functions. Starting from exact projection identities and the spectral representation of the resolvent, we derive a recursive hierarchy for the self-energy in which cross-correlated propagation processes are systematically rewritten in terms of diagonal resolvents. Under a diagonal closure approximation, the hierarchy becomes formally closed yet remains systematically improvable.
The framework combines two nonperturbative mechanisms. First, a Lanczos continued-fraction representation provides a recursive single-resolvent structure that naturally produces non-Lorentzian spectral features beyond self-consistent Born approximations. Second, an exact projected multi-resolvent hierarchy introduces nonlocal frequency couplings through products of resolvents and their Hilbert transforms. These contributions mix parity sectors under energy reflection and generate spectral skewness, which is absent in single-resolvent closures.
To solve the resulting equations, we employ a hierarchy of Lorentzian, Gaussian, and Voigt-type ansätze together with an effective Faddeeva self-energy representation ensuring analyticity and causality. Spectral broadening, distribution tails, and higher-order fluctuations emerge from the interplay between continued-fraction recursion and multi-resolvent correlations.
The framework requires no small expansion parameters or diagrammatic truncations, relying instead on ETH-type statistical assumptions appropriate for dense chaotic spectra. It provides a unified route from microscopic interactions to emergent spectral structure, revealing a progression from single-pole self-consistent dynamics to continued-fraction renormalization and finally to multi-resolvent interference effects.
The spectral weight with which a many-body eigenstate contributes to a given channel is not fixed by the smooth one-point ETH envelope, and its higher-order correlations---the connected cumulants of three and more channel-resolved spectral functions---have lacked a systematic organization. We construct one. In random free fermions, where every channel overlap is a squared Slater minor of a Haar-orthogonal one-body eigenvector matrix, the r-resolvent moment admits an exact replica-sector decomposition: a sum over the double-coset sectors of the orthogonal Weingarten algebra, in which the permutation sectors are the replica contraction classes and each sector value is determined by the overlap pattern of the channels and of the eigenstates. The sector coefficients follow from an exact occupancy-matrix assembly---a signed convolution over the slot permutations, a decomposition into connected blocks, and a finite-state chain transfer matrix whose local transition rules are independent of r and k---so the sector theorem holds at arbitrary finite multiplicity, realized exactly at r=3 (k=2,3,4), r=4 and r=5 (k=2). The structure carries the physical hierarchy: the fully coincident sector is the projection that yields the Selberg moment ratios qr, the levels are nested by exact marginalization, and the sectors are the multi-energy spectral cumulants behind the resolvent cumulants, with the irreducible self-energy vertices of the Feshbach ladder as the connected kernel. The multi-resolvent correlations of random free fermions are therefore organized, level by level, by one replica-sector algebra---an exact microscopic benchmark for the higher-order fluctuation sector of eigenstate thermalization, with its microscopic sector geometry made explicit.
Eigenstate thermalization fixes the smooth energy-resolved envelope of few-body observables but not the microscopic overlap statistics that realize it. In the multi-resolvent hierarchy, this missing information appears in a two-resolvent fluctuation sector. We solve that sector exactly in a tractable benchmark, random free fermions, obtaining a non-Porter-Thomas intensity hierarchy, a structured channel-distance covariance, and an exact overlap kernel that closes the two-resolvent covariance. We then reconstruct the same sector without using its exact solution: a projected self-consistency scheme fixes the irreducible vertex through independently computable projections, while unconstrained residuals identify the sectors missed by a minimal ansatz and are completed exactly by a finite-size Weingarten evaluation, without fitting. For interacting deformations, we prove a model-independent moment identity showing that the deep-plane coefficients are fixed by Hamiltonian-moment covariances. This yields an exact rigidity law: pair-hopping interactions freeze the sector, dense perturbations add only an isotropic layer, and density-density interactions produce an explicit polynomial deformation in channel geometry. By contrast, the bulk sector requires the interacting eigenstate-overlap kernel and exhibits a logarithmically divergent perturbative response. The fluctuation sector therefore splits into a rigid moment-determined part and an intrinsically nonperturbative eigenstate-resolved part.
The eigenstate thermalization hypothesis (ETH) parametrizes off-diagonal matrix elements by a smooth function whose microscopic origin remains largely phenomenological. We develop a multi-resolvent hierarchy that derives this smooth structure from the microscopic Hamiltonian. Starting from exact projection identities, we express the ETH variance as fji2=Dji+gji, where Dji is the diagonal-overlap baseline and gji=∑r≥2gji(r) is a systematically improvable hierarchy of multi-channel interference processes. The leading r=3 sector generates an odd-parity component in the energy difference, inaccessible to parity-preserving single-resolvent closures. The same resolvent construction yields an exact covariance representation of eigenstate fluctuations. Under amplitude isotropy, decorrelation, and regularity, it reduces to the Gaussian limit, with q−3=κ4/⟨c2⟩2; normalization further fixes the cross-channel covariance Cˉi, including its energy-resolved form. Exact diagonalization verifies these relations in random-matrix and structured systems, while the latter exhibit controlled breakdown of the isotropic Gaussian closure. The framework thus provides a microscopic hierarchy for both the smooth and fluctuation sectors of subsystem ETH.
Hidden-free-fermion (HFF) Hamiltonians exhibit free-fermion spectra despite lacking an explicit quadratic representation, yet their algebraic provenance and periodic-boundary behavior remain open questions. Here we establish that a broad class of HFF models arises from distinct Lax operator realizations of a twisted multistate Perk--Schultz R matrix. Its reflected degeneration generates a closed finite hierarchy of mutually commuting transfer matrices, whose functional relations determine the complete many-body spectrum. As a decisive example, we solve the inhomogeneous free-fermions-in-disguise (FFD) chain, which exhibits HFF, nested HFF, or interacting spectra, depending on its coupling. The homogeneous FFD chain with both one-wrap and periodic boundary conditions exhibits an emergent cubic finite-size gap while retaining the exact z=3/2 bulk behavior. The construction extends to higher-state and face-type realizations, providing a common algebraic framework for coupling-dependent spectral organization and multiscale behavior in interacting quantum systems.
We investigate the generic emergence of entropy in isolated quantum systems from the spectral structure of their long-time dephased states. Using a resolvent-based hierarchy, we separate the smooth energy-shell envelope of eigenstate intensities from their intra-shell fluctuations. The single-resolvent sector determines the coarse-grained intensity envelope and yields an entropy contribution that is quantitatively captured by a Lorentzian-Gaussian form. We then show that the difference between this envelope entropy and the von Neumann entropy is governed exactly by the conditional statistics of the residual intensity fluctuations. The two-resolvent covariance sector identifies their leading variance, while higher correlation sectors control the full non-Gaussian correction. A Gaussian-amplitude closure provides a benchmark entropy deficit, whereas exact diagonalization of an interacting Ising chain exhibits substantial finite-size deviations from this benchmark. These results establish a hierarchical spectral framework in which entropy beyond the smooth envelope is directly tied to multi-resolvent correlations, providing a systematic route toward its microscopic characterization.