Quantum Monte Carlo investigation of the 2D Heisenberg model with S=1/2

The two-dimensional (2D) Heisenberg model with anisotropic exchange (Δ = 1−J _x /J _z ) and S=1/2 is investigated by the quantum Monte Carlo method. The energy, susceptibility, specific heat, spin-spin correlation functions, and correlation radius are calculated. The sublattice magnetization (σ) and the Néel temperature of the anisotropic antiferromagnet are logarithmic functions of the exchange anisotropy: 1/σ+1+0.13(1)ln(1/Δ). Crossover of the static magnetic structural factor as a function of temperature from power-law to exponential occurs for T _c /J≈0.4. The correlation radius can be approximated by 1/ξ=2.05T ^1.0(6)/exp(1.0(4)/T). For La₂CuO₄ the sublattice magnetization is calculated as σ=0.45, the exchange is J=(1125–1305) K; for Er₂CuO₄ J∼625 K and the exchange anisotropy Δ∼0.003. The temperature dependence of the static structural magnetic factor and the correlation radius above the Néel temperature in these compounds can be explained by the formation of topological excitations (spinons). Fiz. Tverd. Tela (St. Petersburg) 41, 116–121 (January 1999)