Monte Carlo simulation of ice models

We study a number of Monte Carlo algorithms for the simulation of ice models, and compare their efficiency. One of them, a cluster algorithm for the equivalent three-color model, appears to have a dynamic exponent close to zero, making it particularly useful for simulations of critical ice models. We have performed extensive simulations using our algorithms to determine a number of critical exponents for the square ice and $F$ models.