Thermodynamic Properties of Ice: A Monte Carlo Study

Abstract: In this text, we implement a monte carlo algorithm to study thermodynamic properties of ice. Our program, written in Python, is open-sourced and available at https://github.com/AKnightWing/ColdAsIce. We develop a novel scheme to compute the residual entropy of a two dimensional lattice model of ice at absolute zero. A model of energetic ice is also considered and its thermodynamic properties are studied with temperature. We report a remarkable result for the energetic ice model, the presence of a phase transition at non zero temperatures.