inflation in the early universe @xcite is now an indispensable ingredient of modern cosmology not only to explain the global properties of homogeneous and isotropic space with a vanishingly small spatial curvature but also to account for the origin of the primordial curvature perturbation that seeded cosmic structure formation @xcite . at present , despite the significant progress in the state - of - the - art precise measurements of the cosmic microwave background radiation ( cmb ) by wmap  @xcite and planck  @xcite missions , there is no single observational result in conflict with the single - field inflation paradigm @xcite . in particular , the anti - correlation of the temperature and the e - mode polarization anisotropies on large scales observed by the wmap mission
strongly supports the superhorizon perturbations suggested by inflation @xcite .    in other words , once inflation sets in , virtually all the available cosmological observation data can be explained simultaneously irrespective of the initial condition of the universe
this does not mean that we may be indifferent to the initial condition of the universe before inflation .
on the contrary , in order to achieve complete understanding of the cosmic history , we must work out the very beginning of the universe that may smoothly evolve into the inflationary phase .    as is well known , as long as the null energy condition ( nec ) is satisfied in the expanding phase , the hubble parameter and the energy density of the universe increase backward in cosmic time .
so , it is often claimed that , if one tries to discuss what happened before inflation and/or how inflation started , one needs to know the information of very high energy physics , and challenge the initial singularity problem @xcite in terms of quantum gravity .
but , this is not always the case .
recently , it was recognized that , if an action includes higher derivative terms of a scalar field like the galileon terms , the nec can be violated without encountering ghost nor gradient instabilities .
see , e.g. , ref .
@xcite for a recent review and ref .
@xcite for a subtle issue of nonlinear instabilities .
if the nec is violated , the energy density can grow as time proceeds , contrary to the conventional wisdom . in the nec violating theories , the universe can therefore start from the static zero - energy state described by the minkowski spacetime from infinite past  @xcite , and the universe starts expansion with the increase of the energy density .
such a picture of the emergence of the universe was first proposed by creminelli _
et al . _
@xcite with the name galilean genesis . in their model , however , the hot big bang state was postulated to be realized after the effective field theory description breaks down as the energy density blows up beyond its realm of validity .
therefore , the theory to describe the most important epoch of the early universe is lacking there .
nevertheless , since their original idea is so interesting that a number of extension has been made in a wider class of scalar field theories  @xcite and various aspects of the genesis scenario have been explored in the literature @xcite , such as avoidance of the superluminal propagation of perturbations and absence of primordial tensor perturbations .
they have been unsuccessful , however , to realize transition from the genesis phase to the hot big bang state within their model lagrangians .    in this paper , we take a different approach , namely , to make use of the galilean genesis to explain the initial condition of the universe before inflation and smoothly connect it to the inflationary phase , thereby solving the initial singularity problem @xcite and the trans - planckian problem @xcite ( see also @xcite ) in inflationary cosmology .
in fact , such an approach has also been put forward by pirtskhalava _
et al . _
@xcite recently .
their model lagrangian , however , gives rise to gradient instability as it is , although it has been argued there that higher - order structure of the effective field theory for perturbations possesses enough freedom to cure the gradient instability .
discussion on termination of inflation and reheating is not presernted there , either .    in the present paper ,
we construct a specific model free from any catastrophic instabilities and with subluminal velocities of primordial perturbations . in our setup
the universe starts from the minkowski spacetime from infinite past and is smoothly connected to the inflationary phase followed by the graceful exit . for this purpose , we provide a generic lagrangian capable of describing the background and perturbation evolution in all the above phases instead of choosing the effective field theory approach because the latter can not capture the evolution of the background and perturbations from pre - inflationary genesis to the exit from inflation with the same single lagrangian .
although we start with asymptotically minkowski space at the past infinity for aesthetic beauty , it has been shown that the galilean genesis solution is an attractor for a variety of initial conditions including those with a negative hubble parameter and/or finite curvature , provided that the time derivative of the scalar field has the right sign @xcite .
the horndeski theory  @xcite or the generalized galileon @xcite , whose mutual equivalence was first shown in @xcite , is known to be the most general scalar - tensor tensor theory with the second - order field equations , and thereby avoid ostrogradski instabilities in spite of having higher derivative terms in the action .
the theory can be generalized to have second - order field equations only in a specific gauge while maintaining the number of propagating degrees of freedom .
this possibility was realized recently by gleyzes _
et al . _
@xcite ( see also ref .
@xcite ) and was extended further by gao  @xcite .
the number of propagating degrees of freedom in these theories is indeed shown to be the same as that of the horndeski theory  @xcite . in this paper , we use the subclass of gao s framework as a concrete realization of the unified scenario starting from galilean genesis through inflation to the graceful exit .
this paper is organized as follows . in the next section ,
we give a framework of our model and derive the background equations of motion and the quadratic actions of cosmological perturbations . in sec .
iii , a concrete lagrangian is constructed to describe our scenario beginning from the genesis phase through the inflationary one to the graceful exit , and such a background dynamics is presented explicitly . in sec .
iv , we discuss the stability during each phase based on the quadratic actions of cosmological perturbations . in sec .
v , a concrete realization of our scenario is given .
the final section is devoted to our conclusions and discussion .
let us start with describing the general framework to construct and study our explicit realization of the early - time completion of inflation .
we would like to consider theories composed of a metric @xmath0 and a single scalar field @xmath1 , and hence it will be appropriate to work in the horndeski theory .
the lagrangian of the horndeski theory is of the form @xmath2,\label{horndeski - l}\end{aligned}\ ] ] where @xmath3 , @xmath4 is the four - dimensional ricci scalar , and @xmath5 is the four - dimensional einstein tensor .
we have four arbitrary functions of @xmath1 and @xmath6 in the horndeski theory .
this is the most general lagrangian having second - order field equations .
nevertheless , it will turn out that this framework is insufficient for our purpose , and hence we have to go beyond the horndeski theory .
one can generalize the horndeski theory to possess higher order field equations while maintaining the number of propagating degrees of freedom  @xcite .
the first step to do so is to perform an adm decomposition by taking @xmath7 const hypersurfaces as constant time hypersurfaces . in the adm language ,
the metric is written as @xmath8 by definition @xmath1 is a function of only @xmath9 , @xmath10 , and @xmath11 , where a dot denotes differentiation with respect to @xmath9 , so any function of @xmath1 and @xmath6 can be regarded as a function of @xmath9 and the lapse function @xmath12 , provided that @xmath13 and @xmath14 never vanish .
then , the horndeski lagrangian  ( [ horndeski - l ] ) can be written in terms of the adm variables as @xmath15 with @xmath16 where @xmath17 and @xmath18 are the extrinsic and intrinsic curvature tensors on the constant time hypersurfaces , and @xmath19 , @xmath20 , @xmath21 , and @xmath22 are subject to the relations @xmath23 variation of the above lagrangian with respect to @xmath12 gives a second - class constraint that eliminates only one degree of freedom , as opposed to general relativity .
the key trick to generalize the horndeski theory is to notice that this property remains the same even if one liberates @xmath19 and @xmath20 from the restriction imposed by eq .
( [ ab - constraint ] )  @xcite .
we thus arrive at the so called glpv theory that is more general than horndeski but has the same number of propagating degrees of freedom .
one can move back to a covariant form of the lagrangian by introducing the unit normal to the constant time hypersurfaces as @xmath24 , writing the extrinsic curvature tensor in terms of @xmath25 , and using the gauss - codazzi equations .
since there are six arbitrary functions of @xmath9 and @xmath12 in the adm form , the resultant covariant lagrangian has six arbitrary functions of @xmath1 and @xmath6 .    the above idea has been pushed forward by gao  @xcite , who proposed a unified framework to study single scalar - tensor theories beyond horndeski .
one can write a general lagrangian in the adm form as @xmath26,\end{aligned}\ ] ] where the coefficients @xmath27 , @xmath28 , ... are arbitrary functions of @xmath9 and @xmath12 .
the hamiltonian depends nonlinearly on @xmath12 as in the glpv theory , giving rise to a single scalar degree of freedom on top of the traceless and transverse gravitons  @xcite .    in this paper
, we will employ the lagrangian @xmath29 with @xmath30 where @xmath31 , @xmath32 , and @xmath33 are constant parameters of the theory .
this is a deformation of the glpv lagrangian and belongs to a subclass of gao s framework .
the generalization to this level is sufficient for the purpose of the present paper .
the glpv theory is recovered by taking @xmath34 .
given the lagrangian  ( [ ourlag ] ) in the adm form , one can restore the scalar degree of freedom @xmath1 to write its covariant expression in the same way as in the glpv theory .
however , it will be more convenient for our purpose to use the explicitly time - dependent lagrangian , because by doing so one can easily design the lagrangian so as to admit the desired cosmological evolution .    before specifying the suitable form of @xmath35 , @xmath36 , ... to construct our early universe model , let us derive the general equations governing the background and perturbation dynamics of cosmologies based on the lagrangian  ( [ ourlag ] ) .
the adm variables are given by @xmath37 where @xmath38 is the curvature perturbation in the unitary gauge and @xmath39 is the transverse and traceless tensor perturbation .
a spatially flat background has been assumed and the spatial diffeomorphism invariance was used to write @xmath40 in the above form . in the following , the background value of the lapse function is denoted by @xmath12 where there is no worry about confusion .
substituting eq .
( [ adm - cosmological ] ) to the lagrangian  ( [ ourlag ] ) , we obtain the background part of the lagrangian as @xmath41 where @xmath42 , @xmath43 , and @xmath44 . at the background level , @xmath31 , @xmath32 , and @xmath33 just rescale @xmath19 and @xmath20 . in what follows
we simply consider the case with @xmath45 .
since we are considering a spatially flat universe , we have @xmath46 at zeroth order , and hence @xmath21 and @xmath22 play no role in the background dynamics .
varying eq .
( [ zero - bg ] ) with respect to @xmath12 and @xmath47 , we obtain , respectively , @xmath48 where a prime represents differentiation with respect to @xmath12 .
the background equations contain at most second derivatives of the scale factor and first derivatives of the lapse function .
the quadratic lagrangian for the tensor perturbation is given by @xmath49,\end{aligned}\ ] ] where @xmath50 the equation of motion contains at most second derivatives both in time and space .
the tensor perturbation is stable provided that @xmath51 and @xmath52 .
the quadratic lagrangian for the scalar perturbations is given by @xmath53,\label{lags2}\end{aligned}\ ] ] where the coefficients are defined as @xmath54 and note the relation @xmath55 .
one has @xmath56 in the horndeski and glpv theories , in which @xmath34 .
therefore , the last term in the lagrangian  ( [ lags2 ] ) is the novel consequence of theories beyond glpv .    from @xmath57 and @xmath58
we obtain @xmath59,\label{eldeltan } \\
\frac{\partial^2\chi}{a^2 } & = & \frac{1}{\theta^2+\sigma{\cal c}}\left [ ( 3\theta^2+\sigma{\cal g}_b)\frac{\dot\zeta}{n } -\theta{\cal g}_b\frac{\partial^2\zeta}{a^2 } \right].\label{elchi}\end{aligned}\ ] ] substituting eqs .  ( [ eldeltan ] ) and  ( [ elchi ] ) into eq .
( [ lags2 ] ) , we obtain the reduced lagrangian for the curvature perturbation , @xmath60,\end{aligned}\ ] ] where @xmath61 thus , if @xmath62 , the equation of motion for @xmath38 has the fourth derivative in space , giving the dispersion relation @xmath63 we require that @xmath64 in order to avoid ghost instabilities .
however , we allow for a negative sound speed squared , @xmath65 , for a short period of time . in the absence of the @xmath66 term ( @xmath56 ) , a negative sound speed squared would cause a rapid growth of instabilities for large @xmath67 modes . in this paper , we consider theories with @xmath62 , so that the curvature perturbation with large @xmath67 can be stabilized by requiring that @xmath68 .    as will be seen in the rest of the paper ,
the sound speed squared becomes negative at the transition from one phase to another .
such a behavior should not occur even for a tiny period because high wavenumber modes would grow exponentially rapidly .
however , we could not avoid it not only within the horndesky theory but also the glpv theory despite we analyzed extensive models . on the other hand ,
we have not been successful in proving that this is an inevitable consequence . since our primary purpose is to show an existence proof of the model to realize our intended cosmic evolution without any instabilities , we construct a specific model by going beyond the glpv theory and invoking the @xmath66 term .
the lagrangian we study in this paper is characterized by a single time - dependent function @xmath69 and four functions @xmath70 of @xmath12 : @xmath71 where @xmath72 is a constant parameter .
we have introduced the mass scales @xmath73 ( and the planck mass @xmath74 ) , so that @xmath69 and @xmath75 are dimensionless .
the other two functions , @xmath21 and @xmath22 , are arbitrary at this stage because they have no impacts on the background dynamics .
note that @xmath76 is not a dynamical variable .
specifying the functions @xmath77 and @xmath78 amounts to defining a concrete theory .
the above forms of @xmath79 are chosen so that the theory admits an inflationary universe preceded by the generalized galilean genesis while retaining much of the generality .
other choices could be possible and hence we do not claim that this is the most general description of such scenarios at all . instead , as we mentioned above , we would provide the existence proof of desired models by demonstrating that a sufficiently wide class of healthy models can indeed be constructed .
we design @xmath69 so as to implement the ( generalized ) galilean genesis followed by inflation and a graceful exit from the prolonged inflationary phase .
our choice is @xmath80 well before @xmath81 , and @xmath82 for @xmath83 .
as our time variable starts at @xmath84 with asymptotically minkowski spacetime configuration , @xmath9 is large and negative in the beginning , so we find @xmath85 in eq .
( [ f : gen ] ) . as will be seen shortly , the initial stage described by eq .
( [ f : gen ] ) corresponds to the generalized galilean genesis , while the subsequent stage described by eq .
( [ f : ds ] ) to inflation .
after a sufficiently long period of the inflationary stage , we assume that @xmath86 for @xmath87 , where @xmath88 is the time at the end of inflation . with this the universe exits from inflation .
in what follows we will investigate the background evolution of each stage .
assuming that @xmath89 in the first stage where @xmath76 is given by eq .
( [ f : gen ] ) , let us look for a consistent solution for large @xmath76 .
the background field equations read @xmath90 it can be seen from eq .
( [ fr : gen ] ) that the lapse function @xmath12 is a constant , @xmath91 , satisfying @xmath92 then , @xmath93 is consistently determined from eq .  (
[ pres : gen ] ) , which can be written as @xmath94 where @xmath95 is a constant .
this leads to the generalized galilean genesis solution  @xcite : @xmath96 it is required that @xmath97 to guarantee @xmath98 .
we have thus arrived at the generalized galilean genesis solution starting from the lagrangian written in the adm form rather than in the covariant form .
the original galilean genesis solution found in ref .
@xcite corresponds to @xmath99 .    in deriving the above solution , @xmath100 and @xmath101
are always subdominant due to the assumed scalings @xmath102 and @xmath103 .
therefore , any choices of @xmath104 and @xmath105 will not spoil the above galilean genesis solution . as will be seen in the next section ,
those two terms are also irrelevant to the stability conditions during the genesis phase .
the galilean genesis phase will end at @xmath106 since the function @xmath76 is constant for @xmath83 .
in the subsequent phase we obtain the de sitter solution , @xmath107 const and @xmath108 const , satisfying @xmath109 ( note that @xmath79 is now a function of @xmath12 only and is independent of @xmath9 . )
a @xmath9-independent lagrangian in the adm form can be recast in a covariant lagrangian with the shift symmetry , @xmath110 .
this implies that the above exact de sitter solution corresponds to kinetically driven g - inflation .
if one invokes a weak time - dependence in @xmath76 , one obtains quasi - de sitter inflation instead .
after the prolonged phase of inflation , @xmath76 is given by eq .
( [ f : exit ] ) .
we assume that @xmath9 is sufficiently large , so that @xmath85
. then , we have a consistent solution with @xmath111 const and @xmath112 satisfying @xmath113 thus , one can implement a graceful exit from inflation .
it follows from eq .
( [ exit : eq1 ] ) that @xmath114 it can be shown using eqs .
( [ exit : eq1 ] ) and  ( [ exit : eq2 ] ) that , during this third stage , @xmath115 it is therefore necessary to impose @xmath116 .    in the standard potential - driven inflation models @xcite inflation
is followed by coherent field oscillation of the inflaton scalar field which decays to radiation to reheat the universe . in the present approach the scalar field @xmath1
is used to specify constant time hypersurfaces , so that @xmath13 may not vanish in order to preserve one - to - one correspondence between @xmath1 and the cosmic time @xmath9 .
hence one must switch from the adm language we used to construct the action to the conventional `` @xmath1 language '' at this point in order to apply the standard reheating mechanism , which is all right but looks like sewing a fox s skin to the lion s .    here
instead we consider another reheating mechanism which can take place without breaking the one - to - one correspondence between @xmath1 and @xmath9 , namely , the gravitational reheating due to the change of geometry or the cosmic expansion law  @xcite .    during the transition from the de sitter inflation to a decelerated power - law expansion , conformally non - invariant particles are produced with the initial energy density @xmath117 where @xmath118 is a factor determined by the effective number of conformally noninvariant fields and the change of the geometry .
for example , for @xmath119 or 4 , a single minimally coupled massless scalar field contributes to @xmath118 by @xmath120 respectively @xcite . here
@xmath121 is the time required for the transition . in case
it is nonminimally coupled with a coupling parameter @xmath122 , a factor @xmath123 is multiplied there .    in order for the radiation thus created to dominate the universe , the energy density of the scalar field must dissipate more rapidly , namely , @xmath124 then , the reheating temperature at the radiation domination is given by @xmath125 where @xmath126 is the effective number of relativistic degrees of freedom and we have assumed the universe would evolve in the same way as in the einstein gravity after inflation . if long - lived massive particles are copiously produced at the gravitational particle production , the reheating temperature may be significantly higher then the above value . furthermore , the decay of quasi - flat direction may produce a large amount of entropy to reheat the universe efficiently and create matter particles  @xcite .    at this point one may wonder if standard particle production due to the decay of the inflaton might work in our scenario . in such a reheating process , the inflaton decay ( and the associated energy loss of the inflaton ) must be very efficient in order to avoid the dominance of the inflaton energy density in comparison to that of produced radiation . for this reason the inflaton must oscillate after inflation .
however , in our construction of the action in the adm language , the field @xmath1 is used to specify constant time hypersurfaces and so @xmath1 and the cosmic time @xmath9 have a one - to - one correspondence .
therefore , our approach can not describe the oscillatory phase of the inflaton and we need to resort to gravitational reheating instead of the standard reheating scenario due to the inflaton decay .
we expect that a different action based on the @xmath1 language instead of the adm language can accommodate the oscillatory phase as well as the genesis and inflationary phases .
minimum modification would be such that one moves to the @xmath1 language at the end of inflation to construct a @xmath1-lagrangian that admits the standard reheating scenario after inflation .
however , this is beyond the scope of the present paper and is left for future work .
having obtained the background evolution of our scenario , let us investigate the nature of primordial perturbations and stability , using the result of the generic analysis in sec .
[ sec : cp ] .      during the genesis phase , we have @xmath127 obviously , the kinetic term of the tensor perturbations has the right sign , @xmath51 . for large @xmath76
, we see @xmath128 ( as long as @xmath62 ) , and hence @xmath129 this implies that @xmath130 const , while @xmath131 . the kinetic term of the curvature perturbation has the right sign if @xmath132 thus , it is sufficient to impose @xmath133 ( we are considering only the case with @xmath134 . ) another stability condition , @xmath135 , is equivalent to requiring that @xmath136    since @xmath137 depends on @xmath138 and @xmath139 and these two functions are irrelevant to the background dynamics , the condition @xmath140 can easily be satisfied without spoiling the genesis background .
suppose for simplicity that @xmath141 where @xmath142 is a constant .
then , @xmath143 .
for the scalar perturbations we have @xmath144 \nonumber\\ & = & { \rm const.}\end{aligned}\ ] ] this can also be made positive by an appropriate choice of @xmath145 .
it should be noted that if @xmath146 then we inevitably have @xmath147 ; the @xmath148 term is crucial for the stable violation of the nec .
note also that , if we take sufficiently small @xmath149 , the sound speed @xmath150 can be smaller than unity , which applies also to the other two phases discussed below .
let us move to discuss the nature of the primordial fluctuations in the genesis phase .
since @xmath151 const , the tensor perturbations behave in the same way as in the minkowski spacetime . therefore ,
no large tensor modes are generated during the first stage of our scenario .
plane with @xmath152 decreasing toward the right . in the region below ( above ) the red broken curve , @xmath153 is dominated by the term proportional to @xmath66 @xmath154 .
modes with @xmath155 experience the break down of the wkb approximation around the point crossing the blue solid curve beyond which @xmath38 is frozen , while modes with @xmath156 do not .
, height=415 ]    the behavior of the curvature perturbation turns out to be more nontrivial , as sketched in fig .
[ fig : horizon.eps ] .
recalling that @xmath157 const , @xmath158 const , and @xmath159 , the equation of motion for @xmath38 in the fourier space is of the form @xmath160 where @xmath161 and @xmath162 with @xmath163 and @xmath164 being some constants .
for sufficiently large @xmath152 , we have @xmath165 .
one may define the time at which this approximation breaks down as @xmath166 , and for @xmath167 we have @xmath168 .    with some manipulation , it is found that @xmath169 where @xmath170 .
this implies that for the modes with @xmath171 the wkb approximation is always good in the genesis phase , @xmath172 giving @xmath173 for @xmath167 .
thus , the amplitude of those modes at late times in the genesis phase is given by @xmath174    for the modes with @xmath175 , the wkb approximation breaks down at some time and then the curvature perturbation freezes .
this `` horizon crossing '' occurs at @xmath176 .
it can be seen that @xmath177 for @xmath175 , for the @xmath178 mode .
the genesis phase could end sufficiently early so that @xmath179 .
if this is the case , we only need to care about the modes with @xmath175 . ] which allows us to study the freezing process by using the solution to eq .
( [ zetaeqgen ] ) with @xmath180 .
the exact solution in this case that matches the positive frequency wkb solution for @xmath181 is given by @xmath182 where @xmath183 is the hankel function of the first kind .
the frozen amplitude can thus be evaluated by taking the limit @xmath184 in the solution  ( [ smallksol ] ) , leading to @xmath185 for @xmath186 , @xmath187 is dominated by the @xmath188 term where the solution  ( [ smallksol ] ) is no longer exact .
the frozen amplitude  ( [ frozen ] ) , however , is still valid even in this regime since the solution to eq .
( [ zetaeqgen ] ) with the effective frequency  ( [ omega ] ) does not oscillate any more and remains constant .
hence , the expression of the power spectrum  ( [ frozen ] ) is correct for the entire range of @xmath189 .    to summarize , the power spectrum of the curvature perturbation generated during the genesis phase is blue and hence is suppressed on large scales .      in the ( de sitter )
inflationary phase , @xmath190 , @xmath191 , @xmath192 , @xmath193 , and @xmath194 are time - independent .
we require that all those coefficients are positive during inflation in order to avoid instabilities .
since the quadratic action for the tensor perturbations is essentially the same as that of generalized g - inflation , the power spectrum of the primordial tensor perturbations is given by  @xcite @xmath195    the equation of motion for the canonically normalized variable @xmath196 during inflation is of the form @xmath197 where @xmath198 with @xmath199 and @xmath200 being dimensionless constants . here
, we have introduced the conformal time @xmath201 defined by @xmath202 .
the dispersion relations of this form have been studied in the context of inflation , e.g. , in refs .
the positive frequency modes are given by  @xcite @xmath203 where @xmath204 is the whittaker function .
taking the limit @xmath205 , the power spectrum of the curvature perturbation can be calculated as @xmath206 where @xmath207 even in the presence of the @xmath66 term in the dispersion relation , the power spectrum is scale - invariant in the case of exact de sitter inflation . since we have @xmath208 as @xmath209 , we recover the result of generalized g - inflation  @xcite in the limit @xmath210 .
for @xmath211 we have @xmath212 , so that one can take the limit @xmath213 smoothly to get @xmath214    we have approximated the inflationary phase as exact de sitter .
if we consider a background slightly different from de sitter by incorporating weak time dependence in @xmath76 , we would be able to obtain a tilted spectrum of @xmath38 .
after inflation , we have @xmath215 , @xmath216 , @xmath217 where to simplify the expression we introduced @xmath218 recalling that we have been imposing @xmath219 , all of these coefficients are positive provided that @xmath220
. this condition can be written equivalently as @xmath221
and ( b ) the lapse function @xmath12 around the genesis - de sitter transition . , height=377 ]    , and ( b ) the coefficient of @xmath222 ( divided by @xmath223 ) around the genesis - de sitter transition .
, height=377 ]    let us provide a concrete lagrangian exhibiting the genesis - de sitter transition .
the lagrangian is characterized by @xmath224 where @xmath225 and @xmath226 are constants .
we take @xmath227 , @xmath228 , and @xmath229 .
we also take @xmath230 to guarantee the stability .
this corresponds to the ( @xmath219 generalization of the ) unitary gauge description of the lagrangian considered in ref .
@xcite . in the genesis stage
we have @xmath231<0.\end{aligned}\ ] ] since @xmath230 and @xmath232 , we see that @xmath233 and @xmath234 .
we also see that @xmath235 -1,\end{aligned}\ ] ] and hence it is easy to satisfy @xmath236 during the genesis phase by choosing the parameters appropriately .
a numerical example of the genesis - de sitter transition is illustrated in figs .
[ fig : hn.eps ] and  [ fig : stability.eps ] .
our numerical calculation was performed as follows : we solve the evolution equations @xmath237 and @xmath238 with initial data @xmath239 satisfying @xmath240 , and confirm that the constraint @xmath240 is satisfied at each time step . in the numerical calculation ,
the parameters are given by @xmath241 , @xmath242 , @xmath243 , @xmath244 , and @xmath245 .
the function @xmath69 is taken to be @xmath246+f_{1},\label{f : ex1}\end{aligned}\ ] ] with @xmath247 , @xmath248 , and @xmath249 . the background evolution is shown in fig .  [ fig : hn.eps ] .
the evolution of the sound speed squared , @xmath250 , and the coefficient of @xmath222 in the dispersion relation is shown in fig .
[ fig : stability.eps ] . as pointed out in ref .
@xcite , @xmath251 flips the sign at the transition .
the sound speed squared is positive except in this finite period . during the genesis and subsequent de sitter phases we have @xmath64 and @xmath135 , and
therefore we may conclude that this model is stable .     and ( b ) the lapse function @xmath12 around the end of inflation .
height=377 ]    , and ( b ) the coefficient of @xmath222 ( divided by @xmath223 ) around the end of inflation . , height=377 ]    although we have thus obtained the stable example of the genesis - de sitter transition , the simple example  ( [ example ] ) is not completely satisfactory if one would want successful gravitational reheating .
indeed , the condition  ( [ exit : c1 ] ) implies that @xmath252 , but @xmath253 for such @xmath254 .
this problem can be evaded easily by the following small deformation of @xmath255 : @xmath256 where @xmath257 is a parameter smaller than @xmath258 .
the condition  ( [ exit : c1 ] ) now reads @xmath259 , i.e. , @xmath260 , while @xmath261 is positive for @xmath262 .
the stability condition further restricts the allowed ranges of @xmath254 and @xmath257 .
the necessary condition for stability is @xmath263 [ see eq .
( [ ge : st ] ) ] .
this translates to @xmath264 , leading to @xmath265 = 4.8 .
note that the small deformation of @xmath255 with @xmath266 does not change the background and perturbation dynamics of the genesis and inflationary phases .    to illustrate the final stage of inflation , let us take @xmath267\right\}^{1/(\alpha + 1)},\label{f : ex2}\end{aligned}\ ] ] where the origin of time is shifted so that the end of inflation is given by @xmath268 .
in the numerical plots presented in figs .
[ fig : hn - ds - r.eps ] and  [ fig : stability - ds - r.eps ] , the parameters are given by @xmath269 , @xmath270 , and @xmath271 , while the other parameters are taken to be the same as the previous example of the genesis - de sitter transition .
it is found that @xmath272 .
again , we see that @xmath273 in the finite period around the transition .
however , @xmath192 and @xmath194 remain positive all through the inflation and subsequent stages .
in this paper , we have introduced a generic description of galilean genesis in terms of the adm lagrangian and constructed a concrete realization of inflation preceded by galilean genesis , _
i.e. _ , the scenario in which the universe starts from minkowski spacetime in the asymptotic past and is connected smoothly to the inflationary phase followed by the graceful exit .
our model utilizes the recent extension of the horndeski theory , which has the same number of propagating degrees of freedom as the horndeski theory and thus can avoid ostrogradski instabilities .
this approach allows us to cover the background and perturbation evolution in all the three phases with the same single lagrangian , as opposed to the effective field theory approach . in our scenario ,
the sound speed squared during the transition from the genesis phase to inflation becomes negative for a short period .
however , thanks to the nonlinear dispersion relation arising from the fourth - order derivative term in the quadratic action , modes with higher momenta are stable and the growth rate of perturbations with smaller momenta is finite and under control .
it should also be noted that the sound speed of the primordial perturbations can be smaller than unity by choosing the parameter of the model appropriately .
although we have constructed our inflation model in order to resolve the initial singularity and possible trans - planckian problems by incorporating galilean genesis phase before inflation , we could make use of our model to realize the original galilean genesis scenario , which is an alternative to inflationary cosmology , simply by taking vanishingly short period of inflation there .
as discussed in the appendix , the sound speed squared becomes negative at the transition also in this case , but the instabilities are relevant only for small @xmath67 modes thanks to the @xmath66 term in the dispersion relation .
thus , the transition from the genesis phase to the reheating stage is described in a healthy and controllable manner .
in fact , it would be fair to say that such a cosmology works quite well among the proposed alternatives to inflation , because , in contrast with the bouncing cosmology , in which all the would - be decaying modes in the expanding universe such as vector fluctuations and spatial anisotropy severely increase in an undesirable manner , the genesis solution is an attractor and generation of nearly scale - invariant curvature perturbation is also possible with an appropriate choice of model parameters @xcite .
since no first - order tensor perturbation is generated in this type of scenarios , detection of tensor perturbation with its amplitude larger than @xmath274 would be a smoking gun of inflation .
this work was supported in part by the jsps grant - in - aid for scientific research nos .
24740161 ( t.k . ) , 25287054 and 26610062 ( m.y . ) , 23340058 and 15h02082 ( j.y . ) .
and ( b ) the lapse function @xmath12 around the genesis - reheating transition .
, height=377 ]      in the main text , we consider the scenario in which galilean genesis is followed by inflation . in this appendix
, we will go back to the original motivation of galilean genesis and study how we can match smoothly the genesis phase to the reheating phase .
our approach based on the adm lagrangian is quite useful in analyzing such a situation as well .
it is now obvious that by taking @xmath275 and gluing the two functions smoothly at around @xmath276 , one can describe the genesis - reheating transition . as a concrete example , we glue @xmath277 and @xmath278 smoothly at around @xmath276 and perform a numerical calculation as shown in figs .  [ fig : hn - g - r.eps ] and  [ fig : stability - g - r.eps ] .
the other parameters are the same as those taken in the main text .
as is expected , the numerical result here is much the same as the case where a duration of the intermediate inflationary phase is taken to be very short . in particular , @xmath279 becomes negative at the genesis - reheating transition .
the model is nevertheless stable since the conditions @xmath64 and @xmath135 remain satisfied .
99 a.  a.  starobinsky , `` a new type of isotropic cosmological models without singularity , '' phys .
b * 91 * , 99 ( 1980 ) .
k.  sato , `` first order phase transition of a vacuum and expansion of the universe , '' mon .  not .
soc .   * 195 * , 467 ( 1981 ) .
a.  h.  guth , `` the inflationary universe : a possible solution to the horizon and flatness problems , '' phys .
d * 23 * , 347 ( 1981 ) .
a.  d.  linde , `` a new inflationary universe scenario : a possible solution of the horizon , flatness , homogeneity , isotropy and primordial monopole problems , '' phys .
b * 108 * , 389 ( 1982 ) . for
a review of inflation , see , e.g. , j.  yokoyama , `` inflation : 1980 - 201x , '' ptep * 2014 * , no .
6 , 06b103 ( 2014 ) .
v.  f.  mukhanov and g.  v.  chibisov , `` quantum fluctuation and nonsingular universe .
( in russian ) , '' jetp lett .   *
33 * , 532 ( 1981 ) [ pisma zh .
fiz .   * 33 * , 549 ( 1981 ) ] .
a.  a.  starobinsky , `` dynamics of phase transition in the new inflationary universe scenario and generation of perturbations , '' phys .
b * 117 * , 175 ( 1982 ) .
s.  w.  hawking , `` the development of irregularities in a single bubble inflationary universe , '' phys .
b * 115 * , 295 ( 1982 ) .
a.  h.  guth and s.  y.  pi , `` fluctuations in the new inflationary universe , '' phys .
lett .   * 49 * , 1110 ( 1982 ) .    c.  l.  bennett _ et al .
_ [ wmap collaboration ] , `` nine - year wilkinson microwave anisotropy probe ( wmap ) observations : final maps and results , '' astrophys .  j.  suppl .
* 208 * , 20 ( 2013 ) [ arxiv:1212.5225 [ astro-ph.co ] ] .
g.  hinshaw _ et al . _
[ wmap collaboration ] , `` nine - year wilkinson microwave anisotropy probe ( wmap ) observations : cosmological parameter results , '' astrophys .
j.  suppl .
* 208 * , 19 ( 2013 ) [ arxiv:1212.5226 [ astro-ph.co ] ] .
p.  a.  r.  ade _ et al .
_ [ planck collaboration ] , `` planck 2013 results .
i. overview of products and scientific results , '' astron .
astrophys .   * 571 * , a1 ( 2014 ) [ arxiv:1303.5062 [ astro-ph.co ] ] .
p.  a.  r.  ade _ et al .
_ [ planck collaboration ] , `` planck 2013 results . xxii .
constraints on inflation , '' astron .
astrophys .   * 571 * , a22 ( 2014 ) [ arxiv:1303.5082 [ astro-ph.co ] ] .
v.  a.  rubakov , `` the null energy condition and its violation , '' phys .
usp .   * 57 * , 128 ( 2014 ) [ arxiv:1401.4024 [ hep - th ] ] .
i.  sawicki and a.  vikman , `` hidden negative energies in strongly accelerated universes , '' phys .
d * 87 * , no .
6 , 067301 ( 2013 ) [ arxiv:1209.2961 [ astro-ph.co ] ] .
v.  a.  rubakov , `` consistent null - energy - condition violation : towards creating a universe in the laboratory , '' phys .
d * 88 * , 044015 ( 2013 ) [ arxiv:1305.2614 [ hep - th ] ] .
b.  elder , a.  joyce and j.  khoury , `` from satisfying to violating the null energy condition , '' phys .
d * 89 * , 044027 ( 2014 ) [ arxiv:1311.5889 [ hep - th ] ] .
j.  martin and r.  h.  brandenberger , `` the transplanckian problem of inflationary cosmology , '' phys .
d * 63 * , 123501 ( 2001 ) [ hep - th/0005209 ] .
a.  a.  starobinsky , `` robustness of the inflationary perturbation spectrum to transplanckian physics , '' pisma zh .  eksp .
teor .  fiz .
* 73 * , 415 ( 2001 ) [ jetp lett .   *
73 * , 371 ( 2001 ) ] [ astro - ph/0104043 ] ; t.  tanaka , `` a comment on transplanckian physics in inflationary universe , '' astro - ph/0012431 .
d.  pirtskhalava , l.  santoni , e.  trincherini and p.  uttayarat , `` inflation from minkowski space , '' jhep * 1412 * , 151 ( 2014 ) [ arxiv:1410.0882 [ hep - th ] ] .
g.  w.  horndeski , `` second - order scalar - tensor field equations in a four - dimensional space , '' int .
j.  theor .
* 10 * , 363 ( 1974 ) .
c.  deffayet , x.  gao , d.  a.  steer and g.  zahariade , `` from k - essence to generalised galileons , '' phys .
d * 84 * , 064039 ( 2011 ) [ arxiv:1103.3260 [ hep - th ] ] .
t.  kobayashi , m.  yamaguchi and j.  yokoyama , `` generalized g - inflation : inflation with the most general second - order field equations , '' prog .
phys .   * 126 * , 511 ( 2011 ) [ arxiv:1105.5723 [ hep - th ] ] .
m.  zumalacrregui and j.  garca - bellido , `` transforming gravity : from derivative couplings to matter to second - order scalar - tensor theories beyond the horndeski lagrangian , '' phys .
d * 89 * , 064046 ( 2014 ) [ arxiv:1308.4685 [ gr - qc ] ] .
x.  gao , `` unifying framework for scalar - tensor theories of gravity , '' phys .
d * 90 * , no .
8 , 081501 ( 2014 ) [ arxiv:1406.0822 [ gr - qc ] ] . c.  lin , s.  mukohyama , r.  namba and r.  saitou , `` hamiltonian structure of scalar - tensor theories beyond horndeski , '' jcap * 1410 * , no .
10 , 071 ( 2014 ) [ arxiv:1408.0670 [ hep - th ] ] .
j.  gleyzes , d.  langlois , f.  piazza and f.  vernizzi , `` exploring gravitational theories beyond horndeski , '' jcap * 1502 * , no .
02 , 018 ( 2015 ) [ arxiv:1408.1952 [ astro-ph.co ] ] .
x.  gao , `` hamiltonian analysis of spatially covariant gravity , '' phys .
d * 90 * , 104033 ( 2014 ) [ arxiv:1409.6708 [ gr - qc ] ] .
m.  fasiello and s.  renaux - petel , `` non - gaussian inflationary shapes in @xmath280 theories beyond horndeski , '' jcap * 1410 * , no .
10 , 037 ( 2014 ) [ arxiv:1407.7280 [ astro-ph.co ] ] .
l.  parker , `` particle creation in expanding universes , '' phys .  rev
.  lett .
* 21 * , 562 ( 1968 ) . l.  parker , `` quantized fields and particle creation in expanding universes . 1 .
, '' phys .
rev .   * 183 * , 1057 ( 1969 ) .
y.  b.  zeldovich and a.  a.  starobinsky , `` particle production and vacuum polarization in an anisotropic gravitational field , '' sov .
jetp * 34 * , 1159 ( 1972 ) [ zh .
* 61 * , 2161 ( 1971 ) ] .
n.  d.  birrell , p.  c.  w.  davies and l.  h.  ford , `` effects of field interactions upon particle creation in robertson - walker universes , '' j.  phys .  a * 13 * , 961 ( 1980 ) .
l.  h.  ford , `` gravitational particle creation and inflation , '' phys .
d * 35 * , 2955 ( 1987 ) .
t.  kunimitsu and j.  yokoyama , `` higgs condensation as an unwanted curvaton , '' phys .
d * 86 * , 083541 ( 2012 ) [ arxiv:1208.2316 [ hep - ph ] ] .
k.  enqvist , s.  kasuya and a.  mazumdar , `` mssm higgses as the source of reheating and all matter , '' phys .
lett .   * 93 * , 061301 ( 2004 ) [ hep - ph/0311224 ] .