Counting without sampling: Asymptotics of the log‐partition function for certain statistical physics models

In this article we propose new methods for computing the asymptotic value for the logarithm of the partition function (free energy) for certain statistical physics models on certain type of finite graphs, as the size of the underlying graph goes to infinity. The two models considered are the hard-core (independent set) model when the activity parameter λ is small, and also the Potts (q-coloring) model. We only consider the graphs with large girth. In particular, we prove that asymptotically the logarithm of the number of independent sets of any r-regular graph with large girth when rescaled is approximately constant if r ≤ 5. For example, we show that every 4-regular n-node graph with large girth has approximately (1.494 · · · )n-many independent sets, for large n. Further, we prove that for every r-regular graph with r ≥ 2, with n nodes and large girth, the number of proper q ≥ r +1 colorings is approximately [q(1 − 1 q ) r 2 ]n, for large n. We also show that these results hold for random regular graphs with high probability (w.h.p.) as well. As a byproduct of our method we obtain simple algorithms for the problem of computing approximately the logarithm of the number of independent sets and proper colorings, in low degree Correspondence to: David Gamarnik *Preliminary version of this paper appeared in Proceedings of the Seventeenth ACM-SIAM Symposium on Discrete Algorithms 2006. †ArXive version http://front.math.ucdavis.edu/math.PR/0510471 ‡Supported by Swedish Natural Science Research Council. © 2008 Wiley Periodicals, Inc.

[1]  M. Mézard,et al.  Reconstruction on Trees and Spin Glass Transition , 2005, cond-mat/0512295.

[2]  Mark Jerrum,et al.  The Markov chain Monte Carlo method: an approach to approximate counting and integration , 1996 .

[3]  Mark Jerrum,et al.  Approximating the Permanent , 1989, SIAM J. Comput..

[4]  Peter Winkler,et al.  Graph Homomorphisms and Long Range Action , 2001, Graphs, Morphisms and Statistical Physics.

[5]  Peter Winkler,et al.  A second threshold for the hard‐core model on a Bethe lattice , 2004, Random Struct. Algorithms.

[6]  G. Brightwell,et al.  Random colorings of a cayley tree , 2002 .

[7]  F. Spitzer Markov Random Fields on an Infinite Tree , 1975 .

[8]  Martin E. Dyer,et al.  A random polynomial-time algorithm for approximating the volume of convex bodies , 1991, JACM.

[9]  M. Simonovits,et al.  Random walks and an O * ( n 5 ) volume algorithm for convex bodies , 1997 .

[10]  Johan Wästlund AN EASY PROOF OF THE ζ ( 2 ) LIMIT IN THE RANDOM ASSIGNMENT PROBLEM , 2006 .

[11]  J. Michael Steele,et al.  The Objective Method: Probabilistic Combinatorial Optimization and Local Weak Convergence , 2004 .

[12]  Cristopher Moore,et al.  The Chromatic Number of Random Regular Graphs , 2004, APPROX-RANDOM.

[13]  Dror Weitz,et al.  Counting independent sets up to the tree threshold , 2006, STOC '06.

[14]  Eric Vigoda,et al.  A polynomial-time approximation algorithm for the permanent of a matrix with non-negative entries , 2001, STOC '01.

[15]  J. Jonasson Uniqueness of uniform random colorings of regular trees , 2002 .

[16]  Jon Warren,et al.  Dynamics and endogeny for recursive processes on trees , 2005, math/0506038.

[17]  F. Papangelou GIBBS MEASURES AND PHASE TRANSITIONS (de Gruyter Studies in Mathematics 9) , 1990 .

[18]  Hans-Otto Georgii,et al.  Gibbs Measures and Phase Transitions , 1988 .

[19]  David Gamarnik Linear phase transition in random linear constraint satisfaction problems , 2004, SODA '04.

[20]  Dror Weitz Combinatorial criteria for uniqueness of Gibbs measures , 2005 .

[21]  Elchanan Mossel,et al.  Information flow on trees , 2001, math/0107033.

[22]  James B. Martin Reconstruction Thresholds on Regular Trees , 2003, DRW.

[23]  Elchanan Mossel,et al.  Survey: Information Flow on Trees , 2004 .

[24]  Martin E. Dyer,et al.  Counting and sampling H-colourings? , 2004, Inf. Comput..

[25]  David Gamarnik,et al.  Maximum Weight Independent Sets and Matchings in Sparse Random Graphs. Exact Results Using the Local Weak Convergence Method , 2004, APPROX-RANDOM.

[26]  M. Mézard,et al.  Spin Glass Theory and Beyond , 1987 .

[27]  Leslie G. Valiant,et al.  The Complexity of Computing the Permanent , 1979, Theor. Comput. Sci..

[28]  M. Mézard,et al.  Glass models on Bethe lattices , 2003, cond-mat/0307569.

[29]  Frank Kelly,et al.  Stochastic Models of Computer Communication Systems , 1985 .

[30]  Santosh S. Vempala,et al.  Simulated annealing in convex bodies and an O*(n4) volume algorithm , 2006, J. Comput. Syst. Sci..

[31]  Alan M. Frieze,et al.  Random graphs , 2006, SODA '06.

[32]  D. Aldous,et al.  A survey of max-type recursive distributional equations , 2004, math/0401388.

[33]  M. Mézard,et al.  The Cavity Method at Zero Temperature , 2002, cond-mat/0207121.

[34]  M. Talagrand The parisi formula , 2006 .

[35]  D. Gamarnik,et al.  Invariant probability measures and dynamics of exponential linear type maps , 2008, Ergodic Theory and Dynamical Systems.

[36]  Martin E. Dyer,et al.  Randomly coloring constant degree graphs , 2004, Random Struct. Algorithms.

[37]  Miklós Simonovits,et al.  Random walks and an O*(n5) volume algorithm for convex bodies , 1997, Random Struct. Algorithms.

[38]  Dorit S. Hochbaum,et al.  Approximation Algorithms for NP-Hard Problems , 1996 .