“…We note that, previously, surface structures have already been a target of a number of Monte Carlo studies using a canonical ensemble, which, however, have been severely limited in the scope, accuracy, and comparability with realistic systems. A more general approach employing the concept of a grand potential has been used in thermodynamic analysis of (i) phase segregation for nonstoichiometric surface terminations of various materials (for example, see refs − ), (ii) nonstoichiometric slabs, (iii) interfaces, , and (iv) surface reconstructions. − Following these thermodynamic approaches, a local or global deviation from precise 1:1 stoichiometry at surfaces of ZnO can be trivially incorporated using chemical potentials of O and Zn atoms or ions that are simply added to the total energy of preoptimized surface potential energies. Such an approach is particularly beneficial for modeling ZnO growth that in solutions would involve transient nonstoichiometric (and, possibly, high-energy) configurations, with a general addition (or removal) of Zn x O y at each elementary step cf.…”