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.
The classical Heisenberg model is one of the most fundamental models in statistical and condensed matter physics. Extensive theoretical and numerical studies suggest that, in two dimensions, this model does not exhibit a finite-temperature phase transition but instead manifests asymptotic freedom. However, some research has also proposed the possibility of a Berezinskii-Kosterlitz-Thouless (BKT) phase transition over the years. In this study, we revisit the classical two-dimensional (2D) Heisenberg model through large-scale simulations with linear system sizes up to
L=16384. Our Monte Carlo data, without any extrapolation, clearly reveal an exponential divergence of the correlation length $\ensuremath{\xi}$ as a function of the inverse temperature $\ensuremath{\beta}$, a hallmark of asymptotic freedom. Moreover, extrapolating $\ensuremath{\xi}$ to the thermodynamic limit in the low-temperature regime achieves close agreement with the three-loop perturbative calculations. We further propose a finite-size scaling (FSS) ansatz for $\ensuremath{\xi}$, demonstrating that the pseudocritical point ${\ensuremath{\beta}}_{L}$ diverges logarithmically with
L. The thermodynamic and finite-size scaling behaviors of the magnetic susceptibility $\ensuremath{\chi}$ are also investigated and corroborate the predictions of asymptotic freedom. Our paper provides solid evidence for asymptotic freedom in the 2D Heisenberg model and advances understanding of finite-size scaling in such systems.
The properties of the spin-1/2 Heisenberg model on 1, 2, and 3-dimensional lattices are calculated using the Decoupled Cell Method of Homma et al., and these results are compared with high temperature and spin-wave expansions, and with other numerical approaches. The DCM has advantages over other Monte Carlo methods currently in wide use in that the transition probability is positive definite, there is no need to introduce an additional imaginary time, or Trotter, dimension, and the acceptance rate for transitions is comparable to that of classical lattice models. We find very good agreement between the DCM and the high temperature expansion in the temperature region where the high temperature expansion is valid, and reasonably good agreement at low temperatures with spin wave theory. The DCM fails for temperatures T < Tc which decreases with the size of the cell.
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).
Monte Carlo renormalization-group methods were applied to the classical three-component Heisenberg model on a two-dimensional lattice. Expectation values of local correlations of spins and various sized block spins were computed using traditional Monte Carlo methods. By matching quantities at different length scales generated by different Monte Carlo Hamiltonians we directly determined the renormalization of the nearest-neighbor coupling. Using these data and the results of high-temperature expansions and low-temperature renormalization-group calculations we have determined that this model does not have a phase transition. We have also obtained the amplitude for the low-temperature divergence of the susceptibility and the correlation length.