An overlap fitted chain of spheres exchange method

The "chain of spheres" (COS) algorithm, as part of the RIJCOSX SCF procedure, approximates the exchange term by performing analytic integration with respect to the coordinates of only one of the two electrons, whereas for the remaining coordinates, integration is carried out numerically. In the present work, we attempt to enhance the efficiency of the method by minimizing numerical errors in the COS procedure. The main idea is based on the work of Friesner and consists of finding a fitting matrix, Q, which leads the numerical and analytically evaluated overlap matrices to coincide. Using Q, the evaluation of exchange integrals can indeed be improved. Improved results and timings are obtained with the present default grid setup for both single point calculations and geometry optimizations. The fitting procedure results in a reduction of grid sizes necessary for achieving chemical accuracy. We demonstrate this by testing a number of grids and comparing results to the fully analytic and the earlier COS approximations. This turns out to be favourable for total and reaction energies, for which chemical accuracy can now be reached with a corresponding ~30% speedup over the original RIJCOSX procedure for single point energies. Results are slightly less favourable for the accuracy of geometry optimizations, but the procedure is still shown to yield geometries with errors well below the method inherent errors of the employed theoretical framework.