First principles correction scheme for linear-response time-dependent density functional theory calculations of core electronic states

Augustin Bussy,Jürg Hutter
DOI: https://doi.org/10.1063/5.0058124
2021-05-26
Abstract:Linear-response time-dependent density functional theory (LR-TDDFT) for core level spectroscopy using standard local functionals suffers from self-interaction error and a lack of orbital relaxation upon creation of the core hole. As a result, LR-TDDFT calculated X-ray absorption near edge structure (XANES) spectra need to be shifted along the energy axis to match experimental data. We propose a correction scheme based on many body perturbation theory to calculate the shift from first principles. The ionization potential of the core donor state is first computed and then substituted for the corresponding Kohn--Sham orbital energy, thus emulating Koopmans' condition. Both self-interaction error and orbital relaxation are taken into account. The method exploits the localized nature of core states for efficiency and integrates seamlessly in our previous implementation of core level LR-TDDFT, yielding corrected spectra in a single calculation. We benchmark the correction scheme on molecules at the K- and L-edges as well as for core binding energies and report accuracies comparable to higher order methods. We also demonstrate applicability in large and extended systems and discuss efficient approximations.
Chemical Physics
What problem does this paper attempt to address?
This paper aims to solve two main problems encountered in the application of linear - response time - dependent density functional theory (LR - TDDFT) in core - level spectroscopy: self - interaction error (SIE) and insufficient orbital relaxation. These problems lead to the need to shift the X - ray near - edge absorption structure (XANES) spectra calculated by LR - TDDFT along the energy axis to match the experimental data. To solve these problems, the authors propose an ab - initio correction scheme based on many - body perturbation theory. This scheme can calculate the correction value and apply it to LR - TDDFT calculations, thereby improving the accuracy of the calculation results. Specifically, this method first calculates the ionization potential of the core - donor state, and then replaces the corresponding Kohn - Sham orbital energy with this value to simulate the Koopmans condition. This method not only takes into account self - interaction error and orbital relaxation, but also utilizes the local properties of the core state to improve computational efficiency, and can be seamlessly integrated into the existing core - level LR - TDDFT implementation, so as to obtain the corrected spectra in a single calculation. Through benchmark tests on the K - and L - edges of molecules and core - binding energies, this paper demonstrates the effectiveness of this correction scheme, and its accuracy is comparable to that of higher - order methods. In addition, the applicability of this method in large and extended systems has also been verified, and effective approximation methods are discussed.