Rainbow Gauge Theory (Loop Model)
$\beta H = - h \sum_{e} (\sigma_e > 0)$
Edges represent spins, with colored lines representing the nontrivial states. Each Monte Carlo step flips all four edges around a random plaquette, accepted with the Glauber probability 1/(1+exp(ΔE/T)).