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 $S=1$, 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

Rocking isolation of a typical bridge pier on spread foundation

It has been observed that after some earthquakes a number of structures resting on spread footings responded to seismic excitation by rocking on their foundation and in some cases this enabled them to avoid failure. Through application to a standard bridge supported by direct foundations, this paper discusses the major differences in response when foundation uplift is taken into consideration. Special focus is given on the modifications of rocking response under biaxial and tri-axial excitation with respect to uniaxial excitation. It is found that inelastic rocking has a significant isolation effect. It is also shown that this effect increases under biaxial excitation while it is less sensitive to the vertical component of the earthquake. Finally, parametric analyses show that the isolation effect of foundation rocking increases as the size of the footing and the yield strength of the underlying soil decreases

The Two-Dimensional S=1 Quantum Heisenberg Antiferromagnet at Finite Temperatures

The temperature dependence of the correlation length, susceptibilities and the magnetic structure factor of the two-dimensional spin-1 square lattice quantum Heisenberg antiferromagnet are computed by the quantum Monte Carlo loop algorithm (QMC). In the experimentally relevant temperature regime the theoretically predicted asymptotic low temperature behavior is found to be not valid. The QMC results however, agree reasonably well with the experimental measurements of La2NiO4 even without considering anisotropies in the exchange interactions.Comment: 4 Pages, 1 table, 4 figure

Transition matrix Monte Carlo method for quantum systems

We propose an efficient method for Monte Carlo simulation of quantum lattice models. Unlike most other quantum Monte Carlo methods, a single run of the proposed method yields the free energy and the entropy with high precision for the whole range of temperature. The method is based on several recent findings in Monte Carlo techniques, such as the loop algorithm and the transition matrix Monte Carlo method. In particular, we derive an exact relation between the DOS and the expectation value of the transition probability for quantum systems, which turns out to be useful in reducing the statistical errors in various estimates.Comment: 6 pages, 4 figure

Energetics and geometry of excitations in random systems

Methods for studying droplets in models with quenched disorder are critically examined. Low energy excitations in two dimensional models are investigated by finding minimal energy interior excitations and by computing the effect of bulk perturbations. The numerical data support the assumptions of compact droplets and a single exponent for droplet energy scaling. Analytic calculations show how strong corrections to power laws can result when samples and droplets are averaged over. Such corrections can explain apparent discrepancies in several previous numerical results for spin glasses.Comment: 4 pages, eps files include

On infrared divergences in spin glasses

By studying the structure of infrared divergences in a toy propagator in the replica approach to the Ising spin glass below $T_c$, we suggest a possible cancellation mechanism which could decrease the degree of singularity in the loop expansion.Comment: 13 pages, Latex , revised versio

A Numerical Study of Ultrametricity in Finite Dimensional Spin Glasses

We use a constrained Monte Carlo technique to analyze ultrametric features of a 4 dimensional Edwards-Anderson spin glass with quenched couplings J=\pm 1. We find that in the large volume limit an ultrametric structure emerges quite clearly in the overlap of typical equilibrium configurations.Comment: 8 one column pages, latex, 4 figures with epsfig.st

Quantum Phase Transition of Randomly-Diluted Heisenberg Antiferromagnet on a Square Lattice

Ground-state magnetic properties of the diluted Heisenberg antiferromagnet on a square lattice are investigated by means of the quantum Monte Carlo method with the continuous-time loop algorithm. It is found that the critical concentration of magnetic sites is independent of the spin size S, and equal to the two-dimensional percolation threshold. However, the existence of quantum fluctuations makes the critical exponents deviate from those of the classical percolation transition. Furthermore, we found that the transition is not universal, i.e., the critical exponents significantly depend on S.Comment: RevTeX, 4 pages including 5 EPS figure

Ageing and dynamical scaling in the critical Ising spin glass

The non-equilibrium ageing behaviour of the 3D and 4D critical Ising spin glass is studied for both binary and gaussian disorder. The same phenomenology of the time-dependent scaling as in non-disordered magnets is found but the non-equilibrium exponents and the universal limit fluctuation-dissipation ratio depend on the distribution of the coupling constants.Comment: Latex2e, 7 pages with epl macro, 4 figures included, final for

Non-equilibrium Relations for Spin Glasses with Gauge Symmetry

We study the applications of non-equilibrium relations such as the Jarzynski equality and fluctuation theorem to spin glasses with gauge symmetry. It is shown that the exponentiated free-energy difference appearing in the Jarzynski equality reduces to a simple analytic function written explicitly in terms of the initial and final temperatures if the temperature satisfies a certain condition related to gauge symmetry. This result is used to derive a lower bound on the work done during the non-equilibrium process of temperature change. We also prove identities relating equilibrium and non-equilibrium quantities. These identities suggest a method to evaluate equilibrium quantities from non-equilibrium computations, which may be useful to avoid the problem of slow relaxation in spin glasses.Comment: 8 pages, 2 figures, submitted to JPS