Three-center nuclear attraction integrals over exponential-type functions are required for ab initio molecular structure calculations and density functional theory (DFT). These integrals occur in many millions of terms, even for small molecules, and they require rapid and accurate numerical evaluation. The use of a basis set of B functions to represent atomic orbitals, combined with the Fourier transform method, led to the development of analytic expressions for these molecular integrals. Unfortunately, the numerical evaluation of the analytic expressions obtained turned out to be extremely difficult due to the presence of two-dimensional integral representations, involving spherical Bessel integral functions. The present work concerns the development of an extremely accurate and rapid algorithm for the numerical evaluation of these spherical Bessel integrals. This algorithm, which is based on the nonlinear D transformation and the W algorithm of Sidi, can be computed recursively, allowing the control of the degree of accuracy. Numerical analysis tests were performed to further improve the efficiency of our algorithm. The numerical results section demonstrates the efficiency of this new algorithm for the numerical evaluation of three-center nuclear attraction integrals.

In a previous work, a new Gauss quadrature was introduced with a view to evaluate multicenter integrals over Slater-type functions efficiently. The complexity analysis of the new approach, carried out using the three-center nuclear integral as a case study, has shown that for low-order polynomials its efficiency is comparable to the SD. The latter was developed in connection with multi-center integrals evaluated by means of the Fourier transform of B functions. In this work we investigate the numerical properties of the Gauss-Bessel quadrature and devise strategies for an efficient implementation of the numerical algorithms for the evaluation of multi-center integrals in the framework of the Gaussian transform/Gauss-Bessel approach. The success of these strategies are essential to elaborate a fast and reliable algorithm for the evaluation of multi-center integrals over STFs.