Spectral Independence via Stability and Applications to Holant-Type Problems

by   Zongchen Chen, et al.

This paper formalizes connections between stability of polynomials and convergence rates of Markov Chain Monte Carlo (MCMC) algorithms. We prove that if a (multivariate) partition function is nonzero in a region around a real point λ then spectral independence holds at λ. As a consequence, for Holant-type problems (e.g., spin systems) on bounded-degree graphs, we obtain optimal O(nlog n) mixing time bounds for the single-site update Markov chain known as the Glauber dynamics. Our result significantly improves the running time guarantees obtained via the polynomial interpolation method of Barvinok (2017), refined by Patel and Regts (2017). There are a variety of applications of our results. In this paper, we focus on Holant-type (i.e., edge-coloring) problems, including weighted edge covers and weighted even subgraphs. For the weighted edge cover problem (and several natural generalizations) we obtain an O(nlogn) sampling algorithm on bounded-degree graphs. The even subgraphs problem corresponds to the high-temperature expansion of the ferromagnetic Ising model. We obtain an O(nlogn) sampling algorithm for the ferromagnetic Ising model with a nonzero external field on bounded-degree graphs, which improves upon the classical result of Jerrum and Sinclair (1993) for this class of graphs. We obtain further applications to antiferromagnetic two-spin models on line graphs, weighted graph homomorphisms, tensor networks, and more.


page 1

page 2

page 3

page 4


Fast Sampling of b-Matchings and b-Edge Covers

For an integer b ≥ 1, a b-matching (resp. b-edge cover) of a graph G=(V,...

Optimal Mixing of Glauber Dynamics: Entropy Factorization via High-Dimensional Expansion

We prove an optimal mixing time bound on the single-site update Markov c...

Sampling from the low temperature Potts model through a Markov chain on flows

In this paper we consider the algorithmic problem of sampling from the P...

Polynomial-time approximation algorithms for the antiferromagnetic Ising model on line graphs

We present a polynomial-time Markov chain Monte Carlo algorithm for esti...

Fast and perfect sampling of subgraphs and polymer systems

We give an efficient perfect sampling algorithm for weighted, connected ...

Dimension reduction for maximum matchings and the Fastest Mixing Markov Chain

Let G = (V,E) be an undirected graph with maximum degree Δ and vertex co...

Zeros of ferromagnetic 2-spin systems

We study zeros of the partition functions of ferromagnetic 2-state spin ...

Please sign up or login with your details

Forgot password? Click here to reset