loislane said:
The distinction I was referring to was between Riemannian and pseudo-riemannian geometry.
I don't see any real distinction between the two with regard to the equations I wrote down; both of them are valid regardless of the metric signature.
loislane said:
The point being made in those posts is that, if we look at the Einstein-Hilbert action for gravity alone, it can only be linear in second derivatives of the metric (i.e., in ##R##, the Ricci scalar); otherwise we would get a third or higher order PDE for the field equation for gravity.
But we haven't been looking at the field equation for gravity; we've been looking at the field equation for electromagnetism, which is derived by varying the Lagrangian with respect to the vector potential ##A_{\mu}##, not the metric. Terms of higher order in ##R## aren't functions of the vector potential, so, from the standpoint of Maxwell's Equations, if we're going to include any curvature terms, we should include them all, regardless of order.
We could, however, look at how adding those terms to the Lagrangian would affect the Einstein Field Equation; what we will find is that even the ##\alpha## term would add extra unwanted terms to that equation. The usual EFE is derived from the action (leaving out irrelevant factors in the terms and including the EM Lagrangian--we assume there is no other matter present):
$$
S = \int d^4 x \sqrt{-g} \left( R - \frac{1}{4} F_{\mu \nu} F^{\mu \nu} \right)
$$
Varying this with respect to the metric gives the usual EFE with the EM stress-energy tensor on the RHS (which I won't bother writing out explicitly since its form won't matter here):
$$
R_{\mu \nu} - \frac{1}{2} g_{\mu \nu} R = T_{\mu \nu}
$$
But now suppose we include the extra ##\alpha## term in the Lagrangian, so that varying with respect to the vector potential puts the extra ##\alpha## term in Maxwell's Equations:
$$
S = \int d^4 x \sqrt{-g} \left( R - \frac{1}{4} F_{\mu \nu} F^{\mu \nu} - \frac{1}{4} F_{\mu \nu} F^{\mu \nu} \alpha R \right)
$$
Now, when we vary this with respect to the metric, we get:
$$
\left( R_{\mu \nu} - \frac{1}{2} g_{\mu \nu} R \right) \left( 1 - \alpha \frac{1}{4} F_{\mu \nu} F^{\mu \nu} \right) = T_{\mu \nu} \left( 1 + \alpha R \right)
$$
This mess is clearly not a clean field equation for the metric (i.e., an equation where the LHS has no matter terms and the RHS has no curvature terms. (And if we try to make it into one by algebraically moving terms and factors around, what we end up with will not be a clean second-order PDE anyway, since it will have to divide the Einstein tensor by the factor ##\left( 1 + \alpha R \right)## on the LHS.)
In other words, you can't add any curvature terms to the field equations for matter (such as Maxwell's Equations) without also adding unwanted terms to the field equation for gravity (as above). That's a powerful argument in favor of the minimal coupling principle.