First-principles study of the polar O-terminated ZnO surface in thermodynamic equilibrium with oxygen and hydrogen

Using density-functional theory in combination with a thermodynamic formalism we calculate the relative stability of various structural models of the polar O-terminated (000(1) over bar)-O surface of ZnO. Model surfaces with different concentrations of oxygen vacancies and hydrogen adatoms are considered. Assuming that the surfaces are in thermodynamic equilibrium with an O-2 and H-2 gas phase we determine a phase diagram of the lowest-energy surface structures. For a wide range of temperatures and pressures we find that hydrogen will be adsorbed at the surface, preferentially with a coverage of 1/2 monolayer. At high temperatures and low pressures the hydrogen can be removed and a structure with 1/4 of the surface oxygen atoms missing becomes the most stable one. The clean, defect-free surface can only exist in an oxygen-rich environment with a very low hydrogen partial pressure. However, since we find that the dissociative adsorption of molecular hydrogen and water (if also the Zn-terminated surface is present) is energetically very preferable, it is very unlikely that a clean, defect-free (000(1) over bar)-O surface can be observed in experiment.