ISING MODEL — GLAUBER DYNAMICS

Single-spin flips accepted via Metropolis: P(flip) = min(1, exp(−ΔE/kT)). Tc = 2J/ln(1+√2) ≈ 2.269J