What we usually call "fundamental solution of the Laplace operator" is the following function defined on $\mathbb{R}^n\setminus\{0\}$: $$\tag{1}\Phi(x)=\begin{cases} \frac{-1}{2\pi} \log r & n=2 \\ \frac{-1}{(2-n)n\alpha(n)} \frac{1}{r^{n-2}} & n \ge 3 \end{cases}$$ where $r=\lvert x \rvert$ is the radial coordinate. The main property of this function, which justifies its name, is that it is a distributional solution of the equation $$\tag{2}-\Delta \Phi = \delta, $$ but it certainly is not the unique solution of such equation. Indeed the solutions of (2) are exactly the distributions $T$ such that $T=\Phi + h$ for a harmonic (entire) function $h$. In this sense there are "a lot" of fundamental solutions of the Laplace operator.
Question. Why do we always choose the one given in equation (1)?
A first reason that comes to mind is that the only fundamental solutions having radial symmetry are of the form $$E=\Phi+ C, $$ for a constant $C$. But then, why do we choose to set $C=0$? Is this choice purely cosmetic? Especially in dimension $2$, where the fundamental solution is unbounded both at $0$ and at $\infty$, this seems to me to be the case.
