Inflation with scalar-tensor theory of gravity

The latest released data from Planck in 2018, put up tighter constraints on inflationary parameters. In the present article, the in-built symmetry of the non-minimally coupled scalar-tensor theory of gravity is used to fix the coupling parameter, the functional Brans-Dicke parameter, and the potential of the theory. It is found that all the three different power-law potentials and one exponential, pass these constraints comfortably, and also gracefully exit from inflation.


Introduction
The standard (FLRW) model of cosmology based on the basic assumption of homogeneity and isotropy, known as the 'cosmological principle', has successfully been able to explain several very important issues in connection with the evolution of the universe. First of all, it predicts the observed expansion of the universe being supported by the Hubble's law. It also postulates the existence of cosmic microwave background radiation (CMBR), formed since recombination when the electrons combined to form atoms, allowing photons to free stream, with extreme precession, being verified by Penzias and Wilson for the first time. It also predicts with absolute precession the abundance of the light atomic nuclei ( H ∼ 10 −9 , by mass and not by number) observed in the present universe. Finally, assuming the presence of the seeds of perturbation in the early universe, it can explain the observed present structure of the universe. Despite such tremendous success, the model inevitably suffers from a plethora of pathologies. The problems at a glance are the following. 1. 'The singularity problem': Extrapolating the FLRW solutions back in time one encounters an unavoidable singularity, since all the physical parameters viz. the energy density (ρ), the thermodynamic pressure (p), the Ricci scalar (R ), the Kretschmann scalar (R αβγδ R αβγδ ) etc. diverge. 2a. 'The flatness problem': The model does not provide any explanation to the observed value of the density parameter Ω ≈ 1 , which depicts that the universe is spatially flat. 2b. 'The horizon problem': It also can not provide any reason to the observed tremendous isotropy of the CMBR being split in 1.4 × 10 4 patches of the sky, that were never causally connected before emission of the CMBR. 2c. 'The structure formation problem': It does not also provide any clue to the seeds of perturbation responsible for the structure formation. 3. 'The dark energy problem': Finally, the standard FLRW model does not fit the redshift versus luminositydistance curve plotted in view of the observed SN1a (Supernova type a) data.
In connection with the first problem, viz. the so called 'Big-Bang singularity' and also to understand the underlying physics of 'Black-Hole' being associated with Schwarzschild singularity, it has been realized long ago that 'General Theory of Relativity' (GTR) must have to be replaced by a quantum theory of gravity when and where gravity is strong enough. However, GTR is not renormalizable and a renormalized theory requires to include higher-order curvature invariant terms in the gravitational action. Despite serious and intense research over several decades and formulation of new high energy physical theories like superstring and supergravity theories, a viable quantum theory of gravity is still far from being realized. In connection with last problem, a host of research is in progress over last two decades. It has been realized that to fit the observed redshift versus luminosity-distance curve, it is either required to take into account some form of exotic matter in addition to the barotropic fluid (ordinary plus the cold dark matter) which violates the strong energy condition (ρ + p ≥ 0, ρ + 3p ≥ 0 ), and is dubbed as 'dark energy' (since it interacts none other than with the gravitational field) or to modify the theory of gravity by including additional curvature scalars in the Einstein-Hilbert action, known as 'the modified theory of gravity'. It has been observed that both the possibilities lead to present accelerated expansion of the universe. The problem is thus rephrased as: why the universe undergoes an accelerated expansion at present? The pathology 2, in connection with the flatness, horizon and structure formation problems has however been solved under the hypothesis called 'Inflation', which is our present concern.
Under the purview of cosmological principle, i.e. taking into account Robertson-Walker line element, the co-moving distance (the present-day proper distance) traversed by light between cosmic time t 1 and t 2 in an expanding universe may be expressed as, d 0 (t i , t f ) = a 0 t f ti dt a(t) , where a(t) is called the scale factor. The co-moving size of the particle horizon at the last-scattering surface of CMBR (a f = a lss ) corresponds to d 0 ∼ 100 Mpc, or approximately 1 0 (one degree) on the CMB sky today. In the decelerated radiation dominated era of the standard model of cosmology (FLRW model), for which a ∝ √ t the integrand, (ȧ) −1 ∼ 2 √ t decreases towards the past, and there exists a finite co-moving distance traversed by light since the Big Bang (a i → 0 ), called the particle horizon. The hypothesis of inflation [1,2] postulates a period of accelerated expansion, a > 0 , in the very early universe, prior to the radiation-dominated era, administering certain initial conditions [3,4,5,6,7,8,9]. During a period of inflation e.g. a de-Sitter universe (a ∝ e Λt ) driven by a cosmological constant (say), (ȧ) −1 ∼ (Λe Λt ) −1 increases towards the past, and hence the integral diverges as (a i → 0 ). This allows an arbitrarily large causal horizon dependent only upon the duration of the accelerated expansion. Assuming that the universe inflates with a finite Hubble rate H i , (instead of a constant exponent Λ ) ending with H f < H i , we is measured in terms of the logarithmic expansion (or 'e-folds'), and describes the duration of inflation. It has been found that a 40 − 60 e-folds of inflation can encompass our entire observable universe today, and thus solves the horizon and the flatness problem discussed earlier. In some situations, e-fold may range between 25 ≤ N ≤ 70 , depending on the model under consideration.
A false vacuum state can drive an exponential expansion, corresponding to a de-Sitter space-time with a constant Hubble rate on spatially-flat hypersurfaces. However, a graceful exit from such exponential expansion requires a phase transition to the true vacuum state. A second-order phase transition [10,11], under the slow roll condition of the scalar field (that can also drive the inflation instead of the cosmological constant), potentially leads to a smooth classical exit from the vacuum-dominated phase. Further, the quantum fluctuations of the scalar field, which essentially are the origin of the structures seen in the universe today, provides a source of almost scale-invariant density fluctuations [12,13,14,15,16], as detected in the CMBR. Accelerated expansion and primordial perturbations can also be produced in some modified theories of gravity (e.g., [1,17] and also a host of models presently available in the literature), which introduce additional non-minimally coupled degrees of freedom. Such inflationary models are conveniently studied by transforming variables to the so-called Einstein frame, in which Einsteins equations apply with minimally coupled scalar fields [18,19], which we shall deal with, in the present manuscript.
Non-minimal coupling with the scalar field φ is unavoidable in a quantum theory, since such coupling is generated by quantum corrections, even if it is primarily absent in the classical action. Particularly, it is required by the renormalization properties of the theory in curved space-time background. Recently, in view of a general conserved current, obtained under suitable manipulation of the field equations [20,21,22,23], a non-minimally coupled scalar-tensor theory of gravity has been studied extensively in connection with the cosmological evolution, starting from the very early stage (Inflationary regime) to the late-stage (presently accelerated matter-dominated era) via a radiation dominated era [24]. It has been found that such a theory admits a viable inflationary regime, since the inflationary parameters viz. the scalar-tensor ratio (r ), and the spectral index n s lie well within the limits of the constraints imposed by Planck's data, released in 2014 [25] and in 2016 [26]. Further the model passes through a Friedmann-like radiation era (a ∝ √ t, q = 1 ,) and also an early stage of long Friedmann-like decelerating matter dominated era (a ∝ t a red-shift, z ≈ 0.75 , which is very much at par with recent observations. Further, the present numerical values of the cosmological parameters obtained in the process are also quite absorbing, since the age of the universe (13.86 < t 0 < 14.26 ) Gyr, the present value of the Hubble parameter (69.24 < H 0 < 69.96 ) Km.s −1 Mpc −1 , so that 0.991 < H 0 t 0 < 1.01 fit with the observation with appreciable precision. Numerical analysis also reveals that the state finder {r, s} = {1, 0} , which establishes the correspondence of the present model with the standard Λ CDM universe. Last but not the least important outcome is: considering the CMBR temperature at decoupling (z ∼ 1080 ) to be 3000 , required for recombination, it's present value is found to be 2.7255 , which again fits the observation with extremely high precision. Thus, non-minimally coupled scalar-tensor theory of gravity appears to serve as a reasonably fair candidate for describing the evolution history of our observable universe, beyond quantum domain.
In the mean time new Planck's data is released [27,28], which imposed even tighter constraints on the inflationary parameters. In this manuscript, we therefore pose if the theory [24] admits these new constraints. However, earlier we considered a particular form of coupling parameter along with the potential in the form V (φ) = V 0 φ 4 − Bφ 2 , where, V 0 and B are constants [24]. Here instead, we choose different forms of the coupling parameters and also different potentials to study the inflationary regime. In the following section 2, we describe the model, write down the field equations, find the parameters involved in the theory in view of a general conserved current. We also present the scalar-tensor equivalent form of the action in Einstein's frame to find the inflationary parameters. In section 3, we choose different forms of the coupling parameters and associated potentials to test the viability of the model in view of the latest released data from Planck [27,28]. We conclude in section 4.
2 The model, Conserved current, scalar-tensor equivalence and inflationary parameters: We start with the non-minimally coupled scalar-tensor theory of gravity, for which the action is expressed in the form, where, L m is the matter Lagrangian density, f (φ) is the coupling parameter, while, ω(φ) is the variable Brans-Dicke parameter. The field equations are, where prime denotes derivative with respect to φ, and denotes D'Alembertian, such that, The model involves three parameters viz. the coupling parameter f (φ), the Brans-Dicke parameter ω(φ) and the potential V (φ). It is customary to choose these parameters by hand in order to study the evolution of the universe. However, we have proposed a unique technique to relate the parameters in such a manner, that choosing one of these may fix the rest [20,21,22,23,24]. This follows in view of a general conserved current which is admissible by the above pair of field equations, briefly enunciated below.
The trace of the field equation (3) reads as, Now eliminating the scalar curvature between equations (4) and (5), one obtains, which may then be expressed as, and finally as, Thus there exists a conserved current J µ , where, for trace-less matter field (T µ µ = T = 0 ), provided To study cosmological consequence of such a conserved current, let us turn our attention to the minisuperspace model (1), in which the conserved current (9), reads as in traceless vacuum dominated and also in radiation dominated eras. In the above, C 1 is the integration constant. Note that, fixing the form of the coupling parameter f (φ), the potential V (φ) is fixed in view of (10), once and forever. Further, we use a relation [24] 3f ′2 + 2ωf where, ω 0 is a constant, to fix the Brans-Dicke parameter as well. As a result, we obtain the relation C being yet another constant. In the process, all the coupling parameters f (φ), ω(φ) and the potential V (φ) may be fixed a-priori. Note that the above choice (12) finally leads to the conserved current associated with the canonical momenta conjugate to the scalar field φ, in the absence of a variable Brans-Dicke parameter, i.e. for . We shall work, in the typical unit, M 2 p 2 = c = 1 , and consider different forms of f (φ), that fixes the Brans-Dicke parameter ω(φ) as well as the potential V (φ). In view of these known functional forms of the parameters of the theory, we focus our attention to study inflation, which must have occurred in the very early vacuum dominated universe. We relax the symmetry by adding an useful term in the potential V (φ), so that one of the terms act just as a constant in the effective potential. This ensures a constant value of the potential as the scalar field dies out, and this constant acts as an effective cosmological constant (Λ e ). In the process, the number of parameters increases to three (ω 0 , V 0 , V 1 ), which is essential to administer good fit with observation.

Scalar-tensor equivalence and inflationary parameters:
As mentioned, it is convenient and hence customary to study inflationary evolution in the Einstein's frame under suitable transformation of variables, where possible. Therefore, in order to study inflation, we consider very early vacuum dominated (p = 0 = ρ, for which trace of the matter field identically vanishes and symmetry holds) era, and express the action (2) in the form, where, K(φ) = 2 ω(φ) φ . The above action (14) may be translated to the Einstein's frame under the conformal transformation (g Eµν = f (φ)g µν ) to take the form [29], where, the subscript 'E ' stands for Einstein's frame. The effective potential (V E ) and the field (σ ) in the Einstein's frame may be found from the following expressions, In view of the action (15), it is also possible to cast the field equations, viz. the Klein-Gordon and the ( 0 0 ) equations of Einstein as, where, H =ȧ E aE denotes the expansion rate, commonly known as the Hubble parameter. The slow-roll parameters and the number of e-foldings, then admit the following forms, where, t i , t f stand for the initiation time and the end time, while φ b , φ e stand for the values of the scalar field at the beginning and at the end of inflation respectively. Comparing expression for the primordial curvature perturbation on super-Hubble scales produced by single-field inflation (P ζ (k)) with the primordial gravitational wave power spectrum (P t (k)), one obtains the tensor-to-scalar ratio for single-field slow-roll inflation r = Pt(k) P ζ (k) = 16ǫ , while, the scalar tilt, conventionally defined as n s − 1 may be expressed as n s − 1 = −6ǫ + 2η , or equivalently n s = 1 − 6ǫ + 2η , dubbed as scalar spectral index. According to the latest released results, the scalar to tensor ratio r ≤ 0.16 (TT,TE,EE+lowEB+lensing), while r ≤ 0.07 (TT,TE,EE+lowE+lensing+BK14+BAO) [27,28]. Further, combination of all the data (TT+lowE, EE+lowE, TE+lowE, TT,TE,EE+lowE, TT,TE,EE+lowE+lensing) constrain the scalar spectral index to 0.9569 ≤ n s ≤ 0.9815 [27,28]. It is useful to emphasize that under the present choice of unit M 2 p 2 = c = 1 (which although appears to be a bit unusual but doesn't cause any harm), φ controls the cosmological evolution in the manner φ > 1 corresponds to the inflationary stage, φ ∼ 1 describes the end of inflation while φ < 1 is the low energy regime which triggers matter dominated era.

Inflation with power law and exponential potentials:
In the non-minimal theory, the flat section of the potential V (φ) responsible for slow-roll is usually distorted. However, flat potential is still obtainable if the Einstein's frame potential V E is asymptotically constant [30,31]. Note that, the symmetry explored in section 2, makes the potential (V E ) in the Einstein's frame to be constant, once and forever. Thus, we need to relax the symmetry as required by the condition (10), by taking into account additional term in the potential, viz. a constant term V 0 , or even a functional form, to ensure that the effective potential in the Einstein's frame (V E ) is asymptotically constant. At the end of inflation, the universe becomes cool due to sudden large (exponential, in the present case) expansion. Therefore, in order that the structure we live in are formed, the universe must be reheated and take the state of a hot thick soup of plasma (the so called hot Big-Bang). This phenomena is possible if at the end of inflation, the scalar field starts oscillating rapidly on the Hubble time scale, about the minimum of the potential. In the process, particles are created under standard quantum field theoretic (in curved space-time) approach, which results in the re-heating of the universe. The universe then eventually transits to the radiation dominated era. At that epoch, the additional term may be absorbed in the potential if it is a constant term (V 0 ), or may even be neglected in case it is a function (since φ goes below the Planck's mass), without any loss of generality, to reassure symmetry. The symmetry leads to the first integral of certain combination of the field equations, which helps in solving the field equations leading to a Friedmann-like radiation dominated era, as shown earlier [24]. However, in the present manuscript, we only concentrate upon inflationary regime and of-course study possibility of graceful exit from inflation. In the following subsection, we shall study different power law potentials, while in the next we shall deal with exponential potential. We consider de-Sitter solution in the form a ∝ e Ht , where the Hubble parameter H is slowly varying during inflation. We repeat that according to our current choice of units ( , which although is uncommon, but doesn't create any problem whatsoever, the value of f (φ) at the end of inflation must be a little greater than 1 .

Power law potential
We shall take into account three different values of n, viz. n = 1, 3 2 and 2 , in the following three sub-subsections. For each value of n we shall study different cases taking into account different additive terms. In the first place however, we shall consider an additive constant V 0 in all the three cases, viz.
corresponding to which, one can now find the expression for the Brans-Dicke parameter ω(φ), the potential V E (σ) in the Einstein's frame, the expression for dσ dφ , and the slow-roll parameters ǫ, η along with the number of e-foldings N , in view of the equations (12), (16) and (18) respectively as, Case − 1 : The effect of the constant term V 0 is now clearly noticeable, since when φ is large, second term in the Einstein's frame potential (V E ) becomes insignificantly small, for n ≥ 1 , and it almost becomes (non-zero) constant, assuring slow-roll. On the contrary, if V 0 is set to vanish from the very beginning, the Einstein's frame potential V E = V 1 would remain flat always, and the universe would have been ever-inflating.
We shall also consider a functional additive term in the potential for all the cases under consideration, such the the potential reads as In view of the above potential (21) it is possible to find the expression for the Brans-Dicke parameter ω(φ), the potential V E (σ) in the Einstein's frame, the expression for dσ dφ , and the slow-roll parameters ǫ, η along with the number of e-foldings N , in view of the equations (12), (16) and (18) respectively as, Case − 2 : In the following sub-subsections, we shall take three different values of n, as already mentioned, and present the data set in tabular form along with appropriate plots, to demonstrate the behaviour of the slow-roll parameters in comparison with the latest data set released by Planck [27,28]. Different additive terms, as indicated, will be considered in each subcase separately. In the subsection (3.1.3), we shall consider an additional case with a pair of additive terms in the form of a whole square.
Case-1: Under the choice n = 1 , the potential (19) takes the form V (φ) = V 0 + V 1 φ 2 , and thus the parameters of the theory under consideration (20) read as, In view of the above forms of the slow roll parameters (23), we present table-1 and table-2, underneath, corresponding to two different values of the parameter V 1 > 0 . The wonderful fit with the latest data sets released by Planck [27,28] is appreciable particularly because 0.968 ≈ n s < 0.982 , while r < 0.0278 . Further, the number of e-fold (36 ≤ N ≤ 62 ) is sufficient to alleviate the horizon and flatness problems. Figure 1 and figure 2 are the two plots r versus n s and r versus ω 0 respectively, presented for visualization. For example, the figures clearly depict that the plot which represents data sets corresponding to table 2, appears to be even better.   Table 2: One very interesting feature is that the above data sets remain unaltered even if the sign of V 0 and V 1 are interchanged. Note that, second derivative of the potential has to be positive, since it represents effective mass of the scalar field. In view of the forms of the potentials V (φ) and V E (σ) presented in (19) and (20) the effective mass of the scalar fields φ and σ respectively are, In our data set, we keep V 1 > 0 , since φ is the scalar field under consideration, while translation to σ only amounts to handling the situation with considerable ease. However, as a matter of taste if one favours Einstein's frame over Jordan's frame, it is possible to revert the sign and keep V 0 > 0 , without changing the data set. As mentioned, at the end of inflation, the scalar field must oscillate rapidly so that particles are produced and the universe turns to the phase of: a hot thick soup of plasma, commonly called the 'hot big-bang'. This phenomena is dubbed as graceful exit, which is required for the structure formation together with the formation of CMB. We therefore proceed to check if the present model admits graceful exit from inflation. Here, V E = V 1 + V0 φ 2 , and so one can express (17) as, At large value of the scalar field, which in the present unit φ > 1 , we obtained slow-roll. However, as the scalar field falls below the Planck's mass M p , then the Hubble rate H also decreases, and once it falls below the effective mass V 1 i.e. H ≪ V 1 , then the above equation may be approximated to, Where,σ = ω0 φφ , in view of (23). Thus, finally we get, which is an oscillatory solution, and the field then oscillates many times over a Hubble time. This coherent oscillating field corresponds to a condensate of non-relativistic massive (inflaton) particles, which ensures graceful exit from the inflationary regime, driving a matter-dominated era at the end of inflation. There is a long standing debate regarding the physical frame. It appears that most of the people favour Einstein's frame over Jordan's frame (we have briefly discussed the issue in conclusion). In this regard, it is important to mention that since in view of (23) σ = ω 0 ln φ, therefore σ executes oscillatory behaviour as well.
Case-2: Under the same situation f (φ) = φ, let us now consider, V (φ) = V 1 φ 2 +V 0 φ 4 , where instead of a constant term, we have added a quartic term in the potential. The expression for the Brans-Dicke parameter (ω(φ)), the potential (V E ) in the Einstein's frame, dσ dφ , the slow-roll parameters ǫ, η and the number of e-foldings N , may then be found in view of the equation (22), respectively as, Although, the potential V E does not appear to attend a flat section, the smallness of the value of η confirms that there indeed exists a flat section, admitting slow-roll. In fact, in the Einstein's frame (15), this is just the case of a standard inflation field theory with quadratic potential. Followings tables 3 and 4, for V 1 > 0 , together with the associated plots n s versus r and r versus ω 0 here again depict appreciably good fit with the recent released Planck's data set, particularly because 0.97 ≤ n s ≤ 0.98 while r ≤ 0.098 . The figure 3 depicts that data of  Table 3:  Table 4: graceful exit. Here V E = (V 1 + V 0 φ 2 ), and so from (17) one obtains, As the Hubble rate H ≪ V 1 , the above equation can be approximated as, The oscillatory behaviour of the scalar field clearly ensures graceful exit from inflationary regime, as already discussed, and in view of (28) σ also executes oscillatory behaviour.

Case-1:
Under the choice n = 3 2 , f (φ) = φ 3 2 , and the potential (19) takes the cubic form, V (φ) = V 0 + V 1 φ 3 . Thus the expressions for the parameters of the theory under consideration along with the slow-roll parameters (20) are, As before, we present two sets of data in tables 5 and 6, for two different values of V 1 > 0 . Plots 5 and 6 depict the variations of the spectral index n s with the scalar-tensor ratio r and the scalar-tensor ratio r with the Brans-Dicke parameter ω 0 respectively. Here again we observe that r ≤ 0.0278 , and 0.96 < n s < 0.9815 , which are very much within the stipulated observational range [27,28], while number of e-folding N is sufficient to remove the flatness and the horizon problems.      To check the behaviour of the scalar field, we proceed as before, to find,  under approximation, as the Hubble rate H ≪ V 1 , and using the relationσ = ω 0φ φ 3 2 , in view of (31). The solution reads as, In the above the hash tag (#n) denotes nth argument of a pure function, and c 1 is a constant. Although, the solution is not obtainable in closed form, rather is a complicated inverse elliptic function, nevertheless its oscillatory behaviour is quite apparent, and σ = −2 ω0 √ φ also oscillates as well.
Case-2: Cubic potentials with additive term have important consequence. For example, a potential in the form V = 1 2 mω 2 x 2 − 1 3 bx 3 can be used to model decay of metastable states [32], and it also describes the global flow [33]. Further, the tunnelling rate in real time in the semiclassical limit may be found for arbitrary energy levels, while it's ground state agrees well with the result found by the instanton method [34]. It is therefore worth to continue the present study in view of such an additive form in the cubic potential.
Under the choice n = 3 2 , f (φ) = φ 3 2 , and taking the potential as cubic form added with a quadratic term , i,e, V (φ) = V 1 φ 3 + V 0 φ 2 , the expressions for the parameters of the theory under consideration along with the slow-roll parameters (22) are, We present two sets of data in tables 7 and 8 underneath, for two different values of V 1 > 0 . Plots 7 and 8 depict the variations of the spectral index n s with the scalar-tensor ratio r and the scalar-tensor ratio r with the Brans-Dicke parameter ω 0 respectively. Here again we observe that r ≤ 0.062 , and 0.973 < n s < 0.983 , which are again in excellent agreement of Planck's data [27,28], while the number of e-folding N is also sufficient to remove the flatness and the horizon problems. It is interesting to note that the variation n s with r for the two sets of data almost overlap in figure-7. In order to study the behaviour of the scalar field at the end of inflation,    we start with the Einstein's frame potential as before, V E = (V 1 + V0 φ ), and express the field equation (17) as,  As the Hubble rate falls, and H ≪ V 1 , the above equation may be approximated to,σ 2 = 2i 2 (V 1 + V0 φ ), which in terms of the scalar field φ reads as, which may be solved to find which unfortunately is not oscillatory. Perhaps, due to the asymmetry of the potential, the oscillatory behaviour of the scalar field with an additive quadratic term is not exhibited.
Case-1: Under the choice n = 2 , f (φ) = φ 2 , and the potential (19) is now V (φ) = V 0 + V 1 φ 4 . Therefore the Brans-Dicke parameter, the Einstein's frame potential together with the slow-roll parameters (20) take the following forms, It is quite transparent that for large value of the scalar field φ, a flat Einstein's frame potential is realizable here too. As before, we take two sets of data corresponding to two different values of V 1 > 0 , and tabulate the parametric values in tables 9 and 10. One can see that the scalar tensor ratio r ≤ 0.0202 , and the spectral index lies between 0.9645 ≤ n s ≤ 0.9809 , which are in excellent agreement with Planck's data [27,28]. Further number of e-folding N is also sufficient to alleviate the flatness and the horizon problems. The n s versus r and r versus ω 0 plots are presented in figures 9 and 10, as well. Equation (17) now reads as, which, as H falls below V 1 , i.e. H ≪ V 1 , may be approximated to, 1 2σ 2 + V0 φ 4 + V 1 = 0 . In terms of the scalar field φ it is expressed asφ 2 + 2 ω 2 0 V 1 φ 4 + V 0 = 0 , since,σ = ω 0φ φ 2 , in view of (38). The solution is,  Table 9:  Table 10:   where JacobiSN is a meromorphic function in both arguments, which certain special arguments may automatically be evaluated to exact values. In any case, under numerical simulation the above solution is found to exhibit oscillatory behaviour of the scalar field φ. It is also clear that σ = − ω0 φ oscillates as well, and the universe transits from inflationary regime to the matter dominated era.

Case-2:
Here, for f (φ) = φ 2 , we consider the potential in the form, V = V 1 φ 4 + V 0 φ 2 , i.e. instead of a constant additive term, we consider V 0 φ 2 in addition. This case was earlier studied in [24]. However, as already mentioned, in the years, Planck's data puts up tighter constraints on inflationary parameters, and so it is quite reasonable to check if this form of potential passes the said constraints [27,28]. One can now find the expression for the Brans-Dicke parameter ω(φ), the potential V E in the Einstein's frame, dσ dφ , the slow-roll parameters ǫ, η and the number of e-folding N , in view of the equations (22), respectively as, Note that the Einstein's frame potential now takes the same form as in case-1 for n = 1 , and a flat section of the potential is still realizable at large value of the scalar field φ. We present two tables 11 and 12, as before for different values of V 1 > 0 . The scalar to tensor ratio r ≤ 0.0636 and the spectral index 0.9715 ≤ n s ≤ 0.9831 lie very much within the Planck's data, while the number of e-folding N is again sufficient to alleviate the horizon and flatness problems. The plots (figures 11 and 12) represent n s versus r and r versus ω 0 respectively. In view of the plots, the data for table 12, here appears to be even better.
To check if the scalar field executes oscillatory behaviour at the end of inflation, we note that here  Table 11:  Table 12:    So in view of equation (17) one obtains, As, H falls below V 1 , and H ≪ V 1 , the above equation can be approximated to, It is also possible to solve for φ and express it in the following form, It is now quite apparent that φ executes oscillatory behaviour and therefore graceful exit from inflation is realizable. Since in view of (41) σ = − ω0 φ , therefore σ also executes oscillatory behaviour.

Case-3:
We consider yet another case for n = 2 , i.e. taking f (φ) = φ 2 , with the potential V (φ) being represented by two additional terms apart from φ 4 , as it should be, to make it a perfect square: As before, one can now find the expression for the Brans-Dicke parameter ω(φ), the potential V E in the Einstein's frame, dσ dφ , the slow-roll parameters ǫ, η and the number of e-foldings N , in view of the equations (12), (16) and (18) respectively as, One can clearly see that the flat section of the potential is still attainable for large value of the scalar field φ. Table 13 and Table 14 depict that the scalar to tensor ratio r < 0.1 ), is quite reasonable, while the spectral index 0.9752 ≤ n s ≤ 0.981 fits perfectly with Planck's data [27,28]. Figures 13, and 14 represent n s versus r and r versus ω 0 plots respectively. Interestingly, two r versus n s plots (figure-13) corresponding to the two sets of data (Table-13 and Table- Table 14: The scalar field executes oscillatory behaviour here to, as we demonstrate below. Here, , and so from (17) we find,  As H falls below V 1 , and H ≪ V 1 , the above equation can be approximated as, where,σ = ω 0φ φ 2 . Therefore finally we obtain, Clearly, φ executes oscillatory behaviour, and graceful exit from inflation may be realized hereto. Here again since in view of (45) σ = − ω0 φ , therefore σ executes oscillatory behaviour, as well.

Exponential potential:
Finally, we consider an exponential form of the potential with: f (φ) = e λφ 2 , with V (φ) = V 1 e λφ + V 0 . It is possible to find the expression for the Brans-Dicke parameter ω(φ), the potential V E in the Einstein's frame, dσ dφ , the slow-roll parameters ǫ, η and the number of e-folding N , in view of the equations (12), (16) and (18) respectively as, ; We present our results in the following table 15, under the only choice of the parameter λ = −1 . The data shows good agreement r < 0.06 , and 0.9612 ≤ n s ≤ 0.9836 with Planck's data [27,28]. Figure 15, represents a plot for r versus n s . The scalar field executes oscillatory behaviour here to, as we demonstrate below.
Scalar-tensor theories of gravity are generalizations of the Brans-Dicke theory, in which the coupling parameter is a function of the scalar field, i.e. ω BD = ω(φ), and therefore is a variable. The requirement for such generalization of Brans-Dicke theory generated from the tight constraints on ω BD established by the solar system experiments [35]. There exists various classification of scalar-tensor theory of gravity [36]. In the present manuscript we have considered standard non-minimal coupling, where the coupling parameter f (φ) is arbitrary. It has been noticed earlier that such a theory has an in-built symmetry being associated with a conserved current for trace-free fields, such as vacuum and radiation dominated eras for barotropic fluids. In view of such a symmetry it is possible to fix all the variables of the theory, including the potential function, fixing the form of one of those. In this manuscript, we have chosen different forms of the coupling parameter f (φ), which fixed ω(φ) and V (φ), to study the cosmological evolution of the very early universe in the context of inflation. Inflation is a quantum mechanical phenomenon, and has occurred around Planck's era. However, it has been argued that since the radiative corrections to the potential are negligible, hence the inflationary parameters can be computed using the classical Lagrangian [37]. This argument leads in general, to calculate inflationary parameters in view of the classical Lagrangian, which we have done in the present manuscript. The so called unification programmes, which essentially claim to unify early inflationary regime with late-time cosmic acceleration have no credentials, since none of the models passes through a well behaved radiation and early matter dominated era. However, a history of cosmic evolution starting from inflationary regime, followed by a Friedmann-like radiation (a ∝ √ t) and early matter dominated eras (a ∝ t 2 3 ), that finally ends up to a late-time accelerated unverse (z = 0.75 ), has already been explored in view of the present [24]. In this connection, the present model makes a reasonably viable attempt to unify early inflation with late-time cosmic acceleration. Nevertheless earlier, only a single form of the coupling parameter f (φ) together with a particular form of V (φ) had been treated. Here, we have extended our work considering at least three power potentials together with an exponential potential. We find that quadratic, cubic, quartic and exponential potential pass the tighter constraints on inflationary parameters released by latest Planck's data [27,28] comfortably. Further, all these potentials admit graceful exit from inflation, except one case of cubic potential associated with a square potential, for which unfortunately the scalar field does not show up oscillatory behaviour at the end of inflation.
For the purpose of the present analysis, we have translated the non-minimally coupled Jordan's frame action to the Einstein's frame, under conformal transformation. It is therefore worth to make certain comments in this regard. There is an age old debate regarding physical equivalence between the two: Jordan's and Einstein's frames, which are related under conformal transformation. Now, indeed if the two formulations are not equivalent, the problem arises in selecting the physically preferred frame. It emerges from the work of several authors, in different contexts on Kaluza-Klein and Brans-Dicke theories, that the formulations of a scalar-tensor theory in the two conformal frames are physically inequivalent [38,39,40,41,42,43,44,45]. Also the Jordan frame formulation of a scalar-tensor theory is not viable because the energy density of the gravitational scalar field present in the theory is not bounded from below, which amounts to the violation of the weak energy condition [46]). The system therefore is unstable and decays toward a lower and lower energy state ad infinitum [43,44,45]. Although, a quantum system may have states with negative energy density [46,47,48], such feature is not acceptable for a viable classical theory of gravity. In fact, a classical theory must have a ground state that is stable against small perturbations. The violation of the weak energy condition by scalar-tensor theories formulated in the Jordan conformal frame makes them unviable descriptions of classical gravity, while the Einstein frame formulation of scalar-tensor theories is free from such problem. However, in the Einstein frame also there is a violation of the equivalence principle due to the anomalous coupling of the scalar field to ordinary matter. Nevertheless, this violation is small and compatible with the available tests of the equivalence principle [49]. Further, Einstein's frame is indeed regarded as an important low energy manifestation of compactified theories [49,50,51,52,53,54]. However, in search of Noether symmetries of F (R) theory of gravity, the two frames have been found to be physically equivalent [55]. So although the debate persists, but somehow it is quite relevant to consider Einstein's frame to be the physical frame. Therefore, In view of the above discussions, it is justified to study the physics associated with non-minimally coupled scalar-tensor theory of gravity, after translating it to the Einstein's frame, as we have done in the present article.