Surface energy and surface proton order of the ice Ih basal and prism surfaces

Ding Pan, Li-Min Liu, Gareth A. Tribello, Ben Slater, Angelos Michaelides, Enge Wang

Research output: Contribution to journalArticlepeer-review

43 Citations (Scopus)


Density-functional theory (DFT) is used to examine the basal and prism surfaces of ice Ih. Similar surface energies are obtained for the two surfaces; however, in each case a strong dependence of the surface energy on surface proton order is identified. This dependence, which can be as much as 50% of the absolute surface energy, is significantly larger than the bulk dependence (< 1%) on proton order, suggesting that the thermodynamic ground state of the ice surface will remain proton ordered well above the bulk order-disorder temperature of about 72 K. On the basal surface this suggestion is supported by Monte Carlo simulations with an empirical potential and solution of a 2D Ising model with nearest neighbor interactions taken from DFT. Order parameters that define the surface energy of each surface in terms of nearest neighbor interactions between dangling OH bonds (those which point out of the surface into vacuum) have been identified and are discussed. Overall, these results suggest that proton order-disorder effects have a profound impact on the stability of ice surfaces and will most likely have an effect on ice surface reactivity as well as ice crystal growth and morphology. S Supplementary data are available from

Original languageEnglish
Article number074209
Pages (from-to)-
Number of pages11
JournalJournal of Physics: Condensed Matter
Issue number7
Publication statusPublished - 24 Feb 2010

Cite this