Why Are Polar Surfaces of ZnO Stable?

We probe and rationalize the complex surface chemistry of wurtzite ZnO by employing interatomic potential calculations coupled with a Monte Carlo procedure that sampled over 0.5 million local minima. We analyze the structure and stability of the (0001) and (0001) ZnO surfaces, rationalizing previous patterns found in STM images and explaining the (1 × 1) periodicity reported by LEED analysis. The full range of Zn/O surface occupancies was covered for a (5 × 5) supercell, keeping |mZn – mO|/N ≈ 0.24 where m and N are the numbers of occupied surface sites and total surface sites, respectively. Our calculations explain why the (5 × 5) reconstructions seen in some experiments and highlight the importance of completely canceling the inherent dipole of the unreconstructed polar surfaces. The experimentally observed rich reconstruction patterns can be traced from the lowest occupancy, showing the thermodynamically most stable configurations of both polar surfaces. Triangular and striped reconstructions are seen...