It is well known that the self-interaction error inherent to DFT causes uncertainty in the absolute energies of the Kohn-Sham states. Since the single particle states are crucial for this study when reconciling defect charge transitions and the density of states, all point defect calculations employed the hybrid HSE06 functional. We set the exact exchange mixing parameter to 0.37, which yielded an indirect band gap of 0.76 eV in agreement with the experimentally measured band gap of 0.78 eV. The defect formation energy ∆H D q of a defect D with charge state q is calculated using the equation