Loop updates for variational and projector quantum Monte Carlo simulations in the valence-bond basis

We show how efficient loop updates, originally developed for Monte Carlo simulations of quantum spin systems at finite temperature, can be combined with a ground-state projector scheme and variational calculations in the valence bond basis. The methods are formulated in a combined space of spin z-components and valence bonds. Compared to schemes formulated purely in the valence bond basis, the computational effort is reduced from up to O(N^2) to O(N) for variational calculations, where N is the system size, and from O(m^2) to O(m) for projector simulations, where m>> N is the projection power. These improvements enable access to ground states of significantly larger lattices than previously. We demonstrate the efficiency of the approach by calculating the sublattice magnetization M_s of the two-dimensional Heisenberg model to high precision, using systems with up to 256*256 spins. Extrapolating the results to the thermodynamic limit gives M_s=0.30743(1). We also discuss optimized variational amplitude-product states, which were used as trial states in the projector simulations, and compare results of projecting different types of trial states.Comment: 12 pages, 9 figures. v2: Significantly expanded, to appear in Phys. Rev.

Quantum Monte Carlo simulations of bosonic and fermionic impurities in a two-dimensional hard-core boson system

A two-dimensional lattice hard-core boson system with a small fraction of bosonic or fermionic impurity particles is studied. The impurities have the same hopping and interactions as the dominant bosons and their effects are solely due to quantum statistics. Quantum Monte Carlo simulations are carried out in which paths of the dominant boson species are sampled and a summation is performed over all second-species paths compatible with the permutation cycles. Both kinds of impurities reduce modestly and equally the Kosterliz-Thouless superfluid transition temperature. However, the effective impurity interactions are found to be qualitatively different at lower temperatures; fermions are repulsive and further suppress superfluidity at low temperatures.Comment: 4 pages, 5 figure

Accessing the dynamics of large many-particle systems using Stochastic Series Expansion

The Stochastic Series Expansion method (SSE) is a Quantum Monte Carlo (QMC) technique working directly in the imaginary time continuum and thus avoiding "Trotter discretization" errors. Using a non-local "operator-loop update" it allows treating large quantum mechanical systems of many thousand sites. In this paper we first give a comprehensive review on SSE and present benchmark calculations of SSE's scaling behavior with system size and inverse temperature, and compare it to the loop algorithm, whose scaling is known to be one of the best of all QMC methods. Finally we introduce a new and efficient algorithm to measure Green's functions and thus dynamical properties within SSE.Comment: 11 RevTeX pages including 7 figures and 5 table

Thermodynamics of a gas of deconfined bosonic spinons in two dimensions

We consider the quantum phase transition between a Neel antiferromagnet and a valence-bond solid (VBS) in a two-dimensional system of S=1/2 spins. Assuming that the excitations of the critical ground state are linearly dispersing deconfined spinons obeying Bose statistics, we derive expressions for the specific heat and the magnetic susceptibility at low temperature T. Comparing with quantum Monte Carlo results for the J-Q model, which is a candidate for a deconfined Neel-VBS transition, we find excellent agreement, including a previously noted logarithmic correction in the susceptibility. In our treatment, this is a direct consequence of a confinement length scale Lambda which is proportional to the correlation length xi raised to a non-trivial power; Lambda ~ xi^(1+a) ~1/T^(1+a), with a>0 (with a approximately 0.2 in the model).Comment: 4+ pages, 3 figures. v2: cosmetic changes onl

Susceptibility of the 2D S=1/2 Heisenberg antiferromagnet with an impurity

We use a quantum Monte Carlo method (stochastic series expansion) to study the effects of a magnetic or nonmagnetic impurity on the magnetic susceptibility of the two-dimensional Heisenberg antiferromagnet. At low temperatures, we find a log-divergent contribution to the transverse susceptibility. We also introduce an effective few-spin model that can quantitatively capture the differences between magnetic and nonmagnetic impurities at high and intermediate temperatures.Comment: 5 pages, 4 figures, v2: Updated data in figures, minor changes in text, v3: Final version, cosmetic change

The sign problem in Monte Carlo simulations of frustrated quantum spin systems

We discuss the sign problem arising in Monte Carlo simulations of frustrated quantum spin systems. We show that for a class of ``semi-frustrated'' systems (Heisenberg models with ferromagnetic couplings $J_z(r) < 0$ along the $z$-axis and antiferromagnetic couplings $J_{xy}(r)=-J_z(r)$ in the $xy$-plane, for arbitrary distances $r$) the sign problem present for algorithms operating in the $z$-basis can be solved within a recent ``operator-loop'' formulation of the stochastic series expansion method (a cluster algorithm for sampling the diagonal matrix elements of the power series expansion of ${\rm exp}(-\beta H)$ to all orders). The solution relies on identification of operator-loops which change the configuration sign when updated (``merons'') and is similar to the meron-cluster algorithm recently proposed by Chandrasekharan and Wiese for solving the sign problem for a class of fermion models (Phys. Rev. Lett. {\bf 83}, 3116 (1999)). Some important expectation values, e.g., the internal energy, can be evaluated in the subspace with no merons, where the weight function is positive definite. Calculations of other expectation values require sampling of configurations with only a small number of merons (typically zero or two), with an accompanying sign problem which is not serious. We also discuss problems which arise in applying the meron concept to more general quantum spin models with frustrated interactions.Comment: 13 pages, 16 figure

Stochastic series expansion method with operator-loop update

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 exp(-BH) (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.Comment: 5 pages, 4 figure

Stochastic series expansion algorithm for the S=1/2 XY model with four-site ring exchange

We describe a stochastic series expansion (SSE) quantum Monte Carlo method for a two-dimensional S=1/2 XY-model (or, equivalently, hard-core bosons at half-filling) which in addition to the standard pair interaction J includes a four-particle term K that flips spins on a square plaquette. The model has three ordered ground state phases; for K/J<8 it has long-range xy spin order (superfluid bosons), for K/J>15 it has staggered spin order in the z direction (charge-density-wave), and between these phases it is in a state with columnar order in the bond and plaquette energy densities. We discuss an implementation of directed-loop updates for the SSE simulations of this model and also introduce a "multi-branch" cluster update which significantly reduces the autocorrelation times for large K/J. In addition to the pure J-K model, which in the z basis has only off-diagonal terms, we also discuss modifications of the algorithm needed when various diagonal interactions are included.Comment: 23 pages, 21 figure

Double-layer Heisenberg antiferromagnet at finite temperature: Brueckner Theory and Quantum Monte Carlo simulations

The double-layer Heisenberg antiferromagnet with intra- and inter-layer couplings $J$ and $J_\perp$ exhibits a zero temperature quantum phase transition between a quantum disordered dimer phase for $g>g_c$ and a Neel phase with long range antiferromagnetic order for $g<g_c$, where $g=J_\perp/J$ and $g_c \approx 2.5$. We consider the behavior of the system at finite temperature for $g \ge g_c$ using two different and complementary approaches; an analytical Brueckner approximation and numerically exact quantum Monte Carlo simulations. We calculate the temperature dependent spin excitation spectrum (including the triplet gap), dynamic and static structure factors, the specific heat, and the uniform magnetic susceptibility. The agreement between the analytical and numerical approaches is excellent. For $T \to 0$ and $g \to g_c$, our analytical results for the specific heat and the magnetic susceptibility coincide with those previously obtained within the nonlinear $\sigma$ model approach for $N\to \infty$. Our quantum Monte Carlo simulations extend to significantly lower temperatures than previously, allowing us to obtain accurate results for the asymptotic quantum critical behavior. We also obtain an improved estimate for the critical coupling: $g_c = 2.525 \pm 0.002$.Comment: 23 pages, 12 figure

Quantum Monte Carlo with Directed Loops

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 back-tracking processes in the loop construction) to more restricted loop algorithms that can be constructed only for a limited range of Hamiltonians (where back-tracking 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 2D 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_perp = 0.0659 +- 0.0002.Comment: v2: Revised and expanded discussion of detailed balance, error in algorithmic phase diagram corrected, to appear in Phys. Rev.