Calculation of two-centre two-electron integrals over Slater-type   orbitals revisited. I. Coulomb and hybrid integrals
  In this paper, which constitutes the first part of the series, we consider calculation of two-centre Coulomb and hybrid integrals over Slater-type orbitals (STOs). General formulae for these integrals are derived with no restrictions on the values of the quantum numbers and nonlinear parameters. Direct integration over the coordinates of one of the electrons leaves us with the set of overlap-like integrals which are evaluated by using two distinct methods. The first one is based on the transformation to the ellipsoidal coordinates system and the second utilises a recursive scheme for consecutive increase of the angular momenta in the integrand. In both methods simple one-dimensional numerical integrations are used in order to avoid severe digital erosion connected with the straightforward use of the alternative analytical formulae. It is discussed that the numerical integration does not introduce a large computational overhead since the integrands are well-behaved functions, calculated recursively with decent speed. Special attention is paid to the numerical stability of the algorithms. Applicability of the resulting scheme over a large range of the nonlinear parameters is tested on examples of the most difficult integrals appearing in the actual calculations including at most 7i-type functions (l=6). 
Micha{\l} Lesiuk and Robert Moszynski
quant-ph physics.chem-ph