Fast algorithm for periodic density fitting for Bloch waves

Abstract

We propose an efficient algorithm for density fitting of Bloch waves for Hamiltonian operators with periodic potential. The algorithm is based on column selection and random Fourier projection of the orbital functions. The computational cost of the algorithm scales as $\mathcal{O}(N_{\text{grid}}{N}^2+N_{\text{grid}}NK\log (NK))$, where $N_{\text{grid}}$ is number of spatial grid points, $K$ is the number of sampling $k$-points in first Brillouin zone, and $N$ is the number of bands under consideration. We validate the algorithm by numerical examples in both two and three dimensions.

References in corpus (2)

• Resolution-of-identity approach to Hartree-Fock, hybrid density functionals, RPA, MP2, and \textit{GW} with numeric atom-centered orbital basis functions
• The tensor hypercontracted parametric reduced density matrix algorithm: coupled-cluster accuracy with O(r^4) scaling