## Counting in Lattices: Combinatorial Problems from Statistical Mechanics (1994)

Citations: | 2 - 1 self |

### BibTeX

@TECHREPORT{Randall94countingin,

author = {Dana Randall},

title = {Counting in Lattices: Combinatorial Problems from Statistical Mechanics},

institution = {},

year = {1994}

}

### OpenURL

### Abstract

In this thesis we consider two classical combinatorial problems arising in statistical mechanics: counting matchings and self-avoiding walks in lattice graphs. The first problem arises in the study of the thermodynamical properties of monomers and dimers (diatomic molecules) in crystals. Fisher, Kasteleyn and Temperley discovered an elegant technique to exactly count the number of perfect matchings in two dimensional lattices, but it is not applicable for matchings of arbitrary size, or in higher dimensional lattices. We present the first efficient approximation algorithm for computing the number of matchings of any size in any periodic lattice in arbitrary dimension. The algorithm is based on Monte Carlo simulation of a suitable Markov chain and has rigorously derived performance guarantees that do not rely on any assumptions. In addition, we show that these results generalize to counting matchings in any graph which is the Cayley graph of a finite group. The second problem is countin...

### Citations

3839 |
Introduction to Automata Theory, Languages, and Computation
- Hopcroft, Ullman
- 1979
(Show Context)
Citation Context ...an algorithm takes as input a natural description of the problem instance and must compute an output in time which is bounded by some fixed polynomial in the size of the input description (see, e.g., =-=[31]-=-). For problems in statistical mechanics involving lattices, the input size is taken to be N , the number of vertices in the lattice; a polynomial-time algorithm (e.g., for the counting or generation ... |

2257 | Equation of state calculations by fast computing machines. The journal of chemical physics - Metropolis, Rosenbluth, et al. - 1953 |

508 |
The complexity of computing the permanent
- Valiant
- 1979
(Show Context)
Citation Context ...ional lattices. 1.2.2 Hardness Computational complexity theory offers further evidence for the apparent hardness of these problems. In 1979, Valiant defined the class #P to classify counting problems =-=[64]-=-. A counting problem asking for the number of elements in a set S n representing configurations on a graph G belongs to #P if there is a non-deterministic polynomial-time Turing machine with input n a... |

346 | Self-testing/correcting with applications to numerical problems
- Blum, Luby, et al.
- 1993
(Show Context)
Citation Context ...tions. This notion of self-testing, which either gives us confidence in our results or warns us that they may be erroneous, has been previously explored in the context of program checking (see, e.g., =-=[5]-=-). The conjecture we require is the following: Conjecture 2: For any dimension d , for some fixed polynomial g , c n c msg(n +m) c n+m ; 8n; m: This conjecture says that if we choose random self-avoid... |

319 | Statistical mechanics: rigorous results - Ruelle - 1969 |

314 | A polynomial-time approximation algorithm for the permanent of a matrix with non-negative entries - Jerrum, Sinclair, et al. - 2000 |

299 | Approximating the permanent
- Jerrum, Sinclair
- 1989
(Show Context)
Citation Context ...nly a polynomial number of steps we are close to stationarity. This is called rapid mixing and can be characterized by the expansion properties or conductancesof the graph underlying the Markov chain =-=[59, 34]-=-. We will explain this connection in greater detail in section 3.2.2. 1.4 Summary The remainder of this thesis is divided into two chapters. The first chapter addresses the monomer-dimer problem. We s... |

287 | Group Representations in Probability and Statistics - Diaconis - 1988 |

282 |
Random generation of combinatorial structures from a uniform distribution. Theoret
- Jerrum, Valiant, et al.
- 1986
(Show Context)
Citation Context ...g allows us to examine "typical" configurations of the system. Moreover, it has been recognized that for large classes of combinatorial problems the counting and sampling problems are closel=-=y related [36]-=-. For 10 each of the problems we consider we will be addressing both of these questions. 1.1.2 Free energy As indicated earlier, the partition function Z is the key to relating the microscopic and mac... |

277 | Geometric bounds for eigenvalues of Markov chains
- Diaconis, Stroock
- 1991
(Show Context)
Citation Context ...in terms of the quantity ae = min \Gamma ae(\Gamma) , using a measure known as conductance [60]. However, we can get a slightly sharper bound in this case by following an idea of Diaconis and Stroock =-=[12]-=- and using the alternative measure ae = min \Gamma ae(\Gamma)`(\Gamma) , where `(\Gamma) is the maximum length of a path in \Gamma . The appropriate version of this bound can be found by combining Pro... |

259 |
Approximate counting, uniform generation and rapidly mixing Markov chains, Inform. and Comput
- Sinclair, Jerrum
- 1989
(Show Context)
Citation Context ...kov chain should be rapidly mixing. This intuition can be formalized and a bound obtained on the mixing time in terms of the quantity ae = min \Gamma ae(\Gamma) , using a measure known as conductance =-=[60]-=-. However, we can get a slightly sharper bound in this case by following an idea of Diaconis and Stroock [12] and using the alternative measure ae = min \Gamma ae(\Gamma)`(\Gamma) , where `(\Gamma) is... |

179 |
Algorithms for Random Generation and Counting : A Markov Chain Approach
- Sinclair
- 1993
(Show Context)
Citation Context ...lly, we will be interested in designing algorithms which meet the following specifications. These are the standard definitions of approximation algorithms for counting and sampling (see, for example, =-=[39, 36, 59]-=-). Let S n be the set of configurations of size n on a lattice G of size N (for example, S n might be the set of self avoiding walks of length n on ansN 1=2 \ThetasN 1=2 lattice). Let a n be the numbe... |

178 | Improved bounds for mixing rates of Markov chains and multicommodity flow
- Sinclair
- 1992
(Show Context)
Citation Context ...ty is analyzed rigorously in [34] and shown to be a polynomial function of the ratio ff(G) = am\Gamma1 =am , specifically O(ff(G) 4 jEj 2 ) . (This has since been improved to O(ff(G) 2 jV jjEj) ; see =-=[58]-=-.) The Monte Carlo procedure will therefore be efficient for graphs G in which the ratio ff(G) is small. x Matchings in G of cardinality m are called perfect matchings, and those of cardinalitysm \Gam... |

157 | Intersections of random walks
- Lawler
- 1991
(Show Context)
Citation Context ...ticians and natural scientists for many years and are the subject of an extensive literature; for a state-of-the-art survey, see the recent book of Madras and Slade [48]. (See also the book by Lawler =-=[45]-=- for related topics.) One of the most important applications is as a model for the spatial arrangement of linear polymer molecules in chemical physics. Here the walk represents a molecule composed of ... |

132 |
The statistics of dimers on a lattice I. the number of dimer arrangements on a quadratic lattice
- Kasteleyn
- 1961
(Show Context)
Citation Context ...cal mechanics was achieved in 1961 when Fisher, Kasteleyn and Temperley independently discovered a polynomial-time algorithm for a special case of the monomer-dimer problem known as the dimer problem =-=[15,40, 63]-=- (see section 2.1). The dimer problem asks for the number of perfect matchings, or coverings of the lattice by N=2 dimers (and no monomers). (More generally, a perfect matching in a graph is a matchin... |

120 |
A two-dimensional model with an order-disorder transformation
- Onsager
- 1944
(Show Context)
Citation Context ...tions to the counting (or weighted counting) problems. In 1944, Onsager discovered such a closed-form expression for the partition function of the two-dimensional Ising model with zero external field =-=[53]-=-. It has not been possible to find such a precise solution for most other physical systems, or even for the Ising model in higher dimensions or in the presence of an external field. Consequently most ... |

117 |
The Self-Avoiding Walk
- Madras, Slade
- 1993
(Show Context)
Citation Context ...sition forces the self-avoidance condition. The dilute solution enables one to study a single polymer without worrying about further restrictions caused by interactions with other polymer chains (see =-=[48]-=- for a survey). The more general class of lattice animals arises similarly in the study of branched polymers in dilute solution. As we have already stated, for each of these physical systems, a primar... |

99 |
Graph theory and crystal physics
- Kasteleyn
- 1967
(Show Context)
Citation Context ...ested in counting the number of monomer-dimer arrangements with any fixed density. In this section we present some of the highlights in the history of this problem. For further information see, e.g., =-=[30, 42, 66]-=- and the references given there. The monomer-dimer problem gained prominence in 1937 through the early paper of Fowler and Rushbrooke [17]. A breakthrough was achieved in 1961, when, independently, Fi... |

89 |
Monte-Carlo approximation algorithms for enumeration problems
- Karp, Luby, et al.
- 1989
(Show Context)
Citation Context ...lly, we will be interested in designing algorithms which meet the following specifications. These are the standard definitions of approximation algorithms for counting and sampling (see, for example, =-=[39, 36, 59]-=-). Let S n be the set of configurations of size n on a lattice G of size N (for example, S n might be the set of self avoiding walks of length n on ansN 1=2 \ThetasN 1=2 lattice). Let a n be the numbe... |

85 |
Dimer statistics and phase transitions
- Kasteleyn
- 1963
(Show Context)
Citation Context ... efficient manner. In fact, this technique is more general and allows the number of dimer coverings of any planar graph (or indeed of any family of graphs with fixed genus) to be computed efficiently =-=[41]-=-. Unfortunately, these methods do not extend to two-dimensional lattices with dimer density less than 1, or to lattices in higher dimensions even when the dimer density remains 1. In fact, the three-d... |

80 |
Theory of monomer-dimer systems
- Heilmann, Lieb
- 1972
(Show Context)
Citation Context ...n-overlapping arrangement of dimers and monomers. A dimer is a diatomic molecule which covers two adjacent vertices in the lattice, and a monomer covers each vertex not covered by a dimer (see, e.g., =-=[30]-=- for a 8 t t t t t t t t t t t t t t t t t t t t t t t t t t t t t t t t f f f f f f f f Figure 1.2: A monomer-dimer arrangement. history of the problem). In graph theoretic terms, the monomer-dimer a... |

68 |
hard is it to marry at random? (On the approximation of the permanent
- BRODER, How
- 1986
(Show Context)
Citation Context ...hm fails: namely, counting matchings of arbitrary size, and counting matchings on non-planar lattices. These results are based on an approximation algorithm due to Jerrum and Sinclair [34] and Broder =-=[7]-=-. In each of these cases we also present an almost uniform generator for the corresponding sampling problems. This chapter is based on work with Claire Kenyon and Alistair Sinclair [43]. In chapter 3 ... |

62 | Computing the volume of convex bodies: a case where randomness provably helps
- Dyer, Frieze
- 1991
(Show Context)
Citation Context ...vement in the running time comes from two observations. The first is a standard trick for reducing the number of samples required to get a good estimate of the product of random variables (see, e.g., =-=[13]-=-). Rather than requiring that the variance of all n random variables be extremely small (so that the product will be close to its mean), we bound the variance of the product directly. This allows us t... |

59 | and Evgeny Lifshitz, Statistical physics - Landau - 1980 |

56 | Some inequalities for reversible Markov chains - Aldous - 1982 |

42 |
Two-dimensional monomer-dimer systems are computationally intractable
- Jerrum
- 1987
(Show Context)
Citation Context ... algorithm fails to generalize to all non-planar graphs. Furthermore, Jerrum showed that counting the total number of matchings (of all sizes) in a graph is #P-complete, even when the graph is planar =-=[32]-=-. This gives evidence for why the Fisher, Kasteleyn, Temperley algorithm does not easily generalize to counting matchings of arbitrary size. Similarly, consider the problem #SAW(G;sN \Gamma 1; x 0 ) w... |

37 |
Statistical mechanics of dimers on a plane lattice
- Fisher
- 1961
(Show Context)
Citation Context ...cal mechanics was achieved in 1961 when Fisher, Kasteleyn and Temperley independently discovered a polynomial-time algorithm for a special case of the monomer-dimer problem known as the dimer problem =-=[15,40, 63]-=- (see section 2.1). The dimer problem asks for the number of perfect matchings, or coverings of the lattice by N=2 dimers (and no monomers). (More generally, a perfect matching in a graph is a matchin... |

37 |
Dimer Problems in Statistical Mechanics-an exact result, Philos
- Temperley, Fisher
- 1961
(Show Context)
Citation Context ...cal mechanics was achieved in 1961 when Fisher, Kasteleyn and Temperley independently discovered a polynomial-time algorithm for a special case of the monomer-dimer problem known as the dimer problem =-=[15,40, 63]-=- (see section 2.1). The dimer problem asks for the number of perfect matchings, or coverings of the lattice by N=2 dimers (and no monomers). (More generally, a perfect matching in a graph is a matchin... |

35 |
An introduction to the Ising model
- Cipra
- 1987
(Show Context)
Citation Context ...ar lattice (or chessboard). Example 1: The Ising Model The Ising model was introduced in the 1920's to study ferromagnetism and is one of the most famous models from statistical mechanics (see, e.g., =-=[8]-=- for a history and review). The system is modeled by a finite lattice where the vertices represent atoms of a ferromagnetic 7 material. Each configuration oe consists of an assignment of a +1 or-1 spi... |

29 |
A mildly exponential approximation algorithm for the permanent, Algorithmica 16
- Jerrum, Vazirani
- 1996
(Show Context)
Citation Context ...ff(G) is large can be efficiently decomposed in such a way that the resulting components have a small value of ff , and hence fall within the scope of the Monte Carlo algorithm; this idea was used in =-=[37]-=- to obtain an approximation scheme for general graphs whose running time, though still exponential, improves substantially on naive deterministic methods. 33 The question of whether ff is polynomially... |

27 |
New Monte Carlo method for the self-avoiding walk
- Beretti, Sokal
- 1985
(Show Context)
Citation Context ...ial time algorithms where the statistical errors are rigorously controlled. Our algorithms are based on modifications and extensions of a Monte Carlo approach studied originally by Berretti and Sokal =-=[3]-=-. In the next subsection we sketch this approach and point out its limitations. Then, in section 3.1.3, we summarize our algorithms and explain how they overcome these problems. 36 3.1.2 Monte Carlo s... |

25 |
Poor man's Monte Carlo
- Hammersley, Morton
- 1954
(Show Context)
Citation Context ...r arbitrary dimensions d , they are of greatest interest in the case of low-dimensional lattices with 2sds4 . One key fact that holds in all dimensions was discovered in 1954 by Hammersley and Morton =-=[27]-=-; they observed that lim n!1 c 1=n n =sexists, and thatsnsc n =sn f(n) , where lim n!1 f(n) 1=n = 1 . This is a straightforward consequence of the obvious fact that 35 the sequence ` n = log c n is su... |

25 |
Bounds on the L spectrum for Markov chains and Markov processes: A generalization of Cheegerâ€™s inequality
- LAWLER, SOKAL
- 1988
(Show Context)
Citation Context ... 1, Sokal and Thomas [62] argue that with such values of fi the Markov chain is rapidly mixing, i.e., it gets very close to stationarity after a number of steps that is only polynomial in n (see also =-=[46]-=-). In order to appreciate the role of fi here, consider a truncated version of this Markov chain in which the length of a walk is never allowed to exceed n , so that the stationary distribution is alw... |

21 | The computational complexity of some classical problems from statistical physics
- Welsh
- 1990
(Show Context)
Citation Context ...ested in counting the number of monomer-dimer arrangements with any fixed density. In this section we present some of the highlights in the history of this problem. For further information see, e.g., =-=[30, 42, 66]-=- and the references given there. The monomer-dimer problem gained prominence in 1937 through the early paper of Fowler and Rushbrooke [17]. A breakthrough was achieved in 1961, when, independently, Fi... |

17 |
Upper bounds for the connective constant of self-avoiding walks
- Alm
- 1993
(Show Context)
Citation Context ...h have also been employed to obtain numerical estimates for the constantssand fl . Elementary considerations show thats2 (d; 2d \Gamma 1) . For d = 2 , it has actually been proven thats2 (2:62; 2:70) =-=[2, 10]-=-. (See also [29] for similar bounds in higher dimensions.) However, these rigorous bounds are much weaker than the non-rigorous estimates obtained by empirical methods, which are typically quoted to f... |

17 | On the number of Eulerian orientations of a graph - Mihail, Winkler - 1992 |

17 | Self-testing algorithms for self-avoiding walks
- Randall, Sinclair
(Show Context)
Citation Context ...ations as well. These are the first algorithms for counting and generating self-avoiding walks where the error bounds are rigorously controlled. These results are based on work with Alistair Sinclair =-=[54]-=-. An optimized version of the algorithm was implemented by Rob Pike yielding numerical estimates for the number of self-avoiding walks in 2 and 3 dimensions; these estimates are presented in section 3... |

16 |
Solution of the Dimer Problem on an Hexagonal Lattice with Boundary
- Elser
- 1984
(Show Context)
Citation Context ...e dimer covering problem for this planar lattice has been known for some time [65, 41]. In contrast to the rectangular lattice, the assumption of periodic boundary conditions is important here: Elser =-=[14]-=- has solved the dimer covering problem on a hexagonal lattice with fixed boundaries, and shown that the result depends significantly on the shape of the boundary. 21 has at most jEj neighbors, where j... |

16 |
Further results on the rate of convergence to the connective constant of the hypercubical lattice, Quart
- Hammersley, Welsh
- 1962
(Show Context)
Citation Context ...) , where lim n!1 f(n) 1=n = 1 . This is a straightforward consequence of the obvious fact that 35 the sequence ` n = log c n is subadditive, i.e., ` n+ms` n + ` m for all n; m . Hammersley and Welsh =-=[28]-=- later showed that f(n) = O(a n 1=2 ) for some constant a . It is a celebrated and long-standing conjecture that f(n) is in fact polynomially bounded. More precisely, we have: Conjecture 1: c n =sn e ... |

15 |
Existence theorems and Monte Carlo methods for the monomer-dimer problem
- HAMMERSLEY
- 1966
(Show Context)
Citation Context ...rectangular lattice with dimers, is one of the classical unsolved problems of solid-state chemistry. A few facts are known: for example, ln(f(n))=n 3 tends 19 to a finite limitsas n tends to infinity =-=[21]-=-. Hammersley [22] proved the lower bounds0:418347 , while the early paper by Fowler and Rushbrooke [17] showed the upper bounds0:54931 . It has been conjectured thatslies between 0:43 and 0:45 . In ot... |

13 |
Statistical theory of perfect solutions
- Fowler, Rushbrooke
- 1937
(Show Context)
Citation Context ...y of this problem. For further information see, e.g., [30, 42, 66] and the references given there. The monomer-dimer problem gained prominence in 1937 through the early paper of Fowler and Rushbrooke =-=[17]-=-. A breakthrough was achieved in 1961, when, independently, Fisher, Kasteleyn and Temperley provided an analytic solution for the case of dimer coveringss(i.e., arrangements with dimer density 1) on a... |

12 |
Statistical mechanics of dimers on a plane lattice. II. Dimer correlations and monomers
- Fisher, Stephenson
- 1963
(Show Context)
Citation Context ...y shed more light on other quantities related to monomerdimer systems, and in particular the correlation between monomers at two specified sites, as studied in two dimensions by Fisher and Stephenson =-=[16]-=-. The remainder of the chapter is organized as follows. In the next section, we prove bounds of the above form for rectangular lattices with periodic boundary conditions in any dimension, and with fix... |

12 |
A Monte Carlo algorithm for estimating the permanent
- Karmarkar, Karp, et al.
- 1993
(Show Context)
Citation Context ...y other currently known algorithm, so it is of interest to establish precisely which graphs are amenable to it. (For other simpler, but apparently less widely applicable approximation algorithms, see =-=[38, 33]-=- and [55].) Moreover, it is conceivable that any graph G for which ff(G) is large can be efficiently decomposed in such a way that the resulting components have a small value of ff , and hence fall wi... |

12 |
Approximating the Permanent: A Simple Approach ,Random Structures and Algorithms 5
- Rasmussen
- 1994
(Show Context)
Citation Context ...ntly known algorithm, so it is of interest to establish precisely which graphs are amenable to it. (For other simpler, but apparently less widely applicable approximation algorithms, see [38, 33] and =-=[55]-=-.) Moreover, it is conceivable that any graph G for which ff(G) is large can be efficiently decomposed in such a way that the resulting components have a small value of ff , and hence fall within the ... |

11 | Algorithmic research problems in molecular bioinformatics - Lengauer - 1993 |

9 |
A lower bound for the monomer-dimer problem
- Hammersley, Menon
- 1970
(Show Context)
Citation Context ...rigorous results holds for the problem at dimer densities less than 1, even in two dimensions. Notable exceptions are series expansions valid at low densities [19] and lower bounds on the free energy =-=[6, 26]-=-. 2.1.2 Results We make progress on the monomer-dimer problem in cases where the technique of Fisher, Kasteleyn and Temperley fails. Specifically, we give a polynomial-time algorithm for computing, to... |

8 | The number of polygons on a lattice - Hammersley - 1961 |

8 |
An analysis of a Monte Carlo algorithm for estimating the permanent
- Jerrum
- 1991
(Show Context)
Citation Context ...y other currently known algorithm, so it is of interest to establish precisely which graphs are amenable to it. (For other simpler, but apparently less widely applicable approximation algorithms, see =-=[38, 33]-=- and [55].) Moreover, it is conceivable that any graph G for which ff(G) is large can be efficiently decomposed in such a way that the resulting components have a small value of ff , and hence fall wi... |

7 |
Exact series-expansion study of the monomer-dimer problem, Phys. Rev
- Gaunt
- 1969
(Show Context)
Citation Context ...f(n) to good accuracy. A similar lack of rigorous results holds for the problem at dimer densities less than 1, even in two dimensions. Notable exceptions are series expansions valid at low densities =-=[19]-=- and lower bounds on the free energy [6, 26]. 2.1.2 Results We make progress on the monomer-dimer problem in cases where the technique of Fisher, Kasteleyn and Temperley fails. Specifically, we give a... |

7 |
Calculations of lattice statistics
- Hammersley
- 1970
(Show Context)
Citation Context ...nomially with n for any fixed dimension d . A similar bound holds for arbitrary Cayley graphs. We stress that our Monte Carlo algorithm differs from earlier ones for monomer-dimer systems (see, e.g., =-=[23]-=-) in that it is guaranteed (independent of any heuristic arguments) to provide statistically reliable estimates in a running time that grows only polynomially with the number of lattice sites. Our pro... |