Chaos in a Two-Dimensional Ising Spin Glass

We study chaos in a two dimensional Ising spin glass by finite temperature Monte Carlo simulations. We are able to detect chaos with respect to temperature changes as well as chaos with respect to changing the bonds, and find that the chaos exponents for these two cases are equal. Our value for the exponent appears to be consistent with that obtained in studies at zero temperature.Comment: 4 pages, LaTeX, 4 postscript figures included. The analysis of the data is now done somewhat differently. The results are consistent with the chaos exponent found at zero temperature. Additional papers of PY can be obtained on-line at http://schubert.ucsc.edu/pete

Quantum simulations of the superfluid-insulator transition for two-dimensional, disordered, hard-core bosons

We introduce two novel quantum Monte Carlo methods and employ them to study the superfluid-insulator transition in a two-dimensional system of hard-core bosons. One of the methods is appropriate for zero temperature and is based upon Green's function Monte Carlo; the other is a finite-temperature world-line cluster algorithm. In each case we find that the dynamical exponent is consistent with the theoretical prediction of $z=2$ by Fisher and co-workers.Comment: Revtex, 10 pages, 3 figures (postscript files attached at end, separated by %%%%%% Fig # %%%%%, where # is 1-3). LA-UR-94-270

Quantum Monte Carlo Loop Algorithm for the t-J Model

We propose a generalization of the Quantum Monte Carlo loop algorithm to the t-J model by a mapping to three coupled six-vertex models. The autocorrelation times are reduced by orders of magnitude compared to the conventional local algorithms. The method is completely ergodic and can be formulated directly in continuous time. We introduce improved estimators for simulations with a local sign problem. Some first results of finite temperature simulations are presented for a t-J chain, a frustrated Heisenberg chain, and t-J ladder models.Comment: 22 pages, including 12 figures. RevTex v3.0, uses psf.te

Kosterlitz-Thouless transition of quantum XY model in two dimensions

The two-dimensional $S=1/2$ XY model is investigated with an extensive quantum Monte Carlo simulation. The helicity modulus is precisely estimated through a continuous-time loop algorithm for systems up to $128 \times 128$ near and below the critical temperature. The critical temperature is estimated as $T_{\rm KT} = 0.3427(2)J$. The obtained estimates for the helicity modulus are well fitted by a scaling form derived from the Kosterlitz renormalization group equation. The validity of the Kosterlitz-Thouless theory for this model is confirmed.Comment: 8 pages, 2 tables, 6 figure

Universality and universal finite-size scaling functions in four-dimensional Ising spin glasses

We study the four-dimensional Ising spin glass with Gaussian and bond-diluted bimodal distributed interactions via large-scale Monte Carlo simulations and show via an extensive finite-size scaling analysis that four-dimensional Ising spin glasses obey universality.Comment: 12 pages, 9 figures, 4 table

Accessing the dynamics of large many-particle systems using Stochastic Series Expansion

The Stochastic Series Expansion method (SSE) is a Quantum Monte Carlo (QMC) technique working directly in the imaginary time continuum and thus avoiding "Trotter discretization" errors. Using a non-local "operator-loop update" it allows treating large quantum mechanical systems of many thousand sites. In this paper we first give a comprehensive review on SSE and present benchmark calculations of SSE's scaling behavior with system size and inverse temperature, and compare it to the loop algorithm, whose scaling is known to be one of the best of all QMC methods. Finally we introduce a new and efficient algorithm to measure Green's functions and thus dynamical properties within SSE.Comment: 11 RevTeX pages including 7 figures and 5 table

Quantum Monte Carlo Study on Magnetization Processes

A quantum Monte Carlo method combining update of the loop algorithm with the global flip of the world line is proposed as an efficient method to study the magnetization process in an external field, which has been difficult because of inefficiency of the update of the total magnetization. The method is demonstrated in the one dimensional antiferromagnetic Heisenberg model and the trimer model. We attempted various other Monte Carlo algorithms to study systems in the external field and compared their efficiency.Comment: 5 pages, 9 figures; added references for section 1, corrected typo

Dual Monte Carlo and Cluster Algorithms

We discuss the development of cluster algorithms from the viewpoint of probability theory and not from the usual viewpoint of a particular model. By using the perspective of probability theory, we detail the nature of a cluster algorithm, make explicit the assumptions embodied in all clusters of which we are aware, and define the construction of free cluster algorithms. We also illustrate these procedures by rederiving the Swendsen-Wang algorithm, presenting the details of the loop algorithm for a worldline simulation of a quantum $S=$ 1/2 model, and proposing a free cluster version of the Swendsen-Wang replica method for the random Ising model. How the principle of maximum entropy might be used to aid the construction of cluster algorithms is also discussed.Comment: 25 pages, 4 figures, to appear in Phys.Rev.