Abstract: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 $\lambda$ then spectral independence holds at $\lambda$. As a consequence, for Holant-type problems (e.g., spin systems) on bounded-degree graphs, we obtain optimal $O(n\log 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(n\log{n})$ 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(n\log{n})$ 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.

Numerical Simulation of Polynomial-Speed Convergence Phenomenon

On the polynomial convergence rate to nonequilibrium steady-states

Probabilistic limit behaviors of numerical discretizations for time-homogeneous Markov processes

Numerical Simulation of Interactions Between Rigid Rod-Like Polymers and Coherent Structures in a Mixing Layer

Numerical Simulation of the Melting of Solid Particles in Thermal Convection with a Modified Immersed Boundary Method

Stochastic models of polymeric fluids at small Deborah number

Convergence of a Stochastic Method for the Modeling of Polymeric Fluids

Numerical simulation for non-linear thermal wave

Convergence Speed and Approximation Accuracy of Numerical MCMC

Long-time Convergence of Numerical Approximations for Semilinear Parabolic Equations (Ⅰ)

Spectral Independence via Stability and Applications to Holant-Type Problems

Convergence of a Metropolized Integrator for Stochastic Differential Equations with Variable Diffusion Coefficient

Numerical Simulation of Microflows using Hermite Spectral Methods

Numerical solutions to time-fractional stochastic partial differential equations

Numerical Approximations and Convergence Analysis of Piecewise Diffusion Markov Processes, with Application to Glioma Cell Migration

Numerical simulation for polymer solutions based on a hybrid computation of multi-particle collision dynamics and coarse-grained molecular dynamics

Probability-conservative simulation for Lévy financial model by a mixed finite element method

Multi-frequency averaging and uniform accuracy towards numerical approximations for a Bloch model

Numerical study of condensation in a Fermi-like model of conterflowing particles via Gini coefficient

Numerical analysis of a full discretization for stochastic Cahn--Hilliard equation driven by additive noise

Direct numerical simulation of water-ethanol flows in a T-mixer