Fast and Scalable GPU-Accelerated Quantum Chemistry for Periodic Systems with Gaussian Orbitals: Implementation and Hybrid Density Functional Theory Calculations

Efficient hybrid DFT simulations of solid state materials would be extremely beneficial for computational chemistry and materials science, but is presently bottlenecked by difficulties in computing Hartree-Fock (HF) exchange with plane wave orbital bases. We present a GPU-accelerated, Gaussian orbital based integral algorithm for systems with periodic boundary conditions, which takes advantage of Ewald summation to efficiently compute electrostatic interactions. We have implemented this approach into the TeraChem software package within the $Γ$ point approximation, enabling simulation of unit cells with hundreds or thousands of atoms at the HF or hybrid DFT level on a single GPU card. Our implementation readily parallelizes over multiple GPUs and paves the road to accurate simulation of the properties and dynamics of extended materials in both the ground and excited states.