Factorization of Laplace operators on higher spin representations

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.