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.
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 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.
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.
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.
We analyze thermalization and eigenstate thermalization in the random free fermion model, a quadratic Hamiltonian with a Gaussian Orthogonal Ensemble hopping matrix, by an exact random-matrix calculation. For a single-particle observable A we compute the ensemble average and fluctuations of Tr(Aρ(t)): we derive the relaxation envelope g2(t/τλ) with τλ=ℏ/(2ηN) and its oscillatory ∣t∣−3 tail, both properties of the semicircular limiting density, and we show that the variance of the fluctuations is O(1/N), so that the relaxation profile self-averages at fixed scaled time. Incorporating particle-number conservation, we prove a factorization theorem for the energy-shell statistics of the many-body eigenstates: the energy-independent ETH structure follows from the statistical independence of the eigenvector matrix and the spectrum, and at fixed realization the energy constraint enters only through a covariance subtraction of relative order 1/N. A controlled eigenvector--eigenvalue correlation deformation H(ε)=H0+εH1 then generates, at first order in ε, the energy dependence of the diagonal energy-resolved slope and of the off-diagonal spectral function, a single-particle diagnostic of the ETH structure in the allowed particle-hole sector. Correlation-function fluctuations and the single-particle spectral form factor finally identify the subleading O(1/N) late-time fingerprints of two-level spectral correlations. The model thus separates eigenbasis chaos from spectral chaos and identifies the statistical origin of energy-independent versus energy-dependent ETH.
We study diagnostics of thermalization in quantum many-body systems with global SU(2) symmetry, where the standard eigenstate thermalization hypothesis (ETH) is generalized to its non-Abelian form. As an eigenstate-level probe, we introduce a symmetry-resolved trace distance constructed from the block structure of the reduced density matrix. This block structure separates spin-sector probabilities from configurational fluctuations within each sector, naturally leading to a decomposition into a probability trace distance and a configurational trace distance. The microcanonical average of the former is bounded by fluctuations of the corresponding spin-sector probabilities within a microcanonical energy window, whereas the latter captures finer intra-sector fluctuations. In non-Abelian thermalizing systems, these spin-sector-probability fluctuations are constrained by the non-Abelian ETH and therefore become exponentially suppressed with system size. Numerical studies of the one-dimensional J1--J2 Heisenberg chain are consistent with this picture and suggest that, in the thermal regime, the trace distance is asymptotically dominated by the configurational trace distance.