Resolutions of the Coulomb operator: VI. Computation of auxiliary   integrals
  We discuss the efficient computation of the auxiliary integrals that arise when resolutions of two-electron operators (specifically, the Coulomb and long-range Ewald operators) are employed in quantum chemical calculations. We derive a recurrence relation that facilitates the generation of auxiliary integrals for Gaussian basis functions of arbitrary angular momentum and propose a near-optimal algorithm for its use. 
Taweetham Limpanuparb, Joshua W. Hollett, Peter M. W. Gill
physics.chem-ph