An extensive theoretical investigation of the nonpolar $\left(101\mathrm{}¯0\right)$ and $\left(112\mathrm{}¯0\right)$ surfaces as well as the polar zinc-terminated (0001)-Zn and oxygen-terminated $\left(0001\mathrm{}¯\right)$-O surfaces of ZnO is presented. Particular attention is given to the convergence properties of various parameters such as basis set, k-point mesh, slab thickness, or relaxation constraints within local-density and generalized-gradient approximation pseudopotential calculations using both plane-wave and mixed-basis sets. The pros and cons of different approaches to deal with the stability problem of the polar surfaces are discussed. Reliable results for the structural relaxations and the energetics of these surfaces are presented and compared to previous theoretical and experimental data, which are also concisely reviewed and commented.

• Received 27 June 2002
• Publisher error corrected 15 January 2003

DOI:https://doi.org/10.1103/PhysRevB.67.035403