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.