implementation based on the pseudopotential and numerical-atomic-orbital basis-set framework: Algorithms and benchmarks
Phys. Rev. B 114, 175139 – Published 30 September, 2026
DOI: https://doi.org/10.1103/47j8-8kll
Abstract
The method delivers substantially improved accuracy in electronic band structure calculations over conventional Kohn–Sham density functional theory (KS-DFT) by explicitly incorporating the electron self-energy effect beyond mean-field approximations. In this work, we present an efficient NAO-PP-based computational framework by interfacing the first-principles software package abacus with librpa, a library for performing low-scaling random-phase approximation and calculations based on NAOs. Our approach employs the localized resolution of identity (LRI) technique with a compression scheme, significantly improving both computational efficiency and numerical stability. In addition, an analytic treatment of the small- limit of the microscopic dielectric function reduces the need for dense -point sampling. Furthermore, we quantify how the choice of explicitly treated valence electrons affects the frequency-dependent macroscopic dielectric response and the resulting quasiparticle energies. Systematic benchmarks validate the effectiveness of our compression scheme and real-space tensor filtering strategies, demonstrating both high accuracy and significant computational efficiency gains. Comparisons with established implementations show excellent agreement in band structures and band gaps, confirming as a reliable and efficient platform for large-scale simulations.