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

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 C, where we were originally limited to C due to the high RAM usage of GOSTSHYP.