Next Article in Journal
Minisuperspace Quantization of f(T, B) Cosmology
Next Article in Special Issue
On A. D. Sakharov’s Hypothesis of Cosmological Transitions with Changes in the Signature of the Metric
Previous Article in Journal
Renormalizability of Alternative Theories of Gravity: Differences between Power Counting and Entropy Argument
Previous Article in Special Issue
The Dirac Sea, T and C Symmetry Breaking, and the Spinor Vacuum of the Universe
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

De Sitter Solutions in Einstein–Gauss–Bonnet Gravity †

Skobeltsyn Institute of Nuclear Physics, Lomonosov Moscow State University, Leninskie Gory 1, 119991 Moscow, Russia
*
Author to whom correspondence should be addressed.
This paper is an extended version from the proceeding paper: Sergey Vernov and Ekaterina Pozdeeva. De Sitter Solutions in Models with the Gauss-Bonnet Term. In Proceedings of the 1st Electronic Conference on Universe, online, 22–28 February 2021.
Submission received: 23 April 2021 / Revised: 12 May 2021 / Accepted: 12 May 2021 / Published: 16 May 2021

Abstract

:
De Sitter solutions play an important role in cosmology because the knowledge of unstable de Sitter solutions can be useful to describe inflation, whereas stable de Sitter solutions are often used in models of late-time acceleration of the Universe. The Einstein–Gauss–Bonnet gravity cosmological models are actively used both as inflationary models and as dark energy models. To modify the Einstein equations one can add a nonlinear function of the Gauss–Bonnet term or a function of the scalar field multiplied on the Gauss–Bonnet term. The effective potential method essentially simplifies the search and stability analysis of de Sitter solutions, because the stable de Sitter solutions correspond to minima of the effective potential.
PACS:
98.80.-k; 04.50.Kd

1. Introduction

Cosmological models with scalar fields play a central role in the description of the global evolution of the Universe. In particular, modified gravity models with the Ricci scalar multiplied by a function of the scalar field are very popular [1,2,3,4]. These models are quite natural because quantum corrections to the effective action with a minimally coupled scalar field include nonminimal coupling terms [5,6,7]. Many inflationary models that connect cosmology and particle physics include nonminimally coupled scalar fields [8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24].
It is well-known that one can add the Gauss–Bonnet term to the Hilbert–Einstein Lagrangian of the General Relativity and it does not change the equations of motion. On the other hand, this term multiplied by some nonconstant function of a scalar field modifies the equations of motion. Models with both the Ricci scalar and the Gauss–Bonnet term multiplied by some functions of the scalar field are natural generalizations of the models with a minimal coupling [25,26,27]. Furthermore, models with a nonlinear function of the Gauss–Bonnet term can be rewritten in the equivalent form that includes a scalar field without kinetic term [28,29,30,31].
The cosmological models with the Gauss–Bonnet term are motivated by the string theory [30,31,32,33,34,35,36,37] and are actively used for describing of both the early Universe evolution (inflation) [25,26,27,38,39,40,41,42,43,44,45,46,47,48,49,50,51,52,53,54,55,56,57] and the current dark energy dominated epoch [28,29,30,31,36,37,58,59,60,61,62,63,64,65,66]. Note that both stages of the Universe evolution are characterized by the quasi de Sitter accelerated expansion of the Universe. So, it is important to have an effective method for the searching of de Sitter solutions and the study of their stability. The method proposed in [67] solves this problem for the Gauss–Bonnet model with the standard scalar field. It is a generalization of the effective potential method for models with scalar field nonminimally coupled to the curvature only [68,69,70]. We have shown in Ref. [27] that the effective potential is a useful tool to generalize the known inflationary models with the Gauss–Bonnet term. On the one hand, the effective potential allows to rule out models with a stable de Sitter solution, for which it is difficult to construct an inflationary scenario with a graceful exit. Models with unstable de Sitter solutions or without an exact de Sitter solution are more suitable. On the other hand, the scalar spectral index n s and the amplitude of the scalar perturbations A s as functions of the e-folding number can be expressed via derivatives of the effective potential [27,55,57].
In this paper, we generalize the effective potential method on models with nonlinear functions of the Gauss–Bonnet term G that can be presented as models with a scalar field without kinetic term. We also show that the situation is more difficult in the case of a phantom scalar field.
The paper is organized as follows. In Section 2, we remind the evolution equations of the model considered and the standard way for the search of de Sitter solutions. In Section 3, we analyze the stability of de Sitter solutions due to the effective potential. A way of constructing different models with the same structure of de Sitter solutions is proposed in Section 4. In Section 5, we analyze the stability of de Sitter solutions in the proposed type of L ( R , G ) gravity models. Section 6 is devoted to our conclusions.

2. Models the Gauss–Bonnet Term

Let us consider the model with the Gauss–Bonnet term described by the following action:
S = d 4 x g U R c 2 g μ ν μ ϕ ν ϕ V F G ,
where the functions U ( ϕ ) , V ( ϕ ) , and F ( ϕ ) are double differentiable ones, c is a constant, R is the Ricci scalar and G is the Gauss–Bonnet term,
G = R 2 4 R μ ν R μ ν + R μ ν α β R μ ν α β .
As known [29,31], the action
S = d 4 x g U 0 R + L ( G ) ,
where U 0 is a constant and L ( G ) is a double differentiable function, can be rewritten in the form of action (1) with c = 0 :
S = d 4 x g U 0 R + L ( ϕ ) ( G ϕ ) + L ( ϕ ) ,
where a prime denotes the derivatives with respect to ϕ . Varying action (3) over ϕ , one gets ϕ = G and the initial L ( G ) model.
In the spatially flat Friedmann–Lemaître–Robertson–Walker metric with the interval
d s 2 = d t 2 + a 2 ( t ) d x 1 2 + d x 2 2 + d x 3 2 ,
where a ( t ) is the scale factor, one gets the following evolution equations [56]:
6 H 2 U + 6 H U ϕ ˙ = c 2 ϕ ˙ 2 + V + 24 H 3 F ϕ ˙ ,
4 U 4 H F ˙ H ˙ = c ϕ ˙ 2 2 U ¨ + 2 H U ˙ + 8 H 2 F ¨ H F ˙ ,
c ϕ ¨ + 3 c H ϕ ˙ 6 H ˙ + 2 H 2 U + V + 24 H 2 F H ˙ + H 2 = 0 ,
where H = a ˙ / a is the Hubble parameter, dots and primes denote the derivatives with respect to the cosmic time and the scalar field ϕ , respectively. At c = 1 , these equations have been investigated in many papers (see, for example [25,56,67]).
To find de Sitter solutions with a constant ϕ in the model (2) we substitute ϕ = ϕ d S and H = H d S into Equations (5) and (7). A de Sitter solution does not depend on the value of c, so we obtain the same results as in the case of c = 1 considered in [67]:
H d S 2 = V d S 6 U d S
and
F d S = 3 U d S ( 2 U d S V d S V d S U d S ) 2 V d S 2 ,
where A d S A ( ϕ d S ) for any function A. Therefore, for arbitrary functions U ( ϕ ) and V ( ϕ ) with V d S U d S > 0 , we can choose F ( ϕ ) such that the corresponding point becomes a de Sitter solution, with the Hubble parameter defined by Equation (8). We always choose that H d S > 0 .

3. Stability of de Sitter Solutions

To analyze the stability of a de Sitter solution we transform Equations (6) and (7) to the following dynamical system:
ϕ ˙ = ψ , ψ ˙ = 1 2 B ˜ 4 c F H ψ 2 H 3 B + 4 F V 6 U 2 6 c U ψ 2 V 2 U X + 12 H 2 2 U + 3 c F + 2 U F 96 F F H 4 3 2 U + c U ψ 2 , H ˙ = 1 4 B ˜ 4 c F H ψ 8 c U 4 F H 2 H ψ 2 V 2 U 2 4 F H 2 U X + 8 F H 2 2 U c c ψ 2 ,
where
B ˜ = 3 4 H 2 F U 2 + c U ,
X = U 2 V 2 24 H 4 F 12 H 2 U + V .
In the case c = 0 , the last equation is essentially simplified:
H ˙ = 24 H 4 F 12 H 2 U + V 6 U 4 H 2 F .
At a de Sitter point system (10) is
ϕ ˙ = 0 , ψ ˙ = 0 , H ˙ = 0 ,
that corresponds to X d S = 0 .
In Ref. [67], the effective potential has been proposed for models with the Gauss–Bonnet term:
V e f f = U 2 V + 2 3 F .
Using Equations (8) and (12), we obtain
X d S = 2 3 F d S 2 U d S U d S V d S + V d S U d S 2 V d S 2 = V e f f ( ϕ d S ) = 0 ,
therefore, de Sitter solutions correspond to extremum points of the effective potential V e f f .
To investigate the Lyapunov stability of a de Sitter solution we use the following expansions:
H ( t ) = H d S + ε H 1 ( t ) , ϕ ( t ) = ϕ d S + ε ϕ 1 ( t ) , ψ ( t ) = ε ψ 1 ( t ) ,
where ε is a small parameter. Therefore,
X = ε ( X , H H 1 + X , ϕ ϕ 1 ) + O ( ε 2 ) ,
where
X , H = X H ϕ = ϕ d S = 4 6 V d S 5 / 2 U d S 3 / 2 U d S V d S V d S U d S ,
X , ϕ = X ϕ ϕ = ϕ d S = 1 V d S 2 2 3 F d S V d S 2 2 U d S U d S V d S + V d S U d S 2 .
The functions H 1 ( t ) , ϕ 1 ( t ) , and ψ 1 ( t ) are connected by Equation (5):
H 1 ( t ) = V d S U d S U d S V d S 2 U d S V d S H d S ϕ 1 ( t ) ψ 1 ( t ) .
This expression does not depend on the value of c and coincides with the corresponding expression obtained in Ref. [67].
Substituting (16)–(18) into Equation (10) in the first order of ε , we obtain the following system of two linear differential equations:
ϕ ˙ 1 = ψ 1 ,
ψ ˙ 1 = 2 F d S V d S 3 6 U d S U d S V d S 2 + 3 V d S U d S 2 V d S 6 U d S V d S V d S U d S 2 3 U d S V d S B d S ϕ 1 6 U d S V d S 2 U d S ψ 1 ,
where
B d S = 3 V d S 2 V d S U d S U d S V d S 2 + c U d S .
This system can be rewritten in the matrix form:
ϕ ˙ 1 ψ ˙ 1 = A 11 A 21 A 12 A 22 ϕ 1 ψ 1
where the matrix
A = 0 , 1 V d S 2 V e f f ( ϕ d S ) U d S B d S , 3 H d S
The general solution of system (22) has the following form:
ϕ 1 = c 11 e λ t + c 21 e λ + t ,
ψ 1 = c 21 e λ t + c 22 e λ + t ,
where c i j are some constants. Solving the characteristic equation:
det ( A ˜ λ · I ) = λ 2 3 H d S λ + V d S 2 V e f f ( ϕ d S ) U d S B d S = 0 ,
we get the following roots:
λ ± = 3 2 H d S ± 9 4 H d S 2 V d S 2 U d S B d S V e f f ( ϕ d S ) .
A de Sitter solution is stable if real parts of both λ and λ + are negative. We consider the case
H d S = V d S 6 U d S > 0 ,
hence, e ( λ ) < 0 .
In the case of a positive U d S , we see that B d S > 0 for c 0 and the condition e ( λ + ) < 0 is equivalent to V e f f ( ϕ d S ) > 0 . In the cases c > 0 and c = 0 , a de Sitter solution is stable if V e f f ( ϕ d S ) > 0 and unstable if V e f f ( ϕ d S ) < 0 .
In the case of c < 0 , we see that B d S can be negative. So, in this case de Sitter solution is stable if the V e f f ( ϕ d S ) B d S > 0 . So, the main result of Ref. [67] can be generalized on the case c = 0 without any correction, whereas the condition should be change to V e f f ( ϕ d S ) B d S > 0 in the case of c < 0 that corresponds to a phantom scalar field ϕ .

4. Different Models with the Same Structure of de Sitter Solutions

It is evident that any change of functions U and V can be compensated by the change of function F such that the first derivative of the effective potential does not change. This property can be used for the generalization of the inflationary scenarios [27]. Let us consider a more nontrivial question: how one can change the functions U and V only to obtain the model with the same structure of de Sitter solutions. To be concrete we seek for models with the same values of ϕ d S and the same stability properties of de Sitter solutions, but maybe with different values of H d S .
Considering the stability analysis of de Sitter solutions, we define the effective potential as such a function that its minima correspond to the stable de Sitter solutions and maxima correspond to unstable de Sitter solutions. For this reason, the effective potential is not unique. We can add a constant to it or multiply it on a positive number. If we transform functions U ( ϕ ) and V ( ϕ ) to functions U ˜ ( ϕ ) = f ( ϕ ) U ( ϕ ) and
V ˜ ( ϕ ) = V ( ϕ ) f ( ϕ ) 2 U ( ϕ ) 2 W ( ϕ ) V ( ϕ ) + U ( ϕ ) 2
then the effective potentials of the original and transformed models are connected as follows:
V ˜ e f f = V e f f + W .
Therefore, the structure of de Sitter solutions does not change if the function W = C V e f f + W 0 , where a constant C > 1 and W 0 is an arbitrary constant. Furthermore, one should check that V ˜ ( ϕ d S ) > 0 and U ˜ ( ϕ d S ) > 0 for all de Sitter solutions.
Moreover, if V e f f ( ϕ ) > 0 for any ϕ , then the functions V e f f ( ϕ ) n and 1 / V e f f ( ϕ ) n , where n is a natural number, can be considered as new effective potentials. Different transformations of the effective potential can be combined, for example, functions 1 / ( V e f f + W ) and ( V e f f + W ) 2 can be considered as effective potentials if V e f f > W for any ϕ . So, transformations of the model with V ˜ e f f = ( V e f f + W ) 2 or V ˜ e f f = 1 / ( V e f f + W ) does not change the structure of de Sitter solutions.

5. Examples of L ( R , G ) Models

5.1. Evolution Equations

Several examples of models with an ordinary scalar field coupled with the Gauss–Bonnet term have been considered in Refs. [55,67]. In this paper, we consider examples in the case of c = 0 , namely, the case of L ( R , G ) model, described by the following action:
S = d 4 x g U 0 R + L ( C 1 R + C 2 G ) ,
where L is a double differentiable function, U 0 , C 1 , and C 2 are constants. A linear function L corresponds to the General Relativity, whereas a nonlinear function L corresponds to the modified gravity. Note that L ( R ) and L ( G ) models are particular cases of the model considered.
In the case of a nonlinear function L , action (29) can be rewritten in the following form:
S = d 4 x g U 0 R + L ( ϕ ) ( C 1 R + C 2 G ϕ ) + L ( ϕ ) .
Varying action (30) over ϕ , one gets ϕ = C 1 R + C 2 G and the initial L ( R , G ) model with action (29). Action (30) is a particular case of action (1) with the functions:
U = U 0 + C 1 L , V = ϕ L L , F = C 2 L , c = 0 .
The effective potential is the following combination of the function L and its first derivative:
V e f f = 2 3 C 2 L U 0 + C 1 L 2 ϕ L L .
So,
V e f f = L ϕ L L 2 ( 3 C 1 2 2 C 2 ϕ ) ϕ L 2 L L 2 C 2 L 2 6 C 1 U 0 L 3 U 0 2 ϕ
and a point with L = 0 correspond to de Sitter solutions if ϕ L > L and U 0 + C 1 L > 0 . Furthermore, the second multiplied can be equal to zero that also can correspond to a de Sitter solution. We explorer de Sitter solutions in detail in a few examples of L ( R , G ) models.

5.2. The Function F in a Role of the Effective Potential

If V = C U 2 , where C is a constant, then V e f f = 2 F / 3 = 2 C 2 L / 3 and the function L plays the role of the effective potential. For the considering L ( R , G ) models the condition V = C U 2 is the following first order differential equation
C C 1 2 L 2 + ( 2 C C 1 U 0 ϕ ) L + C U 0 2 + L = 0 .
This equation has two solutions:
L 1 = A ϕ C ( C 1 A + U 0 ) 2 , L 2 = 1 4 C C 1 2 ϕ 2 U 0 C 1 ϕ ,
where A is an integration constant. The function L 1 is a linear one, so this case is the General Relativity model with the cosmological constant. In the case of the function L 2 , the function F is a linear one and this model has no de Sitter solution.

5.3. The Case of a Power Function L of the Gauss–Bonnet Term

Let L = C G α , where C and α are constants. Substituting this function with C 1 = 0 and C 2 = 1 into Equation (33), we get
V e f f = 2 α 2 C 2 ( α 1 ) 2 ϕ 2 α + 3 U 0 2 ϕ 6 ( C ( α 1 ) ϕ α + 2 ) ,
There is no de Sitter solutions at ϕ d S > 0 . So, we consider only such values of α that ϕ α is real and negative at ϕ < 0 and obtain
ϕ d S = 3 U 0 2 2 C 2 ( α 1 ) 2 1 2 α 1 .
Let us consider several examples:
1.
At α = 2 , we get de Sitter solution
ϕ d S = 3 U 0 2 2 C 2 1 / 3
and the potential V = C ϕ 2 , so, C > 0 is a necessary condition for de Sitter solutions existence.
The effective potential
V e f f = 4 C ϕ 3 U 0 C ϕ 2 ,
and its second derivative is negative: V e f f = 6 U 0 C ϕ 4 . Therefore, we can conclude that the de Sitter solution is unique and unstable.
2.
At α = 3 we get
ϕ d S = 3 U 0 2 8 C 2 1 / 5 .
The potential
V = 2 C ϕ 3 .
The condition V ( ϕ d S ) > 0 demands C < 0 for ϕ d S > 0 . Using
V e f f ( ϕ d S ) = 20 C ,
we get that the considered model with α = 3 and C < 0 has an unstable de Sitter solution.
3.
The case of α = 1 / 3 . Similar models have been proposed in [31]. For
ϕ d S = 512 C 6 19683 U 0 6
we obtain
V ( ϕ d S ) = 16 C 3 81 U 0 2 and V e f f ( ϕ d S ) = 3486784401 U 0 16 4194304 C 15 .
So, a de Sitter solution exist if and only if C > 0 and it is unstable.

5.4. The Case of a Quadratic Polynomial L

5.4.1. Equation for ϕ d S

Let us consider the case of
L = b 2 ϕ 2 + b 1 ϕ + b 0 ,
where b i are constants, b 2 0 . In this case,
U = 2 C 1 b 2 ϕ + Q 0 , V = b 2 ϕ 2 b 0 , F = C 2 ( 2 b 2 ϕ + b 1 ) ,
where Q 0 = U 0 + b 1 C 1 . Note that the function F is defined up to a constant, so we can put b 1 = 0 without loss of generality.
De Sitter solution corresponds to ϕ d S that is a solution of the following equation:
2 C 2 b 2 2 ϕ d S 4 + 2 ( 3 C 1 b 2 Q 0 2 C 2 b 0 b 2 ) ϕ d S 2 + 3 ( 4 C 1 2 b 0 b 2 + Q 0 2 ) ϕ d S + 6 C 1 b 0 Q 0 + 2 C 2 b 0 2 = 0 .
Let us consider a few interesting particular cases of this equation.

5.4.2. The Case L ( R ) Gravity

If C 2 = 0 , then the model has no Gauss–Bonnet term and is a L ( R ) gravity model with the effective potential
V e f f = 2 C 1 b 2 ϕ + U ˜ 0 2 b 0 b 2 ϕ 2 .
Equation (46) is a quadratic one and has the following two solutions:
ϕ 1 = 2 C 1 b 0 U ˜ 0 , ϕ 2 = U ˜ 0 2 C 1 b 2 .
The point ϕ = ϕ 2 does not correspond to de Sitter solution because of U ( ϕ 2 ) = 0 . At ϕ = ϕ 1 , we get
V e f f ϕ = ϕ 1 = 2 b 2 U ˜ 0 4 b 0 2 4 C 1 2 b 0 b 2 U ˜ 0 2
At de Sitter point ϕ = ϕ 1 , we obtain
U d S = U ˜ 0 2 4 C 1 2 b 0 b 2 U ˜ 0 , V d S = 4 b 2 C 1 2 b 0 2 U ˜ 0 2 b 0 .
From Equation (8), we get
H d S 2 = b 0 6 U ˜ 0
a de Sitter solution exists only if b 0 U ˜ 0 < 0 . Note that the Starobinsky R 2 inflationary model [71,72,73] does not include the cosmological constant, so, b 0 = 0 . In this case, a de Sitter solution does not exist. The potential V > 0 for all values of ϕ if b 2 > 0 and b 0 < 0 . In this case, the de Sitter solution is stable.

5.4.3. The Case L ( G ) Gravity

If C 1 = 0 , then the model is a L ( G ) gravity models with the effective potential
V e f f = 2 b 2 2 C 2 b 2 2 ϕ 4 4 C 2 b 0 b 2 ϕ 2 + 3 U 0 2 ϕ + 2 C 2 b 0 2 3 ( b 2 ϕ 2 b 0 ) 2
The value of ϕ at a de Sitter point is a real solution of the following equation:
2 C 2 b 0 b 2 ϕ d S 2 2 + 3 U 0 2 ϕ d S = 0 .
In the case of b 0 = 0 , one gets the following real solutions of Equation (53):
ϕ 1 = 0 , ϕ 2 = 12 U 0 2 C 2 2 b 2 1 / 3 2 C 2 b 2 .
Note that ϕ 1 does not correspond to a de Sitter solution, because V ( ϕ 1 ) = 0 . So, we get the unique de Sitter solution ϕ d S = ϕ 2 .
At the de Sitter point, the second derivative of the effective potential is
V d S = 2 3 144 b 2 5 C 2 4 U 0 2 1 / 3 .
Furthermore, we demand that V d S = b 2 ϕ d s 2 > 0 , so, b 2 > 0 and we get a model with one unstable de Sitter solution.

5.4.4. The Case of the Absence of the Cosmological Constant

In the case of b 0 = 0 , the potential V = b 2 ϕ 2 , therefore, ϕ d S 0 and b 2 > 0 are necessary condition for de Sitter solution existence. Assuming that Q 0 0 , we obtain the first derivative of the effective potential in the following form:
V e f f = 2 Q 0 2 2 B 1 ϕ 3 + 6 B 2 ϕ + 3 3 b 2 ϕ 3 ,
where
B 1 b 2 Q 0 C 1 , B 2 = b 2 2 Q 0 2 C 2 .
The effective potential can be multiplied on any positive constant, so we can consider the function
J ( ϕ ) = 2 B 1 ϕ 3 + 3 ( 2 B 2 ϕ + 1 ) ϕ 3
instead of V e f f and the number of its zeros depends on values of parameters B 1 and B 2 only.
The function
U = Q 0 ( 2 B 2 ϕ + 1 )
is not always positive, so it is possible that an extremum of the effective potential does not corresponds to de Sitter solution. The sign of U ( ϕ d S ) depends on the sign of the parameter Q 0 , whereas V e f f does not depend on the sign of the parameter Q 0 . Note that the functions U and J are not equal to zero at the same point, because B 1 0 .
The same values of parameters B 1 and B 2 correspond both to models with one de Sitter solution and to models without de Sitter solution in dependence of the sign of Q 0 . A decreasing behavior of the function J in the neighborhood of J = 0 corresponds to an unstable de Sitter solution. In Figure 1, blue solid curves and green dash-dot curves correspond to models either with one unstable de Sitter solution or without de Sitter solutions in dependence of the sign of U ( ϕ d S ) .
In Figure 2, we present the function J ( ϕ ) that has three roots (red solid curves) and the function U ( ϕ ) for the same values of parameters B 1 and B 2 and the parameter Q 0 = ± 1 . One can see that the model has one unstable de Sitter solution or one stable and one unstable de Sitter solutions in dependence on the sign of Q 0 . The crosses of gray and black lines with axis ϕ correspond to switch of gravity and antigravity regimes in the model considered.

6. Conclusions

In this paper, we consider de Sitter solutions in models with the Gauss–Bonnet term, including L ( R , G ) gravity models. We investigate evolution equations of a scalar field nonminimally coupled both with the curvature and with the Gauss–Bonnet term and look for the fixed points of scalar field dynamics which correspond to de Sitter solutions. We show that, in the case of a positive coupling function U ( ϕ ) , it is possible to introduce an effective potential V e f f which can be expressed through the function U of nonminimal coupling with the curvature, the scalar field potential V, and the coupling function with the Gauss–Bonnet term denoted by F. We show that it is convenient to investigate the structure of fixed points using the effective potential because the stable de Sitter solutions correspond to minima of the effective potential. The existence and stability of de Sitter solutions in the system under consideration can be studied with the help of function V e f f since one can get the structure and stability properties of de Sitter solutions using a graphical representation of the effective potential only. It should be noted that the effective potential is not uniquely defined. We can multiply it on a positive constant or add a constant to it. If the effective potential is a positive definite function, then we can consider V e f f n and 1 / V e f f n , where n is a natural number, like other forms of the effective potential.
In this paper, we show that the effective potential proposed [67] for models with the Gauss–Bonnet term multiplied on a function of the scalar field can be used in L ( R , G ) models as well. To find de Sitter solutions in some L ( R , G ) model, we rewrite the action of this model in the form (30) and construct the corresponding effective potential V e f f . A stable de Sitter solutions corresponds V e f f ( ϕ d S ) > 0 , where the values of the scalar field at de Sitter point ϕ d S is determined by the condition V e f f ( ϕ d S ) = 0 . We have found de Sitter solutions in a few L ( R , G ) models to demonstrate the effective potential method.
Note that the proposed effective potential is a useful tool for the construction of inflationary scenarios in the models with the Gauss–Bonnet term multiplied to a function of the scalar field [27,55,57]. It is interesting that in the slow-roll approximation the scalar spectral index n s and the amplitude of the scalar perturbations A s as functions of the e-folding number can be expressed via derivatives of the effective potential given in the form (14). Furthermore, the knowledge of unstable de Sitter solutions can be useful to describe inflation (see, for example, Ref. [20]). We plan to generalize this approach to inflationary scenarios in L ( R , G ) models.
The search for stable de Sitter solutions is important for f ( R , G ) models that explain the late-time accelerated expansion of the Universe [29,59,60,62,65,66]. A generic f ( R , G ) action can be transformed into one linear in R and G by including of two scalar fields, whereas the proposed special type of such models describing by action (29) can be linearized in R and G by including of one scalar field without kinetic term. It allows to use the effective potential method and to simplify analysis of the stability of de Sitter solutions in distinguish to the traditional approach [62]. The proposed L ( R , G ) models include not only F ( R ) and F ( G ) gravity models, but also more complicated models with the L ( R + G ) function in the action. We plan to investigate a possibility to describe the late-time accelerated expansion of the Universe in such types of models taking into account the observation restrictions.

Author Contributions

The authors contributed equally to this work. Investigation, S.V. and E.P.; Writing original draft, S.V. and E.P. All authors have read and agreed to the published version of the manuscript.

Funding

This work is partially supported by the Russian Foundation for Basic Research grant No. 20-02-00411.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Capozziello, S.; Faraoni, V. Beyond Einstein Gravity: A Survey of Gravitational Theories for Cosmology and Astrophysics; Springer: New York, NY, USA, 2011. [Google Scholar]
  2. Capozziello, S.; De Laurentis, M. Extended Theories of Gravity. Phys. Rep. 2011, 509, 167–321. [Google Scholar] [CrossRef] [Green Version]
  3. Fujii, Y.; Maeda, K. The Scalar—Tensor Theory of Gravitation; Cambridge University Press: Cambridge, UK, 2004. [Google Scholar]
  4. Faraoni, V. Cosmology in Scalar—Tensor Gravity; Kluwer Academic: Dordrecht, The Netherlands, 2004. [Google Scholar]
  5. Chernikov, E.A.; Tagirov, E.A. Quantum theory of scalar fields in de Sitter space-time. Ann. Poincare Phys. Theor. A 1968, 9, 109. [Google Scholar]
  6. Tagirov, E.A. Consequences of field quantization in de Sitter type cosmological models. Ann. Phys. 1973, 76, 561. [Google Scholar] [CrossRef]
  7. Callan, C.G.; Coleman, S.R.; Jackiw, R. A New improved energy—Momentum tensor. Ann. Phys. 1970, 59, 42. [Google Scholar] [CrossRef]
  8. Barvinsky, A.O.; Kamenshchik, A.Y. Quantum scale of inflation and particle physics of the early universe. Phys. Lett. B 1994, 332, 270. [Google Scholar] [CrossRef] [Green Version]
  9. Bezrukov, F.L.; Shaposhnikov, M. The Standard Model Higgs boson as the inflaton. Phys. Lett. B 2008, 659, 703. [Google Scholar] [CrossRef] [Green Version]
  10. Barvinsky, A.O.; Kamenshchik, A.Y.; Starobinsky, A.A. Inflation scenario via the Standard Model Higgs boson and LHC. J. Cosmol. Asropart. Phys. 2008, 2008, 021. [Google Scholar] [CrossRef]
  11. De Simone, A.; Hertzberg, M.P.; Wilczek, F. Running Inflation in the Standard Model. Phys. Lett. B 2009, 678, 1. [Google Scholar] [CrossRef] [Green Version]
  12. Barvinsky, A.O.; Kamenshchik, A.Y.; Kiefer, C.; Starobinsky, A.A.; Steinwachs, C.F. Higgs boson, renormalization group, and cosmology. Eur. Phys. J. C 2012, 72, 2219. [Google Scholar] [CrossRef]
  13. Lerner, R.N.; McDonald, J. Higgs Inflation and Naturalness. J. Cosmol. Astropart. Phys. 2010, 1004, 015. [Google Scholar] [CrossRef] [Green Version]
  14. Bezrukov, F.L.; Magnin, A.; Shaposhnikov, M.; Sibiryakov, S. Higgs inflation: Consistency and generalisations. J. High Energy Phys. 2011, 1101, 016. [Google Scholar] [CrossRef] [Green Version]
  15. Greenwood, R.N.; Kaiser, D.I.; Sfakianakis, E.I. Multifield Dynamics of Higgs Inflation. Phys. Rev. D 2013, 87, 064021. [Google Scholar] [CrossRef] [Green Version]
  16. Bezrukov, F.L. The Higgs field as an inflaton. Class. Quant. Grav. 2013, 30, 214001. [Google Scholar] [CrossRef]
  17. Cerioni, A.; Finelli, F.; Tronconi, A.; Venturi, G. Inflation and Reheating in Induced Gravity. Phys. Lett. B 2009, 681, 383. [Google Scholar] [CrossRef] [Green Version]
  18. Kallosh, R.; Linde, A.; Roest, D. The double attractor behavior of induced inflation. J. High Energy Phys. 2014, 1409, 062. [Google Scholar] [CrossRef] [Green Version]
  19. Rinaldi, M.; Vanzo, L.; Zerbini, S.; Venturi, G. Inflationary quasi-scale invariant attractors. Phys. Rev. D 2016, 93, 024040. [Google Scholar] [CrossRef] [Green Version]
  20. Elizalde, E.; Odintsov, S.D.; Pozdeeva, E.O.; Vernov, S.Y. Renormalization-group inflationary scalar electrodynamics and SU(5) scenarios confronted with Planck2013 and BICEP2 results. Phys. Rev. D 2014, 90, 084001. [Google Scholar] [CrossRef] [Green Version]
  21. Elizalde, E.; Odintsov, S.D.; Pozdeeva, E.O.; Vernov, S.Y. Cosmological attractor inflation from the RG-improved Higgs sector of finite gauge theory. J. Cosmol. Astropart. Phys. 2016, 1602, 025. [Google Scholar] [CrossRef] [Green Version]
  22. Pozdeeva, E.O.; Vernov, S.Y. Renormalization-group improved inflationary scenarios. Phys. Part. Nucl. Lett. 2017, 14, 386. [Google Scholar] [CrossRef] [Green Version]
  23. Dubinin, M.N.; Petrova, E.Y.; Pozdeeva, E.O.; Sumin, M.V.; Vernov, S.Y. MSSM-inspired multifield inflation. J. High Energy Phys. 2017, 1712, 036. [Google Scholar] [CrossRef] [Green Version]
  24. Kamenshchik, A.Y.; Tronconi, A.; Venturi, G. Quantum cosmology and the inflationary spectra from a nonminimally coupled inflaton. Phys. Rev. D 2020, 101, 023534. [Google Scholar] [CrossRef] [Green Version]
  25. van de Bruck, C.; Longden, C. Higgs Inflation with a Gauss–Bonnet term in the Jordan Frame. Phys. Rev. D 2016, 93, 063519. [Google Scholar] [CrossRef] [Green Version]
  26. Mathew, J.; Shankaranarayanan, S. Low scale Higgs inflation with Gauss–Bonnet coupling. Astropart. Phys. 2016, 84, 1. [Google Scholar] [CrossRef] [Green Version]
  27. Pozdeeva, E.O.; Vernov, S.Y. Construction of inflationary scenarios with the Gauss–Bonnet term and nonminimal coupling. arXiv 2021, arXiv:2104.04995. [Google Scholar]
  28. Nojiri, S.; Odintsov, S.D.; Sasaki, M. Gauss–Bonnet dark energy. Phys. Rev. D 2005, 71, 123509. [Google Scholar] [CrossRef] [Green Version]
  29. Cognola, G.; Elizalde, E.; Nojiri, S.; Odintsov, S.D.; Zerbini, S. Dark energy in modified Gauss–Bonnet gravity: Late-time acceleration and the hierarchy problem. Phys. Rev. D 2006, 73, 084007. [Google Scholar] [CrossRef] [Green Version]
  30. Cognola, G.; Elizalde, E.; Nojiri, S.; Odintsov, S.D.; Zerbini, S. String-inspired Gauss–Bonnet gravity reconstructed from the universe expansion history and yielding the transition from matter dominance to dark energy. Phys. Rev. D 2007, 75, 086002. [Google Scholar] [CrossRef] [Green Version]
  31. Nojiri, S.; Odintsov, S.D.; Oikonomou, V.K. Modified Gravity Theories on a Nutshell: Inflation, Bounce and Late-time Evolution. Phys. Rept. 2017, 692, 1–104. [Google Scholar] [CrossRef] [Green Version]
  32. Antoniadis, I.; Rizos, J.; Tamvakis, K. Singularity-free cosmological solutions of the superstring effective action. Nucl. Phys. B 1994, 415, 497. [Google Scholar] [CrossRef] [Green Version]
  33. Kawai, S.; Soda, J. Evolution of fluctuations during graceful exit in string cosmology. Phys. Lett. B 1999, 460, 41. [Google Scholar] [CrossRef] [Green Version]
  34. Cartier, C.; Hwang, J.C.; Copeland, E.J. Evolution of cosmological perturbations in nonsingular string cosmologies. Phys. Rev. D 2001, 64, 103504. [Google Scholar] [CrossRef] [Green Version]
  35. Hwang, J.C.; Noh, H. Classical evolution and quantum generation in generalized gravity theories including string corrections and tachyon: Unified analyses. Phys. Rev. D 2005, 71, 063536. [Google Scholar] [CrossRef] [Green Version]
  36. Calcagni, G.; Tsujikawa, S.; Sami, M. Dark energy and cosmological solutions in second-order string gravity. Class. Quant. Grav. 2005, 22, 3977. [Google Scholar] [CrossRef]
  37. Tsujikawa, S.; Sami, M. String-inspired cosmology: Late time transition from scaling matter era to dark energy universe caused by a Gauss–Bonnet coupling. J. Cosmol. Astropart. Phys. 2007, 0701, 006. [Google Scholar] [CrossRef]
  38. Guo, Z.; Schwarz, D.J. Slow-roll inflation with a Gauss–Bonnet correction. Phys. Rev. D 2010, 81, 123520. [Google Scholar] [CrossRef] [Green Version]
  39. Koh, S.; Lee, B.H.; Lee, W.; Tumurtushaa, G. Observational constraints on slow-roll inflation coupled to a Gauss–Bonnet term. Phys. Rev. D 2014, 90, 063527. [Google Scholar] [CrossRef] [Green Version]
  40. Jiang, P.X.; Hu, J.W.; Guo, Z.K. Inflation coupled to a Gauss–Bonnet term. Phys. Rev. D 2013, 88, 123508. [Google Scholar] [CrossRef] [Green Version]
  41. De Laurentis, M.; Paolella, M.; Capozziello, S. Cosmological inflation in F(R, G ) gravity. Phys. Rev. D 2015, 91, 083531. [Google Scholar] [CrossRef] [Green Version]
  42. Oikonomou, V.K. Autonomous dynamical system approach for inflationary Gauss–Bonnet modified gravity. Int. J. Mod. Phys. D 2018, 27, 1850059. [Google Scholar] [CrossRef] [Green Version]
  43. Wu, Q.; Zhu, T.; Wang, A. Primordial Spectra of slow-roll inflation at second-order with the Gauss–Bonnet correction. Phys. Rev. D 2018, 97, 103502. [Google Scholar] [CrossRef] [Green Version]
  44. Nozari, K.; Rashidi, N. Perturbation, nonGaussianity, and reheating in a Gauss–Bonnet α-attractor model. Phys. Rev. D 2017, 95, 123518. [Google Scholar] [CrossRef] [Green Version]
  45. Koh, S.; Lee, B.H.; Tumurtushaa, G. Reconstruction of the Scalar Field Potential in Inflationary Models with a Gauss–Bonnet term. Phys. Rev. D 2017, 95, 123509. [Google Scholar] [CrossRef] [Green Version]
  46. Chakraborty, S.; Paul, T.; SenGupta, S. Inflation driven by Einstein–Gauss–Bonnet gravity. Phys. Rev. D 2018, 98, 083539. [Google Scholar] [CrossRef] [Green Version]
  47. Yi, Z.; Gong, Y.; Sabir, M. Inflation with Gauss–Bonnet coupling. Phys. Rev. D 2018, 98, 083521. [Google Scholar] [CrossRef] [Green Version]
  48. Odintsov, S.D.; Oikonomou, V.K. Viable Inflation in Scalar–Gauss–Bonnet Gravity and Reconstruction from Observational Indices. Phys. Rev. D 2018, 98, 044039. [Google Scholar] [CrossRef] [Green Version]
  49. Yi, Z.; Gong, Y. Gauss–Bonnet Inflation and the String Swampland. Universe 2019, 5, 200. [Google Scholar] [CrossRef] [Green Version]
  50. Fomin, I.V.; Chervon, S.V. Reconstruction of GR cosmological solutions in modified gravity theories. Phys. Rev. D 2019, 100, 023511. [Google Scholar] [CrossRef] [Green Version]
  51. Odintsov, S.D.; Oikonomou, V.K.; Fronimos, F.P. Rectifying Einstein–Gauss–Bonnet Inflation in View of GW170817. Nucl. Phys. B 2020, 958, 115135. [Google Scholar] [CrossRef]
  52. Odintsov, S.; Oikonomou, V.K. Swampland Implications of GW170817-compatible Einstein–Gauss–Bonnet Gravity. Phys. Lett. B 2020, 805, 135437. [Google Scholar] [CrossRef]
  53. Fomin, I. Gauss–Bonnet term corrections in scalar field cosmology. Eur. Phys. J. C 2020, 80, 1145. [Google Scholar] [CrossRef]
  54. Pozdeeva, E.O. Generalization of cosmological attractor approach to Einstein–Gauss–Bonnet gravity. Eur. Phys. J. C 2020, 80, 612. [Google Scholar] [CrossRef]
  55. Pozdeeva, E.O.; Gangopadhyay, M.R.; Sami, M.; Toporensky, A.V.; Vernov, S.Y. Inflation with a quartic potential in the framework of Einstein–Gauss–Bonnet gravity. Phys. Rev. D 2020, 102, 043525. [Google Scholar] [CrossRef]
  56. Oikonomou, V.K.; Fronimos, F.P. Non-minimally coupled Einstein–Gauss–Bonnet gravity with massless gravitons: The constant-roll case. Eur. Phys. J. Plus 2020, 135, 917. [Google Scholar] [CrossRef]
  57. Pozdeeva, E.O. Violation of the slow-roll regime in the EGB inflationary models with r N e 1 . arXiv 2021, arXiv:2105.02772. [Google Scholar]
  58. Sami, M.; Toporensky, A.V.; Tretjakov, P.V.; Tsujikawa, S. The Fate of (phantom) dark energy universe with string curvature corrections. Phys. Lett. B 2005, 619, 193. [Google Scholar] [CrossRef] [Green Version]
  59. Nojiri, S.; Odintsov, S.D.; Toporensky, A.V.; Tretjakov, P.V. Reconstruction and deceleration-acceleration transitions in modified gravity. Gen. Relativ. Grav. 2010, 42, 1997–2008. [Google Scholar] [CrossRef] [Green Version]
  60. Elizalde, E.; Myrzakulov, R.; Obukhov, V.V.; Saez-Gomez, D. LambdaCDM epoch reconstruction from F(R,G) and modified Gauss–Bonnet gravities. Class. Quant. Grav. 2010, 27, 095007. [Google Scholar] [CrossRef] [Green Version]
  61. Myrzakulov, R.; Saez-Gomez, D.; Tureanu, A. On the ΛCDM Universe in f(G) gravity. Gen. Relativ. Grav. 2011, 43, 1671–1684. [Google Scholar] [CrossRef]
  62. de la Cruz-Dombriz, A.; Saez-Gomez, D. On the stability of the cosmological solutions in f(R,G) gravity. Class. Quant. Grav. 2012, 29, 245014. [Google Scholar] [CrossRef] [Green Version]
  63. Benetti, M.; Santos da Costa, S.; Capozziello, S.; Alcaniz, J.S.; De Laurentis, M. Observational constraints on Gauss–Bonnet cosmology. Int. J. Mod. Phys. D 2018, 27, 1850084. [Google Scholar] [CrossRef] [Green Version]
  64. Navó, G.; Elizalde, E. Stability of hyperbolic and matter-dominated bounce cosmologies from F(R, G ) modified gravity at late evolution stages. Int. J. Geom. Meth. Mod. Phys. 2020, 17, 2050162. [Google Scholar] [CrossRef]
  65. Odintsov, S.D.; Oikonomou, V.K.; Fronimos, F.P.; Fasoulakos, K.V. Unification of a Bounce with a Viable Dark Energy Era in Gauss–Bonnet Gravity. Phys. Rev. D 2020, 102, 104042. [Google Scholar] [CrossRef]
  66. Odintsov, S.D.; Oikonomou, V.K.; Fronimos, F.P. Late-time cosmology of scalar-coupled f(R, G ) gravity. Class. Quant. Grav. 2021, 38, 075009. [Google Scholar] [CrossRef]
  67. Pozdeeva, E.O.; Sami, M.; Toporensky, A.V.; Vernov, S.Y. Stability analysis of de Sitter solutions in models with the Gauss–Bonnet term. Phys. Rev. D 2019, 100, 083527. [Google Scholar] [CrossRef] [Green Version]
  68. Skugoreva, M.A.; Toporensky, A.V.; Vernov, S.Y. Global stability analysis for cosmological models with nonminimally coupled scalar fields. Phys. Rev. D 2014, 90, 064044. [Google Scholar] [CrossRef] [Green Version]
  69. Pozdeeva, E.O.; Skugoreva, M.A.; Toporensky, A.V.; Vernov, S.Y. Possible evolution of a bouncing universe in cosmological models with nonminimally coupled scalar fields. J. Cosmol. Astropart. Phys. 2016, 1612, 006. [Google Scholar] [CrossRef] [Green Version]
  70. Järv, L.; Toporensky, A. Global portraits of nonminimal inflation. arXiv 2021, arXiv:2104.10183. [Google Scholar]
  71. Starobinsky, A.A. A New Type of Isotropic Cosmological Models without Singularity. Phys. Lett. B 1980, 91, 99. [Google Scholar] [CrossRef]
  72. Starobinsky, A.A. Dynamics of phase transition in the new inflationary universe scenario and generation of perturbations. Phys. Lett. B 1982, 117, 175. [Google Scholar] [CrossRef]
  73. Starobinsky, A.A. The Perturbation Spectrum Evolving from a Nonsingular Initially De-Sitter Cosmology and the Microwave Background Anisotropy. Sov. Astron. Lett. 1983, 9, 302. [Google Scholar]
Figure 1. The function J ( ϕ ) for different values of parameters: B 2 = 1 , B 1 = 2 (blue solid curve) and B 1 = 1 (green dash-dot curve) (left); B 2 = 1 , B 1 = 1 / 2 (blue solid curve) and B 1 = 2 (green dash-dot curve) (right).
Figure 1. The function J ( ϕ ) for different values of parameters: B 2 = 1 , B 1 = 2 (blue solid curve) and B 1 = 1 (green dash-dot curve) (left); B 2 = 1 , B 1 = 1 / 2 (blue solid curve) and B 1 = 2 (green dash-dot curve) (right).
Universe 07 00149 g001
Figure 2. The functions J ( ϕ ) and U ( ϕ ) for B 2 = 1 (left) and B 1 = 1 (right). Red solid curves show the function J ( ϕ ) . Black dash curves show the function U ( ϕ ) at Q 0 = 1 and grey dash-dot curves show the function U ( ϕ ) at Q 0 = 1 .
Figure 2. The functions J ( ϕ ) and U ( ϕ ) for B 2 = 1 (left) and B 1 = 1 (right). Red solid curves show the function J ( ϕ ) . Black dash curves show the function U ( ϕ ) at Q 0 = 1 and grey dash-dot curves show the function U ( ϕ ) at Q 0 = 1 .
Universe 07 00149 g002
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Vernov, S.; Pozdeeva, E. De Sitter Solutions in Einstein–Gauss–Bonnet Gravity. Universe 2021, 7, 149. https://0-doi-org.brum.beds.ac.uk/10.3390/universe7050149

AMA Style

Vernov S, Pozdeeva E. De Sitter Solutions in Einstein–Gauss–Bonnet Gravity. Universe. 2021; 7(5):149. https://0-doi-org.brum.beds.ac.uk/10.3390/universe7050149

Chicago/Turabian Style

Vernov, Sergey, and Ekaterina Pozdeeva. 2021. "De Sitter Solutions in Einstein–Gauss–Bonnet Gravity" Universe 7, no. 5: 149. https://0-doi-org.brum.beds.ac.uk/10.3390/universe7050149

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