Ising Spins on Thin Graphs

The Ising model on ``thin'' graphs (standard Feynman diagrams) displays several interesting properties. For ferromagnetic couplings there is a mean field phase transition at the corresponding Bethe lattice transition point. For antiferromagnetic couplings the replica trick gives some evidence for a spin glass phase. In this paper we investigate both the ferromagnetic and antiferromagnetic models with the aid of simulations. We confirm the Bethe lattice values of the critical points for the ferromagnetic model on $\phi^3$ and $\phi^4$ graphs and examine the putative spin glass phase in the antiferromagnetic model by looking at the overlap between replicas in a quenched ensemble of graphs. We also compare the Ising results with those for higher state Potts models and Ising models on ``fat'' graphs, such as those used in 2D gravity simulations.Comment: LaTeX 13 pages + 9 postscript figures, COLO-HEP-340, LPTHE-Orsay-94-6

Spin Glasses on Thin Graphs

In a recent paper we found strong evidence from simulations that the Isingantiferromagnet on ``thin'' random graphs - Feynman diagrams - displayed amean-field spin glass transition. The intrinsic interest of considering such random graphs is that they give mean field results without long range interactions or the drawbacks, arising from boundary problems, of the Bethe lattice. In this paper we reprise the saddle point calculations for the Ising and Potts ferromagnet, antiferromagnet and spin glass on Feynman diagrams. We use standard results from bifurcation theory that enable us to treat an arbitrary number of replicas and any quenched bond distribution. We note the agreement between the ferromagnetic and spin glass transition temperatures thus calculated and those derived by analogy with the Bethe lattice, or in previous replica calculations. We then investigate numerically spin glasses with a plus or minus J bond distribution for the Ising and Q=3,4,10,50 state Potts models, paying particular attention to the independence of the spin glass transition from the fraction of positive and negative bonds in the Ising case and the qualitative form of the overlap distribution in all the models. The parallels with infinite range spin glass models in both the analytical calculations and simulations are pointed out.Comment: 13 pages of LaTex and 11 postscript figures bundled together with uufiles. Discussion of first order transitions for three or more replicas included and similarity to Ising replica magnet pointed out. Some additional reference

Density of states, Potts zeros, and Fisher zeros of the Q-state Potts model for continuous Q

The Q-state Potts model can be extended to noninteger and even complex Q in the FK representation. In the FK representation the partition function,Z(Q,a), is a polynomial in Q and v=a-1(a=e^-T) and the coefficients of this polynomial,Phi(b,c), are the number of graphs on the lattice consisting of b bonds and c connected clusters. We introduce the random-cluster transfer matrix to compute Phi exactly on finite square lattices. Given the FK representation of the partition function we begin by studying the critical Potts model Z_{CP}=Z(Q,a_c), where a_c=1+sqrt{Q}. We find a set of zeros in the complex w=sqrt{Q} plane that map to the Beraha numbers for real positive Q. We also identify tilde{Q}_c(L), the value of Q for a lattice of width L above which the locus of zeros in the complex p=v/sqrt{Q} plane lies on the unit circle. We find that 1/tilde{Q}_c->0 as 1/L->0. We then study zeros of the AF Potts model in the complex Q plane and determine Q_c(a), the largest value of Q for a fixed value of a below which there is AF order. We find excellent agreement with Q_c=(1-a)(a+3). We also investigate the locus of zeros of the FM Potts model in the complex Q plane and confirm that Q_c=(a-1)^2. We show that the edge singularity in the complex Q plane approaches Q_c as Q_c(L)~Q_c+AL^-y_q, and determine the scaling exponent y_q. Finally, by finite size scaling of the Fisher zeros near the AF critical point we determine the thermal exponent y_t as a function of Q in the range 2<Q<3. We find that y_t is a smooth function of Q and is well fit by y_t=(1+Au+Bu^2)/(C+Du) where u=u(Q). For Q=3 we find y_t~0.6; however if we include lattices up to L=12 we find y_t~0.50.Comment: to appear in Physical Review

Counting Complex Disordered States by Efficient Pattern Matching: Chromatic Polynomials and Potts Partition Functions

Counting problems, determining the number of possible states of a large system under certain constraints, play an important role in many areas of science. They naturally arise for complex disordered systems in physics and chemistry, in mathematical graph theory, and in computer science. Counting problems, however, are among the hardest problems to access computationally. Here, we suggest a novel method to access a benchmark counting problem, finding chromatic polynomials of graphs. We develop a vertex-oriented symbolic pattern matching algorithm that exploits the equivalence between the chromatic polynomial and the zero-temperature partition function of the Potts antiferromagnet on the same graph. Implementing this bottom-up algorithm using appropriate computer algebra, the new method outperforms standard top-down methods by several orders of magnitude, already for moderately sized graphs. As a first application, we compute chromatic polynomials of samples of the simple cubic lattice, for the first time computationally accessing three-dimensional lattices of physical relevance. The method offers straightforward generalizations to several other counting problems.Comment: 7 pages, 4 figure

Exact Potts Model Partition Function on Strips of the Triangular Lattice

In this paper we present exact calculations of the partition function $Z$ of the $q$-state Potts model and its generalization to real $q$, for arbitrary temperature on $n$-vertex strip graphs, of width $L_y=2$ and arbitrary length, of the triangular lattice with free, cyclic, and M\"obius longitudinal boundary conditions. These partition functions are equivalent to Tutte/Whitney polynomials for these graphs. The free energy is calculated exactly for the infinite-length limit of the graphs, and the thermodynamics is discussed. Considering the full generalization to arbitrary complex $q$ and temperature, we determine the singular locus ${\cal B}$ in the corresponding ${\mathbb C}^2$ space, arising as the accumulation set of partition function zeros as $n \to \infty$. In particular, we study the connection with the T=0 limit of the Potts antiferromagnet where ${\cal B}$ reduces to the accumulation set of chromatic zeros. Comparisons are made with our previous exact calculation of Potts model partition functions for the corresponding strips of the square lattice. Our present calculations yield, as special cases, several quantities of graph-theoretic interest.Comment: 43 pages, latex, 24 postscript figures, Physica A, in pres

Exact Potts Model Partition Functions on Strips of the Honeycomb Lattice

We present exact calculations of the partition function of the $q$-state Potts model on (i) open, (ii) cyclic, and (iii) M\"obius strips of the honeycomb (brick) lattice of width $L_y=2$ and arbitrarily great length. In the infinite-length limit the thermodynamic properties are discussed. The continuous locus of singularities of the free energy is determined in the $q$ plane for fixed temperature and in the complex temperature plane for fixed $q$ values. We also give exact calculations of the zero-temperature partition function (chromatic polynomial) and $W(q)$, the exponent of the ground-state entropy, for the Potts antiferromagnet for honeycomb strips of type (iv) $L_y=3$, cyclic, (v) $L_y=3$, M\"obius, (vi) $L_y=4$, cylindrical, and (vii) $L_y=4$, open. In the infinite-length limit we calculate $W(q)$ and determine the continuous locus of points where it is nonanalytic. We show that our exact calculation of the entropy for the $L_y=4$ strip with cylindrical boundary conditions provides an extremely accurate approximation, to a few parts in $10^5$ for moderate $q$ values, to the entropy for the full 2D honeycomb lattice (where the latter is determined by Monte Carlo measurements since no exact analytic form is known).Comment: 48 pages, latex, with encapsulated postscript figure

On the non-ergodicity of the Swendsen-Wang-Kotecky algorithm on the kagome lattice

We study the properties of the Wang-Swendsen-Kotecky cluster Monte Carlo algorithm for simulating the 3-state kagome-lattice Potts antiferromagnet at zero temperature. We prove that this algorithm is not ergodic for symmetric subsets of the kagome lattice with fully periodic boundary conditions: given an initial configuration, not all configurations are accessible via Monte Carlo steps. The same conclusion holds for single-site dynamics.Comment: Latex2e. 22 pages. Contains 11 figures using pstricks package. Uses iopart.sty. Final version accepted in journa