I'm asked to prove the following statement in my physics book:
A vector field with covariant components $v^b$, in order to have a vanishing covariant derivative everywhere in a manifold, must satisfy: $$(\partial_b\Gamma^{d}{}_{ac}-\partial_c\Gamma^{d}{}_{ab}+\Gamma^{e}{}_{ac}\Gamma^{d}{}_{eb}-\Gamma^{e}{}_{ab}\Gamma^{d}{}_{ec})v^a=0.$$
Edit: This is what I tried after @PraharMitra's suggestion:
Since $\nabla_a v^b=0$, clearly $[\nabla_b,\nabla_c]v^d=\nabla_b\nabla_cv^d-\nabla_c\nabla_dv^d=0$.
As the covariant derivative of a contravariant component is defined as $\nabla_bv^d=\partial_bv^d+\Gamma^d{}_{eb}v^e$ and the covariant derivative of a covariant component as $\nabla_bv_d=\partial_bv_d+\Gamma^e{}_{db}v_e$, I got to:
$$ \nabla_b(\nabla_cv^d)=\partial_b(\nabla_cv^d)-\Gamma^e{}_{cb}(\nabla_ev^d)+\Gamma^d{}_{eb}(\nabla_cv^e) $$
Now let's plug in the covariant derivatives with respect to c:
$$ \nabla_b(\nabla_cv^d)=\partial_b(\partial_cv^d+\Gamma^d{}_{ac}v^a)-\Gamma^e{}_{cb}(\partial_ev^d+\Gamma^d{}_{ae}v^a)+\Gamma^d{}_{eb}(\partial_cv^e+\Gamma^e{}_{ac}v^a) $$ $$=\partial_b\partial_cv^d+\partial_b\Gamma^d{}_{ac}v^a-\Gamma^e{}_{cb}\partial_ev^d+\Gamma^e{}_{cb}\Gamma^d{}_{ae}v^a+\Gamma^d{}_{eb}\partial_cv^e+\Gamma^d{}_{eb}\Gamma^e{}_{ac}v^a$$
From this expression it's straightforward to obtain the $\nabla_c(\nabla_bv^d)$ term by exchanging the b and c indexes. I got:
$$ \nabla_c(\nabla_bv^d)=\partial_c\partial_bv^d+\partial_c\Gamma^d{}_{ab}v^a-\Gamma^e{}_{bc}\partial_ev^d+\Gamma^e{}_{bc}\Gamma^d{}_{ae}v^a+\Gamma^d{}_{ec}\partial_bv^e+\Gamma^d{}_{ec}\Gamma^e{}_{ab}v^a $$
Now, on substracting:
$$ \nabla_b(\nabla_cv^d)-\nabla_c(\nabla_bv^d)=\partial_b\partial_cv^d+\partial_b\Gamma^d{}_{ac}v^a-\Gamma^e{}_{cb}\partial_ev^d+\Gamma^e{}_{cb}\Gamma^d{}_{ae}v^a+\Gamma^d{}_{eb}\partial_cv^e+\Gamma^d{}_{eb}\Gamma^e{}_{ac}v^a-[\partial_c\partial_bv^d+\partial_c\Gamma^d{}_{ab}v^a-\Gamma^e{}_{bc}\partial_ev^d+\Gamma^e{}_{bc}\Gamma^d{}_{ae}v^a+\Gamma^d{}_{ec}\partial_bv^e+\Gamma^d{}_{ec}\Gamma^e{}_{ab}v^a] $$
The terms with both the partial derivatives vanish, since those can be exchanged:
$$=\partial_b\Gamma^d{}_{ac}v^a-\Gamma^e{}_{cb}\partial_ev^d+\Gamma^e{}_{cb}\Gamma^d{}_{ae}v^a+\Gamma^d{}_{eb}\partial_cv^e+\Gamma^d{}_{eb}\Gamma^e{}_{ac}v^a-[\partial_c\Gamma^d{}_{ab}v^a-\Gamma^e{}_{bc}\partial_ev^d+\Gamma^e{}_{bc}\Gamma^d{}_{ae}v^a+\Gamma^d{}_{ec}\partial_bv^e+\Gamma^d{}_{ec}\Gamma^e{}_{ab}v^a] $$
In this last expression, the first, sixth, eighth and tenth terms can be put together to be $(\partial_b\Gamma^{d}{}_{ac}-\partial_c\Gamma^{d}{}_{ab}+\Gamma^{e}{}_{ac}\Gamma^{d}{}_{eb}-\Gamma^{e}{}_{ab}\Gamma^{d}{}_{ec})v^a$, which is recognised to be the term we wanted to prove to be 0 in order to have vanishing covariant derivative (note that the terms with two connection coefficients commute). We can tell further that $(\partial_b\Gamma^{d}{}_{ac}-\partial_c\Gamma^{d}{}_{ab}+\Gamma^{e}{}_{ac}\Gamma^{d}{}_{eb}-\Gamma^{e}{}_{ab}\Gamma^{d}{}_{ec})v^a=R^d{}_{abc}v^a$, the riemannian curvature. On rewriting, it remains:
$$\nabla_b(\nabla_cv^d)-\nabla_c(\nabla_bv^d)=R^d{}_{abc}v^a+[-\Gamma^e{}_{cb}\partial_ev^d+\Gamma^e{}_{cb}\Gamma^d{}_{ae}v^a+\Gamma^d{}_{eb}\partial_cv^e]-[-\Gamma^e{}_{bc}\partial_ev^d+\Gamma^e{}_{bc}\Gamma^d{}_{ae}v^a+\Gamma^d{}_{ec}\partial_bv^e] $$
Assuming the manifold is torsionless, we would have $\Gamma^a{}_{bc}=\Gamma^a{}_{cb}$, and so both the first and second therm in each bracket would cancel, leaving us with:
$$\nabla_b(\nabla_cv^d)-\nabla_c(\nabla_bv^d)=0=R^d{}_{abc}v^a+\Gamma^d{}_{eb}\partial_cv^e-\Gamma^d{}_{ec}\partial_bv^e$$
In order to prove our initial statement, it seems like those two terms with the connection coefficients must vanish, but I don't see how do they or where did I make a mistake...