An orthonormal basis for the space of harmonic polynomials in $n$ variables is provided by the spherical harmonics on the $n-1$ sphere, see e.g. wiki.
From there, constructing an orthonormal basis for the subspace of harmonic polynomials symmetric in $n$ variables is straightforward but laborious (or is it?). Are there known closed expressions, or simple derivation, for such a basis?