It comes sometimes from Laplace equation.
But if such dependency, there are options :
one could use exp(-r), which is much easier in calculations
1/sqrt(r), which is the solution in fracture mechanics, and then a potential is not singular.
That will imply a geometrical mean and splittin of variables and would be more easy to handle in Quantum field theories.