The Loop Algorithm

Abstract

A review of the Loop Algorithm, its generalizations, and its relation to some other Monte Carlo techniques is given. The loop algorithm is a Quantum Monte Carlo procedure which employs nonlocal changes of worldline configurations, determined by local stochastic decisions. It is based on a formulation of quantum models of any dimension in an extended ensemble of worldlines and graphs, and is related to Swendsen-Wang algorithms. It can be represented directly on an operator level, both with a continuous imaginary time path integral and with the stochastic series expansion (SSE). It overcomes many of the difficulties of traditional worldline simulations. Autocorrelations are reduced by orders of magnitude. Grand-canonical ensembles, off-diagonal operators, and variance reduced estimators are accessible. In some cases, infinite systems can be simulated. For a restricted class of models, the fermion sign problem can be overcome. Transverse magnetic fields are handled efficiently, in contrast to strong diagonal ones. The method has been applied successfully to a variety of models for spin and charge degrees of freedom, including Heisenberg and XYZ spin models, hard-core bosons, Hubbard, and tJ-models. Due to the improved efficiency, precise calculations of asymptotic behavior and of quantum critical exponents have been possible.