Hardness of sampling for the anti-ferromagnetic Ising model on random graphs
Neng Huang, Will Perkins, Aaron Potechin
TL;DR
The paper demonstrates a fundamental gap between search and sampling for the anti-ferromagnetic Ising model on random graphs: although near-optimal search is tractable via mean-field-like methods, stable sampling of the Gibbs distribution becomes provably hard when the inverse temperature exceeds the condensation threshold, as $d$ grows. The authors extend the interpolation framework between the Sherrington–Kirkpatrick (SK) model and diluted sparse graphs to include average energy and overlap, not just free energy, and they employ a Poisson Stein-Chen approach to mirror Gaussian integration by parts in the sparse setting. They establish that the sparse model inherits SK-like overlap structure in the high-degree limit and that disorder chaos transfers from SK to the sparse model, yielding a constant Wasserstein-distance gap for any stable sampler at $eta>1$. This combination of energy/overlap interpolation and disorder chaos provides a robust mechanism to rule out stable sampling algorithms and highlights a sharp separation between optimization and sampling in dilute mean-field systems.
Abstract
We prove a hardness of sampling result for the anti-ferromagnetic Ising model on random graphs of average degree $d$ for large constant $d$, proving that when the normalized inverse temperature satisfies $β>1$ (asymptotically corresponding to the condensation threshold), then w.h.p. over the random graph there is no stable sampling algorithm that can output a sample close in $W_2$ distance to the Gibbs measure. The results also apply to a fixed-magnetization version of the model, showing that there are no stable sampling algorithms for low but positive temperature max and min bisection distributions. These results show a gap in the tractability of search and sampling problems: while there are efficient algorithms to find near optimizers, stable sampling algorithms cannot access the Gibbs distribution concentrated on such solutions. Our techniques involve extensions of the interpolation technique relating behavior of the mean field Sherrington-Kirkpatrick model to behavior of Ising models on random graphs of average degree $d$ for large $d$. While previous interpolation arguments compared the free energies of the two models, our argument compares the average energies and average overlaps in the two models.
