2D Ising Model — Metropolis Monte Carlo

Ferromagnetic Ising model on a square lattice. Metropolis algorithm samples the Boltzmann distribution. Critical temperature T_c = 2/ln(1+√2) ≈ 2.269 (Onsager 1944).

|M|/N
E/N
2.269T_c
T