“…For positive even values of , the expansion contains a finite number of terms (due to the cosine formula for r ij 2 ), but for positive odd it is an infinite series with the radial function R u (r i , r j ), becoming more complicated as increases. Use of this expansion for integer values of less than Ϫ1 is not needed for traditional nonrelativistic energy calculations, but the expansion of r ij Ϫ2 arises in lower-bound calculations [13,15] and in the treatment of relativistic effects [9,18,19]. Sack's formulation includes the case ϭ Ϫ2, but severe numerical and convergence problems have been encountered with its use for that value of , and several workers have been exploring alternatives or modifications [12,21].…”