Next Article in Journal
Diagonals of Rational Functions: From Differential Algebra to Effective Algebraic Geometry
Previous Article in Journal
Study on the Preparation and Performance of a Prefabricated Pervious-Cement-Based Concrete Pavement Slab
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Stable Exponential Cosmological Type Solutions with Three Factor Spaces in EGB Model with a Λ-Term

by
Kubantai K. Ernazarov
1 and
Vladimir D. Ivashchuk
1,2,*
1
Institute of Gravitation and Cosmology, Peoples’ Friendship University of Russia (RUDN University), 6 Miklukho-Maklaya Street, 117198 Moscow, Russia
2
Center for Gravitation and Fundamental Metrology, VNIIMS, 46 Ozyornaya Street, 119361 Moscow, Russia
*
Author to whom correspondence should be addressed.
Symmetry 2022, 14(7), 1296; https://doi.org/10.3390/sym14071296
Submission received: 2 June 2022 / Revised: 17 June 2022 / Accepted: 20 June 2022 / Published: 22 June 2022
(This article belongs to the Section Physics)

Abstract

:
We study a D-dimensional Einstein–Gauss–Bonnet model which includes the Gauss–Bonnet term, the cosmological term Λ and two non-zero constants: α 1 and α 2 . Under imposing the metric to be diagonal one, we find cosmological type solutions with exponential dependence of three scale factors in a variable u, governed by three non-coinciding Hubble-like parameters: H 0 , h 1 and h 2 , obeying m H + k 1 h 1 + k 2 h 2 0 , corresponding to factor spaces of dimensions m > 1 , k 1 > 1 and k 2 > 1 , respectively, and depending upon sign parameter ε = ± 1 , where ε = 1 corresponds to cosmological case and ε = 1 —to static one). We deal with two cases: (i) m < k 1 < k 2 and (ii) 1 < k 1 = k 2 = k , k m . We show that in both cases the solutions exist if ε α = ε α 2 / α 1 > 0 and α Λ > 0 satisfy certain (upper and lower) bounds. The solutions are defined up to solutions of a certain polynomial master equation of order four (or less), which may be solved in radicals. In case (ii), explicit solutions are presented. In both cases we single out stable and non-stable solutions as u ± . The case H = 0 is also considered.

1. Introduction

In this semi-review article, which generalizes our previous work [1], we deal with the so-called Einstein–Gauss–Bonnet (EGB) gravitational model in dimensions D > 7 , which contains the Gauss–Bonnet term and the cosmological term Λ . The model also includes two non-zero constants: α 1 and α 2 , corresponding to Einstein and Gauss–Bonnet terms, respectively. It is well-known that the equations of motion for this model are of the second order (as it appears in General Relativity). The so-called Gauss–Bonnet term has appeared in (super)string theory as a second order correction in curvature to the effective (super)string effective action [2,3].
At present, the EGB gravitational model, e.g., with a cosmological term, and its modifications [4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25], are under intensive studies in cosmology and astrophysics, aimed at a solution of the dark energy problem, i.e., a possible explanation for the accelerating expansion of the Universe, which follows from supernovae (type Ia) observational data [26,27], and the search for a possible local manifestation of dark energy (related to black holes, wormholes etc.).
In this article we start with the so-called cosmological type solutions with “diagonal” metric d s 2 = ε ( d u ) 2 + i = 1 n ε i a i 2 ( u ) ( d y i ) 2 , governed by n > 3 scale factors ( D = n + 1 , ε = ± 1 , ε i = ± 1 ) depending upon one variable u, which is the synchronous time variable for the cosmological case, when ε = ε i = 1 . For the case ε = 1 and ε 1 = 1 , ε j = 1 ( j > 1 ) we get static configurations described by space-like variable (coordinate) u and time-like coordinate y 1 . In the cosmological case the equations of motion are governed by an effective Lagrangian which contains a 2-metric (or minisupermetric) G i j and a finslerian metric G i j k l , see Refs. [13,14] for Λ = 0 and Ref. [28] for Λ 0 .
Here we consider the cosmological type solutions with exponential dependence of scale factors (upon u-variable) and obtain a class of solutions with three scale factors, governed by three non-coinciding Hubble-like parameters—H, h 1 and h 2 —corresponding to factor spaces of dimensions m > 1 , k 1 > 1 and k 2 > 1 , respectively ( D = 1 + m + k 1 + k 2 ). Here we impose the following restriction S 1 = m H + k 1 h 1 + k 2 h 2 0 , excluding the solutions with constant volume factor and addressing a classification theorem which tells us that for generic anisotropic exponential solutions with Hubble-like parameters h 1 , , h n obeying S 1 = i = 1 n h i 0 the number of different (real) numbers among h 1 , , h n may be 1, or 2, or 3 [21]. The main goal of this paper is to extend the results of Ref. [1] to a class of cosmological type solutions, which include static ones (with ε = 1 ).
Here, as in Ref. [1], we consider without loss of generality two cases: (i) m < k 1 < k 2 and (ii) 1 < k 1 = k 2 = k , k m . (In the case m = k 1 = k 2 , the solutions are absent due to our restrictions.) For H 0 in both cases the solutions exist only if α ε = ε α 2 / α 1 > 0 , Λ ε > 0 and multidimensional cosmological term Λ obeys the bounds: 0 < λ ( m , k 1 , k 2 ) Λ α λ + ( m , k 1 , k 2 ) . For H = 0 , the solutions exist only when α ε > 0 , Λ ε > 0 , k 1 k 2 and Λ α = λ ( k 1 , k 2 ) > 0 . We note that here, as in Ref. [1], we use the Chirkov–Pavluchenko–Toporensky scheme of reduction of the set of polynomial equations [17]. As in Ref. [1] we reduce the problem in the generic H 0 case to solutions of a single polynomial master equation of the fourth order or less, which may be solved in radicals for all m > 1 , k 1 > 1 and k 2 > 1 . In the case (ii) 1 < k 1 = k 2 = k , k m ( H 0 ), the solutions for Hubble-like parameters are found explicitly (see Section 4).
We also study (in Section 5) the stability of the solutions for u ± in a class of cosmological type solutions with diagonal metrics by using an extension of the results of Refs. [1,21] (see also the approach of Ref. [18]) and single out the subclasses of stable/non-stable solutions.
We note that the exponential cosmological type solutions with two non-coinciding Hubble-like parameters H 0 and h obeying S 1 = m H + l h 1 0 with m > 2 , l > 2 were studied earlier in Ref. [29]. In that case there were two sets of solutions obeying: (a) ε α > 0 , α Λ < λ + ( m , l ) and (b) ε α < 0 , α Λ < λ ( m , l ) , where λ ± ( m , l ) > 0 and ε = ± 1 .
It should be noted that, recently, EGB models were used for constricting certain 4-dimensional gravitational models (so-called 4DEGB theories, e.g., belonging to Horndeski class) by using ideas of Glavan–Lin rescaling [30] and/or dimensional reductions. These 4D modified models of gravity are (at the moment) under intensive study and have numerous applications in gravitational physics and cosmology, for a review see Ref. [31].

2. The Cosmological Model

We start with the model governed by the action:
S = M d D z | g | { α 1 ( R [ g ] 2 Λ ) + α 2 L 2 [ g ] } .
Here, g = g M N d z M d z N is the metric on a manifold M ( dim M = D ), | g | = | det ( g M N ) | , Λ is the cosmological term, R [ g ] is scalar curvature,
L 2 [ g ] = R M N P Q R M N P Q 4 R M N R M N + R 2
is the Gauss–Bonnet term and α 1 , α 2 are certain nonzero constants.
Our choice of the manifold is following:
M = R × M 1 × × M n .
In what follows we deal with the metric:
g = ε d u d u + i = 1 n B i ε i e 2 v i u d y i d y i .
Here, B i > 0 are arbitrary constants, ε = ± 1 , ε i = ± 1 , i = 1 , , n ( n > 3 ) and M 1 , , M n are chosen to be 1-dimensional manifolds (either non-compact ( R ) or compact ( S 1 ) ones). The cosmological case ( ε = ε i = 1 ) was considered in detail in Ref. [1]. The case ε = 1 may describe certain static configurations.
The action (1) with the ansatz for the metric (1) imposed gives rise to the equations of motion which are of polynomial type [20]:
E = G i j v i v j + 2 Λ ε α ε G i j k l v i v j v k v l = 0 ,
Y i = 2 G i j v j 4 3 α ε G i j k l v j v k v l i = 1 n v i 2 3 G i j v i v j + 8 3 Λ ε = 0 ,
i = 1 , , n . Here we denote α = α 2 / α 1 and
G i j = δ i j 1 , G i j k l = G i j G i k G i l G j k G j l G k l .
Refs. [13,14]. For n > 3 , we have a set of polynomial equations of order 4.
For the case n > 3 , Λ = 0 and α ε < 0 the set of Equations (4) and (5) has a trivial (isotropic) solution: v 1 = = v n = H [13,14], which was generalized in Ref. [16] to the case Λ 0 .
In Refs. [13,14], the following proposition was proved: there are no more than three different numbers among v 1 , , v n if Λ = 0 . This proposition was generalised in Ref. [21] for Λ 0 , when the following condition is imposed i = 1 n v i 0 .
In this paper we study solutions to Equations (4) and (5) by using the following ansatz:
v = ( H , , H m , h 1 , , h 1 k 1 , h 2 , , h 2 k 2 ) .
Here, H is the Hubble-like parameter which corresponds to an m-dimensional factor space with inequality m > 1 imposed, while h 1 is the Hubble-like parameter which is related to an k 1 -dimensional factor space with k 1 > 1 and h 2 is the Hubble-like parameter assigned to an k 2 -dimensional factor space with k 2 > 1 .
In what follows we add additional restrictions to our ansatz (7):
H h 1 , H h 2 , h 1 h 2 , S 1 = m H + k 1 h 1 + k 2 h 2 0 .
It was shown in Ref. [22] that the set of ( n + 1 ) polynomial Equations (4) and (5) under ansatz (7) and restrictions (8) obeyed are equivalent to a set of polynomial equations:
E = 0 ,
Q = 1 2 α ε ,
L = H + h 1 + h 2 S 1 = 0 ,
which are of fourth, second and first orders, respectively. Here, E is defined in (4) and
Q = Q h 1 h 2 = S 1 2 S 2 2 S 1 ( h 1 + h 2 ) + 2 ( h 1 2 + h 1 h 2 + h 2 2 ) ,
where
S k = i = 1 n ( v i ) k .
For more general prescription of a scheme of the reduction of polynomial equations of motion see Ref. [17] (the so-called Chirkov–Pavluchenko–Toporensky trick).
Relation (10) is a special case of more general relations [22]:
Q h i h j = S 1 2 S 2 2 S 1 ( h i + h j ) + 2 ( h i 2 + h i h j + h j 2 ) = 1 2 α ε , i j ,
i , j = 0 , 1 , 2 , with notation h 0 = H used.
Relation (8) excludes the following case H = h 1 = h 2 = 0 . In the main body of the paper we put:
H 0 .
As in Ref. [1] we denote:
x 1 = h 1 / H , x 2 = h 2 / H .
In terms of dimensionless parameters the restrictions (8) may rewritten as follows:
x 1 1 , x 2 1 , x 1 x 2 , m + k 1 x 1 + k 2 x 2 0 .
Equation (11) is equivalent to the following one:
m 1 + ( k 1 1 ) x 1 + ( k 2 1 ) x 2 = 0 .
In what follows we do not consider the case,
m = k 1 = k 2 ,
which lead us to the empty set of solutions, since we find for m = k 1 = k 2 > 1 from restriction (17): 1 + x 1 + x 2 0 , while (18) implies 1 + x 1 + x 2 = 0 .
Due to (10) and (12) we obtain:
2 α ε P H 2 = 1 ,
where
P = P ( x 1 , x 2 ) ( m + k 1 x 1 + k 2 x 2 ) 2 ( m + k 1 x 1 2 + k 2 x 2 2 ) 2 ( m + k 1 x 1 + k 2 x 2 ) ( x 1 + x 2 ) + 2 ( x 1 2 + x 1 x 2 + x 2 2 ) .
The relation (20) is valid for α ε P < 0 . It can be readily proved that [1]:
P < 0
for m > 1 , k 1 > 1 , k 2 > 1 . Indeed [1],
P = 1 m + ( 1 k 1 ) x 1 2 + ( 1 k 2 ) x 2 2 < 0 .
It follows from (22) that:
α ε > 0 .
Equation (9) reads [1]:
2 Λ ε = G i j v i v j + α ε G i j k l v i v j v k v l = H 2 V 1 + α ε H 4 V 2 ,
where
V 1 = V 1 ( x 1 , x 2 ) = m k 1 x 1 2 k 2 x 2 2 + ( m + k 1 x 1 + k 2 x 2 ) 2
and
V 2 = V 2 ( x 1 , x 2 ) = [ m ] 4 + 4 [ m ] 3 ( k 1 x 1 + k 2 x 2 ) + 6 [ m ] 2 [ k 1 ] 2 x 1 2 + 2 k 1 k 2 x 1 x 2 + [ k 2 ] 2 x 2 2 + 4 m [ k 1 ] 3 x 1 3 + 3 [ k 1 ] 2 k 2 x 1 2 x 2 + 3 k 1 [ k 2 ] 2 x 1 x 2 2 + [ k 2 ] 3 x 2 3 + [ k 1 ] 4 x 1 4 + 4 [ k 1 ] 3 k 2 x 1 3 x 2 + 6 [ k 1 ] 2 [ k 2 ] 2 x 1 2 x 2 2 + 4 k 1 [ k 2 ] 3 x 1 x 2 3 + [ k 2 ] 4 x 2 4 .
Here, [ N ] k = N ( N 1 ) . . . ( N k + 1 ) .
Due to (20) we obtain:
λ = α Λ = V 1 4 P + V 2 8 P 2 ,
or
V 2 ( x 1 , x 2 ) 2 P ( x 1 , x 2 ) V 1 ( x 1 , x 2 ) 8 ( P ( x 1 , x 2 ) ) 2 λ = 0 .
Owing to Equation (18) we get:
x 2 = x 2 ( x 1 ) = m 1 k 2 1 k 1 1 k 2 1 x 1 .
Hence, from Equation (29) we get a master equation in x 1 variable:
V 2 ( x 1 , x 2 ( x 1 ) ) 2 P ( x 1 , x 2 ( x 1 ) ) V 1 ( x 1 , x 2 ( x 1 ) ) 8 ( P ( x 1 , x 2 ( x 1 ) ) ) 2 λ = 0 .
This polynomial equation is of fourth order or less (this depends upon the value of λ ). One may solve it in radicals for all m > 1 , k 1 > 1 and k 2 > 1 .
Relations (23) and (30) imply the identity:
( k 2 1 ) P ( x 1 , x 2 ( x 1 ) ) = ( k 1 1 ) ( k 1 + k 2 2 ) x 1 2 + 2 ( m 1 ) ( k 1 1 ) x 1 + ( m 1 ) ( m + k 2 2 ) ,
which will be used below.

3. The Case k 1 k 2

In this section we put k 1 k 2 . We rewrite relation (28) as follows:
λ = f ( x 1 ) V 1 ( x 1 , x 2 ( x 1 ) ) 4 P ( x 1 , x 2 ( x 1 ) ) + V 2 ( x 1 , x 2 ( x 1 ) ) 8 ( P ( x 1 , x 2 ( x 1 ) ) ) 2 .
Due to (30), we present restrictions (17) in the following form [1]:
x 1 X 1 , x 1 X 2 , x 1 X 3 , x 1 X 4 ,
where
X 1 = 1 ,
X 2 = m + k 2 2 k 1 1 ,
X 3 = m 1 k 1 + k 2 2 ,
X 4 = m k 2 k 2 k 1 .

Extremum Points

We obtain [1]:
d f d x 1 = C ( m , k 1 , k 2 ) ( x 1 X 1 ) ( x 1 X 2 ) ( x 1 X 3 ) ( x 1 X 4 ) ( k 2 1 ) P ( x 1 , x 2 ( x 1 ) ) 3 ,
where
C ( m , k 1 , k 2 ) = ( m 1 ) ( k 1 1 ) 2 ( k 2 k 1 ) ( k 1 + k 2 2 )
and X 1 , X 2 , X 3 , X 4 are given by (35)–(38). Thus, the extreme points of the function f ( x 1 ) are excluded from our consideration Due to (8) we are ought to exclude the extreme points of f ( x 1 ) .
For λ i = f ( X i ) , i = 1 , 2 , 3 , 4 , we have [1]:
λ 1 = λ 1 ( m , k 1 , k 2 ) = u ( k 2 , m + k 1 ) 8 ( m + k 1 + k 2 3 ) ( m + k 1 2 ) ( k 2 1 ) ,
λ 2 = λ 2 ( m , k 1 , k 2 ) = u ( k 1 , m + k 2 ) 8 ( m + k 1 + k 2 3 ) ( m + k 2 2 ) ( k 1 1 ) ,
λ 3 = λ 3 ( m , k 1 , k 2 ) = u ( m , k 1 + k 2 ) 8 ( m 1 ) ( k 1 + k 2 2 ) ( m + k 1 + k 2 3 ) ,
λ 4 = λ 4 ( m , k 1 , k 2 ) = v ( m , k 1 , k 2 ) 8 w ( m , k 1 , k 2 ) .
Here,
u ( m , l ) = l m 2 + ( l 2 8 l + 8 ) m + l ( l 1 ) ,
v ( m , l , k ) = ( k + l ) m 2 + ( m + l ) k 2 + ( m + k ) l 2 6 m l k ,
w ( m , l , k ) = ( k + l 2 ) m 2 + ( m + l 2 ) k 2 + ( m + k 2 ) l 2 + 2 m l + 2 m k + 2 l k 6 m l k .
It was verified in Ref. [1] that:
λ i = λ i ( m , k 1 , k 2 ) > 0
for m > 1 , k 1 > 1 , k 2 > 1 , i = 1 , 2 , 3 , 4 .
In the limit x 1 ± we obtain:
λ = lim x 1 f ( x 1 ) = ( k 1 + k 2 6 ) k 1 k 2 + k 1 2 + k 2 2 + k 1 + k 2 8 ( k 1 1 ) ( k 2 1 ) ( k 1 + k 2 2 ) .
Here we obtain [1]:
λ = λ ( k 1 , k 2 ) = λ ( k 2 , k 1 ) > 0 ,
for all k 1 > 1 and k 2 > 1 .
The definitions of X i imply [1]:
X 2 < X 3 < 0 < X 1 = 1 .
Here, m > 1 , k 1 > 1 and k 2 > 1 .
From this point up to Section 4 we impose the following inequality:
1 < m < k 1 < k 2 .
It was shown in Ref. [1] that:
0 < λ 1 < λ 2 < λ 3 ,
0 < λ 1 < λ 4 < λ 3 ,
and
( A + ) X 4 < X 2 , λ 4 > λ 2 , for 2 k 1 m k 2 > 0 ,
( A ) X 4 > X 2 , λ 4 < λ 2 , for 2 k 1 m k 2 < 0 ,
and
( A 0 ) X 4 = X 2 , λ 4 = λ 2 , for 2 k 1 m k 2 = 0 .
For ( m , k 1 , k 2 ) = ( 4 , 6 , 7 ) the function λ = f ( x 1 ) is presented graphically in Figure 1.
It was proved in Ref. [1] that:
λ 1 < λ < λ 3 .
By using (40) and (52) we get:
C ( m , k 1 , k 2 ) > 0 .
It was proved in Ref. [1] that for the function f ( x 1 ) mentioned above X 3 is the point of absolute maximum and X 1 is the point of absolute minimum, i.e.,
λ 1 λ = f ( x 1 ) λ 3 ,
for all x 1 R . We remind that according to (34) the points X 1 , X 2 , X 3 , X 4 are forbidden for our analysis. We obtain:
λ 1 < λ = f ( x 1 ) < λ 3 ,
for all x 1 X 1 , X 2 , X 3 , X 4 . Let us denote the set of definition of the function f for our consideration ( , ) * { x | x R , x X 1 , X 2 , X 3 , X 4 } . Since the function f ( x 1 ) is a continuous one, the image of the function f (due to intermediate value theorem) is: f ( ( , ) * ) = ( λ 1 , λ 3 ) . Thus, we are led to the following proposition.
Proposition 1.
The solutions to Equations (4) and (5) for ansatz (7) with 1 < m < k 1 < k 2 obeying the inequalities H 0 , H h 1 , H h 2 , h 1 h 2 and S 1 = m H + k 1 h 1 + k 2 h 2 0 do exist if and only if α ε > 0 and
0 < λ 1 < α Λ < λ 3 ,
where λ 1 and λ 3 are defined in (41) and (43), respectively. In this case x 1 = h 1 / H X 1 , X 2 , X 3 , X 4 (see (35)–(38)), x 2 = h 2 / H = x 2 ( x 1 ) is given by (30), x 1 obeys the polynomial master Equation (31) (of fourth order or less) and H 2 is given by (20) and (21).
The case H = 0 . In the case H = 0 the solutions under consideration take place only if α ε > 0 , Λ ε > 0 and
α Λ = λ ( k 1 , k 2 ) = ( k 1 + k 2 6 ) k 1 k 2 + k 1 2 + k 2 2 + k 1 + k 2 8 ( k 1 1 ) ( k 2 1 ) ( k 1 + k 2 2 ) > 0 ,
where k 1 k 2 . Indeed, relation (11) reads as ( k 1 1 ) h 1 + ( k 2 1 ) h 2 = 0 , and relation (10) is equivalent to ( k 1 1 ) ( h 1 ) 2 + ( k 2 1 ) ( h 2 ) 2 = 1 / ( 2 α ε ) . From these relations we get α ε > 0 and
h 1 = ± k 2 1 2 α ε ( k 1 1 ) ( k 1 + k 2 2 ) 1 / 2 ,
h 2 = k 1 1 2 α ε ( k 2 1 ) ( k 1 + k 2 2 ) 1 / 2 ,
which imply, due to H = 0 and (9), the relation (63).

4. The Case k 1 = k 2

We will now turn our attention to the case H 0 , m > 1 and k 1 = k 2 = k > 1 . Due to (18) we obtain:
m 1 + ( k 1 ) ( x 1 + x 2 ) = 0 .
It follows from (23) that:
P = 1 m + ( 1 k ) ( x 1 2 + x 2 2 ) .
Since the case of equal factor-space dimensions is excluded from our consideration (see Section 2) we put:
m k
and α ε > 0 .
Denoting
X α ε H 2 ,
α ε > 0 , we obtain from (20) that
X P = 1 2 .
Relation (69) implies
H = ϵ 0 X / α ε , ϵ 0 = ± 1 .
Plugging the relations (66), (67) into (26), (27) we obtain
V 1 = [ ( m 1 ) ( m k ) + P k ( k 1 ) ] / ( k 1 ) 2 ,
V 2 = [ ( m 1 ) ( m k ) ( m + k 2 ) ( m + 2 k 3 ) + 3 P 2 ( k 1 ) 2 k ] / ( k 1 ) 3 .
By virtue relation (70), we present relation (28) as:
2 λ = 2 α Λ = X V 1 + X 2 V 2 ,
or in an equivalent manner as:
A X 2 + B X + C = 0 .
Here,
A = ( m 1 ) ( m k ) ( m + k 2 ) ( m + 2 k 3 ) ,
B = ( m 1 ) ( m k ) ( k 1 ) ,
C = 1 4 k ( k 1 ) 2 + 2 λ ( k 1 ) 3 .
It follows from (68) that A 0 . The calculation of the discriminant D = B 2 4 A C leads us to the following identity:
D = ( m 1 ) ( m k ) ( k 1 ) 2 ( F 8 λ f ) ,
where we denote
F = F ( m , k ) = ( m 1 ) ( m k ) + ( m + k 2 ) ( m + 2 k 3 ) k ,
f = f ( m , k ) = ( m + k 2 ) ( m + 2 k 3 ) ( k 1 ) > 0 .
It was verified in Ref. [1] that F = F ( m , k ) > 0 for all m > 1 , k > 1 and k m .
By solving Equation (75) we get [1]:
X = ( B + ϵ ¯ 1 D ) / ( 2 A ) , ϵ ¯ 1 = ± 1 .
We seek real solutions obeying:
D > 0 ,
X > 0 .
The case D = 0 should be excluded [1]. Indeed, D = 0 implies either x 1 = 1 or x 2 = 1 , which is in contradiction with (17).
Here, we rewrite the inequality (83) as:
λ < λ 1 for m > k ,
λ > λ 1 for m < k ,
where
λ 1 = λ 1 ( m , k , k ) = F ( m , k ) / ( 8 f ( m , k ) ) .
Equations (66) and (67) may be resolved as:
x 1 = ( ϵ 2 E + m 1 ) / ( 2 k 2 ) ,
x 2 = ( ϵ 2 E + m 1 ) / ( 2 k 2 ) ,
where ε 2 = ± 1 and
E = ( m 1 ) ( m + 2 k 3 ) 2 P ( k 1 ) = ( k 1 ) X 1 ( m 1 ) ( m + 2 k 3 ) .
Here, one should consider the case:
E > 0 .
Indeed, E = 0 implies x 1 = x 2 , which is not allowed by (17). Due to (84) and (91) we obtain:
0 < X < k 1 ( m 1 ) ( m + 2 k 3 ) .
It was verified in Ref. [1] that relations (88), (89) and (92) imply all four inequalities in (17).
Now we proceed with inequalities in (92). By introducing the parameter:
ϵ 1 = ϵ ¯ 1 sign ( m k ) ,
we rewrite relation (82) in the following form:
X = k 1 2 ( m + k 2 ) ( m + 2 k 3 ) + ϵ 1 D 2 | A | ,
ϵ 1 = ± 1 .
First, we consider the case ϵ 1 = 1 . The second inequality in (92) X < k 1 ( m 1 ) ( m + 2 k 3 ) is valid due to 2 ( m + k 2 ) > m 1 . As to the first inequality X > 0 , we obtain:
0 < D < ( m 1 ) | m k | ( k 1 ) .
Due to definition of D in (79) we get:
0 < ( m 1 ) ( m k ) ( k 1 ) 2 ( F 8 λ f ) < ( m 1 ) 2 | m k | 2 ( k 1 ) 2 .
Relations (96) may be presented in the following form:
F < 8 λ f < F , for m > k ,
F < 8 λ f < F , for m < k .
Here,
F F ( m 1 ) ( m k ) .
By using relations:
F 8 f = k 8 ( k 1 ) = λ = λ ( k , k ) ,
where λ ( k , l ) is defined in (49), and (87) and (100) one can present relations (97), (98) in the following form:
λ < λ < λ 1 , for m > k ,
λ 1 < λ < λ , for m < k .
Now, we consider the case ϵ 1 = 1 . Since the inequality X > 0 is obeyed in this case, one should verify the inequality X < k 1 ( m 1 ) ( m + 2 k 3 ) . We find:
0 < D < | m k | ( m + 2 k 3 ) ( k 1 ) ,
or
0 < ( m 1 ) ( m k ) ( F 8 λ f ) < | m k | 2 ( m + 2 k 3 ) 2 .
We write relations (104) in the following form:
F + < 8 λ f < F , for m > k ,
F < 8 λ f < F + , for m < k ,
where
F + F ( m 1 ) 1 ( m k ) ( m + 2 k 3 ) 2 .
Here one can verify that:
F + 8 f = λ 3 = λ 3 ( m , k , k ) .
Due to (87) and (108) we rewrite relations (105), (106) in the following form:
λ 3 < λ < λ 1 , for m > k ,
λ 1 < λ < λ 3 , for m < k .
Here,
λ 1 < λ < λ 3
for m < k , while
λ 3 < λ < λ 1
for k < m . The inequalities in (112) just follow from inequalities F + < F < F for k < m .
Thus, we are led to the following generalisation of the Proposition 2 from Ref. [1].
Proposition 2.
The solutions to Equations (4) and (5) for ansatz (7) imposed with 1 < m , 1 < k 1 = k 2 = k , m k , obeying the inequalities H 0 , H h 1 , H h 2 , h 1 h 2 , S 1 = m H + k h 1 + k h 2 0 do exist if and only if α ε > 0 ,
λ 1 < λ = α Λ < λ 3
for m < k and
λ 3 < λ = α Λ < λ 1 ,
where λ 1 = λ 1 ( k , k ) , λ 3 = λ 3 ( k , k ) are defined in (41) and (43). In this case H satisfies the relation (71) with X from (94), x 1 = h 1 / H and x 2 = h 2 / H are given by relations (88) and (89), λ obeys (101), (102) for ϵ 1 = 1 and (109), (110) for ϵ 1 = 1 with λ = k 8 ( k 1 ) .
The case H = 0 . For k 1 = k 2 = k > 1 and H = 0 the solutions under consideration obeying restrictions (8) are absent [1].

5. The Analysis of Stability

Here, we analyse the stability of our solutions along a line as was done in Refs. [20,21,22].
We impose the following restriction:
det ( L i j ( v ) ) 0 ,
where
L = ( L i j ( v ) ) = ( 2 G i j 4 α ε G i j k s v k v s ) .
Here, one should deal with general cosmological type setup with the metric:
g = ε d u d u + i = 1 n e 2 β i ( u ) ε i d y i d y i ,
where ε = ± 1 , ε i = ± 1 , i = 1 , , n . For the equations of motion we obtain [28]:
E = G i j h i h j + 2 Λ ε α ε G i j k l h i h j h k h l = 0 ,
Y i = d L i d t + ( j = 1 n h j ) L i 2 3 ( G s j h s h j 4 Λ ε ) = 0 ,
where h i = β ˙ i = d β i d u ,
L i = L i ( h ) = 2 G i j h j 4 3 α ε G i j k l h j h k h l ,
i = 1 , , n .
According to previous consideration of Ref. [21] the solution ( h i ( t ) ) = ( v i ) ( i = 1 , , n ; n > 3 ) to Equations (118) and (119) which obeys the restrictions (115) is stable under perturbations:
h i ( t ) = v i + δ h i ( t ) ,
i = 1 , , n , as u + , if and only if
S 1 ( v ) = i = 1 n v i > 0
and it is unstable, as u + , if and only if
S 1 ( v ) = i = 1 n v i < 0 .
In the limit u , the stability condition is given by (123) while the instability condition reads as (122). These conditions just follow from solutions for perturbations δ h i ( t ) = C i exp ( S 1 ( v ) u ) ( C i = const 0 ) which are valid in the leading order.
Here, a key point is the verification of the relation (115). It was fulfilled in Ref. [22] by using first three relations in (8) and (14) and k 1 > 1 , k 2 > 1 and m > 1 .
First we consider the case 1 < m < k 1 < k 2 . By using (18) we find that for H > 0 the condition (122) may be written as:
m + k 1 x 1 + k 2 x 2 = 1 + x 1 + x 2 > 0
or, equivalently,
x 1 > X 4 = m k 2 k 2 k 1 .
For H < 0 the stability condition (122) is as follows:
x 1 < X 4 .
The non-stability condition (123) for u + reads as (126) for H > 0 and as (125) for H < 0 . These conditions are reversed in case u .
Proposition 3.
Let us consider cosmological type solutions to Equations (4) and (5) for ansatz (7) with 1 < k 1 < k 2 , obeying the inequalities H 0 , H h 1 , H h 2 , h 1 h 2 , S 1 = m H + k 1 h 1 + k 2 h 2 0 .
(a) 
Let H > 0 . For u + the solutions are stable if x > X 4 and unstable if x < X 4 , while for u they are stable if x < X 4 and unstable if x > X 4 ;
(b) 
Let H < 0 . For u + the solutions are stable if x < X 4 and unstable if x > X 4 , while for u they are stable if x > X 4 and unstable if x < X 4 .
Now we proceed with considering the case H 0 , 1 < m , 1 < k 1 = k 2 = k , m k . Since x 1 1 , x 2 1 and x 1 x 2 the exact solutions under consideration obey the first three relations in (8), which imply the validity of the key restriction (115).
For the stability condition (122) as u + in this case we get:
H ( m + k 1 x 1 + k 2 x 2 ) = H ( 1 + x 1 + x 2 ) = H 1 m 1 k 1 > 0 ,
or, equivalently,
H ( k m ) > 0 .
The non-stability condition (123) for u + may be written as:
H ( k m ) < 0 .
Thus, we have the following proposition.
Proposition 4.
Let us consider cosmological type solutions (4), (5) for ansatz (7) with 1 < m , 1 < k 1 = k 2 = k , m k , obeying the inequalities H 0 , H h 1 , H h 2 , h 1 h 2 , S 1 = m H + k h 1 + k h 2 0 , is stable, as u + , if and only if H ( k m ) > 0 and it is unstable, as u + , if and only if H ( k m ) < 0 .
(c) 
Let H > 0 . For u + the solutions are stable if k > m and unstable if k < m , while for u they are stable if k < m and unstable if k > m .
(d) 
Let H < 0 . For u + the solutions are stable if k < m and unstable if k > m , while for u they are stable if k > m and unstable if k < m
The case H = 0 . For a completeness we consider the solutions with H = 0 and h 1 , h 2 from (64) and (65), where k 1 k 2 , k 1 > 1 , k 2 > 1 , α ε > 0 and Λ is given by (63). We get:
S 1 = k 1 h 1 + k 2 h 2 = ± ( k 2 k 1 ) 2 α ε ( k 1 1 ) ( k 2 1 ) ( k 1 + k 2 2 ) 1 / 2 .
Here, ± is a sign parameter in (64) and (65). By using our analysis presented above we obtain that the solution with ± ( k 2 k 1 ) > 0 is stable, as u + . This occurs if either k 2 > k 1 and the sign + " are selected in (64) and (65), or if k 2 < k 1 and the sign " are chosen. For the case ± ( k 2 k 1 ) < 0 our solution is unstable, as u + . (Here we also assume the restriction m > 1 ). These conditions are reversed in case u .

6. Conclusions

We have studied the D-dimensional Einstein–Gauss–Bonnet (EGB) model with the Λ -term and two non-zero constants α 1 and α 2 . By dealing with diagonal cosmological type metrics, we have considered a class of solutions with exponential dependence of three scale factors (upon u-variable) for any α = α 2 / α 1 0 , signature parameter ε = ± 1 and generic dimensionless parameter Λ α .
More precisely speaking, we have described a class of cosmological type solutions with exponential dependence of three scale factors, governed by three non-coinciding Hubble-like parameters H, h 1 and h 2 . These parameters correspond, respectively, to factor spaces of dimensions m > 1 , k 1 > 1 and k 2 > 1 ( D = 1 + m + k 1 + k 2 ), and obey the following restriction S 1 = m H + k 1 h 1 + k 2 h 2 0 . We have analyzed two cases: (i) m < k 1 < k 2 and (ii) 1 < k 1 = k 2 = k m . This choice does not restrict the generality, since, as it was shown, there are no solutions under consideration for k 1 = k 2 = m ). It was shown that the solutions exist only if λ = α Λ > 0 and the (dimensionless) parameter λ obey certain restrictions, e.g., upper and lower bounds for H 0 , which depend upon dimensions m, k 1 and k 2 (Proposition 1). In case (ii) we have presented explicit solutions for all k > 1 and k m ( Proposition 2).
By using the Chirkov–Pavluchenko–Toporensky splitting trick from Ref. [17], we have reduced the problem for H 0 to a master equation on the dimensionless variable x 1 = h 1 / H . This equation is of the fourth order (in the generic case) or less (depending on λ ), and may be solved in radicals for all m > 1 , k 1 > 1 , k 2 > 1 and λ . The master equation does not depend upon the signature parameter ε = ± 1 which only controls the sign of α according to inequality α ε > 0 . Due to bounds obtained λ = α Λ > 0 . (This is valid also for H = 0 ). Hence the solutions under consideration do exist if Λ ε > 0 , i.e., when Λ > 0 in the cosmological case ( ε = 1 ) and Λ < 0 in the static case ( ε = 1 ). Here there are no solutions under consideration for Λ = 0 —contrary to the case of two factor spaces [29,32].
Here we have analyzed the stability of solutions as u ± in a class of cosmological type solutions with diagonal metrics. In both cases ((i) and (ii)) for H 0 , the “islands” of stability and instability were singled out. (The case H = 0 was also analysed.) We have shown that in case (i) the solutions with H > 0 are stable as u for x 1 = h 1 / H > X 4 = m k 2 k 2 k 1 and unstable as u for x 1 < X 4 (see Proposition 3). These conditions should be reversed when we consider the case H > 0 , u or we deal with H < 0 , u + (see Proposition 3). It was proved that in case (ii), the solutions with H > 0 are stable as u for k > m and unstable as u for k < m (see Proposition 4). For a given choice of asymptotic u ± , the stability condition for H < 0 is equivalent to the instability conditions for H > 0 and vice versa.
We have also found that the solution with H = 0 exists only for k 1 k 2 , α ε > 0 and a fixed value of ε Λ > 0 depending upon k 1 and k 2 . Here we have two opposite in sign solutions for ( h 1 , h 2 ) with one solution being stable ( u ± ) and the second one unstable, depending upon the sign of k 1 k 2 .
Some cosmological applications of the model ( ε = 1 ), e.g., in the context of the variation of the gravitational constant, were considered in Refs. [1,33,34]. For the static case ( ε = 1 ), possible applications of the obtained solutions may be a subject of a further research, aimed towards a search for topological black hole solutions (with flat horizon) or wormhole solutions, which are coinciding asymptotically (for ( u ± )) with our solutions.

Author Contributions

V.D.I. is responsible for mathematical model development and K.K.E.—for simulation model development. All authors have read and agreed to the published version of the manuscript.

Funding

RUDN University Strategic Academic Leadership Program, RFBR grant 19-02-00346.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

This paper has been supported by the RUDN University Strategic Academic Leadership Program (recipients: V.D.I.—mathematical model development and K.K.E.—simulation model development).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Ernazarov, K.K.; Ivashchuk, V.D. Stable exponential cosmological solutions with three different Hubble-like parameters in EGB model with a Λ-term. Eur. Phys. J. C 2020, 80, 543. [Google Scholar] [CrossRef]
  2. Zwiebach, B. Curvature squared terms and string theories. Phys. Lett. B 1985, 156, 315. [Google Scholar] [CrossRef]
  3. Fradkin, E.S.; Tseytlin, A.A. Effective action approach to superstring theory. Phys. Lett. B 1985, 160, 69–76. [Google Scholar] [CrossRef] [Green Version]
  4. Ishihara, H. Cosmological solutions of the extended Einstein gravity with the Gauss–Bonnet term. Phys. Lett. B 1986, 179, 217. [Google Scholar] [CrossRef]
  5. Deruelle, N. On the approach to the cosmological singularity in quadratic theories of gravity: The Kasner regimes. Nucl. Phys. B 1989, 327, 253–266. [Google Scholar] [CrossRef] [Green Version]
  6. Nojiri, S.; Odintsov, S.D. Introduction to modified gravity and gravitational alternative for Dark Energy. Int. J. Geom. Meth. Mod. Phys. 2007, 4, 115–146. [Google Scholar] [CrossRef] [Green Version]
  7. Elizalde, E.; Makarenko, A.N.; Obukhov, V.V.; Osetrin, K.E.; Filippov, A.E. Stationary vs. singular points in an accelerating FRW cosmology derived from six-dimensional Einstein–Gauss–Bonnet gravity. Phys. Lett. B 2007, 644, 1–6. [Google Scholar] [CrossRef] [Green Version]
  8. Bamba, K.; Guo, Z.-K.; Ohta, N. Accelerating Cosmologies in the Einstein–Gauss–Bonnet theory with dilaton. Prog. Theor. Phys. 2007, 118, 879–892. [Google Scholar] [CrossRef] [Green Version]
  9. Toporensky, A.; Tretyakov, P. Power-law anisotropic cosmological solution in 5 + 1 dimensional Gauss–Bonnet gravity. Grav. Cosmol. 2007, 13, 207–210. [Google Scholar]
  10. Pavluchenko, S.A.; Toporensky, A.V. A note on differences between (4 + 1)- and (5 + 1)-dimensional anisotropic cosmology in the presence of the Gauss–Bonnet term. Mod. Phys. Lett. A 2009, 24, 513–521. [Google Scholar] [CrossRef] [Green Version]
  11. Pavluchenko, S.A. On the general features of Bianchi-I cosmological models in Lovelock gravity. Phys. Rev. D 2009, 80, 107501. [Google Scholar] [CrossRef] [Green Version]
  12. Kirnos, I.V.; Makarenko, A.N.; Pavluchenko, S.A.; Toporensky, A.V. The nature of singularity in multidimensional anisotropic Gauss–Bonnet cosmology with a perfect fluid. Gen. Rel. Grav. 2010, 42, 2633–2641. [Google Scholar] [CrossRef] [Green Version]
  13. Ivashchuk, V.D. On anisotropic Gauss–Bonnet cosmologies in (n + 1) dimensions, governed by an n-dimensional Finslerian 4-metric. Grav. Cosmol. 2010, 16, 118–125. [Google Scholar] [CrossRef] [Green Version]
  14. Ivashchuk, V.D. On cosmological-type solutions in multidimensional model with Gauss–Bonnet term. Int. J. Geom. Methords Mod. Phys. 2010, 7, 797–819. [Google Scholar] [CrossRef] [Green Version]
  15. Maeda, K.-I.; Ohta, N. Cosmic acceleration with a negative cosmological constant in higher dimensions. J. High Energy Phys. 2014, 1406, 95. [Google Scholar] [CrossRef] [Green Version]
  16. Chirkov, D.; Pavluchenko, S.; Toporensky, A. Exact exponential solutions in Einstein–Gauss–Bonnet flat anisotropic cosmology. Mod. Phys. Lett. A 2014, 29, 1450093. [Google Scholar] [CrossRef] [Green Version]
  17. Chirkov, D.; Pavluchenko, S.A.; Toporensky, A. Non-constant volume exponential solutions in higher-dimensional Lovelock cosmologies. Gen. Rel. Grav. 2015, 47, 137. [Google Scholar] [CrossRef] [Green Version]
  18. Pavluchenko, S.A. Stability analysis of exponential solutions in Lovelock cosmologies. Phys. Rev. D 2015, 92, 104017. [Google Scholar] [CrossRef] [Green Version]
  19. Pavluchenko, S.A. Cosmological dynamics of spatially flat Einstein–Gauss–Bonnet models in various dimensions: Low-dimensional Λ-term case. Phys. Rev. D 2016, 94, 084019. [Google Scholar] [CrossRef] [Green Version]
  20. Ernazarov, K.K.; Ivashchuk, V.D.; Kobtsev, A.A. On exponential solutions in the Einstein–Gauss–Bonnet cosmology stability and variation of G. Grav. Cosmol. 2016, 22, 245–250. [Google Scholar] [CrossRef]
  21. Ivashchuk, V.D. On stability of exponential cosmological solutions with non-static volume factor in the Einstein–Gauss–Bonnet model. Eur. Phys. J. C 2016, 76, 431. [Google Scholar] [CrossRef] [Green Version]
  22. Ernazarov, K.K.; Ivashchuk, V.D. Stable exponential cosmological solutions with zero variation of G and three different Hubble-like parameters in the Einstein–Gauss–Bonnet model with a Λ-term. Eur. Phys. J. C 2017, 77, 402. [Google Scholar] [CrossRef] [Green Version]
  23. Ivashchuk, V.D.; Kobtsev, A.A. Stable exponential cosmological solutions with two factor spaces in the Einstein–Gauss–Bonnet model with a Λ-term. Gen. Rel. Grav. 2018, 50, 119. [Google Scholar] [CrossRef] [Green Version]
  24. Benetti, M.; da Costa, S.S.; Capozziello, S.; Alcaniz, J.S.; Laurentis, M.D. Observational constraints on Gauss–Bonnet cosmology. Int. J. Mod. Phys. 2018, 27, 1850084. [Google Scholar] [CrossRef] [Green Version]
  25. Nojiri, S.; Odintsov, S.D.; Oikonomou, V.K. Unifying Inflation with Early and Late-time Dark Energy in F(R) Gravity. Phys. Dark Univ. 2020, 29, 100602. [Google Scholar] [CrossRef]
  26. Riess, A.G.; Filippenko, A.V.; Challis, P.; Clocchiatti, A.; Diercks, A.; Garnavich, P.M.; Gilliland, R.L.; Hogan, C.J.; Jha, S.; Kirshner, R.P. Observational evidence from supernovae for an accelerating universe and a cosmological constant. Astron. J. 1998, 116, 1009–1038. [Google Scholar] [CrossRef] [Green Version]
  27. Perlmutter, S.; Aldering, G.; Goldhaber, G.; Knop, R.A.; Nugent, P.; Castro, P.G.; Deustua, S.; Fabbro, S.; Goobar, A.; Groom, D.E.; et al. Measurements of Omega and Lambda from 42 High-Redshift Supernovae. Astrophys. J. 1999, 517, 565–586. [Google Scholar] [CrossRef]
  28. Ivashchuk, V.D.; Kobtsev, A.A. Exact exponential cosmological solutions with two factor spaces of dimension m in EGB model with a Λ-term. Int. J. Geom. Meth. Mod. Phys. 2019, 16, 1950025. [Google Scholar] [CrossRef]
  29. Ivashchuk, V.D. On Stability of Exponential Cosmological Type Solutions with Two Factor Spaces in the Einstein–Gauss–Bonnet Model with a Λ-term. Gravit. Cosmol. 2020, 20, 16–21. [Google Scholar] [CrossRef]
  30. Glavan, D.; Lin, C. Einstein–Gauss–Bonnet Gravity in Four-Dimensional Spacetime. Phys. Rev. Lett. 2020, 124, 081301. [Google Scholar] [CrossRef] [Green Version]
  31. Fernandes, P.G.S.; Carrilho, P.; Clifton, T.; Mulryne, D.J. The 4D Einstein–Gauss–Bonnet theory of gravity: A review. Class. Quant. Grav. 2022, 39, 063001. [Google Scholar] [CrossRef]
  32. Ivashchuk, V.D.; Kobtsev, A.A. Exponential cosmological solutions with two factor spaces in EGB model with Λ = 0 revisited. Eur. Phys. J. C 2019, 79, 824. [Google Scholar] [CrossRef]
  33. Ernazarov, K.K.; Ivashchuk, V.D. Examples of Stable Exponential Cosmological Solutions with Three Factor Spaces in EGB Model with a Λ-Term. Grav. Cosmol. 2019, 25, 164–168. [Google Scholar] [CrossRef] [Green Version]
  34. Ernazarov, K.K.; Ivashchuk, V.D. Stable exponential cosmological solutions with three different Hubble-like parameters in (1 + 3 + k1 + k2)-dimensional EGB model with a Λ-term. Symmetry 2020, 12, 250. [Google Scholar] [CrossRef] [Green Version]
Figure 1. The graphical representation of the function λ = f ( x 1 ) for m = 4 , k 1 = 6 , k 2 = 7 [1].
Figure 1. The graphical representation of the function λ = f ( x 1 ) for m = 4 , k 1 = 6 , k 2 = 7 [1].
Symmetry 14 01296 g001
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Ernazarov, K.K.; Ivashchuk, V.D. Stable Exponential Cosmological Type Solutions with Three Factor Spaces in EGB Model with a Λ-Term. Symmetry 2022, 14, 1296. https://doi.org/10.3390/sym14071296

AMA Style

Ernazarov KK, Ivashchuk VD. Stable Exponential Cosmological Type Solutions with Three Factor Spaces in EGB Model with a Λ-Term. Symmetry. 2022; 14(7):1296. https://doi.org/10.3390/sym14071296

Chicago/Turabian Style

Ernazarov, Kubantai K., and Vladimir D. Ivashchuk. 2022. "Stable Exponential Cosmological Type Solutions with Three Factor Spaces in EGB Model with a Λ-Term" Symmetry 14, no. 7: 1296. https://doi.org/10.3390/sym14071296

APA Style

Ernazarov, K. K., & Ivashchuk, V. D. (2022). Stable Exponential Cosmological Type Solutions with Three Factor Spaces in EGB Model with a Λ-Term. Symmetry, 14(7), 1296. https://doi.org/10.3390/sym14071296

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