Cl-terminated models. The primitive cell was used for the calculations of bulk crystalline Si. The DOS calculations used an energy smearing of 0.15 eV and 100 sampling points per eV. We note that the absolute values of the DFT calculated band gaps may be underestimated compared to experiment which is often observed; however, the trends in the band gap are expected to be modelled well.