The application of the SCAN density functional to color centers in diamond

Detailed characterization of deep-level color centers requires understanding their electronic and atomic structure, which is most commonly investigated utilizing the Kohn-Sham density functional theory. Standard semilocal functionals based on the generalized gradient approximation (GGA) are inclined toward an imprecise quantitative description of defects' electronic structure. Hybrid functionals provide an improved prediction of electronic properties, albeit at a much higher computational cost. In this work, we test the newly developed Strongly Constrained and Appropriately Normed (SCAN) family of meta-GGA density functionals for selected color centers in diamond. In particular, we study nitrogen-, silicon-, germanium-, and tin-vacancy centers that have been recently investigated for their use in quantum technological applications. We show that SCAN and its derivatives, the rSCAN and r2SCAN functionals, significantly improve the calculated energies of optical transitions within the delta-self-consistent-field approach, almost reaching the accuracy of the hybrid Heyd-Scuseria-Ernzerhof (HSE) functional. In the case of the NV- center, we also show that the SCAN family of functionals improves the description of the adiabatic potential energy surfaces compared to both GGA and hybrid functionals, improving calculated luminescence lineshapes. As a result of these findings, we recommend using the SCAN family of functionals as a promising alternative for studying color centers in solids.