Numerical study of the 6-vertex model with domain wall boundary   conditions
  A Markov process is constructed to numerically study the phase separation in the 6-vertex model with domain wall boundary conditions. It is a random walk on the graph where vertices are states and edges are elementary moves. It converges to the Gibbs measure of the 6-vertex model. Our results show clearly that a droplet of "c" vertices is created when Boltzamnn weights are in the antisegnetoelectric region. The droplet is a diamond-like shaped curve with four cusps. 
David Allison and Nicolai Reshetikhin