NMR relaxation rates for the spin-1/2 Heisenberg chain

The spin-lattice relaxation rate $1/T_1$ and the spin echo decay rate $1/T_{2G}$ for the spin-$1\over 2$ antiferromagnetic Heisenberg chain are calculated using quantum Monte Carlo and maximum entropy analytic continuation. The results are compared with recent analytical calculations by Sachdev. If the nuclear hyperfine form factor $A_q$ is strongly peaked around $q=\pi$ the predicted low-temperature behavior [$1/T_1 \sim \ln{^{1/2}(1/T)}$, $1/T_{2G} \sim \ln{^{1/2}(1/T)}/\sqrt{T}$] extends up to temperatures as high as $T/J \approx 0.5$. If $A_q$ has significant weight for $q \approx 0$ there are large contributions from diffusive long-wavelength processes not taken into account in the theory, and very low temperatures are needed in order to observe the asymptotic $T \to 0$ forms.Comment: 9 pages, Revtex 3.0, 5 uuencoded ps figures To appear in Phys. Rev. B, Rapid Com

Critical temperature and the transition from quantum to classical order parameter fluctuations in the three-dimensional Heisenberg antiferromagnet

We present results of extensive quantum Monte Carlo simulations of the three-dimensional (3D) S=1/2 Heisenberg antiferromagnet. Finite-size scaling of the spin stiffness and the sublattice magnetization gives the critical temperature Tc/J = 0.946 +/- 0.001. The critical behavior is consistent with the classical 3D Heisenberg universality class, as expected. We discuss the general nature of the transition from quantum mechanical to classical (thermal) order parameter fluctuations at a continuous Tc > 0 phase transition.Comment: 5 pages, Revtex, 4 PostScript figures include

Numerical Linked-Cluster Algorithms. I. Spin systems on square, triangular, and kagome lattices

We discuss recently introduced numerical linked-cluster (NLC) algorithms that allow one to obtain temperature-dependent properties of quantum lattice models, in the thermodynamic limit, from exact diagonalization of finite clusters. We present studies of thermodynamic observables for spin models on square, triangular, and kagome lattices. Results for several choices of clusters and extrapolations methods, that accelerate the convergence of NLC, are presented. We also include a comparison of NLC results with those obtained from exact analytical expressions (where available), high-temperature expansions (HTE), exact diagonalization (ED) of finite periodic systems, and quantum Monte Carlo simulations.For many models and properties NLC results are substantially more accurate than HTE and ED.Comment: 14 pages, 16 figures, as publishe

Variational quantum Monte Carlo simulations with tensor-network states

We show that the formalism of tensor-network states, such as the matrix product states (MPS), can be used as a basis for variational quantum Monte Carlo simulations. Using a stochastic optimization method, we demonstrate the potential of this approach by explicit MPS calculations for the transverse Ising chain with up to N=256 spins at criticality, using periodic boundary conditions and D*D matrices with D up to 48. The computational cost of our scheme formally scales as ND^3, whereas standard MPS approaches and the related density matrix renromalization group method scale as ND^5 and ND^6, respectively, for periodic systems.Comment: 4+ pages, 2 figures. v2: improved data, comparisons with exact results, to appear in Phys Rev Let

Effects of intrabilayer coupling on the magnetic properties of YBa2_2Cu3_3O6_6

A two-layer Heisenberg antiferromagnet is studied as a model of the bilayer cuprate YBa$_2$Cu$_3$O$_6$. Quantum Monte Carlo results are presented for the temperature dependence of the spin correlation length, the static structure factor, the magnetic susceptibility, and the $^{63}$Cu NMR spin-echo decay rate $1/T_{2G}$. As expected, when the ratio $J_2/J_1$ of the intrabilayer and in-plane coupling strengths is small, increasing $J_2$ pushes the system deeper inside the renormalized classical regime. Even for $J_2/J_1$ as small as $0.1$ the correlations are considerably enhanced at temperatures as high as $T/J_1 \approx 0.4-0.5$. This has a significant effect on $1/T_{2G}$, and it is suggested that measurements of this quantity at high temperatures can reveal the strength of the intrabilayer coupling in YBa$_2$Cu$_3$O$_6$.Comment: 10 pages (Revtex) + 5 uuencoded ps figures. To appear in Phys. Rev. B, Rapid Com

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 Cluster Series expansion for quantum spin systems

In this paper we develop a cluster-variant of the Stochastic Series expansion method (SCSE). For certain systems with longer-range interactions the SCSE is considerably more efficient than the standard implementation of the Stochastic Series Expansion (SSE), at low temperatures. As an application of this method we calculated the T=0-conductance for a linear chain with a (diagonal) next nearest neighbor interaction.Comment: 5 pages, 7 figure

Quantum-critical scaling and temperature-dependent logarithmic corrections in the spin-half Heisenberg chain

Low temperature dynamics of the S=1/2 Heisenberg chain is studied via a simple ansatz generalizing the conformal mapping and analytic continuation procedures to correlation functions with multiplicative logarithmic factors. Closed form expressions for the dynamic susceptibility and the NMR relaxation rates 1/T_1 and 1/T_{2G} are obtained, and are argued to improve the agreement with recent experiments. Scaling in q/T and \omega/T are violated due to these logarithmic terms. Numerical results show that the logarithmic corrections are very robust. While not yet in the asymptotic low temperature regime, they provide striking qualitative confirmation of the theoretical results.Comment: 4 pages, RevTex, 4 postscript figures in one fil

Two-Dimensional Quantum XY Model with Ring Exchange and External Field

We present the zero-temperature phase diagram of a square lattice quantum spin 1/2 XY model with four-site ring exchange in a uniform external magnetic field. Using quantum Monte Carlo techniques, we identify various quantum phase transitions between the XY-order, striped or valence bond solid, staggered Neel antiferromagnet and fully polarized ground states of the model. We find no evidence for a quantum spin liquid phase.Comment: 4 pages, 4 figure