An electrostatic potential associated with some delocalized charge $\int \rho(\mathbf{r}) d{\mathbf{r}}$ is given by:
$$v_H(\mathbf{r}) = \int \frac{\rho(\mathbf{r'})}{|\mathbf{r}-\mathbf{r'}|}d\mathbf{r'}$$
This potential is finite at $\mathbf{r}=0$. Since $\frac{1}{|\mathbf{r}-\mathbf{r'}|}$ is a solution to a singular Poisson's equation, we can show that:
$$\nabla^2v_H(\mathbf{r})=-4\pi\rho(\mathbf{r})$$ where $\rho$ is a smooth function being finite everywhere.
Would a function
$$v(\mathbf{r}) = \int \frac{\rho(\mathbf{r'})}{|\mathbf{r}-\mathbf{r'}|^n}d\mathbf{r'}$$
(where $n$ is some nonnegative integer)
be also finite at $\mathbf{r}=0$?
In this case, $\frac{1}{|\mathbf{r}-\mathbf{r'}|^n}$ does not represent a solution to Poisson's equation ($\mathbf{r} \in \mathbb{R}^3$) and the above analysis is not valid.