The directed-loop algorithm

Abstract

The directed-loop scheme is a framework for generalized loop-type updates in quantum Monte Carlo, applicable both to world-line and stochastic series expansion methods. Here, the directed-loop equations, the solution of which gives the probabilities of the various loop-building steps, are discussed in the context of the anisotropic $S=1/2$ Heisenberg model in a uniform magnetic field. This example shows how the directed-loop concept emerges as a natural generalization of the conventional loop algorithm, where the loops are selfavoiding, to cases where selfintersection must be allowed in order to satisfy detailed balance.

10 pages, for the proceedings of "The Monte Carlo Method in the Physical Sciences: Celebrating the 50th Anniversary of the Metropolis Algorithm", Los Alamos, June 9-11, 2003