A Monte Carlo method applicable to a class of noncommuting spin variables is developed. The technique relies on a stochastic procedure which simulates the dominant vacuum-to-vacuum graphs of quantum many-body theory. Computational speed is only weakly dependent on the dimensionality of the system and is fast enough to handle large lattices. This method is tested on the
S=21 ferromagnetic Heisenberg model in one, two, and three dimensions and on the
S=21 XY model in one dimension where exact results are known.
We study the phase transition of the
±J Heisenberg model in three dimensions. Using a dynamical simulation method that removes a drift of the system, the existence of the spin-glass (SG) phase at low temperatures is suggested. The transition temperature is estimated to be
TSG∼0.18J from both equilibrium and off-equilibrium Monte-Carlo simulations. Our result contradicts the chirality mechanism of the phase transition reported recently by Kawamura which claims that it is not the spins but the chiralities of the spins that are ordered in Heisenberg SG systems.