Positive harmonic function on $\mathbb{R}^n$ is a constant? Is it true that a positive harmonic function on $\mathbb{R}^n$ must be a constant? How might we show this? The mean value property seems not to be the way...for that we would need boundedness.
 A: Since it is not clear whether the Wikipedia proof uses boundedness or not, please allow me to give a detailed proof that only uses nonnegativity.
Let $u$ be a nonnegative harmonic function in $\mathbb{R}^n$, and let $x,y\in\mathbb{R}^n$. Denote by $B_R(y)$ the ball of radius $R>0$ centred at $y$, and similarly by $B_r(x)$ the ball of radius $r>0$ centred at $x$. We consider such balls with $R>0$ large and $r=R-|x-y|$, that is, $B_r(x)\subset B_R(y)$ and the boundary of $B_r(x)$ is tangent to the boundary of $B_R(y)$.
Then we have
$$
u(x) = \frac1{|B_r|}\int_{B_r(x)}u \leq \frac1{|B_r|}\int_{B_R(y)}u = \frac{|B_R|}{|B_r|}u(y),
$$
where we have used the mean value property in the first and last equalities, and the nonnegativity of $u$ in the second inequality. Now we send $R\to\infty$ and get 
$$
u(x) \leq u(y).
$$
As $x$ and $y$ were arbitrary, we conclude that $u$ is constant.
Note that the Harnack inequality is lurking behind this proof. Also, it was crucial to consider balls with different radii in order to make use of nonnegativity. The Wikipedia proof uses balls of equal radii.
A: Taking $n=2$ just to simplify:
Note that, given $y\in\mathbb{R}^2$, let $r=||y||$. Then, $\forall n\in\mathbb{N}$:
$$B(y,n)\subset B(0,2n)$$
Thus
$$\int_{B(y,n)}u(z)dz\leq\int_{B(0,n+r)}u(z)dz$$
$$u(y) = \frac{1}{\pi n^2}\int_{B(y,n)}u(z)dz\leq\frac{(n+r)^2}{n^2}\frac{1}{\pi(n+r)^2}\int_{B(0,n+r)}u(z)dz=\left(\frac{n+r}{n}\right)^2u(0)$$
That is 
$$u(y)\leq\left(\frac{n+r}{n}\right)^2u(0),\qquad\forall n \in\mathbb{N}$$
Takin $lim_n\rightarrow \infty$, we have that $u(y)\leq u(0)$. But since $y\in\mathbb{R}^2$ was arbitrary, we conclude that $u(z)\leq u(0),\forall z\in\mathbb{R}^2$.
But then, given $R>0$, we have that $\max_{z\in\overline{B(0,R)}}u(z) \leq u(0)$, thus, $\max_{z\in\overline{B(0,R)}}u(z) = u(0)$. Since $0$ is an interior point, by the strong maximum principle, $u$ is constant equal to $u(0)$ in the ball. Since $R$ was arbitrary, it follows that $u(z) = u(0),\forall z\in\mathbb{R}^2$.
A: One can use Louisville's Theorem for harmonic functions to prove this. Note that if $u$ is a harmonic function then so is $e^{-u}$ and if $u$ is bounded below, then $e^{-u}$ is bounded and thus constant. This then forces $u$ to be constant. This shows the slightly stronger fact that if $u$ is a harmonic function bounded below, then $u$ is constant. 
