Transfer matrix computation of critical polynomials for two-dimensional Potts models

TL;DR

This paper introduces a probabilistic transfer matrix method to compute critical polynomials for 2D Potts models on various lattices, achieving high-precision critical point estimates and exploring phase diagrams.

Contribution

It presents a new probabilistic definition of the critical polynomial enabling larger basis computations via transfer matrix, improving accuracy over previous contraction-deletion methods.

Findings

Accurate critical temperature estimates matching or surpassing Monte Carlo results.

Extended computational capacity to bases with up to 243 edges.

Detailed analysis of phase diagram structure in the antiferromagnetic region.

Abstract

In our previous work we have shown that critical manifolds of the q-state Potts model can be studied by means of a graph polynomial P_B(q,v), henceforth referred to as the critical polynomial. This polynomial may be defined on any periodic two-dimensional lattice. It depends on a finite subgraph B, called the basis, and the manner in which B is tiled to construct the lattice. The real roots v = e^K - 1 of P_B(q,v) either give the exact critical points for the lattice, or provide approximations that, in principle, can be made arbitrarily accurate by increasing the size of B in an appropriate way. In earlier work, P_B(q,v) was defined by a contraction-deletion identity, similar to that satisfied by the Tutte polynomial. Here, we give a probabilistic definition of P_B(q,v), which facilitates its computation, using the transfer matrix, on much larger B than was previously possible. We present results for the critical polynomial on the (4,8^2), kagome, and (3,12^2) lattices for bases of up to respectively 96, 162, and 243 edges, compared to the limit of 36 edges with contraction-deletion. We discuss in detail the role of the symmetries and the embedding of B. The critical temperatures v_c obtained for ferromagnetic (v>0) Potts models are at least as precise as the best available results from Monte Carlo simulations or series expansions. For instance, with q=3 we obtain v_c(4,8^2) = 3.742 489(4), v_c(kagome) = 1.876 459 7(2), and v_c(3,12^2) = 5.033 078 49(4), the precision being comparable or superior to the best simulation results. More generally, we trace the critical manifolds in the real (q,v) plane and discuss the intricate structure of the phase diagram in the antiferromagnetic (v<0) region.