Using a modified Lanczos method we have studied the two-dimensional Heisenberg antiferromagnetic model (at zero temperature) with lattices up to 24 sites. The ground-state and first-excited-state energies are evaluated. We show that the model is gapless in the thermodynamic limit in agreement with recent Monte Carlo simulations. Accurate results for the ``square ladder'' are also presented showing that this model is massive. The finite-size dependence of the results is discussed. In addition, we calculate the energy of a resonating-valence-bond state on finite lattices and its overlap with the exact ground state.
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.
We present numerical exact results for the ground state and the low-lying excitations for the spin-1/2 J1-J2 Heisenberg antiferromagnet on finite square lattices of up to N = 40 sites. Using finite-size extrapolation we determine the ground-state energy, the magnetic order parameters, the spin gap, the uniform susceptibility, as well as the spin-wave velocity and the spin stiffness as functions of the frustration parameter J2/J1. In agreement with the generally excepted scenario we find semiclassical magnetically ordered phases for J2 < J c1 2 and J2> J c2 2 separated by a gapful quantum paramagnetic phase. We estimate Jc1 2 ≈ 0.35J1 and J c2 2 ≈ 0.66J1. PACS codes: 75.10.Jm Quantized spin models 75.45.+j Macroscopic quantum phenomena in magnetic systems 1
We present results of a numerical renormalization-group study of the isotropic S=1 Heisenberg chain. The density-matrix renormalization-group techniques used allow us to calculate a variety of properties of the chain with unprecedented accuracy. The ground-state energy per site of the infinite chain is found to be
e0\ensuremath{\simeq}-1.401 484 038 971(4). Open-ended S=1 chains have effective S=1/2 spins on each end, with exponential decay of the local spin moment away from the ends, with decay length \ensuremath{\xi}\ensuremath{\simeq}6.03(1). The spin-spin correlation function also decays exponentially, and although the correlation length cannot be measured as accurately as the open-end decay length, it appears that the two lengths are identical. The string correlation function shows long-range order, with g(\ensuremath{\infty})\ensuremath{\simeq}-0.374 325 096(2). The excitation energy of the first excited state, a state with one magnon with momentum q=\ensuremath{\pi}, is the Haldane gap, which we find to be \ensuremath{\Delta}\ensuremath{\simeq}0.410 50(2). We find many low-lying excited states, including one- and two-magnon states, for several different chain lengths. The magnons have spin S=1, so the two-magnon states are singlets (S=0), triplets (S=1), and quintuplets (S=2). For magnons with momenta near \ensuremath{\pi}, the magnon-magnon interaction in the triplet channel is shown to be attractive, while in the singlet and quintuplet channels it is repulsive.
The ground-state energy, staggered magnetization, and correlation function are computed for the two-dimensional spin-(1/2 quantum antiferromagnetic Heisenberg model (believed to be relevant to the high-temperature superconducting materials without doping) on lattices up to 24\ifmmode\times\else\texttimes\fi{}24 by Monte Carlo methods extrapolated to zero temperature. A more accurate value for the ground-state energy is determined by the method of Neumann and Ulam by extrapolating from lattices up to 12\ifmmode\times\else\texttimes\fi{}12.