Next Article in Journal
Bouncing Cosmology in Modified Gravity with Higher-Order Gauss–Bonnet Curvature Term
Next Article in Special Issue
Cooling Process of White Dwarf Stars in Palatini f(R) Gravity
Previous Article in Journal
QCD Axion Kinetic Misalignment without Prejudice
Previous Article in Special Issue
Galactic Wormhole under Lovelock Gravity
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Weak Field Limit for Embedding Gravity

1
Saint Petersburg Department of Steklov Mathematical Institute of RAS, 191023 Saint Petersburg, Russia
2
Department of High Energy and Elementary Particle Physics, Saint Petersburg State University, 199034 Saint Petersburg, Russia
*
Author to whom correspondence should be addressed.
Universe 2022, 8(12), 635; https://doi.org/10.3390/universe8120635
Submission received: 30 October 2022 / Revised: 21 November 2022 / Accepted: 24 November 2022 / Published: 29 November 2022
(This article belongs to the Collection Modified Theories of Gravity and Cosmological Applications)

Abstract

:
We study a perturbation theory for embedding gravity equations in a background for which corrections to the embedding function are linear with respect to corrections to the flat metric. The remaining arbitrariness after solving the linearized field equations is fixed by an assumption that the solution is static in the second order. A nonlinear differential equation is obtained, which allows for finding the gravitational potential for a spherically symmetric case if a background embedding is given. An explicit form of a spherically symmetric background parameterized by one function of radius is proposed. It is shown that this function can be chosen in such a way that the gravitational potential is in a good agreement with the observed distribution of dark matter in a galactic halo.

1. Introduction

Embedding gravity (also called embedding theory) [1] is a modified theory of gravity based on the idea of our spacetime as a 4D surface in a 10D flat space. From a geometric point of view, this approach renders the theory of gravity similar to string theory. In this case, the space–time metric is considered to be induced and has the form
g μ ν = ( μ y a ) ( ν y b ) η a b ,
where η a b is a flat metric of the ambient space, and y a ( x μ ) is an embedding function defining the shape of the surface. Here and further μ , ν , = 0 , , 3 and a , b , = 0 , , 9 . After the pioneering work of [1], the ideas of the embedding approach were repeatedly used in the works of various authors to describe gravity, including the problem of its quantization; see, for example, [2,3,4,5,6,7,8,9,10,11].
The embedding gravity equations of motion, called Regge–Teitelboim equations, are more general in comparison with Einstein’s equations G μ ν = ϰ T μ ν and have the following form:
G μ ν ϰ T μ ν b μ ν a = 0 ,
where b μ ν a is a second fundamental form of a 4D surface defined by an embedding function y a ( x μ ) . All “Einsteinian” (i.e., satisfying Einstein’s equations) solutions satisfy the Regge–Teitelboim equations, but they do not exhaust the full set of solutions. There may be additional solutions for which G μ ν ϰ T μ ν ; nevertheless, (2) is satisfied. Initially, this was seen as a disadvantage of the approach, since embedding gravity was understood as some reformulation of GR, potentially more convenient for quantization due to the presence of a flat ambient space. Therefore, it was proposed to additionally introduce Einstein constraints into the theory, rendering it truly equivalent to GR [1,12]. However, the presence of extra solutions in the theory can be considered to be an advantage.
If we assume that just the Regge–Teitelboim equations describe the gravitational interaction, from the point of view of describing observations in terms of the usual GR, extra solutions manifest themselves as an additional contribution to the right side of Einstein’s equations. It can be interpreted as the presence of some additional fictitious embedding matter that has nothing to do with ordinary matter. This contribution depends solely on the gravitational variables (which, in embedding theory, are y a ( x μ ) ), and its appearance is connected only with an attempt to reformulate a new theory in the old language. Since the direct detection of dark matter does not currently yield results [13,14], it is of interest to try to treat the embedding matter of the embedding theory as dark matter (and possibly as dark energy). Then, the effects associated with dark matter and energy are purely gravitational. In this way, observational problems that have no explanation within the framework of GR can be solved. The Regge–Teitelboim Equation (2) can be rewritten as a pair of equations
G μ ν = ϰ T μ ν + τ μ ν ,
τ μ ν b μ ν a = 0 ,
i.e., in the form of a set of Einstein’s equations with the contribution of an energy–momentum tensor τ μ ν of fictitious embedding matter and its equations of motion. This was first indicated in [15]. There are also modified theories of gravity which are alternative to the embedding gravity, equations of motion of which can be written in the form of Einstein’s equations supplemented by other equations. The most famous is mimetic gravity [16,17].
In order to obtain the properties of the embedding matter generated by embedding gravity (and thereby to understand whether they are similar to the properties of dark matter), it is necessary to investigate solutions of Equation (2). In general, due to their nonlinearity, this is too difficult a mathematical problem. So, we limit ourselves to the most physically interesting case of weak gravity, when metric g μ ν is close to flat metric η μ ν . In this case, the problem arises in choosing a background value y ¯ a ( x μ ) of the embedding function that would correspond to the flat metric. The simplest choice in the form of y ¯ a ( x μ ) , which defines a 4D plane in the ambient space, is unsuitable since, in such a background, Regge–Teitelboim Equations (2) are not linearized [2]. To linearize the equations, it is necessary to choose "unfolded" embedding [18] as the background, which means that the second fundamental form of the surface is nondegenerate in some sense.
In this paper, we use the unfolded embedding of the Minkowski metric, which is the product of a timelike line y I = c o n s t (we use indices I , K , = 1 , , 9 ; i , k , = 1 , 2 , 3 ) on 9D unfolded embedding y ¯ I x i of the euclidean 3D metric, i.e.,
y ¯ a = x 0 y ¯ I x i , i y ¯ I k y ¯ I = δ i k .
With such a choice of background embedding, a nonrelativistic motion of embedding matter is possible [19].
The purpose of this work is to study Regge-Teitelboim Equation (2), linearized in the background of (5). We look for a solution that corresponds to a galaxy that rotates so slowly that the effect of the rotation can be neglected, i.e., the distribution of ordinary matter is static and spherically symmetric on average. At the same time, we assume that metric (1) is also static (and spherically symmetric), which corresponds to the time independence of value τ μ ν describing embedding matter. The resulting solution determines the dependence of the gravitational potential on the distance to the center of the galaxy, and it can be compared with observations of the rotation curves of galaxies.
In Section 2 and Section 3, we obtain linearized equations and find their solution. In Section 4, the influence of the assumption of the exact static nature of the solution on its behavior in a linear approximation is investigated. In Section 5, we study how the problem is simplified in the case of spherical symmetry. We propose an explicit form of a spherically symmetric background embedding in Section 6. In Section 7, we study the possibility of choosing this embedding in such a way that the corresponding gravitational potential is in agreement with the observed rotation curves of galaxies.

2. Linearization of Regge–Teitelboim Equations

Let us recall some formulas of embedding theory. All convolutions by Latin indices are carried out using the flat metric of ambient space η a b . The induced metric is expressed in terms of the embedding function by Formula (1). We use space–time signature ( , + , + , + ) . The signature is changed by changing the sign of η a b , and the induced metric changes the sign as a consequence. The second fundamental form of a 4D surface is expressed in terms of the covariant derivative D μ of the embedding function, consistent with the metric or through projector Π b a on a subspace transverse to the surface (see, for example, [12]):
b μ ν a = D μ ν y a = Π b a μ ν y b .
We mark with a line the values corresponding to background embedding function y ¯ a (5); for example, b ¯ μ ν a is the second fundamental form of the background surface. We raise and lower the 4D indices of values with a line using the background metric. Since it is flat, the background connection is zero, and the covariant derivative in (6) is reduced to the usual one. The second derivatives of y ¯ 0 and the derivatives by x 0 of y ¯ I are zero; hence, the nonzero components of b ¯ μ ν a are:
b ¯ i j I = i j y ¯ I .
Through index I, this value is transverse to the 3D surface in 9D space described by embedding function y ¯ I ( x i ) . Thus, index I at each point effectively runs through 6 values, exactly like multi-index { i j } . Therefore, b ¯ i j I can be represented as a 6 × 6 matrix. Since this matrix is nonsingular for “unfolded” embeddings (see details in [18]), you can introduce symmetric by l , m and transverse by I value α ¯ I l m , which is inverse to b ¯ i j I in the matrix sense:
α ¯ I l m b ¯ i j I = 1 2 δ i l δ j m + δ j l δ i m .
We look for solutions of the equations in the form of Systems (3) and (4) as a perturbation theory series in degrees of the gravitational constant ϰ in the background (5):
y a = y ¯ a + ϰ y ( 1 ) a + ϰ 2 y ( 2 ) a + .
We write similar expansions in a series of ϰ for any other quantities, including τ μ ν and the energy–momentum tensor of ordinary matter T μ ν . Substituting expansions into Systems (3) and (4), in the first nonvanishing order, we obtain
G ( 1 ) μ ν = T ¯ μ ν + τ ¯ μ ν ,
τ ¯ i k b ¯ i k I = 0 .
Due to the nonsingularity of the (7) mentioned above, the last equation corresponds to Condition
τ ¯ i k = 0 .
Of the 10 linearized Einstein Equations (10), only six (for μ = i , ν = k ) are equations for the first correction y ( 1 ) a to the embedding function and hence for the first correction g ( 1 ) μ ν to the metric, as follows from (1). The remaining four equations define the components of τ ¯ 0 μ ; therefore, they do not restrict the embedding function.
Only the transverse contributions to the background surface in y ( 1 ) a change the shape of the surface, so that only six components of y ( 1 ) a make physical sense. The others four components correspond to the longitudinal deformations of the surface, which are equivalent to diffeomorphisms. Therefore, in order to uniquely find all 10 components of y ( 1 ) a , it is necessary to impose some coordinate conditions.

3. Solving of the Linearized Equations

Assume that at least in the zero order the distribution of ordinary matter is static, i.e., 0 T ¯ μ ν = 0 . We look for solutions of linearized Equation (10) that correspond to the static addition to the metric, i.e., 0 g ( 1 ) μ ν = 0 . This corresponds to the fact that, in the main order, the distribution of embedding matter is also static, i.e., 0 τ ¯ μ ν = 0 .
Instead of solving only 6 out of 10 Equations (10) (which can be performed and gives the same result), it is convenient to solve all 10 equations, treating the quantities of τ ¯ 0 μ as independent variables describing embedding matter. Then, the linearized Regge–Teitelboim equations are simply reduced to linearized Einstein Equation (10). We assume that ordinary matter is dustlike (as can be performed for the matter in galaxies) and is at rest in zero order; then,
T ¯ μ ν = ρ ¯ ( x i ) δ μ 0 δ ν 0 .
Equation (10) can be considered to be an equation not directly for y ( 1 ) a , but as an equation for the correction to metric g ( 1 ) μ ν , from which y ( 1 ) a is then found by solving embedding Equation (1).
The equations for g ( 1 ) μ ν take a simple form:
Δ g ( 1 ) 00 = ( ρ ¯ ( x i ) + τ ¯ 00 ) , Δ g ( 1 ) j k = ( ρ ¯ ( x i ) + τ ¯ 00 ) δ j k ,
Δ g ( 1 ) 0 k = 2 τ ¯ 0 k
(where Δ is the Laplacian operator) if we choose harmonic coordinates and assume that g ( 1 ) μ ν does not depend on time. In addition, we assume that g ( 1 ) μ ν is decreasing at spatial infinity. The conditions for the harmonic coordinates are reduced to only one additional equation:
k g ( 1 ) 0 k = 0 .
Introducing gravitational potential φ ( x i ) (we did not include the gravitational constant ϰ in it for convenience) as a solution of Poisson equation
Δ φ ( x i ) = 1 2 ( ρ ¯ ( x i ) + τ ¯ 00 )
and taking into account that the solutions are decreasing at spatial infinity, we obtain
g ( 1 ) 00 = 2 φ ( x i ) , g ( 1 ) j k = 2 φ ( x i ) δ j k .
Now, we find correction y ( 1 ) a to the embedding function that corresponds to an arbitrary correction g ( 1 ) μ ν to the metric. From Formula (1), in the first order by ϰ , we have:
g ( 1 ) μ ν = ( μ y ¯ a ) ( ν y ( 1 ) a ) + ( ν y ¯ a ) ( μ y ( 1 ) a ) .
Let us introduce a parameterization for the correction to the embedding function as the sum of a tangent and orthogonal to the background components:
y ( 1 ) a = y ( 1 ) γ γ y ¯ a + y ( 1 ) a .
The structure of Background Embedding (5) leads to relation y ( 1 ) 0 = y ( 1 ) 0 . Substituting Representation (20) in (19) yields a system of differential equations for y ( 1 ) γ and y ( 1 ) a :
η γ ν μ y ( 1 ) γ + η γ μ ν y ( 1 ) γ 2 y ( 1 ) I b ¯ μ ν I = g ( 1 ) μ ν .
Writing out its components, we obtain a system of equations:
2 0 y ( 1 ) 0 = g ( 1 ) 00 , k y ( 1 ) 0 + 0 y ( 1 ) k = g ( 1 ) 0 k , i y ( 1 ) j + j y ( 1 ) i 2 y ( 1 ) I b ¯ i j I = g ( 1 ) i j .
This can be solved by integrating the equations over time and taking into account the static g ( 1 ) μ ν . Denoting the integration constants as y ˜ a ( x i ) , we write the answer as follows:
y ( 1 ) 0 = 1 2 g ( 1 ) 00 x 0 + y ˜ 0 ( x i ) = φ ( x i ) x 0 + y ˜ 0 ( x i ) , y ( 1 ) k = g ( 1 ) 0 k + k y ( 1 ) 0 d x 0 + y ˜ k ( x i ) , y ( 1 ) I = 1 2 α ¯ I k j g ( 1 ) k j + 2 k y ( 1 ) j = α ¯ I k j φ ( x i ) δ k j + k y ( 1 ) j ,
where (18) was also used.
Substituting these expressions into (20), we lastly obtain
y ( 1 ) a = δ 0 a φ ( x i ) x 0 + y ˜ 0 ( x i ) + + δ I a g ( 1 ) 0 k + k y ( 1 ) 0 d x 0 + y ˜ k ( x i ) k y ¯ I + α ¯ I k j φ ( x i ) δ k j + k y ( 1 ) j .
This expression gives the solution of the linearized Regge–Teitelboim equations. It is parameterized by the gravitational potential φ ( x i ) that, via Equation (17), is associated with the zero-order τ ¯ 00 of the density of embedding matter, and via the components of metric g ( 1 ) 0 k that, via Equation (15), are related to the zero-order τ ¯ 0 k of the flux density of this matter. Due to (15) and (16), the components of τ ¯ 0 k obey condition k τ ¯ 0 k = 0 (it has the meaning of the embedding matter conservation condition in the case under consideration), so only two of these three components can be considered as independent.
It is interesting to compare this result with the one obtained in [19] in the framework of the nonrelativistic approximation. For this purpose, let us write down the time derivative of Solution (24):
0 y ( 1 ) 0 = φ ( x i ) , 0 y ( 1 ) I = g ( 1 ) 0 k + k y ( 1 ) 0 k y ¯ I + α ¯ I i j i g ( 1 ) 0 j + j y ( 1 ) 0 .
The second equation reproduces Equation (62) from [19], taking into account the assumption g ( 1 ) 0 j = 0 there. However, the first equation differs from Equation (61) of [19] through the absence of a quadratic contribution. Of course, this is not surprising, since the solution of (24) is obtained in the linear approximation and hence with the rejection of the quadratic contributions. However, when analyzing the nonrelativistic limit in [19], it was taken into account that for the large intervals of time x 0 (when time x 0 is replaced by nonrelativistic time t = x 0 / c , where c is the speed of light) quadratic corrections can give a comparable contribution to the linear ones.

4. Accounting for Equations in the Next Order

The solution obtained above in the first order of perturbation theory depends on the choice of the distribution of embedding matter in zero order characterized by τ ¯ 0 μ . Let us find restrictions on this distribution, assuming that both τ μ ν characterizing embedding matter and the energy–momentum tensor of ordinary matter T μ ν are static not only in the zero approximation, but also exactly. In particular, it follows that the metric is also exactly static. When considering a galaxy, this means that the configuration of both ordinary and embedding matter at some point passed into a stationary state. We investigate just this particular case.
Let us write down Equation (4) in the next order in ϰ , compared to already studied Equation (11), and use it:
τ ( 1 ) j k b ¯ j k a + τ ¯ 00 b ( 1 ) 00 a + 2 τ ¯ 0 k b ( 1 ) 0 k a = 0 .
Due to (6) and the properties of Background Embedding (5), we have b ( 1 ) 0 μ a = Π ¯ b a 0 μ y ( 1 ) b . This allows for us to express τ ( 1 ) j k from (26) in the form:
τ ( 1 ) j k = τ ¯ 00 0 0 y ( 1 ) I + 2 τ ¯ i 0 i 0 y ( 1 ) I α ¯ I j k .
Therefore, the assumed time independence of τ μ ν leads to additional restrictions on the quantities in the right-hand side. The time derivative of Equation (27) with (25) and condition 0 τ ( 1 ) j k = 0 provide the equation
τ ¯ i 0 δ i j δ l k + δ l j δ i k l φ + 2 i α ¯ I l m l m φ α ¯ I j k = 0 .
Then, we can also write down Equation (3) in the next order, and it provides the next correction g ( 2 ) μ ν to the metric. However, it is possible to write down a consequence of this equation in which the determined correction does not enter. This is the corresponding order of equation D μ τ μ ν = 0 , reflecting the conservation (with covariant corrections) of embedding matter. In the lowest order, this equation has the form μ τ ¯ μ ν = 0 and is satisfied given (12), (15) and (16). In the next order, it looks like
μ τ ( 1 ) μ ν + Γ ( 1 ) μ γ μ τ ¯ γ ν + Γ ( 1 ) μ γ ν τ ¯ μ γ = 0 .
Taking ν = k in this equation, and using Equation (27) and the assumed time independence of τ μ ν in all orders, from condition 0 τ ( 1 ) 0 k = 0 we obtain an equation:
j τ ¯ 00 0 0 y ( 1 ) I + 2 τ ¯ i 0 i 0 y ( 1 ) I α ¯ I j k τ ¯ 0 k η μ α ( α y ¯ a ) μ 0 y ( 1 ) a τ ¯ 00 ( k y ¯ a ) 0 0 y ( 1 ) a 2 τ ¯ j 0 ( k y ¯ a ) j 0 y ( 1 ) a = 0 ,
where (12) and an expression for connectivity are used (see, for example, [12])
Γ μ γ ν = g ν α ( α y a ) μ γ y a ,
from which it can be found that
Γ ( 1 ) μ 0 ν = η ν α ( α y ¯ a ) μ 0 y ( 1 ) a .
Using Formulas (25) in (30), we obtain an equation that, together with (28), restricts the choice of the background embedding y ¯ a and τ ¯ 0 μ , which parameterize the solution of the linear approximation found in the previous section. From the remaining Equation (29) at ν = 0 , variable τ ( 1 ) k 0 is not excluded, so this equation does not provide a new restrictions on the above quantities.
The background embedding of the form (5) is parameterized by three functions, since nine components of y ¯ I ( x i ) are imposed by six conditions that the metric of the surface described by this function is flat Euclidean. At first sight, the resulting number of equations for the background and τ ¯ 0 μ seems too large—six Equation (28) and three Equation (30) for six unknowns (as noted after (24)), only two of the three components τ ¯ 0 k are independent). However, it is easy to see that a solution nevertheless exists, since the six Equations (28) can be satisfied by choosing three quantities τ ¯ 0 k as
τ ¯ 0 k = 0 .
We study the solutions corresponding just to such a choice. Physically, it means that, in the zero approximation, embedding matter is at rest. Using Background (5) leads to the fact that time independence τ ¯ μ ν actually leads to (33) (otherwise, there are fewer unknowns than equations), i.e., to Embedding Matter at rest in the zero approximation. In the analysis of the nonrelativistic limit in [19], such a background was found when discussing the nonrelativistic character of motion of embedding matter. Taking into account (15) and assuming decreasing corrections to the metric at 3D infinity, it follows from (33) that g ( 1 ) 0 k = 0 .
Given that τ ¯ 0 k = 0 , the remaining Equations (30) are simplified. Using (25), we can now rewrite them as
j ( 2 Δ φ ρ ¯ ) α ¯ I j k α ¯ I l m l m φ ( 2 Δ φ ρ ¯ ) k φ = 0 ,
where τ ¯ 00 was expressed from (17).
These are three equations for four quantities, namely, three that parameterize the background embedding function y ¯ I ( x i ) and one describing the density of embedding matter τ ¯ 00 . Thus, linearized Regge–Teitelboim equations, supplemented by the condition of exact statics of solutions, leave an arbitrariness in the choice of one function of coordinates x i . The physical meaning of this arbitrariness lies in the presence of arbitrariness in the general case when specifying the distribution of embedding matter at the initial moment of time.

5. The Case of Spherical Symmetry

Let us consider the situation when ordinary matter is distributed in a spherically symmetric way. Then, the solution can also be sought in the class of spherically symmetric functions. In this case, both gravitational potential φ and the distribution density of ordinary matter ρ ¯ depend only on the radial coordinate r = x i x i .
The question of how to formulate a condition of any kind of symmetry for an embedding function is rather nontrivial; see, for example, [20]. Since the background embedding function y ¯ I ( x i ) enters Equation (34) only as a quantity
A l m , j k = α ¯ I l m α ¯ I j k ,
it is possible to directly impose symmetry constraints on this. Taking into account all index permutation symmetries, in the case of spherical symmetry, quantity A can be written in the form
A l m , j k = f 1 ( r ) δ l m δ j k + 1 2 f 2 ( r ) δ l j δ m k + δ l k δ m j + 1 2 r 2 f 3 ( r ) δ l m x j x k + δ j k x l x m +   + 1 4 r 2 f 4 ( r ) δ l j x m x k + δ l k x m x j + δ j m x k x l + δ m k x j x l + 1 r 4 f 5 ( r ) x l x m x j x k .
Here, functions f 1 f 5 are characteristics of the particular background embedding used. If we substitute this representation into (34), a common multiplier x k appears, and the three equations turn into one scalar equation. Introducing notations for different combinations of functions f 1 f 5 , the resulting equation can be written in the following form:
φ + 2 r φ 1 2 ρ ¯ φ F 1 + 1 r φ F 2 + 2 r φ + 2 r φ 1 2 ρ ¯ φ F 3 + 1 r φ F 4 = 0 ,
where
F 1 = f 1 + f 2 + f 3 + f 4 + f 5 ,
F 2 = 2 f 1 + f 3 ,
F 3 = f 2 + 1 2 f 3 + f 4 + f 5 ,
F 4 = f 2 + f 3 r 2 2 .
Equation (37) allows for us to determine gravitational potential φ ( r ) if the distribution of ordinary matter ρ ¯ ( r ) and the background embedding function defining f 1 f 5 are known. Since (37) is a nonlinear differential equation, in the general case, its solution can only be sought numerically. First of all, it is necessary to find all the functions f 1 f 5 included in the equation using the given embedding function, which requires first finding the value of α ¯ I l m .
Finding it directly with Formula (8) is a nontrivial task for specific embeddings. Let us find a simpler way that allows for us to find functions f 1 f 5 included in the expansion of (36) without explicit knowledge of α ¯ I l m . When writing this expansion, we took into account three types of symmetry by permutations of indices following from the definition of (35): inside the first pair, inside the second, and when replacing the first pair with the second. As a consequence, we have five independent coefficients f 1 f 5 . However, in fact, there are fewer independent variables parameterizing A. To show this, we define quantity B:
B j k , i s = b ¯ j k I b ¯ i s I = ( j k y ¯ I ) i s y ¯ I .
It is easy to see that this quantity is symmetric through the permutation of any two indices due to the orthogonality of the first and second derivatives of y ¯ I :
B j k , i s = ( j k y ¯ I ) i s y ¯ I = ( i j k y ¯ I ) s y ¯ I = ( i k y ¯ I ) j s y ¯ I = B i k , j s .
Hence, in the spherically symmetric case, it can be parameterized not by five but only by three radial functions:
B j k , i s = 1 3 g 1 ( r ) δ j k δ i s + δ j s δ i k + δ j i δ s k + 1 6 r 2 g 2 ( r ) ( δ j k x i x s + δ j i x k x s + + δ j s x i x k + δ i k x j x s + δ i s x j x k + δ s k x j x i ) + 1 r 4 g 3 ( r ) x j x k x i x s .
There is a connection between A and B that follows from (8):
A l m , j k B j k , i s = 1 2 δ i l δ s m + δ s l δ i m .
It can be used to express five functions f 1 f 5 through three functions g 1 g 3 :
f 1 = 24 g 1 2 + ( 24 g 2 36 g 3 ) g 1 + 3 g 2 2 4 ( 20 g 1 2 + ( 20 g 2 + 24 g 3 ) g 1 g 2 2 ) g 1 , f 2 = 3 2 g 1 , f 3 = ( 12 g 2 + 36 g 3 ) g 1 3 g 2 2 2 ( 20 g 1 2 + ( 20 g 2 + 24 g 3 ) g 1 g 2 2 ) g 1 , f 4 = 3 g 2 g 1 ( 2 g 1 + g 2 ) , f 5 = 360 g 1 2 g 3 + ( 138 g 2 2 + 108 g 2 g 3 ) g 1 3 g 2 3 4 ( 20 g 1 2 + ( 20 g 2 + 24 g 3 ) g 1 g 2 2 ) g 1 ( 2 g 1 + g 2 ) .
As a result, to find the functions f 1 f 5 corresponding to a specific background embedding, we first need to calculate the second fundamental form b ¯ i j I = i j y ¯ I . Then, we have to determine quantity B j k , i s according to (42). Comparing the result with the (44) expansion, we can find functions g 1 g 3 , and then find functions f 1 f 5 by (46).

6. Explicit Unfolded Spherically Symmetric Embedding

Let us present an embedding y ¯ I ( x i ) of the flat Euclidean 3D metric into 9D Euclidean space that would be spherically symmetric in the sense discussed in [20] and would also be unfolded in terms of [18]. The latter means that, for such a surface, the second fundamental form (7) is nonsingular as a 6 × 6 matrix; see the text after (7).
We consider basis λ i j A in the space of symmetric traceless matrices, numbered by index A, which runs through values 4 , , 8 . We subject this basis to the orthonormality condition in the space of such matrices:
λ i j A λ l m A = 1 2 δ i l δ j m + δ i m δ j l 1 3 δ i j δ l m , λ i j A λ i j B = δ A B .
We take the embedding function y ¯ I ( x i ) with components
y ¯ k ( x i ) = f ( r ) x k , y ¯ A ( x i ) = g ( r ) λ l m A x l x m , y ¯ 9 ( x i ) = h ( r ) ,
it is parameterized by three radial functions. One can check that this embedding function describes a surface with spherical symmetry. This means that there is such a representation (generally speaking, reducible) of the group of 3D rotations under the action of which in the space y ¯ I the surface transforms to itself (see [20] for details). The three independent functions f , g , h are related to its irreducible parts, and the reducible representation as a whole can be described as “3+5+1”.
Let us select, among all the embeddings of type (48), those that correspond to a flat metric. Calculating for (48) the components of the induced metric, we obtain the equations
f 2 + 2 r 2 g 2 = 1 , r f 2 + 2 3 r 2 g 2 + h 2 = 1 .
Expressing from these equations functions g and h through f we obtain the desired spherically symmetric embedding of the flat Euclidean 3D metric into 9D Euclidean space. It is parameterized by choosing one function of radius f ( r ) . The fact that such an embedding is unfolded by a non-special choice of this function can be checked directly by computing the second fundamental form i j y ¯ I (7) and making sure that it is nonsingular.
Calculating this quantity and Expression (42) constructed from it lead to the following expressions for the coefficients g 1 g 3 of its expansion (44):
g 1 = 6 g 2 ,
g 2 = 6 f 2 + 2 r 2 g 2 + 4 r g g ,
g 3 = r 2 h r 2 + r 2 f 2 + 2 r f f 3 f 2 + 2 3 r 4 g 2 + 4 r 3 g g + 4 r 2 g g 2 r 2 g 2 4 r g g ,
where the functions g and h must be expressed through a single arbitrary function f ( r ) according to (49).

7. Inverse Problem

Equation (37) obtained above in Section 5 allows for searching for the gravitational potential φ ( r ) if the distribution of ordinary matter ρ ¯ ( r ) and the background embedding function y ¯ I ( x i ) are known. It is easy to see that, among the solutions of this equation, there is an Einsteinian solution for which
φ + 2 r φ 1 2 ρ ¯ = 0 .
However, there are other solutions too. In the galactic halo region where the contribution of ρ ¯ can be neglected, we have a Newtonian behavior of the gravitational potential φ 1 / r for the Einsteinian solution. Such a behavior does not agree well with the data of the rotation curves of galaxies. To find other solutions in this area, it is sufficient to know background embedding y ¯ I ( x i ) .
The background embedding constructed in the previous section has ambiguity in the choice of one function f ( r ) . There is no highlighted natural way to choose this function. As a consequence, there is also no possibility to uniquely find the gravitational potential φ ( r ) as a result of solving Equation (37), since changing f ( r ) also changes φ ( r ) . Therefore, we focus on investigating the inverse problem by posing the question: can we find such a background embedding y ¯ I ( x i ) (i.e., such a function f ( r ) ) so that the solution of Equation (37) corresponds to the gravitational potential in good agreement with galactic halo observations?
As a good fitting of observational data, we take the Burkert profile [21] of the dark matter density; see, for example, review [22]:
ρ B ( r ) = ρ 0 r 0 3 ( r + r 0 ) ( r 2 + r 0 2 ) ,
where r 0 is the core radius, ρ 0 is the density at the center of the galaxy. In order to find a behavior of the gravitational potential that agrees well with the observations, one needs to solve the equation
φ + 2 r φ = 1 2 ( ρ ¯ + ρ B ) .
The solution of this equation in the galactic halo has the form:
φ ( r ) = M r + ρ 0 2 1 + r 0 r arctan r r 0 ln 1 + r r 0 + 1 2 1 r 0 r ln 1 + r 2 r 0 2 ,
where M is the mass of ordinary matter in the galaxy. To avoid introducing this quantity as an additional parameter, we further assume that the proportional contribution to M to the gravitational potential in the halo region can be neglected compared with the contribution of dark matter. The trivial substitution of Function (56) with M = 0 into Equation (37) leads to a very complicated (so we do not show it here) third-order nonlinear differential equation for the unknown function f ( r ) parameterizing the background embedding. This equation fails to be solved analytically. We numerically solved the Cauchy problem for Equation obtained, choosing different values of the initial data at the value of the radius r = r 0 / 2 , which is usually about the visual size of the galaxy. We solve Equation (37) with a zero contribution of ordinary matter, so it is applicable only in the region of the galactic halo.
A typical view of the obtained solutions is shown in Figure 1. The existence of such solutions shows the possibility to choose parameterizing embedding (48) function f ( r ) in a special way. For this choice, equations of motion in the framework of linearization of the embedding theory in this background lead to the same gravitational potential as the dark matter with Burkert profile (54) in the framework of ordinary GR. With a significant change in the chosen values of the initial data, the algorithm used for the numerical solution of the nonlinear differential equation no longer made it possible to extend the solution into the region r 0.95 r 0 . This can be explained by the fact that in the differential equation being solved, the coefficient at the major derivative (which itself is a compound quantity that also depends on the minor derivatives) may become too small, which prevents the correct behavior of the algorithm.

8. Conclusions

We considered the problem of constructing a perturbation theory for Regge-Teitelboim Equation (2) for some background. The background was chosen to be (5), which corresponded to a 4D surface that is a direct product of the timelike straight line on the 9D unfolded embedding of the euclidean 3D metric. This formulation of the problem is a special case of the weak gravitational field approximation in the embedding theory, in which the corrections to the embedding function are linear to the corrections to the metric; see (23). The analysis is simplified if the Regge–Teitelboim equations are rewritten as a set of Einstein Equation (3) with the contribution of fictitious embedding matter and Equation (4) describing the behavior of this matter.
First, the Einstein equations were linearized in the usual way and solved using harmonic coordinates. The result depended on the energy-momentum tensor of embedding matter specified arbitrarily in the lower order. Further, we restricted ourselves to the class of solutions for which both ordinary matter and embedding matter are in a stationary state, which corresponded to the description of an already fully formed galaxy. Taking this condition into account in the second order of perturbation theory leads to additional restrictions on the solutions found in the linear one. As a result, Equation (37) arises in the case of spherical symmetry, which allows for finding the gravitational potential φ ( r ) if the background embedding is given. The nonlinear form of this equation is explained by the fact that it arises as a result of accounting for the static character of the metric in the second (following the linear) order.
It is possible to find a class of spherically symmetric surfaces (48) and (49) that have all necessary properties for their use as a background embedding. These embeddings are parameterized by a single function f ( r ) . This function could be chosen, so that the solution of Equation (37) led to a gravitational potential that corresponded well (if we neglected deviations from spherical symmetry for real galaxies) to the observed distribution of dark matter in the galactic halo.
In our analysis we used the following assumptions. We assumed that ordinary matter has, on average, static spherically symmetric distribution. As a consequence, the energy–momentum tensor of ordinary matter T μ ν has exactly spherically symmetric and static character (we use this fact in the first two orders of perturbation theory: for T ¯ μ ν and T ( 1 ) μ ν ). We further assumed that the metric had an exactly spherically symmetric and static character, and we used this for corrections to the metric g ( 1 ) μ ν and g ( 2 ) μ ν . Given the restrictions on T μ ν , imposing such an assumptions on the metric is the same as imposing them on the energy-momentum tensor of embedding matter in the first two orders: for τ ¯ μ ν and τ ( 1 ) μ ν . As a background embedding function y ¯ a ( x μ ) that corresponded to the flat metric η μ ν , we took the Product (5) of a timelike line on 9D unfolded and the spherically-symmetric embedding (48) of the euclidean 3D metric. The assumptions were internally consistent, but for real galaxies, they are satisfied only approximately. The strongest ones are the deviations from the spherical symmetry.
We discussed the behavior of the gravitational potential φ ( r ) only in the galactic halo. In order to study the behavior of φ ( r ) inside the galaxy where the influence of ordinary matter cannot be neglected, one must additionally find an explicit form of its distribution, which is beyond the scope of this paper. At the same time, deviations of embedding gravity predictions from GR arise at scales where energy–momentum tensor τ μ ν of embedding matter becomes significant compared to the energy–momentum tensor T μ ν of ordinary matter. We assumed that, at the scale of the galaxy, these two quantities were comparable, but at smaller distances, the influence of the embedding matter could be neglected, and the laws of GR would start to be satisfied with good accuracy.

Author Contributions

Conceptualization, S.K. and S.P.; investigation, S.K., M.I., S.M. and S.P.; writing, S.K., M.I., S.M. and S.P. All authors have read and agreed to the published version of the manuscript.

Funding

The work of S.S. Kuptsov was supported by Leonhard Euler International Mathematical Institute grant no. 075-15-2019-1620. The work of S.A. Paston was supported by RFBR grant no. 20-01-00081.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Regge, T.; Teitelboim, C. General relativity à la string: A progress report. In Proceedings of the First Marcel Grossmann Meeting, Trieste, Italy, 7–12 July 1975; pp. 77–88. [Google Scholar]
  2. Deser, S.; Pirani, F.A.E.; Robinson, D.C. New embedding model of general relativity. Phys. Rev. D 1976, 14, 3301–3303. [Google Scholar]
  3. Tapia, V. Gravitation a la string. Class. Quant. Grav. 1989, 6, L49. [Google Scholar]
  4. Maia, M.D. On the integrability conditions for extended objects. Class. Quant. Grav. 1989, 6, 173–183. [Google Scholar]
  5. Franke, V.A.; Tapia, V. The ADM Lagrangian in extrinsic gravity. Nuovo Cimento B 1992, 107, 611. [Google Scholar]
  6. Estabrook, F.B.; Robinson, R.S.; Wahlquist, H.R. Constraint-free theories of gravitation. Class. Quant. Grav. 1999, 16, 911–918. [Google Scholar]
  7. Karasik, D.; Davidson, A. Geodetic Brane Gravity. Phys. Rev. D 2003, 67, 064012. [Google Scholar]
  8. Cordero, R.; Molgado, A.; Rojas, E. Ostrogradski approach for the Regge-Teitelboim type cosmology. Phys. Rev. D 2009, 79, 024024. [Google Scholar]
  9. Paston, S.A. Gravity as a field theory in flat space-time. Theor. Math. Phys. 2011, 169, 1611–1619. [Google Scholar]
  10. Faddeev, L.D. New dynamical variables in Einstein’s theory of gravity. Theor. Math. Phys. 2011, 166, 279–290. [Google Scholar]
  11. Paston, S.A.; Sheykin, A.A. From the Embedding Theory to General Relativity in a result of inflation. Int. J. Mod. Phys. D 2012, 21, 1250043. [Google Scholar]
  12. Paston, S.A.; Franke, V.A. Canonical formulation of the embedded theory of gravity equivalent to Einstein’s general relativity. Theor. Math. Phys. 2007, 153, 1582–1596. [Google Scholar]
  13. Undagoitia, T.M.; Rauch, L. Dark matter direct-detection experiments. J. Phys. G Nucl. Part. Phys. 2015, 43, 013001. [Google Scholar]
  14. Gaskins, J.M. A review of indirect searches for particle dark matter. Contemp. Phys. 2016, 57, 496–525. [Google Scholar]
  15. Pavsic, M. On The Quantization Of Gravity By Embedding Space-Time In A Higher Dimensional Space. Class. Quant. Grav. 1985, 2, 869. [Google Scholar]
  16. Chamseddine, A.H.; Mukhanov, V. Mimetic dark matter. J. High Energy Phys. 2013, 2013, 135. [Google Scholar]
  17. Golovnev, A. On the recently proposed mimetic Dark Matter. Phys. Lett. B 2014, 728, 39–40. [Google Scholar]
  18. Paston, S.A.; Zaitseva, T.I. Nontrivial isometric embeddings for flat spaces. Universe 2021, 7, 477. [Google Scholar]
  19. Paston, S.A. Non-relativistic limit of embedding gravity as General Relativity with dark matter. Universe 2020, 6, 163. [Google Scholar]
  20. Paston, S.A.; Sheykin, A.A. Embeddings for Schwarzschild metric: Classification and new results. Class. Quant. Grav. 2012, 29, 095022. [Google Scholar]
  21. Burkert, A. The Structure of Dark Matter Halos in Dwarf Galaxies. Astrophys. J. 1995, 447, L25. [Google Scholar]
  22. Di Paolo, C.; Salucci, P. Fundamental properties of the dark and the luminous matter from Low Surface Brightness discs. arXiv 2020, arXiv:2005.03520. [Google Scholar]
Figure 1. Numerically constructed solutions with different initial data: (a) f = 0 , f = 0.025 , f = 0.1 ; (b) f = 0.1 , f = 0.02 , f = 0.1 .
Figure 1. Numerically constructed solutions with different initial data: (a) f = 0 , f = 0.025 , f = 0.1 ; (b) f = 0.1 , f = 0.02 , f = 0.1 .
Universe 08 00635 g001
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Kuptsov, S.; Ioffe, M.; Manida, S.; Paston, S. Weak Field Limit for Embedding Gravity. Universe 2022, 8, 635. https://doi.org/10.3390/universe8120635

AMA Style

Kuptsov S, Ioffe M, Manida S, Paston S. Weak Field Limit for Embedding Gravity. Universe. 2022; 8(12):635. https://doi.org/10.3390/universe8120635

Chicago/Turabian Style

Kuptsov, Stanislav, Mikhail Ioffe, Sergey Manida, and Sergey Paston. 2022. "Weak Field Limit for Embedding Gravity" Universe 8, no. 12: 635. https://doi.org/10.3390/universe8120635

APA Style

Kuptsov, S., Ioffe, M., Manida, S., & Paston, S. (2022). Weak Field Limit for Embedding Gravity. Universe, 8(12), 635. https://doi.org/10.3390/universe8120635

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop