He estado trabajando durante algún tiempo en la derivación de las ecuaciones de movimiento (MOE) para el de Gauss-Bonnet-de la Gravedad, la cual está dada por la acción:
$$\int d^D x \sqrt{|g|} (R^2-4R_{ab}R^{ab}+R_{abcd}R^{abcd}).$$
He intentado durante algún tiempo para obtener el primer lugar-fin de variación con respecto al $\delta g_{ab}$. Siempre he sido atrapado con la derivación de la $\delta (R_{ab}R^{ab})$ componente. Supongamos que tengo sólo un tensor de Riemann que contiene sólo las constantes, \begin{align*} \delta(R_{ab}R^{ab}) & = (\delta R_{ab})R^{ab} + (\delta R^{ab})R_{ab} \\ & = (\delta R_{ab})R^{ab}+ \delta(R_{ecfd}g^{ac}g^{bd})g^{ef}R_{ab} \\ & = (\delta R_{ab})R^{ab}+ (\delta R_{ecfd})g^{ac}g^{bd}g^{ef}R_{ab} + (\delta g^{ac}) R_{cd}g^{bd}R_{ab}+(\delta g^{bd})R^{a}{}_{d}R_{ab}\\&\quad+ (\delta g^{ef}) R_{e}{}^{a}{}_{f}{}^{b}R_{ab}. \end{align*}
Esto, sin embargo, nos da: \begin{align*} & = - (\delta g^{hi}) R_{ichd}R^{cd} +( \delta g^{hi}) R_{hcid}R^{cd} + 2(\delta g^{ac}) R_{cd}R_{a}{}^{d} \\ & = 2(\delta g^{ac}) R_{cd}R_{a}{}^{d}. \end{align*}
La respuesta correcta tal y como la conocemos, es $2R^{ab}R_{acbd}\delta g^{cd}$.