Numerical Integration of Slater Basis Functions Over Prolate Spheroidal Grids
Abstract
Slater basis functions have desirable properties that can improve electronic structure simulations, but improved numerical integration methods are needed. This work builds upon the SlaterGPU library for evaluation of Hamiltonian matrix elements in the resolution-of-the-identity approximation. In particular, a Prolate Spheroidal grid will provide sufficient integral accuracy to employ larger basis sets (quadruple-zeta and greater) in practical computations involving polyatomics. To integrate 3-center Coulomb and nuclear attraction terms, an improved grid representation around the 3rd center is introduced. The RMSEs of the integral quantities are evaluated and compared to the previous numerical integration method used in SlaterGPU (Becke Partitioning), resulting in a ~3 order of magnitude reduction in the error for 2-center integral quantities. The procedure is generally applicable to polyatomic systems, GPU accelerated for high performance computing, and tested on self-consistent field and full configuration interaction wavefunctions. Results for a number of 3-atom models as well as propanediyl (C3H6) demonstrate the reliability of the new integration scheme.
Keywords
Cite
@article{arxiv.2510.15204,
title = {Numerical Integration of Slater Basis Functions Over Prolate Spheroidal Grids},
author = {Alexander Stark and Nathan Meier and Jeffrey Hatch and Joshua Kammeraad and Duy-Khoi Dang and Paul Zimmerman},
journal= {arXiv preprint arXiv:2510.15204},
year = {2025}
}