A New Approach to Stochastic State selections in Quantum Spin Systems

We propose a new type of Monte Carlo approach in numerical studies of quantum systems. Introducing a probability function which determines whether a state in the vector space survives or not, we can evaluate expectation values of powers of the Hamiltonian from a small portion of the full vector space. This method is free from the negative sign problem because it is not based on importance sampling techniques. In this paper we describe our method and, in order to examine how effective it is, present numerical results on the 4x4, 6x6 and 8x8 Heisenberg spin one-half model. The results indicate that we can perform useful evaluations with limited computer resources. An attempt to estimate the lowest energy eigenvalue is also stated.Comment: 10 pages, 2 figures, 8 table

Time of flight observables and the formation of Mott domains of fermions and bosons on optical lattices

We study, using quantum Monte Carlo simulations, the energetics of the formation of Mott domains of fermions and bosons trapped on one-dimensional lattices. We show that, in both cases, the sum of kinetic and interaction energies exhibits minima when Mott domains appear in the trap. In addition, we examine the derivatives of the kinetic and interaction energies, and of their sum, which display clear signatures of the Mott transition. We discuss the relevance of these findings to time-of-flight experiments that could allow the detection of the metal--Mott-insulator transition in confined fermions on optical lattices, and support established results on the superfluid--Mott-insulator transition in confined bosons on optical lattices.Comment: 5 pages, 6 figures, published versio

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

Universal scaling at field-induced magnetic phase transitions

We study field-induced magnetic order in cubic lattices of dimers with antiferromagnetic Heisenberg interactions. The thermal critical exponents at the quantum phase transition from a spin liquid to a magnetically ordered phase are determined from Stochastic Series Expansion Quantum Monte Carlo simulations. These exponents are independent of the interdimer coupling ratios, and converge to the value obtained by considering the transition as a Bose-Einstein condensation of magnons, alpha_(BEC) = 1.5. The scaling results are of direct relevance to the spin-dimer systems TlCuCl_3 and KCuCl_3, and explain the broad range of exponents reported for field-induced ordering transitions.Comment: 4 pages, 4 eps-figure

Spin dynamics of SrCu2_2O3_3 and the Heisenberg ladder

The $S=1/2$ Heisenberg antiferromagnet in the ladder geometry is studied as a model for the spin degrees of freedom of SrCu$_2$O$_3$. The susceptibility and the spin echo decay rate are calculated using a quantum Monte Carlo technique, and the spin-lattice relaxation rate is obtained by maximum entropy analytic continuation of imaginary time correlation functions. All calculated quantities are in reasonable agreement with experimental results for SrCu$_2$O$_3$ if the exchange coupling $J \approx 850$K, i.e. significantly smaller than in high-T$_c$ cuprates.Comment: 11 pages (Revtex) + 3 uuencoded ps files. To appear in Phys. Rev. B, Rapid Com

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

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

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

Specific heat of quasi-2D antiferromagnetic Heisenberg models with varying inter-planar couplings

We have used the stochastic series expansion (SSE) quantum Monte Carlo (QMC) method to study the three-dimensional (3D) antiferromagnetic Heisenberg model on cubic lattices with in-plane coupling J and varying inter-plane coupling J_perp < J. The specific heat curves exhibit a 3D ordering peak as well as a broad maximum arising from short-range 2D order. For J_perp << J, there is a clear separation of the two peaks. In the simulations, the contributions to the total specific heat from the ordering across and within the layers can be separated, and this enables us to study in detail the 3D peak around T_c (which otherwise typically is dominated by statistical noise). We find that the peak height decreases with decreasing J_perp, becoming nearly linear below J_perp = 0.2J. The relevance of these results to the lack of observed specific heat anomaly at the ordering transition of some quasi-2D antiferromagnets is discussed.Comment: 7 pages, 8 figure