= Solution
Use the paper's curvature convention throughout, and define the <Ricci tensor> and <scalar curvature> by $R_{kn}=R^i{}_{kin}$ and $R=g^{kn}R_{kn}$. In particular, if semicolons are read from left to right, $X^i{}_{;km}=\nabla_m(\nabla_kX)^i$, including the connection on the derivative index. This ordering matters for the <Ricci identity> sign.
Choose <Riemann normal coordinates> at an arbitrary point $P$. The <Christoffel symbols> vanish there, and the <Levi-Civita connection> formula gives
$$
R_{ikmn}(P)=\frac12\left(g_{im,kn}+g_{kn,im}-g_{km,in}-g_{in,km}\right)(P).
$$
Metric symmetry and commuting partial derivatives show directly that exchanging $m,n$ or $i,k$ changes the sign, while interchanging the pairs $(i,k)$ and $(m,n)$ leaves the expression unchanged. Adding its three cyclic versions over $k,m,n$ cancels every second derivative. These are tensorial statements and $P$ was arbitrary, so
$$
\boxed{R^i{}_{k(mn)}=0,\qquad R^i{}_{[kmn]}=0,\qquad R_{(ik)mn}=0,\qquad R_{ikmn}=R_{mnik}.}
$$
The cyclic relation is the <first Bianchi identity>; its reduction to the total antisymmetrization uses last-pair antisymmetry. These arguments use zero <torsion tensor> and <metric compatibility>, not a field equation.
At $P$, the <covariant derivative> of curvature equals its partial derivative. Differentiating the connection expression for curvature gives derivatives of $\partial_n\Gamma^i{}_{km}-\partial_m\Gamma^i{}_{kn}$; derivatives of the quadratic connection terms vanish because $\Gamma(P)=0$. In the cyclic sum over $m,n,p$, every second derivative of a connection coefficient occurs twice with opposite signs. Thus
$$
R^i{}_{kmn;p}+R^i{}_{knp;m}+R^i{}_{kpm;n}=0,
\qquad\boxed{R^i{}_{k[mn;p]}=0.}
$$
This is the <second Bianchi identity>, valid everywhere by tensoriality.
Lower the first curvature index and contract the differential identity with $g^{im}$. <Metric compatibility> allows the metric to pass through the <covariant derivatives>. Using the pair antisymmetries yields
$$
\boxed{\nabla^iR_{ikmn}=\nabla_mR_{kn}-\nabla_nR_{km}.}
$$
Contract again with $g^{km}$. The left side is $-\nabla^iR_{in}$, while the right side is $\nabla^kR_{kn}-\nabla_nR$. Hence the <contracted Bianchi identity> is
$$
\boxed{\nabla^iR_{in}=\frac12\nabla_nR,\qquad
\nabla^i\left(R_{in}-\frac12g_{in}R\right)=0.}
$$
The <Ricci tensor> is symmetric, as follows by contracting the pair-exchange symmetry. The final divergence is that of the <Einstein tensor>.
For the vector commutator, at the same normal-coordinate point,
$$
X^i{}_{;km}=\partial_m\partial_kX^i+(\partial_m\Gamma^i{}_{jk})X^j.
$$
The terms involving undifferentiated connection coefficients vanish at $P$. Subtracting the reversed expression leaves $(\Gamma^i{}_{jk,m}-\Gamma^i{}_{jm,k})X^j$, exactly the stated curvature convention at $P$. Therefore
$$
\boxed{X^i{}_{;km}-X^i{}_{;mk}=R^i{}_{jkm}X^j.}
$$
For a <covector field>, its connection term has the opposite sign: $\omega_{i;k}=\partial_k\omega_i-\Gamma^j{}_{ik}\omega_j$. The same calculation gives
$$
\boxed{\omega_{i;km}-\omega_{i;mk}=-R^j{}_{ikm}\omega_j.}
$$
Equivalently, apply the commutator to the scalar $\omega_iX^i$, whose two <covariant derivatives> commute, and cancel the vector contribution. Each covariant index contributes a negative curvature action and each contravariant index a positive one in this semicolon ordering.
To prove the requested <double divergence of the Riemann tensor>, let
$$
D_{mn}=R^{ik}{}_{mn;ik}=\nabla_k\nabla_iR^{ik}{}_{mn}.
$$
The once-contracted identity gives
$$
D_{mn}=\nabla_k\nabla_mR^k{}_n-\nabla_k\nabla_nR^k{}_m.
$$
Commuting the outer derivative past $\nabla_m$ and $\nabla_n$ leaves a difference of scalar Hessians, $\tfrac12(\nabla_m\nabla_n-\nabla_n\nabla_m)R=0$, plus curvature terms. The operator commutator $[\nabla_k,\nabla_m]$ has the negative of the paper's curvature sign, since it reverses the semicolon order. On the mixed <Ricci tensor> it gives
$$
[\nabla_k,\nabla_m]R^k{}_n
=-R_{am}R^a{}_n+R^a{}_{nkm}R^k{}_a.
$$
The first term is symmetric in $m,n$: it is the metric contraction of two symmetric <Ricci tensors>. The second is $S_{mn}=R^{ak}R_{ankm}$ and is also symmetric, since
$$
S_{mn}=R^{ak}R_{kman}=R^{ak}R_{amkn}=S_{nm},
$$
where the first equality uses curvature pair exchange and the second swaps the dummy indices $a,k$. Thus the two commutator terms cancel on antisymmetrization in $m,n$, proving
$$
\boxed{R^{ik}{}_{mn;ik}=0.}
$$
No assumption such as vacuum, constant curvature or vanishing <Ricci tensor> was used.
Back to article page