The pyrite and marcasite polymorphs of FeS2 have attracted considerable interests for their potential applications in optoelectronic devices because of their appropriate electronic and optical properties. Controversies regarding their fundamental band gaps remain in both experimental and theoretical materials research of FeS2. In this work, we present a systematic theoretical investigation into the electronic band structures of the two polymorphs by using many-body perturbation theory with the GW approximation implemented in the full-potential linearized augmented plane waves (FP-LAPW) framework. By comparing the quasi-particle (QP) band structures computed with the conventional LAPW basis and the one extended by high-energy local orbitals (HLOs), denoted as LAPW + HLOs, we find that one-shot or partially self-consistent GW (G0W0 and GW0, respectively) on top of the Perdew-Burke-Ernzerhof (PBE) generalized gradient approximation with a converged LAPW + HLOs basis is able to remedy the artifact reported in the previous GW calculations, and leads to overall good agreement with experiment for the fundamental band gaps of the two polymorphs. Density of states calculated from G0W0@PBE with the converged LAPW + HLOs basis agrees well with the energy distribution curves from photo-electron spectroscopy for pyrite. We have also investigated the performances of several hybrid functionals, which were previously shown to be able to predict band gaps of many insulating systems with accuracy close or comparable to GW. It is shown that the hybrid functionals considered in general fail badly to describe the band structures of FeS2 polymorphs. This work indicates that accurate prediction of electronic band structure of FeS2 poses a stringent test on state-of-the-art first-principles approaches, and the G0W0 method based on semi-local approximation performs well for this difficult system if it is practiced with well-converged numerical accuracy.
Solubility prediction from first principles: a density of states approach