2D Ising model with Glauber (heat-bath) dynamics. Watch phase transitions, domain coarsening, critical slowing down at T_c, and measure dynamic critical exponent z ≈ 2.17 via autocorrelation decay.
Glauber dynamics: flip spin i with rate W_i = 1/(1+exp(2β·ΔE)). Satisfies detailed balance → samples from Gibbs distribution. Critical temperature T_c = 2J/ln(1+√2) ≈ 2.269 (Onsager 1944). Dynamic critical exponent: τ_corr ~ ξ^z with z ≈ 2.17 (2D Ising). Critical slowing down: correlation time diverges at T_c. Coarsening for T