Relativistic MOND from modified energetics

We begin to investigate the question of what modifications in the energy-momentum tensor can yield the correct MOND regime. As a starting study, we refrain from insisting on an action principle and focus exclusively on the equations of motion. The present work, despite the absence of an explicit action functional, can be regarded to extend Milgrom’s modified inertia approach to relativistic domain. Our results show that a proper MOND limit arises if the energy-momentum tensor is modified to involve the determinant of the metric tensor in reference to the flat metric, where the latter is dynamically generated as in the gravitational Higgs mechanism. This modified energy-momentum tensor is conserved in both Newtonian and MONDian regimes.


Introduction
Observations of several decades, ranging from the initial measurements by Oort (see the discussion in [1]) to the primal ones by Rubin [2,3], have shown that galaxies exhibit flat rotation curves, manifestly violating the Keplerian dynamics. This universal anomalous dynamics has been interpreted in two distinct ways. The first, first proposed by Zwicky [4] in 1933, refers to the dark matter (DM) hypothesis. According to the DM paradigm, there must be a distribution of non-shining matter at the outer skirts of galaxies to yield the approximately constant velocities after particular distances from the center of galaxies. The DM hypothesis provides viable explanations not only for flat rotation curves but also for various cosmological and astronomical observations describing different phases of the evolution of the Universe. Several experimental groups have been searching for DM particles by utilizing various detection methods (see the recent review volumes [5,6]). So far, no signal of DM has been observed. a e-mail: demir@physics.iztech.edu.tr b e-mail: cananduzturk@iyte.edu.tr The second interpretation, first proposed by Milgrom [7][8][9] in 1983, postulates that the observed flat rotation curves result from modifications in the Newtonian laws of motion. In this approach, instead of adding unknown ingredients to the galactic matter, one exercises modifications in the equations of motion which dominate at the skirts of the galaxies. To this end, Newton's law of motion F = ma changes to where F is the net force acting on the material point which has inertia m and acceleration a (with |a| 2 = a · a). This dynamical equation, structuring Milgrom's MOND theory [7][8][9], is characterized by the empirical function μ(|a|/a 0 ) where a 0 1.2 ×10 −10 ms −2 is a constant acceleration scale for all galaxies [10]. It appears in (1) as a critical acceleration scale set galactically by the mass M and radius R of the galaxy as (G N M)/R 2 a 0 and cosmologically by the present-day value H 0 of the Hubble parameter as (cH 0 )/2π a 0 [11]. The heart of the MOND theory is the empirical function μ(|a|/a 0 ). There is yet no dynamical theory for it; however, its asymptotic behavior is not difficult to guess: if all the successes of the Newtonian theory are to be maintained. Here x does not need to be very large or small compared to unity because μ(x) can attain its asymptotics even when x is close to unity. For instance, the empirical form facilitates the asymptotics in (2) almost independently of x provided that n is large. Indeed, taking n = 50 one finds μ(x) = 0.700, 0.900, 0.986, 0.996, 1.000, 1.000 for x = 0.7, 0.9, 0.99, 1.01, 1.1, 2.0, respectively.
The behavior in (2) ensures that matter in the galaxy exhibits flat rotation curves far away from the galactic center. Indeed, in the limit of small accelerations the equation of motion (1) takes the form so that at large radii R corresponding to the outer skirts of the galaxy one does not find the Keplerian law |F| = (mv 2 )/R but |F| = (mv 4 )/(a 0 R 2 ), which yields the constant speed for |F| = (G N m M)/R 2 . This relation accounts for the observed flat rotation curves [2,3,[12][13][14]. The constant speed (5) is the reason for and result from the whole idea of MOND. It depends crucially on the behavior of the empirical function (2) at low accelerations. The empirical MOND relation in (1), supported by (2) and (3), needs be formulated at a more fundamental level. In this regard, there arise two different interpretations. In the first, after setting a = −∇φ g with φ g being the gravitational potential, one formulates MOND as a modification in gravitational laws (see the reviews [15][16][17]). In this case, one is necessarily led to modified Newtonian gravity [18,19] or General Relativity (GR) extended by geometrical scalar and vector fields [20][21][22][23][24][25]. Besides, there are alternative approaches based on f (R) gravity [26][27][28][29], bimetric gravity [30,31], time foliation [32,33], nonlocal metric theories [34,35], Galileons [36], and Horava-Lifshitz gravity [37]. In general, modified gravity theories introduced to replace the DM necessarily lead to a MONDian structure.
In the second interpretation, one conceives of the equation of motion (1) as defining an acceleration-dependent inertia m(|a|) = mμ(|a|/a 0 ). This approach, the modified inertia approach proposed in [38,39], in the non-relativistic limit, keeps gravitational laws unchanged yet lets in nonlinear kinetic terms. In this framework, it is found that the kinetic term of the point mass involves all derivatives of acceleration [38][39][40][41], yet it is stable and respects causality [15][16][17]34,35]. In the present work, we pursue this modified inertia viewpoint to generalize it to general-relativistic domain. The experience from non-relativistic study [38,39] ensures that forming an action functional must be difficult, if not impossible, in the relativistic domain. We thus focus exclusively on the equations of motion without specifying an action principle to derive them.

Modified energetics
As the beginning phase of a study program aiming at finding dynamical alternatives to modified gravity models of relativistic MOND [15][16][17], in this section we study grav-itational field equations where the MOND phase is understood as changes in the matter energy-momentum tensor. This approach, aiming at carrying Milgrom's modified inertia approach [38,39] into the relativistic domain at the level of equations of motion, is based on the matter energy-momentum tensor T (N ) μν in the Newtonian domain and exploits its expected non-conservation in the MOND regime to derive the MONDian dynamics in an empirical way. Having a complete knowledge of the interactions of matter, its energy-momentum tensor T (N ) μν (with energy density T N 00 , pressure T N ii , momentum density T N 0i and shear stress T N i j ) are strictly conserved in the Newtonian regime. However, the same T μν to a conserved energy-momentum tensor T μν which can be approached only empirically in the absence of a complete dynamical model (see [43,44] for a similar approach to a modified gravity framework for MOND). We now give an empirical implementation of this dynamical picture starting with the Einstein field equations, in which T μν is the conserved energy-momentum tensor of matter at all acceleration scales ranging from |a| = 0 to |a| = ∞. In general, T μν is conserved on the equations of motion, and these equations necessarily encode the novel interactions of matter responsible for the MOND. However, those new interactions are not known and our knowledge of T μν is incomplete; we are able to know it only when |a| > a 0 , for which it equals T (N ) μν . Consequently, on an empirical basis we write for T μν where μ(x) is the MOND function in (2), Q is a scalar, and a is yet another scalar which is to be judiciously constructed to have the empirical limit at non-relativistic energies. This correspondence between the relativistic (a) and non-relativistic (|a|) regimes is crucial for the empirical structure in (7) to give a consistent framework. Physically, the grand energy-momentum tensor T μν must correctly reproduce the Newtonian and MONDian regimes. This is analyzed case by case in Table 1 as a function of the divergence of T (N ) μν . As suggested by the table, the underlying dynamics can be revealed after a proper understanding of T μν and this requires T (N ) μν , a(T ) and Q(T ) to be constructed in detail. We detail these physical variables in the three consecutive subsections that follow. Table 1 The acceleration dependence of the energy-momentum tensor T μν of matter. In general, a = a(T (N ) ) and Q = Q(T (N ) ) are functions of the energy-momentum tensor T (N ) μν . These scalars take appropriate values for Newtonian (T (N ) μν is conserved) and MONDian (T (N ) μν is not conserved) regimes. Namely, matter develops novel interactions (such as the higher-derivative kinetic terms, determined in [38,39] in the non-relativistic regime) at small accelerations and its known energymomentum tensor T (N ) μν starts exhibiting non-conservation properties Acceleration MOND function Energy-momentum tensor Matter dynamics This is the 'Newtonian regime'. The acceleration of matter is above a 0 and μ (a) ensures μν is symmetric and divergence-free (∇ μ T (N ) μν = 0) in agreement with (6). In the Newtonian regime thus T (N ) μν qualifies as the well-known conserved energy-momentum tensor of matter This is the 'MONDian regime'. The acceleration of matter is below a 0 and μ (a) leads to μν , which qualifies as the conserved energy-momentum tensor of matter. In this small acceleration regime, matter develops novel interactions that make ∇ μ T (N ) μν = 0 yet the scalars a and Q help T μν be conserved and give the observed flat rotation curves

Physical properties of T (N ) μν
It has been emphasized previously, specifically in Table 1, that T (N ) μν has the same form as the energy-momentum tensor of matter in the Newtonian regime yet it does not qualify as a true energy-momentum tensor in the MOND regime simply because its conservation is spoiled by novel interactions of matter that arise at accelerations below a 0 . The higher-derivative self interactions studied in [38][39][40][41] form a concrete example of such effects. Let us consider, as an illustrative example, dust (pressureless matter having only energy density in the comoving frame) for which where ρ and u μ are the energy density and velocity, respectively. (One recalls that T (N ) μν = dτρu μ u ν for a relativistic particle with trajectory y μ (τ ) and energy density ρ = mc 2 δ 4 (x − y(τ ))). It is divergence-free, ∇ μ T (N ) μν = 0, because densities and flows of dust are all conserved. However, this conservation property holds only in normal circumstances where the Newtonian laws of motion are valid. In the MONDian regime, where dust develops higher-derivative kinetic interactions for instance, conservation breaks down, ∇ μ T (N ) μν = 0. On dimensional grounds, it is likely to have a structure of the form in addition to terms involving derivatives of the acceleration.
In the absence of an invariant action (like the non-relativistic model in [38,39]), this non-conservation can be understood neither in origin nor in structure (ρa 0 u ν in (10) is just an example). Therefore, our goal is not to construct a model of the non-conservation of T (N ) μν , but to determine its consequences for structures and dynamics of a and Q.

Physical Properties of the Acceleration Scalar a
The acceleration scalar a, which must have the nonrelativistic limit a NR given in (8), must be constructed judiciously to correctly cover the Newtonian and MONDian regimes. Hence, besides the crucial relation (8), it must have the following properties. Table 1, a must vary with the divergence of T

By our construction shown in
while ∇ μ T μν = 0 in both cases. 2. Being a scalar field, a involves contractions of the divergences of T (N ) μν . This necessarily brings in the gravitational acceleration ∇φ g through the gravitational potential φ g = −1 − g 00 arising in the Newtonian limit of the metric tensor g μν . However, the presence of ∇φ g must be prohibited for a to yield the kinetic acceleration in (8). It is easy to see that this cannot be accomplished without using an independent source of φ g and the most natural source as such is the determinant g = Det(g μν ) of the metric tensor. However, being a scalar density rather than a scalar, g cannot appear in a by itself; it must be divided by another scalar density to achieve covariance. This other scalar density necessitates a new metric g μν , and naturally leads to a bimetrical picture (whose relevance for MOND has been discussed in [30,31]). Then the acceleration scalar possess the functional form a = a a 0 , ∇ μ T (N ) μν , T (N ) μν , g μν , g μν g μν , g/g (12) where g = Det(g μν ) arises as an additional variable to be dynamically determined.
These two points plus (8) must be taken into account in formulating a. However, the formulation process becomes utterly incomplete unless the additional metric g μν is demystified. In the two subsections that follow, we first study g μν and then construct a model of a.

Construction of g μν
The second metric tensor g μν , required to eliminate the gravitational acceleration ∇φ g from the acceleration scalar a, can be ascribed different structures depending on the underlying dynamics. For instance, one may consider identifying it with T (N ) μν itself but this attempt fails because its determinant vanishes in the case of dust (see Eq. 9 above). Alternatively, one may take g μν as a second metric tensor with its own curvature and dynamics, but this setup, as was already elaborated by Milgrom in [30,31] (see also [42]), gives a modified gravity theory for MOND. This and other possible modified gravity models fall outside the scope of the present work because the goal here is to develop a dynamical approach to relativistic MOND similar in philosophy to Milgrom's modified inertia approach [38,39].
Our approach to g μν is dynamical rather than geometrical. In other words, the dynamics underlying the asymptotics in Table 1 and structures in (12) proceed not only with T (N ) μν but also g μν . Thus, g μν is a low-acceleration dynamical field, maybe one of many as such, which facilitates the MOND regime. In modeling the dynamics, we interpret the coupling g μν g μν between the two metrics as the kinetic term of four real scalars φ m (m = 0, . . . , 3), and we construct the defining relation where η mn is the flat Minkowski metric, and hence the scalar spectrum contains a ghosty (negative kinetic term) mode. We assume that φ m develops the nontrivial backgrounds depending on whether the diffeomorphism invariance is exact ( φ m = 0) or spontaneously broken ( φ m = M 2 x a ) in the vacuum state governed by the vacuum expectation value φ a of the scalars. Here, the scale M is around a 0 . The dynamics leading to (14) can be known only in a setting where all interactions of matter and extra fields like φ a are specified. The diffeomorphism-breaking vacuum here sets the flat Minkowski metric η μν as the background metric about which g μν can be expanded in a perturbation series. This induction mechanism is similar to what happens in the gravitational Higgs mechanism [46][47][48][49][50], in which a second metric tensor g μν is needed for writing a sensible graviton mass term through the kinetic term g μν g μν of scalars and through the ratio of the determinants g/g. Nevertheless, as was thoroughly analyzed in [50], these two contributions, instead of adding, can cancel each other to keep the graviton massless, or equivalently, gravity unmodified. This does not mean that the metric tensors in (14) do not participate in other physical processes. Indeed, they can well generate our targeted structures involving the gravitational acceleration ∇φ. Consequently, we associate the metric tensors in (14) with the two phases of motion as g μν = 0 ⇒ Newtonian regime, keeping in mind that gravity is not necessarily massive. Indeed, the model of [50] offers a wide parameter space to set V 1 (4) = 0 in Eq. (26) and ζ V 1 (4) = 0 in Eq. (27). Moreover, potential terms in Eq. (11) give enough freedom to realize massless and massive gravity phases. Therefore, as will be proven below, the MOND regime can be realized by using the metrics in (14) without the necessity of modifying gravity.

Construction of a
Having fixed all the variables in (12), we now start formulating the acceleration scalar a. The kinetic term g μν g μν of the scalars does not contribute to ∇φ g , and hence, the argument of a in (12) represents the optimal list of dynamical variables. Out of various possibilities, we consider for a a simple structure, where all indices are raised and lowered with g αβ so that αβ is the trace of the matter energy-momentum tensor in the Newtonian domain. Here, the dimensionless constants c 1,2 will be fixed in the weak field limit by imposing (8). The presence of the metric determinants in (16) is crucially important for MOND because gravitational acceleration ∇φ g is generated by the derivatives of g/g (not g μν g μν , for instance).
Having fixed its functional form in (16), we now start checking if a satisfies its defining asymptotics in (2) and Table 1. This requires its evaluation in the two vacua in (14) since they correspond to the Newtonian and MONDian regimes as indicated in (15).
1. g μν = 0 and ∇ μ T (N ) μν = 0. In this vacuum, g vanishes identically and, as follows from (16), a becomes infinitely large thanks to the fact that c 1,2 > 0, as will be proven below. Now, having found a > 1, one gets μ (a) 1 and this gives T μν T (N ) μν from (7). Thus, the Einstein field equations (6) reduce to (17) in which consistency of the Bianchi identity on G μν is maintained by the conservation of T for a dust distribution characterized by the energymomentum tensor in (9). Clearly, this equation holds if the metric tensor takes the form as appropriate for the non-relativistic limit.
In conclusion, as conjectured in Eq. (15), the minimum energy configuration g gives rise to the Newtonian regime for motion. Small perturbations about this vacuum make g = 0 but this determinant is expected to be sufficiently small to secure the Newtonian regime, a > 1. 2. g μν = η μν and ∇ μ T (N ) μν = 0. In this vacuum, in the non-relativistic limit in which the metric tensor is given by (19), the acceleration scalar defined in (16) becomes for dust whose energy-momentum tensor is given partly by (9) and partly by extra interactions occurring in the low-acceleration regime. It is due to this alleged extra piece that T (N ) μν in (9) satisfies ∇ μ T (N ) μν = 0. It is clear that the acceleration scalar exhibits correct nonrelativistic limit if because then the last two terms of (20) drop out to enable the required limit in (8). Thus, the construction in (16) for a does indeed reduce to the acceleration of the point mass rather than the gravitational acceleration −∇φ g . The non-relativistic result in (20), which holds for a < 1 or equivalently |a| < a 0 , entails μ (a) a, so that the Einstein field equations (6) take the form where the scalar field Q is to be chosen judiciously to make the right-hand side have a vanishing divergence. This constraint, ensuring conservation of T μν , can be difficult to satisfy if Q does not involve T (N μν and g/g. As a plausible structure, we set where one can of course consider alternative structures giving similar results in the non-relativistic limit. In g μν = η μν vacuum, in the non-relativistic limit, conservation of T μν gives where the metric tensor is given by (19). This differential equation is too involved to suggest the MOND dynamics. Nevertheless, a closer look reveals that, if (i) the energy density ρ varies slowly in space (|∇ρ| ρ|∇φ g |) and if (ii) acceleration a varies slowly both in space and time (|∇ μ ρ| ρ|∇ μ φ g |), then one gets from (24) which is the desired MOND relation given in Eq. (4). 3. Non-conservation of T μν N . Having obtained motion equations in the two regimes of g μν , we now turn to a discussion of the non-conservation of T μν N . In view of the dis-cussions summarized in Table 1, the energy-momentum tensor T μν , introduced in (6) and defined in (7), is always conserved. This is necessary for the consistency of the gravitational field equations (6). The T μν N tensor, however, is conserved only in the Newtonian regime. To see how these conservation features hold, it proves useful to examine the divergence of T μν N , as follows from (6) with (7). It is obvious that, in the Newtonian regime, μ (a) → 1 and f ν N vanishes identically to ensure conservation of T μν N . In the MONDian regime, however, μ (a) → a = 1, and f ν N stays non-vanishing. This prohibits conservation of T μν N . These features are precisely the ones listed in Table 1. The MONDian force is consistent with (22). Since a is related to g μν as in (16), the second metric g μν turns out to be a fundamental ingredient of the entire formalism. Not surprisingly, effective forces similar to f ν N also arise in modified gravity theories which couple curvature and energy-momentum tensor T μν N directly [51][52][53].
In this section, we have succeeded to get the MONDian dynamics starting from (6) by defining the acceleration scalar a as in (16), the Q scalar as in (23), and the second metric tensor as in (13). Moreover, we have explicitly ensured conservation of the total energy-momentum tensor T μν , while determining effective MOND force associated with the nonconservation of T μν N . The analysis here provides an existence proof.

Conclusion and future prospects
In the present paper, we reported our results on relativistic MOND as derived from modified dynamics rather than modified gravity. Our approach is an empirical one and gives the beginning stage of a general investigation of relativistic MOND. The formalism developed, though lacking an action principle, can be regarded as generalizing Milgrom's modified inertia approach [38][39][40][41] to the relativistic domain. It is based on the energy-momentum tensor of matter. The reason for this is that the energy-momentum tensor of matter in the Newtonian regime, which necessarily loses its conservation property due to extra interactions occurring at sub-Hubble accelerations, seems to provide the correct path way to quadratic acceleration in the MOND regime. In fact, this dynamical structure cannot follow from other sources such as potentials, metric tensor, and curvature tensor. The main observation behind our approach is that matter possesses its usual energy-momentum tensor under the usual circumstances where the Newtonian laws hold. However, the same matter, at exceedingly small accelerations below the Hubble scale, develops novel interactions causing non-conservation of its energy-momentum tensor, and it is with these interactions that MONDian dynamics arises. Our empirical relativistic model is essentially a bimetric theory. However, our approach to the second metric tensor mimics models of the gravitational Higgs mechanism in which the vacuum expectation value of the second metric tensor equals the flat Minkowski metric, and it provides the requisite terms clearing the gravitational acceleration contributions to enable the quadratic acceleration piece needed for MOND.
The present study can be extended in various respects for rectifying and improving the present model.
• In the present work we have taken matter at the skirts of galaxies as dust. For an accurate analysis of the matter distribution, however, one may need to extend it to a perfect fluid and other forms of matter. • In obtaining the MOND equation of motion (25) we have neglected contributions from the spatial variation of ρ.
The situation can be improved by incorporating such terms from (24). The effect can be pronounced especially at the arms of spirals where the dust density changes sharply.
These points are currently under investigation in [54]. Last but not least, the present model would be grossly improved if an invariant action could be written. The alleged action, which must directly generalize Milgrom's modified inertia approach in [38,39] to relativistic velocities could be too complicated to construct due mainly to the presence of the fixed acceleration scale a 0 . It might necessitate a 0 to be included in relativistic transformations.