Entire $p$-harmonic functions with an isolated critical point and failure of $C^1$-regularity of the natural gradient

For every $n\ge 2$ and $1<p<\mathrm{\infty }$, we show an existence of a nonconstant homogeneous entire $p$-harmonic function whose only critical point is the origin. The angular parts of these functions are axially symmetric and even across the equator. To establish the solution we use contraction argument at the pole in conjuction with a shooting argument. It follows that in every dimension $n\ge 3$ and for $p<2$ sufficiently close to $2$, the natural gradient has pointwise Hölder exponent strictly below one at the origin. Such examples settle the remaining scalar case of the conjecture of Balci, Diening, and Weimar by disproving both the $C^1$ assertion and the linear $L^2$ mean oscillation estimate.