Efficient Implementation of the Random Phase Approximation with Domain-based Local Pair Natural Orbitals

Yu Hsuan Liang,Xing Zhang,Garnet Kin-Lic Chan,Timothy C. Berkelbach,Hong-Zhou Ye
2024-11-12
Abstract:We present an efficient implementation of the random phase approximation (RPA) for molecular systems within the domain-based local pair natural orbital (DLPNO) framework. With optimized parameters, DLPNO-RPA achieves approximately 99.9% accuracy in the total correlation energy compared to a canonical implementation, enabling highly accurate reaction energies and potential energy surfaces to be computed while substantially reducing computational costs. As an application, we demonstrate the capability of DLPNO-RPA to efficiently calculate basis set-converged binding energies for a set of large molecules, with results showing excellent agreement with high-level reference data from both coupled cluster and diffusion Monte Carlo. This development paves the way for the routine use of RPA-based methods in molecular quantum chemistry.
Chemical Physics
What problem does this paper attempt to address?
### What problem does this paper attempt to solve? This paper aims to solve the problem of low computational efficiency of the random - phase approximation (RPA) in large - molecule systems. Specifically, RPA is a powerful method for incorporating non - local multi - electron correlation effects within the Kohn - Sham density functional theory (DFT) framework. It is especially good at capturing long - range dispersion interactions and accurately modeling metallic systems. However, for systems containing more than 100 atoms, performing conventional RPA calculations in the complete basis set (CBS) limit still has a very high computational cost. To solve this problem, the authors proposed an efficient RPA implementation method (DLPNO - RPA) based on the domain - based local pair natural orbital (DLPNO) framework. Through this method, the authors were able to significantly reduce the computational cost while maintaining approximately 99.9% of the total correlation energy accuracy, enabling efficient calculation of reaction energies and potential energy surfaces. In addition, the authors also demonstrated the application of DLPNO - RPA in calculating the basis - set - converged binding energies of large molecules and proved that its results are in good agreement with high - precision reference data (such as the coupled - cluster and diffusion Monte Carlo methods). ### Main contributions 1. **Efficient implementation**: By introducing the DLPNO framework, the computational complexity of RPA calculations in large - molecule systems is greatly reduced, from the traditional O(N^4) to near - linear or sub - linear complexity. 2. **High precision**: The DLPNO - RPA method achieves 99.9% accuracy in the total correlation energy and also shows extremely high accuracy in the calculation of reaction energies, potential energy surfaces, and weak interactions. 3. **Practical applications**: Demonstrated the efficiency and accuracy of DLPNO - RPA in calculating the binding energies of large molecules, especially under basis - set - converged conditions, and the results are highly consistent with advanced correlation wave - function methods. 4. **Extended scope of application**: This work paves the way for the wide application of the RPA method in molecular quantum chemistry, especially for applications involving weak interactions (such as molecular crystals, porous materials, and solid surfaces, etc.). ### Summary The main purpose of this paper is to develop an efficient DLPNO - RPA method to overcome the problem of high computational cost in traditional RPA calculations in large - molecule systems while maintaining high precision. By optimizing parameters and introducing the DLPNO framework, the authors successfully achieved this goal and demonstrated the potential of this method in practical applications.