A quantum Monte Carlo algorithm for arbitrary high-spin Hamiltonians

We present a universal parameter-free quantum Monte Carlo algorithm for simulating arbitrary high-spin (spin greater than $1/2$) Hamiltonians. This approach extends a previously developed method by the authors for spin-$1/2$ Hamiltonians [Phys. Rev. Research 6, 013281 (2024)]. To demonstrate its applicability and versatility, we apply our method to the spin-$1$ and spin-$3/2$ quantum Heisenberg models on the square lattice. Additionally, we detail how the approach naturally extends to general Hamiltonians involving mixtures of particle species, including bosons and fermions. We have made our program code freely accessible on GitHub.