Abstract:Over the past decades, a fascinating computational phase transition has been identified in sampling from Gibbs distributions. Though, the computational complexity at the critical point remains poorly understood, as previous algorithmic and hardness results all required a constant slack from this threshold. In this paper, we resolve this open question at the critical phase transition threshold, thus completing the picture of the computational phase transition. We show that for the hardcore model on graphs with maximum degree $\Delta\ge 3$ at the uniqueness threshold $\lambda = \lambda_c(\Delta)$, the mixing time of Glauber dynamics is upper bounded by a polynomial in $n$, but is not nearly linear in the worst case. For the Ising model (either antiferromagnetic or ferromagnetic), we establish similar results. For the Ising model on graphs with maximum degree $\Delta\ge 3$ at the critical temperature $\beta$ where $|\beta| = \beta_c(\Delta)$, with the tree-uniqueness threshold $\beta_c(\Delta)$, we show that the mixing time of Glauber dynamics is upper bounded by $\tilde{O}\left(n^{2 + O(1/\Delta)}\right)$ and lower bounded by $\Omega\left(n^{3/2}\right)$ in the worst case. For the Ising model specified by a critical interaction matrix $J$ with $\left \lVert J \right \rVert_2=1$, we obtain an upper bound $\tilde{O}(n^{3/2})$ for the mixing time, matching the lower bound $\Omega\left(n^{3/2}\right)$ on the complete graph up to a logarithmic factor. Our mixing time upper bounds are derived from a new interpretation and analysis of the localization scheme method introduced by Chen and Eldan (2022), applied to the field dynamics for the hardcore model and the proximal sampler for the Ising model. As key steps in both our upper and lower bounds, we establish sub-linear upper and lower bounds for spectral independence at the critical point for worst-case instances.

Trickle-Down in Localization Schemes and Applications

Local Escape Rates for -Mixing Dynamical Systems

Locally Stationary Distributions: A Framework for Analyzing Slow-Mixing Markov Chains

Rapid Mixing via Coupling Independence for Spin Systems with Unbounded Degree

Sampling, Diffusions, and Stochastic Localization

Rapid mixing of global Markov chains via spectral independence: the unbounded degree case

Rapid Mixing at the Uniqueness Threshold

Weak Poincaré Inequalities, Simulated Annealing, and Sampling from Spherical Spin Glasses

Rapid mixing of the down-up walk on matchings of a fixed size

Spectral Independence via Stability and Applications to Holant-Type Problems

Global Convergence of Langevin Dynamics Based Algorithms for Nonconvex Optimization

Krickeberg mixing for Z extensions of Gibbs Markov semiflows

Rapid Mixing from Spectral Independence Beyond the Boolean Domain

Parallelising Glauber dynamics

Rapid Mixing of Glauber Dynamics Via Spectral Independence for All Degrees

Improved Distributed Algorithms for Random Colorings

Projected Langevin dynamics and a gradient flow for entropic optimal transport

On Mixing of Markov Chains: Coupling, Spectral Independence, and Entropy Factorization

Universality of Spectral Independence with Applications to Fast Mixing in Spin Glasses

Sharper Bounds for Uniformly Stable Algorithms with Stationary Mixing Process

A Langevin-like Sampler for Discrete Distributions