A dataset containing inputs and outputs for all DFT calculations supporting this study is available under the CC-BY-4.0 license from the University of Bath Research Data Archive . Python codes for calculating defect formation energies and defect concentrations as functions of elemental chemical potentials, and Jupyter notebooks used to generate Figures 1 to 8 are available under the MIT license . Our analysis codes use the scipy , numpy , pandas matplotlib and py-sc-Fermi Python packages. Chemical potential limits stability analysis used the fortran package cplap .