This paper deals with the problem of factorizing integer powers of the
Laplace operator acting on functions taking values in higher spin
representations. This is a far-reaching generalization of the well-known fact
that the square of the Dirac operator is equal to the Laplace operator. Using
algebraic properties of projections of Stein-Weiss gradients, i.e. generalized
Rarita-Schwinger and twistor operators, we give a sharp upper bound on the
order of polyharmonicity for functions with values in a given representation
with half-integral highest weight.