Cosmological bounce and some other solutions in exponential gravity

In this work we present some cosmologically relevant solutions using the spatially flat Friedmann-LemaitreRobertson-Walker (FLRW) spacetime in metric f(R) gravity where the form of the gravitational Lagrangian is given by 1 α e. In the low curvature limit this theory reduces to ordinary Einstein-Hilbert Lagrangian together with a cosmological constant term. Precisely because of this cosmological constant term this theory of gravity is able to support nonsingular bouncing solutions in both matter and vacuum background. Since for this theory of gravity f ′ and f ′′ is always positive, this is free of both ghost instability and tachyonic instability. Moreover, because of the existence of the cosmological constant term, this gravity theory also admits a de-Sitter solution. Lastly we hint towards the possibility of a new type of cosmological solution that is possible only in higher derivative theories of gravity like this one.


Introduction
Investigation of non-singular bouncing cosmological solutions to Einstein field equations has a history that dates back to first half of the twentieth century, and can be attributed to various works of Lemaitre, Tolman, Friedmann and even Einstein himself (see, e.g., Ref. [1] for early histories of bouncing cosmology).But it was then largely considered just as an alternative solution to Einstein equations and did not have any physical motivation as such.The most accepted paradigm about the universe evolution then was the big-bang paradigm, which is based upon the existence of a curvature singularity in the past.However, big-bang singularity was plagued with several other problems.Inflationary scenario came into picture in the later half of the nineteenth century as a very promising candidate to solve these problems.This paradigm was pioneered by the works of Guth [2] and Linde [3].There are many models in the literature so far that realize an inflationary scenario (see [4] for a comprehensive review of all the inflationary models).
Although inflationary cosmology is highly successful in explaining various features of the early universe, it is still plagued with the issue of singularity 1 , which, by definition, is a state of physical lawlessness.When we try to describe our universe with the available physical theories, we usually do not want a singularity to come into the picture.This was the physical motivation which refuelled the interest in nonsingular bouncing scenarios.There are various ways to realize a bouncing scenario (see [6] or [7] for a comprehensive review). 1 The reader can consult Ref. [5] for a review on singularities in inflationary cosmology.
Universe 2018, 4, 105; doi:10.3390/universe4100105www.mdpi.com/journal/universe If one tries to realize a bouncing solution for spatially flat FLRW metric in general relativity (GR), one needs to invoke null energy condition (NEC) violating matter components like ghost fields, ghost condensates or Galileons.If one does not wish to invoke such exotic matter components and still wants to realize a bouncing solution in spatially flat FLRW metric, then he/she has to resort to modified gravity.Modifications to general relativity at high curvature regime near a curvature singularity can indeed be expected.When quantum corrections or string theory motivated effects are taken into account, then the effective low energy gravitational action indeed admits higher order curvature invariant terms [8,9].The simplest of such modifications is when the correction terms depend only on the Ricci scalar R. In such cases the Einstein-Hilbert Lagrangian R is modified to f (R), a function of R (see [10][11][12][13] for beautiful reviews on f (R) gravity).
There has been many attempts to realize bouncing scenario in f (R) theories of gravity.It was first pointed out in Ref. [14] that R + αR 2 gravity with a negative α can give rise to a bouncing scenario in spatially flat FLRW metric.Some authors have even used f (R) gravity to tackle the issues related to cosmological bounce and a cyclic universe [15,16].Bouncing cosmology for quadratic and cubic polynomial f (R) gravity theory was more recently worked out in [17,18] respectively, from both Jordan frame and Einstein frame point of view.It was shown in [17] that for spatially flat FLRW metric a bounce in the Jordan frame is never accompanied by a bounce in the Einstein frame when hydrodynamic matter in the Jordan frame satisfies the condition ρ + P ≥ 0.Here ρ an P specifies the energy density and pressure of matter in the Jordan frame.Working solely in the Jordan frame, it was shown in [19] that quadratic gravity theories of the form λ + R + αR 2 with a negative α and monomial gravity theories of the form R 1+δ can also produce bouncing cosmologies.Carloni et al. in [20] presented the bouncing conditions in f (R) gravity and also analyzed the conditions for R 1+δ , R + αR m and exp αR type of gravity.
Two necessary conditions for the physical viability of any f (R) theory are f (R) > 0 and f (R) > 0. As we will see in a later section, all the other f (R) theories except e αR (α > 0) that has commonly been considered in the literature so far to realize a bounce in spatially flat FLRW metric cannot have f (R) > 0 and f (R) > 0 simultaneously for all R.So exponential gravity may be the only physically viable candidate for achieving a bounce.Carloni et al. in [20] concludes that this gravity theory can give rise to a bouncing solution only in a closed FLRW universe, but as we will show later on, this theory can also produce a bounce in the flat FLRW universe.In previous studies exponential gravity has been used extensively used to study cosmological inflation and late time acceleration of the universe [21][22][23].The motivation of the present paper is to present a metric f (R) theory of gravity, which is free of the above mentioned instabilities, and which can produce a successful cosmological bounce in the early universe.The model of bounce which we present here is to be taken as an effective theory in high energy scales as exponential f (R) only describes the system very near the bounce point.Essentially we present a bounce mechanism and not a full description of cosmology which includes how physics much prior to bounce is related to the bouncing period, although our model can have a transition to low energy GR for small values of the Ricci scalar but we presume such a cross-over to low energy theory may require new physics.As specified in Ref. [24] the prebounce cosmology can be connected to a particular bounce mechanism and out of various bounce mechanism f (R) bounce is one.The inclusion of higher order terms in the Ricci scalar R in the action is often motivated by two main observations: first, adding new terms in the space-time curvature could explain the observations typically associated with dark matter and/or dark energy, and second, since the Einstein-Hilbert action is not renormalizable, any consistent theory of quantum gravity is expected to contain higher order curvature terms in the action that become important near the Planck scale.While f (R) theory is not the most general kind of such theory, it is one of the simplest modifications possible on GR and are often viewed as a good first step in understanding the effect of adding additional terms to the Einstein-Hilbert action.Consequently we assume that the prebounce phase of our theory will be some form of cosmological theory as discussed in [24] or in Ref. [25].In [25] the author takes up a case where there is a bounce followed up by inflation, in the present we do not expect inflation after bounce in exponential gravity but exponential gravity do have an exact de Sitter solution and we briefly discuss cosmology of the de Sitter solution later in this article.
In this article, we briefly study the nature of scalar metric perturbations and opine briefly on the tensor perturbations through bounce.The perturbation theory is completely presented in the Jordan frame.It is shown that scalar perturbations do remain non-singular near the bouncing point but there can be new instabilities in the new system which can make some some of the modes to become non-perturbative near the bounce.We also show that we do not expect a higher tensor-to-scalar ratio in our bounce model.
In this paper, we also present a new solution of f (R) gravity theories which admits a cosmological bounce.The new solution is related to more degrees of freedom of f (R) theories.Unlike GR, f (R) theories depend on the second time derivative of the Hubble parameter and one can tune the cosmological development of a model by specifying various values of Ḧ at some specific time.These new solutions can produce interesting new model universes.We present a new solution in bouncing f (R) theories where the universe transits from an decelerated expansion phase to a normal expansion phase vis a contraction phase.The new result presented in the present paper is general and we have given some specific examples using exponential gravity as an example.
The material in this paper is presented in the following way.In the next section, we present the formalism of f (R) theory in the two conformal frames, the Jordan frame and the Einstein frame.The relation between the frames is presented in the this sections.The bouncing scenario is described in Section 3. Section 4 specifies bounce in exponential gravity.In this section, we present the numerical results depicting various kinds of bounces.Discussion on scalar metric perturbation through bounce is presented in Section 5. Two exact solutions of exponential gravity is presented in Section 6.We discuss the new solutions in f (R) gravity theories in Section 7. The next section concludes the article by summarizing the results obtained.

The Cosmological Set Up in the Two Frame
In this section, we present the formal structure of the theory we will pursue in this article.The field equations in the two conformal frames and the formulae connecting the frames are briefly specified in this section.

The Jordan Frame
In the Jordan frame the field equations for f (R) gravity is given in the tensorial form as follows, For the flat FLRW metric given by and for a perfect fluid given by the energy-momentum tensor the modified Friedmann equations for f (R) gravity are: where ρ eff and P eff are related to energy density and pressure arising from curvature.They are defined by the following expressions, For a barotropic fluid, where ω = 0 corresponds to dust and ω = 1 3 corresponds to radiation.It must be noted that the 4-velocity u µ in Equation ( 3) is the normalized 4-velocity of a fluid element.The other relevant dynamical equation is the continuity equation for the hydrodynamic matter component

The Einstein Frame
The equivalent Einstein frame description of f (R) gravity is defined in terms of a conformally related metric where Let us define in Einstein frame a scalar field φ and the potential V(φ) as Gravitational dynamics in Einstein frame is affected by the dynamics of the scalar field which comes into existence in the Einstein frame due to the conformal transformation.Under the conformal transformation in Equation (10) the energy-momentum tensor transforms as Specifically, we have the relations In the Jordan frame there is only one hydrodynamic matter component, where as in the Einstein frame there is also a scalar field defined above, which has a potential given by Equation (11), which couples non-minimally with the hydrodynamic matter component.Since in the Einstein frame the gravitational theory is GR, we can write the field equations in the tensorial form as where Tµν is the energy-momentum tensor of the hydrodynamic fluid in the Einstein frame Tµν = ( ρ + P) ũµ ũν + P gµν , and T µ ν is the energy-momentum tensor of the scalar field φ.Here Rµ ν is the corresponding Ricci tensor in the Einstein frame.We can recast the conformally related Einstein frame metric gµν as an FLRW metric by using the redefinitions The Hubble parameters of the Jordan frame metric and the Einstein frame metric are now related by the equation where H = ã ( t) ã( t) , prime now stands for d/d t.We can now write the Friedmann equations in the Einstein frame where The other relevant dynamical equations in Einstein frame are the Klein-Gordon equation for the scalar field and the continuity equation for the hydrodynamic matter component It is seen from the above equations that the scalar field and hydrodynamic matter components in the Einstein frame satisfy coupled differential equations.
Before we end this section, we want to point out that we expect physics to be the same as observed from both the conformal frames.Apparently the two conformal frames may not look the same but the features that seem to be different in the two frames may be due to the fact that one needs to transform also the units he/she uses, between one frame and the other.

Description of a Bouncing Scenario
In this section, we describe the scenario of a cosmological bounce in Jordan frame from the point of view of both the frames.A cosmological bounce for the homogeneous and isotropic FLRW universe is defined mathematically by the conditions H b = 0 , and Ḣb > 0 . ( where the subscript b on a time dependent quantity denotes the value of the quantity at the time of the bounce.Using Equation (4), the first bouncing condition in the Jordan frame becomes Throughout the article, we will assume that the matter component satisfies the conditions ρ ≥ 0 , ρ + P ≥ 0 , (25) in the Jordan frame 2 .If matter satisfies the above condition, then the second bouncing condition in the Jordan frame becomes Solving the modified dynamical equations, Equations ( 5) and ( 9), in the Jordan frame, while keeping in mind the constraint, Equation (4), requires three conditions H(t = 0), Ḣ(t = 0) and Ḧ(t = 0).If we choose t = 0 to be the bouncing moment, then the bouncing conditions require H(0) = 0 , and Ḣ(0) > 0 . ( However, Ḧ(0) can be positive, negative or zero.If Ḧ(0) = 0, then we have a completely symmetric bounce, i.e., the contracting phase of the bounce is a mirror image of the expanding phase of the bounce.If Ḧ(0) > 0 ( Ḧ(0) < 0) then the evolution of the scale factor is steeper in the expansion (contraction) phase.
Let us now try to visualize how a Jordan frame bounce looks like from the Einstein frame.Any cosmological evolution in the Einstein frame is governed by the nature of the potential V(φ) and how the scalar field moves on the potential.From the definition of the scalar field φ, we see that and Therefore the Jordan frame bouncing conditions applied on Ḣb and Ḧb determines the Einstein frame values of φ b and φ b 3 at t = 0. Note that if the viability conditions f > 0, f > 0 of an f (R) gravity is respected, then the sign of Ḧb determines the sign of φ b , which is in turn related to the time symmetrical or asymmetrical nature of the bounce.The H b = 0 condition, from Equation (17) implies Before we close this section, we would like to present a short discussion on the difficulties of attaining a bouncing solution in f (R) theory of gravity.One may work with f (R) = AR 1+δ where A and δ are constants.For stability From the bouncing condition in the Jordan frame one can easily check that for a successful bounce one requires Aδ ≤ 0 if we assume that ρ ≥ 0 for the hydrodynamic fluid.These conditions make f < 0 near the bounce where we know that R b > 0 and consequently f (R) = AR 1+δ cannot describe a gravitationally stable bounce.In the quadratic model we have f (R) = λ + R + αR 2 where λ , α are constants and for a bounce α < 0. In this case it is obvious that f < 0 throughout the bounce making the theory unstable.Even cubic gravity models, where f (R) = R + αR 2 + βR 3 where α, , β are constants and for bounce α < 0 and β > 0, can accommodate cosmological bounces.One can tune the parameters in such a way that f > 0 for all values of R [26] by fixing the values of the parameters, but then it is seen that 2 We are not calling this condition as the weak energy condition as the energy conditions are generally stated in the Einstein frame.3 We want to remind the reader at this point that a prime on f implies a derivative with respect to R where as a prime on φ implies a derivative with respect to t.
there are two branches corresponding to positive and negative values of f [18].The unstable branch remains a reality in such theories and the theory becomes more complex as there appears a point where f = 0 which is a singular point in f (R) gravity.In this way one can go on to show that even f (R) = R + αR m where α and m are real constants and m is an integer greater than one gives rise to an unstable bounce.It was shown in Ref. [18] that no polynomial f (R) gravity can simultaneously void both ghost and tachyonic instability for all R.All these examples show that it is very difficult to find a f (R) which gives rise to a stable cosmological bounce.In this regard we show that the exponential f (R) gravity theory, as chosen in the present article, can produce perfectly stable cosmological bounces.Our model is explained in the next section.

Bounce in Exponential Gravity
In this section, we will study cosmological bounce in f (R) gravity where the form of f (R) is given by where α > 0. Note that for this gravity theory implying that this theory is free from both ghost instability and tachyonic instability for all values of R.Moreover, so we can see from the bouncing condition in Equation ( 24) that in this case both matter and matter less bounce is possible.For bounce in presence of matter, whereas for a matter less bounce The cosmological constant term plays a crucial role in producing the bounce If we want to remove the cosmological constant term by taking a theory like which is always positive for any positive value of R.This implies that this theory of gravity cannot support a cosmological bounce in either matter or vacuum background.It is precisely the cosmological constant term that helps to achieve a bounce.The Einstein frame scalar field for exponential gravity is whose potential in the Einstein frame is given by The shape of the above potential is shown in Figure 1.The value of α = 10 12 in our system of units where all the dimensional parameters are ultimately expressed in units of the Planck mass M P .For convenience we assume M P = 1 throughout the article.The potential crosses the V(φ) axis at φ = 0 and the φ axis at φ = 3 2κ .As the gravitational part of the theory, describing the bounce, becomes essentially GR in the Einstein frame, it is easier and interesting to start from the Einstein frame description and track the bounce as the movement of the scalar field on the scalar potential.Firstly, note that the values R = 0 and R = 1/α in the Jordan frame corresponds to the values φ = 0 and φ = 3/(2κ) in the Einstein frame.Since R b has a limited range for bounce in presence of matter where the range is 0 < R b < 1/α in the Jordan frame, φ b is limited in the range 0 < φ b < 3/(2κ) in the Einstein frame.This means that if we want to impose the Einstein frame intermediate 4 conditions at the time of bounce in the Jordan frame (t = t = 0), then the scalar field in the Einstein frame can start from only a very small part of the curve in Figure 1, which is in the fourth quadrant.We may recall from our discussion of the previous section that the sign of φ b , i.e., whether the field is rolling up or down the potential when the bounce has happened is related to whether the time evolution of the scale factor is completely symmetric or steeper on one side and flatter on the other.The dependence of the nature of the time evolution on the initial conditions at bounce is elaborated by Figures 2 and 3.
-   We call the conditions as intermediate instead of initial conditions.The reason being that we impose our conditions on the dynamical system at t = t = 0 and look at the system at both negative and positive times. - In both the figures we show a symmetric bounce, where the intermediate conditions are specified on the figure captions.As exponential gravity can accommodate bounces in presence of hydrodynamic matter and vacuum the two figures show two different kind of universes, one filled with radiation and the other devoid of matter.For symmetric bounces we wee that φ (0) = 0 which translates to Ḧ(0) = 0 in the Jordan frame.For purely symmetrical bounces one must have Ḧ(0) = H(0) = 0 and Ḣ(0) > 0 in the Jordan frame.The last statement is true for all kinds of bounces in metric f (R) theory.
One can also model asymmetric bounces in f (R) theories.In our case we show two asymmetric bounces in Figures 4 and 5.The bounce depicted in Figure 4 is assisted by radiation where as the other bounce takes place in vacuum.In the case of bounce in presence of radiation one can easily calculate from the intermediate conditions that Ḧ(0) = (2/(3e 2 )) × 10 −19 and for vacuum bounce one gets Ḧ(0) = (2 √ e/3) × 10 −19 .In these cases one has non-zero values of the second time derivative of the Hubble parameter in the Jordan frame.In all of the bounces we have discussed in this section φ in the Einstein frame remains positive in sign which implies that R > 0 for all the bounces, in the Jordan frame.In the effective theory approach exponential gravity becomes similar to GR near small R values when one can neglect the higher powers (starting from the quadratic one) of R in the exponential f (R).The theory presented in this paper is consistent when we use the theory for positive non-zero values of R in the Jordan frame.During the end stages of the bouncing scenario, in all the cases, R tends to zero by which one is very near the GR limit.In the bouncing scenarios presented in this article we do not specify how the effective f (R) theory can transform to GR near R ∼ 0, we hope some new physics is involved during this phase.
-400 000 -200 000 200 000 400 000 t On the other hand if we do not take the f (R) theory presented in our article as some form of effective theory, which can change cosmology only for high values of the Ricci scalar, then something interesting happens.If the bounce is not due to some effective change in the gravitational part of the Lagrangian, and the f (R) modifications are valid for all values of R including R = 0 then there are new possible dynamical configurations of the universe.These new configurations depend upon the intermediate conditions applied to produce the dynamics of the bounce.

Evolution of Metric Perturbations through the Bounce
In this section, we discuss the evolution of scalar cosmological perturbation through a bounce in f (R) gravity.We also briefly opine on tensor perturbations in bouncing cosmologies guided by f (R) theories before we end this section.We express the scalar perturbation equations in terms of conformal time defined as follows, In the domain of linear perturbations, the scalar perturbed FLRW metric has two gauge invariant degrees of freedom.In the longitudinal gauge this can be expressed as, Here Φ and Ψ are the two gauge invariant perturbation degrees of freedom, also called the Bardeen potentials.If only a single barotropic matter component is present, the perturbation in the matter sector can be assumed to be adiabatic, so that the sound velocity can be defined as where ω is the constant equation of state of the barotropic matter component.In this section, we will be working with conformal time and a derivative with the conformal time will be represented by a prime above.Consequently we will specify the derivatives of f (R) with respect to R as f R and higher derivatives with respect to R as f RR and so on 5 .The 00, ii, 0i, ij − th(i = j) components of the linearized f (R) field equations in the Fourier space are [27] f The new conventions which are at odds with our previous convention becomes necessary as conformal time is involved in the discussions.From the next section we will use the old conventions.
Here we have defined the matter density perturbation as δ = δρ ρ and the perturbed velocity potential as δu i = 1 a ∂ i v.There are a total of four dynamical scalar perturbation quantities, namely, Φ, Ψ, δ, v, and two constraint equations between them, namely, Equations ( 42) and (44).Therefore only two of them are independent.For our convenience, we can take them to be the two metric degrees of freedom Φ, Ψ.
From Equations ( 42) and (43), we can eliminate δ to write, The equations Equations ( 45) and ( 46) can be solved to get the solutions Φ(η) and Ψ(η).Once Φ, Ψ is known, δ and v is determined from Equations ( 42) and (44).We calculate the gauge invariant comoving curvature perturbation as For the case of vacuum the right hand sides of Equations ( 42)-(44) all vanish.In this case we can just solve the two coupled second order differential equations Equations ( 43) and (45) to obtain Φ(η) and Ψ(η).In this case the comoving curvature perturbation is just R = Ψ. (48)

Scalar Perturbation Evolution through Bounce
In this section, we present the numerical solution of the perturbation equations by assuming suitable initial conditions.Let us first discuss some general conclusions regarding the forms of the solutions of the perturbation equations in the context of a nonsingular bounce.First of all note that the Hubble radius diverges at the bounce (see Figures 6 and 7).Therefore as the bounce approaches, all the perturbation modes become sub-Hubble near the bounce.We consider perturbation modes with wavenumber k = 10 −8 (in Planck units) to illustrate our point regarding perturbation growth near the bounce point.
-400 000 -200 000 200 000 400 000 t In general for all models of non-singular bouncing cosmologies, near the bounce point when η ∼ 0, the scale-factor can be well approximated by an even function of the conformal time.In this section we will be working in the vicinity of the η = 0 point and consequently assume the scale-factor to be an an even function.As the scale-factor is an even function of time, the Hubble parameter would be an odd function of conformal time, as it is first order derivative (with respect to η) of the former.Likewise, R(η) and all the functions of it, like f R , f RR etc. would be even functions of time.
Below we constrain ourselves to the case of a bouncing solution in radiation background, although our considerations hold true also for a bounce in vacuum background.Both the perturbation evolution equations (Equations ( 45) and ( 46)) that we use to solve for Φ(η) and Ψ(η) are of the form where φ 0 (η), φ 1 (η), φ 2 (η), ψ 0 (η), ψ 1 (η), ψ 2 (η) are functions defined by the background evolution and D η ≡ d dη .In case of a completely time-symmetric evolution near η = 0, φ 0 (η), φ 2 (η), ψ 0 (η), ψ 2 (η) are even functions whereas φ 1 (η), ψ 1 (η) are odd functions.For example, for Equation (45) we have whereas for Equation (46) we have The above coefficient functions are seen to be non-singular near η = 0 showing that the perturbation evolution equation, as given in Equation (49), do not produce singular solutions near the bounce point.Since D η → −D η when η → −η, the perturbation equations are time-symmetric, and have even solutions Φ e (η), Ψ e (η) and odd solutions Φ o (η), Ψ o (η).Note that, since the perturbation evolution equations are linear differential equations, in general any linear combination of the even and odd solutions c 1 Φ e + c 2 Φ o and c 1 Ψ e + c 2 Ψ o is also a solution of the perturbation evolution equations.The odd solutions Φ o (η), Ψ o (η) vanish at bounce, meaning that if these solutions start with an absolute value less than unity prior to the bounce, then it remains less than unity throughout the bouncing phase.In other words, these odd solutions, if they start at a perturbative level before the bounce, they remain at the perturbative level throughout the bounce.On the contrary the even solutions Φ e (η), Ψ e (η) must have a local extremum at bounce.
Let us discuss these even solutions in more detail.Suppose these perturbations start with an absolute value less than unity prior to the bounce.If we want to guarantee that the perturbations indeed always remain in the perturbative level then one can specify a general condition as follows.
If the relevant perturbation (Φ e (η) or Ψ e (η)) is in the perturbative region (meaning |Φ e (η)| < 1 or |Ψ e (η)| < 1) during the contracting phase then the perturbation will remain perturbative throughout the bounce if the minimum of the function has a value greater than −1 at η = 0 or the maximum of the function is less than +1 at η = 0.A major part of the above condition is encoded in the following inequalities The above inequalities only form a subset of the various cases which may arise from the condition stated earlier.The above inequalities can be used for a simplified consideration of the problem at hand.Solutions for perturbations are found out by solving two coupled differential equations, each having the form as in Equation (49).Since Φ e (0) = Ψ e (0) = 0, the two perturbation evolution equations written at η = 0 relates the four quantities Φ e (0), Ψ e (0), Φ e (0), Ψ e (0) by two simultaneous algebraic equations.For the numerical initial conditions that we have chosen (as mentioned in the caption of Figure 2) and for k = 10 −8 these two simultaneous algebraic equations are as follows: 1.82Φ (0) − 4.11 × 10 −13 Φ(0) + 1.82Ψ (0) + 4.60 × 10 −13 Ψ(0) = 0. ( These two algebraic equations can then be solved to obtain Φ e (0), Ψ e (0) in terms of Φ e (0), Ψ e (0) as follows Ψ e (0) = −0.22 × 10 −13 Φ e (0) − 1.11 × 10 −13 Ψ e (0).
Let us now check for the four conditions mentioned above which guarantee that if the perturbations started with a perturbative value, then they remain perturbative throughout the bounce.Straightforward check gives that the first condition and the third condition can never be satisfied, at least in the simplified cases we are discussing.The second condition and the fourth condition can be satisfied respectively for Φ e (0) > 0 , −0.2Φ e (0) < Ψ e (0) ≤ 0 and Φ e (0) < 0 , 0 ≤ Ψ e (0) < −0.2Φ e (0).For the perturbation theory to be valid throughout the bouncing time period, the perturbations need to remain always at a perturbative level.The work out in this brief section shows that there are various cases where the perturbations in the longitudinal gauge remains perturbative throughout the bounce.We do not rule out cases where the perturbations may become non-perturbative near the bounce but our analysis shows that the perturbations do not diverge near the bounce point.Before we end this section we want to present a discussion on the instabilities present in general bouncing models based on GR with modified matter sector (as ghosts, Galileon) and f (R) theories with standard matter sector.In GR-based theories the non standard matter component violates the null energy condition (NEC) near the bounce point and consequently this kind of matter component is essential for bounce in spatially flat FLRW models.But this NEC breaking phase may produce instabilities as ghost instabilities.One may eradicate the ghost instability by considering bounce in Galileon theory or in general Horndeski theory.The Horndeski-like theories show gradient instabilities due to which the sound speed squared turns out to be negative and there is an exponential growth of perturbations.The following Refs.[28][29][30] give a detailed discussion on the issues of instabilities arising in GR based bounces and methods to eradicate them.In f (R) bounce theories one generally do not employ non-standard matter, the energy density and pressure for the non-standard matter component is produced from the scalar curvature and encoded in ρ eff and P eff , as shown in Equations ( 4) and ( 5).The matter part specified by ρ and P generally always has ρ + P > 0. During bounce the curvature contribution (to energy density and pressure) acts like non-standard matter in f (R) theory.In Ref. [10] the authors show that in f (R) theory, which is a higher derivative theory, ghosts do not appear.Although in exponential gravity one naturally takes care of the conditions f (R) > 0 and f (R) > 0 which eradicates inherent instabilities of f (R) dynamics there may remain some hidden form of instabilities corresponding to the gradient instability which may make some of the scaler modes to become non-perturbative during the time of bounce.The bounce problem in f (R) theories are still in their infancy and we hope in the near future a stability analysis of the perturbation modes will be presented in full detail.

Brief Comments about Tensor Perturbations through Bounce
The calculations on tensor perturbations in f (R) gravity, in the context of inflation, has been done in Ref. [11].From the above reference one can see that in general the tensor perturbation evolution equation, in coordinate time, is given by ḧλ + (a 3 F) for λ = +, ×, which are the polarization modes of the tensor perturbations.Here F = d f /dR.The evolution equation for the tensor perturbations follows the same generic form of perturbation equation presented in Equation ( 49), as the coefficient function multiplying ḣλ is odd and the coefficient function multiplying h λ is even in coordinate time.This shows that in principle one can produce a similar analysis, of the tensor perturbations in the present case, based on the symmetry of the gravitational wave amplitudes in time.We do not expect any singularity of the tensor modes as the above equation does not have any singular point in the bouncing regime.In Ref. [31] the authors showed that the two scalar field matter bounce model in GR produces a tensor-to-scalar ratio(r) that is within the observed upper bound.In f(R) bounce models, however, explicit analytical calculation of tensor-to-scalar ratio is still missing.We expect that the effect of modified f (R) gravity models are valid near the vicinity of the bounce and if the previous history of the universe does not produce large gravitational power spectrum (compared to comoving curvature power spectrum) then exponential gravity will not amplify r (the tensor-to-scalar ratio).In Starobinsky inflationary scenario where an R 2 correction term is incorporated in the Lagrangian, r comes out to be smaller than the observational bound [11].Also, in Ref. [32] the authors consider inflation in exponential gravity and calculate r to be below the observational limit.From these facts, we can expect that inclusion of such correction terms in the Lagrangian will not abruptly increase the value of r (from observational bound).This point needs explicit analytical consideration, which we hope to address in future publications.
In the next section we will discuss some interesting exact solutions in exponential gravity.The solutions plotted in this section are obtained numerically, rarely in f (R) gravity theory we have the privilege of having exact solutions [33,34].The solutions given in the next section do not define the form of f (R), the solutions are exact solutions of exponential gravity.

Two Exact Solutions in Exponential Gravity
In this section, we show that (1/α)e αR gravity has both an exact bouncing solution and a solution with constant Ricci curvature.As we are not using the concept of conformal time explicitly in this section we go back to our old convention where derivatives with R are specified by primes.We will show that the constant curvature(de Sitter point) solution admits an inflationary scale-factor.We have not worked out the cosmology around the de Sitter point and the analogy with the inflationary scale-factor may turn out to be purely formal.Since out of the three equations Equations ( 4), ( 5) and ( 9), only two are independent, we can choose to work with equations Equations ( 4) and ( 9) for the sake of convenience.Using the standard solution ρ = ρ 0 a −3(1+ω) obtained from equation Equation ( 9) where ρ 0 and ω are constants, the constraint equation Equation ( 4), can be written as For (1/α)e αR gravity, this becomes Using the above equation let us now prove the existence of exact exponential bouncing solution and exact de-Sitter solution.

Exact Exponential Bouncing Solution
In this subsection we choose a scale-factor a(t) = e Ct 2 , where C is a positive constant.In such a case we get Using all these information and the modified constraint equation, Equation (68), we get If a(t) = e Ct 2 is an exact solution of (1/α)e αR gravity theory, then the above equation has to be satisfied for all values of t.This condition can be satisfied when αC = 1 48 , ω = − 4 3 .We see that for these values of the parameters α, C and ω, ρ 0 is given by 8κα , which is always positive.Exponential gravity can yield a bouncing universe with exponential scale-factor in presence of hydrodynamic matter which can have negative pressure but whose energy density is positive definite.The causal nature of the universe where the scale-factor is as given in this section is worked out in [35].

Exact de-Sitter Solution
A de-Sitter solution is a vacuum solution of constant positive curvature in GR.We use the same terminology and name our solution as de Sitter solution, as our solution in f (R) gravity is a constant positive curvature solution in presence of positive vacuum energy.We assume that the scale-factor of the universe as a(t) = e Ht , where H is a positive constant and consequently Using the above information in the modified constraint equation, Equation (68), and setting Therefore (1/α)e αR gravity has an exact de Sitter solution in which the constant value of the Ricci scalar is given by R = 2/α.In the Einstein frame this value of R correspond to φ = 2 3/(2κ), the point at which the potential V(φ) assumes its maximum.
As the de-Sitter point lies at the top of the Einstein frame potential one can intuitively conclude that the de Sitter solution is an unstable solution.In the rest of the section we will show that the de Sitter solution is indeed an unstable solution.To analyze the stability of this solution, we resort to a dynamical system analysis in the Jordan frame in terms of normalized, dimensionless dynamical variables as formulated in references [36][37][38][39].Let us define the dimensionless dynamical variables (71) Equation (68) then implies the following constraint between them which implies only two of them are independent.Note that for the de-Sitter solution Therefore for the constraint equation to be satisfied one must have H 2 = 1/(6α), which is consistent with what we previously obtained.Without loss of generality, we can take u 1 and u 2 as the two independent dynamical variables.With the help of the Equation ( 5), we can find the dynamical equations for u 1 , u 2 in terms of the dimensionless logarithmic time variable N ≡ ln a as: It is straightforward to check that the de-Sitter point given by the coordinates (0, 2) in the u 1 u 2 plane is a fixed point, i.e., at this point The first eigenvalue is positive and the second is negative, meaning the de-Sitter point is a saddle point in the space of isotropic vacuum solutions in exponential gravity.Figure 16 shows the flows of the solutions in the phase space.
After discussing the exact solutions in exponential gravity we discuss the extra solutions we get in exponential gravity if we allow exponential gravity to cross R = 0 value.In general we expect the theory to be similar to GR near R = 0 and the f (R) effect becomes effective for high values of R.But if we assume that f (R) gravity can be also used for small R values then we come across new kind of solutions presented in the next section.

New Solutions in Exponential Gravity
If we allow the cosmological model presented in this article to be valid for very small R then we get new class of solutions which, to our understanding, was never reported before in any discussions about f (R) gravity.The important feature of the solutions presented in this section is the generality of the discussion.Most of the results we present in this section is generally true for any f (R) theory which accommodates a cosmological bounce.We present the results for exponential gravity in the present paper.
For a cosmological bounce one requires H(0) = 0 and Ḣ(0) > 0 in the Jordan frame.In f (R) gravity we can give another intermediate condition, and this condition is on Ḧ(0).Suppose we specify Ḧ(0) = where is a real constant.Near the point where t = 0 we can then write approximately Ḣ(t) ≈ t + b where b is a positive, real constant.The approximation remains valid as long as Ḣ varies linearly near the bouncing point.In such a case we can integrate once more and write the expression of the Hubble parameter near the bounce point as the integration constant has been set to zero because the Hubble parameter vanishes at t = 0.The above equation shows that H can be zero at two separate time instants, given by The above statements are in general true for any f (R) theory which accommodates a cosmological bounce.The system behavior, in the particular case of exponential gravity, is shown in the Figure 17, which depicts an universe filled up radiation, and in Figure 18, which depicts the properties of matter less universe.In both the cases we see that the Hubble parameter reaches zero value at two time instants, as predicted from Equation (77).One of the points when the Hubble parameter is zero is placed at t = 0 where as the other time instant is approximately −2 × 10 5 , of the same order as predicted from Equation (77).If one wants to see how the scale-factor of the universe behaves during this period then one can look at Figures 19 and 20.From both the figures one can see that initially the universe was expanding and this expansion slowly stops and a brief period of contraction sets in, the period when the Hubble parameter turns negative, and then this contraction stops slowly and the universe again enters a period of expansion.It is to be noted that if we increase the value of , or Ḧ(0) in Jordan frame, then the temporal separation of the two roots in Equation (77) decreases.This implies that the initial expansion phase ends near to the point where later expansion phase starts.In the asymmetric bouncing solutions, as depicted in Figures 4 and 5, in Section 4 we have used non-zero values of Ḧ(0) in the Jordan frame, still we are getting perfect bounces and not oscillatory behavior as expected from our present analysis.The reason for this is related to the time period used to study the asymmetric bounces in Section 4 and the specific value of Ḧ(0) used there.From the conditions given for the asymmetric bounces one can easily show that the oscillatory behavior should have been evident if we presented the plot for a larger time period.In the time period of interest, only the bouncing region becomes prominent in Figures 4 and 5.
Before we conclude this section we want to show another important property of the new solutions shown.For flat FLRW solutions one can write the Ricci scalar as close to the bounce point.From the above expression one can easily verify that which shows that the Ricci scalar has to change sign in between the two temporal values when the Hubble parameter vanishes.As a result of this fact, the universe has to cross the point when R = 0, between the two expanding regimes.This result is a general result and is true for any f (R) theory which accommodates a cosmological bounce.If we want to work with large, positive values of R and treat f (R) theory as an effective gravitational theory then the new solutions become redundant as only a bounce exists, as analyzed in Section 4. On the other hand if one wants to keep the effect of f (R) bouncing cosmology for very small values of the Ricci scalar, then the new solutions will show up.

Conclusions
In this article we presented some early universe solutions coming out from exponential f (R) gravity.The primary reason for choosing exponential gravity is related to the stability of the theory.For positive cosmological constant one always gets f (R) > 0 and f (R) > 0 for all values of R. We presume the theory to be unstable when one moves away from bounce point as in the relevant time limit the Ricci scalar becomes negative and remains unbounded from below.This feature of the theory is well represented in the Einstein frame potential, as shown in Figure 1.We would like to interpret the early universe results, coming out from exponential f (R) theory, as an outcome from some effective theory of gravity which affects in the ultraviolet end.
In this paper we initially present the basics of cosmological dynamics in the two conformal frames, the Jordan frame where the original problem is posed and the Einstein frame which is used to calculate the dynamical development of the system.Our results related to bouncing cosmologies uses two conformal transformations.We show that it is difficult in f (R) theory to have stable bounces.In this regard exponential gravity is an interesting exception as it satisfies the stability conditions for all values as R.In Section 3 we present the bouncing condition in the Jordan frame and its corresponding condition in the Einstein frame.We want to specify here that in spatially flat FLRW spacetimes one in general does not get a simultaneous bounce in both the conformal frames when matter in Jordan frame satisfies ρ + P ≥ 0 [17].The conditions of bounce in Jordan frame, when written down in the Einstein frame, shows that the conditions depend upon the second time derivative of the Jordan frame Hubble parameter.Consequently, the bouncing intermediate conditions actually involves the values of all the relevant time derivatives of H in the Jordan frame.It is to be noted that unlike f (R) theories, GR-based cosmology does not require the specification of Ḧ at any instant of time.In this article we point out specifically how Ḧ can affect cosmological dynamics in higher derivative gravity theory.The evolution of scalar metric perturbations through bounce is presented in Section 5.The whole analysis is done in the Jordan frame and it is shown that the perturbations do not attain any singularity at the bounce point although there me be some modes which can become non-perturbative very near the bounce point.The reason for such non-perturbative evolution may be related to the behavior of curvature related terms near the bounce point.A thorough analysis of stability of metric perturbations in f (R) cosmology, near bounce, is not yet present we hope the theory will be formulated in the near future.We also briefly opine on the fate of the tensor-to-scalar ratio in our f (R) theory induced bounce.We do not expect the present theory to amplify the ratio compared to its value in the prebounce phase.
In the article we present the bouncing solutions in the Jordan frame.In the bouncing solution calculations we use the Einstein frame as an auxiliary frame where the main calculation is done and then we transport the solutions in the Jordan frame via a conformal transformation.The bouncing solutions presented in this paper involve both bounce in vacuum and bounce in the presence of hydrodynamic matter.The examples involve both symmetric and asymmetric bounces.The asymmetry of the bounces is related to finite values of Ḧ(0) in the Jordan frame.For symmetric bounces Ḧ(0) = 0.The solutions presented in Section 4 are computed numerically and they are well behaved in our time period of interest, −10 5 ≤ t ≤ 10 5 .Outside the time window the solutions can show other features.All the bounces concerned in Section 4 involves variation of positive values of the Ricci scalar in the Jordan frame.Treated f (R) theory as an effective theory of gravity, the exponential gravity bounces become more plausible at higher positive values of R. As exponential gravity becomes similar to GR with positive cosmological constant for low R values one may like to infer that GR effects becomes stronger as R nears zero.The exact turn over from f (R) to GR may involve new physics and is beyond the scope of this paper.
Another interesting property of exponential gravity, as discussed in the paper, is the existence of exact solutions.We present two such solutions in the article.One exact solution is a bouncing solution where the scale-factor of the universe is given by an exponential function of the square of Jordan frame time.This solution can be realized in an universe with matter having positive energy density and equation of state as −4/3.With an equation of state lesser than −1, the matter does not satisfy ρ + P ≥ 0 in the Jordan frame and consequently for such an universe one may expect simultaneous bounces in both the Jordan frame and the Einstein frame.The other exact solution is the exponential expansion solution with constant H at a de Sitter point.In vacuum, exponential gravity allows such a solution to exist.One can easily show that such a de-Sitter point exists because exponential gravity involves a positive cosmological constant.If the form of f (R) is suitably changed such that it does not contain any cosmological constant the de-Sitter point vanishes.Using the techniques of dynamical systems we have shown that a constant Hubble parameter solution at the de-Sitter point in exponential gravity is an unstable solution.In effect it is a saddle point solution in the phase space of suitable defined dimensionless phase space variables.We do not present the dynamical system approach in the other solutions in this paper because all the other solutions involve the values H = 0 and the phase space variables in our analysis always have the Hubble parameter in the denominator.We hope to construct a suitable dynamical systems approach to tackle bouncing problems in the near future.
We present a new solution in f (R) gravity theory in the penultimate section.The new solution in f (R) theories are allowed only if the theory is allowed to be unmodified in the low R regime.As our theory is stable it can safely be extended to the low R regime.The only cost one has to pay to attain these new solutions is that one has to reject the point of view that f (R) gravity is an ultraviolet modification of GR effects.The new results are related to non-zero values of Ḧ(0) in the Jordan frame when the other bouncing conditions hold.In such a case one can a have a solution which represent decelerated expansion of the universe in the past.At some point in the past the decelerated expansion comes to a halt momentarily and contraction of the universe starts.This contraction does not lead to a spacetime singularity.In time this contraction slows down and the universe comes to a static configuration momentarily after which again the universe starts to expand.This solution is practically not a bouncing solution, although one can get this solution with an extra intermediate condition on top of the bouncing conditions at t = 0.The new result which we obtain in this paper is a general result in f (R) gravity which accommodates a bounce.We explicitly show the nature of the solutions in exponential gravity.
Figure 1. Figure showing the nature of the scalar field potential in the Einstein frame where f (R) = (1/α) exp(αR).In the above figure the field and its potential both are expressed in Planck units (where the Planck mass is set as unity) in which α = 10 12 .

Figure 6 .Figure 7 .
Figure 6.Time evolution of the comoving Hubble radius for the time-symmetric bouncing solution of Figure 2 in radiation background.

Figures 8
Figures 8 and 9 correspond to the first of the above conditions and Figures 10 and 11 correspond to the second.As is clearly seen from the plots, for both of these conditions produce perturbative evolution of the even Bardeen potentials through the bounce.The various other cases of perturbation evolution are shown in Figures 12-15.

du 1 dN 1 2
= du 2 dN = 0.The eigenvalues of the Jacobian matrix at this point are

Figure 16 .
Figure16.Phase space plot in the u 1 u 2 plane of isotropic vacuum solutions in exponential gravity, where u 1 is along the x-axis and u 2 is along the y-axis.The de-Sitter solution is given by the point (0, 2), which is represented by the dot at the center of the figure.The arrows show the direction of the solution flow.The phase flows clearly show that the de-Sitter solution is a saddle fixed point.