Curvature collineation and the Killing identity

Using the identity in here, we have $$ [nabla_{[a}, mathcal{L}_X] R_{bc]de} = pi_{d[a}{}^f R_{bc]fe} + pi_{e[a}{}^f R_{bc]df} tag{1}$$ where we used that the first order deformation tensor is symmetric in its first two indices.

The first order deformation tensor is, explicitly $$ pi_{lm}{}^n = frac12 g^{np} left( nabla_l (nabla_m X_p + nabla_p X_m) + nabla_m(nabla_l X_p + nabla_p X_l) - nabla_p (nabla_m X_l + nabla_l X_m) right) $$ Rearranging we get $$ pi_{lm}{}^n = frac12 g^{np}( 2 nabla_l nabla_m X_p + [nabla_m ,nabla_l ] X_p + [nabla_l, nabla_p]X_m + [nabla_m, nabla_p] X_l ) $$ or $$ pi_{lm}{}^n = frac12 g^{np}( 2 nabla_l nabla_m X_p + (R_{spml} + R_{smlp} + R_{slmp})X^s ) $$ From the first Bianchi identity we finally get $$ pi_{lm}{}^n = g^{np} ( nabla_l nabla_m X_p + R_{slmp} X^s) $$

Notice that the terms in the brackets are exactly the "Killing identity" terms. (I think we use different sign convention for Riemann.)

---

An immediate consequences of the above computation:

Killing vectors satisfy the Killing identity

For Killing vectors, the 0th order deformation tensor vanishes (Killing's equation), and since the 1st order deformation tensor is formed through the covariant derivative of the 0th order ones, it must also vanish. And thus Killing's identity must hold.

---

Now let's see what this allows us to say for curvature collineation. By the second Bianchi identity, $nabla_{[a}R_{bc]de} = 0$. Assuming curvature collineation, the left hand side of equation (1) vanishes. This requires $$ pi_{d[a}{}^fR_{bc]fe} + pi_{e[a}{}^fR_{bc]df} = 0 $$

Let $Gamma$ be the vector space of type $(1,2)$-tensors symmetric in the covariant indices, and define the linear mapping $Upsilon$ from $Gamma$ to type $(0,5)$ tensors by

$$ gamma_{ab}{}^c mapsto gamma_{d[a}{}^fR_{bc]fe} + gamma_{e[a}{}^fR_{bc]df} $$

For your assertion to hold (that curvature collineation implies Killing's identity), it is sufficient that the Riemann curvature tensor is such that the mapping $Upsilon$ is injective. Obviously, as Deane noted, this fails when your manifold is flat.

In low dimensions injectivity however can fail in general. In two dimensions the antisymmetry in the $a, b, c$ indices means that $Upsilon$ is by definition the trivial map. In three dimensions the domain $Gamma$ is 18 dimensional, while the target is 3 dimensional (isomorphic to the space of antisymmetric matrices) and so must have a non-trivial Kernel. Even in four dimensions, the dimensional counting argument is not optimistic: let $R^*_{abcd} = frac12 R_{abef} epsilon^{ef}{}_{cd}$ the right Hodge dual of the Riemann curvature tensor. Taking the Hodge dual w.r.t. to the antisymmetric $a,b,c$ indices in $Upsilon(gamma)$ we can study the equivalent linear map $$ Upsilon^*: gamma mapsto gamma_{d[a]}{}^f R^*_{b]fc}{}^d $$ The image is a $(0,3)$ tensor that is antisymmetric in its first two indies, so $Upsilon^*$ is a mapping from a 40 dimensional space to a 24 dimensional one, and still must have a non-trivial kernel in general.

(Remark: in the two dimensional case, however, using that Riemann curvature has special structure, you see that curvature collineation implies that the vector field $X$ is conformal Killing, so there's still some hope of getting the desired inequality through a different route.)

---

In conclusion, I doubt that there is a general statement that curvature collineation implies Killing identity, because there are too many special cases to worry about. At best (for this line of argument) you may have a statement that states that for generic metrics in sufficiently high dimension that curvature collineation implies Killing identity.