REVIEW 4 minor 3 cited by
Expansion creates spin-glass order in finite-connectivity models: a rigorous and intuitive approach from the theory of LDPC codes
T0 review · 0 major / 4 minor · reviewed 2026-08-06 · deepseek-v4-flash
Pith's one-line read Expansion in a code forces the low-temperature Gibbs state of a finite-connectivity spin model to shatter into exponentially many tiny, incongruent components.
desk verdict Rigorous first proof of finite-temperature spin-glass order on finite-connectivity graphs via code expansion; solid theorem with clearly disclosed scope limits. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The load-bearing object is code expansion, Definition III.1: $H\in\mathbb{F}_2^{m\times n}$ is $(\delta,\gamma)$-expanding if $|Hx|>\gamma|x|$ whenever $|x|<\delta n$, meaning every small error violates a proportional number of checks. Because the energy is $E(x)=|Hx|$, the triangle inequality yields Eq. (20), $E(y)-E(x)\ge\gamma|x\oplus y|-2E(x)$, which implies that any state with energy density below $\delta\gamma/2$ is surrounded by an extensive energy barrier. Full rank of $H$ supplies the exact partition function and the count of states at each energy; combining the Hamming-ball upper bound on cluster size (Lemma D.6) with this exact counting gives Eq. (D11), a lower bound $s_{\mathrm{conf}}\ge r+s(T)$ that is positive for $\gamma>\gamma^*(r)$ at small $T$. That positivity is what proves shattering and incongruence.
What would settle it
Run the exact low-temperature calculation for a family of Gallager codes that provably satisfy rank $H_n=m$ and $\gamma>\gamma^*(r)$, for instance the $(5,6)$-ensemble used in Figure 8, and compute the configurational entropy of the explicit Gibbs decomposition from the closed-form partition function; if $s_{\mathrm{conf}}(T)$ fails to exceed $r$ at arbitrarily low temperature as the lower bound predicts, the theorem is false.
Extended reading notes
Core claim
The central claim is Theorem IV.2: for any family of LDPC codes $H_n$ with rate $r$ that is $(\delta,\gamma)$-expanding with $\gamma>\gamma^*(r)$ and has rank $H_n=m$, the associated classical spin model—energy equal to the number of violated parity checks—realizes spin-glass order at sufficiently low temperature. Spin-glass order is defined concretely as shattering: no Gibbs component carries more than an exponentially small fraction of the Gibbs weight, together with incongruence: the configurational entropy density $s_{\mathrm{conf}}(T)$ strictly exceeds its zero-temperature value $s_{\mathrm{conf}}(0)=r$. The proof constructs an explicit Gibbs-state decomposition, first showing that the energy landscape below a cutoff shatters into clusters separated by extensive Hamming distance and extensive energy barriers, and then converting those barriers into free-energy bottlenecks. The no-redundancy condition gives the exact partition function $Z=2^{n-m}(1+e^{-\beta})^m$, which makes the counting argument precise enough to show $s_{\mathrm{conf}}(T)>r$ for small $T$. Rigorous instantiations are Gallager-type diluted ferromagnetic $p$-spin models for sufficiently large but finite $p$.
Load-bearing premise
The load-bearing premise is that the code satisfies both a sufficiently strong expansion condition ($\gamma>\gamma^*(r)$, meaning every small error creates a large enough number of violated constraints) and full rank of $H_n$ (no redundant constraints); if either fails, the counting estimates that force shattering and incongruence collapse.
Editorial extensions
If this is right
- Rigorous low-temperature spin-glass order for diluted ferromagnetic $p$-spin models on random regular graphs, for sufficiently large but finite $p$.
- The argument is independent of local tree-likeness, so the same conclusion applies in principle to loop-rich expander geometries such as hyperbolic tessellations, where cavity methods become uncontrolled.
- Every Gibbs component is protected by an extensive free-energy barrier, so any local, detailed-balance stochastic dynamics is a passive memory: initialized in a typical low-temperature component, it remains there for a time growing exponentially with system size.
- The proof establishes a nontrivial overlap distribution consistent with one-step replica symmetry breaking, but it does not prove that the cavity or RSB solution is exact, and it does not rigorously establish the intermediate weak-ergodicity-breaking phase seen numerically.
- The exact partition function is analytic at all nonzero temperatures, so these models provide examples where a trivial partition function coexists with a nontrivial low-temperature Gibbs-state decomposition.
Reading between the lines
- Beyond the paper: improved coding-theory lower bounds on $\gamma$ for Sipser-Spielman and hyperbolic Tanner-Ising codes would immediately convert the numerical evidence in Section V into rigorous examples, since the remaining obstruction is the expansion bound rather than the physics.
- Beyond the paper: the no-redundancy assumption is likely replaceable by subextensive redundancies; the partition function then admits a lower bound via a Kramers-Wannier-type duality, which should preserve the entropy inequality.
- Beyond the paper: a plausible testable extension is that the intermediate phase $T_G<T<T_{\mathrm{mem}}$ has exponentially many components but one dominant component, with only weak dynamical signatures, mirroring the Bethe-lattice Ising model; the present theorem proves only the low-temperature shattered phase.
- Beyond the paper: a weaker 'linear confinement' condition with $\delta(n)$ growing only logarithmically may suffice, which would extend the result to codes with subextensive confinement and hence to a broader family of locally constrained spin models.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper proves that families of classical spin models associated with sufficiently expanding, redundancy-free LDPC codes display spin-glass order at low temperatures, in the precise sense of a Gibbs-state decomposition into exponentially many components separated by extensive free-energy barriers, with no component carrying more than an exponentially small fraction of the weight and with most components not containing ground states. The proof proceeds by first showing that code expansion creates extensive energy barriers around all low-energy-density states, then, under strong expansion and full rank of the parity-check matrix, proving shattering and incongruence of the low-energy landscape by a counting argument, and finally lifting these properties to the Gibbs state using the exact partition function and microcanonical concentration. The rigorous results apply to Gallager-type codes (diluted p-spin models) with sufficiently large bit degree. The paper also presents Monte Carlo and cavity-method studies of two Sipser-Spielman-type Tanner-Ising models on high-girth random regular graphs and hyperbolic tessellations, providing evidence for two separate transitions associated with weak ergodicity breaking and spin-glass order, while explicitly noting that these numerical models are not covered by the rigorous theorem.
Significance. If correct, this is the first rigorous demonstration of finite-temperature spin-glass order on closed finite-degree interaction graphs, a problem that has remained open beyond mean-field models. The strength of the paper is that the central argument is elementary and self-contained: the proof in Appendix D states precise assumptions, gives explicit counting bounds, and derives the exact partition function rather than assuming it. The paper is also honest about its scope: the numerical Tanner-Ising models do not satisfy the available rigorous expansion bounds, and the authors say so explicitly. The proof is not machine-checked, but the counting and bottleneck steps are checkable by hand. The work connects coding theory and statistical mechanics in a way that is likely to be influential and offers a concrete path toward rigorous results on loopy non-Euclidean graphs if sharper expansion bounds become available.
minor comments (4)
- [Sec. V A / Table I] The construction of the numerical model says that the authors add all seven nonzero linear combinations of the three local Hamming checks, a set that is linearly dependent, yet Table I states that the parity-check matrix H_G has no redundancies. Please clarify that the no-redundancy statement refers to the minimal three-check representation of the global code, while the simulated Hamiltonian is the symmetrized (locally redundant) form.
- [App. D, Theorem D.8 and Eq. (D11)] The theorem states that γ* depends only on the rate r, but the proof does not give an explicit threshold. From the small-ε expansion of the bound r + (1−r)h(ε/(1−r)) − h(2ε/γ), any fixed γ > 2 suffices to make the bound strictly exceed r for sufficiently small ε; stating this explicitly would make the theorem easier to verify.
- [Sec. II D, Eq. (10)] The replica identity displays no disorder average, although the text says that the overline denotes the disorder average. Adding the missing average to the equation would avoid ambiguity.
- [Sec. III C and Sec. IV] There are several typographical errors that should be corrected in a final version, including 'whcih', 'dam', 'the the system', and the reference to 'Theorem D.1' where 'Definition D.1' is meant.
Circularity Check
No significant circularity: the spin-glass conclusion is derived from explicit axioms (expansion, full rank) and an explicit counting/partition-function argument; the only self-citations are not load-bearing.
full rationale
The central derivation is self-contained. The load-bearing inputs are Definition III.1 ((δ,γ)-expansion) and the explicit assumption rank H = m (no redundancies); neither is derived from the target conclusion. The exact partition function in Lemma D.10, Z = 2^{n-m}(1+e^{-β})^m, is derived from full rank rather than imposed. The counting identity |Ω(ϵ)| = 2^{rn} B_m(ϵn) in Eq. (D9) follows from linearity and full rank, and together with the cluster-diameter bound (Lemma D.6) it yields the lower bound s_conf ≥ r + (1−r) log_2 Υ(ϵ/(1−r)) − log_2 Υ(2ϵ/γ) in Eq. (D11), which is strictly larger than r for γ > γ* and sufficiently small ϵ. The Gibbs-state step uses concentration on the microcanonical shell and the bottleneck condition; no fitted parameter is renamed as a prediction. The numerical T_G and T_mem in Sec. V are simulation outputs, and the authors explicitly state that the rigorous theorem does not cover those models because the available γ bounds are trivial (Sec. IV C 2). The self-citations that appear, notably Ref. [61] for context and Ref. [111] for a version of the bottleneck theorem, are not load-bearing: the bottleneck theorem is also credited to the external reference [63], and the central uniqueness/counting arguments do not reduce to the authors' prior work. Thus the derivation does not reduce to its own inputs.
Assumptions & free parameters
assumptions (6)
- domain assumption The energy of any configuration equals the Hamming weight of its syndrome, H(σ) = |Hx| (Eq. 14).
- domain assumption The parity-check matrix H has full rank (no redundancies) for every n in the family.
- domain assumption The code family is (δ,γ)-expanding with γ > γ*(r).
- standard math Gallager (w_bit, w_check) random code ensembles have full rank with high probability (Lemma C.1) and are unique-neighbor expanders with sufficiently large γ for large enough w_bit (Theorem C.3).
- standard math Hamming ball volume bounds (Ref. 110) and concentration of the energy under full rank via Hoeffding's inequality (Lemma D.10).
- standard math The classical bottleneck theorem for Markov chains (Theorem D.13, citing Refs. 63 and 111).
Cite this review
Pith. "Pith review of Expansion creates spin-glass order in finite-connectivity models: a rigorous and intuitive approach from the theory of LDPC codes." pith.science (2026). https://pith.science/paper/WBP2RZNY
@misc{pith2026250713342,
author = {Pith},
title = {Pith review of: Expansion creates spin-glass order in finite-connectivity models: a rigorous and intuitive approach from the theory of LDPC codes},
year = {2026},
howpublished = {\url{https://pith.science/paper/WBP2RZNY}},
note = {Machine review of arXiv:2507.13342}
}
read the original abstract
Complex free-energy landscapes with many local minima separated by large barriers are believed to underlie glassy behavior across diverse physical systems. This is the heuristic picture associated with replica symmetry breaking (RSB) in spin glasses, but RSB has only been rigorously verified for certain mean-field models with all-to-all connectivity. In this work, we give a rigorous proof of finite temperature spin glass order for a family of models with local interactions on finite-connectivity, non-Euclidean expander graphs. To this end, we bypass the RSB formalism entirely, and instead exploit the mathematical equivalence of such models to certain low-density parity check (LDPC) codes. We use code expansion, a property of LDPC codes which guarantees extensive energy barriers around ground states. Together with mild additional assumptions, this allows us to construct an explicit decomposition of the low-temperature Gibbs state into disjoint components, each hosting an asymptotically long-lived state associated with a local minimum of the landscape. Each component carries at most an exponentially small fraction of the total weight, and almost all components do not contain ground states -- which we take together to define spin-glass order. The proof is elementary, and treats various expanding graph topologies on the same footing, including those with short loops where existing approaches such as the cavity method fail. Our results apply rigorously to diluted p-spin glasses for sufficiently large p, and while unproven, we also expect our assumptions to hold in a broader family of codes. Motivated by this, we numerically study two simple models, on random regular graphs and a regular tesselation of hyperbolic space. We show that both models undergo two transitions as a function of temperature, corresponding to the onset of weak ergodicity breaking and spin glass order, respectively.
Figures
Figures from the paper (24 more)
Forward citations
Cited by 3 Pith papers
-
Topological states and flat bands in exactly solvable decorated Cayley trees
Flat bands on decorated Cayley trees map exactly onto topological edge states of 1D SSH chains, and persist on infinite Bethe lattices by a covering construction.
-
Lifting Lifted Product Codes
Group-extension lifts systematically enlarge any LP code, transfer logical gadgets via chain maps (often with less surgery overhead), improve some code parameters, and give candidate thermodynamic families with cohere...
-
Superconductivity in hyperbolic spaces: Regular hyperbolic lattices and Ginzburg-Landau theory
On hyperbolic lattices, superconductivity can appear only at the boundary above the bulk critical temperature, and rough boundaries with dangling bonds can raise Tc severalfold.
Reference graph
Works this paper leans on
-
[1]
Choose a vertexvrandomly among those vertices inGwith degrees min
-
[2]
Table of expander graphs referenced in this work
Choose a second vertexuamong those vertices inG with degrees min and graph distance dist(u, v)> g 36 HGRRG (s= 7) Hyperbolic{3,7}LPS girthN v Ne λ2 systoleN v Ne λ2 p q N v Ne λ2 3 38 133 3.92 4 24 84 2.65 7 11 1320 5280 5.12 4 218 763 4.67 6 72 252 4.76 7 13 2184 8736 5.00 5 1298 4543 4.80 7 156 546 5.30 7 17 4896 19584 5.11 6 7778 27223 4.87 7 192 672 5...
-
[3]
Constructing expander graphs There are many known constructions of expander graphs, both random and explicit. We will here review three of them which are relevant for this work which are (1) a family of random-regular graphs with guaranteed girth (that is the size of the shortest loop) (2) regular tesselations of the hyperbolic plane, and (3) a symmet- ri...
-
[4]
Add the edge (u, v) toG The procedure described above can get stuck in princi- ple, for example if at any point all vertices with minimum degree left in the graph are close. However, as shown in Ref. 94, the procedure terminates with high probability. Of course, the family of RRGs constructed here is not exactly the same as the one considered by Friedman ...
-
[5]
Clustering and Configurational Entropy of the Energy Landscape We start by characterizing the complexity of the low energy landscape. We do this by analyzing the structure of the configuration space below a given cutoff energy densityϵ: Ω(ϵ) ={x;|Hx|< ϵn}.(D1) We will show that this set, for certain models at suffi- ciently lowϵ, can be decomposed into ex...
-
[6]
Spin Glass Order from Expansion In the previous section, we showed that certain non- redundant expander codes display a complex energy landscape, i.e., the set of configurations below a given en- ergy density cutoff has a complex cluster decomposition with shattering and incongruence. In this subsection, we will be interested in properties of theGibbs dis...
-
[7]
EachΩ j is surrounded by a bottleneck: there exists η >0, and a function∆(n)with∆(n)− − − − → n→∞ 0 such that pG[∂ηΩj,n] pG[Ωj,n] ≤∆(n) (D29) where theη-boundary of a subsetΩis defined as ∂ηΩ≡ {x∈Ω c; dist(x,Ω)≤ηn}
-
[8]
The setΛcontains only a vanishing fraction of the weight pG,n(Λn)− − − − → n→∞ 0.(D30) We callΛthe “junk” set. Gibbs state components are then defined as the (normal- ized) restrictions ofp G toΩ j: p(Ωj,n) G,n = 1Ωj,n ·p G,n pG,n(Ωj,n) (D31) where1 Ω is the indicator function of the setΩ. The Gibbs distribution can now be written the convex sum of its co...
Show all 24 references
-
[9]
AllΩ j are typical subsets (Theorem D.11), 2.diam[Ω j ∩Ξ ω(β)]≤2n(⟨ε⟩ β +ω)/γwithΞ ω(β)the microcanonical shell [Eq. (D13)]
-
[10]
Energy density Configuration Space T typical cluster atypical cluster 0 FIG
The{Ω j}andΛ≡F n 2 /(⊎M j=1Ωj)define a decom- position of the Gibbs state in the sense of Theo- rem D.14. Energy density Configuration Space T typical cluster atypical cluster 0 FIG. 22. Sketch of the various energy scales and distances involved in the Gibbs state decompositio...
-
[11]
Definition E.1.We call a binary matrixH∈F m×n 2 (w, b)-sparseif each row has Hamming weight≤wand each column has Hamming weight≤b
General results We will always assume that the number of checks is proportional to the number of bitsm= Θ(n). Definition E.1.We call a binary matrixH∈F m×n 2 (w, b)-sparseif each row has Hamming weight≤wand each column has Hamming weight≤b. Theorem E.2.Consider an expander cod...
-
[12]
Gallager codes We now show that random LDPC codes satisfy the assumptions of Theorem E.2 and Theorem E.4. By??C.1??C.3 choosing the left-degreeℓto be odd and large enough we obtain a Gallager code withH G be- ing full-rank, (r, ℓ)-sparse and (δ >0, γ >2)-expanding and thus sat...
-
[13]
symmetrized
No redundancies Given a parity check matrixH∈F m×n 2 with full rank m, the energy of a statex∈F n 2 is given byE:=|s|, whereHx=s∈F m 2 is called the syndrome ofx. Sinces i ∈ {0,1}and hence|s|= P i si, we can write the Gibbs distribution as pG ∝e −βE(σ) =e −β|s| = Y i e−βsi .(G...
-
[14]
symmetrized
Local redundancies We now discuss the case of how to sample from the Gibbs distribution of the “symmetrized” models con- sidered in the main text. The symmetrized model is based on the Tanner code on a graphG,T G, H(sym) L , where given a local codeH L ∈F m0 2 without redundan...
-
[15]
redundant
Constant number of redundancies Consider now the general case of a parity check matrix H∈F m×n 2 with rankm 0. As before, we want to sample from the Gibbs distributionp G ∝e −β|s|. In the case of no redundancies, we can simply choose a randomsby sampling the syndrome component...
-
[16]
We can represent this conveniently as a tensor network (see also Refs
T ensor network formulation The unnormalized Gibbs state,e −βE(σ) , of a Tanner code factorizes into Q v e−βEv(σv) whereσ v denotes the spin configuration on the edges adjacent tov. We can represent this conveniently as a tensor network (see also Refs. [112, 113] for similar t...
-
[17]
inconsistent
Recursion relation for general boundary conditions Consider forming a rooted tree by joining together s−1 branches, where theith branch has boundary con- figurationσ (i) ∂ , partition functionZ ∂(σ(i) ∂ ), and condi- tional magnetizationm i on the root spin. Letm= (m2, ..., ms...
-
[18]
ferromagnetic
Memory transition An exceptional case where the fixed point distribution isanalytically solvable isα=∞, where the only al- lowed boundary configurations are those consistent with a global codeword. As a result, the distribution at depth rconcentrates onto a pair of delta funct...
-
[19]
ForT < Tmem,h opp bulk(T)>0: up to this point, a strong enough negative boundary field can overwhelm the positive boundary field
Turning on a bulk field, we leth opp bulk(T) denote the magnitude of the largest bulk field at which negatively magnetized and positively magnetized stable fixed points coexist. ForT < Tmem,h opp bulk(T)>0: up to this point, a strong enough negative boundary field can overwhel...
-
[20]
annealed
Spin glass transitions As we have seen, the tensor network formulation natu- rally allows us to calculate bulk observables conditioned on a particular boundary configurationσ ∂. In partic- ular, lettingσdenote the configuration of spins on the interior, we can consider the mag...
-
[21]
Sample the spin configurationσ 2, ..., σs with Boltz- mann weightT v(σ, σ2, ..., σs)
-
[22]
For eachσ i, independently samplem i from popu- lation ˜Q(r) σi
-
[23]
spin glass
Set thejth element of the population ˜Q(r+1) σ equal toF(m 2, ..., ms). This method was used in Ref. [52] to analyze theq- state Potts model with spin degrees of freedom at the vertices of a tree. For those models, theσ i in step (1) can be sampled independently, a simplificat...
-
[24]
We measure the entropy density in bits, so thats= β(ε−f)/ln(2) whereε, fare the energy density and free- energy density, respectively, at a given temperature
Configurational entropy As discussed in the main text, the configurational entropy is the difference of two terms,s β ands α=1. We measure the entropy density in bits, so thats= β(ε−f)/ln(2) whereε, fare the energy density and free- energy density, respectively, at a given tem...
Reviewed August 6, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.