Eq. 52 suggests that the spherical averaging can be performed before generating the pairs with higher azimuthal quantum numbers through differentiation with respect to the center of charge coordinates. Through the introduction of Hermite Gaussians and the associated recursion scheme for the Hermite-to-Cartesian coefficients, computation of the total coefficients multiplying derivatives of the s-type Hermite Gaussian becomes much simpler in comparison to the previous scheme for the evaluation of the spherically averaged density, .