We study a percolation problem based on critical loop configurations of the O(n) loop model on the honeycomb lattice. We define dual clusters as groups of sites on the dual triangular lattice that are not separated by a loop, and investigate the bond-percolation properties of these dual clusters. The universal properties at the percolation threshold are argued to match those of Kasteleyn-Fortuin random clusters in the critical Potts model. This relation is checked numerically by means of cluster simulations of several O(n) models in the range 1<or=n<or=2. The simulation results include the percolation threshold for several values of n, as well as the universal exponents associated with bond dilution and the size distribution of the diluted clusters at the percolation threshold. Our numerical results for the exponents are in agreement with existing Coulomb-gas results for the random-cluster model, which confirms the relation between both models. We discuss the renormalization flow of the bond-dilution parameter p as a function of n, and provide an expression that accurately describes a line of unstable fixed points as a function of n, corresponding with the percolation threshold. Furthermore, the renormalization scenario indicates the existence, in a p versus n diagram, of another line of fixed points at p=1, which is stable with respect to p.
We study a percolation problem based on critical loop configurations of the O(n) loop model on the honeycomb lattice. We define dual clusters as groups of sites on the dual triangular lattice that are not separated by a loop, and investigate the bond-percolation properties of these dual clusters. The universal properties at the percolation threshold are argued to match those of Kasteleyn-Fortuin random clusters in the critical Potts model. This relation is checked numerically by means of cluster simulations of several O(n) models in the range 1<or=n<or=2. The simulation results include the percolation threshold for several values of n, as well as the universal exponents associated with bond dilution and the size distribution of the diluted clusters at the percolation threshold. Our numerical results for the exponents are in agreement with existing Coulomb-gas results for the random-cluster model, which confirms the relation between both models. We discuss the renormalization flow of the bond-dilution parameter p as a function of n, and provide an expression that accurately describes a line of unstable fixed points as a function of n, corresponding with the percolation threshold. Furthermore, the renormalization scenario indicates the existence, in a p versus n diagram, of another line of fixed points at p=1, which is stable with respect to p.
We study bond percolation in two dimensions between random site variables having one out of q colors, using transfer-matrix and Monte Carlo techniques. We determine the percolation threshold as a function of the Potts temperature T in the disordered Potts range ${T}_{c}\ensuremath{\leqslant}T<\ensuremath{\infty}$ for several q-state Potts Hamiltonians. For high T, these transitions fit, irrespective of q, in the universality class of the ordinary percolation transitions. However, for $T\ensuremath{\downarrow}{T}_{c}$, q-dependent crossover phenomena appear. The topology of the phase diagram changes in a qualitative sense at q=2. For q<2 the Potts critical state appears to enhance percolation, for q>2 it appears to suppress it. Remarkably, for q=2 the percolation line coincides with the only flow line extending to T>Tc from the critical fixed point associated with Potts clusters.
The author proposes a new type of cluster in addition to those of Reynolds, Klein and Stanley for the square lattice and another in addition to those of Chao for the triangular lattice for a real-space renormalisation group approach for studying critical phenomena. Thus each lattice is associated with two types of clusters-the longitudinal and the diagonal. The results for the bond percolation and Ising fixed points obtained from the diagonal square and honeycomb clusters agree with the exact values to within 6-11% error; whereas those obtained from the longitudinal triangular and honeycomb clusters agree with the exact values to within 2-4.5% error. The critical exponents corresponding to the percolation and the Ising problems also show reasonably good results: in particular those calculated by using the longitudinal triangular cluster agree with the exact values to within 6.7% and 3.2% error respectively.
We study the O(N) loop model on the Honeycomb lattice with real value N≥1 by means of a cluster algorithm. The formulation of the algorithm is based on the equivalence of the O(N) loop model and the low-temperature graphical representation of a N-color Ashkin-Teller model on the triangular lattice. The latter model with integer N can be simulated by means of an embedding Swendsen-Wang-type cluster method. By taking into account the symmetry among loops of different colors, we develop another version of the Swendsen-Wang-type method. This version allows the number of colors N to take any real value N≥1. As an application, we investigate the N=1.25,1.50,1.75, and 2 loop model at criticality. The determined values of various critical exponents are in excellent agreement with theoretical predictions. In particular, from quantities associated with half of the loops, we determine some critical exponents that corresponds to those for the tricritical q=N2 Potts model but have not been observed yet. Dynamic scaling behavior of the algorithm is also analyzed. The numerical data strongly suggest that our cluster algorithm {\it hardly} suffers from critical slowing down.
A critical dilute O(n) model on the kagome lattice is investigated analytically and numerically. We employ a number of exact equivalences which, in a few steps, link the critical O(n) spin model on the kagome lattice to the exactly solvable critical q-state Potts model on the honeycomb lattice with q=(n+1)2. The intermediate steps involve the random-cluster model on the honeycomb lattice and a fully packed loop model with loop weight n'=sqrt(q) and a dilute loop model with loop weight n , both on the kagome lattice. This mapping enables the determination of a branch of critical points of the dilute O(n) model, as well as some of its critical properties. These properties differ from those of the generic O(n) critical points. For n=0, our model reproduces the known universal properties of the theta point describing the collapse of a polymer. For n not equal 0 it displays a line of multicritical points, with the same universal behavior as a branch of critical points that was found earlier in a dilute O(n) model on the square lattice. These findings are supported by a finite-size-scaling analysis in combination with transfer-matrix calculations.
In the preceding paper, one of us (F. Y. Wu) considered the Potts model and bond and site percolation on two general classes of two-dimensional lattices, the triangular-type and kagome-type lattices, and obtained closed-form expressions for the critical frontier with applications to various lattice models. For the triangular-type lattices Wu's result is exact, and for the kagome-type lattices Wu's expression is under a homogeneity assumption. The purpose of the present paper is twofold: First, an essential step in Wu's analysis is the derivation of lattice-dependent constants A,B,C for various lattice models, a process which can be tedious. We present here a derivation of these constants for subnet networks using a computer algorithm. Second, by means of a finite-size scaling analysis based on numerical transfer matrix calculations, we deduce critical properties and critical thresholds of various models and assess the accuracy of the homogeneity assumption. Specifically, we analyze the q -state Potts model and the bond percolation on the 3-12 and kagome-type subnet lattices (n×n):(n×n) , n≤4 , for which the exact solution is not known. Our numerical determination of critical properties such as conformal anomaly and magnetic correlation length verifies that the universality principle holds. To calibrate the accuracy of the finite-size procedure, we apply the same numerical analysis to models for which the exact critical frontiers are known. The comparison of numerical and exact results shows that our numerical values are correct within errors of our finite-size analysis, which correspond to 7 or 8 significant digits. This in turn infers that the homogeneity assumption determines critical frontiers with an accuracy of 5 decimal places or higher. Finally, we also obtained the exact percolation thresholds for site percolation on kagome-type subnet lattices (1×1):(n×n) for 1≤n≤6 .
We investigate bond- and site-percolation models on several two-dimensional lattices numerically, by means of transfer-matrix calculations and Monte Carlo simulations. The lattices include the square, triangular, honeycomb kagome, and diced lattices with nearest-neighbor bonds, and the square lattice with nearest- and next-nearest-neighbor bonds. Results are presented for the bond-percolation thresholds of the kagome and diced lattices, and the site-percolation thresholds of the square, honeycomb, and diced lattices. We also include the bond- and site-percolation thresholds for the square lattice with nearest- and next-nearest-neighbor bonds. We find that corrections to scaling behave according to the second temperature dimension X_{t2}=4 predicted by the Coulomb gas theory and the theory of conformal invariance. In several cases there is evidence for an additional term with the same exponent, but modified by a logarithmic factor. Only for the site-percolation problem on the triangular lattice does such a logarithmic term appear to be small or absent. The amplitude of the power-law correction associated with X_{t2}=4 is found to be dependent on the orientation of the lattice with respect to the cylindrical geometry of the finite systems.