Skip to content

Eigenstate Thermalization Hypothesis

The eigenstate thermalization hypothesis (ETH) is an ansatz for matrix elements of suitably simple observables within a resolved chaotic symmetry sector. Its diagonal part makes nearby energy eigenstates locally indistinguishable; its off-diagonal part suppresses late-time fluctuations. ETH can explain why a narrow-energy initial state relaxes to microcanonical values, but it is not a theorem for arbitrary Hamiltonians, eigenstates, or operators.

Required background. Equilibration and dephasing derives the diagonal ensemble whose thermal value ETH must explain. Statistical ensembles define the microcanonical comparison. Helpful background. Prethermalization supplies an important regime in which approximate charges delay the ETH prediction.

Let n\lvert n\rangle be energy eigenstates within one irreducible symmetry sector, and let AA be a few-body or local operator. With

Eˉ=Em+En2,ω=EmEn,\bar E=\frac{E_m+E_n}{2}, \qquad \omega=E_m-E_n,

the standard ansatz is

Amn=A(Eˉ)δmn+eS(Eˉ)/2fA(Eˉ,ω)Rmn.A_{mn} =A(\bar E)\delta_{mn} +e^{-S(\bar E)/2}f_A(\bar E,\omega)R_{mn}.

S(E)S(E) is the thermodynamic entropy of the sector and energy window, A(E)A(E) and fA(E,ω)f_A(E,\omega) are smooth on scales large compared with the level spacing, and RmnR_{mn} is a zero-mean unit-variance fluctuating quantity with Hermiticity correlations. This formulation originated in the random-matrix arguments of Deutsch 1991 and the statistical ansatz of Srednicki 1994.

The formula contains two logically different claims:

  • diagonal ETH: AnnA_{nn} approaches a smooth microcanonical function with shrinking fluctuations;
  • off-diagonal ETH: AmnA_{mn} has entropy-suppressed magnitude and a frequency envelope tied to equilibrium dynamics.

One can hold while a proposed detailed distribution for the other fails. The entropy factor also depends on whether discrete symmetries and conserved charges have been resolved; mixing sectors produces exact zeros, degeneracies, and artificial multimodality.

From eigenstate matrix elements to a quench prediction

Section titled “From eigenstate matrix elements to a quench prediction”

For ψ0=ncnn\lvert\psi_0\rangle=\sum_nc_n\lvert n\rangle, the long-time mean is

A=npnAnn,pn=cn2.\overline{\langle A\rangle} =\sum_n p_nA_{nn}, \qquad p_n=\lvert c_n\rvert^2.

Suppose the energy distribution is centered at E0E_0, has width ΔE\Delta E, lies in one sector, and A(E)A(E) is smooth across that width. Expanding gives

A=A(E0)+A(E0)npn(EnE0)+O ⁣(A(E0)ΔE2)+δETH.\overline{\langle A\rangle} =A(E_0) +A'(E_0)\sum_np_n(E_n-E_0) +O\!\left(A''(E_0)\Delta E^2\right) +\delta_{\mathrm{ETH}}.

The linear term vanishes by the definition of E0E_0. If diagonal ETH holds and ΔE\Delta E is thermodynamically narrow, A(E0)A(E_0) agrees with the microcanonical value up to ensemble and finite-size corrections.

For nondegenerate gaps, the temporal variance is

δA(t)2=mnpmpnAmn2.\overline{\lvert\delta A(t)\rvert^2} =\sum_{m\ne n}p_mp_n\lvert A_{mn}\rvert^2.

Off-diagonal ETH makes typical terms exponentially small in entropy. This explains small fluctuations, but the conclusion still depends on effective dimension, gap degeneracies, the support of the initial state, and possible rare matrix elements.

A reproducible test uses the following order.

  1. Resolve every exact symmetry and fix boundary conditions.
  2. Choose energy windows whose thermodynamic width shrinks while their level count grows.
  3. Compare AnnA_{nn} with a locally smoothed A(En)A(E_n) and report residual distributions, not only their mean.
  4. Test the residual width versus entropy or Hilbert-space dimension at fixed energy density.
  5. Bin eS(Eˉ)/2Amne^{S(\bar E)/2}A_{mn} in both Eˉ\bar E and ω\omega; test its variance and correlations.
  6. Propagate the measured diagonal residuals and initial energy width to the predicted quench value.
  7. Repeat for more than one local operator and inspect rare states separately.

Adjacent-level ratios can help locate an integrability-breaking crossover without unfolding, but Wigner–Dyson statistics are not a substitute for the matrix-element test. Conversely, ETH-like diagonal data for one operator do not establish spectral chaos.

This table is the canonical structured record for the chapter’s thermalization claims. Every row identifies a distinct observable and its strongest defensible conclusion; later pages link here rather than inventing a second standard.

ClaimObservable and sectorLimit or contourApproximation and mandatory checksFinite-size and reproducibility testStrongest supported conclusion
Local equilibrationA(t)\langle A(t)\rangle for declared local AA and initial statetime window before recurrence; thermodynamic order stateddiagonal ensemble, tolerance, conserved chargesvary volume, window, state, and observableselected observables spend most times near a stationary value
ETHAnnA_{nn} and AmnA_{mn} in one irreducible sectorfixed energy density, growing volumesmooth window, entropy normalization, rare-state searchcollapse residual distributions and propagate uncertainty to a quenchfinite-size evidence that the tested operator obeys ETH in that sector
Prethermalizationdrift of HH_* or approximate QjQ_jtfastttheatt_{\mathrm{fast}}\ll t\ll t_{\mathrm{heat}}controlled small parameter, resonance and remainder testsscale frequency/coupling, volume, and truncation orderdynamics is governed by an approximate generator during the stated window
Nonthermal scalingf(t,p)f(t,p), conserved integral, and fluxexpanding time and momentum windowcovariance-aware collapse, exponent relation, cutoff separationheld-out times and multiple initial conditionsself-similar transport compatible with the stated cascade
OTOC growthone fully specified regularized four-point functiondeclared thermal contour and pre-saturation intervaldisconnected subtraction, fit alternatives, operator dependencevary size, contour, operator, and fit windowthe chosen correlator has a resolved growth regime; chaos needs triangulation
Chaos-bound testanalytic normalized correlator and fitted λ\lambdadissipation–scrambling separation under theorem hypothesesstrip analyticity, boundedness, small connected correctiondemonstrate a parametrically stable window and hypothesis residualsλ2π/β\lambda\le2\pi/\beta only for the theorem’s correlator and regime
Operator spreadingsquared commutator or operator weight at separation xxcausal/front scaling limitfront marker, broadening model, saturation floorvary threshold, size, regulator, and operator paira butterfly/front velocity and broadening law for the stated system
Spectral chaosspacings and connected form factor in one sectorunfolded/filter-resolved bulk and late-time windowremove disconnected density, handle degeneracies and averaginggrow size and vary filter/unfolding and symmetry resolutionspectral correlations agree with a stated random-matrix class over resolved scales

No row alone proves universal continuum-QFT chaos. Closed unitary rows have no complete-positivity requirement; open-system guarantees are recorded separately on the master-equation evidence matrix.

Integrability breaking is evidence, not proof

Section titled “Integrability breaking is evidence, not proof”

Breaking an integrable Hamiltonian generally removes exact charges and can produce level repulsion, but the crossover scale depends on size, energy density, symmetry, perturbation, and observable. Constrained sectors, emergent conservation, localization, and scar towers can survive. The review by D’Alessio et al. 2016, §§4–7 collects numerical evidence and its finite-size limitations; it does not turn ETH into a general theorem.

The strongest current statement is conditional: many nonintegrable finite many-body models show convergent ETH signatures for broad classes of local observables. Whether a particular continuum limit, gauge sector, or unbounded operator satisfies the required scaling must be established in that theory.

Mixing symmetry sectors. Exact crossings and degeneracies then mimic nonchaotic statistics and distort matrix-element distributions.

Testing only diagonal means. ETH concerns fluctuations, off-diagonal elements, energy dependence, and operator class.

Ignoring rare states. A vanishing fraction can dominate a specially prepared quench. Report their number, overlap, and stability.

Calling ETH a proof of dynamics. The ansatz predicts outcomes when the initial energy distribution and gap structure cooperate; it does not fix the relaxation time.

Use the figure to place the ETH ansatz beside its principal exceptions and beside diagnostics that answer different questions. On this page, inspect the separation between matrix-element evidence for ETH and dynamical or spectral evidence for chaos.

The upper row places an ETH-compatible sector beside dephasing, prethermal, nonthermal-scaling, scar, and fragmented outcomes without arrows between them; the lower row separately lists operator-front, OTOC, Lyapunov-fit, and spectral evidence.

ETH is one possible account of late-time local observables, not a synonym for dephasing or chaos. A convincing ETH test combines diagonal smoothness, off-diagonal scaling, symmetry resolution, and an explicit search for exceptional states. The diagram is schematic: its upper alternatives and lower diagnostics are two independent, arrow-free rows with no universal time ordering.

In text-equivalent form, ETH concerns the energy dependence and entropy scaling of matrix elements. Operator spreading, regularized OTOC growth, level statistics, form factors, and recurrences instead test propagation or spectral correlations, so none substitutes for the matrix-element test.

Show that the off-diagonal ETH scale suppresses the temporal variance for a state spread over eSe^S comparable levels.

Solution

With pneSp_n\sim e^{-S} and Amn2eS\lvert A_{mn}\rvert^2\sim e^{-S}, there are e2Se^{2S} terms, each of size e3Se^{-3S}. The total is of order eSe^{-S}, up to the frequency envelope and correlations. Degenerate gaps or anomalously large elements invalidate this counting and must be tested.

A numerical study finds Gaussian off-diagonal elements but uses the full Hilbert space containing even and odd parity. Diagnose the test.

Solution

It is not an ETH test in a single irreducible sector. Parity selection makes some elements exactly zero and mixes two spectra. Resolve parity, normalize entropy within each sector, and repeat the diagonal and off-diagonal analysis separately.

Continue to exceptions and independent diagnostics

Section titled “Continue to exceptions and independent diagnostics”

Scars and fragmentation classify bounded failures of the typical-state expectation. Spectral statistics and OTOCs provide independent spectral and dynamical tests. Model-specific scars, fragmentation, and localization continue in Volume XII.