104 research outputs found
Fast Algorithms at Low Temperatures via Markov Chains
For spin systems, such as the hard-core model on independent sets weighted by fugacity lambda>0, efficient algorithms for the associated approximate counting/sampling problems typically apply in the high-temperature region, corresponding to low fugacity. Recent work of Jenssen, Keevash and Perkins (2019) yields an FPTAS for approximating the partition function (and an efficient sampling algorithm) on bounded-degree (bipartite) expander graphs for the hard-core model at sufficiently high fugacity, and also the ferromagnetic Potts model at sufficiently low temperatures. Their method is based on using the cluster expansion to obtain a complex zero-free region for the partition function of a polymer model, and then approximating this partition function using the polynomial interpolation method of Barvinok. We present a simple discrete-time Markov chain for abstract polymer models, and present an elementary proof of rapid mixing of this new chain under sufficient decay of the polymer weights. Applying these general polymer results to the hard-core and ferromagnetic Potts models on bounded-degree (bipartite) expander graphs yields fast algorithms with running time O(n log n) for the Potts model and O(n^2 log n) for the hard-core model, in contrast to typical running times of n^{O(log Delta)} for algorithms based on Barvinok\u27s polynomial interpolation method on graphs of maximum degree Delta. In addition, our approach via our polymer model Markov chain is conceptually simpler as it circumvents the zero-free analysis and the generalization to complex parameters. Finally, we combine our results for the hard-core and ferromagnetic Potts models with standard Markov chain comparison tools to obtain polynomial mixing time for the usual spin system Glauber dynamics restricted to even and odd or "red" dominant portions of the respective state spaces
Algorithms for #BIS-hard problems on expander graphs
We give an FPTAS and an efficient sampling algorithm for the high-fugacity hard-core model on bounded-degree bipartite expander graphs and the low-temperature ferromagnetic Potts model on bounded-degree expander graphs. The results apply, for example, to random (bipartite) Δ-regular graphs, for which no efficient algorithms were known for these problems (with the exception of the Ising model) in the non-uniqueness regime of the infinite Δ-regular tree
A Spectral Independence View on Hard Spheres via Block Dynamics
The hard-sphere model is one of the most extensively studied models in statistical physics. It describes the continuous distribution of spherical particles, governed by hard-core interactions. An important quantity of this model is the normalizing factor of this distribution, called the partition function. We propose a Markov chain Monte Carlo algorithm for approximating the grand-canonical partition function of the hard-sphere model in d dimensions. Up to a fugacity of ? < e/2^d, the runtime of our algorithm is polynomial in the volume of the system. This covers the entire known real-valued regime for the uniqueness of the Gibbs measure.
Key to our approach is to define a discretization that closely approximates the partition function of the continuous model. This results in a discrete hard-core instance that is exponential in the size of the initial hard-sphere model. Our approximation bound follows directly from the correlation decay threshold of an infinite regular tree with degree equal to the maximum degree of our discretization. To cope with the exponential blow-up of the discrete instance we use clique dynamics, a Markov chain that was recently introduced in the setting of abstract polymer models. We prove rapid mixing of clique dynamics up to the tree threshold of the univariate hard-core model. This is achieved by relating clique dynamics to block dynamics and adapting the spectral expansion method, which was recently used to bound the mixing time of Glauber dynamics within the same parameter regime
Ferromagnetic Potts Model: Refined #BIS-hardness and Related Results
Recent results establish for 2-spin antiferromagnetic systems that the
computational complexity of approximating the partition function on graphs of
maximum degree D undergoes a phase transition that coincides with the
uniqueness phase transition on the infinite D-regular tree. For the
ferromagnetic Potts model we investigate whether analogous hardness results
hold. Goldberg and Jerrum showed that approximating the partition function of
the ferromagnetic Potts model is at least as hard as approximating the number
of independent sets in bipartite graphs (#BIS-hardness). We improve this
hardness result by establishing it for bipartite graphs of maximum degree D. We
first present a detailed picture for the phase diagram for the infinite
D-regular tree, giving a refined picture of its first-order phase transition
and establishing the critical temperature for the coexistence of the disordered
and ordered phases. We then prove for all temperatures below this critical
temperature that it is #BIS-hard to approximate the partition function on
bipartite graphs of maximum degree D. As a corollary, it is #BIS-hard to
approximate the number of k-colorings on bipartite graphs of maximum degree D
when k <= D/(2 ln D).
The #BIS-hardness result for the ferromagnetic Potts model uses random
bipartite regular graphs as a gadget in the reduction. The analysis of these
random graphs relies on recent connections between the maxima of the
expectation of their partition function, attractive fixpoints of the associated
tree recursions, and induced matrix norms. We extend these connections to
random regular graphs for all ferromagnetic models and establish the Bethe
prediction for every ferromagnetic spin system on random regular graphs. We
also prove for the ferromagnetic Potts model that the Swendsen-Wang algorithm
is torpidly mixing on random D-regular graphs at the critical temperature for
large q.Comment: To appear in SIAM J. Computin
Delay performance in random-access grid networks
We examine the impact of torpid mixing and meta-stability issues on the delay
performance in wireless random-access networks. Focusing on regular meshes as
prototypical scenarios, we show that the mean delays in an toric
grid with normalized load are of the order . This
superlinear delay scaling is to be contrasted with the usual linear growth of
the order in conventional queueing networks. The intuitive
explanation for the poor delay characteristics is that (i) high load requires a
high activity factor, (ii) a high activity factor implies extremely slow
transitions between dominant activity states, and (iii) slow transitions cause
starvation and hence excessively long queues and delays. Our proof method
combines both renewal and conductance arguments. A critical ingredient in
quantifying the long transition times is the derivation of the communication
height of the uniformized Markov chain associated with the activity process. We
also discuss connections with Glauber dynamics, conductance and mixing times.
Our proof framework can be applied to other topologies as well, and is also
relevant for the hard-core model in statistical physics and the sampling from
independent sets using single-site update Markov chains
Ferromagnetic Potts Model: Refined #BIS-hardness and Related Results
Recent results establish for the hard-core model (and more generally for 2-spin antiferromagnetic systems) that the computational complexity of approximating the partition function on graphs of maximum degree D undergoes a phase transition that coincides with the uniqueness/non-uniqueness phase transition on the infinite D-regular tree. For the ferromagnetic Potts model we investigate whether analogous hardness results hold. Goldberg and Jerrum showed that approximating the partition function of the ferromagnetic Potts model is at least as hard as approximating the number of independent sets in bipartite graphs, so-called #BIS-hardness. We improve this hardness result by establishing it for bipartite graphs of maximum degree D. To this end, we first present a detailed picture for the phase diagram for the infinite D-regular tree, giving a refined picture of its first-order phase transition and establishing the critical temperature for the coexistence of the disordered and ordered phases. We then prove for all temperatures below this critical temperature (corresponding to the region where the ordered phase "dominates") that it is #BIS-hard to approximate the partition function on bipartite graphs of maximum degree D.
The #BIS-hardness result uses random bipartite regular graphs as a gadget in the reduction. The analysis of these random graphs relies on recent results establishing connections between the maxima of the expectation of their partition function, attractive fixpoints of the associated tree recursions, and induced matrix norms. In this paper we extend these connections to random regular graphs for all ferromagnetic models. Using these connections, we establish the Bethe prediction for every ferromagnetic spin system on random regular graphs, which says roughly that the expectation of the log of the partition function Z is the same as the log of the expectation of Z. As a further consequence of our results, we prove for the ferromagnetic Potts model that the Swendsen-Wang algorithm is torpidly mixing (i.e., exponentially slow convergence to its stationary distribution) on random D-regular graphs at the critical temperature for sufficiently large q
Algorithmic Pirogov-Sinai theory
We develop an efficient algorithmic approach for approximate counting and
sampling in the low-temperature regime of a broad class of statistical physics
models on finite subsets of the lattice and on the torus
. Our approach is based on combining contour
representations from Pirogov-Sinai theory with Barvinok's approach to
approximate counting using truncated Taylor series. Some consequences of our
main results include an FPTAS for approximating the partition function of the
hard-core model at sufficiently high fugacity on subsets of with
appropriate boundary conditions and an efficient sampling algorithm for the
ferromagnetic Potts model on the discrete torus at
sufficiently low temperature
- …