---
title: Diluted p-Spin Glasses Overview
url: https://www.emergentmind.com/topics/diluted-p-spin-glasses
type: topic
---

# Diluted p-Spin Glasses Overview

Searching arXiv for relevant papers on diluted p-spin glasses and closely related results.
arxiv_search.query({"search_query":"all:\"diluted p-spin\" OR ti:\"diluted p-spin\"","start":0,"max_results":10})
arxiv_search.search({"query":"diluted p-spin glass", "max_results": 10})
Diluted \(p\)-spin glasses are spin-glass models in which the interaction set is not the complete graph of the fully connected mean-field theory, but a random hypergraph or a power-law diluted long-range network. In the cited literature, this class includes Poisson random \(p\)-uniform hypergraphs with Hamiltonian
\[
H_N(\sigma) = \sum_{(i_1\ldots i_p)\in E_N} J_{i_1\ldots i_p}\,\sigma_{i_1}\cdots \sigma_{i_p} + h \sum_{i=1}^N \sigma_i,
\]
general diluted models with random interaction functions \(\theta_k\) on uniformly chosen \(p\)-tuples, mixed diluted \(p\)-spin systems on random hypergraphs \(E=\bigcup_{p=2}^\Delta E_p\), and one-dimensional long-range proxies with fixed average coordination and power-law dilution [1308.1944] [2410.15599] [2404.09409] [2604.15433]. The subject combines replica and cavity methods, asymptotic Gibbs-measure theory, recursive distributional equations, finite-size scaling, and spectral methods for disorder perturbations; its main themes are replica symmetry, one-step and full replica-symmetry breaking, and the existence or absence of thermodynamic glass transitions.

## 1. Model classes and Hamiltonian structure

A standard diluted \(p\)-spin model is defined on \(N\) Ising spins \(\sigma_i\in\{-1,+1\}\). In one formulation, a \(\mathrm{Poisson}(\alpha N)\) number of \(p\)-tuples is sampled i.i.d. uniformly from \(\{1,\dots,N\}^p\), each sampled tuple receives an independent standard Gaussian coupling, and the Hamiltonian is
\[
H_N(\sigma) = \sum_{(i_1\ldots i_p)\in E_N} J_{i_1\ldots i_p}\,\sigma_{i_1}\cdots \sigma_{i_p} + h \sum_{i=1}^N \sigma_i.
\]
The free energy density is
\[
F_N=\frac1N E\log \sum_{\sigma\in\{\pm1\}^N}\exp H_N(\sigma),
\]
and \(F=\lim_{N\to\infty}F_N\) exists [1308.1944].

A broader formulation places the model on a random \(p\)-uniform hypergraph with sparsity parameter \(\alpha\). One draws \(\pi\sim\mathrm{Pois}(\alpha N)\), chooses \(I(k)\subset [N]\) uniformly of size \(p\), assigns i.i.d. interaction functions \(\theta_k\), and writes
\[
H_N(\sigma)=\sum_{k\le\pi(\alpha N)} \theta_k(\sigma_{I(k,1)},\dots,\sigma_{I(k,p)}) + \sum_{i=1}^N \psi(\sigma_i).
\]
This formulation accommodates Ising \(p\)-spin, Potts, XY, and continuous hardcore examples within a common cavity framework [2410.15599].

A mixed diluted \(p\)-spin model allows hyperedges of several cardinalities. For each \(p\in\{2,3,\dots,\Delta\}\), every \(p\)-subset \(e\subseteq[N]\) is included independently with probability
\[
\Pr(e\in E_p)=\frac{\alpha_p\,N}{\binom Np},
\]
and the Hamiltonian is
\[
H_J(\sigma)=\sum_{p=2}^\Delta \sum_{e\in E_p}\rho(J_e)\,\sigma_e,\qquad \sigma_e=\prod_{i\in e}\sigma_i,
\]
with \(\rho\) an odd measurable function and \(J_e\) independent standard Gaussians [2404.09409].

A different dilution mechanism appears in the balanced \(M=4\), \(p=4\) long-range proxy. There are \(M=4\) Ising spins \(\{S_i,T_i,U_i,V_i\}\) on each of \(L\) sites of a ring, and the Hamiltonian is
\[
H = - \sum_{i<j}\sum_{(a<b),(c<d)} J_{ij}^{(ab,cd)}\,S_i^a S_i^b S_j^c S_j^d.
\]
In the diluted fixed-coordination-number \(z\) version, each rung is connected on average to \(z\) others, with
\[
\mathrm{Prob}(J_{ij}\neq 0)\propto 1-e^{-A/r_{ij}^{2\sigma}},
\]
where \(r_{ij}=(L/\pi)\sin[\pi|i-j|/L]\). In both the fully connected and diluted long-range constructions, the limit \(\sigma\to 0\) recovers an infinite-range model, while larger \(\sigma\) moves the system toward a short-range-like regime [2604.15433].

## 2. Overlaps, order parameters, and mean-field organization

The basic order parameter is the overlap. For two configurations \(\sigma,\tau\in\{\pm1\}^N\),
\[
R_N(\sigma,\tau)=\frac1N\sum_{i=1}^N \sigma_i\tau_i.
\]
In replicated formulations of \(p\)-spin glasses, one also writes
\[
q_{ab}=\frac1L\sum_i S_i^a S_i^b.
\]
In the asymptotic 1-RSB setting, the overlap takes exactly two values \(q_*<q^*\): the self-overlap is \(R^{\ell,\ell}=q^*\), whereas for \(\ell\neq \ell'\) one has \(R^{\ell,\ell'}=q_*\) [1308.1944].

Mean-field \(p\)-spin glasses with \(p\ge 3\) exhibit a dynamical transition at \(T_d\), described as mode-coupling-like nonergodicity, and a static Kauzmann temperature \(T_K\) at which the configurational entropy vanishes. Their mean-field low-temperature phase is one-step replica-symmetry breaking: the order-parameter overlaps jump discontinuously at \(T_K\). At still lower temperature, a Gardner transition induces a full replica-symmetry-breaking hierarchy inside the 1-RSB phase [2604.15433].

To distinguish weakly first-order 1-RSB behavior from a continuous FRSB transition, the replicated Gibbs free-energy expansion introduces
\[
\lambda\equiv \omega_2/\omega_1,
\]
where \(\omega_1\) and \(\omega_2\) are linear combinations of third-order cumulants of overlap fluctuations. For the balanced \(M=4\), \(p=4\) model, mean field gives \(\lambda=2\). In that setting, \(\lambda>1\) is associated with a discontinuous 1-RSB transition, whereas \(\lambda<1\) yields a continuous FRSB transition [2604.15433].

The overlap distribution \(P(q)\) is correspondingly diagnostic. A two-peak structure is the characteristic mean-field signature of 1-RSB, while a broad continuous \(P(q)\) extending from \(q\approx 0\) to a maximum is characteristic of a direct transition to an FRSB-like phase in the numerical long-range proxy [2604.15433]. A frequent misconception is that dilution by itself enforces replica symmetry; the cited work does not support that conclusion. Rather, replica symmetry, 1-RSB, and FRSB appear as distinct regimes whose realization depends on temperature, sparsity, and the geometry of the underlying interaction structure.

## 3. Asymptotic Gibbs measures and the rigorous 1-RSB description

The rigorous theory of diluted \(p\)-spin glasses studies infinite-volume Gibbs measures via exchangeability and cavity identities. In Panchenko’s framework, one draws \(w\sim\mathrm{Unif}[0,1]\), constructs a random probability measure \(G_w\) on
\[
H=\{u:[0,1]\to[-1,1]\mid \|u\|_2\le 1\},
\]
and samples i.i.d. \(u^1,u^2,\dots\sim G_w\). The support of \(G_w\) then yields pure-state functions with weights \(V_\alpha\), and independent uniform variables are inserted into these functions to generate spin variables in each pure state [1308.1944].

Under the 1-RSB hypothesis, the Mézard–Parisi ansatz predicts two key structures. First, the atomic masses \((V_\alpha)_{\alpha\ge 1}\) are independent of the array of spin magnetizations and have the Poisson–Dirichlet law \(\mathrm{PD}(m)\), \(0<m<1\). Second, there is a measurable function \(f:[0,1]^3\to[-1,1]\) and an auxiliary uniform random variable \(\theta\) such that, conditionally on \((V_\alpha)\) and independent uniforms \((\theta_\alpha)\), \((\phi_i)\), the magnetization of spin \(i\) in pure state \(\alpha\) is
\[
f(\theta,\theta_\alpha,\phi_i).
\]
This gives the factorized 1-RSB description of pure-state magnetizations in the diluted setting [1308.1944].

The associated cavity formalism produces a distributional fixed-point equation for cavity fields. In the 1-RSB setting, a single cavity field \(x\) satisfies in law
\[
x \;\overset{d}{=}\; f(x_1,\dots,x_K;J),
\]
with \(K\sim\mathrm{Poisson}(\alpha_p)\), \(J\sim N(0,1)\), and
\[
f(x_1,\dots,x_K;J)=\mathrm{atanh}\!\Bigl[\tanh(\beta J)\prod_{\ell=1}^{p-1}\tanh(x_\ell)\Bigr]+h.
\]
This is the usual population-dynamics fixed-point equation in the 1-RSB cavity method [1308.1944].

The rigorous identification of the ansatz proceeds through cavity equations and hierarchical exchangeability. A hierarchical version of Aldous–Hoover shows that any such law factors through nested i.i.d. randomizations at each RSB level, while the cavity equations remove an extra Ghirlanda–Guerra coordinate and reduce the magnetization function from \(f:[0,1]^4\to[-1,1]\) to \(f:[0,1]^3\to[-1,1]\) [1308.1944].

The main theorem is conditional on being in the 1-RSB regime. If \(h=0\) and \(q_*\neq 0\), any 1-RSB asymptotic Gibbs measure satisfies the Mézard–Parisi ansatz. If \(h\neq 0\), then necessarily \(q_*>0\), and every 1-RSB Gibbs measure satisfies the full 1-RSB ansatz. The exceptional case \(h=0\) and \(q_*=0\) is only partially characterized: one proves a symmetry property of the conditional law of the spin-function, but not the full factorization formula [1308.1944]. This establishes the rigorous structure of diluted 1-RSB Gibbs measures without asserting that every diluted \(p\)-spin model actually realizes a 1-RSB phase.

## 4. Replica symmetry and recursive distributional equations

Replica symmetry in general diluted spin glasses can also be formulated as a fixed-point problem on laws of single-site log-densities. Let \(X\in C(\Sigma)\) represent the log-density of the single-site marginal. If \(X_1,X_2,\dots\) are i.i.d. with law \(\lambda\), \(K\sim\mathrm{Pois}(\alpha p)\), and \(\theta_1,\dots,\theta_K\) are i.i.d. interaction functions, then the cavity operator is
\[
T_{\nu,\beta,K}(X_1,\dots,X_{(p-1)K})(t)
=\frac1\beta \ln
\frac{\int e^{\beta\sum_{k=1}^K \theta_k(\sigma_{k,1},\dots,\sigma_{k,p-1},t)+\beta\psi(t)}
\prod_{j=1}^{(p-1)K} e^{\beta X_j(\sigma_j)}\,\nu(d\sigma)}
{\int e^{\beta\sum_{k=1}^K \theta_k(\sigma_{k,1},\dots,\sigma_{k,p-1},s)+\beta\psi(s)}
\prod_{j=1}^{(p-1)K} e^{\beta X_j(\sigma_j)}\,\nu(d\sigma)\,\nu(ds)}.
\]
Replica symmetry means that the law \(\mu\) of \(X\) satisfies the recursive distributional equation
\[
\mu=\mathcal F_{\beta,\alpha}(\mu)
\]
obtained by averaging over \(K\) and the \(\theta\)’s [2410.15599].

The corresponding RS free-energy functional is
\[
\Phi_{RS}(\lambda)
=
E\Bigl[
\ln\int e^{\beta\sum_{k=1}^K \theta_k(\sigma_{k,1},\dots,\sigma_{k,p-1},\varepsilon)+\beta\psi(\varepsilon)}\,\nu(d\varepsilon)
\Bigr]
-\alpha(p-1)\,
E\bigl[\ln E_\beta(\theta,X_1,\dots,X_p)\bigr],
\]
where
\[
E_\beta(\theta,X_1,\dots,X_p)
=
\int e^{\beta\theta(\sigma_1,\dots,\sigma_p)}\prod_{i=1}^p e^{\beta X_i(\sigma_i)}\,\nu(d\sigma_i).
\]
Under the RS ansatz,
\[
f_{RS}(\beta,\alpha)=\min_{\lambda:\lambda=\mathcal F_{\beta,\alpha}(\lambda)} \Phi_{RS}(\lambda),
\]
and if the fixed point \(\lambda^*\) is unique then
\[
\lim_{N\to\infty}F_N(\beta)=\Phi_{RS}(\lambda^*).
\]
The same framework admits a zero-temperature operator \(T_\infty\) and a variational formula for the ground-state energy [2410.15599].

Two replica-symmetric regimes are isolated. The first is a high-temperature regime in which the cavity operator is contractive in \(W_1\); one sufficient condition is
\[
6\,\beta\,e^{4\beta\|\psi\|_\infty}\,E\|\theta\|_\infty\,e^{4\beta\|\theta\|_\infty}\cdot \alpha\,p(p-1)\le 1.
\]
The second is the subcritical regime
\[
\alpha\,p(p-1)\le 1,
\]
where the underlying Galton–Watson tree of branching mean \(\alpha p(p-1)\) is almost surely finite, giving uniqueness of the fixed point independently of temperature [2410.15599].

For the symmetric Ising \(p\)-spin case with \(\theta(\sigma)=J\sigma_1\cdots \sigma_p\), \(\psi=0\), and uniform \(\nu\) on \(\{\pm1\}\), the unique RS solution is \(\lambda^*=\delta_0\), and the explicit formulas are
\[
f_{RS}(\beta,\alpha)=\alpha\,E\log 2\cosh(\beta J)-\alpha p\log 2,
\qquad
e_0(\alpha,p)=\alpha\,E|J|.
\]
These results delineate where dilution is compatible with a complete RS description and where more complicated Gibbs structures must be expected [2410.15599].

## 5. Numerical phase structure in a power-law diluted long-range proxy

A detailed numerical study of phase transitions was carried out for the balanced \(M=4\), \(p=4\) spin-glass model used as a one-dimensional long-range proxy for finite-dimensional short-range \(p\)-spin glasses. The simulations employed Metropolis single-spin updates augmented by parallel tempering across \(N_T\) temperatures between \(T_{\min}\) and \(T_{\max}\). Four real replicas per disorder sample were simulated to compute overlap distributions and third-order cumulants, and equilibration was checked with the energy–link-overlap identity
\[
U=-(36/T)\,[N_b/L\,(1-q_l)]_{\rm av},
\]
where \(N_b=Lz/2\) in the diluted model [2604.15433].

The spin-glass susceptibility at zero momentum is
\[
\chi_{\rm SG}
=
\frac{1}{L M^4}\sum_{ij}\sum_{r=1}^{36}
\bigl[
\langle O_{r,i}O_{r,j}\rangle-\langle O_{r,i}\rangle\langle O_{r,j}\rangle
\bigr]^2_{\rm av},
\]
where \(O_{r,i}\) runs over the 6 pairwise spin products on each rung, giving 36 channels in total. The finite-size scaling forms are known exactly in the long-range mean-field regime \((\sigma<2/3)\),
\[
\chi_{\rm SG}\sim L^{1/3}\,\tilde\chi[(T-T_c)L^{1/3}],
\]
and in the non-mean-field regime \((\sigma>2/3)\),
\[
\chi_{\rm SG}\sim L^{2-\eta}\,\tilde\chi[(T-T_c)L^{1/\nu}],
\]
with crossings drifting according to
\[
T^*(L,2L)=T_c+A/L^{2/3}+\cdots .
\]
This supports extrapolation of \(T_c\) from size-pair crossings [2604.15433].

The extracted critical temperatures are as follows.

| Model and regime | \(\sigma\) | Critical observation |
|---|---:|---|
| Fully connected (variance-scaled) | \(0\) | \(T_c=0.320(7)\) |
| Fully connected (variance-scaled) | \(0.25\) | \(T_c=0.322(6)\) |
| Diluted \(z=6\) | \(0\) | \(T_c=0.701(4)\) |
| Diluted \(z=6\) | \(0.25\) | \(T_c=0.695(4)\) |
| Diluted \(z=6\) | \(0.55\) | \(T_c=0.673(4)\) |
| Diluted non-MF | \(0.85\) | no crossing of \(\chi_{\rm SG}/L^{2-\eta}\rightarrow\) no finite-\(T_c\) |

Below their respective \(T_c\), both the fully connected and diluted models exhibit broad, continuous \(P(q)\) distributions extending from \(q\approx 0\) to a maximum, rather than the two-peak structure expected for a 1-RSB transition. The estimated \(\lambda\)-parameter, obtained from three- and four-replica estimators, remains below unity near \(T_c\) for all accessible sizes and all studied \(\sigma\). The numerical interpretation is therefore a direct transition from the paramagnetic state to an FRSB-like phase with renormalized \(\lambda<1\), despite the mean-field value \(\lambda=2\) for the balanced model [2604.15433].

For \(\sigma=0.85\), which is described as roughly corresponding to a three-dimensional system, the results are stronger: \(P(q)\) shows no bimodality, \(\lambda\ll 1\) throughout, and there are no signs of either a 1-RSB transition or a continuous FRSB transition at finite temperature. The authors argue that strong finite-size effects and closely spaced transition temperatures remove the expected 1-RSB transition for the accessible system sizes [2604.15433]. Taken together with the rigorous 1-RSB theory, this separates two questions that are often conflated: how a 1-RSB Gibbs measure is structured if such a phase occurs, and whether a specific diluted model actually realizes that phase.

## 6. Disorder chaos, dimensional implications, and open problems

Disorder chaos in diluted \(p\)-spin glasses concerns the instability of equilibrium states under perturbations of the couplings. In the mixed diluted model, one compares the original disorder \(J\) with a perturbed copy \(J(t)\), either through an Ornstein–Uhlenbeck interpolation,
\[
J_e(t)=e^{-t}J_e+\sqrt{1-e^{-2t}}\,J'_e,
\]
or through independent replacement,
\[
J_e(t)=B_e J_e+(1-B_e)J'_e,\qquad B_e\sim\mathrm{Bernoulli}(e^{-t}),
\]
with \(J'_e\) i.i.d. standard Gaussians. Sampling \(\sigma\sim G_{J,\beta}\) and \(\tau\sim G_{J(t),\beta}\), one defines the site overlap
\[
R(\sigma,\tau)=\frac1N\sum_{i=1}^N \sigma_i\tau_i.
\]
If
\[
\lambda=\sum_{p=2}^\Delta p(p-1)\,\alpha_p>1,
\]
then for either perturbation, every \(\beta\in[0,\infty]\), and every \(t\ge 0\), there is a constant \(C=C(\alpha_2,\dots,\alpha_\Delta)\) such that
\[
R(\sigma,\tau)^2_t\le \frac{C}{N^{\,t/(t+2\ln\lambda)}}.
\]
For any fixed \(t>0\), the squared overlap therefore vanishes at a polynomial rate as \(N\to\infty\) [2404.09409].

The proof uses the Hermite spectral method. Writing the two-spin correlation functions \(\phi_{ij}(J)=\langle \sigma_i\sigma_j\rangle_{G_{J,\beta}}\) in the Hermite basis gives a semigroup identity under perturbation. A sign-flip lemma then eliminates low-frequency coefficients whenever a suitable sign vector can be chosen, and a local hypertree lemma shows that on tree-like balls \(B_r(i)\), coefficients with \(|E(n)|<\min\{d(i,j),r\}\) must vanish. Combining this with an exploration process on the sparse hypergraph yields the overlap bound above [2404.09409]. The result applies to any mixture, including odd \(p\), and removes the large-connectivity restriction of earlier work discussed in the same paper.

The dimensional implications of dilution remain contested. In the numerical long-range proxy, the mapping \(\sigma\approx 0.85\) to short-range \(d\approx 3\) is used to suggest the absence of any thermodynamic glass transition \(T_K>0\) in three dimensions; the data are described as consistent with \(T_K=0\), and the paper states that this supports the view that structural glasses do not undergo a finite-temperature Kauzmann transition in \(d\le 3\) [2604.15433]. This is a numerical interpretation of a proxy model rather than a theorem, but it directly addresses a central question in the theory of structural glasses.

Several open problems are explicit in the rigorous literature. These include extension from 1-RSB to full \(r\)-RSB, treatment of the random \(K\)-SAT model, identification of the exact 1-RSB region, and proving absence of condensation transitions beyond one step [1308.1944]. Another unresolved issue is how the rigorous RS regimes of subcritical or sufficiently high-temperature dilution interface with the non-RS phases observed numerically in long-range diluted proxies and with the universal disorder-chaos behavior on supercritical sparse hypergraphs. The current body of work shows that diluted \(p\)-spin glasses do not admit a single universal phase description: replica symmetry is exact in identifiable regimes, 1-RSB has a rigorous asymptotic structure when it occurs, FRSB-like behavior can dominate accessible numerical regimes, and sparse geometry produces strong sensitivity to perturbations of the disorder.

Source: https://www.emergentmind.com/topics/diluted-p-spin-glasses