2,594 research outputs found
Monte Carlo computation of correlation times of independent relaxation modes at criticality
We investigate aspects of universality of Glauber critical dynamics in two
dimensions. We compute the critical exponent and numerically corroborate
its universality for three different models in the static Ising universality
class and for five independent relaxation modes. We also present evidence for
universality of amplitude ratios, which shows that, as far as dynamic behavior
is concerned, each model in a given universality class is characterized by a
single non-universal metric factor which determines the overall time scale.
This paper also discusses in detail the variational and projection methods that
are used to compute relaxation times with high accuracy
Universal Dynamics of Independent Critical Relaxation Modes
Scaling behavior is studied of several dominant eigenvalues of spectra of
Markov matrices and the associated correlation times governing critical slowing
down in models in the universality class of the two-dimensional Ising model. A
scheme is developed to optimize variational approximants of progressively
rapid, independent relaxation modes. These approximants are used to reduce the
variance of results obtained by means of an adaptation of a quantum Monte Carlo
method to compute eigenvalues subject to errors predominantly of statistical
nature. The resulting spectra and correlation times are found to be universal
up to a single, non-universal time scale for each model
Transfer-matrix approach to the three-dimensional bond percolation: An application of Novotny's formalism
A transfer-matrix simulation scheme for the three-dimensional (d=3) bond
percolation is presented. Our scheme is based on Novotny's transfer-matrix
formalism, which enables us to consider arbitrary (integral) number of sites N
constituting a unit of the transfer-matrix slice even for d=3. Such an
arbitrariness allows us to perform systematic finite-size-scaling analysis of
the criticality at the percolation threshold. Diagonalizing the transfer matrix
for N =4,5,...,10, we obtain an estimate for the correlation-length critical
exponent nu = 0.81(5)
Optimization of ground and excited state wavefunctions and van der Waals clusters
A quantum Monte Carlo method is introduced to optimize excited state trial
wavefunctions. The method is applied in a correlation function Monte Carlo
calculation to compute ground and excited state energies of bosonic van der
Waals clusters of upto seven particles. The calculations are performed using
trial wavefunctions with general three-body correlations
Numerical Studies of the Two Dimensional XY Model with Symmetry Breaking Fields
We present results of numerical studies of the two dimensional XY model with
four and eight fold symmetry breaking fields. This model has recently been
shown to describe hydrogen induced reconstruction on the W(100) surface. Based
on mean-field and renormalization group arguments,we first show how the
interplay between the anisotropy fields can give rise to different phase
transitions in the model. When the fields are compatible with each other there
is a continuous phase transition when the fourth order field is varied from
negative to positive values. This transition becomes discontinuous at low
temperatures. These two regimes are separated by a multicritical point. In the
case of competing four and eight fold fields, the first order transition at low
temperatures opens up into two Ising transitions. We then use numerical methods
to accurately locate the position of the multicritical point, and to verify the
nature of the transitions. The different techniques used include Monte Carlo
histogram methods combined with finite size scaling analysis, the real space
Monte Carlo Renormalization Group method, and the Monte Carlo Transfer Matrix
method. Our numerical results are in good agreement with the theoretical
arguments.Comment: 29 pages, HU-TFT-94-36, to appear in Phys. Rev. B, Vol 50, November
1, 1994. A LaTeX file with no figure
Potts and percolation models on bowtie lattices
We give the exact critical frontier of the Potts model on bowtie lattices.
For the case of , the critical frontier yields the thresholds of bond
percolation on these lattices, which are exactly consistent with the results
given by Ziff et al [J. Phys. A 39, 15083 (2006)]. For the Potts model on
the bowtie-A lattice, the critical point is in agreement with that of the Ising
model on this lattice, which has been exactly solved. Furthermore, we do
extensive Monte Carlo simulations of Potts model on the bowtie-A lattice with
noninteger . Our numerical results, which are accurate up to 7 significant
digits, are consistent with the theoretical predictions. We also simulate the
site percolation on the bowtie-A lattice, and the threshold is
. In the simulations of bond percolation and site
percolation, we find that the shape-dependent properties of the percolation
model on the bowtie-A lattice are somewhat different from those of an isotropic
lattice, which may be caused by the anisotropy of the lattice.Comment: 18 pages, 9 figures and 3 table
The Dynamic Exponent of the Two-Dimensional Ising Model and Monte Carlo Computation of the Sub-Dominant Eigenvalue of the Stochastic Matrix
We introduce a novel variance-reducing Monte Carlo algorithm for accurate
determination of autocorrelation times. We apply this method to two-dimensional
Ising systems with sizes up to , using single-spin flip dynamics,
random site selection and transition probabilities according to the heat-bath
method. From a finite-size scaling analysis of these autocorrelation times, the
dynamical critical exponent is determined as (12)
Anomalous dynamics in two- and three- dimensional Heisenberg-Mattis spin glasses
We investigate the spectral and localization properties of unmagnetized
Heisenberg-Mattis spin glasses, in space dimensionalities and 3, at T=0.
We use numerical transfer-matrix methods combined with finite-size scaling to
calculate Lyapunov exponents, and eigenvalue-counting theorems, coupled with
Gaussian elimination algorithms, to evaluate densities of states. In we
find that all states are localized, with the localization length diverging as
, as energy . Logarithmic corrections to density of
states behave in accordance with theoretical predictions. In the
density-of-states dependence on energy is the same as for spin waves in pure
antiferromagnets, again in agreement with theoretical predictions, though the
corresponding amplitudes differ.Comment: RevTeX4, 9 pages, 9 .eps figure
Decoupling in the 1D frustrated quantum XY model and Josephson junction ladders: Ising critical behavior
A generalization of the one-dimensional frustrated quantum XY model is
considered in which the inter and intra-chain coupling constants of the two
infinite XY (planar rotor) chains have different strengths. The model can
describe the superconductor to insulator transition due to charging effects in
a ladder of Josephson junctions in a magnetic field with half a flux quantum
per plaquette. From a fluctuation-effective action, this transition is expected
to be in the universality class of the two-dimensional classical XY-Ising
model. The critical behavior is studied using a Monte Carlo transfer matrix
applied to the path-integral representation of the model and a
finite-size-scaling analysis of data on small system sizes. It is found that,
unlike the previous studied case of equal inter and intra-chain coupling
constants, the XY and Ising-like excitations of the quantum model decouple for
large interchain coupling, giving rise to pure Ising model critical behavior
for the chirality order parameter and a superconductor-insulator transition in
the universality class of the 2D classical XY model.Comment: 15 pages with figures, RevTex 3.0, INPE-93/00
- …