I am learning the basics of reformulating vector field equations in terms of differential forms. In particular, I am trying to express the equations for the Riemannian connection using differential forms. The Riemannian connection equations for an affine connection $\nabla$ on a Riemannian manifold $(M,g)$ are
$$X (g(Y,Z)) - g(\nabla_X Y, Z) - g(Y,\nabla_XZ) = 0, \ \forall X,Y,Z \in \Gamma(TM),$$ $$ \nabla_XY - \nabla_YX - [X,Y] = 0, \ \forall X,Y \in \Gamma(TM). $$
However, I'm finding it difficult to convert each to one-form equations and would appreciate any pointers or references on this.
Work (resolved)
We work in an orthonormal framing $\{ E_i \}$ on the tangent bundle, with dual basis $\{ \omega^i \}$. We define the one-form $\omega^j_i$ by $\omega^j_i (X) := \langle \nabla_X E_i, E_j \rangle$. Then, working on the first equation, and letting $\langle X,Y \rangle := g(X,Y),$ we get
$$\begin{align} \langle \nabla_XY, Z \rangle &= \ldots = Z^j(X^iE_i(Y^j) + Y^i\omega^{j}_{i}(X)), \\ \langle Y,\nabla_X Z \rangle &= \ldots = Y^j(X^iE_i(Z^j) + Z^i \omega^j_i (X)), \\ X \langle Y,Z \rangle &= \ldots = X^i Y^j E_i (Z^j) + X^i Z^j E_i(Y^j), \\ \implies 0 &= (Z^jY^i\omega^j_i + Y^jZ^i\omega^j_i)(X), \ \forall X,Y,Z \in \Gamma(TM) \\ \implies 0 &= (Z^jY^i)(\omega^j_i + \omega^i_j)(X), \ \forall X,Y,Z \in \Gamma(TM) \\ \implies 0 &= \omega^j_i + \omega^i_j. \end{align}$$
This part is fine. Now, on setting $[E_i,E_j]=c_{ij}^k E_k$ and slotting in the coordinates to the torsion condition, I get $$ \Gamma^k_{ij} - \Gamma^k_{ji} - c^k_{ij} =0, $$ where $\Gamma^i_{jk}E_i:=\nabla_{E_j}E_k$. Now clearly $\Gamma^{k}_{ij} = \omega^k_j(E_i)$ and $-c^{k}_{ij}=d\omega^k(E_i,E_j)$, so we have $$ \begin{align} \omega^k_j(E_i) - \omega^k_i (E_j) + d\omega^k(E_i,E_j) &=0, \ \forall i,j. \end{align}$$ However, $$ (\omega^p \wedge \omega^{k}_p)(E_i,E_j) = \begin{vmatrix} \omega^p(E_i) & \omega^p(E_j) \\ \omega^k_p(E_i) & \omega^k_p(E_j) \end{vmatrix} = \ldots = \omega^k_i(E_j) - \omega^k_j(E_i).$$ Putting this all together, we get $$ \begin{align} d\omega^k(E_i,E_j) &= \omega^p \wedge \omega^k_p (E_i,E_j), \ \forall i,j \\ \implies d\omega^k &= \omega^p \wedge \omega^k_p. \end{align} $$
So our corresponding differential system is
$$ \begin{align} \omega^j_i + \omega^i_j &= 0, \\ d\omega^i - \omega^j \wedge \omega^i_j &= 0.\end{align} $$