As Sabine Hossenfelder has pointed out, MOND appears to be a more successful explanation for the phenomena called dark matter than the Einstein field equations with axions. She has a rather complicated theory called Covariant Extended Gravity. Instead, I would like to try a merely Lorentz-covariant hypothesis which uses just the gravitational force field, Christoffel symbols of the second kind for the Levi-Civita connection, and the Lagrange multiplier λμ for a coordinate condition instead of what she calls the "imposter field".
which was chosen because it is Lorentz-covariant and it ignores the scale of the metric. Consider the covariant derivative
The last term emerges because is not an invariant scalar, and so its covariant derivative is not the same as its ordinary derivative. Rather, .
Contracting ν with ρ and applying the coordinate condition, we get:
Thus, we get that an alternative way of expressing the coordinate condition is:
.
.
Note that this coordinate condition must be satisfied at all times, not just the present, but also the past and the future. Only then, is the theory consistent with special relativity. And only then can one expect it to yield MOND as a possible result.
The equations of motion will include our coordinate condition (from varying the Lagrange multiplier) and a modified form of the Einstein field equations which will incorporate the result of varying the new terms with respect to the inverse metric tensor. Variation of the new terms gives
.
where parentheses indicate symmetrization of the indices. So each of the ten Einstein field equations will have these four terms added to it:
If we turn this around, we can express the partial derivative in terms of the covariant derivative
which will help us to convert more of our equation of motion into an explicitly invariant form. Using our coordinate condition, we get
and thus
.
also
.
so the correction to the Einstein field equations becomes
where the first two terms are the stress-energy tensor for an invariant symmetric divergence of λμ (with zero trace and thus zero mass, similar to electromagnetism) and the last two terms are a kind of source term for λμ which is required to comply with our coordinate condition.
Assuming that is a scalar field (so its gradient is a covariant vector) and observing that when in one of our preferred reference frames
we can infer that our entire correction to each of the Einstein field equations is equal to
which is an invariant since it is built entirely from tensors. Thus we have rendered our theory into a generally invariant form so that we can use spherical coordinates or cylindrical coordinates or whatever coordinate system we may need.
However, to properly interpret this correction, we must remember that
This would be the case, if one is considering a static spherically-symmetric black hole or collapsar in an asymptotically Minkowskian space-time where:
there is no ordinary matter or radiation outside rmin;