Non-relativistic limit of embedding gravity as General Relativity with dark matter

Regge-Teitelboim embedding gravity is the modified gravity based on a simple string-inspired geometrical principle: our spacetime is considered here as a 4-dimensional surface in a flat bulk. This theory is similar to the recently popular theory of mimetic gravity: the modification of gravity appears in both theories as a result of the change of variables in the action of General Relativity. Embedding gravity, as well as mimetic gravity, can be used in explaining the dark matter mystery since, in both cases, the modified theory can be presented as General Relativity with additional fictitious matter (embedding matter or mimetic matter). For the general case, we obtain the equations of motion of embedding matter in terms of embedding function as a set of first-order dynamical equations and constraints consistent with them. Then we construct a non-relativistic limit of these equations, in which the motion of embedding matter turns out to be slow enough so that it can play the role of cold dark matter. The non-relativistic embedding matter turns out to have a certain self-interaction, which could be useful in the context of solving the core-cusp problem that appears in the LambdaCDM model.


Introduction
In the attempts to explain some astronomical observations in the framework of General Relativity (GR), certain discrepancies arise, which can be solved by an assumption of the existence of an additional matter. This matter is called Dark Matter (DM) since it can not be directly observed. Such discrepancies appear at a wide range of scales: from galactic to Hubble ones, see, e.g., [1]. In the case of galactic, it is the deviation of the stars' motion from the GR predictions in the assumption that stars are moving only in the gravitational field of visible matter ("rotation curves"). In order to explain this discrepancy, an assumption is made that there is some amount of DM in galaxies besides ordinary matter. The density of DM decreases slower than ordinary matter density, so it forms a halo around a galaxy. At the larger scales, the observations of gravitational lensing, interpreted in the GR framework, indicate a significant amount of matter in the regions with a low amount of visible matter. This fact can also be explained using the assumption of the presence of DM in such regions. At the largest cosmological scale, the assumption of the homogeneous distribution of DM (in the amount of 26% of the Universe total mass-energy), alongside with the existence of dark energy (69%) whose properties are close to the cosmological constant, ensures the matching of average mass-energy density with the critical density, as well as provides the solution of the structure formation problem.
All currently available observational data (including the very important data related to the cosmic microwave background anisotropy) can be satisfactorily described by Λ-Cold Dark Matter (ΛCDM) model, which is a Standard Model of modern cosmology. In its framework, DM can be considered as non-relativistic dust matter, which produces the same gravitational field as ordinary matter, but non-gravitational interaction of dark matter with ordinary one is either absent or very weak and thus undetectable [2].
Numerous attempts to detect DM have been fruitless so far [3,4], which makes a mystery of the DM nature one of the most intriguing problems of modern physics. There are many hypotheses (see, e.g., [5][6][7]) related to the nature of DM, the choice of which determines a way to detect it. The failures of detecting can be connected to the weakness of the coupling between the dark and the ordinary matter. However, a more radical explanation exists. Since DM undoubtedly interacts with gravity, the absence of its participation in other interactions can be interpreted as an indication of the fact that DM is not a real matter coupled with Einsteinian gravity, but rather is an effect of a description of the gravitational interaction in the framework of a certain modified theory of gravity. In such an approach, DM turns out to be absent from the point of view of the fundamental theory, but it reveals its existence through interaction with gravity (via a contribution to Einstein equations) if the latter is described in the framework of GR instead of its description in the context of the modified theory.
One possible candidate for the role of this modified theory is MOND [8]. However, modifications that have additional degrees of freedom in comparison with GR seem more promising (see the review of modifications in [9], see also [10]). In consideration of modified theory in the context of GR, these degrees of freedom turn out to correspond to DM. In recent years strong attention has been drawn to the mimetic gravity [11], which allows one to introduce "dark matter without dark matter" [11] in the described manner. This model appears as a result of the change of variables, which isolates the conformal mode of the metric. Since this change of variables involves differentiation (being a particular case of the differential field transformations in GR [12]), it leads to the extension of dynamics of the theory, and some additional degrees of freedom w.r.t. GR appear. These degrees of freedom describe the mimetic matter, which might be interpreted as DM. The possibility to reformulate the mimetic gravity as the GR with an additional contribution to the action [13] shows that the mimetic matter is a potentially moving dust effectively described by two scalar fields. The restriction of the potentiality of motion can be lifted [14], the introduction of a pressure is also possible [15], but the theory nevertheless requires significant extensions to resolve the discrepancies mentioned above in the framework of mimetic gravity, see, e.g., [16] and also review [17].
Another example of modified theory, whose origins lie in the procedure of differential field transformations in the GR, is Regge-Teitelboim embedding gravity (also known as embedding theory) [18]. It appears as a result of the following change of variables: (here and hereafter µ, ν, . . . = 0, 1, 2, 3) in the Einstein-Hilbert action of GR. This change of variables, in contrast with the mimetic one, has a deep geometric nature, since (1) is a formula for induced metric on the 4-dimensional surface described by embedding function y a (x µ ) in an ambient spacetime with Minkowski metric η ab . Therefore the embedding gravity is based on a simple string-inspired assumption that our spacetime is not merely an abstract pseudo-Riemannian space but rather a surface in a flat ambient spacetime of higher dimension. In accordance with the Janet-Cartan-Friedman theorem [19], the minimal number of dimensions of a spacetime in which an arbitrary 4-dimensional surface can be locally embedded, is equal to ten, so the ambient spacetime in embedding gravity is usually taken to be 10-dimensional.
After the appearance of paper [18] the ideas of embedding gravity were discussed in the paper [20] and were used lately by many authors with the aim of description of gravity, including its quantization (see, e.g., [21][22][23][24][25][26][27][28][29][30]). In this framework, one can try to explain the effects of DM in the assumption of a strict spatial homogeneity and isotropy of the universe [29,31]. To do it at the smaller scales, one might try to reformulate the theory as the GR with an additional contribution corresponding to some fictitious matter [32] and to examine the properties of this embedding matter. The number of variables that describe embedding matter turns out to be sufficiently large, and the properties of this matter are sufficiently rich to hope for a successful description of observational properties of DM in the framework of embedding gravity, without the introduction of additional extensions (as it is done, e.g., in the framework of mimetic gravity).
The fact that embedding gravity is based on the simple geometric idea without any ad hoc extensions and arbitrary constants , as well as the potential quantization advantages connected with the presence of flat ambient spacetime in the formulation of the theory (see, e.g., [21,24]), makes this modification of gravity, in our opinion, an especially promising in the search for an explanation of the mystery of DM.
In the present paper, we will analyze the equations of motion of embedding matter in the non-relativistic limit, whose definition was proposed in the paper [33]. In this analysis, we will follow the paper [33] in part, but the calculations will be performed in a more straightforward and geometrically clear manner. In the [33] a non-square vielbein e a µ ≡ ∂ µ y a was used instead of embedding function y a (such an approach was also used in another variant of embedding gravity [28]), and such vielbein was required to satisfy the integrability condition ∂ µ e a ν = ∂ ν e a µ . Due to this choice of independent variables, the embedding function y a , which describes 4dimensional surface, appeared in the calculations only in an implicit way. In the present paper, we show a straightforward way to obtain the same result, using the embedding function y a as one of the variables describing embedding matter.
In Section 2, we remind the formulation of embedding gravity, both in its original form, based on the differential field transformations in GR, as well as in the form of GR with additional embedding matter contribution in the action, and write down its equations of motion. In Section 3, we analyze and rewrite these equations in the form of dynamical equations (with first-order time derivatives only) for the independent variables describing embedding matter and constraints that restrict the choice of these variables at the initial moment of time. In Section 4, we consider the equations of motion in the limiting case of non-relativistic approximation, whereas in Section 5 the deviations from this limit are considered.
We succeeded in showing that the embedding matter in the studied limit turns out to be a non-relativistic matter with certain self-interaction. If this self-interaction is not very strong, this behavior quite well corresponds to the known properties of DM at all scales, but the very presence of self-interaction could be useful. For example, the models in which DM is the Self-Interacting Dark Matter (SIDM) [7] can be considered as solutions to the well-known core-cusp problem. This problem manifest itself in the fact that according to calculations the density of DM without self-interaction must increase rapidly towards the center of a galaxy, forming "cuspy" profiles. However, the observations do not confirm this, giving "cored" profiles.

Embedding gravity
Regge-Teitelboim embedding gravity (or embedding theory) is the theory of modified gravity appearing as a result of the change of variables (1) in the action of GR: where S m is a contribution of an ordinary matter. In this approach the independent variable describing gravity is an embedding function y a (x µ ) (we remind that µ, ν, . . . = 0, 1, 2, 3; a, b, . . . = ∅, 1, . . . , 9; we also assume that the signature of the flat ambient spacetime is (+ − . . . −), so the metric g µν has the signature (+ − −−)). The variation of (2) w.r.t. independent variable y a leads to the Regge-Teitelboim (RT) equations [18] where D µ is a covariant derivative, G µν is the Einstein tensor and T µν is the EMT of the ordinary matter. In order to simplify the analysis of the embedding gravity from the familiar general-relativistic point of view, one can rewrite the RT equations (3) as a set of equations [21] G µν = κ (T µν + τ µν ) , The first one can be interpreted as the Einstein equation with an additional contribution τ µν to EMT from a certain embedding matter, whereas the second one plays the role of the equation of motion for this fictitious matter. The equations (5) resembles a local conservation law for a set of some currents: There is a possibility to reformulate the embedding gravity as GR with some additional matter not only at the level of equations of motion but also at the level of action. To do this, one needs to add certain action of embedding matter S add to the action (2): The metric here is assumed to be an independent variable, and the condition (1) appears as a consequence of the equations of motion. In the paper [32] several ways to choose S add were considered. Let us introduce two of them, which seem the most promising in the context of the analysis of embedding matter properties. First of all, one can choose the symmetric tensor τ µν alongside with embedding function y a as independent variables describing embedding matter [14]: so the independent quantity τ µν , which is playing the role of a Lagrange multiplier, turns out to be an EMT of embedding matter. As an alternative, one can assume the independence of the set of currents j µ a in the equations (6) instead of τ µν (y a is still independent) [32]: Here it is assumed that one needs to take the square root of matrix with the indices µ and α and subsequently take the trace of the result (the operation tr). In terms of j µ a the EMT of embedding matter can be written as where the square root of matrix is again present, see details in [32]. Note that (11) is consistent with (7). The set of currents j µ a can be interpreted as four vectors (indexed by µ) of the 10-dimensional flat bulk. In Sections 4 and 5 we will analyze the behavior of embedding matter in such case that all these four vectors are non-relativistic in the bulk [33], i.e.
As we said in the Introduction, the analysis will be performed partially in spite of the paper [33], but with more straightforward and geometrically clear way of calculations. If the expression (9) is used as a contribution of embedding matter, the full set of its equations of motion has the form: and the metric g µν needs to be found through a separately considered Einstein equation (4) (the connection Γ α µν that is present in the covariant derivative D µ , is represented by Christoffel symbols which are expressed through the metric). The form of equations of motion of embedding matter is different in the case where one takes the set of currents j µ a as independent variables instead of τ µν , and (10) is chosen as an action of the embedding matter, see [32].

Equations of motion of embedding matter
Let us note that the satisfaction of (13) leads to the relation that is well known in the formalism of the embedding theory (see, e.g., [25]): It means that the second fundamental form of the surface is transverse w.r.t. index a, i.e. orthogonal to the surface. Using (15), one can rewrite (14) as a set of two equations: Let us consider the equations of motion of embedding matter (13), (17),(18) as a time x 0 dynamics of the independent variables y a , τ µν . Let us also recall another known formula of the theory of embedded surfaces, which connects the curvature tensor R αβµν and second fundamental form of the surface b a µν (see, e.g., [25]): Taking α = µ = 0, β = i, ν = k in it (here and hereafter i, k, . . . = 1, 2, 3) and using (16), we obtain the corollary of (13): The components of the second fundamental form of the surface b a ik can be treated as a 6 × 6 matrix, since the pair of indices {ik}, w.r.t. which this quantity is symmetric, have six possible values, and the index a also can take only six values (due to the transversality condition (15)) without the vanishing of the quantity. In the generic case, where this matrix is invertible, one can use an "inverse" matrix, introducing the quantity α ik a , which is uniquely defined by the relations Using this quantity, one can express the transverse part of D 0 ∂ 0 y a from (20), and since the tangent part of it vanished due to (15), one finds that We used only one equation of motion in the derivation of the relation (22), namely (13). Using (22) in the another equations of motion (18), we obtain new equation which now can be used instead of (18). Using the quantity α ik a again, one can extract the components of τ ik from this equation: We are therefore able to express the components of τ ik through other variables, excluding it from the set of independent variables. As a result, only y a and τ 0ν remain independent variables describing embedding matter, and the equations of motions of this matter are reduced to (13) and (17) with (24), i.e. to the system of equations The equation (25), in contrast with remaining equations (26)-(29), does not contain time derivatives. In the canonical formulation (different forms of it for embedding theory was studied in the works [18,22,25,26,[34][35][36][37]) such equations are called constraints. In the examination of constraints, after the discovery of primary constraints one usually needs to check their conservation in time in order to find secondary constraints. Let us find the conservation condition for the constraint (25), using the covariant time derivative D 0 instead of ordinary one ∂ 0 for the sake of convenience. These are equivalent, since change from ∂ 0 to D 0 leads only to the appearance of contribution proportional to a constraint. We obtain: where we used that D 0 ∂ i y a = D i ∂ 0 y a . Using the product rule in (30) and taking (26) into account, we obtain the condition of the conservation of constraint (25): (note that this equation is a particular case of the abovementioned relation (15)). Using the expression for a covariant derivative as well as relation (25), (26), one obtains In the generic case, when the three vectors (at different values of k) ∂ k y a together with six vectors (at different values of symmetric indices i, k) ∂ i ∂ k y a form a basis of nine vectors of the ambient space, the equations (26) and (32) together with normalization equation (27) allow to uniquely determine the quantity ∂ 0 y a if the function y a (x 0 , x i ) at the given moment is known. Such expression defines the dynamics of the embedding function y a , providing the satisfaction of (26) and (27), as well as the conservation of the constraint (25) that must be imposed when one chooses y a at the initial moment of time x 0 . The remaining equations (28), (29) define the dynamics of the remaining independent variables τ 00 and τ 0m , so the equations of motion of embedding matter are reduced to the constraint (25) and dynamical equations for independent variables y a and τ 0ν which expresses their first time derivative through their values at the given moment of time.
In the general case dynamical equation for y a turns out to be quite bulky, but it can nevertheless be written in an explicit form. To do that, let us notice that at the fixed x 0 the embedding function y a (x 0 , x i ) defines a certain 3-dimensional surface in the ambient spacetime, so one can define a second fundamental form of it: Using it, let us introduce a quantity 3 α ik a , defined by relations close to (21) [25]: Contracting the quantities 3 α ik b and 3 b a lm the other way, we obtain a projector on a corresponding 6-dimensional space: where the second and the third terms at the r.h.s. are projectors on 3-dimensional and 1dimensional spaces, correspondingly.
Using (26), let us rewrite (32) in the form Multiplying it by 3 α ik b and using (36), we obtain Using it together with (26), we can write where we denote ζ = ω c ω c ω a ∂ 0 y a . Taking the properties (34), (35) into account, one can find the quantity ζ by substitution of (39) into (27): As a result, we obtain the explicit form of the dynamical equation for the independent variable y a in the form (39) with (40). Together with dynamical equations (28)

The limiting case of non-relativistic approximation
Let us analyze the equations of motion of embedding gravity in the non-relativistic in the bulk limit (12). We will assume that j µ a and τ µν (as well as the EMT of matter T µν ) are small enough, so due to the small value of the gravitational constant κ the r.h.s. of Einstein equations (4) is small and gravitational field turns out to be weak: where η µν is a flat metric. Let us firstly find a solution in a limiting case when δj µ a = 0 (see (12)) and h µν = 0. The corresponding quantities will be marked by a bar.
The limiting value of the set of currents can be written as where the arbitrary vector j µ that appeared in (12) (it must be timelike for the action (10) to be real) is presented in the form Substituting (42) in (11), we obtain that Substituting (42) and (45) in (7), we notice that where we have assumed thatρ τ = 0 (the caseρ τ = 0 corresponds toτ µν = 0, i.e. the contribution of embedding matter is absent, so this case is of no interest in the context of the analysis of its properties). One can always choose coordinates x µ on a surface defined by embedding functionȳ a (x) in such a way thatȳ 0 = x 0 . The relation (46) means that the vector δ ∅ a is tangent to the surface at any point of it so that one can impose an additional restriction on a coordinate system: (here and hereafter I, K, . . . = 1, . . . , 9). Therefore we can conclude that there is a certain coordinate system such that the embedding functionȳ a (x) in the considered limit has the form On the other side, there is a certain coordinate system (possibly not the same as above, generally speaking) in which the relation (44) is satisfied, so the curvature tensor of the surface vanishes. If a 3-dimensional submanifold x 0 = 0 of the surface (48) is considered (ȳ I (x i ) is its embedding function to 9-dimensional space), then, using the Gauss equation (it connects 3-dimensional components of 4-dimensional curvature tensor with 3-dimensional curvature tensor and the second fundamental form of the considered 3-submanifold of Riemannian 4-space, which is equal to zero in the present case), it can be noticed that its curvature is also equal to zero. One can therefore transform the metric of the discussed 3-submanifold to the formḡ ik = −δ ik by choosing the coordinates x i . After that the metric of 4-surface becomes to coincide with (44), which can be easily seen using (1). Therefore there are coordinates in which the relations (44) and (48) are simultaneously satisfied, and we will use such coordinates. Substituting (48) in (46) at a = 0, we obtain thatū 0 = 1. Using this fact in (43) together with (44), we find thatū According to (45), this leads to the expression for EMT: which corresponds to the resting dust embedding matter with densityρ τ . As it was noticed after (48), the quantityȳ I (x i ) turns out to be an embedding of flat 3-metric in a 9-dimensional Euclidean space. Since generic 3-metric is parametrized by six independent components, such embedding can be parametrized by three arbitrary functions. In a general case (which, for example, has to be realized if we consider the equations of motion of embedding matter as a dynamics w.r.t. time x 0 and do not perform a fine-tuning of initial values) this embedding has such form that there is a quantityᾱ ik a defined by (21) which corresponds to (48), andᾱ ik ∅ = 0. Such embedding structure corresponds to the notion of free embedding introduced in [38] in context of 3-dimensional manifold, and to the notion of spatially free embedding in context of 4-dimensional one. Note that in the considered limit, when the embedding function has the form (48), the quantities defined by relations (21) and (35) which means that the embedding matter density does not change with time. In this caseρ τ can arbitrarily depend on spatial coordinates x i ; this dependence is determined by the choice of initial values. It can also be noticed that according to (40) we haveζ = 1.

Equations of motion in the nonrelativistic limit
Now let us analyze the equations of motion in such regime that the independent variables are close to their limiting values considered in the previous Section. In this case the metric has the form (41). We will additionally assume that ordinary matter is a dust with density ρ and non-relativistic motion, so its EMT in the leading order has the same form as (50): and ∂ 0 ρ is small. Taking into account that EMT of embedding matter τ µν is close to (50) and using harmonic coordinates, as it is usually done in such cases, one can write down linearize Einstein equations (4) as follows: It leads to the fact that in the leading order h µν is diagonal and their components are proportional to κ, whereas its non-diagonal elements are of higher order, so, in particular, h 0i is small in comparison with κ. Note also that the small value of ∂ 0 ρ together with (52) leads to the slow change of gravitational field, so ∂ 0 h µν is also small in comparison with κ. If we consider a finite time interval ∆x 0 only, then the embedding function y a can be obtained as a small deviation from its limiting value (48): where δy a is small. Note that the expression (48) gives an explicit expression of the component y 0 only, while for remaining components, it merely means that they are time-independent according to the condition (47). In the case that is close to the limiting one, this condition is replaced by the smallness of ∂ 0 y I . Therefore, if the behavior of independent variables is studied not only at finite, but also at large time intervals; in particular, at those intervals, ∆x 0 whose inverse is proportional to the smallness parameter (we will do it below in the transition to the non-relativistic limit), then y I could start to change noticeably, and it turns out to be impossible to write it down as (55) with a small δy I . For that reason, we will assume that the only component that can be written as its limiting value plus a small deviation is y 0 , and the remaining components are restricted only by the condition that their derivatives w.r.t. x 0 are small: where w is small. For finite intervals ∆x 0 , such formulation of conditions corresponds to the expression (55) with small δy a and, at the same time, allows components y I to change noticeably at the large (i.e., of the order of inverse values of ∂ 0 y I components) intervals ∆x 0 . Note that we do not choose any particular coordinate system when writing the embedding function in the form (56), so the metric g µν can be defined in any coordinates. This fact was used above (54) when harmonic coordinates had been chosen.
In the considered equations of motion of embedding matter (namely, in (29) and (39)) the quantities α ik a and 3 α ik are present, which are defined by relations (21) and (35). As was noted above (51), in the limiting case (and hence in the leading approximation, when we consider the deviations from the limiting case) these quantities are equal to each other and vanish at a = 0. We will denote their values in the leading approximations asα ik a , using the following exact definition of it: which is the limit of (21) and (35). Defined as such, the leading approximation of the quantities α ik a and 3 α ik a coincides with time-independent quantityᾱ ik a for the embedding function (48). However, after its replacement by (56) this coincidence is no longer exact due to the appearance of a slight dependence of y I on x 0 . As a result,α ik a can also be slightly time-dependent (by slightness we mean that ∂ 0α ik a is small due to the satisfaction of the condition from (56)). Let us start the analysis of the equations of motion of embedding matter from the equation (39) which defined a dynamics of embedding function. To do that, let us write down the formula (36) at a = 0 and b = I using (33): where we took into account that ω a ω a = 1 (59) according to (34),(51) and the small value of deviations. It follows from (58) that with the use of (41), (51) and (56), in the leading approximation one can find: Now let us find the leading approximation of (39), assuming that κ, as well as w, is small, and we can neglect higher powers of them. Using (40), (41), (56) and (59) together with the fact that h 0i and ∂ 0 h µν are small in comparison with κ (see below (54)), in such approximation we obtain where the quantity ω I is given bt (60) and has the same order of magnitude as w.
Writing down the remaining dynamical equations (28), (29), at the leading order in a similar manner, we obtain where we denote ρ τ ≡ τ 00 . We used (56),(57),(61),(62) and took into account that h 0i , ∂ 0 h µν ∼ o(κ). Let us examine the time evolution of independent variables governed by equations (61)-(64). Consider the situation when at the initial moment the variables w a τ 0i (as well as ω I according to (60)) are either vanished or very small. While this condition is satisfied well enough, the contributions at the r.h.s. of (61) and (64) containing these quantities can be neglected in comparison with remaining contributions of order κ. The time dependence of the latter can also be neglected, so with these assumptions, the r.h.s. of (61) and (64) will also be timeindependent in the leading order, and w, as well as τ 0i , will increase linearly w.r.t. x 0 . In a finite time interval ∆x 0 they will necessarily reach the order of magnitude of κ, but despite this the contributions in r.h.s of (61), (64), in which they are present, will be nevertheless small in comparison with leading contributions of order κ, so the linear increasing of w and τ 0i will be continued. Consequently, according to (62) and(63), the quantities y I and ρ τ will increase quadratically with time x 0 , i.e. will receive an additive perturbation proportional to x 0 2 .
Such trivial (i.e. corresponding to a linear or quadratic time dependence) dynamics would take place during an arbitrary finite interval of time ∆x 0 , while the quantities w and τ 0i have order of κ. However, for a time ∆x 0 ∼ 1/ √ κ (which is not assumed to be finite anymore, since κ is small), these quantities would increase to the order √ κ. After that, it would become impossible to drop out the contributions at the r.h.s. of equations (61) and (64), which were neglected earlier. Note also that after non-finite time intervals ∆x 0 , it is no longer possible to assume time independence of the contributions to the r.h.s. of (61) and (64), that were initially of the order of κ.
Therefore the dynamics defined by equations (61)-(64) becomes nontrivial at time scales of order ∆x 0 ∼ 1/ √ κ. In order to analyze this nontrivial dynamics (which can easily be recognized as a standard non-relativistic one) let us perform a time change, introducing a large quantity c: One can assume that c is the speed of light, then, taking into account that it must be assumed that Newtonian gravitational constant G is a finite quantity, whereas the small value of κ corresponds to the large value of c according to (66). Using the fact that ∂ 0 h µν is small in comparison with κ (see below (54)), one can obtain from (54) that at the leading order where ϕ is the Newtonian potential corresponding to matter distribution with density ρ + ρ τ : As it was mentioned above, in the nontrivial regime the independent variables w and τ 0i have order of √ κ ∼ 1/c. Due to this, in order to analyze non-relativistic limit, let us introduce some new quantities instead of them: which remain finite at κ → 0, c → ∞. Let us also introduce instead of ω I , where (60) was used. In terms of these new variables and derivatives w.
the equations (61)-(64) can be rewritten up to terms of order 1/c as follows: where we used (67), (69) and (70). The resulting dynamical equations (72)-(75) describe the motion of the embedding matter in the non-relativistic limit c → ∞, and, together with (67),(68) (to which the Einstein equations are reduced in the considered case) and the equations of motion of ordinary matter fully describe the non-relativistic dynamic of a system. The possibility to construct consistent non-relativistic equations means, in particular, that non-relativistic approximation, appearing as a result of non-relativistic character of currents (12) in the ambient space, turns out to be stable: the embedding matter remains non-relativistic.
The dynamic equations (72)-(75) need to be supplemented by imposition of the constraint (25) at the initial moment, which (if (41),(56),(69) are taken into account), up to terms of order 1/c, takes the form This equation means that in the non-relativistic limit (we remind that the time conservation of the constraint is provided by dynamical equations, see below (32)) y I resembles an embedding function of flat 3-dimensional metric in 9-dimensional ambient space at any moment of time. As we mentioned above (see the text below (50)), such embedding is parametrized by three arbitrary functions. Unfortunately, we did not succeed in the construction of an explicit form of such functions. These functions evolve with time according to the dynamical equation (73), which can be interpreted as dynamical isometric bending of a 3-dimensional surface with a flat Euclidean metric. It can be straightforwardly checked that the deformation governed by (73) is indeed isometric bending. We obtain that the embedding matter in the non-relativistic limit is described by the following eight parameters: • three functions which parametrize an embedding function of a flat 3-dimensional Euclidean metric in a 9-dimensional Euclidean space, • the quantity ψ (see (69) and (56)), • the density ρ τ of the embedding matter, • the velocity v i τ of the embedding matter.
To obtain an equation that describes the velocity v i τ dynamics, (75) needs to be rewritten using (74). Let us do that, together with the formulation of this equation in a familiar way of dust matter description, where it is not just ∂ t v m τ that is present at the l.h.s. of the equation, but rather some combination resembling an acceleration of a "single particle" of the embedding matter (although we assume that in the fundamental theory such particles do not exist, and all variables describing the embedding matter are, in fact, gravitational variables in the framework of embedding gravity). As a result, we obtain: The non-relativistic equations of motion of the embedding matter in the form (72)-(74),(77) reproduce the ones that were obtained by the more complicated (and less geometrically clear) way in the work [33]. The equation (74) is a usual conservation law of this matter, and (77) defines the force (per unit volume) acting on its "single particles". Its first term is a well-known gravitational force corresponding to Newtonian approximation. The remaining terms can be interpreted as a certain self-interaction force of the embedding matter, which depend not only on familiar characteristics of this matter, namely its density ρ τ and velocity v i τ , but also on its additional characteristics ψ and y I , whose dynamics is governed by equations (72) and (73).

Embedding function in the non-relativistic limit
Let us find the approximate expression for the embedding function y a (x) corresponding to the non-relativistic limit that was found in the previous Section. According to (72) and (73), the quantities ψ and y I have certain dependence on t = x 0 /c governed by these equations. Therefore, according to (56) together with (69), the components of the embedding functions can be written as (note that in this notation the condition (56) is satisfied), where the quantities ψ(t, x i ), y I (t, x i ) satisfy the equations that are reduced to (72), (73) and (76) in the limit c → ∞. The limiting forms of these quantities that exactly satisfy the mentioned equations, are denoted asψ(t, x i ), y I (t, x i ) in this Section, so one can write where δψ(t, x i ) and δy I (t, x i ) are small corrections. According to (76),ŷ I (t, x i ) at each value of t is an embedding function of a flat 3-dimensional Euclidean metric in 9-dimensional Euclidean space, so the corresponding induced metric exactly coincides with δ ik . Note that an arbitrary perturbation of this induced metric can be obtained by choosing the embedding function de-theory in the form of GR with an additional contribution of the embedding matter, we find its equations of motion in the general case. These equations are reduced to a set of first-order dynamical equations (28), (29),(39) and the constraint (25) that is in involution with them. One can try to explain various effects in terms of embedding matter by analyzing these equations in different regimes. Such effects could be related not only to DM (which is the topic of the present work) but also probably to the dark energy and the inflation. Since observational discrepancies, which exist in the framework of GR, can be quite well explained by an assumption that DM is a cold dust matter, we consider the class of solutions of embedding gravity, which corresponds to the non-relativistic character of motion of the embedding matter. This class of solutions naturally appears due to the assumption that conserved currents j a µ describing the embedding matter are non-relativistic in the bulk (12). The obtained non-relativistic regime of motion of fictitious embedding matter turns out to be stable in such sense that it remains non-relativistic during the evolution if the choice of initial values corresponds to the non-relativistic regime. Especially interesting is the fact that in the non-relativistic regime, besides usual coupling with gravity, this embedding matter possesses certain self-interaction, as in the models of Self-Interacting Dark Matter [7], which could be perspective in the context of solving the core-cusp problem appearing in the ΛCDM model.
In the present work, we use the method of obtaining non-relativistic equations of motion, which has a simpler geometric meaning than the one that was used in the work [33], since from the start to the end of the calculations, we use 4-dimensional embedding function y a (x µ ). It should be noted that in the present work, as well as in [33], the choice of the spatial part of an embedding function is essentially different from the one that was used at the cosmological scales in the works [29,31], where it has minimal codimension (i.e., the 3-dimensional surface was embedded in 4-dimensional flat bulk). On the contrary, in the present work, we suppose that for a random choice of initial values, this 3-dimensional surface turns out to be "unfolded" in 9-dimensional flat bulk, which corresponds to the existence of the quantityα ik I defined in (57), so the 4-dimensional embedding is spatially free [38].
In the non-relativistic limit, the embedding matter is described by eight fields: besides its density ρ τ and velocity v i τ an additional field ψ is present alongside with an embedding function y I of a flat 3-dimensional Euclidean metric in 9-dimensional Euclidean space, which form is described by three functions. These four additional degrees of freedom have their own dynamics and affect the process of self-interaction of the embedding matter. An important role in the embedding matter's behavior is played by the initial values of these variables, which correspond to their values at the beginning of the non-relativistic regime. These values are formed during the preceding relativistic period, which possibly corresponds to the inflationary era.