An efficient implementation of the GOSTSHYP pressure model by applying shell-bounding Gaussian 1-electron-3-center integral screening

Abstract

We implemented a screening algorithm for one-electron-three-center overlap integrals over contracted Gaussian-type orbitals into the Q-Chem program package. The respective bounds were derived using shell-bounding Gaussians and the Obara–Saika recurrence relations. Using integral screening, we reduced the computational scaling of the Gaussians On Surface Tesserae Simulate HYdrostatic Pressure (GOSTSHYP) model in terms of calculation time and memory usage to a linear relationship with the tesserae used to discretize the surface area. Further code improvements allowed for additional performance boosts. To demonstrate the algorithm’s better performance, we calculated the compressibility of fullerenes up to C180, where we were originally limited to C40 due to the high RAM usage of GOSTSHYP.