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".
where: LO is the Lagrangian in classical physics of ordinary matter and radiation, excluding gravity and dark matter. (Here, LO includes quantum effects. That is, this is classical in the sense of excluding gravity, not quantum.) κG and κD are proportionality constants to match the units of gravity and dark matter to those of classical physics.
If dark matter is described as a merely Lorentz-covariant field, then LD might take the form
or
where s is a scale factor, the classical Newtonian gravitational field. M and λ are Lagrange multipliers. (See also Nordström's theory of gravitation and the Weyl curvature tensor.)
Really special relativity
Let us now assume that there is a non-empty set of preferred coordinate systems in which
holds everywhere and always where ηαβ is the Minkowski metric (not a tensor, but just a mathematical constant).
Moving towards general invariance
We need to transform to an arbitrary coordinate system rather than confining ourselves to the preferred coordinates. This is necessary because: we may not initially know which are the preferred coordinates; we may want to test whether there are in fact such preferred coordinates; and the symmetries of the systems we want to study may require the use of different coordinate systems. The main difficulty with such a change is caused by the need to calculate the covariant derivatives with respect to the new coordinates. That, in turn, requires knowledge of
for arbitrary coordinates. To obtain that, we can first calculate it for our preferred coordinates and then transform it to arbitrary coordinates.
where δεα is the Kronecker delta. Notice that the last line is meaningful for arbitrary coordinates.
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.
If we turn this around, we can express the partial derivative in terms of the covariant derivative
Simplification
In one of our preferred coordinate systems, the metric tensor is simply the product of the Minkowski metric and a scale factor
.
From this we get that the gravitational force field (in such a preferred coordinate system) is
.
So the non-invariant terms in the correction to the Einstein field equations become
.
Establishing general invariance
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
so that
when the metric is diagonal
Black hole
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;