“…At the start, we repeated the original local optimizations of the pure ZnO polytypes on ab initio level, using the CRYSTAL14 code, in order to establish quantitative compatibility with the calculations performed on ZnS and the mixed ZnO/ZnS compounds, for both LDA and HSE functionals. As in the earlier work, we found that the wurtzite polymorph of zinc oxide is the energetically lowest and thermodynamically most stable one, which is in agreement with experimental observations (Liu et al, 2013;Sowa & Ahsbahs, 2006;Kisi & Elcombe, 1989;Bates et al, 1962;Ves et al, 1990;Ö zgü r et al, 2005) and previous calculations (Catlow et al, 2008;Ostapenko et al, 2016;Boulfelfel & Leoni, 2008;Mora-Fonz et al, 2017;Sponza et al, 2015;Jaffe et al, 2000;Wang & Zhang, 2016;Demiroglu et al, 2014;Rasoulkhani et al, 2017;Lacerda & de Lazaro, 2017). As shown in Table S2 in the supporting information, when using LDA, all of the predicted ZnO polytypes (4H, 5H, 6H, 8H, 9R, 12R and 15R), as well as the sphalerite structure, exhibit essentially the same energy as the wurtzite modification.…”