q-state Potts model: H = -J Σ δ(s_i, s_j). Each spin takes q values. Metropolis dynamics. For q≤4: continuous (2nd order) transition. For q>4: first-order (discontinuous) transition with latent heat. Domain walls form between different-state regions.