We study two types of generalized Baxter-Wu models, by means of transfer-matrix and Monte Carlo techniques. The first generalization allows for different couplings in the up-and down-triangles, and the second generalization is to a q-state spin model with three-spin interactions. Both generalizations lead to self-dual models, so that the probable locations of the phase transitions follow. Our numerical analysis confirms that phase transitions occur at the self-dual points. For both generalizations of the Baxter-Wu model, the phase transitions appear to be discontinuous.

We formulate a single-cluster Monte Carlo algorithm for the simulation of the random-cluster model. This algorithm is a generalization of the Wolff single-cluster method for the q-state Potts model to non-integer values q>1. Its results for static quantities are in a satisfactory agreement with those of the existing Swendsen-Wang-Chayes-Machta (SWCM) algorithm, which involves a full cluster decomposition of random-cluster configurations. We explore the critical dynamics of this algorithm for several two-dimensional Potts and random-cluster models. For integer q, the single-cluster algorithm can be reduced to the Wolff algorithm, for which case we find that the autocorrelation functions decay almost purely exponentially, with dynamic exponents zexp=0.07 (1), 0.521 (7), and 1.007 (9) for q=2,3, and 4 respectively. For non-integer q, the dynamical behavior of the single-cluster algorithm appears to be very dissimilar to that of the SWCM algorithm. For large critical systems, the autocorrelation function displays a range of power-law behavior as a function of time. The dynamic exponents are relatively large. We provide an explanation for this peculiar dynamic behavior.

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≤n≤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.

In this article, the effects of irreversible detachments of particles in totally asymmetric simple exclusion processes (TASEPs) with extended particles which occupy more than one lattice site, are investigated. First, an approximate meanfield theory is used to calculate phase diagrams and density profiles. The results show that the detachment and the size of the particles have distinct effects on the stationary phases in the two sublattices divided by the detachment, especially in the mc/hd and the hd/hd phase. Here, symbols"mc", "hd", and"ld" are initials of maximal current, high density, and low density, respectively, and"mc/hd"represents the stationary state that the left sublattice is in the maximum-current (mc) phase while the right sublattice is in the high density (hd) phase. When the detachment rate is very large, there are four stationary phases, including ld/ld, ld/hd, mc/ld, and mc/hd phases. When the value of the detachment rate is in the middle range, the hd/hd phase occurs, and hence there exist five stationary phases. When the rate is very small, the mc/hd phase disappears, and there are only four phases again. These theoretical results qualitatively and quantitatively agree with computer Monte Carlo simulations especially in the case of large value of the detachment rate.

We derive the scaling dimension associated with crossing bonds in the random-cluster representation of the two-dimensional Potts model, by means of a mapping on the Coulomb gas. The scaling field associated with crossing bonds appears to be irrelevant, on the critical as well as on the tricritical branch. The latter result stands in a remarkable contrast with the existing result for the tricritical O(n) model that crossing bonds are relevant. In order to obtain independent confirmation of the Coulomb gas result for the crossing-bond exponent, we perform a finite-size-scaling analysis based on numerical transfer-matrix calculations.

We introduce several infinite families of new critical exponents for the random-cluster model and present scaling arguments relating them to the k-arm exponents. We then present Monte Carlo simulations confirming these predictions. These new exponents provide a convenient way to determine k-arm exponents from Monte Carlo simulations. An understanding of these exponents also leads to a radically improved implementation of the Sweeny Monte Carlo algorithm. In addition, our Monte Carlo data allow us to conjecture an exact expression for the shortest-path fractal dimension dmin in two dimensions: dmin?=(g+2)(g+18)/(32g) where g is the Coulomb-gas coupling, related to the cluster fugacity q via q=2+2 cos(gπ/2) with 2≤g≤4.

We present a Markov-chain Monte Carlo algorithm of worm type that correctly simulates the fully-packed loop model with n = 1 on the honeycomb lattice, and we prove that it is ergodic and has uniform stationary distribution. The honeycomblattice fully-packed loop model with n=1 is equivalent to the zero-temperature triangular-lattice antiferromagnetic Ising model, which is fully frustrated and notoriously difficult to simulate. We test this worm algorithm numerically and estimate the dynamic exponent zexp=0.515(8). We also measure several static quantities of interest, including loop-length and face-size moments. It appears numerically that the face-size moments are governed by the magnetic dimension for percolation.

We present Monte Carlo simulations of the spanning-forest model (q→0 limit of the ferromagnetic Potts model) in spatial dimensions d=3,4,5. We show that, in contrast to the two-dimensional case, the model has a ferromagnetic second-order phase transition at a finite positive value wc. We present numerical estimates of wc and of the thermal and magnetic critical exponents. We conjecture that the upper critical dimension is 6.

We investigate several tricritical models on the square lattice by means of Monte Carlo simulations. These include the Blume-Capel model, Baxter's hard-square model, and the q=1,3, and 4 Potts models with vacancies. We use a combination of the Wolff and geometric cluster methods, which conserves the total number of vacancies or lattice-gas particles and suppresses critical slowing down. Several quantities are sampled, such as the specific heat C and the structure factor Cs, which accounts for the large-scale spatial inhomogeneity of the energy fluctuations. We find that the constraint strongly modifies some of the critical singularities. For instance, the specific heat C reaches a finite value at tricriticality, while Cs remains divergent as in the unconstrained system. We are able to explain these observed constrained phenomena on the basis of the Fisher renormalization mechanism generalized to include a subleading relevant thermal scaling field. In this context, we find that, under the constraint, the leading thermal exponent yt1 is renormalized to 2−yt1, while the subleading exponent yt2 remains unchanged.

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 bondpercolation 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 Xt2=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 such a logarithmic term appears to be small or absent. The amplitude of the power-law correction associated with Xt2=4 is found to be dependent on the orientation of the lattice with respect to the cylindrical geometry of the finite systems.