First-principles calculations of defects in metal halide perovskites: A performance comparison of density functionals

Metal halide perovskite semiconductors have outstanding optoelectronic properties. Although these perovskites are defect-tolerant electronically, defects hamper their long-term stability and cause degradation. Density functional theory (DFT) calculations are an important tool to unravel the microscopic structures of defects, but results suffer from the different approximations used in the DFT functionals. In the case of metal halide perovskites, qualitatively different results have been reported with different functionals, either predicting vacancy or interstitial point defects to be most dominant. Here, we conduct a comprehensive comparison of a wide range of functionals for calculating the equilibrium defect formation energies and concentrations of point defects in the archetype metal halide perovskite, ${\mathrm{MAPbI}}_3$. We find that it is essential to include long-range van der Waals interactions in the functional, and that it is vital to self-consistently optimize structure and volume of all compounds involved in the defect formation. For calculating equilibrium formation energies of point defects in ${\mathrm{MAPbI}}_3$ and similar metal halide perovskites, we argue that the exact values of the chemical potentials of the species involved, or of the intrinsic Fermi level, are not important. In contrast to the simple Schottky or Frenkel pictures, we find that the dominant defects are MA and I interstitials, and Pb vacancies.