Pawpyseed: Perturbation-extrapolation band shifting corrections for point defect calculations

Abstract: Significant progress has been made recently in the automation and standardization of ab initio point defect calculations. However, the task of developing, implementing, and benchmarking charge corrections for density functional theory (DFT) point defect calculations is still an open challenge. Here we present a high-performance Python package called pawpyseed, which can read PAW DFT wave functions and calculate the overlap between wavefunctions from different structures. Using pawpyseed, we implement a new band shifting correction derived from first order perturbation theory. We benchmark this method by calculating the transition levels of several point defects in silicon and comparing to experimental and hybrid functional results. The new band shifting method can shift single-particle energies to improve transition level predictions and can be automated and parallelized using pawpyseed, suggesting it could be a useful method for high-throughput point defect calculations.