Additionally, a quasiparticle calculation was performed using the "hybrid quasiparticle self-consistent GW" (hQSGW) method as implemented in the questaal code using a linearised muffin-tin orbital (LMTO) basis set. In this hybrid approach a converged self-energy is obtained by the quasiparticle self-consistent GW (QSGW) approach, and an empirical 80% of the self-energy is combined with the LDA self-energy to obtain a prediction of the semiconductor bandgap. This correction is routinely used to account for neglected interactions and improve the accuracy of bandgap estimations. The lattice parameters and atomic positions were used from the results of the PBE0 geometry optimisation. The LMTO basis set was generated with the recommended parameters in questaal, using a 10 Ry 1/2 cut-off for the interstitial mesh in the one-particle Hamiltonian steps, 2.7 Ry 1/2 cut-off for the interstitial mesh of the two-particle objects and 3.3 Ry 1/2 cut-off for the basis envelope functions in GW steps. The 7 × 7 × 4 k-point grid used for the DFT calculations was reduced to 3 × 3 × 2 for the GW steps. Effective masses were determined by quadratic fitting to band structure data, using tools in the questaal package.