The problem: I am trying to prove that, for a Riemannian manifold $(M,\langle\ , \ \rangle)$, $X \in \Gamma(TM)$ is a Killing field (i.e. one for which $\langle u,v \rangle_{p \in M} = \langle (d \phi^X_{t})_{\phi^{X}(t, p)} (u), (d \phi^X_{t})_{\phi^{X}(t, p)} (v) \rangle_{\phi^{X}(t, p)}$, where $\phi^{X}_{t} = \phi^{X}(t,) $ is the local flow of $X$) if $\forall Y,Z \in \Gamma(TM)$ $$ \langle \nabla_Y X, Z \rangle + \langle \nabla_Z X, Y \rangle = 0, $$ where $\nabla$ is the Levi-Civita connection.
My work so far:
By using the compatibility and torsion-freeness of the connection, it is easy to reduce this assumption to $$ X \langle Y,Z \rangle = \langle [X,Z],Y \rangle + \langle [X,Y],Z \rangle. $$ Then, if we choose $Y,Z$ such that $Y(p) = u$, $Z(p) = v$, and remembering that $X(p) = \frac{d}{dt}\Big|_{t=0} \phi^{X}(t,p)$, our LHS evaluated at $p$ becomes $$ \frac{d}{dt} \Big|_{t=0} \langle Y(\phi^{X}(t, p)), Z(\phi^{X}(t, p)) \rangle_{\phi^{X}(t, p)}. $$ Our RHS at $p$ becomes $$ \langle [X,Z](p),Y(p) \rangle_p + \langle [X,Y](p),Z(p) \rangle_p . $$ At this point I start looking at the definition of the Lie bracket: $$ [X,Y](p) = \frac{d}{dt}\Big|_{t=0} d \phi^X_{-t}(Y(\phi^X(t,p))), $$ so our RHS becomes $$ \langle \frac{d}{dt}\Big|_{t=0} d \phi^X_{-t}(Z(\phi^X(t,p))), Y(p) \rangle_p + \langle \frac{d}{dt}\Big|_{t=0} d \phi^X_{-t}(Y(\phi^X(t,p))) , Z(p) \rangle_p,$$ but I am having a hard time completing the proof. My hunch is that we want to show the RHS at $p$ is zero (I could be wrong here). Any hints or pointers would be very welcome. Apologies if my notation is a bit untidy. There is also probably a much easier way to do this, so any hints away from my solution are also welcome.