We study the computational complexity of approximating the partition function of a q-state spin system with an external field. There are just three possible levels of computational difficulty, depending on the interaction strengths between adjacent spins: (i) efficiently exactly computable, (ii) equivalent to the ferromagnetic Ising model, and (iii) equivalent to the antiferromagnetic Ising model. Thus, every nontrivial q-state spin system, irrespective of the number q of spins, is computationally equivalent to one of two fundamental two-state spin systems.computational complexity | partition function | spin system L et Q = f1,2, . . . , qg denote a set of spins. A q-state spin system is specified by a nonnegative real symmetric interaction matrix A ∈ ðR ≥0 Þ q×q . The entries ða ij Þ of A represent "interaction strengths" between spins in Q. An instance of such a spin system is a graph G = ðV , EÞ, where V is a set of vertices (sites) and E is a set of edges (bonds) together with a collection h = fh w : Q → R ≥0 jw ∈ V g of functions, representing the action of an external field. The function h w represents the effect of the field on vertex w. The partition function of the system is then a weighted sum over configurations σ : V → Q.The above setting encompasses all spin systems with uniform interactions between pairs of sites and includes many familiar models. For example, the interaction matricescapture the Ising, independent set (hard-core), three-state Potts, and Widom-Rowlinson models, respectively. In the case of the Ising and Potts models, λ > 1 corresponds to a ferromagnetic system and λ < 1 to an antiferromagnetic one.
The Problem and What Is KnownWe study the following computational problem. Fix an interaction matrix A. Given an instance consisting of a graph G and field h, approximately evaluate the partition function Z A ðG, hÞ. An important point to note is that the interaction matrix A does not form part of the problem instance. By fixing A, we fix a particular model, say the hardcore model or the q-state ferromagnetic Potts model. We then ask, for that model, What is the computational complexity of approximately evaluating the partition function given the input pair ðG, hÞ?We are interested in determining how that complexity depends on A. The purpose of this paper is to map the space of interaction matrices, delineating which case arises for each interaction matrix A. Before stating our result, we fill in some details about the three possibilities above (and approximation complexity in general), we say a little bit more about the result of ref. 1, and we present some matrix preliminaries.The first possibility, that the partition function can be evaluated exactly in polynomial time, is straightforward. Unfortunately, as we shall see, it rarely arises. Let us turn to the second possibility. It will be helpful to introduce the problem #BIS, which plays an important role in the complexity theory of approximate counting problems (2). #BIS is the problem of counting the independent sets of a bipartite g...