Quartic scaling MP2 for solids: A highly parallelized algorithm in the plane wave basis

TL;DR

This paper introduces a highly parallelizable MP2 algorithm for solids with quartic scaling, utilizing Laplace transforms and plane wave basis sets to efficiently compute correlation energies in periodic systems.

Contribution

The paper presents a novel quartic-scaling MP2 algorithm for solids that leverages Laplace transforms and plane wave basis sets, improving computational efficiency and parallelization.

Findings

Achieves $ ext{O}(N^4)$ scaling for MP2 calculations in solids.

Demonstrates high parallelization efficiency.

Provides a basis set extrapolation method to improve convergence.

Abstract

We present a low-complexity algorithm to calculate the correlation energy of periodic systems in second-order M\o ller-Plesset perturbation theory (MP2). In contrast to previous approximation-free MP2 codes, our implementation possesses a quartic scaling, $\mathcal O (N^4)$, with respect to the system size $N$ and offers an almost ideal parallelization efficiency. The general issue that the correlation energy converges slowly with the number of basis functions is eased by an internal basis set extrapolation. The key concept to reduce the scaling is to eliminate all summations over virtual orbitals which can be elegantly achieved in the Laplace transformed MP2 (LTMP2) formulation using plane wave basis sets and Fast Fourier transforms. Analogously, this approach could allow to calculate second order screened exchange (SOSEX) as well as particle-hole ladder diagrams with a similar low complexity. Hence, the presented method can be considered as a step towards systematically improved correlation energies.