For calculation of the band alignment of surfaces, the slabs were cleaved from the HSE06-relaxed bulk conventional unit cell. The atomic positions were then relaxed with HSE06 using the same force criterion as before, with constraints on unit cell shape and volume. This was done as hybrid functionals offer a better description of the electronic structure compared to GGA functionals, so the band alignment presented is more accurate.