This is my first idea for how the swaps should work, with P = min(1, Exp(Beta*hstep*(m2-m1))). It doesn’t work, basically as beta increases, it flattens out the magnetization until everything is zero magnetized. I’m going to try something else next.