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

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

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

High-energy magnon dispersion in the half-filled Hubbard model: A comparison with La2_2CuO4_4

We use quantum Monte Carlo methods and single-mode approximation to study the magnon dispersion in the 2D half-filled Hubbard and phonon-coupled Heisenberg models. We find that in the Hubbard model with $U/t< 8$, high-energy magnon dispersion is similar to those observed in inelastic neutron scattering experiments in ${La}_2{CuO}_4$. On the other hand, our studies of a 2D Heisenberg model coupled to dynamic optical bond phonons, fails to reproduce the experimental dispersion. These results can be interpreted as evidence for intermediate $U/t$ and charge fluctuations in the cuprate materials