the global evolution of the observable universe can be separated on four epochs : inflation , radiation domination era , matter domination era and dark energy domination era .
all these epochs can be described by general relativity models with minimally coupled scalar fields , because the assumption that the hubble parameter is a monotonically decreasing function does not contradict with the observation data . in such type of models the initial period of the universe evolution with energies above the planck scale
should be described by quantum gravity because the classical evolution includes the initial singularity .
important question of theoretical cosmology is whether the entire universe evolution can remain classical and has no singularity .
it is possible to avoid this singularity considering modified gravity .
bouncing universe smoothly transits from a period of contraction to a period of expansion and its evolution remains classical .
models of bouncing universes attract a lot of attention  @xcite .
we can mention , in particular , bouncing models with galileon fields  @xcite , @xmath0 gravity models  @xcite , non - local gravity models  @xcite and models with non - minimal coupling  @xcite .
the goal of this paper is to explore new mathematical features of cosmological models with non - minimally coupled scalar fields that admit bounce solutions .
we consider models with the spatially flat friedmann
lematre  robertson  walker ( flrw ) metric . at the bounce point
the period of universe contraction changes to a period of universe expansion .
thereby , a bounce point is characterized by two condition : at this point the hubble parameter @xmath1 is equal to zero and its cosmic time derivative is positive .
note that the general relativity models with minimally coupled standard scalar fields does not have such a type of solutions , because in these models the hubble parameter is a monotonically decreasing function .
the simplest way to get a bounce solution is to consider a model with a phantom scalar field . note that similar models is used to describe dark energy with the state parameter less then minus one  @xcite .
to get non - monotonic behaviour of the hubble parameter quintom models have been developed  @xcite . in phantom field models
the null energy condition is violate and instability problems arise  @xcite .
at the same time it is easy to get increasing or non - monotonic behavior of the hubble parameter in models with a standard scalar field non - minimally coupled to gravity , for example , in induced gravity models  @xcite .    to obtain bounce solutions in cosmological models with standard scalar fields one can use non - minimal coupling between the scalar fields and gravity and consider the following action : @xmath2 where @xmath3 is the determinant consisting of metric tensor components @xmath4
, @xmath5 is the ricci scalar , @xmath6 and @xmath7 are differentiable functions of a scalar field @xmath8 .    from the condition that at the bounce point the hubble parameter @xmath9
, the potential should be negative @xmath10 .
the characteristic property of models with non - positive definite potential is the existence of the unreachable domain on the phase plane , which corresponds to non - real hubble parameter  @xcite .    in this paper
we consider models with bounce solutions that generalize the integrable model proposed in  @xcite with the ricci scalar @xmath5 as an integral of motion .
this model has the potential in the form @xmath11 , with @xmath12 and @xmath13 , and the standard quadratic coupling @xmath14 with @xmath15 .
as known in the spatially flat flrw metric the ricci scalar @xmath5 is a function of the hubble parameter @xmath1 and its derivative .
so , the condition that @xmath5 is a constant defines the hubble parameter as a solution of the first order differential equation .
thereby , we get only monotonic behavior of the hubble parameter ( formulae are presented in section 3 ) . in particular
, this integrable model contains bouncing cosmological solutions with smooth future behavior tending to de sitter solution ( while past behavior shows emerging of a universe at a point with zero effective newtonian constant ) .
it has been shown in  @xcite that the integrable cosmological model  @xcite with a constant @xmath5 belongs to one - parameter set of integrable models with one and the same function @xmath16 and different potentials @xmath17 . however , the ricci scalar @xmath5 is not an integral to motion for these integrable models . in this paper
we show that the model , proposed in  @xcite , is unique in this sense that the property of having @xmath5 to be a constant disappears , when one changes coupling function or scalar field potential ( or both ) .    to construct more realistic models with both a bounce and non - monotonic behavior of the hubble parameter one can consider non - integrable models that are close to the above - mention integrable model with a constant @xmath5 . in
@xcite the authors got such types of the hubble parameter considering non - integrable models with the same potential and other values of @xmath18 . in the present paper
we mostly vary the potential in searching for conditions of bouncing behavior similar to studied in  @xcite to exist in these more general models . in section  4
we consider the model with @xmath15 and the higgs - type potential multiplying by a negative constant .
in section  5 we consider models with different positive values of @xmath18 , a cosmological constant and different even degrees of monomial potentials .
using numerical calculations we demonstrate the possibility of different evolutions of the hubble parameters in the cases of potentials of the forth and sixth degrees .
we analyze the reasons of such behavior and find the corresponding conditions of the model parameters .
the paper is organized as follows . in sect .
2 we remind a reader the general bounce conditions in non - minimally coupled scalar field theory as well as the conception of the effective potential . in sect .
3 we show that the choice of potential in the paper  @xcite leading to the property that the curvature @xmath5 is the integral of motion is the unique one . on the contrary , in sect .  4 and 5 we show that bounce behavior admits much larger class of potentials . in sect .
4 bounce behavior is studied for potentials with quadratic term in addition to quartic and constant terms and @xmath15 .
the models with @xmath19 are considered in sect .
6 contains a summary of the result obtained .
in the spatially flat flrw metric with the interval : @xmath20 we obtain from action ( [ s1 ] ) the following equations : @xmath21 @xmath22 + 2u'\left[\ddot{\varphi}+2h\dot{\varphi}\right ] = v-\left[2u''+\frac12\right ] \dot{\varphi}^2,\ ] ] @xmath23+v ' = 0 , \label{kgoc}\ ] ] where a `` dot '' means a derivative with respect to the cosmic time @xmath24 and a `` prime '' means a derivative with respect to the scalar field @xmath8 .
the function @xmath25 is the scale factor , its logarithmic derivative @xmath26 is the hubble parameter .
we note that an effective gravitation constant in the model considered is @xmath27 the dynamics of flrw universe can be prolonged smoothly into the region of @xmath28 ( see , for example  @xcite ) , however , since any anisotropic corrections are expected to diverge while @xmath29 tends to infinity  @xcite . in this paper , we analyze such bounce solutions that @xmath30 for any @xmath31 and conditions of their existence .    if a solution of eqs .
( [ frequoc00])([kgoc ] ) has such a point @xmath32 that @xmath33 and @xmath34 , then it is a bounce solution .
let us find conditions that are necessary for the existence of a bounce in models with action  ( [ s1 ] ) .
using eq .
( [ frequoc00 ] ) we get that from @xmath33 it follows @xmath35 .
subtracting equation  ( [ frequoc00 ] ) from equation  ( [ frequocii ] ) , we obtain @xmath36    therefore , at the bounce point we also get @xmath37v,\ ] ] where functions @xmath16 and @xmath17 and their derivatives are taken at the point @xmath38 . the condition @xmath34 gives the restriction on functions @xmath16 and @xmath17 at the bounce point .    to analyze the dynamic
it is suitable to introduce a new variable and present equations ( [ frequoc00 ] ) and ( [ frequ21 ] ) in the form similar to the friedmann equations in the einstein frame .
following  @xcite , we introduce the effective potential @xmath39 if the constant @xmath40 , then @xmath41 coincides with the potential of the corresponding model in the einstein frame .
note that we do not transform the metric and consider the model in the jordan frame only .
for our purposes any positive value of @xmath42 is suitable ( in  @xcite it was chosen @xmath43 ) .
also in  @xcite functions @xmath44 have been introduced . if @xmath45 , then @xmath46 as well .    in terms of these functions eqs .
( [ frequoc00 ] ) and ( [ frequ21 ] ) take the following form : @xmath47 @xmath48    as known the hubble parameter is a monotonically decreasing function in models with a standard scalar field minimally coupled to gravity ( @xmath6 is a constant ) . for a model with an arbitrary positive @xmath6 function @xmath49
has the same property .
equation  ( [ frequoc00 ] ) is a quadratic equation in @xmath50 that has the following solutions : @xmath51    the values of the function @xmath49 that correspond to @xmath52 are @xmath53 }      = \pm\sqrt{\frac{a}{3}{\dot\varphi}^2+\frac{2}{3}k^2v_{eff } } \ , .\ ] ]    if the potential @xmath17 is negative at some values of @xmath8 , whereas the function @xmath16 is positive at these points , then one should restrict the domain of absolute values of  @xmath54 from below to get a real hubble parameter .
in other words , there exists the unreachable domain on the phase plane  @xcite .
the boundary of this domain is defined by the condition @xmath55 . if @xmath56 then the sign of the function @xmath57 can not change .
if in some domain @xmath58 is negative , then the sign of @xmath57 can be changed from plus to minus , but not vice verse .
so , the sign of @xmath57 can not be changed twice .
equations  ( [ frequoc00])([kgoc ] ) can be transformed into the following system of the first order equations which is useful for numerical calculations and analysis of stability : @xmath59u'+2uv'}{2\left(3 { u'}^2 + u\right)},\\ \dot h&={}-\frac{2u''+1}{4(3{u'}^2+u)}\psi^2+\frac{2u'}{{3{u'}^2+u}}h\psi -\frac{6{u'}^2}{3{u'}^2+u}h^2+\frac{u'v'}{2(3{u'}^2+u)}\ , .
\end{split } \right . \label{fosequ}\ ] ] if eq .
( [ frequoc00 ] ) is satisfied in the initial moment of time , then from system ( [ fosequ ] ) it follows that eq .
( [ frequoc00 ] ) is satisfied at any moment of time .
so , one can use eq .
( [ frequoc00 ] ) to fix initial conditions of system  ( [ fosequ ] ) . from eqs .
( [ equp ] ) and ( [ fr21qm ] ) it is easy to get the following system  @xcite : @xmath60 system of equations ( [ fr21qm ] ) and ( [ system2 ] ) is equivalent to system ( [ fosequ ] ) .
the de sitter solutions corresponds to @xmath61 , and hence , @xmath62 .
the corresponding hubble parameter is @xmath63    we analyze the stability of de sitter solutions with @xmath64 and @xmath65 only . from ( [ hf ] )
it follows that @xmath66 . in this case
the hubble parameter in the neighborhood of @xmath67 is uniquely defined by eq .
( [ frequoc00 ] ) , so it is enough to consider system ( [ system2 ] ) to analyze the stability of de sitter solution . substituting the de sitter solution with the first order perturbations : @xmath68 into ( [ system2 ] ) , we get the following linear system on @xmath69 and @xmath70 : @xmath71    we find eigenvalues for the system ( [ linsystem ] ) : @xmath72    the real part of @xmath73 always , the condition that the real part of @xmath74 is equivalent to @xmath75 . therefore , for arbitrary differentiable functions @xmath17 and @xmath76 , the model has a stable de sitter solution with @xmath64 only if the potential @xmath41 has a minimum  @xcite and @xmath77 .    using the obtained values of @xmath78 ,
it is easy to find  @xcite that the de sitter point is a stable node ( the scalar field decreases monotonically ) at @xmath79 and a stable focus ( the scalar field oscillations exist ) at @xmath80
in the paper  @xcite the cosmological model with a constant ricci scalar @xmath5 has been considered . in this section
we demonstrate how this form of function @xmath16 and the corresponding potential can be found from the requirement that the ricci scalar is an integral of motion .    in the spatially flat flrw metric @xmath81 .
from eqs .
( [ frequoc00])([kgoc ] ) we get @xmath82 @xmath83 combining eqs .
( [ equtrace ] ) and ( [ kgequ ] ) we obtain    @xmath84    from the structure of eq .  ( [ equr ] ) it is easy to see that the simplest way to get a constant @xmath5 is to choose such @xmath6 that @xmath85 where @xmath86 is a constant .
the first two equations of system  ( [ equ ] ) has the following solution @xmath87 where @xmath88 is an integration constant . without loss of generality we put @xmath89 .
for such a choice of @xmath6 we get eq .
( [ equr ] ) as follows : @xmath90 considering eq .
( [ equr2 ] ) as a differential equation for @xmath7 , for a constant @xmath5 we get the following solution : @xmath91 where @xmath92 is an integration constant , for other constants we choose @xmath93    thus , requiring that the ricci scalar is a constant one can define both functions @xmath6 and @xmath7 . to get a positive @xmath29 for some values of @xmath8 we choose @xmath94 .
substituting @xmath95 and @xmath96 into eq .
( [ bouncecond2 ] ) , we get that the condition @xmath97 is equivalent to @xmath98 , hence , from @xmath99 it follows @xmath13 .
this integrable cosmological model has been considered in  @xcite , where the behavior of bounce solutions has been studied in detail .
equation @xmath100 is a differential equation for the hubble parameter : @xmath101 multiplying this equation by @xmath102 , we get @xmath103= 0.\ ] ] therefore , there exists the following integral of motion : @xmath104    equation  ( [ equhr ] ) with a positive @xmath105 has two possible real solutions in dependence of the initial conditions : @xmath106 where @xmath107 is an integration constant . for convenience
we locate the bounce at @xmath108 , so hereafter @xmath109 .
note that the behavior of the hubble parameter does not depend on the specific dynamics of the scalar field @xmath8 , because two - parametric set of functions @xmath110 corresponds to one - parametric set of @xmath50 .
the monotonically increasing hubble parameter @xmath111 is not suitable for construction of a realistic cosmological scenario , whereas @xmath112 is not a bounce solution . in  @xcite
the authors slightly modify the function @xmath16 to get a bounce solution with non - monotonic behaviour of the hubble parameters . in this section
we modify the potential instead of the function @xmath16 and analyze the obtained bounce solutions . in other words ,
we choose the function @xmath113 , given by ( [ uc ] ) , with @xmath114 . as a minimal generalization of the potential @xmath96 we consider the following potential @xmath115 where @xmath116 are constants .
we consider the case @xmath13 only , because the integrable model with the potential @xmath96 has a bounce solution and a stable de sitter solution only if @xmath98 and @xmath13 .
we plan to consider the evolution of bounce solutions fixing the initial conditions at bounce point .
without loss of generality we can consider as initial value @xmath117 only .
the value of @xmath118 is defined by the condition @xmath119 , hence , @xmath120 .
our first goal is to find conditions on the coefficients of the potential @xmath121 that corresponds to the existence of bounce solutions . for @xmath113 ,
a bounce solution exists only if there exists such a point @xmath122 that the following conditions are satisfied  @xcite : @xmath123 so , @xmath124 , and we obtain the following conditions on parameters @xmath116 : @xmath125 thus , at least one of the constants @xmath126 or @xmath127 should be positive .
a bounce solution has the physical sense only if @xmath128 .
the condition @xmath129 means @xmath130 .
the effective potential ( [ veff ] ) is @xmath131 the even potential @xmath41 has an extremum at @xmath132 and at points @xmath133 we consider such values of parameter that points @xmath134 are real and @xmath135 .
note that for integrable model ( @xmath13 , @xmath136 , @xmath12 ) @xmath137 are non - zero real numbers and the condition @xmath138 is equivalent to @xmath139 .    from  ( [ bouncecond ] ) , @xmath13 and @xmath130 we get @xmath140 so , the model with a bounce solution has real @xmath134 only at @xmath141 using these conditions , we get @xmath142 so , the potential @xmath41 has a minimum at @xmath132 and maxima at @xmath143 .
a few examples of the effective potential @xmath41 with two maxima at nonzero points @xmath134 and a minimum at @xmath132 are presented in fig .
[ veffc0 ] .     at different values of parameters . in all pictures
we choose @xmath144 .
the values of parameters are @xmath145 , @xmath146 ( left picture ) , @xmath147 , @xmath146 ( middle picture ) , @xmath147 , @xmath148 ( right picture ) .
the parameter @xmath149 ( black curve ) , @xmath150 ( red curve ) , @xmath151 ( blue curve ) , @xmath152 ( green curve ) , and @xmath153 ( cyan curve ) .
the yellow curve in the middle picture corresponds to @xmath154.,title="fig:",width=185 ]          at different values of parameters . in all pictures we choose @xmath144 .
the values of parameters are @xmath145 , @xmath146 ( left picture ) , @xmath147 , @xmath146 ( middle picture ) , @xmath147 , @xmath148 ( right picture ) .
the parameter @xmath149 ( black curve ) , @xmath150 ( red curve ) , @xmath151 ( blue curve ) , @xmath152 ( green curve ) , and @xmath153 ( cyan curve ) .
the yellow curve in the middle picture corresponds to @xmath154.,title="fig:",width=185 ]         at different values of parameters . in all pictures we choose @xmath144 .
the values of parameters are @xmath145 , @xmath146 ( left picture ) , @xmath147 , @xmath146 ( middle picture ) , @xmath147 , @xmath148 ( right picture ) .
the parameter @xmath149 ( black curve ) , @xmath150 ( red curve ) , @xmath151 ( blue curve ) , @xmath152 ( green curve ) , and @xmath153 ( cyan curve ) .
the yellow curve in the middle picture corresponds to @xmath154.,title="fig:",width=185 ]    stable de sitter solution at @xmath155 corresponds to @xmath156 .
therefore , such a solution exists for @xmath157 only .
it means that the stable de sitter solution exists only if @xmath158 . in this case @xmath159 .
the potential @xmath41 has the following zeros : @xmath160    at @xmath13 , @xmath12 one gets @xmath161 .
therefore , @xmath162 and there exist only two real roots : @xmath163 the bounce point @xmath164 belongs to the interval @xmath165 .
we get the condition @xmath166
.      for @xmath113 and an arbitrary potential system ( [ fosequ ] ) has the following form  @xcite : @xmath167 . \end{split } \label{fosequc}\ ] ] we integrate this system with @xmath168 numerically .
we consider a positive @xmath122 such that @xmath169 .
the evolution of the scalar field starts at the bounce point with a negative velocity , defined by relation @xmath170 the field @xmath8 can come to zero passing the maximum of the potential .
so , we keep in mind that the following subsequence of inequalities : @xmath171    in the case @xmath12 there are three possible evolutions of the bounce solutions , depending on whether the solution passes the points of the maximum of @xmath41 or not . in the left and middle pictures of fig .
[ threesolutions ] we present an example of three possible behaviors of the bounce solutions .
the solution with less initial value of @xmath8 tends to infinity ( blue curve ) , whereas the bounce solutions with greater initial values tends to zero ( cyan curve ) or to minus infinity as a monotonically decreasing function ( green curve ) .
all trajectories start at bounce points .
the corresponding behaviors of the hubble parameter are presented in the right picture of fig .
[ threesolutions ] ( colors the hubble parameter evolutions on this picture coincide to the colors of the corresponds phase trajectories ) .     and
the values of constants are @xmath144 , @xmath147 , @xmath172 , and @xmath173 .
the initial conditions are @xmath174 , @xmath175 ( blue line ) , @xmath176 , @xmath177 ( cyan line ) , and @xmath178 , @xmath179 ( green line ) .
the black curves are the lines of the points that correspond to @xmath9 .
the red dash lines are the boundary of unreachable domain defined by @xmath55 .
the blue point lines correspond to @xmath180 .
a zoom of the central part of phase plane is presented in the middle picture .
the corresponding hubble parameters as functions of cosmic time are presented on the right picture.,title="fig:",width=178 ]         and @xmath121 .
the values of constants are @xmath144 , @xmath147 , @xmath172 , and @xmath173 .
the initial conditions are @xmath174 , @xmath175 ( blue line ) , @xmath176 , @xmath177 ( cyan line ) , and @xmath178 , @xmath179 ( green line ) .
the black curves are the lines of the points that correspond to @xmath9 .
the red dash lines are the boundary of unreachable domain defined by @xmath55 .
the blue point lines correspond to @xmath180 .
a zoom of the central part of phase plane is presented in the middle picture .
the corresponding hubble parameters as functions of cosmic time are presented on the right picture.,title="fig:",width=178 ]         and @xmath121 .
the values of constants are @xmath144 , @xmath147 , @xmath172 , and @xmath173 .
the initial conditions are @xmath174 , @xmath175 ( blue line ) , @xmath176 , @xmath177 ( cyan line ) , and @xmath178 , @xmath179 ( green line ) .
the black curves are the lines of the points that correspond to @xmath9 .
the red dash lines are the boundary of unreachable domain defined by @xmath55 .
the blue point lines correspond to @xmath180 .
a zoom of the central part of phase plane is presented in the middle picture .
the corresponding hubble parameters as functions of cosmic time are presented on the right picture.,title="fig:",width=178 ]    let us compare two bounce solutions with initial conditions : @xmath181 and @xmath182 and negative initial values of @xmath118 .
the corresponding hubble parameter that is equal to zero at a bounce point is given by ( [ hpm ] ) : @xmath183 .
the initial value of @xmath118 for the first solution we denote as @xmath184 .
the potential @xmath10 for all @xmath185 $ ] , therefore , @xmath186 for any solution that pass @xmath187 .
thus , at some moment of time @xmath188 the solution @xmath189 comes at the point @xmath187 : @xmath190 .
the value of its velocity at the bounce point @xmath191 has a minimal absolute value by comparison with any solutions that pass via @xmath187 with a negative velocity .
it follows from ( [ frequoc00 ] ) , because both @xmath192 , and @xmath193 are positive .
it proves that if a solution with the bounce point @xmath187 passes through the maximum of the potential and tends to zero , then any solution with the bounce point @xmath194 tends to zero as well .
all bounce solutions with negative initial value of @xmath118 come in the domain where @xmath195 , so , @xmath196 . if a solution does not pass the maximum of the effective potential , then after some moment @xmath8 starts to grow and this solution comes to antigravity domain with @xmath197 ( the blue curve in fig .
[ threesolutions ] denotes an example of such a solution ) .
solutions passing the maximum of the potential can be different and need considering in detail .    if @xmath10 for all @xmath8 , then @xmath110 is a monotonic function , because at any point @xmath198 .
similar dynamics is possible even if @xmath195 ( see blue curve in fig .
[ threesolutions ] )
.    for @xmath12 there exists the stable de sitter solution @xmath155 and @xmath199 .
it is a stable node at @xmath200 and a stable focus in the opposite case @xmath201 . in fig .
[ threesolutions ] ( cyan curves ) , in the left picture of fig .
[ phaseportnode ] , and fig .
[ phaseportv4c0 ] solutions in the case of a stable focus are presented .
an example with a stable node at @xmath132 is given in the middle picture of fig .
[ phaseportnode ] .
one can see that in the case of a stable node the presented hubble parameter ( the right picture of fig .  [ phaseportnode ] ) is close to the hubble parameter obtained in the paper  @xcite that has a maximum .
note that in the integrable case @xmath136 and there exists a stable node at @xmath132 .    , @xmath147 , @xmath202 and @xmath203 .
initial values are @xmath204 and @xmath205 ( cyan curve ) , @xmath206 and @xmath207 ( green curve ) .
the cyan curve is an example of a stable focus . on the middle and right pictures
the example of a stable node at @xmath132 is presented .
the field @xmath8 ( blue line ) and the hubble parameter ( red line ) as functions of the cosmic time are presented in the right picture .
the values of parameters are @xmath208 , @xmath209 , @xmath172 and @xmath210 .
the initial conditions of the bounce solution are @xmath211 , @xmath212 , and @xmath119.,title="fig:",width=178 ]       , @xmath147 , @xmath202 and @xmath203 .
initial values are @xmath204 and @xmath205 ( cyan curve ) , @xmath206 and @xmath207 ( green curve ) .
the cyan curve is an example of a stable focus . on the middle and right pictures the example of a stable node at @xmath132
is presented .
the field @xmath8 ( blue line ) and the hubble parameter ( red line ) as functions of the cosmic time are presented in the right picture .
the values of parameters are @xmath208 , @xmath209 , @xmath172 and @xmath210 .
the initial conditions of the bounce solution are @xmath211 , @xmath212 , and @xmath119.,title="fig:",width=178 ]        , @xmath147 , @xmath202 and @xmath203 .
initial values are @xmath204 and @xmath205 ( cyan curve ) , @xmath206 and @xmath207 ( green curve ) .
the cyan curve is an example of a stable focus . on the middle and right pictures
the example of a stable node at @xmath132 is presented .
the field @xmath8 ( blue line ) and the hubble parameter ( red line ) as functions of the cosmic time are presented in the right picture .
the values of parameters are @xmath208 , @xmath209 , @xmath172 and @xmath210 .
the initial conditions of the bounce solution are @xmath211 , @xmath212 , and @xmath119.,title="fig:",width=178 ]     and @xmath121 is presented in the left picture .
the values of constants are @xmath144 , @xmath147 , @xmath202 , @xmath173 .
the initial conditions are @xmath206 and @xmath213 .
a zoom of the central part of phase plane is presented in the middle picture .
the hubble parameter ( red ) and the scalar field ( blue ) of functions of cosmic time are presented in the right picture.,title="fig:",width=178 ]         and @xmath121 is presented in the left picture .
the values of constants are @xmath144 , @xmath147 , @xmath202 , @xmath173 .
the initial conditions are @xmath206 and @xmath213 .
a zoom of the central part of phase plane is presented in the middle picture .
the hubble parameter ( red ) and the scalar field ( blue ) of functions of cosmic time are presented in the right picture.,title="fig:",width=178 ]         and @xmath121 is presented in the left picture .
the values of constants are @xmath144 , @xmath147 , @xmath202 , @xmath173 .
the initial conditions are @xmath206 and @xmath213 . a zoom of the central part of phase plane is presented in the middle picture .
the hubble parameter ( red ) and the scalar field ( blue ) of functions of cosmic time are presented in the right picture.,title="fig:",width=178 ]    in the left picture of fig .
[ phaseportnode ] the phase trajectories have been constructed for @xmath144 , @xmath147 , @xmath202 and @xmath203 .
let us now change the value of @xmath127 only and consider the model with @xmath173 . in fig .
[ phaseportv4c0 ] the corresponding phase trajectory is presented .
we see that now the bounce solution that starts at @xmath206 tends to zero and finishes at the point @xmath132 .
we can see that the trajectories that revolve around @xmath214 points looks similar at @xmath203 and at @xmath173 .
let us consider now the phase trajectory at @xmath215 that is presented in fig .
[ phaseportv4cnega ] .
we see that trajectories are similar at the beginning only .
the scalar field tends to infinity and the system comes to antigravity domain with @xmath216 .
the form of the effective potential does not depend essentially from the sign value of @xmath127 ( see fig .
[ veffc0 ] ) .
but the sign of @xmath217 are different . at the same time the behavior of bounce solutions are essentially different . in the right pictures of fig .
[ phaseportv4c0 ] and fig .
[ phaseportv4cnega ] we show that the behavior of the hubble parameter also essentially depends on the sign of  @xmath127 .     and
@xmath121 is presented in the left picture .
the values of constants are @xmath144 , @xmath147 , @xmath202 , @xmath215 .
the initial conditions are @xmath206 and @xmath218 . a zoom of the central part of phase plane is presented in the middle picture .
the hubble parameter ( red ) and the scalar field ( blue ) of functions of cosmic time are presented in the right picture.,title="fig:",width=178 ]         and @xmath121 is presented in the left picture .
the values of constants are @xmath144 , @xmath147 , @xmath202 , @xmath215 .
the initial conditions are @xmath206 and @xmath218 . a zoom of the central part of phase plane is presented in the middle picture .
the hubble parameter ( red ) and the scalar field ( blue ) of functions of cosmic time are presented in the right picture.,title="fig:",width=178 ]         and @xmath121 is presented in the left picture .
the values of constants are @xmath144 , @xmath147 , @xmath202 , @xmath215 .
the initial conditions are @xmath206 and @xmath218 . a zoom of the central part of phase plane is presented in the middle picture .
the hubble parameter ( red ) and the scalar field ( blue ) of functions of cosmic time are presented in the right picture.,title="fig:",width=178 ]    the difference between the solutions of system ( [ fosequc ] ) with a positive and a negative @xmath127 is demonstrated in fig .
[ c0pm ] as well .
the cyan curves correspond to @xmath203 , whereas the red curves correspond to @xmath149 .
we see that the phase trajectories of the field @xmath8 and the hubble parameter are similar in the beginning , but stand essentially different in the future .
we see that the behavior of solutions essentially depends on the sign of @xmath127 .
let us understand the reason of this dependence .
, @xmath147 , @xmath219 . the parameter @xmath203 for cyan curves and @xmath149 for red curves .
the initial conditions of bounce solution are @xmath204 , @xmath220 ( cyan curve ) and @xmath221 ( red curve).,title="fig:",width=272 ]        , @xmath147 , @xmath219 .
the parameter @xmath203 for cyan curves and @xmath149 for red curves .
the initial conditions of bounce solution are @xmath204 , @xmath220 ( cyan curve ) and @xmath221 ( red curve).,title="fig:",width=272 ]    let us consider the domain @xmath222 . from ( [ hpm ] )
it follows that @xmath223    if @xmath224 for all @xmath225 , then this condition is always satisfied and @xmath8 tends to a minimum of @xmath41 .
we see such evolution in the right picture of fig .
[ phaseportv4c0 ] .
the constants @xmath116 are such that @xmath195 for all @xmath226 $ ] .
for suitable initial @xmath54 the evolution is finished at @xmath132 .
note that in this case if @xmath227 in the moment when potential stands positive and the function @xmath8 tends to zero , then @xmath227 at any moment in future .    if the constants @xmath116 are such that the potential change the sign and @xmath228 then the evolution is different ( see the right picture in fig .  [ phaseportv4cnega ] ) .
the hubble parameter becomes negative and positive again , so , there are two bounce points . after that the hubble parameter tends to infinity .
let us consider this case in detail .
if @xmath229 , then there is a restricted domain in the neighborhood of @xmath214 point on phase plane such that the values of the scalar field @xmath8 and its derivative correspond to non - real value of the hubble parameter .
the boundary of this domain is defines by equation @xmath55 . on fig .
[ phaseportv4cnega ] we see that the phase trajectory rotates around this domain .
the trajectory can not cross the boundary , but can touch it .
let us show that all such trajectories touch the boundary at some finite moment of time .
let for some moments of time @xmath230 and @xmath231 we have @xmath232 , then , using @xmath233 and formula ( [ fr21qm ] ) , we get @xmath234 where @xmath235 is a negative number .
therefore , this integral has a finite negative value .
for any circle value of @xmath57 decreases on some positive value , which does nt tend to zero during evolution , when number of circles increase .
we come to conclusion that only a finite number of circles is necessary to get the value @xmath55 . at this point @xmath236 as well , so the function @xmath57 changes the sign .
when @xmath55 two possible values of the hubble parameters : @xmath237 and @xmath238 coincide . at this moment the value of the hubble parameter changes from @xmath237 to @xmath238 .
the value of the function @xmath57 continues to decrease , so , the distance between the trajectory and the boundary of unreachable domain increases .
we do not say that its increase monotonically but @xmath57 that is a characteristic of this distance increases on a finite quantity after any circle .
so , after some finite number of circles the absolute value of @xmath8 becomes more then @xmath134 .
after this moment @xmath8 monotonically tends to infinity and at some finite moment we get @xmath239 .
thus , the final of trajectory will be in antigravity domain always .
note that in the domain with @xmath195 the hubble parameter is uniquely defined as a function @xmath8 and @xmath118 by  ( [ hpm ] ) .
if @xmath240 , then whole evolution of bounce solutions is evolution a solution of the second order system , whereas for @xmath241 the third order system  ( [ fosequc ] ) with the additional condition  ( [ frequoc00 ] ) is not equivalent to any second order system .
in this section we consider the scalar field potential of the form @xmath242 where @xmath243 is an even natural number , and following non - minimally coupled function @xmath244 where @xmath18 is positive .
we are interested in cosmological scenarios in the physical region @xmath245 , where the bounce occurs at first and after that stable de sitter solution @xmath246 is realized .
so , the potential should change the sign .
the potential @xmath41 is an even function , hence , it has an extremum at @xmath132 . if we additionally suppose that @xmath155 , then we get the condition @xmath12 .
therefore , the bounce condition @xmath247 for the chosen potential @xmath17 and the function @xmath16 gives ( see ): @xmath248    from the condition @xmath249 it follows @xmath250    taking into account the inequalities @xmath251 and @xmath12 we obtain from conditions   and that the bounce exists for    * \(i ) @xmath252 * * ( 1 ) . if @xmath253 ,   then  @xmath254 . * * ( 2 ) . if @xmath255 ,   then  @xmath256 . * \(ii ) @xmath257 + @xmath256 .    here
@xmath122 is the value of the scalar field at the bounce .
we see that there is an additional restriction on the location of points of bounce in comparison with the cases studied earlier in  @xcite , appearing when @xmath252 ( the case ( i1 ) ) .
this happens due to @xmath258 condition of the bounce , which can be violated in the case of small enough @xmath18 , causing a recollapse instead of bounce .
the bounce point should corresponds to @xmath245 that gives the following condition @xmath259    let us study properties of the effective potential @xmath41 to analyse the lyapunov stability of de sitter solutions .
the effective potential @xmath58 have zeros for @xmath260 they are located in the region @xmath245 only if @xmath261 thus , we have received a restriction from above for @xmath18 .    the fist derivative @xmath41 with respect to @xmath8 for the chosen functions @xmath17 and @xmath16 is given by @xmath262 therefore , the effective potential has the extremum at the point @xmath132 for any even @xmath243 .
non - zero extrema @xmath134 for @xmath263 are @xmath264 for @xmath265 values
@xmath266 are roots of the following cubic equation @xmath267 which has three real roots for @xmath268 .
we consider only @xmath269 ( see ) , hence , there are three real roots in this case .
only positive ones have the physical sense ( @xmath270 ) .
we calculate @xmath271 for de sitter solution @xmath155 , @xmath272 .
then it is stable for @xmath273 , namely , @xmath274 from this it follows that there is the restriction from below for the coupling constant @xmath275 in the case of @xmath263 . on the other hand , for @xmath276 the interval of possible @xmath18 allowing evolution towards de sitter solution after the bounce is not restricted from below ( we remind that @xmath277 and we consider only positive values of @xmath18 in this paper ) .
when de sitter solution is stable we have ( see and ) a monotonic decreasing of the scalar field for @xmath278 or scalar field oscillations for @xmath279 this result agrees with the one received in  @xcite for @xmath280 .    let us summarize our requirements to the effective potential @xmath281 and the parameter @xmath18 :    *  @xmath12 ,  @xmath251 , *  @xmath282   for  @xmath276 , *  @xmath283   for  @xmath263 .      the case of @xmath280 and @xmath15 is an integrable model with monotonic behavior of the hubble parameter only .
the change of parameter @xmath18 can not only give a hubble parameter with one maximum ( see  @xcite ) , but also oscillating hubble parameter that tends to a constant .
such solution has been found at @xmath284 and are presented in fig .
[ n4xi20 ] .
, @xmath285 .
the parameters are @xmath284 , @xmath286 , @xmath287 , @xmath288 .
the initial conditions are @xmath289 and @xmath290 ( green line ) , @xmath291 and @xmath292 ( blue line ) .
the black curves are the lines of the points that correspond to @xmath9 .
the blue point lines correspond to @xmath180 .
the hubble parameter @xmath50 is presented in the right picture .
, title="fig:",width=181 ] , @xmath285 .
the parameters are @xmath284 , @xmath286 , @xmath287 , @xmath288 .
the initial conditions are @xmath289 and @xmath290 ( green line ) , @xmath291 and @xmath292 ( blue line ) .
the black curves are the lines of the points that correspond to @xmath9 .
the blue point lines correspond to @xmath180 .
the hubble parameter @xmath50 is presented in the right picture .
, title="fig:",width=181 ] , @xmath285 .
the parameters are @xmath284 , @xmath286 , @xmath287 , @xmath288 .
the initial conditions are @xmath289 and @xmath290 ( green line ) , @xmath291 and @xmath292 ( blue line ) .
the black curves are the lines of the points that correspond to @xmath9 .
the blue point lines correspond to @xmath180 .
the hubble parameter @xmath50 is presented in the right picture .
, title="fig:",width=181 ]    we explore numerically cases with @xmath280 and @xmath265 applying the first order system of differential equations .
initial data are chosen @xmath293 ,  @xmath294  and   @xmath295 is real and calculated using @xmath237 ( see ) .
values of the scalar field @xmath122 and its time derivative @xmath296 are at the bounce line @xmath297 ( we denote @xmath298 ) . in figs .
[ n4xi01 ] ,  [ n4xi01767 ] and [ v6ksi086 ] parameters @xmath286 , @xmath299 are chosen for all these plots and , therefore , stable de sitter solution is @xmath155 , @xmath300 .    what is more interesting
, numerical integration reveal another possible fate of a trajectory with initial conditions in the region suitable for bounce . in fig .
[ n4xi01767 ] the outcome of the evolution with given initial value have been marked by different colors .
black and green zones correspond to a bounce followed by a smooth evolution towards de sitter attractor ( black ) or @xmath8-turn then growing @xmath8 , leading to crossing @xmath180 boundary ( green ) .
however , the blue zone exists starting from which a trajectory does not experience any bounce and fall into a singularity .
this happens because a trajectory can touch the unreachable part of phase space and jump to another branch of two possible solutions for @xmath1 ( see  [ hpm ] ) .
we have described such a phenomenon in the previous section , where this happens near @xmath132 due to negative cosmological constant @xmath127 .
here it happens in the bounce zone , where we also have negative potential and unreachable part of the phase space .
after branch changing , the trajectory can not go through a bounce because bounce region for other branch is located in different phase space region .
one can see in fig .
[ n4xi01 ] that de sitter attractors do not exist for sufficiently small positive @xmath18 .    ,
@xmath285 : singularities with @xmath301 , @xmath302 ( green points ) and with @xmath303 , @xmath304 ( blue points ) .
the parameters are @xmath305 , @xmath286 , @xmath306 , @xmath299 . if the bounce occurs then we mark initial data with circles and in otherwise  points .
the red domain corresponds to unreachable domain.,width=453 ]    , @xmath285 in the left picture : de sitter solution ( black circles ) , the analog of ruzmaikin one ( green circles ) and the singularity with @xmath303 , @xmath304 ( blue points ) .
if the bounce occurs , then we mark initial data with circles and in otherwise  points .
the parameters are @xmath307 , @xmath286 , @xmath306 , @xmath299 .
three phase trajectories ( middle picture ) and dependencies @xmath50 ( right picture ) corresponding these attractors are plotted for initial data @xmath308 and @xmath309 ( blue curve ) , @xmath310 ( green curve ) , @xmath311 ( black curve ) .
the gold trajectory starts at the bounce point @xmath312 .
the black dash line corresponds to @xmath9 .
the blue point line corresponds to @xmath180 .
the red color in the left picture indicates the unreachable domain .
the red dash line in the middle picture is the boundary of unreachable domain defined by @xmath55 .
, title="fig:",width=185 ] , @xmath285 in the left picture : de sitter solution ( black circles ) , the analog of ruzmaikin one ( green circles ) and the singularity with @xmath303 , @xmath304 ( blue points ) .
if the bounce occurs , then we mark initial data with circles and in otherwise  points .
the parameters are @xmath307 , @xmath286 , @xmath306 , @xmath299 .
three phase trajectories ( middle picture ) and dependencies @xmath50 ( right picture ) corresponding these attractors are plotted for initial data @xmath308 and @xmath309 ( blue curve ) , @xmath310 ( green curve ) , @xmath311 ( black curve ) .
the gold trajectory starts at the bounce point @xmath312 .
the black dash line corresponds to @xmath9 .
the blue point line corresponds to @xmath180 .
the red color in the left picture indicates the unreachable domain .
the red dash line in the middle picture is the boundary of unreachable domain defined by @xmath55 . ,
title="fig:",width=181 ] , @xmath285 in the left picture : de sitter solution ( black circles ) , the analog of ruzmaikin one ( green circles ) and the singularity with @xmath303 , @xmath304 ( blue points ) .
if the bounce occurs , then we mark initial data with circles and in otherwise  points .
the parameters are @xmath307 , @xmath286 , @xmath306 , @xmath299 .
three phase trajectories ( middle picture ) and dependencies @xmath50 ( right picture ) corresponding these attractors are plotted for initial data @xmath308 and @xmath309 ( blue curve ) , @xmath310 ( green curve ) , @xmath311 ( black curve ) .
the gold trajectory starts at the bounce point @xmath312 .
the black dash line corresponds to @xmath9 .
the blue point line corresponds to @xmath180 .
the red color in the left picture indicates the unreachable domain .
the red dash line in the middle picture is the boundary of unreachable domain defined by @xmath55 .
, title="fig:",width=181 ]     in the case of the sixth degree potential and @xmath313 .
dependencies @xmath110 and @xmath50 for initial data : @xmath314 and @xmath315 ( gold curve ) , @xmath314 and @xmath316 ( blue curve ) , @xmath317 and @xmath316 ( green curve ) .
parameters are chosen @xmath265 , @xmath313 , @xmath286 , @xmath318 , @xmath299 .
the critical value of @xmath319 .
, title="fig:",width=181 ]   in the case of the sixth degree potential and @xmath313 .
dependencies @xmath110 and @xmath50 for initial data : @xmath314 and @xmath315 ( gold curve ) , @xmath314 and @xmath316 ( blue curve ) , @xmath317 and @xmath316 ( green curve ) .
parameters are chosen @xmath265 , @xmath313 , @xmath286 , @xmath318 , @xmath299 .
the critical value of @xmath319 .
, title="fig:",width=181 ]   in the case of the sixth degree potential and @xmath313 .
dependencies @xmath110 and @xmath50 for initial data : @xmath314 and @xmath315 ( gold curve ) , @xmath314 and @xmath316 ( blue curve ) , @xmath317 and @xmath316 ( green curve ) .
parameters are chosen @xmath265 , @xmath313 , @xmath286 , @xmath318 , @xmath299 .
the critical value of @xmath319 .
, title="fig:",width=181 ]    finally let us mention the case of sixth degree potentials .
we consider @xmath320 . in this case
@xmath41 has a minimum at @xmath321 and two maxima in the antigravity domain ( see the left picture of fig .
[ v6ksi086 ] ) . both the minimum of @xmath41 at @xmath76 and a maxima of @xmath41 at @xmath322 corresponds to a stable de sitter solutions that are presented in fig .  [ v6ksi086 ] . in fig .
[ v6ksi086 ] we see oscillations in behavior the hubble parameter .
solutions that do not leave the domain @xmath76 are similar to the case of the fourth order potential .
there is the difference between solutions that cross the line @xmath180 .
namely , such a solution tends to a de sitter solution instead of the infinity point .
in this paper we have considered a class of spatially flat flrw models with bounce solutions , generalizing those found in @xcite .
we show that despite the model with the potential @xmath323 , studied in   @xcite is the unique model with the curvature scalar @xmath5 being the integral of motion , a much wider class of potentials admits bounce solutions followed by smooth future asymptotic behavior .
one class of models studied in the present paper generalizes the above mentioned potential by adding a term quadratic in @xmath8 .
this generalization allows us to find bounce solutions with non - monotonic behavior of the hubble parameter , in particular we found solutions with essential decreasing hubble parameter , including the solutions with the hubble parameter tending to zero .
moreover , we present a list of possible types of the bounce solutions and show that bounce behavior is still present even in the case when @xmath324 being negative , though in this case the nice future asymptotic is evidently lost .    another generalization replaces the quartic term in the potential by other even index .
we studied the quartic and sixth - power cases , changing also the coupling constant @xmath18 .
we show that despite the model becomes non - integrable , the qualitative features of solution found in  @xcite are still valid for a not so restricted conditions on the potential . on the other hand , some new interesting features like possible branch changing during the cosmological evolution appears , which can decrease the measure of initial conditions good for bounce , though
does not eliminate the possibility for the bounce completely . on the other hand , for sufficient small positive @xmath18
the bounce solution that tends to de sitter point does not exist .
note that inflationary models with the ricci scalar multiplied by a function of the scalar field are very popular  @xcite .
they not only do not contradict the recent observational data  @xcite , but also connect cosmology and particle physics .
note that models with non - minimally coupled scalar fields are quite natural because quantum corrections to the effective action with minimal coupling include induced gravity terms  @xcite .
also , the einstein frame model that corresponds to the bouncing model  @xcite describes inflation  @xcite ( see also  @xcite ) . note that non - local model  @xcite and its generalization  @xcite not only has bounce solution , but also ( at other values of parameters ) can describe inflation  @xcite .
the observation data  @xcite are consistent with ekpyrotic cyclic models in which the universe is smoothed and flattened during a period of slow contraction followed by a bounce ( see  @xcite and references therein ) .
it would be very interesting to construct cosmological model with a non - minimal coupling standard scalar field , a bounce solution of which is suitable for inflationary scenario .
the present paper can be considered as a step in this direction .
we plan to consider a possibility to construct such a solution in future publications .
research of e.p . is supported in part by grant mk-7835.2016.2 of the president of russian federation .
is supported by rfbr grant 14 - 02 - 00894 and by the russian government program of competitive growth of kazan federal university .
research of s.v . is supported in part by grant nsh-7989.2016.2 of the president of russian federation .
researches of e.p . and
s.v . are supported in part by the rfbr grant 14 - 01 - 00707 .
m.  novello and s.e.p .
bergliaffa , _ bouncing cosmologies _ , phys
.  rept .
* 463 * ( 2008 ) 127 ( arxiv:0802.1634 ) d.  battefeld and p.  peter , _ a critical review of classical bouncing cosmologies _ ,
* 571 * ( 2015 ) 1 ( arxiv:1406.2790 ) r.  brandenberger and p.  peter , _ bouncing cosmologies : progress and problems _ , arxiv:1603.05834 [ hep - th ] .
m.  koehn , j.  l.  lehners and b.  a.  ovrut , _ cosmological super - bounce , _ phys .
d * 90 * ( 2014 ) 025005 ( arxiv:1310.7577 ) m.  libanov , s.  mironov , and v.  rubakov , _ generalized galileons : instabilities of bouncing and genesis cosmologies and modified genesis _ , arxiv:1605.05992 [ hep - th ] .
a.  ijjas and p.j .
steinhardt , _ classically stable non - singular cosmological bounces _ , arxiv:1606.08880 [ gr - qc ] .
k.  bamba , a.n .
makarenko , a.n .
myagky , s.  nojiri , and s.d .
odintsov , _ bounce cosmology from @xmath325 gravity and @xmath325 bigravity _ , jcap * 1401 * ( 2014 ) 008 ( arxiv:1309.3748 ) s.  nojiri , s.d .
odintsov , and v.k .
oikonomou , _ bounce universe history from unimodular @xmath325 gravity _ , phys .
d * 93 * ( 2016 ) 084050 ( arxiv:1601.04112 ) .
tretyakova , a.a .
shatskiy , i.d .
novikov and s.  alexeyev , _ non - singular brans - dicke cosmology with cosmological constant _ , phys .
d * 85 * ( 2012 ) 124059 ( arxiv:1112.3770 ) b.  boisseau , h.  giacomini , d.  polarski , and a.a .
starobinsky , _ bouncing universes in scalar - tensor gravity models admitting negative potentials _ , j. cosmol .
* 1507 * ( 2015 ) 002 ( arxiv:1504.07927 )    a.yu .
kamenshchik , e.o .
pozdeeva , a.  tronconi , g.  venturi and s.yu .
vernov , _ interdependence between integrable cosmological models with minimal and non - minimal coupling _ , class . quantum grav . * 33 * ( 2016 ) 015004 ( arxiv:1509.00590 )            m.p .
dabrowski , t.  stachowiak and m.  szydlowski , _ phantom cosmologies _ ,
d * 68 * ( 2003 ) 103519 ( hep - th/0307128 ) i.ya .
arefeva , _ nonlocal string tachyon as a model for cosmological dark energy _ , aip conf .
* 826 * ( 2006 ) 301 ( astro - ph/0410443 ) ; + i.ya .
arefeva , a.s .
koshelev , and s.yu .
vernov , _ exact solution in a string cosmological model _ ,
* 148 * ( 2006 ) 895 ( arxiv : astro - ph/0412619 ) ; + m.l .
pucheu , f.a.p .
alves - junior , a.b .
barreto and c.  romero , _ cosmological models in weyl geometrical scalar - tensor theory _ , arxiv:1602.06966 [ gr - qc ] .
yi - fu cai , e.n .
saridakis , m.r .
setare , and jun - qing xia , _
quintom cosmology : theoretical implications and observations _ , phys .
rept . * 493 * ( 2010 ) 160 , arxiv:0909.2776 ; + zong - kuan guo , yun - song piao , xinmin zhang , yuan - zhong zhang , _ cosmological evolution of a quintom model of dark energy _ ,
b * 608 * ( 2005 ) 177182 ( astro - ph/0410654 ) ; + r.  lazkoz , g.  len , and i. quiros , _ quintom cosmologies with arbitrary potentials _ , phys . lett .
b * 649 * ( 2007 ) 103110 , astro - ph/0701353 ; + m.r .
setare , e.n .
saridakis , _ quintom cosmology with general potentials _ , int.j.mod.phys .
d * 18 * ( 2009 ) 549557 ( arxiv:0807.3807 ) ; + yi - fu cai , taotao qiu , yun - song piao , mingzhe li , xinmin zhang , _ bouncing universe with quintom matter _ , jhep * 0710 * ( 2007 ) 071 ( arxiv:0704.1090 )    i.ya .
arefeva , a.s .
koshelev , and s.yu .
vernov , _ crossing the @xmath326 barrier in the d3-brane dark energy model _ , phys .
d * 72 * ( 2005 ) 064017 ( arxiv : astro - ph/0507067 ) ; + s.yu .
vernov , _ construction of exact solutions in two - field cosmological models _ , theor .
* 155 * ( 2008 ) 544 ( arxiv : astro - ph/0612487 ) ; + i.ya .
arefeva , n.v .
bulatov , s.yu .
vernov , _ stable exact solutions in cosmological models with two scalar fields _ , theor .
* 163 * ( 2010 ) 788 ( arxiv:0911.5105 )      a.yu .
kamenshchik , a.  tronconi , g.  venturi , and s.yu .
vernov , _ reconstruction of scalar potentials in modified gravity models _ , phys . rev .
d * 87 * ( 2013 ) 063503 ( arxiv:1211.6272 ) ; + e.o .
pozdeeva and s.yu .
vernov , _ stable exact cosmological solutions in induced gravity models _ // aip conf
* 1606 * ( 2014 ) 4858 ( arxiv:1401.7550 )          i.ya .
arefeva , n.v .
bulatov , r.v .
gorbachev , s.yu .
vernov , _ non - minimally coupled cosmological models with the higgs - like potentials and negative cosmological constant _ , class .
quantum grav .
* 31 * ( 2014 ) 065007 ( arxiv:1206.2801 )      m.  sami , m ,  shahalam , m.  skugoreva and a.  toporensky , _ cosmological dynamics of a nonminimally coupled scalar field system and its late time cosmic relevance _
d * 86 * ( 2012 ) 103532 ( arxiv:1207.6691 )          b.l .
spokoiny , _ inflation and generation of perturbations in broken symmetric theory of gravity _ , phys .
b * 147 * ( 1984 ) 3943 ; + t.  futamase , k .-
chaotic inflationary scenario in models having nonminimal coupling with curvature _ , phys .
d * 39 * ( 1989 ) 399404 ; + r.  fakir , w.g .
unruh , _ improvement on cosmological chaotic inflation through nonminimal coupling , _ phys .
d * 41 * ( 1990 ) 17831791 ; + m.v .
libanov , v.a .
rubakov , p.g .
tinyakov , _ cosmology with nonminimal scalar field : graceful entrance into inflation _ ,
b * 442 * ( 1998 ) 63 ( arxiv : hep - ph/9807553 ) ; + a.  cerioni , f.  finelli , a.  tronconi , g.  venturi , _ inflation and reheating in induced gravity _ , phys .
b * 681 * ( 2009 ) 383386 ( arxiv:0906.1902 ) ;   r.  kallosh , a.  linde , d.  roest , _ the double attractor behavior of induced inflation _ , j. high energy phys . * 1409 * ( 2014 ) 062 ( arxiv:1407.4471 ) ; + m.  rinaldi , l.  vanzo , s.  zerbini , g.  venturi , _ inflationary quasi - scale invariant attractors _ , phys . rev .
d * 93 * ( 2016 ) 024040 ( arxiv:1505.03386 )    f.l .
bezrukov , m.  shaposhnikov , _ the standard model higgs boson as the inflaton , _ phys .
b * 659 * ( 2008 ) 703 ( arxiv:0710.3755 ) ; + a.o .
barvinsky , a.y .
kamenshchik , a.a .
starobinsky , _ inflation scenario via the standard model higgs boson and lhc _ , j. cosmol .
. phys . * 0811 * ( 2008 ) 021 ( arxiv:0809.2104 ) ; + f.  bezrukov , d.  gorbunov , m.  shaposhnikov , _ on initial conditions for the hot big bang _ , j. cosmol .
* 0906 * ( 2009 ) 029 ( arxiv:0812.3622 ) ; + a.o .
barvinsky , a.y .
kamenshchik , c.  kiefer , a.a .
starobinsky , c.f .
steinwachs , _ asymptotic freedom in inflationary cosmology with a nonminimally coupled higgs field _
, j. cosmol .
phys . * 0912 * ( 2009 ) 003 ( arxiv:0904.1698 ) ; + a.  de simone , m.p .
hertzberg and f.  wilczek , _ running inflation in the standard model _ , phys .
b * 678 * ( 2009 ) 1 ( arxiv:0812.4946 ) ; + f.l .
bezrukov , a.  magnin , m.  shaposhnikov , s.  sibiryakov , _
higgs inflation : consistency and generalisations _ , j. high energy phys . * 1101 * ( 2011 ) 016 ( arxiv:1008.5157 ) ; + a.o .
barvinsky , a.yu .
kamenshchik , c.  kiefer , a.a .  starobinsky , c.f .
steinwachs , _ higgs boson , renormalization group , and cosmology _ , eur .
j. c * 72 * ( 2012 ) 2219 ( arxiv:0910.1041 ) ; + f.  bezrukov , _ the higgs field as an inflaton _ , class .
quant . grav .
* 30 * ( 2013 ) 214001 ( arxiv:1307.0708 ) ; + j.  ren , z .- z .
xianyu , h.j .
higgs gravitational interaction , weak boson scattering , and higgs inflation in jordan and einstein frames _ , j. cosmol . astropart
* 1406 * ( 2014 ) 032 ( arxiv:1404.4627 )    d.i .
kaiser , _ induced - gravity inflation and the density perturbation spectrum _ , phys .
b * 340 * ( 1994 ) 2328 ( arxiv : astro - ph/9405029 ) ; + d.i .
kaiser , _ primordial spectral indices from generalized einstein theories _ , phys .
d * 52 * ( 1995 ) 4295 ( arxiv : astro - ph/9408044 ) r.  kallosh and a.  linde , _ universality class in conformal inflation _ , j. cosmol .
* 1307 * ( 2013 ) 002 ( arxiv:1306.5220 ) .
k.  bamba , sh .
nojiri and s.d .
odintsov , _ reconstruction of scalar field theories realizing inflation consistent with the planck and bicep2 results _ , phys
b * 737 * ( 2014 ) 374378 ( arxiv:1406.2417 )    e. elizalde , s.d .
odintsov , e.o .
pozdeeva , and s.yu .
vernov , _ renormalization - group inflationary scalar electrodynamics and @xmath327 scenarios confronted with planck2013 and bicep2 results _ , phys .
d * 90 * ( 2014 ) 084001 ( arxiv:1408.1285 ) ; + e.  elizalde , s.d .
odintsov , e.o .
pozdeeva , s.yu .
vernov , _ cosmological attractor inflation from the rg - improved higgs sector of finite gauge theory _ , j. cosmol .
* 1602 * ( 2016 ) 025 ( arxiv:1509.08817 )        n.a .
chernikov , e.a .
tagirov , _ quantum theory of scalar fields in de sitter space - time _ , annales poincare phys .
theor .  a * 9 * ( 1968 ) 109 ; e.a .
tagirov , _ consequences of field quantization in de sitter type cosmological models _ , annals phys . * 76 * ( 1973 )  561 c.g .
callan , s.r .
coleman , r.  jackiw , _ a new improved energy  momentum tensor _ , annals phys . * 59 * ( 1970 ) 42