Next Article in Journal
Path Integral Action for a Resonant Detector of Gravitational Waves in the Generalized Uncertainty Principle Framework
Previous Article in Journal
Generating CP Violation from a Modified Fridberg-Lee Model
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Calculation of the Cosmological Constant for the Planetary System in Schwarzschild’s Cosmological Model

by
Alvaro Humberto Salas Salas
1,
Jairo Ernesto Castillo Hernandez
2,* and
Jorge Enrique Pinzon Quintero
3
1
FIZMAKO Research Group, Department of Mathematics and Statistics, Universidad Nacional de Colombia Sede Manizales, Manizales 500001, Colombia
2
FIZMAKO Research Group, Universidad Distrital Francisco José de Caldas, Bogotá 11021, Colombia
3
FIZMAKO Research Group, Universidad del Tolima, Ibagué 730001, Colombia
*
Author to whom correspondence should be addressed.
Universe 2022, 8(9), 449; https://doi.org/10.3390/universe8090449
Submission received: 12 July 2022 / Revised: 18 August 2022 / Accepted: 22 August 2022 / Published: 28 August 2022
(This article belongs to the Section Cosmology)

Abstract

:
In this work, the static cosmological model of the Schwarzschild solution for the solar system is proposed taking into account the cosmological constant in the equation of the general theory of relativity (GTR) proposed by A. Einstein. We found the nonlinear differential equation that describes the behavior of the planets around the Sun; this is solved exactly by the Jacobi and Weierstrass elliptic functions. The obtained solution allows for us to estimate the value of the cosmological constant knowing the perihelion of the different planets and from different mathematical approaches; that is, the inverse problem is solved. From the obtained results, the Schwarzschild static cosmological model for the solar system is proposed, establishing the Schwarzschild cosmological radius and the curvature limit of the solar system. From the curvature limit, different regions are proposed for the planets, exoplanets, and a region is predicted where the existence of new planets and exoplanets belonging to the solar system is possible. The proposed theory of the static Schwarzschild cosmological model may be of great interest to astronomers, cosmologists, and all those interested in the study of the universe.

1. Introduction

Albert Einstein in 1915 [1,2] formulated the general theory of gravitation (TGR), synthesizing it in the equation
R α β 1 2 R g α β = 8 π G c 4 T α β
He verified that his theory explained the anomaly of the perihelion of Mercury and the curvature of light in the gravitational field of the Sun. The first exact static solution was found by Schwarzschild [3]. This solution led to deducing two great physical phenomena in our planetary system: the perihelion of the planet Mercury, the deviation of light rays in the gravitational field of the Sun, establishing that spacetime in the vicinity of a gravitational mass is curved, later experimentally demonstrated [4,5], this being a great triumph for GTG.
In the following years and recently, nonstatic solutions were found, such as those of: E. Kasner [6,7,8,9], H.P Robertson [10], A.G Walker [11], A. Friedman [12], and G. Lemaitre [13], which are called the Friedmann–Lemaitre–Robertson–Walker metric or FLRW cosmological model.
In 1917, Albert Einstein [14] could not find a static universe of his original equation. This led him to reformulate the original equation of the GTG, introducing a term called the lambda Λ -cosmological constant, which encompass the physical structure of the static universe. Later, in 1937, Albert Einstein rejected this idea once De sitter [15] had established that there is a solution of Equation (1) for empty space, which represents a model of the expanding universe, and it was observed by Edwin Hubbles [16] in the redshift of galaxies, suggesting that the universe was not static, and that Eddington [17] demonstrated the expansion of the universe and that the static universe with cosmological constant was unstable [18,19].
R α β 1 2 R g α β Λ g α β = 8 π G c 4 T α β
Edwin Hubbles’ discovery of the redshift, and subsequent observations of different galaxies and supernovae yielded new results that confirmed the theory of redshift, the discoveries of the radiation of the cosmic background, the acceleration of the universe, the large-scale structure of the universe, dark energy, and dark matter, leading to retaking the importance of the cosmological constant and the search for solutions to the equations proposed by A. Einstein with the cosmological constant.
Many works studied the perihelion of Mercury [20,21] and of other planets with constant zero and nonzero cosmology [22,23,24,25].
This work is organized as follows: In Section 1, Einstein’s equation is solved with the cosmological constant in vacuum with the aim that the work is complete and independent for the reader. In Section 2 are the algebraic equations of geodesics, and the nonlinear differential equation that describes the behavior of a planet around the Sun. In Section 3 and Section 4, the nonlinear differential equation is solved by means of the elliptic Jacobi functions and Weierstrass elliptic function. The cosmological constant is found as a function of the perihelion and aphelion of the planets around the Sun from different perspectives. Lastly, an exhaustive analysis of the obtained results and their consequences in Schwarzschild’s cosmological static model was carried out, for which A. Einstein so yearned and sought.

2. Solution of Einstein’s Equation with the Cosmological Constant

In this section, we propose to derive the Schwarzschild metric while taking into account the cosmological constant. For this purpose, we solve Einstein’s equation in a vacuum while taking into account the cosmological constant
R α β Λ g α β = 0
where R α β is the Ricci tensor, g α β the metric tensor, and Λ the cosmological constant.
Since the Schwarzschild solution is static and isotropic, let us find the solution of (3) of the form:
d s 2 = D ( r ) d t 2 E ( r ) d r 2 F ( r ) r 2 ( d θ 2 + sin 2 θ d ϕ 2 )
Now, we perform a coordinate transformation ( c t , r , θ , ϕ ) > ( c t , ρ , θ , ϕ ) where the new coordinate ρ is defined by the equation
ρ 2 F ( r ) r 2 :
d s 2 = A ( ρ ) c 2 d t 2 B ( ρ ) d ρ 2 ρ 2 ( d θ 2 + sin 2 θ d ϕ 2 )
The metric tensor that corresponds to the linear differential element is clearly
g α β = A ( ρ ) 0 0 0 0 B ( ρ ) 0 0 0 0 ρ 2 0 0 0 0 ρ 2 sin 2 θ
This is the metric of the more general isotropic static field problem. Functions A ( ρ ) and B ( ρ ) are still indeterminate, as is the relationship between coordinates ρ and r. When calculating the components of Ricci tensor R α β = ρ Γ β α ρ β Γ ρ α ρ + Γ ρ λ ρ Γ β α λ Γ β λ ρ Γ ρ α λ , the only nonzero components are:
R 00 = A ( ρ ) 2 B ( ρ ) ρ A ( ρ ) + 2 A ( ρ ) ρ A ( ρ ) B ( ρ ) ρ B ( ρ ) A ( ρ ) 2 4 ρ A ( ρ ) B ( ρ ) 2 R 11 = A ( ρ ) ρ A ( ρ ) + 4 A ( ρ ) B ( ρ ) + ρ B ( ρ ) A ( ρ ) 2 2 A ( ρ ) A ( ρ ) 4 ρ A ( ρ ) 2 B ( ρ ) R 22 = ρ B ( ρ ) A ( ρ ) + ρ A ( ρ ) B ( ρ ) 2 2 A ( ρ ) B ( ρ ) 2 A ( ρ ) B ( ρ ) 2 R 33 = sin 2 ( ϕ ) A ( ρ ) ρ B ( ρ ) 2 2 B ( ρ ) ρ B ( ρ ) A ( ρ ) 2 A ( ρ ) B ( ρ ) 2
Replacing (7) and (6) in (3), we obtain the following system of equations for A ( ρ ) end B ( ρ ) :
2 ρ A ( ρ ) B ( ρ ) A ( ρ ) ρ A ( ρ ) A ( ρ ) B ( ρ ) ρ B ( ρ ) A ( ρ ) 2 + 4 A ( ρ ) B ( ρ ) A ( ρ ) 2 B ( ρ ) 2 + Λ A ( ρ ) = 0 2 ρ A ( ρ ) B ( ρ ) A ( ρ ) + ρ A ( ρ ) A ( ρ ) B ( ρ ) + ρ B ( ρ ) A ( ρ ) 2 + 4 A ( ρ ) 2 B 2 B ( ρ ) 2 Λ B ( ρ ) = 0 ρ B ( ρ ) A ( ρ ) + ρ A ( ρ ) B ( ρ ) 2 A ( ρ ) B ( ρ ) 2 + 2 A ( ρ ) B ( ρ ) 2 2 A ( ρ ) B ( ρ ) Λ ρ 2 = 0 sin 2 ( ϕ ) ρ B ( ρ ) A ( ρ ) ρ A ( ρ ) B ( ρ ) 2 A ( ρ ) B ( ρ ) 2 2 A ( ρ ) B ( ρ ) 2 + 2 A ( ρ ) B ( ρ ) Λ ρ 2 sin 2 ( θ ) = 0
Adding the first and second equations of System (8), we have A ( ρ ) B ( ρ ) + A ( ρ ) B ( ρ ) = ( A B ) = 0 . This means that:
A ( ρ ) B ( ρ ) = k 1 ,
so that B ( ρ ) = k 1 A ( ρ ) . By placing this equation into the third equation of System (8), we obtain:
ρ A ( ρ ) + A ( ρ ) + k 1 Λ ρ 2 k 1 = 0 ,
from where:
A ( ρ ) = k 1 1 3 k 1 ρ 2 Λ + c 1 ρ
In order to determine the values of constants k 1 and c 1 , some physical condition must be used, which is achieved by studying the situation at the Newtonian limit. When Λ = 0 , the Newtonian potential is Φ = G M r ; therefore,
g 00 = A ( ρ ) 1 2 G M c 2 r
Comparing with (11), we see that k 1 = 1 and c 1 ρ = 2 G M c 2 r The set of two variables ( c 1 , ρ ) is controlled by a single equation. This means that the values of variables ( c 1 , ρ ) are undetermined. We then need to arbitrarily determine one of the variables. The most comfortable choice is ρ = r ; consequently, c 1 is equal to 2 G M c 2 . Therefore, we lastly have:
A = 1 B = 1 1 3 r 2 Λ 2 G M c 2 r
Linear Differential Element (4) reads [15,26].
d s 2 = ( 1 1 3 r 2 Λ 2 G M c 2 r ) c 2 d t 2 ( 1 1 3 r 2 Λ 2 G M c 2 r ) 1 d r 2 r 2 ( d θ 2 + sin 2 θ d ϕ 2 ) ,
where
g α β = ( 1 1 3 r 2 Λ 2 G M c 2 r ) 0 0 0 0 ( 1 1 3 r 2 Λ 2 G M c 2 r ) 1 0 0 0 0 r 2 0 0 0 0 r 2 sin 2 θ

3. Geodesic Equations and Equation of Motion

Let us consider the motion of a planet in the gravitational field of a much heavier body (the Sun). The gravitational field of the center in spherical coordinates is given by Linear Differential Element (14).
d s 2 = ( 1 1 3 r 2 Λ s r ) c 2 d t 2 ( 1 1 3 r 2 Λ s r ) 1 d r 2 r 2 ( d θ 2 + sin 2 θ d ϕ 2 ) ,
where s = 2 G M c 2 . With Metric (16), we calculate the Christoffel symbols.
Γ μ ν λ = 1 2 g λ σ ( μ g σ ν + ν g σ μ σ g μ ν )
The nonzero quantities are:
Γ 01 0 = Γ 10 0 = 3 s 2 r 3 Λ 2 Λ r 4 + 6 r 2 6 s r Γ 00 1 = c 2 ( Λ r 3 3 r + 3 s ) ( 2 r 3 Λ 3 s ) 18 r 3 Γ 11 1 = Γ 01 0 Γ 22 1 = r 3 Λ 3 r + s Γ 33 1 = 1 3 ( r 3 Λ 3 r + 3 s ) sin 2 θ Γ 12 2 = Γ 21 2 = 1 r Γ 13 3 = Γ 31 3 = 1 r Γ 33 2 = cos θ sin θ Γ 32 3 = Γ 23 3 = cot θ
Using the equations of geodesic lines d 2 x σ d τ 2 + Γ μ ν σ d x μ d τ d x ν d τ = 0 and taking into account the central symmetry of our problem, any plane through the center can, however, be chosen as the plane θ = π 2 , that is, the orbit can be located in any plane through the center; we obtain the algebraic equations from the geodesics:
2 r t 3 s 2 Λ r 3 2 Λ r 4 + 6 r 2 6 r s + t = 0 . c 2 t 2 Λ r 3 3 r + 3 s 2 Λ r 3 3 s 18 r 3 + ϕ 2 Λ r 3 3 r + s + r ( k ) r 2 3 s 2 Λ r 3 2 Λ r 4 + 6 r 2 6 r s = 0 . 2 r ϕ r + ϕ = 0 .
From the first and third equations of (19), we obtain that:
t = E r 3 s 3 r + r 3 Λ , ϕ = L r 2 ,
where E and L are the constants of integration that represent the effective energy of the system and the angular momentum.
Dividing (14) by d s 2 = c 2 d τ 2 gives
1 = r 2 c 2 1 3 Λ r ( τ ) 2 s r ( τ ) + 1 r ( τ ) 2 ϕ 2 c 2 + t 2 1 3 Λ r ( τ ) 2 s r ( τ ) + 1
Replacing (20) into (21):
L 2 c 2 r ( τ ) 2 r 2 c 2 1 3 Λ r ( τ ) 2 s r ( τ ) + 1 + E 2 r ( τ ) 2 1 3 Λ r ( τ ) 2 s r ( τ ) + 1 Λ r ( τ ) 3 3 r ( τ ) + 3 s 2 1 = 0
Applying chain rule r ( τ ) = d r d φ d φ d τ , r ( τ ) 2 = ( d r d φ ) 2 L 2 r 4 ,
r ( ϕ ) 4 c 2 E 2 9 + 3 Λ L 2 3 c 2 Λ r ( ϕ ) 6 9 c 2 s r ( ϕ ) 3 + 9 L 2 r 2 9 L 2 s r ( ϕ ) + 9 L 2 r ( ϕ ) 2 = 0
Using relation r ( φ ) = 1 u ( φ ) and taking into account that s = 2 G M c 2 , we obtain that
c 2 E 2 9 L 2 2 G M u ( ϕ ) 3 c 2 c 2 Λ 3 L 2 u ( ϕ ) 2 + c 2 L 2 2 G M u ( ϕ ) L 2 Λ 3 + u 2 + u ( ϕ ) 2 = 0
or
2 G M u ( ϕ ) 3 c 2 c 2 Λ 3 L 2 u ( ϕ ) 2 2 G M u ( ϕ ) L 2 + u 2 + u ( ϕ ) 2 = c 2 E 2 9 L 2 c 2 L 2 + Λ 3 = C
Derivating Equation (24) after simplifying (24), we obtain the equation of motion taking into account the cosmological constant:
d 2 u d ϕ 2 + u G M L 2 3 G M c 2 u 2 + c 2 Λ 3 L 2 u 3 = 0

4. Einstein Equation for Planetary Motion

Let us consider the equation
u ( ϕ ) = G M L 2 u ( ϕ ) + 3 G M c 2 u 2 ( ϕ ) , u ( 0 ) = 1 P and u ( 0 ) = 0 .
If we take into account cosmological constant Λ , the equation is modified by adding another term as follows:
u ( ϕ ) = G M L 2 c 2 Λ 3 L 2 u 3 ( ϕ ) u ( ϕ ) + 3 G M c 2 u 2 ( ϕ ) = 0 , u ( 0 ) = 1 P and u ( 0 ) = 0 .
Equations (27) and (28) may both be solved in closed form. However, Equation (28) demands inverting some hyper elliptic Abelian integral and its solution is expressed in terms of a generalized Weierstrass function, which is a difficult task. Let us solve the easier Equation (27). Let
u ( ϕ ) = 1 P S + S cn 2 λ + μ ϕ , μ 2 ( λ + μ ) .
Define the residual function
R ( ϕ ) = α u ( ϕ ) + γ u 2 ( ϕ ) u ( ϕ ) , α = G M L 2 = 3 A P ( A + P ) 2 γ A 2 + A P + P 2 6 A 2 P 2 , γ = 3 G M c 2 .
In these formulas, A stands for the aphelion and P stands for the perihelion of the planet. Number L corresponds to the angular momentum of the planet. We have
R ( ϕ ) = P + P 2 S + P 2 α + γ 2 P S γ + P 2 S 2 γ 2 P 2 S λ P 2 S μ P 2 ( S ( P 2 γ + 2 P S γ 4 P λ ) ) P cn 2 + S ( S γ + 3 μ ) cn 4 , cn = cn λ + μ ϕ , μ 2 ( λ + μ ) .
Equating the coefficients of cn 0 , cn 2 , and cn 4 to zero gives an algebraic system whose solution is easily obtained, as follows:
λ = 2 γ + 2 γ P S + P 4 P , μ = γ S 3 . S = 6 γ 3 P + 12 γ 2 + 9 P 2 + 12 γ P ( 1 4 α P ) 4 γ P .
Then,
u ( ϕ ) = 1 P S + S cn 2 ω ϕ , m ,
where
m = 1 2 3 ( ρ + 2 γ P ) 2 ( ρ 6 γ + 3 P ) = 2 G M ( A P ) 2 G M ( 2 A + P ) A c 2 P . ω = 1 2 P 2 γ P ( 1 + m )   a n d   S = 1 A 3 2 γ + 2 P .
The perihelion shift is given by
Δ G R T = 2 K ( m ) ω 2 π · 23668612128 π · Siederal arc-sec / cy .
Here, Siederal stands for the orbit period of the planet. Since, for small m ,
4 K ( m ) P 2 γ ( m + 1 ) P 2 π 2 π 3 m P 2 γ P + 9 m 16 P 2 γ P + 16 ( 16 3 m ) P 2 γ P
number Δ G R T may be approximated as follows:
Δ ˜ G R T = 47337224256 8 A c 2 P ( 1 υ ) + G M ( 29 A υ 41 A + 19 P υ 7 P ) υ 8 A c 2 P G M ( 29 A + 19 P ) · Siederal arcsec / cy , w h e r e   υ = 1 2 γ P = 1 6 G M c 2 P
Another formula:
Δ ^ G T R = 71005836384 μ c 2 P ( Q + 1 ) + 2 μ Q ( Q + 5 ) c 4 P 2 Siederal arcsec / century , μ = G M , Q = P A .

Exact Value for the Perihelion Shifts

We used the following data:
Sun’s mass: M = 1.9885469241387177397120698128813396785059276218789 × 10 30
Gravitational constant: G = 6.6738399987961360272969834 × 10 11 .
Speed of light in vacuum: c = 299792458
The results are depicted in Table 1.
Using NASA data, we obtained the following values for perihelion shifts:
Figure 1 illustrates the way in which Mercury’s perihelion moves.
Remark 1.
Equation (27) may also be written as
d u d ϕ = 2 γ 3 u ( ϕ ) 1 A 1 P u ( ϕ ) 3 2 γ 1 A 1 P u ( ϕ ) .
Integrating it gives
T = 2 ( ϕ ( 1 / P ) ϕ ( 1 / A ) ) = 4 K 2 G M ( A P ) 2 G M ( 2 A + P ) A c 2 P P 2 γ P 1 + 2 G M ( A P ) 2 G M ( 2 A + P ) A c 2 P = 2 6 A K A P P ( A δ 0 1 ) γ A δ 0 1 ,
δ 0 = 3 2 γ 1 A 1 P .
This formula also allows for us to evaluate the perihelion shift of the planet.
Remark 2.
The solution to Equation (27) is also expressed in terms of Weierstrass elliptic function as follows:
( ϕ ) = 1 P ( A P ) ( [ 3 P 4 γ ] A 2 γ P ) A 2 P ( P 2 γ ) 1 + 12 P P 2 γ ( ϕ ; g 2 , g 3 ) ,
where
g 2 = 4 γ 2 A 2 + A P + P 2 + 3 A 2 P 2 6 A γ P ( A + P ) 36 A 2 P 2 . g 3 = ( A 2 γ ) ( 2 γ P ) ( γ ( A + P ) A P ) 216 A 2 P 2 .
The solution is periodic with period
T = 2 r d x 4 x 3 g 2 x g 3 ,
where r is the greatest real root to the cubic 4 x 3 g 2 x g 3 = 0 . This cubic has the following real roots:
x 1 = γ 6 A 1 12 , x 2 = γ 6 P 1 12 , x 3 = 1 6 γ ( A + P ) 6 A P .
The greatest real root is r = x 3 ; so,
T = 2 x 3 d x 4 x 3 g 2 x g 3 = 2 6 A P F sin 1 2 ( A P ) γ 3 A P 2 γ P 4 A γ , 2 ( 2 A + P ) γ 3 A P 2 ( A P ) γ γ P A .
For Mercury (see Table 2), we have the following value from (46):
T 2 π = 5.018176594262513 × 10 7
Remark 3.
Solution (29) may be trigonometrically approximated in as follows. We have the approximation
cn ( t , m ) c o s m ( t ) : = 1 + λ + μ cos 1 + λ + μ 1 μ t 1 + λ cos 2 1 + λ + μ 1 μ t + μ cos 4 1 + λ + μ 1 μ t , λ = 1 μ + 2 m + 7 μ 2 m μ μ 2 + m μ 2 . μ = 40 m 2 27 m 2 128 m + 128 2409 m 4 29600 m 3 + 111520 m 2 163840 m + 81920 .
Then, the approximate trigonometric solution to I.V.P. (27) reads
u trigo ( ϕ ) = 3 2 γ 1 P 1 A + 1 A 3 2 γ + 2 P cos m 2 ( ω ϕ ) . m = 2 G M ( A P ) 2 G M ( 2 A + P ) A c 2 P , ω = 1 2 P 2 γ P ( 1 + m ) , γ = 3 G M c 2 .
Solution (49) is periodic with period
T trigo = π ω 1 μ 1 + λ + μ
For the Mercury data (see Table 2), we have
Δ G R T ( T trigo 2 π ) · 23668612128 π · Siederal = 42.98145640188336 arcsec / cy
The trigonometric solution is highly accurate.

5. Contribution of the Cosmological Constant

Let us consider I.V.P. (28). The equation may be written in the form
u ( ϕ ) α + β u ( ϕ ) 3 + u ( ϕ ) γ u ( ϕ ) 2 = 0 , u ( 0 ) = 1 / P   and   u ( 0 ) = 0 . α = G M L 2 , β = c 2 Λ 3 L 2 , γ = 3 G M c 2 .
Next, we multiply the equation u ( ϕ ) α + β u ( ϕ ) 3 + u ( ϕ ) γ u ( ϕ ) 2 = 0 by u ( ϕ ) and then we integrate it with regard to ϕ to obtain
1 2 u ( ϕ ) 2 α u ( ϕ ) β 2 u ( ϕ ) 2 1 3 γ u ( ϕ ) 3 + u ( ϕ ) 2 2 = C ,
C is the constant of integration. Letting ϕ = 0 and taking into account conditions u ( 0 ) = 1 / P and u ( 0 ) = 0 gives
C = γ 3 P 3 β P 2 2 + 1 2 P 2 α P .
On the other hand, let ϕ A be the angle for aphelion position. Then, u ( ϕ A ) = 1 / A and u ( ϕ A ) = 0 so that
C = γ 3 P 3 β P 2 2 + 1 2 P 2 α P = γ 3 A 3 A 2 β 2 + 1 2 A 2 α A .
From (55),
α = G M L 2 = 3 A P ( A + P ) A 2 β P 2 + 1 2 γ A 2 + A P + P 2 6 A 2 P 2 .
The angular momentum of the planet now depends on the contribution of the cosmological constant:
L 2 = A 2 P 2 6 G M A c 2 Λ P ( A + P ) 3 A P ( A + P ) 2 γ A 2 + A P + P 2 .
From (53)–(56),
d u d ϕ 2 = 2 γ u 1 A 1 P u u 3 + 3 A P 2 A γ 2 P γ 2 A P γ u 2 3 A P ( A + P ) β 2 γ u 3 A P β 2 γ 3 u 2
The solution to (58) is a periodic function with period
T Λ = 2 3 2 γ 1 / A 1 / P u d u u 1 A 1 P u u 3 + 3 A P 2 A γ 2 P γ 2 A P γ u 2 3 A P ( A + P ) β 2 γ u 3 A P β 2 γ .

6. Estimation of the Value of the Cosmological Constant

6.1. First Approach

Our aim is to find a value for Λ , such that the value of T Λ is as close as possible to value T given by (40). To this end, assuming a positive cosmological constant, we obtain a range for Λ , say 10 k Λ 10 m . Then, we evaluate (59) for m = 1 , 2 , . . . , k . We obtain a list of pairs ( j , T 10 j ) that allows us to construct an interpolating function object. This function is denoted with Ψ . The next step consists of minimizing the quantity ( Ψ ( x ) T ) 2 on the interval 10 k x 10 m . Let k = 56 and m = 30 , so that 10 56 Λ 10 30 . For Mercury data our calculations produced the optimal value Λ + = 4.42355 × 10 44 . On the other hand, assuming a negative cosmological constant value, we obtained the optimal value Λ = 7.41311 × 10 42 ; see Table 3.
The confidence interval for Table 4 is:
Λ Λ + Λ Δ Λ { 1.4 × 10 37 , 2.3 × 10 37 } { 1.4 × 10 37 , 2.3 × 10 37 }

6.2. Second Approach

Integral (59) is hard to evaluate, so we give an approximate analytical expression. Assume a positive cosmological constant Λ = Λ + . We have
u 3 + 3 A P 2 A γ 2 P γ 2 A P γ u 2 3 A P ( A + P ) β 2 γ u 3 A P β 2 γ < u 3 + 3 A P 2 A γ 2 P γ 2 A P γ u 2 3 A P ( A + P ) β 2 γ u < u 3 + 3 A P 2 A γ 2 P γ 2 A P γ u 2 3 A P ( A + P ) β 2 γ u 2 = u 2 ( δ 1 u ) , where δ 1 = 3 A P 2 A γ 2 P γ 2 A P γ 3 A P ( A + P ) β 2 γ . Then u u 3 + 3 A P 2 A γ 2 P γ 2 A P γ u 2 3 A P ( A + P ) β 2 γ u 3 A P β 2 γ > 1 δ 1 u .
On the other hand, since u 1 / A , u 1 A , u 1 A and u 2 A 2 , and then
u 3 + 3 A P 2 A γ 2 P γ 2 A P γ u 2 3 A P ( A + P ) β 2 γ u 3 A P β 2 γ = u 2 u + 3 A P 2 A γ 2 P γ 2 A P γ 3 A P ( A + P ) β 2 γ u 1 3 A P β 2 γ u 2 u 2 u + 3 A P 2 A γ 2 P γ 2 A P γ 3 A P ( A + P ) β 2 γ A 3 A P β 2 γ A 2 = u 2 δ 2 u , where δ 2 = 3 A P 2 A γ 2 P γ 2 A P γ 3 A P ( A + P ) β 2 γ A 3 A P β 2 γ A 2 . Then u u 3 + 3 A P 2 A γ 2 P γ 2 A P γ u 2 3 A P ( A + P ) β 2 γ u 3 A P β 2 γ 1 δ 2 u .
From (60) and (61), we obtain
1 δ 1 u u u 3 + 3 A P 2 A γ 2 P γ 2 A P γ u 2 3 A P ( A + P ) β 2 γ u 3 A P β 2 γ < 1 δ 2 u ,
Then, taking into account (59), we have estimates
F ( δ 1 ) T Λ + < F ( δ 2 ) ,
where
F ( δ ) : = 2 3 2 γ 1 / A 1 / P d u u 1 A 1 P u ( δ u ) = 2 6 A K A P P ( A δ 1 ) γ A δ 1 .
Then,
T Λ + 1 2 ( F ( δ 1 ) + F ( δ 2 ) ) .
Equating (40) and (65), we obtain the required equation to determine Λ + .
On the other hand, suppose that Λ = Λ is negative. Then, since u > u 2 ,
u 3 + 3 A P 2 A γ 2 P γ 2 A P γ u 2 3 A P ( A + P ) β 2 γ u 3 A P β 2 γ > u 3 + 3 A P 2 A γ 2 P γ 2 A P γ u 2 3 A P ( A + P ) β 2 γ u > u 3 + 3 A P 2 A γ 2 P γ 2 A P γ u 2 3 A P ( A + P ) β 2 γ u 2 = u 2 ( 3 A P 2 A γ 2 P γ 2 A P γ 3 A P ( A + P ) β 2 γ u ) = u 2 ( δ 4 u ) , δ 4 = 3 A P 2 A γ 2 P γ 2 A P γ 3 A P ( A + P ) β 2 γ .
so that
u u 3 + 3 A P 2 A γ 2 P γ 2 A P γ u 2 3 A P ( A + P ) β 2 γ u 3 A P β 2 γ < 1 δ 4 u
Now, since u 1 A ,
u 3 + 3 A P 2 A γ 2 P γ 2 A P γ u 2 3 A P ( A + P ) β 2 γ u 3 A P β 2 γ = u 2 u + 3 A P 2 A γ 2 P γ 2 A P γ 3 A P ( A + P ) β 2 γ u 1 3 A P β 2 γ u 2 u 2 u + 3 A P 2 A γ 2 P γ 2 A P γ 3 A P ( A + P ) β 2 γ A 3 A P β 2 γ A 2 = u 2 ( δ 3 u ) , w h e r e δ 3 = 3 A P 2 A γ 2 P γ 2 A P γ 3 A P ( A + P ) β 2 γ A 3 A P β 2 γ A 2 .
We obtained the following estimates:
1 δ 3 u u u 3 + 3 A P 2 A γ 2 P γ 2 A P γ u 2 3 A P ( A + P ) β 2 γ u 3 A P β 2 γ
Using (68) and (69), and taking into account (59)–(64) we have:
F ( δ 3 ) T Λ < F ( δ 4 ) .
The equation for Λ is 1 2 ( F ( δ 3 ) + F ( δ 4 ) ) = T , where T is found from (40).

7. Approximate Analytical Solution to the Motion Equation with a Nonzero Cosmological Constant

7.1. First Solution Method

Since Λ is small,
u 3 + 3 A P 2 A γ 2 P γ 2 A P γ u 2 3 A P ( A + P ) β 2 γ u 3 A P β 2 γ u 3 + 3 A P 2 A γ 2 P γ 2 A P γ u 2 3 A P ( A + P ) β 2 γ u 2 .
Then, Hard Ode (58) may be approximated by the ode
d u d ϕ 2 = 2 γ 3 u 1 A 1 P u δ u , where δ = 3 A P 2 A γ 2 P γ 2 A P γ 3 A P ( A + P ) β 2 γ , β = c 2 Λ 3 L 2 . L 2 = A 2 P 2 6 G M A c 2 Λ P ( A + P ) 3 A P ( A + P ) 2 γ A 2 + A P + P 2 .
The exact solution to Ode (72) given initial conditions u ( 0 ) = 1 / P and u ( 0 ) = 0 is given by
u ( ϕ ) = δ + 1 P δ nd γ ( A δ 1 ) 6 A ϕ , A P P ( A δ 1 ) 2
Let us compare this solution with the numerical solution to equivalent Ode (28) for Mercury data assuming the value of Λ = 10 46 , see Figure 2.
Our semianalytical solution was highly accurate and involved the cosmological constant. This also offers another way to estimate the value of the cosmological constant. As a matter of fact, the period of the semianalytical solution reads
T Λ = 2 K A P P ( A δ 1 ) γ ( A δ 1 ) 6 A
Then, the equation for Λ is
T Λ = 2 6 A γ ( A δ 1 ) K A P P ( A δ 1 ) = T : = 4 K 2 G M ( A P ) A c 2 P 2 G M ( A + 2 P ) 1 2 G M ( A + 2 P ) A c 2 P
Since the value of m = A P P ( A δ 1 ) is small, we may use the following approximation:
2 K ( m ) 1 5 m 16 1 9 m 16 π .
Then,
T Λ = 6 π A ( A ( 16 δ P 5 ) 11 P ) γ A δ 1 ( A ( 16 δ P 9 ) 7 P ) = T ,
from where
256 A 3 P 2 T 2 γ δ 3 96 A 2 P 16 A P π 2 + 3 A T 2 γ + 5 P T 2 γ δ + 3 A 320 A 2 P π 2 + 704 A P 2 π 2 + 27 A 2 T 2 γ + 138 A P T 2 γ + 91 P 2 T 2 γ δ 2 150 A 3 π 2 660 A 2 P π 2 726 A P 2 π 2 81 A 2 T 2 γ 126 A P T 2 γ 49 P 2 T 2 γ = 0 .
Solving this cubic, we obtain the value for δ . Then, from (72),
Λ = 6 μ Q 2 ( 2 γ δ P 3 P + 2 γ ( Q + 1 ) ) c 2 P 2 ( Q + 1 ) 2 γ + 2 γ δ P 2 + 2 γ ( Q + 1 ) ( P + Q ) 3 P ( P + Q + 1 ) , Q : = P A .
We may go further using other approximates for K ( m ) . For example,
K ( m ) π 409 m 2 3984 m + 4864 1025 m 2 5200 m + 4864
This approximation gives us the following quintic for δ :
C 0 δ 5 + C 1 δ 4 + C 2 δ 3 + C 3 δ 2 + C 4 δ + C 5 = 0 ,
where
C 0 = 23658496 A 5 γ P 4 T 2 .
C 1 = 155648 A 4 P 3 912 π 2 A P + 325 A γ T 2 + 435 γ P T 2 .
C 2 = 256 A 3 P 2 908352 π 2 A 2 P + 144575 A 2 γ T 2 + 1309632 π 2 A P 2 + 501250 A γ P T 2 + 278335 γ P 2 T 2 .
C 3 = 32 A 2 P 3722064 π 2 A 3 P + 333125 A 3 γ T 2 + 14356320 π 2 A 2 P 2 + 2470425 A 2 γ P T 2 + 8537424 π 2 A P 3 + 3544575 A γ P 2 T 2 + 1045155 γ P 3 T 2 .
C 4 = A 19553472 π 2 A 4 P + 1050625 A 4 γ T 2 + 179551680 π 2 A 3 P 2 + 17117500 A 3 γ P T 2 + 279850560 π 2 A 2 P 3 + 53377350 A 2 γ P 2 T 2 + 88848192 π 2 A P 4 + 40032700 A γ P 3 T 2 + 6714305 γ P 4 T 2 .
C 5 = 1003686 π 2 A 5 15538728 π 2 A 4 P 1050625 A 4 γ T 2 66467748 π 2 A 3 P 2 6457500 A 3 γ P T 2 48971688 π 2 A 2 P 3 11334950 A 2 γ P 2 T 2 9969126 π 2 A P 4 4340700 A γ P 3 T 2 474721 γ P 4 T 2 .
Solving the quintic, we obtain one or more estimates for Λ . Let us evaluate the Λ using Formula (79) for Mercury. The cubic resolvent reads
3.224134057504410567808405 × 10 61 δ 3 + 1.0917240150709210311767171311214 × 10 58 δ 2 3.6333609324359012885422051569104 × 10 47 δ + 3.0230424838943767016875726381236 × 10 36 = 0 .
Its roots are:
δ 1 = 1.66403214678195905033219994 × 10 11 . δ 2 = 1.66406282054722562545891615 × 10 11 . δ 3 = 0.0003386099610925325989588409 .
The values of Λ from (79) are:
Λ 1 = 2.3818552552068309770172818 × 10 29 .
Λ 2 = 2.3818552552068309770172040 × 10 29
Λ 3 = 2.574519419707803 × 10 43
Using the quintic, we obtain the negative value
Λ = 1.0241093318102231833808 × 10 47
The algebraic approach may give either positive or negative estimates for the cosmological constant. Letting Λ = 2.5 × 10 43 , we obtain the following polar plot (Figure 3) for Mercury’s trajectory using the numerical solution to Ode (28).
After many calculations to obtain Λ estimates, we found the following interesting fact. Number δ in Equation (72) was very close to δ = 1 / 2953 for all planets in the solar system. The mean of these values was
δ = 0.000338609988014895212199968489565 .
This implies that
6 K A P P ( A δ 1 ) γ δ 1 A = 2 K 2 ( A P ) μ 2 ( 2 A + P ) μ A c 2 P 1 2 μ ( 2 A + P ) A c 2 P , μ = G M .
That is, the planets are related to Curve (93) in the A P plane. This allows for us to predict the value of the aphelion given the the perihelion and vice versa. Indeed, if we are given the aphelion, we simply solve transcendental Equation (93) for the perihelion. On the other hand, the fact that δ is almost a constant tells us that Λ defined by (79) must indeed be a constant !. This justifies the name ’cosmological constant’; see Figure 4.

7.2. Second Solution Method

Let
F ( u ) = u 3 + 3 A P 2 A γ 2 P γ 2 A P γ u 2 3 A P ( A + P ) β 2 γ u 3 A P β 2 γ , β = c 2 Λ 3 L 2 .
Let us introduce the following notations:
d 0 = 3 A P β 2 γ , d 1 = 3 A P ( A + P ) β 2 γ , d 2 = 3 A P 2 A γ 2 P γ 2 A P γ
For a positive Λ , the cubic d 0 + d 1 u + d 2 u 2 u 3 = 0 has two small real roots of opposite signs. For a negative Λ , we have only one positive real root, and the two other roots are imaginary and small in magnitude. This allows for us to perform the approximation
F ( u ) u 2 d 2 4 + 3 d 1 d 2 2 + 2 d 0 d 2 + d 1 2 d 2 3 + 2 d 1 d 2 + d 0 u .
Then, hard Ode (58) is replaced with easy ode
d u d ϕ 2 = 2 γ 3 u 1 A 1 P u d 2 4 + 3 d 1 d 2 2 + 2 d 0 d 2 + d 1 2 d 2 3 + 2 d 1 d 2 + d 0 u
so that, in view of (64),
T Λ 2 3 2 γ 1 / A 1 / P d u u 1 A 1 P u ( δ u ) , δ = d 2 4 + 3 d 1 d 2 2 + 2 d 0 d 2 + d 1 2 d 2 3 + 2 d 1 d 2 + d 0
The required equation to estimate the Λ is then
2 6 A γ ( A δ 1 ) K A P P ( A δ 1 ) = 4 K 2 G M ( A P ) A c 2 P 2 G M ( A + 2 P ) 1 2 G M ( A + 2 P ) A c 2 P .
For Mercury data, the estimated value using (99) is Λ + = 6.51 × 10 52 . For Pluto, Λ = 2.83 × 10 52 . For data [25] P = 46005766800 , A = 69818681600 , the estimated value using Formula (99) is Λ = 3.52 × 10 51 .

7.3. Third Approach

We may proceed quite differently in order to estimate the Λ value. Let u app ( ϕ ) be some approximate solution to I.V.P. (27), for example, u = u trigo . Then, replacing u with u app in (58) gives
Λ ( ϕ ) = f ( ϕ ) : = L 0 2 u app ( ϕ ) 2 3 A P A u app ( ϕ ) 1 P u app ( ϕ ) 1 A P u app ( ϕ ) 2 2 γ A u app ( ϕ ) 1 P u app ( ϕ ) 1 A P u app ( ϕ ) + A + P A 2 c 2 P 2 A u app ( ϕ ) 1 P u app ( ϕ ) 1 ( A + P ) u app ( ϕ ) + 1 .
Then, we take
Λ integral = Λ ( A , P ) = 1 T 0 T 0 / 2 T 0 / 2 f ( ϕ ) d ϕ .
See Table 5 for the different Λ estimates.
Predicted value for Λ + :
Arithmetic   Mean : Λ + = 2.87 × 10 44 . Geometric   Mean : Λ + = 1.25 × 10 52
Predicted value for Λ :
Arithmetic   Mean : Λ = 7.25 × 10 44 . Geometric   Mean : Λ = 1.01 × 10 52
Confidence interval for Table 6:
Λ Λ + Λ Δ Λ { 1.3 × 10 43 , 1.9 × 10 43 } { 4.8 × 10 43 , 3.3 × 10 43 }

8. Analysis and Discussion

We obtained the estimated values for the cosmological constant, and we compared our results with that obtained in [25]. The nonlinear differential equation in [25] taking into account the cosmological constant reads
d u d ϕ 2 = 2 G M c 2 u ( ϕ ) 2 u ( ϕ ) 1 A 1 P u ( ϕ ) d 0 + d 1 u ( ϕ ) + d 2 u ( ϕ ) 2 u 3 ( ϕ ) ,
where
d 0 = c 2 A 2 c 2 P 2 G M + A P c 2 P 2 G M 2 G M P 2 2 A G M P A c 2 Λ P ( A + P ) + 6 G M Λ . d 1 = c 2 ( A + P ) A 2 c 2 P 2 G M + A P c 2 P 2 G M 2 G M P 2 2 A G M P A c 2 Λ P ( A + P ) + 6 G M Λ . d 2 = c 2 2 G M 1 A 1 P .
In the case when Λ = 0 the exact solution is
u ( ϕ ) = c 2 6 G M + 2 c 2 G M t + t 0 ; g 2 , g 2 ,
where
g 2 = g 2 = A 2 c 4 P 2 6 A 2 c 2 G M P + 12 A 2 G 2 M 2 6 A c 2 G M P 2 + 12 A G 2 M 2 P + 12 G 2 M 2 P 2 12 A 2 c 4 P 2 . g 3 = g 3 = A c 2 6 G M c 2 P 6 G M A c 2 P 3 A G M 3 G M P 216 A 2 c 6 P 2 . t 0 = 1 G M 2 c 2 1 P c 2 6 G M ; g 2 , g 2
The solution is periodic, and its period is given by
T = 4 K 2 G M ( A P ) A c 2 P 2 G M ( A + 2 P ) 1 2 G M ( A + 2 P ) A c 2 P .
In [25] authors use the following data:
P = 45997620600 , A = 69820729300 .
The perihelion shift for these data is 42.9825 arc-sec/century. The value in [25] was =42.9817 arc-sec/century. On the other hand, for the data in [25],
P = 46001260500 , A = 69817089400 .
The perihelion shift for these data is 42.9814 arc-sec/century. The value in [25] is 42.9805 arc-sec/century. The authors in [25] used the data
P = 46010448900 , A = 69807901000 .
For these data, the perihelion shift is 42.9784 arc-sec/century. The value in [25] is =42.9776 arc-sec/century.
Using Mercury data in Table 2, we obtain a perihelion shift of 42.9815 arc-sec/century.
In the case when Λ 0 , cubic d 0 + d 1 u + d 2 u 2 u 3 = 0 has two roots very close to zero. Then, we may use the approximation
d 2 u 2 + d 1 u + d 0 u 3 d 0 d 2 2 + d 1 d 2 3 + 2 d 1 d 2 + d 0 + d 1 2 + d 2 2 d 1 + d 0 d 2 d 2 3 + 2 d 1 d 2 + d 0 u + u 2 d 2 4 + 3 d 1 d 2 2 + 2 d 0 d 2 + d 1 2 d 2 3 + 2 d 1 d 2 + d 0 u
For example, let Λ = 10 56 , P = 46005766800 , A = 69818681600 as in [25]. We have:
d 0 + d 1 u + d 2 u 2 u 3 = ( u 3.6156749990681577 × 10 25 ) ( u + 3.6156749990680067 × 10 25 ) ( 0.0003386099606939169 u ) .
On the other hand, the right-hand side of (109) is written as
( u 3.6156749990681577 × 10 25 ) ( u + 3.6156749990680067 × 10 25 ) ( 0.00033860996069391686 u )
Using these facts and taking into account that
1 u 2 d 0 d 2 2 + d 1 d 2 3 + 2 d 1 d 2 + d 0 + d 1 2 + d 2 2 d 1 + d 0 d 2 d 2 3 + 2 d 1 d 2 + d 0 u + u 2 1 A 2 d 0 d 2 2 + d 1 d 2 3 + 2 d 1 d 2 + d 0 + d 1 2 + d 2 2 d 1 + d 0 d 2 d 2 3 + 2 d 1 d 2 + d 0 · 1 A + 1 A 2 : = N
hard Ode (104) may be replaced with easy ode
d u d ϕ 2 = D u ( ϕ ) 1 A 1 P u ( ϕ ) d 2 4 + 3 d 1 d 2 2 + 2 d 0 d 2 + d 1 2 d 2 3 + 2 d 1 d 2 + d 0 u ( ϕ ) .
The exact solution to Ode (113) is
u ( ϕ ) = 1 P + 3 ( P A ) ( δ P 1 ) P ( A δ P 2 A + P ) 1 + 12 A P D ( A δ P 2 A + P ) ( t ; g 2 , g 3 ) ,
where
D = 2 G M N c 2 , δ = d 2 4 + 3 d 1 d 2 2 + 2 d 0 d 2 + d 1 2 d 2 3 + 2 d 1 d 2 + d 0 . g 2 = D 2 A 2 δ 2 P 2 A 2 δ P + A 2 A δ P 2 A P + P 2 12 A 2 P 2 . g 3 = D 3 ( A δ P + A 2 P ) ( A δ P 2 A + P ) ( 2 A δ P A P ) 432 A 3 P 3 .
Solution (113) is periodic, and its period equals
T Λ = 2 r 1 4 x 3 g 2 x g 3 = 2 P D ( δ P 1 ) K A P A ( P δ 1 ) ,
where r is the greatest real root to cubic 4 x 3 g 2 x g 3 = 0 . Observe that
4 x 3 g 2 x g 3 = x + D ( A ( 1 + P δ ) 2 P ) 12 A P x + D ( P ( 1 + A δ ) 2 A ) 12 A P x D ( 2 A P δ A P ) 12 A P .
The required equation to determine Λ is
2 P D ( δ P 1 ) K A P A ( 1 P δ ) = T : = 4 K 2 G M ( A P ) 2 G M ( 2 A + P ) A c 2 P 1 2 G M ( 2 A + P ) A c 2 P
The authors in [25] gave the predicted values Λ = 10 55 and Λ + = 10 56 for the data
P = 46005766800 , A = 69818681600 .
On the other hand, the predicted values for Λ solving (118) for data (119) are
Λ + = 3.0849194019988974 × 10 45 and Λ = 1.8381582529562642 × 10 45
Now, using the method in previous sections, we obtained the following estimates for Data (119) in [25]
Λ + = 1.024 × 10 47 and Λ = 2.58 × 10 43
Lastly, for the other planets, the estimates for Λ by solving (118) are summarized in Table 7:
Confidence interval for Table 7:
Λ Λ + Λ Δ Λ { 1.1 × 10 37 , 1.6 × 10 37 } { 5.5 × 10 32 , 3.8 × 10 32 }
The different previous approaches show that the cosmological constant takes different values, both positive and negative, which leads to establishing that the static cosmological model is unstable and antisymmetric. Due to the difference in the perihelion for each planet, the cosmological constant is different for each one, which implies that the trajectory of the planets must be antisymmetric, that is, elliptical.
The Schwarzschild radius for the Sun is r S = 2 G M c 2 = 2.7 km. Taking into account the cosmological constant, the Schwarzschild radius is obtained by solving the following cubic equation for r.
1 1 3 r 2 Λ 2 G M c 2 r = 0 .
Solving (122), we obtain that r S Λ = 4.1 km. Both radii are inside the Sun.
The region of spacetime due to the gravitational field of the Sun is obtained from
1 1 3 r 2 Λ 2 G M c 2 r = 1 ,
from where r lim = 5.1 × 10 10 Km = 340 UA.
Between the planet Pluto and the edge of the curved space, an area of the curved spacetime region originates that allows for, in Schwarzschild’s cosmological model, predicting the existence of new planets (planetoids X). Among those catalogued planets are: Kuiper, whose distance to the Sun is 30 AU; Quaoar, 43.4 AU; Makameke, 45.8 AU; Sedna, 88.5 AU; Eris, 98.3 AU; and V774104, 103 AU.
This limit was calculated for the cosmological constant obtained for Mercury Λ = 6.56145 × 10 38 m 2 .
After the boundary of the curved spacetime, a region of Minkowski plane spacetime originates.
From Schwarzschild’s stationary cosmological model (Figure 5), it follows that the planets are in the region of spacetime curved due to the Sun; therefore, each planet follows a stationary geodesic trajectory. All of the above are represented in the following graph.
The static cosmological model of Schwarzschild, the cosmological constant according to the results obtained in Table 4 and Table 6, which were obtained by different approaches, shows that, in the planetary system, the cosmological constant is unstable. In order to understand this effect of instability in the static cosmological model, it is necessary to take into account the effects of the acceleration due to the Big Bang until the current static state is reached.
Lastly, the instability effect of the cosmological constant leads to the fluctuation of the limit radius of the curved spacetime, affecting the stability of the curvature region.

9. Conclusions

In this work, the equation of Einstein’s general theory of relativity for vacuum was solved. Taking into account the cosmological constant, the nonlinear differential equation that describes the movement of the planets was constructed and exactly solved.
Solving the inverse problem, different theoretical estimates were obtained to calculate the value of the cosmological constant. The obtained results were compared with those of other authors, both theoretical and experimental. In the Schwarzschild cosmological radius r S Λ = 4.1 km, and the radius limit of the curvature of the gravitational field of the solar system r lim = 5.1 × 10 10 Km = 340 UA were obtained, and the existence of several regions was discovered, which were classified into: planetary curvature region, exoplanet curvature region, and the prediction region of the existence of new exoplanets.
The proposed methodology can be of great interest to astronomers, cosmologists, nonlinear physics researchers, and all those interested in the study of the universe.

Author Contributions

Conceptualization: J.E.P.Q., J.E.C.H. methodology: J.E.P.Q., J.E.C.H., A.H.S.S. Software: A.H.S.S. Formal analysis: J.E.P.Q., J.E.C.H., A.H.S.S. writing—revision and edition: J.E.P.Q., J.E.C.H., A.H.S.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

Not applicable.

Acknowledgments

The authors thank the University of Francisco Jose de Caldas, University of Tolima, and National university of Colombia for the support to carry out this work.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Einstein, A. Die Feldgeichungen der Gravitation. Sitzungsber. Preuss. Akad. Wiss. 1915, 844–847. Available online: https://ui.adsabs.harvard.edu/abs/1915SPAW.......844E/abstract (accessed on 11 July 2022).
  2. Einstein, A. Die Grundlage der allgeminen Relativitatstheorie. Ann. Phys. 1916, 49, 769. [Google Scholar] [CrossRef]
  3. Schwarzschild, K. Ueber das Gravitattionsfeld eines Massenpunktes nach der Einsteinschen theorie. Preuss. Ak. Wiss. Sitzungsber 1916, 189–196. [Google Scholar]
  4. Urbain J -J Le verrier. Annales De L´observatoire de paris; Nabu Press: Warsaw, Poland, 1859. [Google Scholar]
  5. Eddington, A.S. The total eclipse of 1919 May 29 and the influence of gravitation on light. Observatory 1919, 42, 119–122. [Google Scholar]
  6. Kasner, K. Geometrical theorems on Einstein’s cosmological equations. Am. J. Math. 1921, 43, 217–221. [Google Scholar] [CrossRef]
  7. Pinzón, J.E.; Salas, A.H.; Castillo, J.E. Non-static solutions of the Einstein equation in the vacuum. Int. J. Math. Comput. Sci. 2022, 17, 439–445. [Google Scholar]
  8. Pinson, K.H.E. Kasner solution and the causality principle in the relativistic theory of gravitation. Theor. Mat. Phys. 1991, 89, 1236–1238. [Google Scholar] [CrossRef]
  9. Pinson, K.H.E. Nonstatic solutions of Kazner type in the relativistic theory of gravitation. Theor. Mat. Phys. 1992, 91, 449–451. [Google Scholar] [CrossRef]
  10. Robertson, H.P. Relativistic Cosmology. Rev. Mod. Phys. 1933, 5, 62–90. [Google Scholar] [CrossRef]
  11. Walker, A.G. On Milne´s Theory of world structure. Proc. Lond. Math. Soc. 1936, 42, 90–127. [Google Scholar]
  12. Friedmann, A. Uber die Krummung des Raumes. Zeitschrift für Physik 1922, 10, 377–386. [Google Scholar] [CrossRef]
  13. Georges, L. Un univers homogène de masse constante et de rayon croissant rendant compte de la vitesse radiale des nébuleuses. Ann. Soc. Sci. Bruxelles 1927, A47, 49–59. [Google Scholar]
  14. Einstein, A. Cosmological Considerations in the General Theory of Relativity. Sitzungsber. Preuss. Akad. Wiss 1917, 1, 142–152. [Google Scholar]
  15. De Sitter, W. On Einstein’s theory of gravitation and its astronomical consequences: Third paper. Mon. Not. R. Astron. Soc. 1917, 78, 3–28. [Google Scholar] [CrossRef]
  16. Edwin, H. A Relation Betwewn Distance and radial velocity Among Extra-Galatic Nebulae. Proc. Natl. Acad. Sci. USA 1929, 15, 168–173. [Google Scholar]
  17. Eddington, A.S. The Expanding Universe; Cambridge University Press: Cambridge, UK, 1993. [Google Scholar]
  18. Eddington, A.S. On the Instability of Einstein’s Spherical World. Mon. Not. R. Astron. Soc. 1930, 90, 668–678. [Google Scholar] [CrossRef]
  19. Einstein, A. Prinzipielles zur allgemeinen Relativitätstheorie. Ann. Phys. 1918, 55, 241–244. [Google Scholar] [CrossRef]
  20. Tartaglia, A. Testing General Relativity. arXiv 2015, arXiv:1503.00920. [Google Scholar]
  21. Wilhelm, K.; Dwivedi, B.N. Secular perihelion advances of the inner planets and asteroid Icarus. New Astron. 2014, 31, 51–55. [Google Scholar] [CrossRef]
  22. Salas Salas, A.H.; Castillo Hernandez, J.E.; Restrepo R., C.E. The Jacobi elliptic functions and their applications in the advance of mercury’s perihelion. Rev. Mex. Fis. 2015, 61, 392–399. [Google Scholar]
  23. Roig, F.; Nesvorný, D.; DeSouza, S.R. Jumping jupiter can explain mercury’s orbit. Astrophys. J. Lett. 2016, 820, L30. [Google Scholar] [CrossRef]
  24. Lorenzo, I. The recently determined anomalous perihelion precession of saturn. Astron. J. 2009, 137, 3615–3618. [Google Scholar]
  25. Kraniotis, G.V.; Whitehouse, S.B. Compact calculation of the perihelion precession of Mercury in general relativity, the cosmological constant and Jacobi’s inversion problem. Class. Quantum Gravity 2003, 20, 22. [Google Scholar] [CrossRef]
  26. Gibbons, G.W.; Hawking, S.W. Cosmological event horizons, thermodynamics, and particle creation. Phys. Rev. D 1977, 15, 2378. [Google Scholar] [CrossRef] [Green Version]
  27. Available online: https://calgary.rasc.ca/orbits.htm (accessed on 11 July 2022).
Figure 1. Advance of Mercury’s perihelion.
Figure 1. Advance of Mercury’s perihelion.
Universe 08 00449 g001
Figure 2. Comparison between semianalytical (73) and numerical solutions to the ode.
Figure 2. Comparison between semianalytical (73) and numerical solutions to the ode.
Universe 08 00449 g002
Figure 3. Mercury’s trajectory using Ode (28) with Λ = Λ + = 2.57 × 10 43 .
Figure 3. Mercury’s trajectory using Ode (28) with Λ = Λ + = 2.57 × 10 43 .
Universe 08 00449 g003
Figure 4. Curve (93).
Figure 4. Curve (93).
Universe 08 00449 g004
Figure 5. Schwarzschild’s cosmological model.
Figure 5. Schwarzschild’s cosmological model.
Universe 08 00449 g005
Table 1. Exact and approximate values for perihelion shifts.
Table 1. Exact and approximate values for perihelion shifts.
Planet Δ GTR (Exact Value) Δ ^ GTR Δ GTR Δ ˜ GTR
Mercury 42.981518799592365 42.98151983984209 1.04 × 10 6
Venus 8.625078587720143 8.625078762663101 1.75 × 10 7
Earth 3.8387718031722122 3.8387718587113353 5.55 × 10 8
Mars 1.3508767421610324 1.3508767534274582 1.13 × 10 8
Jupiter 0.062311625478272306 0.06231162564333644 1.65 × 10 10
Saturn 0.013688932210705418 0.013688932230249766 1.95 × 10 11
Uranus 0.0023848021538887348 0.002384802155606265 1.72 × 10 12
Neptune 0.000774085957956637 0.0007740859583331821 3.77 × 10 13
Pluto 0.0004175446900080054 0.0004175446900902417 8.22 × 10 14
Table 2. Planets data. From [27].
Table 2. Planets data. From [27].
PlanetA:Apelion (m)P:Perihelion (m)Siederal Period Q = P A
Mercury6981,7332,00046,000,870,000 87.969089 0.658875
Venus108,939,198,000107,475,372,000 224.69562 0.986563
Earth152,100,915,000147,094,882,000 365.256622 0.967087
Mars249,226,166,000206,662,107,000 686.99329 0.829215
Jupiter816,054,481,000740,505,444,000 4334.24677 0.907422
Saturn1,506,619,721,0001,348,156,111,00010,765.21936 0.894822
Uranus3,004,984,160,0002,735,977,617,00030,700.24558 0.91048
Neptune4,538,617,181,0004,458,057,447,00060,226.53638 0.98225
Pluto7,377,158,662,0007,437,141,859,00090,631.02406 0.60147
Table 3. NASA data.
Table 3. NASA data.
Universe 08 00449 i001A:Apelion (m)P:Perihelion (m)Siederal Period Δ ^ GTR
Mercury69.817 × 10 9 46.002 × 10 9 87.96842.9814956752253
Venus107.476 × 10 9 108.939 × 10 9 224.7018.62485450672991
Earth147.092 × 10 9 152.099 × 10 9 365.2423.838987490384246
Mars249.229 × 10 9 206.617 × 10 9 686.9801.351057133139973
Jupiter816.618 × 10 9 740.522 × 10 9 4332.5890.062314275340345
Saturn1514.504 × 10 9 1352.555 × 10 9 10,759.220.012639384845858
Table 4. Estimation of the cosmological constant using the first approach.
Table 4. Estimation of the cosmological constant using the first approach.
Planet Λ + Λ
Mercury 6.56145 × 10 38 7.94328 × 10 45
Venus 3.03249 × 10 37 10 46
Earth 4.22357 × 10 38 1.2 × 10 46
Mars 3.8824 × 10 47 3.98 × 10 39
Jupiter 2.02148 × 10 40 1.2 × 10 48
Saturn 3.11889 × 10 41 1.2 × 10 49
Uranus 3.98107 × 10 42 10 50
Neptune 1.58489 × 10 42 1.26 × 10 51
Pluto 5.12861 × 10 44 2.5 × 10 51
Table 5. δ value in Equation (72).
Table 5. δ value in Equation (72).
Mercury δ = 0.000338609961092532513876607325187
Venus δ = 0.000338609978670453457775396666207
Earth δ = 0.000338609983781427330631896843727
Mars δ = 0.000338609988303106856661478518689
Jupiter δ = 0.000338609994578505351780356447122
Saturn δ = 0.000338609995748851561048198277604
Uranus δ = 0.000338609996456062380901924635523
Neptune δ = 0.000338609996709698414518729947531
Pluto δ = 0.000338609996793419042605127744494
Table 6. Estimation of the cosmological constant using the second and third approaches.
Table 6. Estimation of the cosmological constant using the second and third approaches.
Planet Λ + Λ Λ integral
Mercury 2.58 × 10 43 1.02 × 10 47 4.18351 × 10 37
Venus 6.34 × 10 49 6.51 × 10 43 5.54028 × 10 38
Earth 1.05 × 10 49 4.26 × 10 45 1.45838 × 10 38
Mars 4.0 × 10 51 5.35 × 10 45 2.22584 × 10 39
Jupiter 3.33 × 10 53 2.46 × 10 55 2.62516 × 10 41
Saturn 5.08 × 10 55 3.45 × 10 60 2.53073 × 10 42
Uranus 2.8 × 10 55 5.33 × 10 59 1.42783 × 10 43
Neptune 4.13 × 10 56 2.00 × 10 59 2.13098 × 10 44
Pluto 1.38 × 10 56 8.303 × 10 60 3.78788 × 10 45
Table 7. Theoretical values for Λ .
Table 7. Theoretical values for Λ .
Planet Λ + Λ
Mercury 3.09 × 10 45 1.83 × 10 45
Venus 2.64 × 10 48 2.56 × 10 47
Earth 1.46 × 10 47 2.12 × 10 46
Mars 8.71 × 10 48 7.56 × 10 32
Jupiter 9.98 × 10 49 1.885 × 10 33
Saturn 8.01 × 10 50 2.74 × 10 51
Uranus 2.23 × 10 51 1.99 × 10 50
Neptune 2.18 × 10 37 9.1 × 10 51
Pluto 3.68 × 10 51 4.66 × 10 36
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Salas Salas, A.H.; Castillo Hernandez, J.E.; Pinzon Quintero, J.E. Calculation of the Cosmological Constant for the Planetary System in Schwarzschild’s Cosmological Model. Universe 2022, 8, 449. https://doi.org/10.3390/universe8090449

AMA Style

Salas Salas AH, Castillo Hernandez JE, Pinzon Quintero JE. Calculation of the Cosmological Constant for the Planetary System in Schwarzschild’s Cosmological Model. Universe. 2022; 8(9):449. https://doi.org/10.3390/universe8090449

Chicago/Turabian Style

Salas Salas, Alvaro Humberto, Jairo Ernesto Castillo Hernandez, and Jorge Enrique Pinzon Quintero. 2022. "Calculation of the Cosmological Constant for the Planetary System in Schwarzschild’s Cosmological Model" Universe 8, no. 9: 449. https://doi.org/10.3390/universe8090449

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