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

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 stacks.iop.org/JPhysCM/22/074209/mmedia

Surface Energy and Surface Proton Order of Ice Ih

Ice Ih is comprised of orientationally disordered water molecules giving rise to positional disorder of the hydrogen atoms in the hydrogen bonded network of the lattice. Here we arrive at a first principles determination of the surface energy of ice Ih and suggest that the surface of ice is significantly more proton ordered than the bulk. We predict that the proton order-disorder transition, which occurs in the bulk at similar to 72 K, will not occur at the surface at any temperature below surface melting. An order parameter which defines the surface energy of ice Ih surfaces is also identified.