Plz Help: How to calculate magnetization in SSE method?? I need help from someone experienced with Stochastic Series Expansion algorithm. I'm working on SSE algorithm for lattice of isotropic Heisenberg spins and it works fine, well, mostly. I have nice values for energy and heat capacity, but i have problem with magnetization. Since magnetization is diagonal in basis used in SSE method I have simply summed up spins for every configuration generated by algorithm but this gives very, very wrong results. Obviously something different must be done. Can anyone help, please? I really need it.