Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:cond-mat/9707221 · Advances In Physics

The loop algorithm

Hans Gerd Evertz

Published Jan 1, 2003• 295 citations
arXivHTMLPDF

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 that employs non-local 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. 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 fields. 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 t - J -models. Owing to the improved efficiency, precise calculations of asymptotic behaviour and of quantum critical exponents have been possible.

DOI
10.1080/0001873021000049195

Related papers

Quantum Monte Carlo with directed loops

Olav F. Syljuåsen, Anders W. Sandvik

We introduce the concept of directed loops in stochastic series expansion and path-integral quantum Monte Carlo methods. Using the detailed balance rules for directed loops, we show that it is possible to smoothly connect generally applicable simulation schemes (in which it is necessary to include backtracking processes in the loop construction) to more restricted loop algorithms that can be constructed only for a limited range of Hamiltonians (where backtracking can be avoided). The "algorithmic discontinuities" between general and special points (or regions) in parameter space can hence be eliminated. As a specific example, we consider the anisotropic S=1/2 Heisenberg antiferromagnet in an external magnetic field. We show that directed-loop simulations are very efficient for the full range of magnetic fields (zero to the saturation point) and anisotropies. In particular, for weak fields and anisotropies, the autocorrelations are significantly reduced relative to those of previous approaches. The back-tracking probability vanishes continuously as the isotropic Heisenberg point is approached. For the XY model, we show that back tracking can be avoided for all fields extending up to the saturation field. The method is hence particularly efficient in this case. We use directed-loop simulations to study the magnetization process in the two-dimensional Heisenberg model at very low temperatures. For LxL lattices with L up to 64, we utilize the step structure in the magnetization curve to extract gaps between different spin sectors. Finite-size scaling of the gaps gives an accurate estimate of the transverse susceptibility in the thermodynamic limit: chi( perpendicular )=0.0659+/-0.0002.

arXiv:cond-mat/0202316•Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics•Oct 2, 2002•825 citations

Loop algorithms for quantum simulations of fermion models on lattices

Naoki Kawashima, James E. Gubernatis, Hans Gerd Evertz

Two cluster algorithms, based on constructing and flipping loops, are presented for world-line quantum Monte Carlo simulations of fermions and are tested on the one-dimensional repulsive Hubbard model. We call these algorithms the loop-flip and loop-exchange algorithms. For these two algorithms and the standard world-line algorithm, we calculated the autocorrelation times for various physical quantities and found that the ordinary world-line algorithm, which uses only local moves, suffers from very long correlation times that makes not only the estimate of the error difficult but also the estimate of the average values themselves difficult. These difficulties are especially severe in the low-temperature, large-U regime. In contrast, we find that new algorithms, when used alone or in combinations with themselves and the standard algorithm, can have significantly smaller autocorrelation times, in some cases being smaller by three orders of magnitude. The new algorithms, which use nonlocal moves, are discussed from the point of view of a general prescription for developing cluster algorithms. The loop-flip algorithm is also shown to be ergodic and to belong to the grand canonical ensemble. Extensions to other models and higher dimensions are briefly discussed.

arXiv:cond-mat/9403082•Physical review. B, Condensed matter•Jul 1, 1994•56 citations

Quantum Monte Carlo loop algorithm for the t−J model

Beat Ammon, Hans Gerd Evertz, Naoki Kawashima +2 more

We propose a generalization of the quantum Monte Carlo loop algorithm to the $t\ensuremath{-}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\ensuremath{-}J$ chain, a frustrated Heisenberg chain, and $t\ensuremath{-}J$ ladder models.

arXiv:cond-mat/9711022•Physical review. B, Condensed matter•Aug 15, 1998•14 citations

Stochastic series expansion method with operator-loop update

Anders W. Sandvik

A cluster update (the ``operator loop'') is developed within the framework of a numerically exact quantum Monte Carlo method based on the power series expansion of $\mathrm{exp}(\ensuremath{-}\ensuremath{\beta}H)$ (stochastic series expansion). The method is generally applicable to a wide class of lattice Hamiltonians for which the expansion is positive definite. For some important models the operator-loop algorithm is more efficient than loop updates previously developed for ``worldline'' simulations. The method is here tested on a two-dimensional anisotropic Heisenberg antiferromagnet in a magnetic field.

arXiv:cond-mat/9902226•Physical review. B, Condensed matter•Jun 1, 1999•747 citations

Directed-Loop Quantum Monte Carlo Method for Retarded Interactions

Manuel Weber, Fakher F. Assaad, Martin Hohenadler

The directed-loop quantum Monte Carlo method is generalized to the case of retarded interactions. Using the path integral, fermion-boson or spin-boson models are mapped to actions with retarded interactions by analytically integrating out the bosons. This yields an exact algorithm that combines the highly efficient loop updates available in the stochastic series expansion representation with the advantages of avoiding a direct sampling of the bosons. The application to electron-phonon models reveals that the method overcomes the previously detrimental issues of long autocorrelation times and exponentially decreasing acceptance rates. For example, the resulting dramatic speedup allows us to investigate the Peierls quantum phase transition on chains of up to 1282 sites.

arXiv:1704.07913•Physical Review Letters•Aug 31, 2017•19 citations

Decoupling the Magnetic Field Operator as an Independent Operator Class in Stochastic Series Expansion Quantum Monte Carlo

Shijie Jin, Lu Liu

The Stochastic Series Expansion (SSE) quantum Monte Carlo method with loop updates is among the most powerful approaches for quantum spin and boson systems. In the standard formulation, the magnetic field term is routinely absorbed into the Heisenberg interactions---a strategy that has proven highly efficient across a wide range of field strengths. In this work, we propose a more general SSE framework in which the magnetic field operator is treated as an independent operator class, enabling it to participate in Monte Carlo updates on an equal footing with all other operators. Importantly, the field operator and other interaction operators can transform into one another during the update process. We demonstrate this algorithm using the two-dimensional antiferromagnetic Heisenberg model in a magnetic field as a concrete example. The simulation results agree with those of the conventional algorithm, confirming the correctness of the new formulation. A comparison of the integrated autocorrelation time shows that the new algorithm is competitive. This flexibility not only facilitates measurements of observables associated with operators but also offers potential computational advantages for a broader class of problems.

arXiv:2609.23978•Sep 21, 2026•0 citations
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 that employs non-local 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. 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 fields. 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 t - J -models. Owing to the improved efficiency, precise calculations of asymptotic behaviour and of quantum critical exponents have been possible.
We introduce the concept of directed loops in stochastic series expansion and path-integral quantum Monte Carlo methods. Using the detailed balance rules for directed loops, we show that it is possible to smoothly connect generally applicable simulation schemes (in which it is necessary to include backtracking processes in the loop construction) to more restricted loop algorithms that can be constructed only for a limited range of Hamiltonians (where backtracking can be avoided). The "algorithmic discontinuities" between general and special points (or regions) in parameter space can hence be eliminated. As a specific example, we consider the anisotropic S=1/2 Heisenberg antiferromagnet in an external magnetic field. We show that directed-loop simulations are very efficient for the full range of magnetic fields (zero to the saturation point) and anisotropies. In particular, for weak fields and anisotropies, the autocorrelations are significantly reduced relative to those of previous approaches. The back-tracking probability vanishes continuously as the isotropic Heisenberg point is approached. For the XY model, we show that back tracking can be avoided for all fields extending up to the saturation field. The method is hence particularly efficient in this case. We use directed-loop simulations to study the magnetization process in the two-dimensional Heisenberg model at very low temperatures. For LxL lattices with L up to 64, we utilize the step structure in the magnetization curve to extract gaps between different spin sectors. Finite-size scaling of the gaps gives an accurate estimate of the transverse susceptibility in the thermodynamic limit: chi( perpendicular )=0.0659+/-0.0002.
Two cluster algorithms, based on constructing and flipping loops, are presented for world-line quantum Monte Carlo simulations of fermions and are tested on the one-dimensional repulsive Hubbard model. We call these algorithms the loop-flip and loop-exchange algorithms. For these two algorithms and the standard world-line algorithm, we calculated the autocorrelation times for various physical quantities and found that the ordinary world-line algorithm, which uses only local moves, suffers from very long correlation times that makes not only the estimate of the error difficult but also the estimate of the average values themselves difficult. These difficulties are especially severe in the low-temperature, large-U regime. In contrast, we find that new algorithms, when used alone or in combinations with themselves and the standard algorithm, can have significantly smaller autocorrelation times, in some cases being smaller by three orders of magnitude. The new algorithms, which use nonlocal moves, are discussed from the point of view of a general prescription for developing cluster algorithms. The loop-flip algorithm is also shown to be ergodic and to belong to the grand canonical ensemble. Extensions to other models and higher dimensions are briefly discussed.
We propose a generalization of the quantum Monte Carlo loop algorithm to the $t\ensuremath{-}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\ensuremath{-}J$ chain, a frustrated Heisenberg chain, and $t\ensuremath{-}J$ ladder models.
A cluster update (the ``operator loop'') is developed within the framework of a numerically exact quantum Monte Carlo method based on the power series expansion of $\mathrm{exp}(\ensuremath{-}\ensuremath{\beta}H)$ (stochastic series expansion). The method is generally applicable to a wide class of lattice Hamiltonians for which the expansion is positive definite. For some important models the operator-loop algorithm is more efficient than loop updates previously developed for ``worldline'' simulations. The method is here tested on a two-dimensional anisotropic Heisenberg antiferromagnet in a magnetic field.
The directed-loop quantum Monte Carlo method is generalized to the case of retarded interactions. Using the path integral, fermion-boson or spin-boson models are mapped to actions with retarded interactions by analytically integrating out the bosons. This yields an exact algorithm that combines the highly efficient loop updates available in the stochastic series expansion representation with the advantages of avoiding a direct sampling of the bosons. The application to electron-phonon models reveals that the method overcomes the previously detrimental issues of long autocorrelation times and exponentially decreasing acceptance rates. For example, the resulting dramatic speedup allows us to investigate the Peierls quantum phase transition on chains of up to 1282 sites.
The Stochastic Series Expansion (SSE) quantum Monte Carlo method with loop updates is among the most powerful approaches for quantum spin and boson systems. In the standard formulation, the magnetic field term is routinely absorbed into the Heisenberg interactions---a strategy that has proven highly efficient across a wide range of field strengths. In this work, we propose a more general SSE framework in which the magnetic field operator is treated as an independent operator class, enabling it to participate in Monte Carlo updates on an equal footing with all other operators. Importantly, the field operator and other interaction operators can transform into one another during the update process. We demonstrate this algorithm using the two-dimensional antiferromagnetic Heisenberg model in a magnetic field as a concrete example. The simulation results agree with those of the conventional algorithm, confirming the correctness of the new formulation. A comparison of the integrated autocorrelation time shows that the new algorithm is competitive. This flexibility not only facilitates measurements of observables associated with operators but also offers potential computational advantages for a broader class of problems.