400 research outputs found
Monte Carlo Renormalization of the 3-D Ising model: Analyticity and Convergence
We review the assumptions on which the Monte Carlo renormalization technique
is based, in particular the analyticity of the block spin transformations. On
this basis, we select an optimized Kadanoff blocking rule in combination with
the simulation of a d=3 Ising model with reduced corrections to scaling. This
is achieved by including interactions with second and third neighbors. As a
consequence of the improved analyticity properties, this Monte Carlo
renormalization method yields a fast convergence and a high accuracy. The
results for the critical exponents are y_H=2.481(1) and y_T=1.585(3).Comment: RevTeX, 4 PostScript file
Generalization of the Fortuin-Kasteleyn transformation and its application to quantum spin simulations,
We generalize the Fortuin-Kasteleyn (FK) cluster representation of the
partition function of the Ising model to represent the partition function of
quantum spin models with an arbitrary spin magnitude in arbitrary dimensions.
This generalized representation enables us to develop a new cluster algorithm
for the simulation of quantum spin systems by the worldline Monte Carlo method.
Because the Swendsen-Wang algorithm is based on the FK representation, the new
cluster algorithm naturally includes it as a special case. As well as the
general description of the new representation, we present an illustration of
our new algorithm for some special interesting cases: the Ising model, the
antiferromagnetic Heisenberg model with , and a general Heisenberg model.
The new algorithm is applicable to models with any range of the exchange
interaction, any lattice geometry, and any dimensions.Comment: 46 pages, 10 figures, to appear in J.Stat.Phy
Graphical representations and cluster algorithms for critical points with fields
A two-replica graphical representation and associated cluster algorithm is
described that is applicable to ferromagnetic Ising systems with arbitrary
fields. Critical points are associated with the percolation threshold of the
graphical representation. Results from numerical simulations of the Ising model
in a staggered field are presented. The dynamic exponent for the algorithm is
measured to be less than 0.5.Comment: Revtex, 12 pages with 2 figure
Transition Matrix Monte Carlo Reweighting and Dynamics
We study an induced dynamics in the space of energy of single-spin-flip Monte
Carlo algorithm. The method gives an efficient reweighting technique. This
dynamics is shown to have relaxation times proportional to the specific heat.
Thus, it is plausible for a logarithmic factor in the correlation time of the
standard 2D Ising local dynamics.Comment: RevTeX, 5 pages, 3 figure
The Block Spin Renormalization Group Approach and Two-Dimensional Quantum Gravity
A block spin renormalization group approach is proposed for the dynamical
triangulation formulation of two-dimensional quantum gravity. The idea is to
update link flips on the block lattice in response to link flips on the
original lattice. Just as the connectivity of the original lattice is meant to
be a lattice representation of the metric, the block links are determined in
such a way that the connectivity of the block lattice represents a block
metric. As an illustration, this approach is applied to the Ising model coupled
to two-dimensional quantum gravity. The correct critical coupling is
reproduced, but the critical exponent is obscured by unusually large finite
size effects.Comment: 10 page
Reconstructing the Density of States by History-Dependent Metadynamics
We present a novel method for the calculation of the energy density of states
D(E) for systems described by classical statistical mechanics. The method
builds on an extension of a recently proposed strategy that allows the free
energy profile of a canonical system to be recovered within a pre-assigned
accuracy,[A. Laio and M. Parrinello, PNAS 2002]. The method allows a good
control over the error on the recovered system entropy. This fact is exploited
to obtain D(E) more efficiently by combining measurements at different
temperatures. The accuracy and efficiency of the method are tested for the
two-dimensional Ising model (up to size 50x50) by comparison with both exact
results and previous studies. This method is a general one and should be
applicable to more realistic model systems
An Almost Perfect Quantum Lattice Action for Low-energy SU(2) Gluodynamics
We study various representations of infrared effective theory of SU(2)
Gluodynamics as a (quantum) perfect lattice action. In particular we derive a
monopole action and a string model of hadrons from SU(2) Gluodynamics. These
are lattice actions which give almost cut-off independent physical quantities
even on coarse lattices. The monopole action is determined by numerical
simulations in the infrared region of SU(2) Gluodynamics. The string model of
hadrons is derived from the monopole action by using BKT transformation. We
illustrate the method and evaluate physical quantities such as the string
tension and the mass of the lowest state of the glueball analytically using the
string model of hadrons. It turns out that the classical results in the string
model is near to the one in quantum SU(2) Gluodynamics.Comment: 39 pages, 10 figure
Cluster Algorithm for a Solid-On-Solid Model with Constraints
We adapt the VMR (valleys-to-mountains reflections) algorithm, originally
devised by us for simulations of SOS models, to the BCSOS model. It is the
first time that a cluster algorithm is used for a model with constraints. The
performance of this new algorithm is studied in detail in both phases of the
model, including a finite size scaling analysis of the autocorrelations.Comment: 10 pages, 3 figures appended as ps-file
- …