Three-center expansion of electron repulsion integrals with linear combination of atomic electron distributions

A new systematic way of constructing auxiliary basis functions for approximating the evaluation of electron repulsion integrals is proposed and applied to SCF and MCSCF wavefunction calculations. In the approximation, the one-electron density is expanded in terms of a linear combination of atomic electron distributions (LCAD), and the four-center two-electron repulsion integrals are reduced to the three- and two-center quantities. This results in a high-accuracy approximation as well as a large reduction in disk storage and input/output requirement, proportional to N3 rather than N4, N being the number of basis functions. Numerical results indicate that the error from the present approximation decreases as the size of molecular basis functions increases and that the LCAD version of MCSCF calculations requires only a fractional amount of the CPU time required in the conventional procedure without loss of accuracy.