to account for Pulay stress. Because VASP cannot restrain the periodicity into just two dimensions, the slab model was employed in this study to study the surface properties. The surfaxe package was used to cleave the slabs, calculate the parameters used in convergence testing and perform the relevant analyses. All nonpolar, symmetric slabs up to a maximum Miller index of two were cleaved from the relaxed conventional unit cell.