Spin-wave expansions are carried out for the Hamiltonian version of the O(2) Heisenberg spin model, and evaluated for the cases of the (1+1)-dimensional linear chain and the (2+1)-dimensional square and triangular lattices. Expressions are obtained for the bulk limits and finite-size corrections to the ground-state energy, the one-boson excitation energy, the magnetization and susceptibility at low 'temperatures'. Results are derived for the critical parameters in the (1+1)-dimensional model, and the Goldstone mode parameters in the (2+1)-dimensional models. Comparison is made with the finite-size scaling predictions of Fisher and Privman (1985).
Spin-wave expansions are carried out for the Hamiltonian version of the O(2) Heisenberg spin model, and evaluated for the cases of the (1+1)-dimensional linear chain and the (2+1)-dimensional square and triangular lattices. Expressions are obtained for the bulk limits and finite-size corrections to the ground-state energy, the one-boson excitation energy, the magnetization and susceptibility at low 'temperatures'. Results are derived for the critical parameters in the (1+1)-dimensional model, and the Goldstone mode parameters in the (2+1)-dimensional models. Comparison is made with the finite-size scaling predictions of Fisher and Privman (1985).
A linear chain of S =1/2 spins with the isotropic Hamiltonian H = J Σ( S i x S i +1 x + S i y S i +1 y + S i z S i +1 z -1/4), ( J <0) , at low temperatures is studied numerically by solving a set of nonlinear integral equations which are derived from the Bethe Ansatz. We find that the free energy and susceptibility are both expanded in T/∣J∣ about T =0. The critical exponents are a =-1/2 and y =2. The coefficient of T 3/2 term in the free energy agrees with the spin wave theory with the accuracy of 0.1%. The coefficient of T -2 term in the zero-field susceptibility agrees with Fisher's solution of the classical Heisenberg model with the accuracy of 0.3%. The contribution from the second and higher order terms of the expansions in these thermodynamic quantities is not negligible in the region T /| J |>0.002.
The ground-state parameters of the two-dimensional S=1/2 antiferromagnetic Heisenberg model are calculated using the stochastic series expansion quantum Monte Carlo method for $L\ifmmode\times\else\texttimes\fi{}L$ lattices with L up to 16. The finite-size results for the energy E, the sublattice magnetization M, the long-wavelength susceptibility ${\ensuremath{\chi}}_{\ensuremath{\perp}}(q=2\ensuremath{\pi}/L)$, and the spin stiffness ${\ensuremath{\rho}}_{s},$ are extrapolated to the thermodynamic limit using fits to polynomials in 1/L, constrained by scaling forms previously obtained from renormalization-group calculations for the nonlinear $\ensuremath{\sigma}$ model and chiral perturbation theory. The results are fully consistent with the predicted leading finite-size corrections, and are of sufficient accuracy for extracting also subleading terms. The subleading energy correction $(\ensuremath{\sim}{1/L}^{4})$ agrees with chiral perturbation theory to within a statistical error of a few percent, thus providing numerical confirmation of the finite-size scaling forms to this order. The extrapolated ground- state energy per spin is $E=\ensuremath{-}0.669437(5)$. The result from previous Green's function Monte Carlo (GFMC) calculations is slightly higher than this value, most likely due to a small systematic error originating from ``population control'' bias in GFMC. The other extrapolated parameters are M=0.3070(3), ${\ensuremath{\rho}}_{s}=0.175(2)$, ${\ensuremath{\chi}}_{\ensuremath{\perp}}=0.0625(9)$, and the spin-wave velocity c=1.673(7). The statistical errors are comparable with those of previous estimates obtained by fitting loop algorithm quantum Monte Carlo data to finite-temperature scaling forms. Both M and ${\ensuremath{\rho}}_{s}$ obtained from the finite-T data are, however, a few error bars higher than the present estimates. It is argued that the T=0 extrapolations performed here are less sensitive to effects of neglected higher-order corrections, and therefore should be more reliable.
Abstract Motivated by the need for high-precision benchmark results for basic quantum many-body models, several ground state properties of the square lattice S = 1 / 2 Heisenberg antiferromagnet are computed (the energy, order parameter, spin stiffness, spinwave velocity, long-wavelength susceptibility and staggered susceptibility) using extensive quantum Monte Carlo simulations with the stochastic series expansion method. Moderately sized lattices are studied at temperatures T sufficiently low to realize the T → 0 limit, with significantly reduced statistical errors relative to previous works. Results for periodic L × L lattices with L ∈ [ 6 , 96 ] are tabulated versus L and extrapolations to infinite system size are carried out. The relative statistical errors of finite-size energies are below 10 −7 . The extrapolated ground state energy density is e 0 = − 0.669 441 857 ( 7 ) , where the number within parentheses is the statistical error (one standard deviation) of the preceding digit. This value represents an improvement in precision of three orders of magnitude over the previous best result. The leading and subleading finite-size corrections to e 0 are in full quantitative agreement with predictions from chiral perturbation theory, thus further supporting the soundness of both the extrapolations and the theory. The extrapolated sublattice magnetization is m s = 0.307 447 ( 2 ) , which agrees well with previous estimates but with a much smaller statistical error. The coefficient of the linear in L −1 correction to m s 2 agrees with the value from chiral perturbation theory and the presence of a factor ln γ ( L ) in the second-order correction is also confirmed, with the previously unknown value of the exponent being γ = 0.82 ( 4 ) . The finite-size corrections to the staggered susceptibility point to logarithmic corrections in this quantity as well. To facilitate benchmarking of methods for which periodic boundary conditions are challenging, results for systems with open and cylindrical boundaries are also listed and their spatially inhomogeneous order parameters are analyzed.
The Green's-function Monte Carlo (GFMC) method is used to calculate very accurate ground-state energies of the two-dimensional, spin-1/2 Heisenberg antiferromagnet. The computations are performed on L\ifmmode\times\else\texttimes\fi{}L square lattices up to L=16 with varying uniform magnetization, which allows the extraction of the perpendicular susceptibility (\ensuremath{\chi}) and spin-wave velocity (c). These two quantities are the lattice- or cutoff-dependent parameters that allow one to map the long-wavelength properties of the antiferromagnet onto the nonlinear \ensuremath{\sigma} model and so are of general interest. Systematic errors present in previous GFMC calculations are addressed and corrected to yield results in excellent agreement with other numerical methods. I find, for the ground-state energy per site, -0.669 34(3); the susceptibility renormalization factor, ${\mathit{Z}}_{\mathrm{\ensuremath{\chi}}}$=0.535(5); and the spin-wave velocity renormalization factor, Zc=1.10(3). Finite-size effects in the extraction of Zc and ${\mathit{Z}}_{\mathrm{\ensuremath{\chi}}}$ are discussed. The value of ${\mathit{Z}}_{\mathrm{\ensuremath{\chi}}}$ computed here is in agreement with the series-expansion results of Singh and of Zheng, Oitmaa, and Hamer, thereby clearing up a previous inconsistency between the series-expansion and quantum Monte Carlo predictions.