The Schwarzschild-de Sitter Metric of Nonlocal dS Gravity

It is already known that a simple nonlocal de Sitter gravity model, which we denote as 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 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 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 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.

Keywords:

Subject: Physical Sciences - Astronomy and Astrophysics

1. 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 [13,39]. In the nonlocal gravity models, besides R in the EH action there is a nonlocal term with some invariants usually composed of R and □, where □ is the d’Alembert-Beltrami operator.

Depending on how the □ is built into the nonlocal term, there are mainly two typical examples of nonlocal gravity models: (i) non-polynomial analytic expansion of □, see e.g. [14,15,16,17,18,19,20], i.e.

F (□) = \sum_{n = 1}^{\infty} f_{n} □^{n}

(see various examples [21,22,23,24,25,26]), and (ii) a polynomial of

□^{- 1}

, see e.g. [27,28,29,30,39]. 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

□^{- 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 (

\sqrt{d S}

gravity):

\begin{matrix} S = \frac{1}{16 π G} \int d^{4} x \sqrt{- g} (R - 2 Λ + \sqrt{R - 2 Λ} F (□) \sqrt{R - 2 Λ}), \end{matrix}

(1)

where

Λ

is the cosmological constant and nonlocal operator

F (□)

has the following general form:

\begin{matrix} F (□) = \sum_{n = 1}^{+ \infty} (f_{n} □^{n} + f_{- n} □^{- n}) . \end{matrix}

(2)

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) = A t^{\frac{2}{3}} e^{\frac{Λ}{14} t^{2}}

, which mimics an interplay of dark matter (

t^{\frac{2}{3}}

) and dark energy (

e^{\frac{Λ}{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.

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

\sqrt{d S}

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.

2. $\sqrt{dS}$ Nonlocal Gravity

Our nonlocal gravity model is given by its action (1). It can be rewritten in the compact form

\begin{matrix} S = \frac{1}{16 π G} \int d^{4} x \sqrt{- g} \sqrt{R - 2 Λ} F (□) \sqrt{R - 2 Λ}, \end{matrix}

(3)

where

F (□)

\begin{matrix} F (□) = 1 + F (□) = 1 + \sum_{n = 1}^{+ \infty} (f_{n} □^{n} + f_{- n} □^{- n}) \end{matrix}

(4)

with the general form of the d’Alembert-Beltrami operator

\begin{matrix} □ = \nabla_{μ} \nabla^{μ} = \frac{1}{\sqrt{- g}} \partial_{μ} (\sqrt{- g} g^{μ ν} \partial_{ν}) . \end{matrix}

(5)

F (□) = 1

, i.e.

F (□) = 0

, then (3) becomes local de Sitter gravity with action

\begin{matrix} S_{0} = \frac{1}{16 π G} \int d^{4} x \sqrt{- g} \sqrt{R - 2 Λ} \sqrt{R - 2 Λ} = \frac{1}{16 π G} \int d^{4} x \sqrt{- g} (R - 2 Λ) . \end{matrix}

(6)

Note that action (3) can be easily obtained from (6) by embedding operator (4) inside the product

\sqrt{R - 2 Λ} \sqrt{R - 2 Λ} .

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 (□)

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

\sqrt{R - 2 Λ} \to - \sqrt{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.

2.1. Equations of Motion

To obtain equations of motion for

\sqrt{d S}

gravity given by action (1), it is useful to start from more general nonlocal de Sitter model

\begin{matrix} S = \frac{1}{16 π G} \int d^{4} x \sqrt{- g} (R - 2 Λ + P (R) F (□) Q (R)), \end{matrix}

(7)

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:

\begin{matrix} G_{μ ν} + Λ g_{μ ν} - \frac{1}{2} g_{μ ν} P (R) F (□) Q (R) + R_{μ ν} W - K_{μ ν} W + \frac{1}{2} Ω_{μ ν} = 0, \end{matrix}

(8)

where

G_{μ ν} = R_{μ ν} - \frac{1}{2} R g_{μ ν}

is the Einstein tensor,

R_{μ ν}

is the Ricci tensor and

\begin{matrix} W = P^{'} (R) F (□) Q (R) + Q^{'} (R) F (□) P (R), K_{μ ν} = \nabla_{μ} \nabla_{ν} - g_{μ ν} □, \end{matrix}

(9)

\begin{matrix} Ω_{μ ν} = \sum_{n = 1}^{+ \infty} f_{n} \sum_{ℓ = 0}^{n - 1} S_{μ ν} (□^{ℓ} P, □^{n - 1 - ℓ} Q) - \sum_{n = 1}^{+ \infty} f_{- n} \sum_{ℓ = 0}^{n - 1} S_{μ ν} (□^{- (ℓ + 1)} P, □^{- (n - ℓ)} Q), \end{matrix}

(10)

where

\begin{matrix} S_{μ ν} (A, B) = g_{μ ν} (\nabla^{α} A \nabla_{α} B + A □ B) - 2 \nabla_{μ} A \nabla_{ν} B . \end{matrix}

(11)

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

\sqrt{d S}

, we have to take

Q (R) = P (R) = \sqrt{R - 2 Λ}

. To this end, let us first consider the case

Q (R) = P (R)

. Consequently, equations (8)–(11) reduce to

\begin{matrix} G_{μ ν} + Λ g_{μ ν} - \frac{g_{μ ν}}{2} P (R) F (□) P (R) + R_{μ ν} W - K_{μ ν} W + \frac{1}{2} Ω_{μ ν} = 0, \end{matrix}

(12)

\begin{matrix} W = 2 P^{'} (R) F (□) P (R), K_{μ ν} = \nabla_{μ} \nabla_{ν} - g_{μ ν} □, \end{matrix}

(13)

\begin{matrix} Ω_{μ ν} = \sum_{n = 1}^{+ \infty} f_{n} \sum_{ℓ = 0}^{n - 1} S_{μ ν} (□^{ℓ} P, □^{n - 1 - ℓ} P) - \sum_{n = 1}^{+ \infty} f_{- n} \sum_{ℓ = 0}^{n - 1} S_{μ ν} (□^{- (ℓ + 1)} P, □^{- (n - ℓ)} P) . \end{matrix}

(14)

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 □ the following equations (eigenvalue problem) is satisfied:

\begin{matrix} □ P (R) = q P (R), □^{- 1} P (R) = q^{- 1} P (R), F (□) P (R) = F (q) P (R), q \neq 0, \end{matrix}

(15)

where q is a parameter of the same dimensionality as □. Applying (15) to equations (12)–(14), we have

\begin{matrix} W = 2 F (q) P^{'} (R) P (R), F (q) = \sum_{n = 1}^{+ \infty} (f_{n} q^{n} + f_{- n} q^{- n}), \end{matrix}

(16)

\begin{matrix} Ω_{μ ν} = F^{'} (q) S_{μ ν} (P, P), F^{'} (q) = \sum_{n = 1}^{+ \infty} n f_{n} q^{n - 1} - \sum_{n = 1}^{+ \infty} n f_{- n} q^{- n - 1}, \end{matrix}

(17)

\begin{matrix} G_{μ ν} + Λ g_{μ ν} + F (q) (2 (R_{μ ν} - K_{μ ν}) P P^{'} - \frac{g_{μ ν}}{2} P^{2} (R)) + \frac{1}{2} F^{'} (q) S_{μ ν} (P, P) = 0 . \end{matrix}

(18)

Now, let us take

P (R) = \sqrt{R - 2 Λ}

. Then

P^{'} (R) P (R) = \frac{1}{2}

and

P^{2} (R) = R - 2 Λ .

Finally, EoM (18) become

(G_{μ ν} + Λ g_{μ ν}) (1 + F (q)) + \frac{1}{2} F^{'} (q) S_{μ ν} (\sqrt{R - 2 Λ}, \sqrt{R - 2 Λ}) = 0 .

(19)

If the nonlocal operator satisfies

\begin{matrix} F (q) = - 1, F^{'} (q) = 0, \end{matrix}

(20)

then equations of motion (19) are also satisfied.

According to the above consideration, the main problem is to solve equation

□ \sqrt{R - 2 Λ} = q \sqrt{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

\sqrt{d S}

gravity. In the sequel, we will investigate the corresponding Schwarzschild-de Sitter metric around static spherically symmetric massive body.

3. The Schwarzschild-de Sitter Metric

We are going to explore the Schwarzschild-de Sitter space-time metric in the case of

\sqrt{d S}

gravity given by its action (1).

3.1. General Consideration

First, we want to consider some general aspects of the corresponding Schwarzschild-de Sitter space-time metric. To this end, we start from the usual expression for the Schwarzschild metric in the pseudo-Rimannian manifold

d s^{2} = - A (r) d t^{2} + B (r) d r^{2} + r^{2} d θ^{2} + r^{2} {sin}^{2} θ d φ^{2}, (c = 1) .

(21)

Non-zero elements of metric tensor

g_{μ ν} (r)

are:

\begin{matrix} g_{00} (r) = - A (r), g_{11} (r) = B (r), g_{22} (r) = r^{2}, g_{33} (r) = r^{2} {sin}^{2} θ . \end{matrix}

(22)

Non-zero components of the corresponding Christoffel symbol

Γ_{μ ν}^{α} = \frac{1}{2} g^{α β} (\partial_{μ} g_{ν β} + \partial_{ν} g_{ν β} - \partial_{β} g_{μ ν})

are:

\begin{matrix} Γ_{10}^{0} = Γ_{01}^{0} = \frac{1}{2} \frac{A^{'}}{A}, Γ_{00}^{1} = \frac{1}{2} \frac{A^{'}}{B}, Γ_{11}^{1} = \frac{1}{2} \frac{B^{'}}{B}, Γ_{22}^{1} = - \frac{r}{B}, \end{matrix}

(23)

\begin{matrix} Γ_{33}^{1} = - \frac{r}{B} {sin}^{2} θ, Γ_{12}^{2} = Γ_{21}^{2} = \frac{1}{r}, Γ_{33}^{2} = - sin θ cos θ, \end{matrix}

(24)

\begin{matrix} Γ_{13}^{3} = Γ_{31}^{3} = \frac{1}{r}, Γ_{23}^{3} = Γ_{32}^{3} = \frac{cos θ}{sin θ}, \end{matrix}

(25)

where ’ denotes derivative with respect to radius

r .

Non-zero components of the Ricci tensor

R_{μ ν} = R_{μ α ν}^{α} = \partial_{α} Γ_{μ ν}^{α} - \partial_{ν} Γ_{μ α}^{α} + Γ_{α ρ}^{α} Γ_{μ ν}^{ρ} - Γ_{ν ρ}^{α} Γ_{μ α}^{ρ}

are:

\begin{matrix} R_{00} = \frac{A^{''}}{2 B} - \frac{A^{'} B^{'}}{4 B^{2}} - \frac{A^{' 2}}{4 A B} + \frac{A^{'}}{r B}, R_{11} = - \frac{A^{''}}{2 A} + \frac{A^{'} B^{'}}{4 A B} + \frac{A^{' 2}}{4 A^{2}} + \frac{B^{'}}{r B}, \end{matrix}

(26)

\begin{matrix} R_{22} = - \frac{r A^{'}}{2 A B} + \frac{r B^{'}}{2 B^{2}} - \frac{1}{B} + 1, R_{33} = (- \frac{r A^{'}}{2 A B} + \frac{r B^{'}}{2 B^{2}} - \frac{1}{B} + 1) {sin}^{2} θ, \end{matrix}

(27)

where as a result of spherical symmetry

R_{33} = R_{22} {sin}^{2} θ .

The corresponding Ricci curvature

R = g^{μ ν} R_{μ ν}

\begin{matrix} R = \frac{2}{r^{2}} - \frac{1}{B (r)} (\frac{2}{r^{2}} + \frac{2 A^{'} (r)}{r A (r)} - \frac{A^{'} {(r)}^{2}}{2 A {(r)}^{2}} - \frac{2 B^{'} (r)}{r B (r)} - \frac{A^{'} (r) B^{'} (r)}{2 A (r) B (r)} + \frac{A^{''} (r)}{A (r)}) . \end{matrix}

(28)

Equation that should be solved is

\begin{matrix} □ u (r) = \frac{1}{B (r)} (▵ u (r) + \frac{1}{2} (\frac{A^{'} (r)}{A (r)} - \frac{B^{'} (r)}{B (r)}) u^{'} (r)) = q u (r), u (r) = \sqrt{R - 2 Λ}, \end{matrix}

(29)

where

\begin{matrix} ▵ = \frac{1}{r^{2}} \frac{\partial}{\partial r} [r^{2} \frac{\partial}{\partial r}] = \frac{\partial^{2}}{\partial r^{2}} + \frac{2}{r} \frac{\partial}{\partial r} \end{matrix}

(30)

is the Laplace operator in spherical coordinate r.

Note that in eigenvalue problem (29) parameter q has the same dimension as operator □, 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 (□)

, defined in (2), satisfies conditions

F (q) = - 1

and

F^{'} (q) = 0

, introduced in (20), if we take it as

\begin{matrix} F (□) = e [a \frac{□}{ζ Λ} e^{(- \frac{□}{ζ Λ})} + b \frac{ζ Λ}{□} e^{(- \frac{ζ Λ}{□})}], where a + b = - 1 . \end{matrix}

(31)

3.2. 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

\begin{matrix} A (r) = A_{0} (r) = 1 - \frac{μ}{r} - \frac{Λ r^{2}}{3}, B (r) = B_{0} (r) = \frac{1}{A_{0} (r)} = \frac{1}{1 - \frac{μ}{r} - \frac{Λ r^{2}}{3}}, μ = \frac{2 G M}{c^{2}} . \end{matrix}

(32)

It makes sense to suppose that solution of equation (29) is of the form

\begin{matrix} A (r) = A_{0} (r) - α (r), B (r) = \frac{1}{A_{0} (r) - β (r)}, \end{matrix}

(33)

where

α (r)

and

β (r)

are some dimensionless functions. When

ζ \to 0

then nonlocal operator (31) tends to zero and consequently nonlocal de Sitter

\sqrt{d S}

gravity model (1) becomes local. Hence, it must be that

A (r) \to A_{0} (r)

and

B (r) \to B_{0} (r)

when

ζ \to 0,

that is

α (r) \to 0

and

β (r) \to 0

ζ \to 0

Replacing

A = A_{0} - α (r)

and

B = \frac{1}{A_{0} - β (r)}

in scalar curvature R (28) and in operator □ of equation (29), we obtain

\begin{matrix} R & = \frac{2}{r^{2}} (1 - A_{0} + β) + 2 \frac{A_{0} - β}{A_{0} - α} (A_{0}^{'} - α^{'}) (\frac{1}{4} \frac{A_{0}^{'} - α^{'}}{A_{0} - α} - \frac{1}{r}) \end{matrix}

\begin{matrix} - 2 (A_{0}^{'} - β^{'}) (\frac{1}{4} \frac{A_{0}^{'} - α^{'}}{A_{0} - α} + \frac{1}{r}) - \frac{A_{0} - β}{A_{0} - α} (A_{0}^{''} - α^{''}), \end{matrix}

(34)

\begin{matrix} □ u & = (A_{0} - β) ▵ u + \frac{1}{2} [\frac{A_{0} - β}{A_{0} - α} (A_{0}^{'} - α^{'}) + A_{0}^{'} - β^{'}] u^{'} = q u, u = \sqrt{R - 2 Λ} \end{matrix}

(35)

If we substitute expressions (34) and (35) in equation (29) then we will get a differential equation with two unknown functions:

α (r \sqrt{ζ Λ})

and

β (r \sqrt{ζ Λ})

. 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) = \frac{1}{A_{0} (r)}

. Hence, there is a sense to take also

B (r) = \frac{1}{A (r)}

in the nonlocal case and it yields

\begin{matrix} R (r) & = \frac{1}{r^{2}} [2 - 2 A (r) - 4 r A^{'} (r) - r^{2} A^{''} (r)] = \frac{1}{r^{2}} \frac{\partial^{2}}{\partial r^{2}} [r^{2} (1 - A (r))], \end{matrix}

(36)

\begin{matrix} □ u (r) & = A (r) u^{''} (r) + (A^{'} (r) + \frac{2}{r} A (r)) u^{'} (r) = \frac{1}{r^{2}} \frac{\partial}{\partial r} [r^{2} A (r) \frac{\partial u}{\partial r}] . \end{matrix}

(37)

In fact, it means that we take

\begin{matrix} β (r) = α (r) . \end{matrix}

(38)

Employing (38) in (34) and (35), we get

\begin{matrix} R & = 4 Λ + \frac{2 α}{r^{2}} + \frac{4 α^{'}}{r} + α^{''}, \end{matrix}

(39)

\begin{matrix} □ u & = (A_{0} - α) ▵ u + (A_{0}^{'} - α^{'}) u^{'} = q u, u = \sqrt{R - 2 Λ} . \end{matrix}

(40)

We are interested in finding function

α (r)

and the next step should be substitution of

\begin{matrix} u = \sqrt{R - 2 Λ} = \sqrt{2 Λ + \frac{2 α}{r^{2}} + \frac{4 α^{'}}{r} + α^{″}} \end{matrix}

(41)

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 □ can be replaced by the Laplacian ▵ in equation (40). In such case we will take

A (r) \approx 1

in (40), that is

\begin{matrix} A (r) = A_{0} (r) - α (r) = 1 - \frac{μ}{r} - \frac{Λ r^{2}}{3} - α (r) \approx 1, \end{matrix}

(42)

what makes sense if the following is satisfied:

\begin{matrix} \frac{μ}{r} ≪ 1, \frac{Λ r^{2}}{3} ≪ 1, | α (r) | ≪ 1 . \end{matrix}

(43)

Applying approximation (42) in (40), we get the following simple equation linear in

u (r)

\begin{matrix} ▵ u = q u, that is \frac{\partial^{2} u}{\partial r^{2}} + \frac{2}{r} \frac{\partial u}{\partial r} = q u, u = \sqrt{R - 2 Λ} . \end{matrix}

(44)

Note that requested function

α (r)

is contained in

u (r) = \sqrt{R - 2 Λ}

, through (39):

R = 4 Λ + \frac{2 α}{r^{2}} + \frac{4 α^{'}}{r} + α^{''}

. Hence, the next step is linearization of

\sqrt{R - 2 Λ}

, and we get

\begin{matrix} \sqrt{R - 2 Λ} = i \sqrt{2 Λ - R} \approx i \sqrt{2 Λ} (1 - \frac{R}{4 Λ}) = - \frac{i}{2 \sqrt{2 Λ}} (R - 4 Λ), if | \frac{R}{2 Λ} | ≪ 1 . \end{matrix}

(45)

It is worth noting that approximation (45) was successfully applied in [24,25].

Now, replacing

\sqrt{R - 2 Λ}

in (44) by

- \frac{i}{2 \sqrt{2 Λ}} (R - 4 Λ)

, we get

\begin{matrix} {(R - 4 Λ)}^{''} + \frac{2}{r} {(R - 4 Λ)}^{'} = q (R - 4 Λ) . \end{matrix}

(46)

The next step is to replace scalar curvature R by

R = 4 Λ + \frac{2 α}{r^{2}} + \frac{4 α^{'}}{r} + α^{''}

(see (39)) in (46) and we obtain

\begin{matrix} α^{''''} + \frac{6}{r} α^{'''} + \frac{2}{r^{2}} α^{''} - \frac{4}{r^{3}} α^{'} + \frac{4}{r^{4}} α = q (α^{''} + \frac{4}{r} α^{'} + \frac{2}{r^{2}} α), \end{matrix}

(47)

that is a linear differential equation of the fourth order.

Equation (47) has a general solution of the form

\begin{matrix} α (r) = \frac{C_{1}}{r} + \frac{C_{2}}{r^{2}} + C_{3} e^{- \sqrt{q} r} (\frac{1}{q r} + \frac{2}{q^{\frac{3}{2}} r^{2}}) + C_{4} e^{\sqrt{q} r} (\frac{1}{q r} - \frac{2}{q^{\frac{3}{2}} r^{2}}), q = ζ Λ . \end{matrix}

(48)

There are four constants (

C_{1}, C_{2}, C_{3}, C_{4}

) and we have chosen them so that the appropriate particular solution

α (r) \to 0

when

ζ \to 0

and that it has some physical meaning.

To exclude term with

e^{\sqrt{q} r}

in (48), we took

C_{4} = 0

, since this exponential function increases indefinitely for very large values of r.

3.2.1. Case $C_{1} = - \frac{δ}{\sqrt{q}}, C_{2} = \frac{2 δ}{q}, C_{3} = - δ \sqrt{q}; C_{4} = 0$ .

In this case solution for

α (r)

\begin{matrix} α (r) = - \frac{δ}{\sqrt{q} r} (1 + e^{- \sqrt{q} r}) + \frac{2 δ}{q r^{2}} (1 - e^{- \sqrt{q} r}), q = ζ Λ, \end{matrix}

(49)

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

ζ \to 0

, it is useful to expand exponential function

e^{- \sqrt{q} r}

into the Taylor series. We have

\begin{matrix} α (r) & = - \frac{δ}{\sqrt{q} r} (1 + e^{- \sqrt{q} r}) + \frac{2 δ}{q r^{2}} (1 - e^{- \sqrt{q} r}) \\ = δ (- \frac{2}{\sqrt{q} r} + 1 - \frac{\sqrt{q} r}{2} + \frac{q r^{2}}{6} - \frac{\sqrt{q} q r^{3}}{24} + \dots) \end{matrix}

(50)

\begin{matrix} + δ (\frac{2}{\sqrt{q} r} - 1 + \frac{\sqrt{q} r}{3} - \frac{q r^{2}}{12} + \dots) \end{matrix}

(51)

\begin{matrix} = - \frac{δ \sqrt{q} r}{6} + \frac{δ q r^{2}}{12} - \dots = - \frac{δ \sqrt{ζ Λ} r}{6} + \frac{δ ζ Λ r^{2}}{12} - \dots \end{matrix}

(52)

As it follows from (52), we conclude that

α (r) \to 0

when

ζ \to 0

According to (49), we get

\begin{matrix} A (r) = 1 - \frac{μ}{r} - \frac{Λ r^{2}}{3} + \frac{δ}{\sqrt{q} r} (1 + e^{- \sqrt{q} r}) - \frac{2 δ}{q r^{2}} (1 - e^{- \sqrt{q} r}), q = ζ Λ, \end{matrix}

(53)

where

μ = \frac{2 G M}{c^{2}} .

When

ζ \to 0

, obtained expression (53) for

A (r)

tends to

A_{0} (r)

, as necessary.

3.3. 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

\sqrt{d S}

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

\begin{matrix} Φ (r) & = \frac{c^{2}}{2} (1 - A (r)) = \frac{G M}{r} + \frac{Λ c^{2} r^{2}}{6} + \frac{c^{2}}{2} α (r) \\ = \frac{G M}{r} + \frac{Λ c^{2} r^{2}}{6} - \frac{δ c^{2}}{2 \sqrt{q} r} (1 + e^{- \sqrt{q} r}) + \frac{δ c^{2}}{q r^{2}} (1 - e^{- \sqrt{q} r}) . \end{matrix}

(54)

Note that here

Φ (r)

is the intensity of the gravitational potential.

The corresponding gravitational acceleration for potential (54) is

\begin{matrix} a_{g} (r) = - \frac{\partial Φ}{\partial r} = \frac{G M}{r^{2}} - \frac{Λ c^{2} r}{3} + \frac{δ c^{2}}{\sqrt{q} r^{2}} (\frac{2}{\sqrt{q} r} - \frac{1}{2}) - \frac{δ c^{2}}{r} (\frac{1}{2} + \frac{3}{2 \sqrt{q} r} + \frac{2}{q r^{2}}) e^{- \sqrt{q} r} . \end{matrix}

(55)

Velocity of the rotation curve

\bar{v} (r)

follows from equality

\frac{{\bar{v}}^{2} (r)}{r} = a_{g} (r)

and it is

\begin{matrix} \bar{v} (r) = \sqrt{a_{g} (r) r} = c \sqrt{\frac{G M}{c^{2} r} - \frac{Λ r^{2}}{3} + \frac{δ}{\sqrt{q} r} (\frac{2}{\sqrt{q} r} - \frac{1}{2}) - δ (\frac{1}{2} + \frac{3}{2 \sqrt{q} r} + \frac{2}{q r^{2}}) e^{- \sqrt{q} r}} . \end{matrix}

(56)

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

\bar{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 \sim 10^{8} M_{⊙}

and there will be no significant changes in the calculated circular velocity

\bar{v} (r)

at very large r.

3.3.1. 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

\bar{v}

using (56) and relative error

δ v = \frac{| v - \bar{v} |}{v}

. A pictorial comparison of measured and calculated velocities is presented in Figure 2.

3.3.2. Spiral Galaxy M33 Case

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

\bar{v} (r)

using formula (56). Measured and computed data are presented in Table 2 and illustrated in Figure 3.

4. Discussion and Concluding Remarks

This paper presents the results of our research regarding the Schwarzschild-de Sitter metric of the nonlocal

\sqrt{d S}

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

\sqrt{d S}

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

\bar{v} (r)

depends on three terms: (i)

\frac{G M}{c^{2} r}

, (ii)

- \frac{Λ r^{2}}{3}

and (iii)

\frac{δ}{\sqrt{q} r} (\frac{2}{\sqrt{q} r} - \frac{1}{2}) - δ (\frac{1}{2} + \frac{3}{2 \sqrt{q} r} + \frac{2}{q r^{2}}) e^{- \sqrt{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 \times 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 general solution (48) of equation (47) was found.
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 Figure 2 and Figure 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

\sqrt{d S}

, we will continue with the further study of this model of nonlocal de Sitter gravity.

Author Contributions

All authors have equally contributed to conceptualization and methodology, formal analysis and investigation, written and editing of the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This research was partially funded by the Ministry of Science, Technological Development and Innovation of the Republic of Serbia, grant numbers: 451-03-66/2024-03/ 200104 and 451-03-65/2024-03/ 200138. This paper is also based on work partially related to the COST Action CA21136, “Addressing observational tensions in cosmology with systematics and fundamental physics” (CosmoVerse) supported by COST (European Cooperation in Science and Technology).

Acknowledgments

The authors are grateful to D. Stojkovic, P. Jovanovic and D. Borka for helpful discussions.

Conflicts of Interest

The authors declare no conflicts of interest.

References

Wald, R.M. General Relativity; University of Chicago Press: Chicago, IL, USA, 1984. [Google Scholar]
Ellis, G.F.R. 100 years of general relativity. In General Relativity and Gravitation: A Centennial Perspective; Cambridge University Press: Cambridge, UK, 2015; pp. 10–19. [Google Scholar]
Robson, B.A. Introductory chapter: Standard Model of Cosmology. In Redefining Standard Model Cosmology; 2019; pp. 1–4. [Google Scholar] [CrossRef]
Oks, E. Brief review of recent advances in understanding dark matter and dark energy. New Ast. Rev. 2021, 93, 101632. [Google Scholar] [CrossRef]
Novello, M.; Bergliaffa, S.E.P. Bouncing cosmologies. Phys. Rep. 2008, 463, 127–213. [Google Scholar] [CrossRef]
Stelle, K.S. Renormalization of higher derivative quantum gravity. Phys. Rev. D 1977, 16, 953. [Google Scholar] [CrossRef]
Clifton, T.; Ferreira, P.G.; Padilla, A.; Skordis, C. Modified gravity and cosmology. Phys. Rep. 2012, 513, 1–189. [Google Scholar] [CrossRef]
Nojiri, S.; Odintsov, S.D. Unified cosmic history in modified gravity: from F(R) theory to Lorentz non-invariant models. Phys. Rep. 2011, 505, 59–144. [Google Scholar] [CrossRef]
Nojiri, S.; Odintsov, S.D.; Oikonomou, V.K. Modified gravity theories on a nutshell: Inflation, bounce and late-time evolution. Phys. Rep. 2017, 692, 1–104. [Google Scholar] [CrossRef]
Dimitrijevic, I.; Dragovich, B.; Grujic, J.; Rakic, Z. On modified gravity. Springer Proc. Math. Stat. 2013, 36, 251–259. [Google Scholar]
Sotiriou, T.P.; Faraoni, V. f(R) theories of gravity. Rev. Mod. Phys. 2010, 82, 451–497. [Google Scholar] [CrossRef]
Capozziello, S.; Bajardi, F. Nonlocal gravity cosmology: An overview. Int. J. Mod. Phys. D 2022, 2022 31, 2230009. [Google Scholar] [CrossRef]
Dragovich, B. On Nonlocal modified gravity and cosmology. Springer Proc. Math. Stat. 2014, 111, 251–262. [Google Scholar]
Biswas, T.; Koivisto, T.; Mazumdar, A. Towards a resolution of the cosmological singularity in non-local higher derivative theories of gravity. J. Cosmol. Astropart. Phys. 2010, 1011, 008. [Google Scholar] [CrossRef]
Biswas, T.; Gerwick, E.; Koivisto, T.; Mazumdar, A. Towards singularity and ghost free theories of gravity. Phys. Rev. Lett. 2012, 108, 031–101. [Google Scholar] [CrossRef] [PubMed]
Biswas, T.; Koshelev, A.S.; Mazumdar, A.; Vernov, S.Y. Stable bounce and inflation in non-local higher derivative cosmology. J. Cosmol. Astropart. Phys. 2012, 08, 024. [Google Scholar] [CrossRef]
Biswas, T.; Conroy, A.; Koshelev, A.S.; Mazumdar, A. Generalized gost-free quadratic curvature gravity. Class. Quantum Grav. 2014, 31, 159501. [Google Scholar] [CrossRef]
Biswas, T.; Koshelev, A.S.; Mazumdar, A. Consistent higher derivative gravitational theories with stable de Sitter and Anti-de Sitter backgrounds. Phys. Rev. 2017, D 95, 043533. [Google Scholar] [CrossRef]
Dimitrijevic, I.; Dragovich, B.; Koshelev, A.S.; Rakic, Z.; Stankovic, J. Cosmological solutions of a nonlocal square root gravity. Phys. Lett. B 2019, 797, 134848. [Google Scholar] [CrossRef]
Kol’ar, I.; Torralba, F.J.M.; Mazumdar, A. New non-singular cosmological solution of non-local gravity. Phys. Rev. D 2022, 105, 044045. [Google Scholar] [CrossRef]
Koshelev, A.S.; Kumar, K.S.; Starobinsky, A.A. Analytic infinite derivative gravity, R2-like inflation, quantum gravity and CMB. Int. J. Mod. Phys. D 2020, 29, 2043018. [Google Scholar] [CrossRef]
Modesto, L.; Rachwal, L. Super-renormalizable and finite gravitational theories. Nucl. Phys. B 2014, 889, 228. [Google Scholar] [CrossRef]
Dimitrijevic, I.; Dragovich, B.; Grujic, J.; Koshelev, A. S.; Rakic, Z. Cosmology of non-local f(R) gravity. Filomat 2019, 33, 1163–1178. [Google Scholar] [CrossRef]
Dimitrijevic, I.; Dragovich, B.; Koshelev, A.S.; Rakic, Z.; Stankovic, J. Some cosmological solutions of a new nonlocal gravity model. Symmetry 2020, 12, 917. [Google Scholar] [CrossRef]
Dimitrijevic, I.; Dragovich, B.; Rakic, Z.; Stankovic, J. New cosmological solutions of a nonlocal gravity model. Symmetry 2021, 1, 0. [Google Scholar] [CrossRef]
Dimitrijevic, I.; Dragovich, B.; Stankovic, J.; Koshelev, A. S.; Rakic, Z. On nonlocal modified gravity and its cosmological solutions. Springer Proc. Math. Stat. 2016, 191, 35–51. [Google Scholar]
Deser, S.; Woodard, R. Nonlocal cosmology. Phys. Rev. Lett. 2007, 99, 111–301. [Google Scholar] [CrossRef] [PubMed]
Woodard, R.P. Nonlocal models of cosmic acceleration. Found. Phys. 2014, 44, 213–233. [Google Scholar] [CrossRef]
Belgacem, E.; Dirian, Y.; Foffa, S.; Maggiore, M. Nonlocal gravity. Conceptual aspects and cosmological predictions. J. Cosmol. Astropart. Phys. 2018, 2018, 002. [Google Scholar] [CrossRef]
Dimitrijevic, I.; Dragovich, B.; Rakic, Z.; Stankovic, J. Nonlocal de Sitter gravity and its exact cosmological solutions. JHEP 2022, 2022, 054–12. [Google Scholar] [CrossRef]
Dragovich, B.; Khrennikov, A.Y.; Kozyrev, S.V.; Volovich, I.V.; Zelenov, E.I. p-Adic mathematical physics: the first 30 years. p-Adic Num. Ultrametr. Anal. Appl. 2017, 9, 87–121. [Google Scholar] [CrossRef]
Dimitrijevic, I.; Dragovich, B.; Rakic, Z.; Stankovic, J. On the Schwarzschild-de Sitter metric of nonlocal de Sitter gravity. Filomat 2023, 37, 25–8641. [Google Scholar] [CrossRef]
Dimitrijevic, I.; Dragovich, B.; Rakic, Z.; Stankovic, J. Variations of infinite derivative modified gravity. Springer Proc. Math. Stat. 2018, 263, 91–111. [Google Scholar]
Corbelli1, E.; Salucci, P. The extended rotation curve and the dark matter halo of M33. Mon. Not. R. Astron. Soc. 2000, 311, 441–447. [Google Scholar] [CrossRef]
Bottema, R.; Pestaña, J.L.G. The distribution of dark and luminous matter inferred from extended rotation curves. MNRAS 2015, 448, 2566–2593. [Google Scholar] [CrossRef]
Genzel, R. et al. Strongly baryon-dominated disk galaxies at the peak of galaxy formation ten billion years ago. Nature 2017, 543, 397–401. [Google Scholar] [CrossRef]
Banik, I.; Thies, I.; Candlish, G.; Famaey, B.; Ibata, R.; Kroupa, P. The global stability of M33 in MOND. ApJ 2020, 2020 905, 135. [Google Scholar] [CrossRef]
Dai, D.-C.; Glenn Starkman, G.; Stojkovic, D. Milky Way and M31 rotation curves: ΛCDM versus MOND. Phys. Rev. D 2022, 105, 104067. [Google Scholar] [CrossRef]
D’Agostino, R.; Jusufi, K.; Capozziello, S. Testing Yukawa cosmology at the Milky Way and M31 galactic scales. arXiv 2024, arXiv:2404.01846. [Google Scholar] [CrossRef]
Jiao, Y.; Hammer, F.; Wang, H.; Wang, J.; Amram, P.; Chemin, L.; Yang, Y. Detection of the Keplerian decline in the Milky Way rotation curve. A&A 2023, 678, A208. [Google Scholar]
Kam, S.Z.; Carignan, C.; Chemin, L.; Foster, T.; Elson, E.; Jarrett, T.H. HI kinematics and mass distribution of Messier 33. AJ 2017, 154, 41. [Google Scholar] [CrossRef]
Granados, A.; Torres, D.; Castañeda, L.; Henao, O.L.-J.; Vanegas, S. GalRotpy: an educational tool to understand and parametrize the rotation curve and gravitational potential of disk-like galaxies. New Astr. 2021, 82, 101456. [Google Scholar] [CrossRef]

Figure 1. We consider the Schwarzschild-de Sitter metric of nonlocal

\sqrt{d S}

gravity at the distances far from a spherically symmetric massive body.

Figure 1. We consider the Schwarzschild-de Sitter metric of nonlocal

\sqrt{d S}

gravity at the distances far from a spherically symmetric massive body.

Figure 2. Rotation curve for the Milky Way galaxy. Red points are measured observational values [40] and blue line is computed

\bar{v} (r)

by formula (56), where

δ = 1.9 \times 10^{- 5}

ζ = 4.4 \times 10^{10}

Λ = 10^{- 52} m^{- 2}

and

M = 4.28 \times 10^{6} M_{⊙}

Figure 2. Rotation curve for the Milky Way galaxy. Red points are measured observational values [40] and blue line is computed

\bar{v} (r)

by formula (56), where

δ = 1.9 \times 10^{- 5}

ζ = 4.4 \times 10^{10}

Λ = 10^{- 52} m^{- 2}

and

M = 4.28 \times 10^{6} M_{⊙}

Figure 3. Rotation curve for spiral galaxy M33. Red points are measured observational values and blue line is computed

\bar{v} (r)

by formula (56), where

δ = 5.7 \times 10^{- 6}

ζ = 3.62 \times 10^{10}

Λ = 10^{- 52} m^{- 2}

and

M = 1.5 \times 10^{3} M_{⊙}

Figure 3. Rotation curve for spiral galaxy M33. Red points are measured observational values and blue line is computed

\bar{v} (r)

by formula (56), where

δ = 5.7 \times 10^{- 6}

ζ = 3.62 \times 10^{10}

Λ = 10^{- 52} m^{- 2}

and

M = 1.5 \times 10^{3} M_{⊙}

Table 1. Milky Way rotation curve data, from [40] and this work.

$r$ [kpc]	$v$ [km/s]	$Δ v$ [km/s]	$\bar{v}$ [km/s]	relative error $δ v$ [%]
9.5	221.75	3.17	217.36	1.98
10.5	223.32	3.02	220.19	1.40
11.5	220.72	3.47	221.93	0.55
12.5	222.92	3.19	222.72	0.09
13.5	224.16	3.48	222.66	0.67
14.5	221.60	4.20	221.85	0.11
15.5	218.79	4.75	220.37	0.72
16.5	216.38	4.96	218.28	0.88
17.5	213.48	6.13	215.63	1.01
18.5	209.17	4.42	212.47	1.58
19.5	206.25	4.63	208.83	1.25
20.5	202.54	4.40	204.77	1.10
21.5	197.56	4.62	200.29	1.38
22.5	197.00	3.81	195.42	0.80
23.5	191.62	12.95	190.17	0.75
24.5	187.12	8.06	184.57	1.36
25.5	181.44	19.58	178.62	1.55
26.5	175.68	24.68	172.32	1.91

Table 2.

M 33

galaxy data, from [41] and this work.

Table 2.

M 33

galaxy data, from [41] and this work.

$r$ [kpc]	$v$ [km/s]	$Δ v$ [km/s]	$\bar{v}$ [km/s]	relative error $δ v$ [%]	$r$ [kpc]	$v$ [km/s]	$Δ v$ [km/s]	$\bar{v}$ [km/s]	relative error $δ v$ [%]
0.5	42.0	2.4	35.62	15.18	12.2	115.7	9.6	120.69	4.31
1.0	58.8	1.5	49.61	15.63	12.7	115.1	7.7	121.05	5.17
1.5	69.4	0.4	59.83	13.79	13.2	117.1	5.1	121.30	3.58
2.0	79.3	4.0	68.02	14.22	13.7	118.2	3.2	121.45	2.75
2.4	86.7	1.8	73.59	15.12	14.2	118.4	1.4	121.50	2.62
2.9	91.4	3.1	79.64	12.86	14.7	118.2	1.8	121.47	2.76
3.4	94.2	4.8	84.90	9.88	15.1	117.5	2.4	121.38	3.30
3.9	96.5	5.5	89.51	7.25	15.6	119.6	0.8	121.19	1.33
4.4	99.8	3.9	93.58	6.23	16.1	118.6	1.5	120.93	1.96
4.9	102.1	1.7	97.21	4.80	16.6	122.6	0.5	120.59	1.64
5.4	103.6	0.4	100.44	3.05	17.1	124.1	2.9	120.17	3.16
5.9	105.9	0.7	103.32	2.44	17.6	125.0	2.2	119.69	4.24
6.4	105.7	1.7	105.90	0.19	18.1	125.5	2.5	119.15	5.06
6.8	106.8	2.2	107.76	0.90	18.6	125.2	8.1	118.54	5.32
7.3	107.3	3.0	109.86	2.39	19.1	122.0	9.8	117.87	3.38
7.8	108.3	4.0	111.73	3.17	19.5	120.4	8.5	117.29	2.58
8.3	109.7	4.0	113.34	3.37	20.0	114.0	26.6	116.52	2.21
8.8	112.0	4.8	114.86	2.55	20.5	110.0	34.6	115.70	5.18
9.3	116.1	2.2	116.15	0.04	21.0	98.7	27.4	114.82	16.33
9.8	117.2	2.5	117.27	0.06	21.5	100.1	33.4	113.89	13.77
10.3	116.5	6.5	118.24	1.49	22.0	104.3	35.2	112.91	8.25
10.8	115.7	8.1	119.07	2.91	22.5	101.2	27.4	111.88	10.56
11.2	117.4	8.2	119.63	1.90	23.0	123.5	39.1	110.81	10.27
11.7	116.8	8.9	120.22	2.93	23.5	115.3	26.7	109.69	4.86

Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

© 2024 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).

Copyright: This open access article is published under a Creative Commons CC BY 4.0 license, which permit the free download, distribution, and reuse, provided that the author and preprint are cited in any reuse.

MDPI Initiatives

Important Links

Choose an area of interest and we will send you notifications of new preprints at your preferred frequency.

Disclaimer

The Schwarzschild-de Sitter Metric of Nonlocal dS Gravity

Abstract

1. Introduction

2. dS Nonlocal Gravity

2.1. Equations of Motion

3. The Schwarzschild-de Sitter Metric

3.1. General Consideration

3.2. Solutions

3.2.1. Case C 1 = − δ q , C 2 = 2 δ q , C 3 = − δ q ; C 4 = 0 .

3.3. The Rotation Curves of Spiral Galaxies

3.3.1. Milky Way Case

3.3.2. Spiral Galaxy M33 Case

4. Discussion and Concluding Remarks

Author Contributions

Funding

Acknowledgments

Conflicts of Interest

References

MDPI Initiatives

Important Links

Subscribe

2. $\sqrt{dS}$ Nonlocal Gravity

3.2.1. Case $C_{1} = - \frac{δ}{\sqrt{q}}, C_{2} = \frac{2 δ}{q}, C_{3} = - δ \sqrt{q}; C_{4} = 0$ .