If we’d consider a general shape, then the possibility of the vectors corresponding to r_< and r_> could swap at any point, meaning an expansion would be very complicated. Also, the interactions between the nucleus and the electron cloud would be very hard to compute analytically (and/or numerically) since we aren’t able to split up the interaction in different ‘levels’ (monopole, dipole etc.)