Electronic and optical properties of Fe, Zn and Pb sulfides

Abstract

Ab initio quantum-chemical calculations of the spatial and electronic structures of sphalerite (ZnS), pyrite (FeS2) and galena (PbS), using the density functional theory (DFT) local density approximation (LDA) and generalized gradient approximation (GGA), the Hartree–Fock (HF) method and the hybrid functional B3LYP, have been carried out. For galena, the DFT LDA and GGA functionals provided the best estimate of the band gap, from within −0.1 eV to +0.4 eV of the measured value. B3LYP and RHF gave rise to errors of +1.3 and +5.4 eV, respectively. The unit cell parameter error varied from between −1.1% and +2.3% for all the functionals examined. For sphalerite the B3LYP functional provided the best estimate of the band gap (error +0.3 eV). The unit cell parameter error varied between −2.1% and +2.0% for the various DFT functionals and B3LYP. RHF gave rise to an error of +3.8%. For FeS2, the DFT-GGA approach provides the best results for both the unit cell and the band gap. This may be due to mutual cancellation of the crystal field splitting and band separation force, which are of equal but opposite magnitudes. The calculated density of states (DOS) for the conduction band is used to interpret the experimental features of the S 1s XANES (X-ray absorption near-edge structure) spectra obtained using synchrotron radiation. Because of the Δl = ±1 selection rule for electron excitation, the S K-edge XANES spectra represent a transition of the S 1s electron to conduction band S p-like orbitals. The near-edge region, up to 15 eV past the edge is approximated well by the DOS. Individual peaks in the DOS correlate with peaks in the XANES spectra. In addition, the imaginary part of the dielectric function, which reflects the transitions from occupied to unoccupied levels, is used to model the near-edge region of the XANES, using the DFT-GGA formalism. Individual peaks in the XANES spectrum are moderately well resolved using the dielectric function, especially for ZnS and FeS2, while the DOS for the conduction band is more successful in predicting the shape of the XANES spectra for all three minerals.