Comment on ``Quantum Phase Transition of the Randomly Diluted Heisenberg Antiferromagnet on a Square Lattice''

In Phys. Rev. Lett. 84, 4204 (2000) (cond-mat/9905379), Kato et al. presented quantum Monte Carlo results indicating that the critical concentration of random non-magnetic sites in the two-dimensional antiferromagnetic Heisenberg model equals the classical percolation density; pc=0.407254. The data also suggested a surprising dependence of the critical exponents on the spin S of the magnetic sites, with a gradual approach to the classical percolation exponents as S goes to infinity. I here argue that the exponents in fact are S-independent and equal to those of classical percolation. The apparent S-dependent behavior found by Kato et al. is due to temperature effects in the simulations as well as a quantum effect that masks the true asymptotic scaling behavior for small lattices.Comment: Comment on Phys. Rev. Lett. 84, 4204 (2000), by K. Kato et al.; 1 page, 1 figur

Spin nematic ground state of the triangular lattice S=1 biquadratic model

Motivated by the spate of recent experimental and theoretical interest in Mott insulating S=1 triangular lattice magnets, we consider a model S=1 Hamiltonian on a triangular lattice interacting with rotationally symmetric biquadratic interactions. We show that the partition function of this model can be expressed in terms of configurations of three colors of tightly-packed, closed loops with {\em non-negative} weights, which allows for efficient quantum Monte Carlo sampling on large lattices. We find the ground state has spin nematic order, i.e. it spontaneously breaks spin rotation symmetry but preserves time reversal symmetry. We present accurate results for the parameters of the low energy field theory, as well as finite-temperature thermodynamic functions

Spin-Peierls transition in the Heisenberg chain with finite-frequency phonons

We study the spin-Peierls transition in a Heisenberg spin chain coupled to optical bond phonons. Quantum Monte Carlo results for systems with up to N=256 spins show unambiguously that the transition occurs only when the spin-phonon coupling α exceeds a critical value α_c. Using sum rules, we show that the phonon spectral function has divergent (for N→∞) weight extending to zero frequency for α<α_c. The phonon correlations decay with distance r as 1/r. This behavior is characteristic for all 0<α<α_c and the q=π phonon does not soften (to zero frequency) at α=α_c.First author draf

Comment on ``Ground State Phase Diagram of a Half-Filled One-Dimensional Extended Hubbard Model''

In Phys. Rev. Lett. 89, 236401 (2002), Jeckelmann argued that the recently discovered bond-order-wave (BOW) phase of the 1D extended Hubbard model does not have a finite extent in the (U,V) plane, but exists only on a segment of a first-order SDW-CDW phase boundary. We here present quantum Monte Carlo result of higher precision and for larger system sizes than previously and reconfirm that the BOW phase does exist a finite distance away from the phase boundary, which hence is a BOW-CDW transition curve.Comment: 1 page, 1 figure, v2: final published versio

Monte Carlo Simulations of Quantum Spin Systems in the Valence Bond Basis

We discuss a projector Monte Carlo method for quantum spin models formulated in the valence bond basis, using the S=1/2 Heisenberg antiferromagnet as an example. Its singlet ground state can be projected out of an arbitrary basis state as the trial state, but a more rapid convergence can be obtained using a good variational state. As an alternative to first carrying out a time consuming variational Monte Carlo calculation, we show that a very good trial state can be generated in an iterative fashion in the course of the simulation itself. We also show how the properties of the valence bond basis enable calculations of quantities that are difficult to obtain with the standard basis of Sz eigenstates. In particular, we discuss quantities involving finite-momentum states in the triplet sector, such as the dispersion relation and the spectral weight of the lowest triplet.Comment: 15 pages, 7 figures, for the proceedings of "Computer Simulation Studies in Condensed Matter Physics XX

Ground state of the random-bond spin-1 Heisenberg chain

Stochastic series expansion quantum Monte Carlo is used to study the ground state of the antiferromagnetic spin-1 Heisenberg chain with bond disorder. Typical spin- and string-correlations functions behave in accordance with real-space renormalization group predictions for the random-singlet phase. The average string-correlation function decays algebraically with an exponent of -0.378(6), in very good agreement with the prediction of $-(3-\sqrt{5})/2\simeq -0.382$, while the average spin-correlation function is found to decay with an exponent of about -1, quite different from the expected value of -2. By implementing the concept of directed loops for the spin-1 chain we show that autocorrelation times can be reduced by up to two orders of magnitude.Comment: 9 pages, 10 figure