Monte Carlo Determination of the Low-Energy Constants of a Spin 1/2 Heisenberg Model with Spatial Anisotropy
Motivated by the possible mechanism for the pinning of the electronic liquid crystal direction in YBCO as proposed in \cite{Pardini08}, we use the first principles Monte Carlo method to study the spin 1/2 Heisenberg model with antiferromagnetic couplings and on the square lattice. The corresponding low-energy constants, namely the spin stiffness , the staggered magnetization density , the spin wave velocity , as well as the ground state energy density are determined by fitting the Monte Carlo data to the predictions of magnon chiral perturbation theory. In particular, the spin stiffnesses and are investigated as a function of the ratio of the couplings. Although we find a good agreement between our results with those obtained by the series expansion method in the weakly anisotropic regime, for strong anisotropy we observe discrepancies.