The S and G transformations for computing three‐center nuclear attraction integrals

Abstract

AbstractIt is now well established that nonlinear transformations can be extremely useful in the case of oscillatory integrals. In previous work, we could show that the G transformation, which is not so well known among those interested in the numerical evaluation of highly oscillatory integrals, works very well for the extremely challenging integral called Twisted Tail. In this work, we demonstrate that these techniques also apply to three‐center nuclear attraction integrals over exponential type functions. The accurate and rapid evaluation of these integrals is required in ab initio molecular structure calculations and density functional theory. Using a basis set of B functions and profiting from their relatively simple Fourier representation, these integrals are formulated as analytical expressions involving highly oscillatory spherical Bessel integral functions. In the present work, we implement two highly accurate algorithms for three‐center nuclear attraction integrals. The first algorithm is based on the G transformation and the second is based on a combination of the S and G transformations. The application of these transformations is largely due to the properties of special functions that allow the computation of higher order derivatives of the integrands with exceptional simplicity. The numerical results illustrate the accuracy of these algorithms applied to three‐center nuclear attraction integrals over exponential type functions with a miscellany of different parameters. © 2009 Wiley Periodicals, Inc. Int J Quantum Chem, 2009