An improvement of the resolution of the identity approximation for the formation of the Coulomb matrix

A straightforward modification of the resolution of the identity (RI) approximation to the Coulomb interaction is described. In the limit of basis sets that are dominated by high angular momentum functions the observed speedups in realistic test systems reach a factor of 2 compared to the standard RI algorithm, and a factor of up to 300 compared to the standard algorithm to form the Coulomb matrix. More moderate savings on the order of 0-20% are obtained for the more commonly used smaller basis sets. A series of test calculations is reported to illustrate the efficiency of the algorithm.