I've a charge distribution, whose potential will be given as $$\begin{aligned} V(\mathbf{r})=& \frac{1}{4 \pi \epsilon_{0}}\left[\frac{1}{r} \int \rho\left(\mathbf{r}^{\prime}\right) d \tau^{\prime}+\frac{1}{r^{2}} \int r^{\prime} \cos \alpha \rho\left(\mathbf{r}^{\prime}\right) d \tau^{\prime}\right.\\ &+\frac{1}{r^{3}} \int\left(r^{\prime}\right)^{2}\left(\frac{3}{2} \cos ^{2} \alpha-\frac{1}{2}\right) \rho\left(\mathbf{r}^{\prime}\right) d \tau^{\prime}+\ldots \end{aligned}$$
If the above charge distribution is replaced by a dipole having two discrete charges distanced by infinitesimal length the potential reduces in the infinitesimal length limit to $$V_{\mathrm{dip}}(\mathbf{r})=\frac{1}{4 \pi \epsilon_{0}} \frac{\mathbf{p} \cdot \hat{\mathbf{r}}}{r^{2}}$$
Is there a way in which I can have a charge distribution however small, whose potential also reduces to the above limit? Can an infinitesimal volume of charge ( neutral overall so that the monople term is zero) also have the same potential given by the above equation by arranging the distribution such that all the higher terms in the multipole expansion vanish and all that's left is the dipole term?