Next Article in Journal
Ion-Acoustic Rogue Waves in Double Pair Plasma Having Non-Extensive Particles
Next Article in Special Issue
Bouncing Quantum Cosmology
Previous Article in Journal
What Is Matter According to Particle Physics, and Why Try to Observe Its Creation in a Lab?
Previous Article in Special Issue
Bounce Cosmology in Generalized Modified Gravities
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Stability of a Viable Non-Minimal Bounce

Department of Physics, Indian Institute of Technology Madras, Chennai 600036, India
Current address: Department of Physical Sciences, IISER Mohali, Manauli 140306, India.
Universe 2021, 7(3), 62; https://doi.org/10.3390/universe7030062
Submission received: 10 February 2021 / Revised: 6 March 2021 / Accepted: 7 March 2021 / Published: 10 March 2021
(This article belongs to the Special Issue Bounce Cosmology)

Abstract

:
The main difficulties in constructing a viable early Universe bouncing model are: to bypass the observational and theoretical no-go theorem, to construct a stable non-singular bouncing phase, and perhaps, the major concern of it is to construct a stable attractor solution which can evade the Belinsky–Khalatnikov–Lifshitz (BKL) instability as well. In this article, in the homogeneous and isotropic background, we extensively study the stability analysis of the recently proposed viable non-minimal bouncing theory in the presence of an additional barotropic fluid and show that, the bouncing solution remains stable and can evade BKL instability for a wide range of the model parameter. We provide the expressions that explain the behavior of the Universe in the vicinity of the required fixed point i.e., the bouncing solution and compare our results with the minimal theory and show that ekpyrosis is the most stable solution in any scenario.

1. Introduction

In solving non-linear differential equations, initial conditions play a crucial role in determining the dynamical behavior of the system as different initial conditions may lead to different trajectories that may behave completely differently. Since gravity is highly non-linear in nature, solving early Universe evolution depends highly on the initial conditions as well. However, the already highly acceptable paradigm of the early Universe—the inflationary paradigm [1,2,3,4,5,6,7,8,9,10,11,12,13], perhaps succeeds based on the fact that, not only it is in line with the recent tighter observational constraints [14,15] but it also provides an attractor solution [16,17]. This implies that the early Universe inflationary solution does not depend on the initial conditions and because of this reason, even if the Universe contains additional matter(s), inflation makes them decay exponentially and the Universe solely dominates by the inflaton field responsible for inflationary evolution.
However, the greatest achievement comes from the fact that the inflationary solution can also evade Belinsky–Khalatnikov–Lifshitz (BKL) instability: the anisotropic energy density can grow faster than the energy density responsible for the evolution of the Universe [18] which can cause a highly unstable system. However, even with this much success, the inflationary paradigm also suffers many crucial difficulties. First, despite the ever-tightening observational constraints, there seem to exist many inflationary models that continue to remain consistent with the data [19,20,21,22], even leading to the concern whether inflation can be falsified at all [23]. Inflation being insensitive to initial conditions is still debatable. For instance, in Ref. [24], using non-perturbative simulations, authors showed that the inflationary expansion starts under very specific circumstances. In Ref. [25], authors pointed out the fact that small field potential fails to start the inflationary expansion under most initial conditions. It suffers from the trans-Planckian problem: the cosmological scales that we nowadays observe may correspond to length scales smaller than the Planck length at the onset of inflation, which in contradiction to the inherent assumption that the scales are classical even at the initial stage of inflation [26].
These issues lead to in search for alternatives to the inflationary paradigm and the most popular one is the classical non-singular bouncing scenario where, the Universe undergoes a phase of contraction until the scale factor reaches a minimum value, before it enters the expanding phase [27,28,29,30,31,32]. However, the problem with most bouncing models is that, it is extremely difficult to construct. While only known attractors are the ekpyrotic models [33], in fact, in Ref. [34], authors showed that the ekpyrotic contraction is a ‘super-smoother’, i.e., it is robust to a very wide range of initial conditions and avoid Kasner/mixmaster chaos and a similar statement could never be proven about any other primordial scenario, it fails to be in line with the observations. On the other hand, while matter bounce models may provide correct theoretical predictions consistent with the observations, the solution fails to be stable and thus, any tiny initial deviation can grow exponentially to cause an instability in the system. Another difficulty is, while constructing the contracting phase is rather easy to achieve, obtaining the non-singular bouncing phase is extremely difficult as it requires to violate the null energy condition: often called as the theoretical no-go theorem [35,36,37,38,39,40,41,42,43,44]. Most importantly, even if one may construct a model evading the instability and theoretical no-go theorem, perturbations suffer from many difficulties such as gradient instability, and these models fail to be in line with the observational constraints: a small tensor-to-scalar ratio ( r 0.002 0.06 ) and simultaneously, very small scalar non-Gaussianity parameter f NL O ( 1 ) : referred to as the observational no-go theorem [45,46] (This has also been addressed in the past in the context of the so-called pre-big bang scenario based on the string cosmology equations, as reported for instance in a recent review in Ref. [47]). Apart from these main three issues, in general, bouncing models possess another, rather weak, difficulty as well: a natural exit mechanism from the bouncing phase to enter into the conventional reheating phase.
In solving these problems, it is soon realized that one needs to go beyond the canonical theories and consider non-minimal couplings, viz. the Horndeski theories or even beyond Horndeski theories [35,36,37,38,39,40,41,42,45,46,48,49,50,51,52,53,54,54]. However, it remains an open problem until in Ref. [55], we showed that a simple non-minimal coupling can solve all those issues.
In Refs. [56,57], we have shown that, the issue of stability can easily be resolved by using a simple non-minimal coupling. This has been achieved with the help of conformal transformation and we also have shown that, it may lead to viable tensor spectra as well [58]. Recently, we extended the work to construct the first viable non-minimal bouncing model which evades all the issues and is consistent with the constraints [55]. This model is constructed in such a way that the scalar sector of the action is conformal to that of the slow-roll inflationary model. Since, under conformal transformation, perturbations remain invariant and the constraints are obtained from the correlation functions of the scalar and tensor perturbations, we achieved to bypass observational no-go theorem. The choice of the coupling function also helps us to achieve the non-singular bouncing phase without any instability and therefore, it evades the theoretical no-go theorem and leads to the conventional reheating scenario. It has been pointed out that since the work is similar to works in Refs. [56,57], it can also resolve the issue of stability, mainly the BKL instability as well. However, the model contains a free model parameter α (which needs to be greater than zero to achieve the bouncing solution) and there is no bound on it. It is not clear how the system behaves in the presence of an additional matter.
At this moment, we need to stress that, when a model is studied for perturbations and is compared with the observations, it is assumed that, either the initial conditions are chosen with extreme measure, or the model solution is stable enough so that there is no need to fine-tune for choosing the initial conditions. Otherwise, even without the external matter, the perturbations themselves can cause instability as it also possesses a small amount of energy. This is the main reason why a stable solution is always preferred (one may even argue that only stable solutions are allowed). Therefore, even our model is conformal to inflationary solutions and obeys all relevant observational constraints, one needs to verify the stability of the model, especially in the presence of an external matter. This leads to the main motivation of this article.
In this work, we analyze the stability of the model in presence of an additional barotropic matter. The aim is to find the general attractor behavior of the bouncing solution that is consistent with the observations. There is mainly one free model parameter α (for minimal theory, α = 1 , since for this value, the coupling function becomes unity) along with the equation of state of the additional barotropic fluid w m . In addition to that, since we consider conformal single field model that leads to slow-roll inflation, we also consider the potential slow-roll parameter in the minimal Einstein frame to be (nearly) constant (in case of exponential potential is strictly constant, as examined in Refs. [16,17]). With the help of these model parameters, we find the condition for stability for the bouncing solution and show that, to avoid the BKL instability leads to α < 2 . Since bouncing solution requires α > 0 , the bound becomes 0 < α < 2 . To understand the behavior of the stability, we study the phase space in the vicinity of the fixed point corresponding to the bouncing solution and find that the system acts like the Universe contains three different types of matter and due to the attractor behavior, only the leading bouncing matter component remains intact, whereas, other two decay very fast. This helps us realize how the external barotropic fluid influences the system and how the attracting nature of the bouncing solution gets rid of the initial deviations.
A few words on our conventions and notations are in order at this stage of our discussion. In this work, we work with the natural units such that = c = 1 , and we define the Planck mass to be M Pl = ( 8 π G ) 1 / 2 . We adopt the metric signature of ( , + , + , + ) . We should mention that, while the Greek indices are contracted with the metric tensor g μ ν , the Latin indices contracted with the Kronecker delta δ i j . Moreover, we shall denote the partial and the covariant derivatives as ∂ and ∇. The overdots and overprimes, as usual, denote derivatives with respect to the cosmic time t and the conformal time η associated with the Friedmann-Lemaître-Robertson-Walker (FLRW) line-element, respectively. The sub(super)script ‘I’ and b denote the quantity in the minimal Einstein theory and non-minimal theory, respectively.

2. A Brief Introduction of the Bouncing Model

The non-minimal bouncing model is constructed in such a way that the scalar sector of the non-minimal action transforms conformally to that of a minimal canonical Einstein model that leads to slow-roll inflation. The slow-roll inflationary model can easily be constructed by using a single canonical scalar field ϕ minimally coupled to the gravity, i.e., the Einstein gravity as
S I = 1 2 d 4 x g I M Pl 2 R I g I μ ν μ ϕ ν ϕ 2 V I ( ϕ ) + S M I ( g μ ν , Ψ M ) .
R I is the Ricci scalar for the metric g μ ν I , and V I ( ϕ ) is the potential responsible for the inflationary solution. The part of the action S M I ( g μ ν , Ψ M ) indicates the presence of an additional fluid. As mentioned above, since inflation is an attractor, the evolution of the field does not depend on the additional fluid and is solely governed by the scalar sector of the system. Assuming the above theory is responsible for the slow-roll inflationary dynamics, the scale factor solution, during inflation, can approximately be written as a function of the scalar field ϕ as
a I ( ϕ ) exp ϕ d ϕ M Pl 2 V I V I , ϕ ,
with V I , ϕ V I ϕ . Now we shall construct a model which is conformal to the scalar part of the above inflationary action (1) in such a way that the new scale factor behaves as a bouncing solution. The transformation can be written as
g μ ν I = f 2 ( ϕ ) g μ ν b a I ( η ) = f ( ϕ ) a b ( η ) .
g μ ν b and a b ( η ) are the required bouncing metric and scale factor solutions, respectively, along with f ( ϕ ) being the coupling function. Using the above conformal transformation along with the action (1), we can construct the action responsible for such bouncing solution as
S b = 1 2 d 4 x g b M Pl 2 f 2 ( ϕ ) R b ω ( ϕ ) g b μ ν μ ϕ ν ϕ 2 V b ( ϕ ) + S M b .
ω ( ϕ ) and the potential V b ( ϕ ) depend on the coupling function f ( ϕ ) and the inflationary potential V I ( ϕ ) as
ω ( ϕ ) = f 2 ( ϕ ) 1 6 M Pl 2 f , ϕ 2 f 2 , V b ( ϕ ) = f 4 ( ϕ ) V I ( ϕ ) ,
where f , ϕ f ϕ . S M b is the external fluid present in the system in the non-minimal theory. Please note that, the above action resembles very similar to the gravi-dilaton string effective action, with the inclusion of non-perturbative string-coupling corrections [47]. Therefore, one can explain the existence of such kind of coupling functions in the domain of string theory. Note that, while the scalar part of the above action (4) is conformal to that of (1), these two actions in reality are not conformal due to the presence of the additional fluid. However, if the scalar field dominated solution in the non-minimal theory also becomes stable, similar to minimal theory, then, for a given inflationary model (1) with its solution (2), one can obtain the bouncing scale factor solution a b ( η ) in the non-minimal theory by setting f ( ϕ ) in a certain manner. For simplicity, if we desire the bouncing solution of the form
a b ( η ) ( η ) α exp α ϕ d ϕ M Pl 2 V I V I , ϕ , α > 0 ,
where η is the comoving time, then by using Equations (2) and (3), one can obtain the required solution of the coupling function f ( ϕ ) as
f ( ϕ ) = f 0 exp α + 1 M Pl 2 ϕ d ϕ V I V I , ϕ .
f 0 is normalized in such a way that at the minima of the potential, f ( ϕ ) becomes unity. The non-minimal theory (4) is the desired bouncing model: given an inflationary potential V I ( ϕ ) , we can completely determine the coupling function f ( ϕ ) , which in turn, provides non-minimal theory (4) responsible for the bouncing scale factor solution (6). Note that, since Equation (2) is an approximated form, by using the above expressions (7) and (3), the obtained scale factor solution a b ( η ) is also not exact, but close to (6) (we will evaluate it later). We need to stress that the form of the scale factor (6) is valid only in the contracting phase, and during and after the bounce, it is extremely difficult to obtain an analytical solution. Therefore, the bouncing model, in our case, is also asymmetrical. The underlying assumption of the theory is that, either the bouncing solution is stable or we must choose a precise initial condition that leads to the bounce. Last of all, while α > 0 provides a bouncing solution, there are actually no bound on α as it can take any value and negative values of it will lead to different Universe with different scale factor solution (6). For instance, α = 1 brings us the conventional minimal Einstein scenario as, for this value, the coupling function becomes unity. In this work, we will concentrate only on the bouncing solutions where α > 0 and obtain the condition for stability and the behavior of the system in the vicinity of the contracting solution, which we will explore in the next section.
Few things to note before we move in the next section. In case of scalar perturbation, in the Einstein frame (1), the perturbed action is
d η d x 3 z I 2 ζ I 2 ζ I 2 , z I ( η ) a I ϕ H I ,
where, ζ is the curvature perturbation and ∇ is the divergence operator. In the non-minimal frame (4), the expression of the action for the scalar perturbation changes by replacing z I ( η ) with z b ( η ) . However, z b ( η ) is simply the conformally transformed z I ( η ) , i.e.,
z b ( η ) = f ( ϕ ) a b ( η ) ϕ H b + ϕ f , ϕ ( ϕ ) f ( ϕ )
This implies that, z b ( η ) = z I ( η ) which leads to the action being identical in both frames and hence, ζ I = ζ b at linear order. As the speed of sound is conformally invariant, c s is unity in both minimal inflationary as well as in non-minimal bouncing scenarios. Similarly, one can evaluate the perturbed interaction Hamiltonian (for detailed evaluation, see Refs. [59,60,61,62,63]) at any order and show that it also remains invariant in all conformal frames. This implies that the curvature perturbations at any order in both the frames are identical, which is not surprising, as we know, under conformal transformation, curvature perturbation remains invariant. The argument extends to tensor perturbation as well. Therefore, given a viable inflationary model which is consistent with the observations, in our bouncing model, there appear no divergences or instabilities (e.g., gradient instability) in the perturbations anywhere, even at the bounce and, in addition to that, the model also satisfies all observational constraints: thus evading the no-go theorem. The model leads to stable non-singular bounce as well as, shortly after the bounce, it leads to conventional reheating phase (for detailed information, see Figure 1 and Figure 2, along with Ref. [55] as well as [64], to compare reheating in minimal theory with the same in non-minimal theory).
In short, for α > 0 , the model (4) leads to a stable, asymmetrical bouncing Universe, where the scale factor in the contracting region behaves as a ( η ) ( η ) α . If the slow-roll inflationary potential V I ( ϕ ) satisfies the observational constraints (e.g., the Starobinsky potential), the corresponding bouncing model, irrespective of the value of α , also satisfied all observational constraints, as the perturbations remain invariant under conformal transformation, thus evading the observational no-go theorem.

3. Stability Analysis

Let us now analyze the stability condition for the above system (4). Assuming the additional fluid is barotropic, in nature, with the equation of state w m , the governing equations become
f 2 R μ ν 1 2 g μ ν R 2 μ f ν f 2 f μ ν f + 2 g μ ν λ f λ f + f f = 1 M Pl 2 ω ( ϕ ) μ ϕ ν ϕ 1 2 g μ ν λ ϕ λ ϕ g μ ν V b + T μ ν ( M ) ,
ϕ + 1 2 ω ( ϕ ) ω , ϕ λ ϕ λ ϕ 2 V , ϕ + 2 M Pl 2 f f , ϕ R = 0 ,
where, T μ ν ( M ) is the energy-momentum tensor of the additional barotropic fluid in the non-minimal theory (4), satisfying
μ T μ ν ( M ) = 0 .
We have not used any subscript b since, from now on, all concerning equations will correspond only to the non-minimal theory. These equations, using Friedmann–Robertson–Walker (FRW) homogeneous and isotropic Universe with the line element
d s 2 = d t 2 + a 2 ( t ) d x 2 ,
turn into
3 f 2 H 2 = 6 H f f , ϕ + 1 M Pl 2 1 2 ω ϕ ˙ 2 + V ( ϕ ) + ρ M = 0 ,
ϕ ¨ + 3 H ϕ ˙ + 1 2 ω ω , ϕ ϕ ˙ 2 + V , ϕ 2 M Pl 2 f f , ϕ R = 0 ,
ρ ˙ M + 3 H ρ M + P M = 0 ,
where ρ M T ( M ) 0 0 and T j i P M δ j i = w m ρ M δ j i are the energy density and the pressure of the barotropic fluid with w m being the equation of state. To study the stability analysis, instead of using quantities with dimensions, we can express and re-write these field equations in terms of dimensionless quantities. These are defined as
x ϕ ˙ 6 M Pl H , y V 3 M Pl f H
Using these dimensionless variables, one can quickly express the evolution equation of them as
d x d N 1 H d x d t = 1 2 γ e 4 ( 36 ( 3 w m 1 ) ( 1 + α ) 4 x 3 18 6 ( 3 w m 1 ) ( 1 + α ) 3 γ e x 2 6 ( 1 + α ) 2 γ e 2 ( 3 4 x 2 y 2 3 w m ( 3 2 x 2 y 2 ) ) x + 6 ( 1 7 x 2 y 2 3 w m ( 1 3 x 2 y 2 ) ) ( 1 + α ) γ e 3 3 x ( ( w m 1 ) ( 1 x 2 y 2 ) + 2 y 2 α ) γ e 4 + 6 y 2 γ e 5 ) ,
d y d N 1 H d y d t = 1 2 γ e 4 ( ( 6 ( 1 3 w m ) ( y 2 1 ) ( 1 + α ) 2 γ e 2 y + 3 ( 1 + w m + y 2 w m y 2 + 2 y 2 α ) γ e 4 + 6 γ e ( 12 ( 3 w m 1 ) ( 1 + α ) 3 x 6 w m ( 1 + α ) γ e 2 + γ e 4 ) + 3 ( 12 ( 1 + α ) 4 8 ( 1 + α ) 2 γ e 2 + γ e 4 w m ( 6 ( 1 + α ) 2 γ e 2 ) 2 ) x 2 ) ) ,
along with the constraint equation
Ω M ρ M 3 M Pl 2 f 2 H 2 = 1 2 6 ( 1 + α ) γ e x + 6 ( 1 + α ) 2 γ e 2 1 x 2 y 2 , γ e M Pl V I , ϕ V I .
N is defined as the e-folding number in the non-minimal frame ln a / a 0 , a 0 is the initial value of the scale factor and V I ( ϕ ) is the potential in the minimal Einstein frame and is related to V ( ϕ ) through Equation (5). Note that, the model parameter γ e actually represents a parameter in the Einstein frame which is completely arbitrary and can be fixed from the observations. For example, while near the pivot scale, for chaotic inflation, γ e 0.01 , for Starobinsky model, at the same scale, it becomes 0.003 . However, chaotic inflation has already been ruled out from the observations while, the Starobinsky model remains consistent with it and one of the most important models in the paradigm. We should stress that the phase during the stability is analyzed is the contraction in a bouncing model and therefore, instead of the e-N-fold time convention, here we are using the e-fold time convention.
One can also express the slow-roll parameter ϵ b in terms of these dimensionless parameters as
ϵ b H ˙ H 2 = 1 2 γ e 4 ( 6 ( 1 3 w m ) ( 1 y 2 ) ( 1 + α ) 2 γ e 2 + 3 ( 1 + w m + y 2 w m y 2 + 2 y 2 α ) γ e 4 + 2 6 γ e ( 3 w m 1 ) ( 1 + α ) ( 6 ( 1 + α ) 2 γ e 2 ) x + 3 ( 12 ( 1 + α ) 4 8 ( 1 + α ) 2 γ e 2 + γ e 4 w m ( 6 ( 1 + α ) 2 γ e 2 ) 2 ) x 2 ) .
Using the above expression, one can define the effective equation of state as
w eff = 1 + 2 3 ϵ b .
Although γ e is a function of ϕ and thus, it is dynamical, the minimal theory that we are considering leads to slow-roll inflation and thus one can consider it to be nearly constant, similar to exponential potential.
In this system, as one can see, there exist three model parameters: { α , w m , γ e } . w m is arbitrary as it describes the nature of the additional field. For example, in the case of anisotropic stress fluid that corresponds to the stiff matter, w m = 1 . This field results in the BKL instability that, later, we will investigate. As mentioned above, γ e is nearly constant and hence, we are left with one model parameter α that can be constrained. In the next section, we will try to constrain it using the Lyapunov exponents.

3.1. Fixed Points and the Bouncing Solution

Using the evolution equation, one can find the fixed/critical point of the system and these points correspond to the solutions of the system, in this case, the Universe depicted by the non-minimal theory (4). These can be found by setting Equations (18) and (19) to be equal to zero, i.e., the velocities of x and y vanishes at these points. There are seven such fixed points:
1 . x 1 * = γ e 6 α , y 1 * = 6 γ e 2 6 α
2 . x 2 * = γ e 6 α , y 2 * = 6 γ e 2 6 α
The sign of the y * signifies whether the Universe is expanding or contracting. Since, H appears in the denominator in the expression of y in (17), the negative sign of y leads to contraction (bouncing), while positive y corresponds to expanding Universe. Therefore, { x 1 * , y 1 * } leads to contraction, and { x 2 * , y 2 * } implies the expanding Universe. In fact, { x 1 * , y 1 * } is our desired solution and the scalar field dominated solution as the fractional energy density of the additional fluid and the effective equation of state (see Equation (22)) of corresponds to this point are
Ω M ( 1 , 2 ) = 0
w eff ( 1 , 2 ) = 2 α γ e 2 3 α
which corresponds to the the scale factor solution as
a b ( η ) ( η ) β , β 2 α 2 γ e 2 .
As you can see, the exponent is not exactly α but β , which is contrary to the Equation (6). However, as mentioned before, γ e 1 for slow-roll models and only by using such model, we achieved such non-minimal bounce, and therefore, β α . Therefore, { x 1 * , y 1 * } corresponds to our bouncing solution that we will investigate thoroughly in the next section. Other fixed points are
3 . x 3 * = 6 γ e 6 ( 1 + α ) 6 γ e , y 3 * = 0
4 . x 4 * = 6 γ e 6 ( 1 + α ) 6 γ e , y 4 * = 0
5 . x 5 * = 3 2 ( 1 + w m ) γ e 4 ( 1 + α ) + γ e 2 , y 5 * = 2 ( 1 3 w m ) 2 α ( 1 + α ) 2 + α ( 1 3 w m 2 2 α + 6 w m ( 1 + α ) ) γ e 2 ( 2 α ( 4 ( 1 + α ) γ e 2 ) )
6 . x 6 * = 3 2 ( 1 + w m ) γ e 4 ( 1 + α ) + γ e 2 , y 6 * = 2 ( 1 3 w m ) 2 α ( 1 + α ) 2 + α ( 1 3 w m 2 2 α + 6 w m ( 1 + α ) ) γ e 2 ( 2 α ( 4 ( 1 + α ) γ e 2 ) )
7 . x 7 * = 2 3 ( 3 w m 1 ) ( 1 + α ) γ e 2 ( 1 + 3 w m ) ( 1 + α ) 2 ( 1 + w m ) γ e 2 , y 7 * = 0
Since our desired fixed point is { x 1 * , y 1 * } , in this work, we will focus only on this point. In the next subsection, we will explore the Lyapunov exponents that determine the stability of the solution, and then later, we will focus on the behavior of the Universe close to the fixed point.

3.2. Lyapunov Exponent and the Stability of the Bouncing Solution

Now, in order to study the stability of these fixed points, we need to linearize the Equations (18) and (19) as
d δ x d N d δ y d N = A ( x , y ) x | * A ( x , y ) y | * B ( x , y ) x | * B ( x , y ) y | * δ x δ y ,
where, A ( x , y ) and B ( x , y ) are the right-hand side of (18) and (19), respectively. | * denotes the value at the fixed point. δ x and δ y are the deviations of the fixed points. By linearizing equations, we assume that we are studying the stability condition in the vicinity of the fixed points, i.e., the deviation from the fixed points are very small. We are considering homogeneous and isotropic background and assume that the anisotropies and inhomogeneities caused by the deviations do not impact on the background. This can be achieved by assuming that the energy fractions due to the deviations are extremely small.
In order to check the stability of the fixed points, we need to evaluate the eigenvalues of this matrix. These eigenvalues are the Lyapunov exponents that arise in the evolution of both δ x and δ y and thus, only these values determine the stability of the fixed points as
δ x = C 11 e λ 1 N + C 12 e λ 2 N δ y = C 21 e λ 1 N + C 22 e λ 2 N
As defined earlier, N is the e-folding number; in an expanding Universe, while it increases, it decreases for contracting solution, and therefore, if both the eigenvalues are negative (positive) in an expanding (contracting) Universe, then δ x and δ y approach zero as N approaches ( ) . This implies that the deviation from the actual trajectory reduces with time and the new trajectory returns to the original trajectory asymptotically in time. In other words, both eigenvalues being negative (positive) in an expanding (contracting) Universe implies that the fixed point is stable and the corresponding solution is an attractor solution.
Let us now calculate the Lyapunov exponents for our desired fixed point (23). It takes the form
λ 1 ( 1 ) = 4 + α 3 w m α γ e 2 α , λ 2 ( 1 ) = 6 γ e 2 2 α
The above expression bears one of the main results of this work. As one can see properly since the desired fixed point is the bouncing solution, we require both of the λ ’s to be positive. Few things to note: since α = 1 represents the minimal solution as mentioned before (coupling function (7) becomes unity), one can easily obtain the results in Refs. [16,17] by substituting α = 1 . Since α , w m and γ e are positive, increasing values of w m and γ e narrow the permissible regime of α for being an attractor. Anisotropic fluid or the stress fluid resembles similar to a stiff matter with the equation of state w m = 1 and the energy density ρ M a 6 , where a is the scale factor. Therefore, since BKL instability is caused by it, the condition to evade the BKL instability can easily be obtained by substituting w m = 1 :
γ e 2 < 6 , α < 2 γ e 2 2 .
Therefore, for exponential potential, γ e is constant and positive domain of α can only be obtained if γ e < 2 . In the case of slow-roll, γ e 1 and therefore, for slow-roll inflationary models which are conformal to the non-minimal bounce, the condition becomes
0 < α < 2 .
as the bouncing solution can only be obtained for α > 0 . This is the main result of this article as, now, α cannot possess arbitrary positive values for a stable viable solution. In the next section, in order to understand how the deviations impact the system near the fixed point, we will study the phase space in the concerned region.

4. Evolution of the Deviations

Using the eigenvalues and the eigenvectors, one can easily solve (33). For the fixed point (23), it becomes
δ x ( 1 ) ( N ) = A 11 ( α , w m , γ e ) e λ 1 ( 1 ) N + A 12 ( α , w m , γ e ) e λ 2 ( 1 ) N
δ y ( 1 ) ( N ) = A 21 ( α , w m , γ e ) e λ 1 ( 1 ) N + A 22 ( α , w m , γ e ) e λ 2 ( 1 ) N ,
where, λ 1 ( 1 ) and λ 2 ( 1 ) are given in (34). The A’s can be evaluated as
A 11 = ( 2 ( 3 w m 1 ) ( 1 + α ) ( w m 1 ) γ e 2 ) ( ( 6 + 6 α γ e 2 ) δ x 0 + γ e 6 γ e 2 δ y 0 ) α γ e 2 ( γ e 2 2 2 ( 1 3 w m ) α )
A 12 = 1 α γ e 2 ( γ e 2 2 2 ( 1 3 w m ) α ) ( γ e 2 6 ) ( 2 ( 3 w m 1 ) ( 1 + α ) 2 + ( 1 w m + α ) γ e 2 ) δ x 0 + γ e 6 γ e 2 ( 2 ( 1 3 w m ) ( 1 + α ) + ( w m 1 ) γ e 2 ) δ y 0
A 21 = 1 α γ e 3 ( γ e 2 2 2 ( 1 3 w m ) α ) ( 2 ( 3 w m 1 ) ( 1 + α ) 2 + ( 1 w m + α ) γ e 2 ) ( 6 γ e 2 ( γ e 2 6 6 α ) δ x 0 + γ e ( γ e 2 6 ) δ y 0 )
A 22 = 1 α γ e 3 ( γ e 2 2 2 ( 1 3 w m ) α ) ( 6 + 6 α γ e 2 ) ( 6 γ e 2 ( 2 ( 3 w m 1 ) ( 1 + α ) 2 + ( 1 w m + α ) γ e 2 ) δ x 0 + γ e ( 2 ( 3 w m 1 ) ( 1 + α ) ( w m 1 ) γ e 2 ) δ y 0 )
These solutions are obtained by using the initial conditions δ x ( 0 ) = δ x 0 , δ y ( 0 ) = δ y 0 . The slow-roll parameter can also be obtained in the vicinity of the fixed point (23) by perturbing the Equation (21) and this becomes
ϵ b = ϵ b ( 1 ) + ϵ b x δ x ( N ) + ϵ b y δ y ( N ) , ϵ b ( 1 ) = 1 + 1 α γ e 2 2 α ,
where,
ϵ b x ( 1 ) = 3 2 ( 6 + 12 α + 6 α 2 γ e 2 ) ( 2 + 2 α γ e 2 w m ( 6 + 6 α γ e 2 ) ) α γ e 3
ϵ b y ( 1 ) = 9 3 γ e 2 2 ( 2 + 2 α 2 γ e 2 + 2 α ( 2 γ e 2 ) w m ( 6 + 12 α + 6 α 2 γ e 2 ) ) α γ e 2 .
Once we obtain the solution expression for the slow-roll parameter, we can solve the energy density equation as ϵ ( N ) = d N ln H ( N ) and by performing the integral, one can easily solve for H. It becomes, H 2 ( N ) = H 0 2 exp 2 d N ϵ ( N ) , H 0 is the initial value of the Hubble parameter. Using (37) and (43), and by using the fact that, δ x 0 and δ y 0 are tiny and act like perturbations, we obtain the evolution of the Hubble parameter as
H 2 H 0 2 = ( 1 Ω m 1 Ω m 2 ) a a 0 2 ϵ b ( 1 ) + Ω m 1 a a 0 λ 1 ( 1 ) 2 ϵ b ( 1 ) + Ω m 2 a a 0 λ 2 ( 1 ) 2 ϵ b ( 1 )
where,
Ω m 1 = 2 ϵ b x ( 1 ) A 11 λ 1 ( 1 ) 2 ϵ b y ( 1 ) A 21 λ 1 ( 1 )
Ω m 2 = 2 ϵ b x ( 1 ) A 12 λ 2 ( 1 ) 2 ϵ b y ( 1 ) A 22 λ 2 ( 1 )
One can easily obtain the the above expressions in terms of model parameters { α , w m , γ e } by using relations (39)–(42), (44) and (45). It is clear from the above equation that the system behaves as the Universe contains three fluids, where the leading solution corresponds the scalar field dominated solution, in the case of additional fluids having the initial energy density fraction Ω m 1 and Ω m 2 , the energy density of them grow/decay as ρ m 1 a a 0 λ 1 ( 1 ) 2 ϵ b ( 1 ) and ρ m 2 a a 0 λ 2 ( 1 ) 2 ϵ b ( 1 ) . It is apparent that, relative to the leading order term with a a 0 2 ϵ b ( 1 ) , the other two decay/grow with the factor a a 0 λ 1 ( 1 ) and a a 0 λ 2 ( 1 ) , respectively. Therefore, for a contracting solution, λ ’s being positive signifies that the those additional fluids represented by the Ω m 1 and Ω m 2 decay. This can even be seen from the expression of the energy fraction of the additional fluid (20). This becomes
Ω M = 2 3 ( ( 6 + 6 α γ e 2 ) δ x 0 + γ e 6 γ e 2 δ y 0 ) ( α γ e ) a a 0 λ 1 N .
From the above expression, it becomes obvious that, for a bouncing solution, if λ ’s are positive, the energy fraction vanishes after a while. To understand the correct nature of it, below, we presented three different scenarios with different choices of model parameters that try to help us realize the system.

4.1. Radiation Bounce

As we know from (36), α must be positive and less than 2 in order to satisfy bouncing solution as well as to evade the BKL instability. Therefore, in the first example, we consider radiation bounce where α = 1 , i.e., the scale factor solution approximately becomes a b ( η ) ( η ) . We consider the additional fluid as the anisotropic stress fluid with stiff equation of state w m = 1 . Assuming typical conformal slow-roll model (e.g., chaotic inflation with V I ( ϕ ) = 1 2 m 2 ϕ 2 ) with γ e 0.01 , the Hubble solution takes the form
H 2 H 0 2 = ( 1 Ω m 1 Ω m 2 ) a a 0 4 + Ω m 1 a a 0 2 + Ω m 1 a a 0 1
with
Ω m 1 Ω m 2 = 1.125 δ x 0 + 0.002 δ y 0 δ x 0 + 0.002 δ y 0 .
The equation represents the attractor nature of the solution as the leading order solution, i.e., the effective radiation matter solution dominates over the additional fluid evolution. The actual energy density fraction of the additional fluid (49) takes the form
Ω M = 980 δ x 0 + 2 δ y 0 a a 0 2 .
For bouncing solution, as it is obvious, it quickly vanishes after a while.

4.2. Comparing with the Minimal Scenario

As it has been mentioned before, negative α represents the expanding solution, and α = 1 reduces to the minimal Einstein inflationary scenario as for this value, the coupling function becomes unity. In this subsection, we will compare the above result with the minimal inflationary case. Again using the similar value of the additional fluid equation of state as well the γ e , i.e., w m = 1 , γ e 0.01 , the Hubble solution becomes
H 2 H 0 2 = ( 1 Ω m 1 Ω m 2 ) + Ω m 1 a a 0 6 + Ω m 1 a a 0 3 ,
with
Ω m 1 Ω m 2 = 0.05 + 122.47 δ y 0 δ x 0 .
Notice that, in this case, the additional fluids with Ω m 1 and Ω m 2 initial mass fractions decay faster compared to the radiation bounce scenario as the Lyapunov exponents, i.e., the λ ’s are having high values. This can even be seen from the expression of energy fraction of the additional fluid (49):
Ω M = 0.008 δ x 0 2 δ y 0 a a 0 6
It vanishes faster than that of the radiation bounce case. Therefore, one may consider, by judging the capability of stability solution, minimal inflationary models are preferred compared to the non-minimal bouncing models. In the next subsection, we will show that, this is clearly not the case.

4.3. Comparing with the Viable Ekpyrosis

In Ref. [34], it has been studied extensively and showed that the ekpyrosis is a super-attractor (although it has been studied for minimal models, one can easily extend it for non-minimal solution as well. See Ref. [56] in this context). It can also been seen from (34) as α appears in the denominator in the expressions. This implies that the smaller values of α leads to higher values of Lyapunov exponents and this, in return, implies the solutions are, in comparison, more stable. Consider the example of α = 0.1 : the solution leads to the ekpyrosis solution. By using similar values of w m and γ e , one can obtain the Hubble solution for the ekpyrosis as
H 2 H 0 2 = ( 1 Ω m 1 Ω m 2 ) a a 0 22 + Ω m 1 a a 0 16 + Ω m 1 a a 0 8
with
Ω m 1 Ω m 2 = 0.95 δ x 0 0.003 δ y 0 δ x 0 + 0.003 δ y 0 .
Again, compared to the leading solution, the externals fluids decay much faster than the minimal case, shown above. This is because the Lyapunov exponents are, in this case, 38 and 30, respectively. For instance, the energy fraction of the additional fluid becomes
Ω M = 5389 δ x 0 + 20 δ y 0 a a 0 38 .
Therefore, even if all the examples, that are presented in this work, are in line with the observations, there is a difference in the behavior of the stability and among them, ekpyrosis is the most stable, and in this sense, the ekpyrosis is the most preferred, which again, has been re-established.

5. Conclusions

In this article, we construct a non-minimal bouncing model that is conformal to the scalar sector of the minimal inflationary model. However, along with the presence of the additional matter in the system, the minimal and non-minimal models do not conform and therefore, the stability of the two systems is not related by the conformal factor. α is the main model parameter and by setting it to be equal to 1 , one can arrive at the stability conditions of the minimal model as well as it had been studied before in the literature. While, the scalar part of the non-minimal bounce is conformal to that of the minimal inflationary model, choosing a viable inflationary model that satisfies the observational constraints, the non-minimal bouncing model also satisfies the observation constraints with the inherent assumption that the model solution is an attractor. If the solution is not stable, then even the quantum fluctuations can grow and diverge and the system becomes unstable. In this article, therefore, we study the stability of the non-minimal bouncing model with the presence of an additional fluid with the equation of state w m . We show that, while α being positive leads to a viable bouncing solution, not all positive values of it leads to a stable solution, and to avoid BKL instability, we obtained the constraint in α as the domain for stability as well as bounce as 0 < α < 2 . We also showed that the ekpyrosis is the best attractor among them and thus is always preferred.
While we studied the system for homogeneous and isotropic background, one can easily extend the work for homogeneous and anisotropic systems as well and may arrive at a similar condition. One can even go beyond the homogeneous system and study the stability using numerical relativity. However, this is beyond the scope of this article. We did not consider the negative values of α as well as α = 0 case. The last one signifies the emergent gravity scenario. Apart from that, we also did not study the behavior of the other fixed points (28)–(32). While in the minimal scenario, most of them are unstable, in our case, it may not be the same and may impact the system heavily. We reserve these works for our future prospects.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The author thanks the Indian Institute of Technology Madras (IIT Madras), Chennai, India for support through the institute postdoctoral fellowship.

Conflicts of Interest

The author declares no conflict of interest.

References

  1. Starobinsky, A. A new type of isotropic cosmological models without singularity. Phys. Lett. B 1980, 91, 99–102. [Google Scholar] [CrossRef]
  2. Sato, K. First-order phase transition of a vacuum and the expansion of the Universe. Mon. Not. R. Astron. Soc. 1981, 195, 467–479. [Google Scholar] [CrossRef]
  3. Guth, A.H. Inflationary universe: A possible solution to the horizon and flatness problems. Phys. Rev. D 1981, 23, 347–356. [Google Scholar] [CrossRef] [Green Version]
  4. Linde, A. A new inflationary universe scenario: A possible solution of the horizon, flatness, homogeneity, isotropy and primordial monopole problems. Phys. Lett. B 1982, 108, 389–393. [Google Scholar] [CrossRef]
  5. Albrecht, A.; Steinhardt, P.J. Cosmology for Grand Unified Theories with Radiatively Induced Symmetry Breaking. Phys. Rev. Lett. 1982, 48, 1220–1223. [Google Scholar] [CrossRef]
  6. Linde, A.D. Chaotic Inflation. Phys. Lett. 1983, 129B, 177–181. [Google Scholar] [CrossRef]
  7. Mukhanov, V.F.; Chibisov, G.V. Quantum Fluctuations and a Nonsingular Universe. JETP Lett. 1981, 33, 532–535. [Google Scholar]
  8. Hawking, S. The development of irregularities in a single bubble inflationary universe. Phys. Lett. B 1982, 115, 295–297. [Google Scholar] [CrossRef]
  9. Starobinsky, A. Dynamics of phase transition in the new inflationary universe scenario and generation of perturbations. Phys. Lett. B 1982, 117, 175–178. [Google Scholar] [CrossRef]
  10. Guth, A.H.; Pi, S.Y. Fluctuations in the New Inflationary Universe. Phys. Rev. Lett. 1982, 49, 1110–1113. [Google Scholar] [CrossRef]
  11. Vilenkin, A. Quantum fluctuations in the new inflationary universe. Nucl. Phys. B 1983, 226, 527–546. [Google Scholar] [CrossRef]
  12. Bardeen, J.M.; Steinhardt, P.J.; Turner, M.S. Spontaneous creation of almost scale-free density perturbations in an inflationary universe. Phys. Rev. D 1983, 28, 679–693. [Google Scholar] [CrossRef]
  13. Starobinsky, A.A. Spectrum of relict gravitational radiation and the early state of the universe. JETP Lett. 1979, 30, 682–685. [Google Scholar]
  14. Akrami, Y.; Arroja, F.; Ashdown, M.; Aumont, J.; Baccigalupi, C.; Ballardini, M.; Banday, A.J.; Barreiro, R.B.; Bartolo, N.; Basak, S.; et al. Planck 2018 results. X. Constraints on inflation. arXiv 2018, arXiv:1807.06211. [Google Scholar]
  15. Aghanim, N.; Akrami, Y.; Ashdown, M.; Aumont, J.; Baccigalupi, C.; Ballardini, M.; Banday, A.J.; Barreiro, R.B.; Bartolo, N.; Basak, S.; et al. Planck 2018 results. VI. Cosmological parameters. arXiv 2018, arXiv:1807.06209. [Google Scholar]
  16. Copeland, E.J.; Liddle, A.R.; Wands, D. Exponential potentials and cosmological scaling solutions. Phys. Rev. D 1998, 57, 4686–4690. [Google Scholar] [CrossRef] [Green Version]
  17. Ng, S.; Nunes, N.; Rosati, F. Applications of scalar attractor solutions to cosmology. Phys. Rev. D 2001, 64, 083510. [Google Scholar] [CrossRef] [Green Version]
  18. Belinskii, V.; Khalatnikov, I.; Lifshitz, E. Oscillatory approach to a singular point in the relativistic cosmology. Adv. Phys. 1970, 19, 525–573. [Google Scholar] [CrossRef]
  19. Martin, J.; Ringeval, C.; Trotta, R. Hunting Down the Best Model of Inflation with Bayesian Evidence. Phys. Rev. 2011, D83, 063524. [Google Scholar] [CrossRef] [Green Version]
  20. Martin, J.; Ringeval, C.; Vennin, V. Encyclopædia Inflationaris. Phys. Dark Univ. 2014, 5–6, 75–235. [Google Scholar] [CrossRef] [Green Version]
  21. Martin, J.; Ringeval, C.; Trotta, R.; Vennin, V. The Best Inflationary Models After Planck. J. Cosmol. Astropart. Phys. 2014, 1403, 039. [Google Scholar] [CrossRef]
  22. Martin, J.; Ringeval, C.; Vennin, V. How Well Can Future CMB Missions Constrain Cosmic Inflation? J. Cosmol. Astropart. Phys. 2014, 1410, 038. [Google Scholar] [CrossRef] [Green Version]
  23. Gubitosi, G.; Lagos, M.; Magueijo, J.; Allison, R. Bayesian evidence and predictivity of the inflationary paradigm. J. Cosmol. Astropart. Phys. 2016, 1606, 002. [Google Scholar] [CrossRef]
  24. East, W.E.; Kleban, M.; Linde, A.; Senatore, L. Beginning inflation in an inhomogeneous universe. J. Cosmol. Astropart. Phys. 2016, 1609, 010. [Google Scholar] [CrossRef] [Green Version]
  25. Clough, K.; Lim, E.A.; DiNunno, B.S.; Fischler, W.; Flauger, R.; Paban, S. Robustness of Inflation to Inhomogeneous Initial Conditions. J. Cosmol. Astropart. Phys. 2017, 1709, 025. [Google Scholar] [CrossRef] [Green Version]
  26. Martin, J.; Brandenberger, R.H. The TransPlanckian problem of inflationary cosmology. Phys. Rev. D 2001, 63, 123501. [Google Scholar] [CrossRef] [Green Version]
  27. Novello, M.; Bergliaffa, S.E.P. Bouncing Cosmologies. Phys. Rept. 2008, 463, 127–213. [Google Scholar] [CrossRef]
  28. Cai, Y.F. Exploring Bouncing Cosmologies with Cosmological Surveys. Sci. China Phys. Mech. Astron. 2014, 57, 1414–1430. [Google Scholar] [CrossRef] [Green Version]
  29. Battefeld, D.; Peter, P. A Critical Review of Classical Bouncing Cosmologies. Phys. Rept. 2015, 571, 1–66. [Google Scholar] [CrossRef] [Green Version]
  30. Lilley, M.; Peter, P. Bouncing alternatives to inflation. Comptes Rendus Phys. 2015, 16, 1038–1047. [Google Scholar] [CrossRef]
  31. Ijjas, A.; Steinhardt, P.J. Implications of Planck2015 for inflationary, ekpyrotic and anamorphic bouncing cosmologies. Class. Quant. Grav. 2016, 33, 044001. [Google Scholar] [CrossRef]
  32. Brandenberger, R.; Peter, P. Bouncing Cosmologies: Progress and Problems. Found. Phys. 2017, 47, 797–850. [Google Scholar] [CrossRef] [Green Version]
  33. Levy, A.M.; Ijjas, A.; Steinhardt, P.J. Scale-invariant perturbations in ekpyrotic cosmologies without fine-tuning of initial conditions. Phys. Rev. 2015, D92, 063524. [Google Scholar] [CrossRef] [Green Version]
  34. Garfinkle, D.; Lim, W.C.; Pretorius, F.; Steinhardt, P.J. Evolution to a smooth universe in an ekpyrotic contracting phase with w > 1. Phys. Rev. 2008, D78, 083537. [Google Scholar] [CrossRef] [Green Version]
  35. Kobayashi, T. Generic instabilities of nonsingular cosmologies in Horndeski theory: A no-go theorem. Phys. Rev. 2016, D94, 043511. [Google Scholar] [CrossRef] [Green Version]
  36. Libanov, M.; Mironov, S.; Rubakov, V. Generalized Galileons: Instabilities of bouncing and Genesis cosmologies and modified Genesis. J. Cosmol. Astropart. Phys. 2016, 1608, 037. [Google Scholar] [CrossRef]
  37. Ijjas, A.; Steinhardt, P.J. Fully stable cosmological solutions with a non-singular classical bounce. Phys. Lett. 2017, B764, 289–294. [Google Scholar] [CrossRef] [Green Version]
  38. Banerjee, S.; Cai, Y.F.; Saridakis, E.N. Evading the theoretical no-go theorem for nonsingular bounces in Horndeski/Galileon cosmology. Class. Quant. Grav. 2019, 36, 135009. [Google Scholar] [CrossRef] [Green Version]
  39. Cai, Y.; Wan, Y.; Li, H.G.; Qiu, T.; Piao, Y.S. The Effective Field Theory of nonsingular cosmology. Eur. Phys. J. C 2017, 01, 090. [Google Scholar] [CrossRef] [Green Version]
  40. Cai, Y.; Piao, Y.S. A covariant Lagrangian for stable nonsingular bounce. J. High Energy Phys. 2017, 09, 027. [Google Scholar] [CrossRef]
  41. Kolevatov, R.; Mironov, S.; Sukhov, N.; Volkova, V. Cosmological bounce and Genesis beyond Horndeski. J. Cosmol. Astropart. Phys. 2017, 1708, 038. [Google Scholar] [CrossRef]
  42. Mironov, S.; Rubakov, V.; Volkova, V. Bounce beyond Horndeski with GR asymptotics and γ-crossing. J. Cosmol. Astropart. Phys. 2018, 1810, 050. [Google Scholar] [CrossRef] [Green Version]
  43. Easson, D.A.; Sawicki, I.; Vikman, A. G-Bounce. J. Cosmol. Astropart. Phys. 2011, 11, 021. [Google Scholar] [CrossRef]
  44. Sawicki, I.; Vikman, A. Hidden Negative Energies in Strongly Accelerated Universes. Phys. Rev. D 2013, 87, 067301. [Google Scholar] [CrossRef] [Green Version]
  45. Quintin, J.; Sherkatghanad, Z.; Cai, Y.F.; Brandenberger, R.H. Evolution of cosmological perturbations and the production of non-Gaussianities through a nonsingular bounce: Indications for a no-go theorem in single field matter bounce cosmologies. Phys. Rev. 2015, D92, 063532. [Google Scholar] [CrossRef] [Green Version]
  46. Li, Y.B.; Quintin, J.; Wang, D.G.; Cai, Y.F. Matter bounce cosmology with a generalized single field: Non-Gaussianity and an extended no-go theorem. J. Cosmol. Astropart. Phys. 2017, 1703, 031. [Google Scholar] [CrossRef] [Green Version]
  47. Gasperini, M.; Veneziano, G. String Theory and Pre-big bang Cosmology. Nuovo Cim. C 2016, 38, 160. [Google Scholar] [CrossRef]
  48. Horndeski, G.W. Second-order scalar-tensor field equations in a four-dimensional space. Int. J. Theor. Phys. 1974, 10, 363–384. [Google Scholar] [CrossRef]
  49. Gleyzes, J.; Langlois, D.; Piazza, F.; Vernizzi, F. Healthy theories beyond Horndeski. Phys. Rev. Lett. 2015, 114, 211101. [Google Scholar] [CrossRef] [PubMed]
  50. Kobayashi, T. Horndeski theory and beyond: A review. arXiv 2019, arXiv:1901.07183. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  51. Cai, Y.F.; Brandenberger, R.; Peter, P. Anisotropy in a Nonsingular Bounce. Class. Quant. Grav. 2013, 30, 075019. [Google Scholar] [CrossRef] [Green Version]
  52. Cai, Y.F.; Easson, D.A.; Brandenberger, R. Towards a Nonsingular Bouncing Cosmology. J. Cosmol. Astropart. Phys. 2012, 08, 020. [Google Scholar] [CrossRef]
  53. Ilyas, A.; Zhu, M.; Zheng, Y.; Cai, Y.F.; Saridakis, E.N. DHOST Bounce. arXiv 2020, arXiv:2002.08269. [Google Scholar] [CrossRef]
  54. Dobre, D.A.; Frolov, A.V.; Ghersi, J.T.G.; Ramazanov, S.; Vikman, A. Unbraiding the Bounce: Superluminality around the Corner. J. Cosmol. Astropart. Phys. 2018, 03, 020. [Google Scholar] [CrossRef] [Green Version]
  55. Nandi, D. Bounce from Inflation. Phys. Lett. B 2020, 809, 135695. [Google Scholar] [CrossRef]
  56. Nandi, D. Stable contraction in Brans-Dicke cosmology. J. Cosmol. Astropart. Phys. 2018, 1905, 040. [Google Scholar] [CrossRef] [Green Version]
  57. Nandi, D. Note on stability in conformally connected frames. Phys. Rev. 2019, D99, 103532. [Google Scholar] [CrossRef] [Green Version]
  58. Nandi, D.; Sriramkumar, L. Can a nonminimal coupling restore the consistency condition in bouncing universes? Phys. Rev. 2020, D101, 043506. [Google Scholar] [CrossRef] [Green Version]
  59. Chen, X.; Huang, M.x.; Kachru, S.; Shiu, G. Observational signatures and non-Gaussianities of general single field inflation. J. Cosmol. Astropart. Phys. 2007, 0701, 002. [Google Scholar] [CrossRef] [Green Version]
  60. Maldacena, J.M. Non-Gaussian features of primordial fluctuations in single field inflationary models. JHEP 2003, 0305, 013. [Google Scholar] [CrossRef]
  61. Nandi, D.; Shankaranarayanan, S. Complete Hamiltonian analysis of cosmological perturbations at all orders. J. Cosmol. Astropart. Phys. 2016, 1606, 038. [Google Scholar] [CrossRef] [Green Version]
  62. Nandi, D.; Shankaranarayanan, S. Complete Hamiltonian analysis of cosmological perturbations at all orders II: Non-canonical scalar field. J. Cosmol. Astropart. Phys. 2016, 1610, 008. [Google Scholar] [CrossRef] [Green Version]
  63. Nandi, D. Hamiltonian formalism of cosmological perturbations and higher derivative theories. arXiv 2017, arXiv:1707.02976. [Google Scholar]
  64. Nandi, D.; Saha, P. Einstein or Jordan: Seeking answers from the reheating constraints. arXiv 2019, arXiv:1907.10295. [Google Scholar]
Figure 1. The slow-roll parameter in the contracting region (top-left), the slow-roll parameter in the bouncing region (top-right), the coupling function (bottom-left), and ω ( ϕ ) (bottom-right) in the non-minimal theory (4) are plotted for the chaotic bouncing model with α = 1 with respect to the e-N-fold time convention N ( N 2 2 ln a / a 0 for contraction, N is negative, whereas, at bounce, it becomes zero and then during expansion, it becomes positive). It is apparent that, during the contraction period, the scale factor behaves as a ( η ) ( η ) , and thus, ϵ b 2 . One also finds that the bouncing solution is asymmetrical. The coupling function f ( ϕ ) as well as ω ( ϕ ) become nearly unity in the reheating epoch, which arises shortly after the bounce.
Figure 1. The slow-roll parameter in the contracting region (top-left), the slow-roll parameter in the bouncing region (top-right), the coupling function (bottom-left), and ω ( ϕ ) (bottom-right) in the non-minimal theory (4) are plotted for the chaotic bouncing model with α = 1 with respect to the e-N-fold time convention N ( N 2 2 ln a / a 0 for contraction, N is negative, whereas, at bounce, it becomes zero and then during expansion, it becomes positive). It is apparent that, during the contraction period, the scale factor behaves as a ( η ) ( η ) , and thus, ϵ b 2 . One also finds that the bouncing solution is asymmetrical. The coupling function f ( ϕ ) as well as ω ( ϕ ) become nearly unity in the reheating epoch, which arises shortly after the bounce.
Universe 07 00062 g001
Figure 2. The slow-roll parameter ϵ b (top-left), Hubble parameter H b (top-right), the coupling function f ( ϕ ) (bottom-left) and ω ( ϕ ) (bottom-right) in the non-minimal theory (4) are plotted for the chaotic bouncing model with α = 2 with respect to the e-N-fold time convention N ( N 2 2 ln a / a 0 , for contraction, N is negative, whereas, at bounce, it becomes zero and then during expansion, it becomes positive) during and after the bounce. One can see that the bouncing solution is asymmetrical. The coupling function f ( ϕ ) as well as ω ( ϕ ) become nearly unity in the reheating epoch, which appears shortly after the bounce.
Figure 2. The slow-roll parameter ϵ b (top-left), Hubble parameter H b (top-right), the coupling function f ( ϕ ) (bottom-left) and ω ( ϕ ) (bottom-right) in the non-minimal theory (4) are plotted for the chaotic bouncing model with α = 2 with respect to the e-N-fold time convention N ( N 2 2 ln a / a 0 , for contraction, N is negative, whereas, at bounce, it becomes zero and then during expansion, it becomes positive) during and after the bounce. One can see that the bouncing solution is asymmetrical. The coupling function f ( ϕ ) as well as ω ( ϕ ) become nearly unity in the reheating epoch, which appears shortly after the bounce.
Universe 07 00062 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

Nandi, D. Stability of a Viable Non-Minimal Bounce. Universe 2021, 7, 62. https://doi.org/10.3390/universe7030062

AMA Style

Nandi D. Stability of a Viable Non-Minimal Bounce. Universe. 2021; 7(3):62. https://doi.org/10.3390/universe7030062

Chicago/Turabian Style

Nandi, Debottam. 2021. "Stability of a Viable Non-Minimal Bounce" Universe 7, no. 3: 62. https://doi.org/10.3390/universe7030062

APA Style

Nandi, D. (2021). Stability of a Viable Non-Minimal Bounce. Universe, 7(3), 62. https://doi.org/10.3390/universe7030062

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