General background conditions for K-bounce and adiabaticity

We study the background conditions for a bounce uniquely driven by a single scalar field model with a generalized kinetic term K(X), without any additional matter field. At the background level we impose the existence of two turning points where the derivative of the Hubble parameter H changes sign and of a bounce point where the Hubble parameter vanishes. We find the conditions for K(X) and the potential which ensure the above requirements. We then give the examples of two models constructed according to these conditions. One is based on a quadratic K(X), and the other on a K(X) which is avoiding divergences of the second time derivative of the scalar field, which may otherwise occur. An appropriate choice of the initial conditions can lead to a sequence of consecutive bounces, or oscillations of H. In the region where these models have a constant potential they are adiabatic on any scale and because of this they may not conserve curvature perturbations on super-horizon scales. While at the perturbation level one class of models is free from ghosts and singularities of the classical equations of motion, in general gradient instabilities are present around the bounce time, because the sign of the squared speed of sound is opposite to the sign of the time derivative of H. We discuss how this kind of instabilities could be avoided by modifying the Lagrangian by introducing Galilean terms in order to prevent a negative squared speed of sound around the bounce.


Introduction
The inflationary scenario provides a model able to explain many of the observational features of the observed Universe, such as the cosmic microwave background radiation isotropy or the spatial flatness. An alternative scenario to solve the horizon problem is the possibility that the Universe underwent a contraction phase before the big-bang, i.e. a bounce. a e-mail: antonio.romano@udea.edu.co Several attempts have been made to understand if the bounce could be of a fully consistent cosmological model [1][2][3][4][5][6][7][8][9][10][11][12][13][14][15][16], and could indeed lead to viable alternatives to inflation.
In this paper we study the general conditions to realize the bounce with a K-essence type scalar field, without introducing any additional field. At the background level we impose the existence of two turning points where the derivative of the Hubble parameter H changes sign and of a bounce point where the Hubble parameter vanishes. We first obtain some general conditions on the function K (X ) and the potential V (φ) and then give the examples of two models constructed according to these conditions. One is based on a quadratic K (X ), and the other on a K (X ) which is always avoiding divergences of the second time derivative of the scalar field, which may otherwise occur. We integrate numerically the classical equations of motion for these models, verifying that they indeed have the expected evolution and can produce a bounce.
In the region where these models have a constant potential they became adiabatic on any scale, i.e. globally adiabatic, and because of this they may not conserve curvature perturbations on super-horizon scales. We then study the stability of the perturbations, focusing attention on the presence of ghosts and gradient instabilities. We confirm that gradient instabilities can arise around the bounce because the ghost free condition implies that the sign of the squared speed of sound is opposite to the sign of the time derivative of H . Finally we discuss how this kind of instabilities could be avoided by modifying the Lagrangian by introducing Galilean terms in order to prevent a negative squared speed of sound around the bounce.

K (X)-bounce
We will consider a scalar field with Lagrangian We denote the Planck mass with m p = (8π G) −1/2 , and we adopt units in which 8π G/3 = 1, so that the Einstein equation takes the form where the energy density and pressure are In the above equations we have introduced the function T (X ) because it is convenient to express the equations of motion in terms of it, and it plays a crucial role in defining the conditions for avoiding singularities in the solutions of the classical equations of motion. The energy conservation equations can then give the equation of motion of the scalar field, The minimally coupled scalar field equation of motion corresponds to the case when K (X ) = T (X ) = X . Another important quantity to determine the dynamics of the Universe is the second derivative of H , which after substituting the equation of motion forφ gives We now define turning points the local extrema of H (t), whereḢ = 0. From the Einstein equations we know thatḢ can be zero either because K = 0 or X = 0, and the corresponding value of the second order time derivativeḦ will play a crucial role in determining if these are actual extrema or flex points. From the previous equation after substituting, respectively, K = 0 or X = 0 we get This implies that X = 0 is a flex point and consequently can never be a turning point, while points at which K = 0 can be either local maxima or minima depending on the relative sign of V andφ. For a bounce with dynamics similar to the one shown in Fig. 1 to occur, the following three critical points are necessary: where we have defined and {t 1 , t 2 } correspond to the first and second turning point times, while t b is the bounce time. The construction of models satisfying the above conditions will be studied in the following sections.

Phase space dynamic of the bounce
In the previous section we have derived the conditions which have to be satisfied at the turning and bounce points, in Eqs. (11)- (13). In order to give a qualitative description of the phase space behavior of the models is convenient to define H 2 0 as the phase space locus of the points where the energy density is zero, i.e. the set of points satisfying the equation The bounce occurs when the phase space trajectory reaches H 2 0 , and as shown in Appendix A, the trajectory is always tangent to H 2 0 . The other important points are the ones where K = 0, which can correspond to a local minimum or maximum of H , depending on the sign of V φ , according to Eqs. (11) and (13).
It should be noted that even if the trajectory does not reach H 2 0 , the conditions for the local extrema of H could be satisfied. This can lead to an oscillating H as shown in Fig. 9, with oscillations that can continue even after the bounce. In phase space this behavior corresponds to oscillations around the turning point X T where K (X T ) = 0. The conditions in Eqs. (11) and (13) by them self do not guarantee the bounce, but only the existence of local extrema of H , which can produce oscillations of H . If somewhere between the turning points also the condition in Eq. (12) is satisfied, than there can be a bounce, or a sequence of bounces as shown in Fig. 9. The particular choice of K (X ), V (φ) and the initial conditions can determine different phase space behaviors, but we can use the conditions in Eqs. (11)-(13) in order to design models appropriately. In order to have turning points we need a kinetic term with at least one local extrema where K = 0, while for the choice of the potential there is more freedom as long as V φ has the right sign at the turning points according to Eqs. (11)- (13).
A more detailed analytical study of the phase space trajectories and their dependence on the initial conditions goes beyond the scope of the present paper, which is mainly focusing on deriving the general conditions for the bounce and turning points contained in the previous section. We devote the rest of the paper to the construction of models specifically designed to satisfy those conditions and show by numerical integration of the differential equations that the expected phase space behavior is realized. It should be noted that the conditions we obtained in Eqs. (11)-(13) are completely general since we do not assume any specific form of K (X ) contrary to the analysis performed in previous work such as [17,18], where it was assumed K (X ) = X n . Another important difference is that we consider models with only one scalar field, without any additional matter field or perfect fluid, and find the most general conditions for a bounce uniquely driven by the scalar field with a dynamics depicted schematically in Fig. 1.

An example: quadratic bounce
Based on the background conditions derived in the previous section we can introduce a quadratic kinetic term model and an appropriate potential which will give the correct dynamic for H (t): In the rest of the paper we will consider a model with The phase space behavior of the model defined in Eq. (17) is shown in Figs. 2, 3, 4, 5, 6, and 7. As can be seen in Fig.  4 the phase space trajectory is tangent to the zero energy locus H 2 0 at the bounce time, as proved in the appendix. The turning points correspond to X = X T = X (t 1 ) = X (t 2 ), and the sign of V iφ i determine, according to Eqs. (11) and (13), if the turning point is a local minimum or maximum for H (t).

Possible divergences of the second derivative of φ
In this section we will try to construct some other model satisfying the conditions in Eq. (14). Let us consider models of the type where X T is the value of X at which the turning point for H should occur, i.e. where K (X T ) = 0, based on the considerations in the previous sections. Note that in this case we have X T = X (t 1 ) = X (t 2 ), but V 1φ1 = V 2φ2 in order to satisfy the conditions for the turning points given in Eqs. (11) and (13). The problem with this kind of model is that there is a line in phase space plane where T (X c ) = 0, which corresponds to a divergence of the second derivative of φ. At that point the equation of motion becomes a first order equation, which correspond to a constraint of the value of the field φ c when X = X c . This critical point in phase space is defined by the following conditions: These equations can be solved to give To avoid this kind of divergence we will introduce a new class of models in the next section.

Models avoiding the divergence of the second derivative
In order to avoid the degenerate solutions encountered in the model studied in the previous section we consider a class of models for which T (X ) is never zero where a is a constant. Solving for K (X ) the above differential equation we get where b and c are two arbitrary integration constants. The socalled cuscuton model [19] corresponds to the case in which a = 0. In this case both T and T vanish. In order to obtain a K (X ) with a local minimum we need a negative b < 0, leading to a minimum at Note that an oscillatory H (t) can be obtained for an appropriate choice of the initial conditions, leading to a sequence of consecutive bounces, or to an oscillating H (t). The linear part of K (X ) can give an ekpyrotic phase, and it arises naturally by imposing the condition T (X ) = a > 0. The phase space behavior of the model defined in Eq. (17) is shown in Figs. 9, 8, 10, and 11. As can be seen in Figs. 9 and 8 the phase space trajectory is tangent to the zero energy locus H2 0 at the bounce time, as proved in the appendix. As in the case of the model defined in Eq. (17) the turning points correspond to X = X T = X (t 1 ) = X (t 2 ), and the sign of V iφ i determine, according to Eqs. (11) and (13), if the turning point is a local minimum or maximum for H (t). For the purpose of integrating numerically the equations of motion we consider this specific model

Stability of cosmological perturbations and adiabaticity
So far we have focused on the conditions to realize the bounce at the background level, but in order for these models to be viable cosmological models which could provide an alternative to inflation, we need to check if the cosmological perturbations are stable. In particular there could be problems with the growth of anisotropies [14] and gradient instabilities or ghosts [20] could arise. The quadratic action for scalar perturbations in the comoving slice gauge for a general Lagrangian of the form K (X, φ) + G(φ) can be written as [20] S (2) where we have defined and the squared sound speed for the scalar perturbations is The condition for absence of ghosts [20] is A > 0, while gradient instabilities are avoided when the square of speed of propagation of perturbations c 2 s is always positive. For the class of models we considered, corresponding to G = 0, these conditions imply We can see that an additional advantage of the class of models satisfying T (X ) = a > 0 given in Eq. (28) is that they are ghost free by construction, beside avoiding a singularity of the classical equations since as discussed previously T (X ) is also the coefficient ofφ in Eq. (6). Nevertheless around the bounce time, whenȧ(t b ) = 0, gradient instabilities could arise [21] becauseȧ is negative before and positive after t b , which impliesḢ > 0 in some neighborhood of t b . From the Einstein equation (2) we can then immediately deduce that gradient instabilities will necessarily arise in any K (X ) model which is also ghost free, i.e. also satisfy T (X ) > 0, because in this case the sign ofḢ is the opposite of the sign of c 2 s , These kinds of instabilities could be healed by adding Galilean terms of the form G(X, φ) φ, to the Lagrangian [20], because in this case the above relation would not be valid anymore. Here we can note that a modification consisting in adding terms of the type G(φ) would not be sufficient because in this case we have The last equation implies that also in this case c 2 s must be negative around the bounce, sinceḢ is positive around the bounce and D > 0 in order to avoid ghosts. Consequently also for these models, the ghosts free condition implies the relation sgnḢ = − sgn c 2 s , and gradient instabilities are unavoidable around the bounce. Only introducing a general G(X, φ), which is also a function of the kinetic term X , the gradient instabilities can be avoided, but even in this case a careful analysis is necessary [24][25][26].
Another general problem of K (X ) models related to the sign of c 2 s is that around the bounce the Hamiltonian could be not positive definite when the gradient term dominates the kinetic term. Also this problem could be avoided by adding a Galilean term of the form G(X, φ) φ to the Lagrangian in order to avoid a negative sound speed, according to Eq. (33), while still realizing the conditions for the bounce. We will leave to a future upcoming work a detailed study of the stability of cosmological perturbation for K (X ) models with the addition of general Galilean terms, since the purpose of the present paper is mainly the derivation of general conditions to realize the bounce at the background level with a K (X ) model. It should be noted nevertheless that these types of instabilities are general for any K (X ) driven bounce model [21], and can arise also for the other models which were previously studied in [17,18].

Adiabaticity
Due to the form of the Lagrangian we can anticipate that the behavior of curvature perturbations for these models requires a careful treatment. The non-adiabatic pressure perturbation for these kind of scalar field models is in fact zero on any scales [23] in the region where the potential is constant, because in that regime they are equivalent to a barotropic fluid with equation of state P(ρ) so that For this reason they are called globally adiabatic(GA) models, to distinguish them from attractor models which are adiabatic only on super-horizon scales (SHA).
After defining the perturbed metric as in [24] it was shown that in general relativity, in absence of anisotropies, the conservation of the perturbed energy momentum tensor gives which is true on any scale, and is not based on neglecting any gradient term. For GA models the adiabatic sound speed c s = δ P ∂ρ c is equal to the phase speed c w =Ṗρ and for this reason, despite δ P nad = 0, nothing can be inferred about the behavior ofṘ c , otherwise R c should be constant on any scale. Since in ekpyrotic models the potential is not constant during the bouncing phase, c s = c w . The exponential growth of R c [14] then implies that δ P nad = 0, i.e. the bounce is strongly non-adiabatic. This large non-adiabaticity is in fact the cause of the growth of anisotropies during the bouncing phase.
On the contrary, in the region where the potential asymptotes a constant value, these models could enter a phase of ultra-slow inflation in which curvature perturbations may not be conserved on super-horizon scales [23]. This is due to the fact that when the potential is constant c s = c w and conservation of curvature perturbation on super-horizon scales is not guaranteed anymore [24] despite adiabaticity. A similar ultra-slow phase could occur for the model defined in Eq. (17), in the region where the potential tends to a constant value.

Conclusion
We have derived the general conditions to obtain a bounce with a scalar field model with generalized kinetic term of the form K (X ) in general relativity. The requirement of the existence of a bounce point and of two turning points for H (t) give some conditions which K (X ) and the potential V (φ) have to satisfy. We have given the examples of two models constructed according to these conditions. One is based on a quadratic K (X ), and the other on a K (X ) which is always avoiding divergences of the second time derivative of the scalar field, which may otherwise occur. An appropriate choice of the initial conditions can lead to a sequence of consecutive bounces.
While at the perturbation level one class of models is free from ghosts, in general gradient instabilities are present around the bounce time, and they could be avoided only by modifying the Lagrangian. In the future it will be interesting to add Galilean type terms to the Lagrangian and study in detail the cosmological perturbations during and after the bounce, and the effects of global adiabaticity in the region where the potential tends to a constant value.