We study the low-temperature $(2+1)$D solid-on-solid model on with zero boundary conditions and nonnegative heights (a floor at height $0$). Caputo et al. (2016) established that this random surface typically admits either $\mathfrak h $ or $\mathfrak h+1$ many nested macroscopic level line loops $\{\mathcal L_i\}_{i\geq 0}$ for an explicit $\mathfrak h\asymp \log L$, and its top loop $\mathcal L_0$ has cube-root fluctuations: For example, if $\rho (x)$ is the vertical displacement of $\mathcal L_0$ from the bottom boundary point $(x,0)$, then $\max \rho (x) = L^{1/3+o(1)}$ over . It is believed that rescaling $\rho $ by $L^{1/3}$ and $I_0$ by $L^{2/3}$ would yield a limit law of a diffusion on $[-1,1]$. However, no nontrivial lower bound was known on $\rho (x)$ for a fixed $x\in I_0$ (e.g., $x=\frac L2$), let alone on $\min \rho (x)$ in $I_0$, to complement the bound on $\max \rho (x)$. Here, we show a lower bound of the predicted order $L^{1/3}$: For every $\epsilon>0$, there exists $\delta>0$ such that $\min _{x\in I_0} \rho (x) \geq \delta L^{1/3}$ with probability at least $1-\epsilon $. The proof relies on the Ornstein–Zernike machinery due to Campanino–Ioffe–Velenik and a result of Ioffe, Shlosman and Toninelli (2015) that rules out pinning in Ising polymers with modified interactions along the boundary. En route, we refine the latter result into a Brownian excursion limit law, which may be of independent interest. We further show that in a $ K L^{2/3}\times K L^{2/3}$ box with boundary conditions $\mathfrak h-1,\mathfrak h,\mathfrak h,\mathfrak h$ (i.e., $\mathfrak h-1$ on the bottom side and $\mathfrak h$ elsewhere), the limit of $\rho (x)$ as $K,L\to \infty $ is a Ferrari–Spohn diffusion.

Many connections and dualities in representation theory and Lie theory can be explained using quasi-hereditary covers in the sense of Rouquier. Recent work by the first-named author shows that relative dominant (and codominant) dimensions are natural tools to classify and distinguish distinct quasi-hereditary covers of a finite-dimensional algebra. In this paper, we prove that the relative dominant dimension of a quasi-hereditary algebra, possessing a simple preserving duality, with respect to a direct summand of the characteristic tilting module is always an even number or infinite and that this homological invariant controls the quality of quasi-hereditary covers that possess a simple preserving duality. To resolve the Temperley–Lieb algebras, we apply this result to the class of Schur algebras $S(2, d)$ and their $q$-analogues. Our second main result completely determines the relative dominant dimension of $S(2, d)$ with respect to $Q=V^{\otimes d}$, the $d$-th tensor power of the natural two-dimensional module. As a byproduct, we deduce that Ringel duals of $q$-Schur algebras $S(2,d)$ give rise to quasi-hereditary covers of Temperley–Lieb algebras. Further, we obtain precisely when the Temperley–Lieb algebra is Morita equivalent to the Ringel dual of the $q$-Schur algebra $S(2, d)$ and precisely how far these two algebras are from being Morita equivalent, when they are not. These results are compatible with the integral setup, and we use them to deduce that the Ringel dual of a $q$-Schur algebra over the ring of Laurent polynomials over the integers together with some projective module is the best quasi-hereditary cover of the integral Temperley–Lieb algebra.

We give algorithms for approximating the partition function of the ferromagnetic $q$-color Potts model on graphs of maximum degree $d$. Our primary contribution is a fully polynomial-time approximation scheme for $d$-regular graphs with an expansion condition at low temperatures (that is, bounded away from the order-disorder threshold). The expansion condition is much weaker than in previous works; for example, the expansion exhibited by the hypercube suffices. The main improvements come from a significantly sharper analysis of standard polymer models; we use extremal graph theory and applications of Karger’s algorithm to count cuts that may be of independent interest. It is #BIS-hard to approximate the partition function at low temperatures on bounded-degree graphs, so our algorithm can be seen as evidence that hard instances of #BIS are rare. We also obtain efficient algorithms in the Gibbs uniqueness region for bounded-degree graphs. While our high-temperature proof follows more standard polymer model analysis, our result holds in the largest-known range of parameters $d$ and $q$.

Given the full shift over a countable state space on a countable amenable group, we develop its thermodynamic formalism. First, we introduce the concept of pressure and, using tiling techniques, prove its existence and further properties, such as an infimum rule. Next, we extend the definitions of different notions of Gibbs measures and prove their existence and equivalence, given some regularity and normalization criteria on the potential. Finally, we provide a family of potentials that nontrivially satisfy the conditions for having this equivalence and a nonempty range of inverse temperatures where uniqueness holds.

We study the asymptotic limit of random pure dimer coverings on rail yard graphs when the mesh sizes of the graphs go to 0. Each pure dimer covering corresponds to a sequence of interlacing partitions starting with an empty partition and ending in an empty partition. Under the assumption that the probability of each dimer covering is proportional to the product of weights of present edges, we obtain the limit shape (law of large numbers) of the rescaled height functions and the convergence of the unrescaled height fluctuations to a diffeomorphic image of the Gaussian free field (Central Limit Theorem), answering a question in [7]. Applications include the limit shape and height fluctuations for pure steep tilings [9] and pyramid partitions [20; 36; 39; 38]. The technique to obtain these results is to analyze a class of Macdonald processes which involve dual partitions as well.

We study the locations of complex zeroes of independence polynomials of bounded-degree hypergraphs. For graphs, this is a long-studied subject with applications to statistical physics, algorithms, and combinatorics. Results on zero-free regions for bounded-degree graphs include Shearer’s result on the optimal zero-free disc, along with several recent results on other zero-free regions. Much less is known for hypergraphs. We make some steps towards an understanding of zero-free regions for bounded-degree hypergaphs by proving that all hypergraphs of maximum degree $\Delta$ have a zero-free disc almost as large as the optimal disc for graphs of maximum degree $\Delta$ established by Shearer (of radius $\sim 1/(e \Delta )$). Up to logarithmic factors in $\Delta$ this is optimal, even for hypergraphs with all edge sizes strictly greater than $2$. We conjecture that for $k\ge 3$, $k$-uniform linear hypergraphs have a much larger zero-free disc of radius $\Omega (\Delta ^{- \frac{1}{k-1}} )$. We establish this in the case of linear hypertrees.

We consider a class of weakly interacting particle systems of mean-field type. The interactions between the particles are encoded in a graph sequence, i.e. two particles are interacting if and only if they are connected in the underlying graph. We establish a law of large numbers for the empirical measure of the system that holds whenever the graph sequence is convergent to a graphon. The limit is the solution of a non-linear Fokker–Planck equation weighted by the (possibly random) graphon limit. In contrast with the existing literature, our analysis focuses on both deterministic and random graphons: no regularity assumptions are made on the graph limit and we are able to include general graph sequences such as exchangeable random graphs. Finally, we identify the sequences of graphs, both random and deterministic, for which the associated empirical measure converges to the classical McKean–Vlasov mean-field limit.

Given a locally finite graph $\Gamma $, an amenable subgroup G of graph automorphisms acting freely and almost transitively on its vertices, and a G-invariant activity function $\unicode{x3bb} $, consider the free energy $f_G(\Gamma ,\unicode{x3bb} )$ of the hardcore model defined on the set of independent sets in $\Gamma $ weighted by $\unicode{x3bb} $. Under the assumption that G is finitely generated and its word problem can be solved in exponential time, we define suitable ensembles of hardcore models and prove the following: if $\|\unicode{x3bb} \|_\infty < \unicode{x3bb} _c(\Delta )$, there exists a randomized $\epsilon $-additive approximation scheme for $f_G(\Gamma ,\unicode{x3bb} )$ that runs in time $\mathrm {poly}((1+\epsilon ^{-1})\lvert \Gamma /G \rvert )$, where $\unicode{x3bb} _c(\Delta )$ denotes the critical activity on the $\Delta $-regular tree. In addition, if G has a finite index linearly ordered subgroup such that its algebraic past can be decided in exponential time, we show that the algorithm can be chosen to be deterministic. However, we observe that if $\|\unicode{x3bb} \|_\infty> \unicode{x3bb} _c(\Delta )$, there is no efficient approximation scheme, unless $\mathrm {NP} = \mathrm {RP}$. This recovers the computational phase transition for the partition function of the hardcore model on finite graphs and provides an extension to the infinite setting. As an application in symbolic dynamics, we use these results to develop efficient approximation algorithms for the topological entropy of subshifts of finite type with enough safe symbols, we obtain a representation formula of pressure in terms of random trees of self-avoiding walks, and we provide new conditions for the uniqueness of the measure of maximal entropy based on the connective constant of a particular associated graph.

We show that the $4$ -state anti-ferromagnetic Potts model with interaction parameter $w\in (0,1)$ on the infinite $(d+1)$ -regular tree has a unique Gibbs measure if $w\geq 1-\dfrac{4}{d+1_{_{\;}}}$ for all $d\geq 4$ . This is tight since it is known that there are multiple Gibbs measures when $0\leq w\lt 1-\dfrac{4}{d+1}$ and $d\geq 4$ . We moreover give a new proof of the uniqueness of the Gibbs measure for the $3$ -state Potts model on the $(d+1)$ -regular tree for $w\geq 1-\dfrac{3}{d+1}$ when $d\geq 3$ and for $w\in (0,1)$ when $d=2$ .

In this text, we provide a fully rigorous and complete proof of E.H. Lieb’s statement that (topological) entropy of square ice (or six-vertex model, XXZ spin chain for anisotropy parameter $\Delta =1/2$ ) is equal to $\tfrac 32\log _{2} (4/3)$ .

We determine the asymptotics of the number of independent sets of size $\lfloor \beta 2^{d-1} \rfloor$ in the discrete hypercube $Q_d = \{0,1\}^d$ for any fixed $\beta \in (0,1)$ as $d \to \infty$ , extending a result of Galvin for $\beta \in (1-1/\sqrt{2},1)$ . Moreover, we prove a multivariate local central limit theorem for structural features of independent sets in $Q_d$ drawn according to the hard-core model at any fixed fugacity $\lambda>0$ . In proving these results we develop several general tools for performing combinatorial enumeration using polymer models and the cluster expansion from statistical physics along with local central limit theorems.

We consider the problem of computing the partition function $\sum _x e^{f(x)}$ , where $f: \{-1, 1\}^n \longrightarrow {\mathbb R}$ is a quadratic or cubic polynomial on the Boolean cube $\{-1, 1\}^n$ . In the case of a quadratic polynomial f, we show that the partition function can be approximated within relative error $0 < \epsilon < 1$ in quasi-polynomial $n^{O(\ln n - \ln \epsilon )}$ time if the Lipschitz constant of the non-linear part of f with respect to the $\ell ^1$ metric on the Boolean cube does not exceed $1-\delta $ , for any $\delta>0$ , fixed in advance. For a cubic polynomial f, we get the same result under a somewhat stronger condition. We apply the method of polynomial interpolation, for which we prove that $\sum _x e^{\tilde {f}(x)} \ne 0$ for complex-valued polynomials $\tilde {f}$ in a neighborhood of a real-valued f satisfying the above mentioned conditions. The bounds are asymptotically optimal. Results on the zero-free region are interpreted as the absence of a phase transition in the Lee–Yang sense in the corresponding Ising model. The novel feature of the bounds is that they control the total interaction of each vertex but not every single interaction of sets of vertices.

We present a polynomial-time Markov chain Monte Carlo algorithm for estimating the partition function of the antiferromagnetic Ising model on any line graph. The analysis of the algorithm exploits the ‘winding’ technology devised by McQuillan [CoRR abs/1301.2880 (2013)] and developed by Huang, Lu and Zhang [Proc. 27th Symp. on Disc. Algorithms (SODA16), 514–527]. We show that exact computation of the partition function is #P-hard, even for line graphs, indicating that an approximation algorithm is the best that can be expected. We also show that Glauber dynamics for the Ising model is rapidly mixing on line graphs, an example being the kagome lattice.

Random constraint satisfaction problems play an important role in computer science and combinatorics. For example, they provide challenging benchmark examples for algorithms, and they have been harnessed in probabilistic constructions of combinatorial structures with peculiar features. In an important contribution (Krzakala et al. 2007, Proc. Nat. Acad. Sci.), physicists made several predictions on the precise location and nature of phase transitions in random constraint satisfaction problems. Specifically, they predicted that their satisfiability thresholds are quite generally preceded by several other thresholds that have a substantial impact both combinatorially and computationally. These include the condensation phase transition, where long-range correlations between variables emerge, and the reconstruction threshold. In this paper we prove these physics predictions for a broad class of random constraint satisfaction problems. Additionally, we obtain contiguity results that have implications for Bayesian inference tasks, a subject that has received a great deal of interest recently (e.g. Banks et al. 2016, Proc. 29th COLT).

Given complex numbers w1,…,wn, we define the weight w(X) of a set X of 0–1 vectors as the sum of $w_1^{x_1} \cdots w_n^{x_n}$ over all vectors (x1,…,xn) in X. We present an algorithm which, for a set X defined by a system of homogeneous linear equations with at most r variables per equation and at most c equations per variable, computes w(X) within relative error ∊ > 0 in (rc)O(lnn-ln∊) time provided $|w_j| \leq \beta (r \sqrt{c})^{-1}$ for an absolute constant β > 0 and all j = 1,…,n. A similar algorithm is constructed for computing the weight of a linear code over ${\mathbb F}_p$. Applications include counting weighted perfect matchings in hypergraphs, counting weighted graph homomorphisms, computing weight enumerators of linear codes with sparse code generating matrices, and computing the partition functions of the ferromagnetic Potts model at low temperatures and of the hard-core model at high fugacity on biregular bipartite graphs.

A self-avoiding polygon is a lattice polygon consisting of a closed self-avoiding walk on a square lattice. Surprisingly little is known rigorously about the enumeration of self-avoiding polygons, although there are numerous conjectures that are believed to be true and strongly supported by numerical simulations. As an analogous problemto this study, we considermultiple self-avoiding polygons in a confined region as a model for multiple ring polymers in physics. We find rigorous lower and upper bounds for the number ${{p}_{m\times n}}$ of distinct multiple self-avoiding polygons in the $m\,\times \,n$ rectangular grid on the square lattice. For $m\,=\,2,\,{{p}_{2\times n}}\,=\,{{2}^{n-1}}\,-1$. And for integers $m,\,n\,\ge \,3$,

$${{2}^{m+n-3}}\left( \frac{17}{10} \right){{\,}^{\left( m-2 \right)\left( n-2 \right)}}\,\le \,{{p}_{m\times n}}\,\le \,{{2}^{m+n-3}}\left( \frac{31}{16} \right){{\,}^{\left( m-2 \right)\left( n-2 \right)}}.$$

A cutset is a non-empty finite subset of ℤd which is both connected and co-connected. A cutset is odd if its vertex boundary lies in the odd bipartition class of ℤd. Peled [18] suggested that the number of odd cutsets which contain the origin and have n boundary edges may be of order eΘ(n/d) as d → ∞, much smaller than the number of general cutsets, which was shown by Lebowitz and Mazel [15] to be of order dΘ(n/d). In this paper, we verify this by showing that the number of such odd cutsets is (2+o(1))n/2d.

We consider the Widom–Rowlinson model of two types of interacting particles on d-regular graphs. We prove a tight upper bound on the occupancy fraction, the expected fraction of vertices occupied by a particle under a random configuration from the model. The upper bound is achieved uniquely by unions of complete graphs on d + 1 vertices, Kd+1. As a corollary we find that Kd+1 also maximizes the normalized partition function of the Widom–Rowlinson model over the class of d-regular graphs. A special case of this shows that the normalized number of homomorphisms from any d-regular graph G to the graph HWR, a path on three vertices with a loop on each vertex, is maximized by Kd+1. This proves a conjecture of Galvin.

Correlation functions in the (n) models below the critical temperature are considered. Based on Monte Carlo (MC) data, we confirm the fact stated earlier by Engels and Vogt, that the transverse two-plane correlation function of the (4) model for lattice sizes about L = 120 and small external fields h is very well described by a Gaussian approximation. However, we show that fits of not lower quality are provided by certain non-Gaussian approximation. We have also tested larger lattice sizes, up to L = 512. The Fourier-transformed transverse and longitudinal two-point correlation functions have Goldstone mode singularities in the thermodynamic limit at k → 0 and h = +0, i.e., G⊥ (k) ≃ ak–λ⊥ and G‖(k)≃bk–λ‖, respectively. Here a and b are the amplitudes, k = |k| is the magnitude of the wave vector k. The exponents λᚆ, λ‖ and the ratio bM2/a2, where M is the spontaneous magnetization, are universal according to the GFD (grouping of Feynman diagrams) approach. Here we find that the universality follows also from the standard (Gaussian) theory, yielding bM2/a2=(n−1)/16. Our MC estimates of this ratio are 0.06±0.01 for n=2, 0.17±0.01 for n = 4 and 0.498±0.010 for n = 10. According to these and our earlier MC results, the asymptotic behavior and Goldstone mode singularities are not exactly described by the standard theory. This is expected from the GFD theory. We have found appropriate analytic approximations for G⊥(k) and G‖(k), well fitting the simulation data for small k. We have used them to test the Patashinski-Pokrovski relation and have found that it holds approximately.

The singularity of specific heat CV of the three-dimensional Ising model is studied based on Monte Carlo data for lattice sizes L≤1536. Fits of two data sets, one corresponding to certain value of the Binder cumulant and the other — to the maximum of CV, provide consistent values of C0 in the ansatz CV(L)=C0+ALα/ν at large L, if α/ν=0.196(6). However, a direct estimation from our data suggests that α/ν, most probably, has a smaller value (e.g., α/ν= 0.113(30)). Thus, the conventional power-law scaling ansatz can be questioned because of this inconsistency. We have found that the data are well described by certain logarithmic ansatz.