AbstractClassical Monte Carlo simulation of the Heisenberg model poorly describes many thermodynamic phenomena due to its neglect of the quantum nature of spins. Alternatively, we discuss how to semiclassically approach the quantum problem and demonstrate a simple method for introducing a locally approximate form of spin quantization. While the procedure underestimates magnetic short-range order, our results suggest a simple correction for recovering realistic spin–spin correlations above the critical temperature. Moreover, ensemble fluctuations are found to provide reasonably accurate thermodynamics, largely reproducing quantum mechanically calculated heat capacities and experimental magnetometry for ferromagnetic Fe and antiferromagnetic RbMnF₃. Extensions of the method are proposed to address remaining inaccuracies.