Point defects at the ice (0001) surface

Using density functional theory we investigate whether intrinsic defects in ice surface segregate. We predict that hydronium, hydroxide, and the Bjerrum L- and D-defects are all more stable at the surface. However, the energetic cost to create a D-defect at the surface and migrate it into the bulk crystal is smaller than its bulk formation energy. Absolute and relative segregation energies are sensitive to the surface structure of ice, especially the spatial distribution of protons associated with dangling hydrogen bonds. It is found that the basal plane surface of hexagonal ice increases the bulk concentration of Bjerrum defects, strongly favoring D-defects over L- defects. Dangling protons associated with undercoordinated water molecules are preferentially injected into the crystal bulk as Bjerrum D-defects, leading to a surface dipole that attracts hydronium ions. Aside from the disparity in segregation energies for the Bjerrum defects, we find the interactions between defect species to be very finely balanced; surface segregation energies for hydronium and hydroxide species and trapping energies of these ionic species with Bjerrum defects are equal within the accuracy of our calculations. The mobility of the ionic hydronium and hydroxide species is greatly reduced at the surface in comparison to the bulk due to surface sites with high trapping affinities. We suggest that, in pure ice samples, the surface of ice will have an acidic character due to the presence of hydronium ions. This may be important in understanding the reactivity of ice particulates in the upper atmosphere and at the boundary layer.