Harmonic cubic homogeneous polynomials such that the norm-squared of the Hessian is a multiple of the Euclidean quadratic form

There is considered the problem of describing up to linear conformal equivalence those harmonic cubic homogeneous polynomials for which the squared-norm of the Hessian is a nonzero multiple of the quadratic form defining the Euclidean metric. Solutions are constructed in all dimensions and solutions are classified in dimension at most $4$. Techniques are given for determining when two solutions are linearly conformally inequivalent.