I am currently inversitgating a system consists of positive and negative point particles satisfying charge neutrality condition. They have a unit charge +1 and -1. In most case, we deal with coulomb interaction potential energy as
$$U = \sum\limits_{i,j} \frac{q_{i}q_{j}}{r_{ij}},$$
where $r_{ij} = |\vec{r}_{i}-\vec{r}_{j}|$. The above equation just considers monopole interaction between charged particles. If I expand $\frac{1}{r_{ij}}$ in terms of spherical harmonics, do you think the spherical harmonics consider multipolar interaction such as dipole and quadrupole between charged particles?