Critical Behavior of the Ferromagnetic Ising Model on a Sierpinski Carpet: Monte Carlo Renormalization Group Study

We perform a Monte Carlo Renormalization Group analysis of the critical behavior of the ferromagnetic Ising model on a Sierpi\'nski fractal with Hausdorff dimension $d_f\simeq 1.8928$. This method is shown to be relevant to the calculation of the critical temperature $T_c$ and the magnetic eigen-exponent $y_h$ on such structures. On the other hand, scaling corrections hinder the calculation of the temperature eigen-exponent $y_t$. At last, the results are shown to be consistent with a finite size scaling analysis.Comment: 16 pages, 7 figure

Monopole Condensation in Lattice SU(2) QCD

This is the short review of Monte-Carlo studies of quark confinement in lattice QCD. After abelian projections both in the maximally abelian and Polyakov gauges, it is seen that the monopole part alone is responsible for confinement. A block spin transformation on the dual lattice suggests that lattice $SU(2)$ QCD is always ( for all $\beta$) in the monopole condensed phase and so in the confinement phase in the infinite volume limit.Comment: Contribution to Confinement '95, March 1995, Osaka, Japan. Names of figure files are corrected. 8 page uuencoded latex file and 10 ps figure

On generalized cluster algorithms for frustrated spin models

Standard Monte Carlo cluster algorithms have proven to be very effective for many different spin models, however they fail for frustrated spin systems. Recently a generalized cluster algorithm was introduced that works extremely well for the fully frustrated Ising model on a square lattice, by placing bonds between sites based on information from plaquettes rather than links of the lattice. Here we study some properties of this algorithm and some variants of it. We introduce a practical methodology for constructing a generalized cluster algorithm for a given spin model, and investigate apply this method to some other frustrated Ising models. We find that such algorithms work well for simple fully frustrated Ising models in two dimensions, but appear to work poorly or not at all for more complex models such as spin glasses.Comment: 34 pages in RevTeX. No figures included. A compressed postscript file for the paper with figures can be obtained via anonymous ftp to minerva.npac.syr.edu in users/paulc/papers/SCCS-527.ps.Z. Syracuse University NPAC technical report SCCS-52

Loop Algorithms for Asymmetric Hamiltonians

Generalized rules for building and flipping clusters in the quantum Monte Carlo loop algorithm are presented for the XXZ-model in a uniform magnetic field along the Z-axis. As is demonstrated for the Heisenberg antiferromagnet it is possible from these rules to select a new algorithm which performs significantly better than the standard loop algorithm in strong magnetic fields at low temperatures.Comment: Replaced measurement of helicity modulus at H=2J with a measurement at H=3.95J + other small changes in the section on numerical result

The Percolation Signature of the Spin Glass Transition

Magnetic ordering at low temperature for Ising ferromagnets manifests itself within the associated Fortuin-Kasteleyn (FK) random cluster representation as the occurrence of a single positive density percolating network. In this paper we investigate the percolation signature for Ising spin glass ordering -- both in short-range (EA) and infinite-range (SK) models -- within a two-replica FK representation and also within the different Chayes-Machta-Redner two-replica graphical representation. Based on numerical studies of the $\pm J$ EA model in three dimensions and on rigorous results for the SK model, we conclude that the spin glass transition corresponds to the appearance of {\it two} percolating clusters of {\it unequal} densities.Comment: 13 pages, 6 figure

A Cluster Method for the Ashkin--Teller Model

A cluster Monte Carlo algorithm for the Ashkin-Teller (AT) model is constructed according to the guidelines of a general scheme for such algorithms. Its dynamical behaviour is tested for the square lattice AT model. We perform simulations on the line of critical points along which the exponents vary continuously, and find that critical slowing down is significantly reduced. We find continuous variation of the dynamical exponent $z$ along the line, following the variation of the ratio $\alpha/\nu$, in a manner which satisfies the Li-Sokal bound $z_{cluster}\geq\alpha/\nu$, that was so far proved only for Potts models.Comment: 18 pages, Revtex, figures include

Scaling functions for O(4) in three dimensions

Monte Carlo simulation using a cluster algorithm is used to compute the scaling part of the free energy for a three dimensional O(4) spin model. The results are relevant for analysis of lattice studies of high temperature QCD.Comment: 12 pages, 6 figures, uses epsf.st

Three-dimensional Ising model in the fixed-magnetization ensemble: a Monte Carlo study

We study the three-dimensional Ising model at the critical point in the fixed-magnetization ensemble, by means of the recently developed geometric cluster Monte Carlo algorithm. We define a magnetic-field-like quantity in terms of microscopic spin-up and spin-down probabilities in a given configuration of neighbors. In the thermodynamic limit, the relation between this field and the magnetization reduces to the canonical relation M(h). However, for finite systems, the relation is different. We establish a close connection between this relation and the probability distribution of the magnetization of a finite-size system in the canonical ensemble.Comment: 8 pages, 2 Postscript figures, uses RevTe

Nonperturbative bound on high multiplicity cross sections in phi^4_3 from lattice simulation

We have looked for evidence of large cross sections at large multiplicities in weakly coupled scalar field theory in three dimensions. We use spectral function sum rules to derive bounds on total cross sections where the sum can be expresed in terms of a quantity which can be measured by Monte Carlo simulation in Euclidean space. We find that high multiplicity cross sections remain small for energies and multiplicities for which large effects had been suggested.Comment: 23 pages, revtex, seven eps figures revised version: typos corrected, some rewriting of discusion, same resul

Cluster Monte Carlo Simulations of the Nematic--Isotropic Transition

We report the results of simulations of the Lebwohl-Lasher model of the nematic-isotropic transition using a new cluster Monte Carlo algorithm. The algorithm is a modification of the Wolff algorithm for spin systems, and greatly reduces critical slowing down. We calculate the free energy in the neighborhood of the transition for systems up to linear size 70. We find a double well structure with a barrier that grows with increasing system size, obeying finite size scaling for systems of size greater than 35. We thus obtain an estimate of the value of the transition temperature in the thermodynamic limit.Comment: 4 figure