Abstract:We consider the problem of sampling from the Ising model when the underlying interaction matrix has eigenvalues lying within an interval of length $\gamma$. Recent work in this setting has shown various algorithmic results that apply roughly when $\gamma< 1$, notably with nearly-linear running times based on the classical Glauber dynamics. However, the optimality of the range of $\gamma$ was not clear since previous inapproximability results developed for the antiferromagnetic case (where the matrix has entries $\leq 0$) apply only for $\gamma>2$. To this end, Kunisky (SODA'24) recently provided evidence that the problem becomes hard already when $\gamma>1$ based on the low-degree hardness for an inference problem on random matrices. Based on this, he conjectured that sampling from the Ising model in the same range of $\gamma$ is NP-hard. Here we confirm this conjecture, complementing in particular the known algorithmic results by showing NP-hardness results for approximately counting and sampling when $\gamma>1$, with strong inapproximability guarantees; we also obtain a more refined hardness result for matrices where only a constant number of entries per row are allowed to be non-zero. The main observation in our reductions is that, for $\gamma>1$, Glauber dynamics mixes slowly when the interactions are all positive (ferromagnetic) for the complete and random regular graphs, due to a bimodality in the underlying distribution. While ferromagnetic interactions typically preclude NP-hardness results, here we work around this by introducing in an appropriate way mild antiferromagnetism, keeping the spectrum roughly within the same range. This allows us to exploit the bimodality of the aforementioned graphs and show the target NP-hardness by adapting suitably previous inapproximability techniques developed for antiferromagnetic systems.

Approximating the Spectral Gap of the Pólya-Gamma Gibbs Sampler

A spectral theorem for the semigroup generated by a class of infinitely many master equations

Spectral gap bounds for reversible hybrid Gibbs chains

The Spectral Gap of a Gaussian Quantum Markovian Generator

Improved Estimation of Relaxation Time in Non-reversible Markov Chains

Optimality of Glauber dynamics for general-purpose Ising model sampling and free energy approximation

On the properties of variational approximations of Gibbs posteriors

Controlling Uncertainty of Empirical First-Passage Times in the Small-Sample Regime

On randomized estimators of the Hafnian of a nonnegative matrix

Approximation of Pufferfish Privacy for Gaussian Priors

Asymptotic Gap Probability Distributions of the Gaussian Unitary Ensembles and Jacobi Unitary Ensembles

Adapting The Gibbs Sampler

Entropy contraction of the Gibbs sampler under log-concavity

On Sampling from Ising Models with Spectral Constraints

The Inverse Gamma-Gamma Prior for Optimal Posterior Contraction and Multiple Hypothesis Testing

Normal approximation for statistics of Gibbsian input in geometric probability

Ordered and size-biased frequencies in GEM and Gibbs models for species sampling

Comparison Theorems for the Mixing Times of Systematic and Random Scan Dynamics

Spectral Gaps via Imaginary Time

Perfect sampling from rapidly mixing Markov chains

Geometric ergodicity of Gibbs samplers for Bayesian error-in-variable regression