Universal SSE algorithm for Heisenberg model and Bose Hubbard model with interaction

We propose universal SSE method for simulation of Heisenberg model with arbitrary spin and Bose Hubbard model with interaction. We report on the first calculations of soft-core bosons with interaction by the SSE method. Moreover we develop a simple procedure for increase efficiency of the algorithm. From calculation of integrated autocorrelation times we conclude that the method is efficient for both models and essentially eliminates the critical slowing down problem.