 Solution

The conclusion holds in dimension two. Shift and change sign to obtain a nonnegative harmonic function $v$. Its circular mean satisfiesHere the angular second derivative in the polar-coordinate Laplacian integrates to zero. Thus $m(r)=a\log r+b$; nonnegativity for every $r>0$ forces $a=0$. For $|x|=r$, the mean value property for harmonic functions and nonnegativity giveThe removable singularity for a bounded harmonic function fills in the origin, and the harmonic Liouville theorem makes the extension constant.

The conclusion fails for $n\ge 3$: $\boxed{u(x)=|x{|}^{2-n}}$ is positive and nonconstant, and its radial Laplacian $u^{\prime \prime }+(n-1)u^{\prime }/r$ vanishes away from zero. For $n=1$, $\boxed{u(x)=|x|}$ is a nonconstant function bounded below whose second derivative vanishes on both components of the punctured line.