---
title: Fock Space Prethermalization (FSP)
url: https://www.emergentmind.com/topics/fock-space-prethermalization-fsp
type: topic
---

# Fock Space Prethermalization (FSP)

Fock Space Prethermalization (FSP) denotes a nonequilibrium regime in which a many-body system develops a long-lived quasi-steady state because the wave function spreads only within a restricted, highly structured region of Fock space, while full ergodic exploration of the Hilbert space is delayed. In this sense, observables can appear nearly stationary even though the state remains far from the uniform Fock-space distribution expected in a fully thermalizing system. In the recent literature, FSP has been formulated in several complementary ways: as constrained spreading near a chaos–integrable crossover in the complex SYK model, as dephasing at fixed mode occupations in an integrable Luttinger liquid, and as sparse-subnetwork dynamics in disorder-free Floquet circuits with approximately conserved domain-wall number [2303.16019] [1212.4645] [2510.24059].

## 1. Definition and operational meaning

In the SYK\(_4\)+SYK\(_2\) setting, prethermalization is defined operationally by an intermediate-time plateau of the infinite-temperature on-site density–density correlator,
\[
\mathcal{C}_i(t)=\left\langle \big| \langle n_i(t)n_i(0)\rangle_\infty - n_{\rm fill}^2 \big| \right\rangle_{\rm dis},
\]
whose value is close to the saturation value of the integrable SYK\(_2\) limit, followed only at much later times by decay toward the small ergodic SYK\(_4\) value [2303.16019]. In that formulation, the hallmark of FSP is not merely slow relaxation, but the coexistence of quasi-stationary observables with a wave function that remains strongly nonuniform in Fock space.

A central consequence is that FSP is not restricted to the conventional setting of nearly integrable local lattices. In the complex SYK model, it arises in an all-to-all quantum dot with no spatial locality and no appeal to quasi-local conserved quantities; the controlling structure is instead the geometry and connectivity of Fock space itself [2303.16019]. This directly separates FSP from the standard near-integrable prethermal scenario in which relaxation is delayed by approximate conservation laws in real space.

A broader operational picture emerges from other works. In the Luttinger-liquid description of a rapidly split one-dimensional Bose gas, prethermalization occurs because the system dephases at approximately fixed mode occupations \(n_k\), so the long-time observables are governed by a diagonal ensemble rather than by true redistribution of occupation numbers [1212.4645]. In the disorder-free Floquet Ising experiment on a superconducting processor, FSP is defined explicitly as a mechanism that “divides the Fock-space network into linearly many sparse sub-networks,” thereby suppressing heating and prolonging the thermalization timescale even at high energy density [2510.24059]. Across these formulations, the shared content is quasi-equilibration inside a restricted Fock-space manifold.

## 2. Canonical finite-size realization in the complex SYK model

The clearest analytic and numerical realization of finite-size FSP is the complex-fermion SYK quantum dot with both quartic and quadratic all-to-all random couplings,
\[
\mathcal{H}=\lambda\,\mathcal{H}_{\text{SYK}_4}+(1-\lambda)\,\mathcal{H}_{\text{SYK}_2}, \qquad \lambda\in[0,1],
\]
with
\[
\mathcal{H}_{\text{SYK}_4}=\sum_{i,j,k,l=1}^{N} J_{ijkl}\,c_i^\dagger c_j^\dagger c_k c_l,
\qquad
\mathcal{H}_{\text{SYK}_2}=\sum_{i,j=1}^N t_{ij}\,c_i^\dagger c_j .
\]
Here the random couplings are Gaussian with zero mean and variances
\[
\langle t_{ij}^2 \rangle_{\rm dis}=\frac{J^2}{64N},
\qquad
\langle J_{ijkl}^2 \rangle_{\rm dis}=\frac{J^2}{2N^3},
\]
which ensure an extensive many-body bandwidth [2303.16019].

The two limiting theories define the relevant endpoints of the crossover. At \(\lambda=1\), the pure SYK\(_4\) model has Wigner–Dyson statistics, ETH eigenstates, and fully ergodic infinite-temperature dynamics. At \(\lambda=0\), the pure SYK\(_2\) model is integrable, has Poisson level statistics, and its eigenstates are Slater determinants of single-particle orbitals [2303.16019]. Between these limits, a finite-size crossover scale
\[
\lambda_c(N)\propto N^{-5/2}
\]
(with logarithmic corrections \(\ln N\)) separates regimes that are dynamically close to Fock-space localization from those that are extended and RMT-like in the SYK\(_2\) basis [2303.16019].

The dynamics is studied at infinite temperature and half filling,
\[
n_{\rm fill}=\frac12,
\]
using the correlator \(\mathcal C_i(t)\) above, with time evolution computed by a Chebyshev expansion of \(e^{-i\mathcal H t}\) and stochastic typicality, and long-time limits obtained by exact diagonalization [2303.16019]. Operationally, the thermalization time \(t_{\rm th}\) is defined by a threshold such as
\[
\mathcal C_i(t_{\rm th}) = 1.5\,\mathcal C_i(\infty).
\]
For small but nonzero \(\lambda\) on the ergodic side \(\lambda>\lambda_c(N)\), \(\mathcal C_i(t)\) exhibits a pronounced prethermal plateau before eventual decay to the ergodic saturation value [2303.16019].

## 3. Fock-space geometry, distance, and constrained spreading

The Fock-space formulation is essential. In the SYK analysis, the basis states are occupation-number configurations
\[
|\psi_a\rangle=|n_1^a,n_2^a,\dots,n_N^a\rangle,\qquad n_i^a\in\{0,1\},
\]
restricted to half filling. The natural metric is the Fock-space or Hamming distance
\[
d_H(\psi_a,\psi_b)=\frac12\sum_{i=1}^{N}|n_i^a-n_i^b|,
\]
with maximal distance \(d_H^{\max}=N/2\) [2303.16019].

For an initial Fock state \(|\psi_0\rangle\), the probability to be at distance \(d\) at time \(t\) is
\[
\mathcal P_d(t)=\sum_{\psi_a:\, d_H(\psi_a,\psi_0)=d} |\langle \psi_a|\psi(t)\rangle|^2.
\]
To remove trivial shell combinatorics, this is compared with the “thermal” distribution \(\mathcal P_d^{\rm th}\), corresponding to a wave function uniformly random over all basis states at fixed filling [2303.16019]. The ratio \(\mathcal P_d(t)/\mathcal P_d^{\rm th}\) then measures deviation from uniform Fock-space spreading.

The two limits are sharply distinct. In the integrable SYK\(_2\) limit, the infinite-time profile is strongly biased toward small \(d\) and follows a stretched exponential,
\[
\mathcal P_d(t\to\infty)\sim \exp\!\big[-\sqrt{d/\xi}\,\big]\;\mathcal P_d^{\rm th},
\]
with \(\xi\) a Fock-space correlation length [2303.16019]. In the chaotic SYK\(_4\) limit, by contrast,
\[
\mathcal P_d(t\to\infty)\approx \mathcal P_d^{\rm th},
\]
which is the ETH/RMT expectation for a structureless wave function in the occupation basis [2303.16019].

FSP is the intermediate regime in which the wave function has relaxed within a restricted Fock-space region but not yet over the full Fock graph. In the SYK crossover, intermediate-time profiles retain the stretched-exponential SYK\(_2\)-like shape inside a shell \(d\lesssim \xi(\lambda)\), while probabilities outside that shell remain suppressed relative to \(\mathcal P_d^{\rm th}\). Only for \(t\gtrsim t_{\rm th}(\lambda)\) does the profile flatten toward the ergodic distribution [2303.16019]. This motivates the formulation of FSP as a Fock-space constrained ergodic regime.

A closely related but dynamically opposite benchmark is provided by dual-unitary Floquet dynamics. In the self-dual kicked Ising model, the generalized inverse participation ratios
\[
I_q(t)=\sum_z |\langle z|\psi(t)\rangle|^{2q}
\]
approach their Haar-random values exponentially fast in time, and the overlap distribution converges to Porter–Thomas on a timescale independent of system size [2408.02732]. That result functions as an “anti-FSP” reference point: FSP corresponds precisely to the failure of such rapid, system-size-independent delocalization.

## 4. Timescales, ETH, and the finite-size character of FSP

The central scaling result in the SYK crossover is
\[
t_{\rm th}(\lambda)\propto 2^{a/\lambda^{2/5}},
\]
with \(a\approx 2.2\) from numerical fits [2303.16019]. The origin of this form is a matching argument between the crossover scale \(\lambda_c\propto N^{-5/2}\) and the Heisenberg time
\[
t_H\sim \binom{N}{N/2}\propto 2^{bN}.
\]
Substituting \(N\sim \lambda_c^{-2/5}\) gives the observed exponential dependence on \(1/\lambda^{2/5}\) [2303.16019]. The prethermal plateau therefore survives essentially up to Heisenberg time as \(\lambda\to \lambda_c^+\).

This is why the phenomenon is termed *finite-size prethermalization*. The plateau exists for finite \(N\), is bounded above by \(t_H(N)\), and disappears straightforwardly in the thermodynamic limit at fixed \(\lambda>0\), where the asymptotic state is ergodic [2303.16019]. The long-time Fock-space extension is also captured by the first moment
\[
\Delta x(t)=\sum_d \mathcal P_d(t)\,d-\frac{4}{N},
\]
whose plateau at large nonzero values near \(\lambda_c\) signals localized-like intermediate dynamics, followed by decay toward zero once full ergodicity is restored [2303.16019].

The ETH connection can be made more explicit through the inverse participation ratio in a different clean interacting fermion problem. There, exact diagonalization shows that eigenstate-to-eigenstate fluctuations of the momentum occupation \(\hat f_k\) are proportional to the IPR
\[
\chi_\alpha=\sum_i (p_i^\alpha)^2
\]
in the Fock basis of the noninteracting Hamiltonian, so that Fock-space delocalization is the microscopic mechanism for the onset of eigenstate thermalization [1007.5306]. In that language, FSP corresponds to an intermediate regime in which \(\chi\) is already decreasing, so local observables can become quasi-stationary, but the wave function is still far from the \(1/D_K\)-type scaling associated with a fully chaotic state [1007.5306].

A recurrent misconception is that FSP implies permanent nonergodicity. The SYK analysis shows the opposite for any fixed \(\lambda>0\): despite an exponentially long delay, the system ultimately thermalizes, and Fock-space observables approach their ergodic SYK\(_4\) values [2303.16019]. A second misconception is that FSP must rely on spatial locality. The SYK realization shows that it can instead arise from finite-size hybridization structure on Fock space itself [2303.16019].

## 5. Related paradigms: dephasing, hidden thermal structure, transport, and fragmentation

One broad paradigm realizes prethermalization through dephasing at fixed Fock-space weights. In the experiment on a coherently split one-dimensional Bose gas, the relative sector is described by a harmonic Luttinger liquid,
\[
\hat H=\sum_k \hbar c|k|\left(\hat b_k^\dagger \hat b_k+\frac12\right),
\]
so the mode occupations \(n_k=\hat b_k^\dagger \hat b_k\) are conserved. Rapid splitting prepares a nonthermal distribution \(n_k\sim 1/|k|\), and time evolution dephases phases between Fock components \(|\{n_k\}\rangle\) without redistributing the occupations. The resulting quasi-steady state is described by an effective temperature
\[
k_B T_{\rm eff}=\frac{\mu}{2}=\frac{g\rho}{2},
\]
yet it is not the true final thermal state [1212.4645]. In FSP language, this is prethermalization generated by phase scrambling at fixed Fock-space populations.

A complementary static line of work identifies a “hidden thermal structure” in Fock space. For an ideal Fermi gas, coarse-grained occupations
\[
N_m=\sum_{\nu\in \mathcal G_m} n_\nu
\]
maximize the combinatorial entropy
\[
S[\Lambda]=\ln W[\Lambda]
\]
under particle-number and energy constraints, producing the typical occupation fraction
\[
\frac{N_m^*}{G_m}=\frac{1}{e^{(\varepsilon_m-\mu)/T}+1}.
\]
This shows that an overwhelming number of Fock states share an observable-resolved Fermi–Dirac limit shape \(\Lambda^*\) [1808.00784]. A plausible implication is that some FSP plateaus can be interpreted as dynamics constrained to, or dephasing within, coarse-grained Fock-space manifolds that already possess thermal structure for appropriate observables.

In disordered systems, the language of probability transport and fragmentation sharpens the dynamical picture. For a disordered transverse-field Ising chain written as a tight-binding model on a Fock-space hypercube, the probability distribution
\[
P_{IJ}(t)=|\langle J|e^{-i\mathcal H t}|I\rangle|^2
\]
shows strongly inhomogeneous, multifractal intermediate-time structure even in the ergodic phase, while in the MBL phase the inhomogeneity persists indefinitely [2212.14333]. In a distinct disordered interacting-fermion quench problem, the accessible phase space is organized into “potential-energy shells,” and at strong disorder these shells decay into disconnected fragments; the fragment containing the initial state then controls long-lived non-ergodic relaxation [2510.19510]. These results suggest that FSP can arise not only from approximate integrability, but also from multifractal transport or shell fragmentation on the Fock-space graph.

## 6. Experimental realizations and diagnostics

The most direct explicit realization of FSP on hardware is the periodically driven Ising chain implemented on 72 superconducting qubits. There the Floquet unitary is engineered so that strong Ising interactions and small perturbations make the total domain-wall number \(w\) approximately conserved. The Fock-space network then separates into \(O(L)\) sparse sub-networks labeled by even \(w\), and only spin flips that preserve local domain-wall structure are resonant at leading order [2510.24059]. This yields long-lived Fock-space confinement, suppressed heating, and FSP-based discrete time-crystalline order persisting over 120 cycles for generic initial Fock states [2510.24059].

The principal observables in that experiment are the Fock-space wave-packet distribution
\[
\Pi(d,nT)=\sum_{\boldsymbol s:\,D(\boldsymbol s,\boldsymbol s_0)=d}
|\langle \boldsymbol s|U_F^n|\boldsymbol s_0\rangle|^2,
\]
its mean position
\[
\langle x\rangle(t)=\sum_d d\,\Pi(d,t),
\]
its width
\[
\Delta x(t)=\sqrt{\sum_d (d-\langle x\rangle)^2\Pi(d,t)},
\]
and the domain-wall distribution
\[
\mathcal D(w,nT)=\sum_{\boldsymbol s:W(\boldsymbol s)=w}
|\langle \boldsymbol s|U_F^n|\boldsymbol s_0\rangle|^2.
\]
Finite-size scaling across \(L=24,40,56,72\) reveals size-independent crossover behavior in both domain-wall and Fock-space dynamics, and the eigenstates of the Floquet unitary cluster by average domain-wall number \(\langle w\rangle_n\) in the FSP regime [2510.24059].

A broader experimental toolkit for FSP-like questions is provided by the 24-qubit two-dimensional Bose–Hubbard/XY experiment that directly measures Fock-space wave-packet propagation. There the many-body problem is mapped to an Anderson-like model on a graph of \(\mathcal N=\binom{24}{12}=2{,}704{,}156\) Fock states, and the radial distribution
\[
\Pi(d,t)=\sum_{\mathbf s:\,D(\mathbf s,\mathbf s_0)=d}
|\langle \mathbf s|e^{-iH_F t/\hbar}|\mathbf s_0\rangle|^2
\]
is tracked experimentally, together with the displacement \(\mathrm x(t)\), width \(\Delta \mathrm x(t)\), and Bhattacharyya coefficient \(\mathcal B(t)\) relative to the ergodic radial distribution \(\Pi^{\rm Erg}(d)\) [2211.05803]. Although that work does not use the term FSP, it suggests a practical diagnostic program: FSP corresponds to long-lived confinement of \(\Pi(d,t)\) away from \(\Pi^{\rm Erg}(d)\), saturation of \(\mathrm x(t)\) and \(\Delta \mathrm x(t)\) at non-ergodic values, and eventual drift only on much longer timescales.

Taken together, these developments establish FSP as a unifying language for several nonequilibrium mechanisms. In one class, the wave function is temporarily trapped in a stretched-exponential or multifractal region of Fock space near a localization–delocalization crossover; in another, it dephases at fixed mode occupations; in another still, Floquet dynamics creates sparse sub-networks labeled by approximately conserved domain-wall number. The common structure is a separation between local or coarse-grained equilibration and global Fock-space ergodization. Whether the resulting plateau is transient, exponentially long, or asymptotically stable depends on the underlying mechanism: finite-size hybridization in the SYK crossover, exact or approximate integrability in harmonic mode problems, or fragmentation-like constraints in driven and disordered systems [2303.16019]

Source: https://www.emergentmind.com/topics/fock-space-prethermalization-fsp