Dynamical Scaling from Multi-Scale Measurements

We present a new measure of the Dynamical Critical behavior: the "Multi-scale Dynamical Exponent (MDE)"Comment: 9 pages,Latex, Request figures from [email protected]

Idealized Multigrid Algorithm for Staggered Fermions

An idealized multigrid algorithm for the computation of propagators of staggered fermions is investigated. Exemplified in four-dimensional $SU(2)$ gauge fields, it is shown that the idealized algorithm preserves criticality under coarsening. The same is not true when the coarse grid operator is defined by the Galerkin prescription. Relaxation times in computations of propagators are small, and critical slowing is strongly reduced (or eliminated) in the idealized algorithm. Unfortunately, this algorithm is not practical for production runs, but the investigations presented here answer important questions of principle.Comment: 11 pages, no figures, DESY 93-046; can be formatted with plain LaTeX article styl

Measuring the Decorrelation Times of Fourier Modes in Simulations

We describe a method to study the rate at which modes decorrelate in numerical simulations. We study the XY model updated with the Metropolis and Wolff dynamics respectively and compute the rate at which each eigenvector of the dynamics decorrelates. Our method allows us to identify the decorrelation time for each mode separately. We find that the autocorrelation function of the various modes is markedly different for the `local' Metropolis compared to the `non-local' Wolff dynamics. Equipped with this new insight, it may be possible to devise highly efficient algorithms.Comment: 16 pp (LaTeX), PUPT-1378 , IASSNS-HEP-93/

Ising Model Coupled to Three-Dimensional Quantum Gravity

We have performed Monte Carlo simulations of the Ising model coupled to three-dimensional quantum gravity based on a summation over dynamical triangulations. These were done both in the microcanonical ensemble, with the number of points in the triangulation and the number of Ising spins fixed, and in the grand canoncal ensemble. We have investigated the two possible cases of the spins living on the vertices of the triangulation (``diect'' case) and the spins living in the middle of the tetrahedra (``dual'' case). We observed phase transitions which are probably second order, and found that the dual implementation more effectively couples the spins to the quantum gravity.Comment: 11 page

Block Spin Effective Action for 4d SU(2) Finite Temperature Lattice Gauge Theory

The Svetitsky-Yaffe conjecture for finite temperature 4d SU(2) lattice gauge theory is confirmed by observing matching of block spin effective actions of the gauge model with those of the 3d Ising model. The effective action for the gauge model is defined by blocking the signs of the Polyakov loops with the majority rule. To compute it numerically, we apply a variant of the IMCRG method of Gupta and Cordery.Comment: LaTeX2e, 22 pages, 8 Figure

Noncomputability Arising In Dynamical Triangulation Model Of Four-Dimensional Quantum Gravity

Computations in Dynamical Triangulation Models of Four-Dimensional Quantum Gravity involve weighted averaging over sets of all distinct triangulations of compact four-dimensional manifolds. In order to be able to perform such computations one needs an algorithm which for any given $N$ and a given compact four-dimensional manifold $M$ constructs all possible triangulations of $M$ with $\leq N$ simplices. Our first result is that such algorithm does not exist. Then we discuss recursion-theoretic limitations of any algorithm designed to perform approximate calculations of sums over all possible triangulations of a compact four-dimensional manifold.Comment: 8 Pages, LaTex, PUPT-132

Spectrum of the Dirac Operator and Multigrid Algorithm with Dynamical Staggered Fermions

Complete spectra of the staggered Dirac operator \Dirac are determined in quenched four-dimensional $SU(2)$ gauge fields, and also in the presence of dynamical fermions. Periodic as well as antiperiodic boundary conditions are used. An attempt is made to relate the performance of multigrid (MG) and conjugate gradient (CG) algorithms for propagators with the distribution of the eigenvalues of~\Dirac. The convergence of the CG algorithm is determined only by the condition number~$\kappa$ and by the lattice size. Since~$\kappa$'s do not vary significantly when quarks become dynamic, CG convergence in unquenched fields can be predicted from quenched simulations. On the other hand, MG convergence is not affected by~$\kappa$ but depends on the spectrum in a more subtle way.Comment: 19 pages, 8 figures, HUB-IEP-94/12 and KL-TH 19/94; comes as a uuencoded tar-compressed .ps-fil

Mean Field Behavior of Cluster Dynamics

The dynamic behavior of cluster algorithms is analyzed in the classical mean field limit. Rigorous analytical results below $T_c$ establish that the dynamic exponent has the value $z_{sw}=1$ for the Swendsen-Wang algorithm and $z_{uw}=0$ for the Wolff algorithm. An efficient Monte Carlo implementation is introduced, adapted for using these algorithms for fully connected graphs. Extensive simulations both above and below $T_c$ demonstrate scaling and evaluate the finite-size scaling function by means of a rather impressive collapse of the data.Comment: Revtex, 9 pages with 7 figure

Kinematics of Multigrid Monte Carlo

We study the kinematics of multigrid Monte Carlo algorithms by means of acceptance rates for nonlocal Metropolis update proposals. An approximation formula for acceptance rates is derived. We present a comparison of different coarse-to-fine interpolation schemes in free field theory, where the formula is exact. The predictions of the approximation formula for several interacting models are well confirmed by Monte Carlo simulations. The following rule is found: For a critical model with fundamental Hamiltonian H(phi), absence of critical slowing down can only be expected if the expansion of in terms of the shift psi contains no relevant (mass) term. We also introduce a multigrid update procedure for nonabelian lattice gauge theory and study the acceptance rates for gauge group SU(2) in four dimensions.Comment: 28 pages, 8 ps-figures, DESY 92-09