Density-functional study of the structure and stability of ZnO surfaces

An extensive theoretical investigation of the nonpolar (10 (1) over bar0) and (11 (2) over bar0) surfaces as well as the polar zinc-terminated (0001)-Zn and oxygen-terminated (000 (1) over bar)-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.