Multipole expansion of integral powers of cosine theta

Legendre polynomials form the basis for multipole expansion of spatially varying functions. The technique allows for decomposition of the function into two separate parts with one depending on the radial coordinates only and the other depending on the angular variables. In this work, the angular function [Formula: see text] is expanded in the Legendre polynomial basis and the algorithm for determining the corresponding coefficients of the Legendre polynomials is generated. This expansion together with the algorithm can be generalized to any case in which a dot product of any two vectors appears. Two alternative multipole expansions for the electron-electron Coulomb repulsion term are obtained. It is shown that the conventional multipole expansion of the Coulomb repulsion term is a special case for one of the expansions generated in this work.