By using a Morse function and a Witten deformation argument, we obtain an upper bound for the dimensions of the space of divergence-free symmetric Killing p-tensors on a closed Riemannian manifold and explicitly calculate it for p=2.
Choi et al. (Tue,) studied this question.