Optimal phase space sampling for Monte Carlo simulations of Heisenberg spin systems