Abstract
An exact Monte Carlo method for calculating thermodynamic properties of quantum spin systems is described. Results for one-dimensional ferromagnetic and antiferromagnetic systems and for three-dimensional ferromagnetic systems show that the method can be used to study quantum spin systems as extensively as classical spin systems are studied with conventional Monte Carlo methods.