All surface slabs studied were nonpolar and contained inversion symmetry. The latter requirement is directly related to the slab model employed as it assumes the two surfaces created are equivalent. The importance of non-polarity is twofold; (1) polar surfaces are typically higher in energy, so less likely to spontaneously form and (2) charged surfaces are difficult to accurately simulate due to the additional dipole corrections required in VASP. When cleaving the surfaces, no restrictions were placed on the bonds broken as none of the bonds in the compound were identified to be strong covalent bonds. In total, one non-polar symmetric termination was identified for each of the eleven Miller indices under investigation. The constructed slabs have bulk-like surfaces, left in a high energetic state due to the large number of dangling bonds produced by the cleavage. To account for this, all atoms were relaxed to find a more favourable, lower energy configuration. All bulk-like and relaxed slabs can be found in the Supporting Information data repository.