General formula evaluation of the electron-repulsion integrals and the first and second integral derivatives over Gaussian-type orbitals

Recurrence relations between auxiliary functions to calculate the electron-repulsion integrals (ERIs) over Cartesian Gaussian-type orbitals (GTOs) can be derived. With the use of the Rys polynomials, this serves as a general formula method for evaluating ERIs and the ERI analytical derivatives. The computer code of the present method is quite general (i.e., it has no limitation for the shell type of GTOs), very compact (i.e., it has no coding difficulty), and extremely rapid (e.g., 2.1 μs per one ERI, 0.51 μs per one ERI first derivative, and 0.40 μs per one ERI second derivative for the [II‖II] shell block). It is found that the present method is especially efficient for evaluating ERIs including g-type or higher GTOs, the first derivatives including f-type or higher GTOs, and the second derivatives including d-type or higher GTOs.