The Schwarzschild-de Sitter Metric of Nonlocal $\sqrt{dS}$ Gravity

It is already known that a simple nonlocal de Sitter gravity model, which we denote as $\sqrt{dS}$ gravity, contains an exact vacuum cosmological solution which mimics dark energy and dark matter and is in very good agreement with the standard model of cosmology. This success of $\sqrt{dS}$ gravity motivated us to investigate how it works at lower than cosmic scale -- galactic and the solar system. This paper contains our investigation of the corresponding Schwarzschild-de Sitter metric of the $\sqrt{dS}$ gravity model. To get exact solution, it is necessary to solve the corresponding nonlinear differential equation, what is a very complicated and difficult problem. What we obtained is a solution of linearized equation, which is related to space metric far from the massive body, where gravitational field is weak. The obtained approximate solution is of particular interest for examining the possible role of non-local de Sitter gravity $\sqrt{dS}$ in describing the effects in galactic dynamics that are usually attributed to dark matter. The solution has been tested on the Milky Way and the spiral galaxy M33 and is in good agreement with observational measurements.


Introduction
The general theory of relativity (GR) [1] is considered as one of the most successful and beautiful physical theories.It is worth mentioning the following main predictions and successful confirmations: deflection of light near the Sun, black holes, gravitational light redshift, lensing, and gravitational waves, see e.g.[2].
In the Standard Model of Cosmology (SMC) [3], also known as ΛCDM model, GR is adopted as theory of gravitation at all space-time scales, from the solar system to the galactic and cosmic ones.To describe galactic rotational curves and accelerated expansion of the universe by GR, it was supposed existence of dark matter (DM) and dark energy (DE), respectively.According to the current SMC, the universe matter/energy budget consists approximately of 68 % of dark energy, 27 % of dark matter and only 5 % of standard visible matter.However, despite many experimental and theoretical investigations (as a review, see e.g.[4]), the existence of DM and DE has not been proven, thus they are still hypothetical constituents of the dark side of the universe.
It should be also mentioned that GR suffers from the singularities -the black hole and Big Bang singularity [5].As we know, if theory contains a singularity it means its inapplicability when approaching to singularity and that a more general appropriate theory should be invented.Also, it should be mentioning that GR can not be consistently quantized [6].As we know, other physical theories have their own domain of validity, usually limited by the space-time scale and some parameters, or by complexity of the system.In this sense, GR should not be an exception and serve as a theory of gravity from the Planck scale to the universe as a whole [7].
Based on all the aforementioned shortcomings, it can be concluded that general relativity is not a final theory of gravitation and that there is a sense to look for a more general theory than GR [8].In principle, there is a huge number of possibilities to extend the Einstein-Hilbert (EH) action and for now there is no rule on how to choose the right path [9,10].Hence, in practice there are many phenomenological approaches and the most elaborated is f (R) theory [11], where in the EH action the scalar curvature R is replaced by some function f (R).Among other interesting approaches is the nonlocal one [39,13].In the nonlocal gravity models, besides R in the EH action there is a nonlocal term with some invariants usually composed of R and 2, where 2 is the d'Alembert-Beltrami operator.
Depending on how the 2 is built into the nonlocal term, there are mainly two typical examples of nonlocal gravity models: (i) non-polynomial analytic expansion of 2, see e.g.[14,15,17,16,18,19,20], i.e.F(2) = ∞ n=1 f n 2 n (see various examples [21,22,23,24,25,26]), and (ii) a polynomial of 2 −1 , see e.g.[39,27,28,29,30].The motivation for using a local operator of the form (i) is found in string theory -ordinary and p-adic one [31].It is obvious that in the case (i) dynamics depends not only on the first and second space-time derivatives but also on the all higher ones.Nonlocal operator 2 −1 in (ii) has its origin in (one-loop) quantum corrections to some classical field Lagrangians and is used in investigation of the late time cosmic acceleration without dark energy [39].
In several papers, see [19,30] and references therein, we investigated the following nonlocal de Sitter gravity model ( √ dS gravity): where Λ is the cosmological constant and nonlocal operator F(2) has the following general form: This nonlocal model is unique compared to other non-local models and its properties will be described in the next section.
It is worth mentioning that (1) applied to the homogeneous and isotropic universe gives several exact cosmological solutions [19,30].One of them is a(t) = At The fact that the model works well on a cosmological scale was motivation to test it on galactic and planetary systems.To this end, it is necessary to get the corresponding Schwarzschild-de Sitter metric.In the paper [32] we presented an initial research with an approximate solution.In this paper, we provide a much wider and more detailed investigation with some new solutions.It also contains preliminary test in our galaxy Milky Way and the spiral galaxy M33 with very satisfactory agreement of obtained theoretical results and observational measurements.
This paper is organized as follows.In Section 2, nonlocal de Sitter gravity √ dS is introduced and the equations of motion for the gravitational field are derived.Various aspects of the Schwarzschild-de Sitter metric are presented in Section 3, which in particular contains the solutions and their comparison with observations of the rotation curves of spiral galaxies.Some discussion and presentation of the main results is contained in the last section.

√ dS Nonlocal Gravity
Our nonlocal gravity model is given by its action (1).It can be rewritten in the compact form where F (2) is with the general form of the d'Alembert-Beltrami operator If F (2) = 1, i.e.F(2) = 0, then (3) becomes local de Sitter gravity with action Note that action (3) can be easily obtained from (6) by embedding operator (4) inside the product It should be noted that the degree of R − 2Λ remains unchanged when we go from local (6) to non-local action (3), as well as that F (2) operator is dimensionless.It is also worth noting that above local and nonlocal action has the same discrete symmetry, i.e. remains unchanged under transformation √ R − 2Λ → − √ R − 2Λ.In this paper, we will not consider the extension of action (1) with the matter sector, since we are looking for the Schwarzschild-de Sitter metric outside the spherically symmetric massive body.

Equations of Motion
To obtain equations of motion for √ dS gravity given by action (1), it is useful to start from more general nonlocal de Sitter model where P (R) and Q(R) are some differentiable functions of the Ricci scalar R. Variation of ( 7) with respect to g µν yields the corresponding equations of motion (EoM) derived in [33], see also [30].
According to [33], the EoM for nonlocal de Sitter gravity model ( 7) are as follows: where G µν = R µν − 1 2 Rg µν is the Einstein tensor, R µν is the Ricci tensor and where P ′ (R) and Q ′ (R) denote the derivative of P (R) and Q(R) with respect to R.
Comparing the above equations of motion with respect to their local Einstein counterpart G µν + Λg µν = 0 they look very complex and finding some exact solutions may be a hard problem.
Since we are interested in the EoM of √ dS, we have to take To this end, let us first consider the case Q(R) = P (R).Consequently, equations ( 8)-( 11) reduce to According to our experience, the above equations of motion ( 12)-( 14) can be significantly simplified and easily solved if there exists a metric tensor g µν such that for the corresponding d'Alembert-Beltrami operator 2 the following equations (eigenvalue problem) is satisfied: where q is a parameter of the same dimensionality as 2. Applying (15) to equations ( 12)-( 14), we have Ω µν = F ′ (q)S µν P, P , Now, let us take If the nonlocal operator satisfies then equations of motion (19) are also satisfied.According to the above consideration, the main problem is to solve equation R − 2Λ for an appropriate metric tensor g µν .This my be a hard problem, and it is the case with the corresponding Schwarzschild-de Sitter metric in the nonlocal √ dS gravity.In the sequel, we will investigate the corresponding Schwarzschild-de Sitter metric around static spherically symmetric massive body.

The Schwarzschild-de Sitter Metric
We are going to explore the Schwarzschild-de Sitter space-time metric in the case of √ dS gravity given by its action (1).

General Consideration
First, we want to consider some general aspects of the corresponding Schwarzschildde Sitter space-time metric.To this end, we start from the usual expression for the Schwarzschild metric in the pseudo-Rimannian manifold Non-zero elements of metric tensor g µν (r) are: Non-zero components of the corresponding Christoffel symbol where ' denotes derivative with respect to radius r.Non-zero components of the Ricci tensor where as a result of spherical symmetry Equation that should be solved is where is the Laplace operator in spherical coordinate r.
Note that in eigenvalue problem (29) parameter q has the same dimension as operator 2, i.e. dimension is L −2 .Since the cosmological constant Λ also has dimension L −2 , it is useful to write q = ζΛ, where ζ is a dimensionless parameter.Note also that nonlocal operator F(2), defined in (2), satisfies conditions F(q) = −1 and F ′ (q) = 0, introduced in (20), if we take it as , where a

Solutions
Recall that in the local de Sitter case (6), with static spherically symmetric body of mass M , the Schwarzschild-de Sitter metric (21) is It makes sense to suppose that solution of equation ( 29) is of the form where α(r) and β(r) are some dimensionless functions.When ζ → 0 then nonlocal operator (31) tends to zero and consequently nonlocal de Sitter √ dS gravity model (1) becomes local.Hence, it must be that A(r) → A 0 (r) and B(r r) in scalar curvature R (28) and in operator 2 of equation ( 29), we obtain If we substitute expressions (34) and (35) in equation ( 29) then we will get a differential equation with two unknown functions: α(r √ ζΛ) and β(r √ ζΛ).It is obvious that we must have another equation which is a relation between functions α and β.Recall that in the local case holds B 0 (r) = 1 A 0 (r) .Hence, there is a sense to take also B(r in the nonlocal case and it yields In fact, it means that we take Employing (38) in ( 34) and ( 35), we get We are interested in finding function α(r) and the next step should be substitution of into equation (40).In that case, one gets an ordinary nonlinear differential equation of the fourth order.Because of nonlinearity, it is a very difficult task to find the corresponding exact solution.After many attempts, we did not succeed to find a reasonable exact solution and concluded that a much more sophisticated approach is required.In the sequel of this paper we will turn our attention to the corresponding linear differential equation.It means we will limit ourselves to studying the Schwarzshild-de Sitter metric in weak gravity field approximation.Practically, it is like considering gravity field far from a massive body (see Figure 1) so that the d'Alembertian 2 can be replaced by the Laplacian △ in equation (40).In such case we will take A(r) ≈ 1 in (40), that is what makes sense if the following is satisfied: Figure 1: We consider the Schwarzschild-de Sitter metric of nonlocal √ dS gravity at the distances far from a spherically symmetric massive body.
The next step is to replace scalar curvature R by R = 4Λ + 2α r 2 + 4α ′ r + α ′′ (see (39)) in ( 46) and we obtain that is a linear differential equation of the fourth order.Equation ( 47) has a general solution of the form There are four constants (C 1 , C 2 , C 3 , C 4 ) and we have chosen them so that the appropriate particular solution α(r) → 0 when ζ → 0 and that it has some physical meaning.
To exclude term with e √ qr in (48), we took C 4 = 0, since this exponential function increases indefinitely for very large values of r.

Case
In this case solution for α(r) is where δ is dimensionless parameter.Since integration constants C 1 , C 2 , C 3 are proportional to δ and C 4 = 0, by this way we reduced the number of parameters from 4 to 1. Altogether, we hove two free parameters (δ and ζ) which should be determined from measurements.
To see how α(r) behaves when ζ → 0, it is useful to expand exponential function e − √ qr into the Taylor series.We have As it follows from (52), we conclude that α(r) → 0 when ζ → 0.

The Rotation Curves of Spiral Galaxies
The rotation curves of spiral galaxies have been the subject of intensive research motivated by the need to determine the amount and distribution of dark matter comparing to visible matter, see e.g.[34,35,36,37,38,39] and references therein.
It is interesting to examine whether this √ dS gravitational model gives the possibility of describing the rotation curves of spiral galaxies.To this end, we should start with A(r) given by (53) and present the corresponding gravitational potential Φ(r), which is Note that here Φ(r) is the intensity of the gravitational potential.
The corresponding gravitational acceleration for potential (54) is Velocity of the rotation curve v(r) follows from equality We checked the validity of the obtained formula for the circular velocity (56) on two cases: the Milky Way galaxy and the spiral galaxy M33.The values of parameters δ and ζ in(56) are estimated by best fitting of measured data using the least-squares method.Since at large distances r, the velocity v weakly depends on mass variation, we employed only the mass of the black hole in the center of the galaxy.Namely, the central mass can be taken up to M ∼ 10 8 M ⊙ and there will be no significant changes in the calculated circular velocity v(r) at very large r.

Milky Way case.
The Milky Way rotation curve data have taken from recent paper [40], where the Keplerian decline in the rotation curve is detected.Measured data for distance r, velocity v and velocity error ∆v are obtained by Gaia telescope and they are presented in Table 1, see [40].In this table also are presented computed velocity v using (56) and relative error δv = |v−v| v .A pictorial comparison of measured and calculated velocities is presented in Figure 2.
Table 1: Milky Way rotation curve data, from [40] and this work.We have used data for the galaxy Messier 33, based on observations obtained at the Dominion Radio Astrophysical Observatory and presented in [41].We have compared measured rotation velocity v(r) [41] with the computed v(r) using formula (56).Measured and computed data are presented in Table 2 and illustrated in Figure 3.

Discussion and Concluding Remarks
This paper presents the results of our research regarding the Schwarzschild-de Sitter metric of the nonlocal √ dS gravity model (1).We found the Schwarzschild-de Sitter metric in the form of A(r) (53), what corresponds to the weak gravity approximation and the linearization of nonlinear differential equation (40).The obtained results were tested on the rotation curves of the Milky Way and the spiral galaxy M33.The calculated and measured values of circular velocities are in good agreement.
Some additional explanations should be given to some parts of these investigations.First, we need to clarify why the weak gravitational field approximation works well here.On the one hand, we derived the Schwarzschild-de Sitter metric away from the massive spherically symmetric body.And on the other hand, we applied the obtained formula for the circular motion of the test body to the circular velocities in spiral galaxies far from their centers where the black hole is located.Recall that the rotation curves were observed in the domain: 9.5 -26.5 kpc for the Milky Way galaxy [40] and 0.5 -23.5 kpc for the M33 galaxy [41].In the Lambda Cold Dark Matter model, it is assumed that dark matter plays an important role in the mentioned domains.However, there is no dark matter in our nonlocal model.The good agreement between observational measurements and theoretical predictions tells us that the role of dark matter can be played by the nonlocality in the presence of the cosmological constant Λ in the √ dS gravity model.
Regarding the applicability of the obtained formula for the circular velocity (56) at smaller distances, such as the solar system, the following should be noted.The circular velocity v(r) depends on three terms: (i) GM c 2 r , (ii) − Λr 2 3 and (iii) δ √ qr 2 qr 2 e − √ q r .The third term depends linearly on δ and with a fixed q(= ζΛ) its value can be controlled by choosing the appropriate value of δ.One can always take a small enough value of δ, e.g.δ < |8.2 × 10 −14 |, so that the first term has a dominant role, since the second term has an important meaning only at distances of the size of the visible universe.Therefore, the velocity formula (56) is also valid for the solar system.
The main new results presented in this article can be summarized as follows.
• In the approximation of the weak gravitational field, a fourth-order linear differential equation for the Schwarzschild-de Sitter metric was obtained (47).
• A particular solution of α(r) was found (49) such that it satisfies the necessary condition that it tends to zero when the nonlocality vanishes.
• The obtained theoretical formula for circular velocity (56) was tested on the rotation curves of two spiral galaxies: the Milky Way and M33.The agreement between the calculated and measured circular velocities is good, especially for the Milky Way, see Figures 2 and 3 and the corresponding tables.To our knowledge, this is the first good description of "the Keplerian decline in the Milky Way rotation curve" by some modified gravity model.
In summary, it can be said that the presented results in this paper are encouraging and deserve further research, especially taking into account the mass distribution in spiral galaxies using [42].Bearing in mind also previously obtained results on the evolution of the universe [19,30], where the effects that are usually attributed to dark energy and dark matter can be described by the nonlocality of the gravity model √ dS, we will continue with the further study of this model of nonlocal de Sitter gravity.

2 3
which mimics an interplay of dark matter (t ) and dark energy (e Λ 14 t 2 ) in very good agreement with the standard model of cosmology.There are also nonsingular bounce solutions in the flat, closed and open universe as well as singular and cyclic solutions.