inflation is a successful paradigm for the physics of the early universe @xcite . besides solving the classical problems of modern cosmology ( e.g. the flatness , entropy , and horizon problems )
, it provides primordial perturbations in perfect agreement with the observations @xcite .
a challenge for inflationary models is to protect the required flatness of the inflaton potential against radiative corrections .
this protection can be provided by an approximate shift symmetry as in models of natural inflation @xcite .
the symmetry can be broken by a controllably small amount , the most known example of this is the case of an axion field acquiring a potential from instantons .
the application of this to inflation , however , requires a greater than planckian axion decay constant @xmath1 @xcite , which may not be stable against gravitational corrections @xcite and which may be impossible to realize in string theory @xcite .
proposed solutions to this problem include using two @xcite or more @xcite axions , which provide an effective large scale evolution even if the decay constants of the original axions are sub - planckian , requiring nontrivial compactifications in string theory @xcite , suitably coupling the axion to a @xmath7-form @xcite , modifying the axion kinetic term @xcite , and slowing down the axion evolution through particle production @xcite as in warm inflation @xcite .
in particular , in the mechanism of @xcite the dissipation occurs through the production of a u(1 ) field coupled to the inflaton @xmath0 through the interaction @xmath8 ( where @xmath9 is the u(1 ) field strength , and @xmath10 its dual ) .
@xcite showed that this coupling can also affect the background evolution if the u(1 ) field is replaced by a su(2 ) field with a nonvanishing vacuum expectation value ( vev ) .
specifically , due to the interaction with the vev of the vector multiplet , the inflaton can be in slow roll even if its potential would otherwise ( i.e. , in absence of this interaction ) be too steep to give inflation .
such a model has been dubbed _ chromo - natural inflation _ in @xcite .
term besides the usual @xmath11 yang - mills term - in the limit in which the axion is close to the bottom of its potential @xcite .
perturbations of gauge - flation were studied in the last of @xcite . ]
the model is characterized by the action @xmath12 , \label{action}\end{aligned}\ ] ] where @xmath0 is the axion inflaton and @xmath13 .
we use the convention @xmath14 for the levi - civita tensor .
the vector field has the vev for space indices , and @xmath15 for internal su(2 ) indices .
the index @xmath16 should not be confused with the scale factor @xmath17 , that enters in the line element as @xmath18 $ ] .
we denote by dot a derivative with respect to physical time @xmath19 , and by prime a derivative with respect to conformal time @xmath20 . ]
@xmath21 which is chosen to give isotropic expansion .
we note that , for a generic theory with vector fields , an isotropic solution may be unstable against anisotropic perturbations .
for instance , this can be the case when a dilaton - like coupling @xmath22 of a scalar inflaton @xmath23 to a su(2 ) field with the vev ( [ a - vev ] ) or to an orthogonal u(1 ) triplet is arranged so to produce a scale invariant spectrum for the vectors .
isotropy is instead preserved if the vector fields are massive @xcite , as in the current model ( if the mass arises from an explicit breaking , one should also check that the theory has no ghosts @xcite ) .    ref .
@xcite performs a thorough analysis of the background evolution of the model .
the study of the chromo - natural inflation theory at the perturbative level was recently performed in @xcite , for @xmath24 where @xmath3 is the mass of the vector field fluctuations in this limit @xcite ( as we also discuss below ) . when this condition is realized , the vector field can be integrated out while , at the same time , leaving its mark on the inflationary dynamics ( this is an explicit realization of the so called _ gelaton _
@xcite mechanism ) . in this limit
, chromo - natural inflation is equivalent to a single scalar field @xmath25 theory ( where @xmath26 and @xmath27 ) .
the non - canonical kinetic term precisely encodes the effect of the gauge fields .
in @xcite it is also shown that an effective field theory equivalence of chromo - natural inflation to a non - canonical @xmath28 lagrangian holds as long as @xmath29 . beyond this limit
, a general perturbative study of the dynamics of the gauge field and the axion becomes necessary in order to test the stability of the theory and formulate its predictions .
this is precisely the scope of the present work .
we perform a full linear order study of scalar , tensor and vector perturbations and we show that the inflationary background solution of the model is stable for @xmath30 , and it is otherwise unstable .
a full phenomenological study of the model is beyond our purposes .
nonetheless , we explore the scalar and tensor modes production for a given choice of @xmath31 and @xmath32 for which the coupling with the vector fields is crucial to ensure slow - roll . specifically , we choose @xmath33 , and @xmath34 , while @xmath35 is a free parameter controlled by the value of @xmath36 .
the stability condition @xmath37 provides a lower bound for this ratio .
too small values are also excluded because they lead to a too red spectrum of the scalar modes ( we obtained this numerically ; this behavior is also seen in the analytic study of @xcite in the regime of validity of their analysis ) . on the other hand ,
the amount of gravity wave signal increases with @xmath35 , and a level which can be observed in the current or forthcoming experiments @xcite can be obtained even if the inflaton in the model spans a range of some orders of magnitude below @xmath38 , evading the so called lyth bound @xcite ( in contrast to what was expected in @xcite ) .
therefore , requiring a stable solution , with a sufficiently flat scalar power spectrum , and sufficiently small tensor modes provides constraints on the parameters of the model that go in opposite directions . for our choice @xmath33 , and @xmath34 we could not find any acceptable solution .
we expect that the situation should improve at larger @xmath1 , where the inflaton potential becomes flatter ( and the model approaches conventional slow roll inflation ) .
the paper is organized as follows : in sec .
[ sec : model ] we review the background evolution of the model for different ( stable or unstable ) ranges of the theory ; in sec .  [
sec : eqs - perturbations ] we start the perturbation analysis by introducing the most general decomposition of the metric , the inflaton , and the gauge field fluctuations , identifying the physical degrees of freedom and verifying that scalar , vector and tensor perturbations decouple at linear order ; we also quantize the system and define the cosmological correlators ; in secs .
[ sec : tensors ] , [ sec : vectors ] and [ sec : scalars ] we study the fluctuations for , respectively , tensor , vector and scalar modes ; finally in sec .
[ sec : conclusions ] we draw our conclusions .
+ we supplement our perturbation analysis in appendix  [ app - dg - scalar ] , where we show that neglecting scalar metric fluctuations does not affect the stability analysis of the model ( while it considerably simplifies the computations ) .
chromo - natural inflation is described by the action ( [ action ] ) of a pseudo - scalar field @xmath0 coupled to an @xmath39 gauge multiplet .
the vector multiplet has the vev ( [ a - vev ] ) .
the _ 00 _ component of einstein equations reads @xmath40 inflation can only occur if the potential @xmath41 is the major contribution to the total energy density , i.e. @xmath42 ; the parameters @xmath43 , @xmath44 , @xmath45 and @xmath46 are all much smaller than unity .
+ the equations of motion for the inflaton and the gauge field are : @xmath47 if we neglect @xmath48 , @xmath49 and @xmath50 , one can solve eqs .
( [ eqone]-[eqtwo ] ) for @xmath51 and @xmath52 @xcite @xmath53 in these equations we then assume @xmath54 and we obtain @xmath55 eq .
( [ eqthree ] ) can be rewritten as an equation of motion for the gauge field in terms of an effective potential @xmath56 where @xmath57 is minimized by @xmath58 in the minimum , the two terms of @xmath57 are parametrically equal to each other .
therefore , the same is true for the last two terms in ( [ eqfour ] ) .
as a consequence , the first term on the right hand side of ( [ eqfour ] ) dominates over the other two in the regime @xmath59 studied in @xcite .
this implies that @xcite @xmath60    for general values of the parameters , inserting @xmath61 , eq .
( [ eqfour ] ) becomes @xmath62 by looking at eq .
( [ eqone ] ) , one realizes that , when the gauge field settles in its minimum , it originates an actual damping term for the motion of the axion .
notice also that , when @xmath63 , the kinetic energy of the fields can be disregarded ( see e.g. eq .  ( [ eqthree ] ) ) ; therefore the slow - roll parameter ( which would normally receive contributions also from @xmath64 and @xmath65 ) will be , instead , mostly due to @xmath66 and @xmath67 @xmath68 finally we have @xmath69 from eq .
( [ chid ] ) we can compute the number of e - foldings @xmath70 where we defined @xmath71 .
+ let us introduce the parameter @xmath72 @xmath73 in terms of which the expression ( [ efolds ] ) rewrites @xmath74 an upper bound on @xmath75 can be obtained by `` pretending '' that the slow roll approximation holds at all values of @xmath0 , corresponding to @xmath76 .
the resulting upper limit is maximized for @xmath77 , where it evaluates to @xmath78 . in the @xmath79 limit
, the upper limit acquires the simple form @xmath80 .    in terms of @xmath72 ,
the slow roll solutions read @xmath81 where we have defined @xmath82    figure [ fig : bck ] allows to appreciate the accuracy of these slow roll approximations . in the figure
we show the evolution of @xmath51 for four choices of @xmath72 , comparing the exact evolution with the slow roll approximation ( [ slow ] ) . ; we then employ ( [ slow ] ) to have the initial value for @xmath51 and @xmath83 .
we set @xmath84 , and we then obtain the initial value for @xmath85 by solving the @xmath86 einstein equations exactly at the initial time .
we make sure that the initial value of @xmath0 leads to more than @xmath87 e - folds of inflation ( so that the slow - roll solution can be achieved ; we note that the system starts slightly displaced from the slow roll solution , since we set @xmath84 ; the displacement is however very small , since @xmath88 in the slow roll solution , and the background evolution quickly approaches the slow roll solution ) .
the evolutions shown in figure [ fig : bck ] are restricted to the final @xmath87 e - folds of inflation . ] in all cases , we take @xmath89 , as the main purpose of this model is to allow slow roll for @xmath31 . this is `` compensated '' by @xmath32 , and in all cases we fixed for definiteness @xmath34 .
we note that fixing the value of @xmath72 does specify a relation between @xmath90 and @xmath36 , but does not fix these two values .
as typical in inflationary models , the parameters can be specified only from fixing the normalization of the power spectrum of the scalar modes to the observed value @xmath91 @xcite .
the evolutions shown in the figure cover @xmath87 e - folds of inflation .
we note from the figure that @xmath51 performs very small oscillations around the slow roll solution .
we believe that they are due to the fact that @xmath92 is tracking a time dependent minimum ( [ qmin ] ) ( we do not show this here , but also this tracking is extremely accurate during the slow roll phase ) .
the oscillations appear somewhat large in the logarithmic scale chosen , but we see that they do not lead to any net departure from the slow roll solution .
moreover , in section [ sec : scalars ] we compare full numerical solutions of the scalar perturbations , obtained using the full numerical background solutions , with analytical solutions , for which the slow roll approximations are used , and we also find excellent agreement . as we discussed after eq .
( [ efolds2 ] ) , decreasing @xmath72 in the @xmath93 region , while keeping the other parameters ( including the initial value of @xmath0 ) fixed , decreases the amount of inflation . as a consequence
, @xmath94 can result in sufficient inflation only if the inflaton is initially close to the top of the potential .
another way to express this is to note that the inflaton rolls faster as @xmath72 decreases ( with the other parameters kept fixed ) , as it is clear both from the slow roll expression ( [ slow ] ) and from the figure .
perhaps the most surprising feature of the model is that , although all the solutions shown in figure [ fig : bck ] appear to be acceptable inflationary solutions ( in all cases the slow roll solution appears to be an attractor ; we note however that the background dynamics only probes homogeneous departures from the slow roll solution ) , and although they only differ from each other by the value of @xmath72 , the background solution with @xmath95 is unstable , while the other ones are stable .
this emerges from the study of the scalar perturbations around these solutions that we perform below .
in this section we discuss at the formal level how we compute , we quantize , and we solve the linearized theory for the perturbations of the model around the background solution discussed in the previous section .
the discussion is divided in two subsections . in the first subsection
we discuss how the perturbations can be divided into three groups , decoupled from each other at the linearized level . in the second subsection
we give the form of the quadratic action and we discuss how we compute the corresponding linearized theory for the perturbations .
there are @xmath96 perturbations in the system , one of the inflaton , @xmath97 of the su(2 ) vector field , and @xmath98 of the metric , that we decompose as : @xmath99 \nonumber\\ g_{00 } & = & - a^2 \left ( 1 - 2 \phi \right ) \nonumber\\ g_{0i } & = & a^2 \left ( b_i + \partial_i b \right ) \nonumber\\ g_{ij } & = & a^2 \left [ \left ( 1 + 2 \psi \right ) \delta_{ij } + 2 \partial_i \partial_j e + \partial_i e_j + \partial_j e_i + h_{ij } \right ] \nonumber\\ \label{deco1}\end{aligned}\ ] ] in this expression , @xmath100 is the su(2 ) index ( we also denote by @xmath16 the scale factor , as there is no ambiguity between it and an su(2 ) index ) , while @xmath101 ranges over the spatial coordinates .
we denote as `` tensor modes '' the perturbations @xmath102 and @xmath103 , which we impose to be transverse ( @xmath104 ) and traceless ( @xmath105 ) ; due to these properties , the tensor sector contains @xmath7 perturbations .
we denote as `` vector modes '' the perturbations @xmath106 , which we impose to be transverse ( @xmath107 ) ; due to this , the vector sector contains @xmath98 perturbations .
we denote as `` scalar modes '' the remaining @xmath108 perturbations .
we point out that the terms  tensor / vector / scalar " are appropriate for the perturbations of the metric and of the inflaton , as they indicate how these modes transform under a spatial rotation .
we extend this terminology also to the perturbations of the vector field , following the notation of @xcite , even if , strictly speaking , these terms are inappropriate ( given that the su(2 ) index has been used in the decomposition ) .
we nonetheless adopt it since the fact that the vector vev is diagonal ( @xmath109 ) plus the transversality properties that we have imposed guarantee that the tensor / vector / sectors that we have defined above remain decoupled from each other at the linearized level @xcite .
we fourier transform these modes @xmath110 where @xmath111 denotes any of the perturbations , and we study the theory in fourier space . in our stability study , we solve for the perturbations at the linearized level , and therefore we study a mode with a given momentum @xmath112 ( modes of different momenta are coupled to each other at the nonlinear level ) .
we can actually fix the orientation of @xmath112 along the @xmath113axis without loss of generality . starting from a general direction for @xmath112 , we rotate the coordinates so that @xmath114 , where @xmath115 is the unit vector along the @xmath113axis , and @xmath116 . after a general rotation ,
@xmath117 is no longer proportional to @xmath118 ; however we can re - obtain @xmath119 through a global su(2 ) rotation .
therefore , we can set @xmath120 without loss of generality .
this choice simplifies our algebra .
we need to remove the redundancy associated to general coordinate and su(2 ) transformations . under an infinitesimal coordinate transformation with parameter @xmath121 , @xmath122 and we remove the freedom of infinitesimal coordinate transformations by setting @xmath123 . consider instead an su(2 ) transformation with infinitesimal parameter @xmath124 ( with @xmath125 transverse ) . under this transformation ,
@xmath126 and we can fix the su(2 ) freedom by setting @xmath127 . clearly , also other modes of the metric and of the gauge field change under these transformations , and different gauge choices can be made .
our choices are motivated by the fact that ( i ) they completely fix the freedom , and ( ii ) they preserve all the @xmath128 and @xmath129 modes .
these perturbations are nondynamical , as they enter in the quadratic action of the perturbations without time derivative , and can be immediately integrated out .    with our gauge choices , and with @xmath130 , the decomposition ( [ deco1 ] )
acquires the explicit form @xmath131 and @xmath132 we verified explicitly that the scalar modes ( @xmath133 ) , the vector modes ( @xmath134 ) and the tensor modes ( @xmath135 ) are decoupled from each other at the linearized level .
namely , the quadratic action for the perturbations splits into three decoupled parts @xmath136      we have seen that the total action for the perturbations splits in a sum of three decoupled quadratic actions .
let us denote by @xmath137 the vector formed by the perturbations in one of these three systems .
we can perform a transformation @xmath138 so that the action for the array @xmath139 is of the type @xmath140 \label{formal}\ ] ] hermitianity of the action implies that @xmath141 is an anti - hermitian matrix , and @xmath142 a hermitian matrix ; the matrices obtained in the current model are actually real . in the following
we quantize the system ( [ formal ] ) , so to obtain the initial conditions for the modes , and we give an expression for the correlators between the modes .
this discussion summarizes the one of @xcite .
we first `` rotate '' @xmath143 where @xmath144 is a unitary matrix ( so that the kinetic term in ( [ formal ] ) remains canonical ) , satisfying @xmath145 as @xmath141 is real , @xmath144 can also be taken real , and , therefore , orthogonal .
we note that @xmath144 is not uniquely determined by this condition , and we fix it by setting @xmath146 at the initial time @xmath147 ; the goal of the present discussion is also to understand when this initial time can be set .
as we will see , the explicit solution for @xmath144 is not needed . in terms of the vector @xmath148 ,
the action becomes @xmath149 , \nonumber\\   { \tilde \omega}^2 & \equiv & r \left ( \omega^2 + k^t \ , k \right ) r^t.\end{aligned}\ ] ]    we then introduce the orthogonal matrix @xmath150 satisfying @xmath151 and we decompose @xmath152 , \label{deco - psi}\ ] ] where @xmath153 destroys / creates a quantum with the frequency @xmath154 .
these operators satisfy the algebra @xmath155 = \delta^{(3 ) } \left ( \vec{k } - \vec{p } \right ) \ , \delta_{ij}.\ ] ]    for the systems that we study , in the deep sub - horizon regime @xmath156 @xmath157 with @xmath158 constant at leading order in slow roll .
therefore , provided that @xmath159 , we can set the initial conditions in the adiabatic vacuum @xmath160 where @xmath161 is a constant arbitrary orthogonal matrix ( at the practical level , eq .
( [ h - in ] ) is an approximate solution of the equations of motion at early times ; the closer @xmath162 is to @xmath163 , the better this approximation is .
this determines how early @xmath147 needs to be taken ) .
the matrix @xmath161 is unphysical , as it drops from the equations of motion for the modes , and from the physical correlators ( which , as we will see , are given in terms of @xmath164 ) .
the freedom associated to @xmath161 is the generalization to @xmath75 fields of the freedom of changing by a constant phase the wave function in the single field case . in the final expression in ( [ h - in ] )
we have used the freedom associated to @xmath161 to set the wave functions to be real at the initial time .    combining ( [ psi - delta ] ) and ( [ deco - psi ] ) , we have @xmath165    using the fact that @xmath166 at @xmath147 , so that initially @xmath167 , we arrive at @xmath168 we start from these initial condition and evolve the equations of motion for the modes following from ( [ formal ] ) @xmath169    in this discussion , all the expressions are given in fourier space .
let us denote by @xmath170 the original fields in real space , @xmath171    we have the two point correlation functions @xmath172 where the power spectra are given by @xmath173 \nonumber\\ & = &   \frac{k^3}{2 \pi^2 } \ , { \rm re } \left [ \left ( { \cal m } \ , { \cal d } \ , { \cal d}^\dagger \ , { \cal m}^t \right)_{ij } \right ] . \label{power}\end{aligned}\ ] ] these correlators are the theoretical prediction , to be confronted with the statistical average of the corresponding quantities .
the power spectra ( [ power ] ) are the generalization to a system of @xmath75 fields @xmath170 of the standard power spectrum of single field inflation .
to summarize , starting from the original fields @xmath174 in momentum space , we perform ( [ y - delta ] ) to have a canonical kinetic term for @xmath175 .
we then decompose this field in terms of annihilation / creation operators of the physical quanta in the system ( the particles of frequencies @xmath154 ) .
we work in terms of the coefficients @xmath176 of this decomposition .
we can set the initial conditions ( [ delta - in ] ) for these coefficients , provided that the initial time is chosen sufficiently early such that @xmath177 .
starting from these initial conditions , the coefficients @xmath176 evolve according to ( [ delta - eom ] ) .
we will see that for some choice of parameters some of the perturbations become exponentially large on a timescale @xmath178 , signaling an instability of the linearized theory . for values of parameters leading to stable solutions ,
the coefficients enter in observable quantities through the power spectra ( [ power ] ) .
we introduce the two doublets @xmath179 for the left and right helicities , where @xmath180 the action for the tensor modes splits into two separate actions , one for the left and one for the right helicity doublet , which are formally of the type ( [ formal ] ) .
the action for the left - helicity doublet is characterized by @xmath181 and @xmath182 \nonumber\\ \omega_{12}^2 & = & a k \ , \frac{2 g q^2}{m_p } + \frac{\cal h}{m_p }   \left ( q ' + { \cal h } \ , q \right ) - \frac{g \lambda q^2 a \chi'}{f m_p }   \nonumber\\ & & \simeq a k h   { \rm o } \left ( \frac{1}{\sqrt{\lambda } } \right ) +     a^2 h^2   { \rm o } \left ( \frac{1}{\sqrt{\lambda } } \right ) \nonumber\\   \omega_{22}^2 & = &   k^2 - a k \left ( 2 g q + \frac{\lambda}{f } a \chi ' \right ) + g \frac {   \lambda } { f }   q a \chi '    \nonumber\\   & &   \simeq k^2 - a k h \ , { \cal a}+ a^2 h^2 { \cal b } \label{left - action}\end{aligned}\ ] ] where the first expression for each coefficient is exact , while the second one is obtained from the slow roll approximation ( [ slow ] ) ( we note that @xmath183 ) . the quantities @xmath184 and @xmath185 are both of @xmath186 . using the slow roll result ( [ eqthree ] ) with @xmath83 given by ( [ qmin ] ) , they evaluate to @xmath187 the action for the right - helicity doublet is related to ( [ left - action ] ) by @xmath188 , signaling the breaking of parity invariance induced by the evolution of the pseudo - scalar inflaton .
the interactions between the gravity wave and vector field tensor perturbations are slow - roll suppressed . however , the effective frequency squared @xmath189 of @xmath190 turns negative for an intermediate interval of time next to horizon crossing.dependence becomes negligible , and the mass term for the two vector perturbations @xmath191 reproduces the value @xmath3 given in @xcite . ]
this leads to a tachyonic growth of @xmath192 in this interval of time , and , correspondingly , to a growth of @xmath193 .
the same growth does not occur in the right - helicity sector due to the opposite sign of the linear term in @xmath194 .
the situation is analogous to that first studied in @xcite , where the interaction @xmath195 between a vector field and a pseudo - scalar rolling inflaton results in a tachyonic growth of the vector modes of a given helicity .
we note from ( [ ab ] ) that the length of the tachyonic region increases with increasing @xmath196 ( the last expression is the slow roll approximation ) .
this corresponds to a larger tensor mode production with growing @xmath197 .
this corresponds to decreasing @xmath72 in the numerical examples that we show below .
the vector modes @xmath198 are non - dynamical , and can be integrated out .
namely , they enter in the quadratic action of the perturbations without time derivatives , and therefore their equations of motion are algebraic equations in these variables ( recall that we are in momentum space ) .
when we solve these equations , we obtain an expression for the non - dynamical modes in terms of the two dynamical modes @xmath199 and @xmath200 .
we then insert this expression back in to the quadratic action for the vector modes , and obtain an action for @xmath201 only .
in other words , the non - dynamical modes do not introduce additional degrees of freedom , but are completely determined by the dynamical ones . after integrating out the non dynamical modes , we define @xmath202 where @xmath203 ( with the upper @xmath204 sign corresponding to @xmath205 and the lower @xmath206 sign corresponding to @xmath207 ) .
the modes @xmath208 are the canonical modes of the system , and are decoupled from each other at the linearized level : @xmath209 .
\end{aligned}\ ] ]    an explicit computation gives @xmath210 } {   m_p^2 p^2 \pm 2 g m_p^2 q p + 2 g^2 q^2 \left ( m_p^2 + q^2 \right ) } \label{om2vpm}\ ] ] where @xmath211 is the physical momentum , and the coefficients @xmath212 are functions of background quantities .
we note that @xmath213 we also note that the numerator in ( [ om2vpm ] ) is always positive .
the exact expressions for @xmath212 are readily obtained , but they are not illuminating .
we report here the leading expressions in the slow roll @xmath32 expansion : @xmath214   \nonumber\\ c_3 & = & h^3   \left [ \frac { g^2 q^2}{h^2 } \left ( 6 g \frac{q}{h } + \frac{\lambda}{f } \frac{\dot{\chi}}{h } \right ) + { \rm o } \left ( \frac { 1 } {   \lambda } \right )   \right ]   \nonumber\\ c_2 & = & h^4   \left [ \frac { g^2 q^2}{h^2 } \left ( 3 +   4 g   \frac{\lambda}{f } \frac{\dot{\chi}}{h } + 8 g^2 \frac{q^2}{h^2 }   \right ) + { \rm o } \left ( \frac { 1 } {   \lambda } \right )   \right ]   \nonumber\\ c_1 & = & h^5   \left [ \frac { g^4 q^4}{h^4 } \left ( 4 g \frac{q}{h } + 6   \frac{\lambda}{f } \frac{\dot{\chi}}{h } \right ) + { \rm o } \left ( \frac { 1 } {   \lambda } \right )   \right ]   \nonumber\\ c_0 & = & h^6   \left [   4 g^5 \frac{q^5}{h^5 } \frac{\lambda}{f } \frac{\dot{\chi}}{h }      + { \rm o } \left ( \frac { 1 } {   \lambda } \right )   \right ]   \nonumber\\\end{aligned}\ ] ]    using ( [ slow ] ) , it is immediate to see that , in each of the above expression , the dominant term in the square parenthesis that we have written explicitly is a @xmath215 coefficient .
moreover , the second term on the right hand side of ( [ om2vpm ] ) evaluates to @xmath216 .
therefore , at leading order in slow roll , @xmath217 in the deep sub - horizon and super - horizon regime , respectively . we conclude that the vector sector is stable .
we note that the super horizon limit of ( [ omega - v ] ) is precisely the slow roll expression of @xmath218 in the @xmath219 limit .
the two dynamical vector modes originate from the perturbations of the vector field , and the expression for the mass in this limit coincides with what found in @xcite .
after gauge fixing , the system of scalar perturbations comprises of @xmath220 mode from the inflaton , @xmath221 , @xmath222 from the su(2 ) field , @xmath223 , and @xmath224 from the metric , @xmath225 and @xmath226 . among , these @xmath227 modes , @xmath228 and @xmath229 are dynamical , while the other three modes are nondynamical ; namely they enter in the quadratic action without time derivatives , and they can be integrated out ( see the discussion at the start of the previous section ) . integrating out the metric perturbations makes the algebra extremely involved ; for this reason , in the study presented in the main text we make the approximation of setting @xmath230 from the start .
we then integrate out @xmath231 , and we denote the resulting quadratic action by @xmath232 .
we also performed the full computation , which we present in appendix [ app - dg - scalar ] .
we denote by @xmath233 the quadratic action obtained by including all the modes and by integrating out the @xmath222 nondynamical ones .
both @xmath232 and @xmath233 are functionals of the @xmath222 dynamical modes .
we expand all entries in these actions in slow roll . as we show in appendix [ app - dg - scalar ]
, all entries of the matrices of the two actions agree at the leading order in this expansion at all scales ( namely , for all values of @xmath234 ) , with a single exception .
the exception is the the @xmath235 coefficient of the frequency matrix @xmath142 , for which the agreement is excellent only up to @xmath236 .
the discrepancy that takes place afterwards is surely completely irrelevant for the stability study that we perform here ( as we will see , the instability , if present , manifests itself deeply inside the horizon ) .
not surprisingly , the metric perturbations do not affect the stability of the background solution .
very likely , this disagreement has also no significant consequence for the power spectra that we show below , since it manifests itself only after the power has frozen ( see the appendix for a detailed study ) .
we therefore set @xmath230 , and integrate out @xmath231 .
we then define @xmath237 in terms of the multiplet @xmath238 , the quadratic action for the scalar modes acquires , up to a total derivative , the form ( [ formal ] ) .
we denote as @xmath239 and @xmath240 the @xmath241 matrices entering in this action .
these matrices have the following entries @xmath242 and @xmath243 \nonumber\end{aligned}\ ] ] @xmath244 we stress that these expressions are obtained by disregarding the scalar metric perturbations , but that they are otherwise exact .
we can also verify that the eigenvalues of the @xmath245 part of @xmath246 are @xmath247 in the super horizon regime and for @xmath219 , in agreement with @xcite .
we solve the theory specified by this action following the steps outlined in subsection [ subsec : quantization ] .
we assume that , after inflation , only the inflaton field provides a sizable contribution to reheating ( we note that the energy in @xmath248 is much smaller than the inflaton energy during inflation ) . in this limit
, we have the curvature perturbation @xmath249 , with the power spectrum @xmath250 where @xmath251 is given in ( [ power ] ) , with @xmath174 being the three fields on the left hand side of ( [ canonical ] ) .    in figure
[ fig : pkt ] , we then present the time evolution of @xmath252 for a single mode ( a given @xmath194 ) and for the same choices of background parameters that we used in figure [ fig : bck ] for the background evolution .
for definiteness , we considered in all cases the mode that leaves the horizon @xmath87 e - folds before the end of inflation , and we denoted the corresponding comoving momentum by @xmath253 .
we choose the initial time of the evolution such that the mode is deeply inside the horizon at the start , and the last term in ( [ omega ] ) is negligible .
we observe that the choice @xmath95 leads to an instability of the linearized theory , while the other three cases are stable , and are characterized by a power that freezes outside the horizon , as in the standard inflationary models .
the instability of the @xmath95 choice manifests itself at some scale inside the horizon , namely @xmath254 .
even the stable solutions show a different evolution at this scale ( note in the figure the small change in the slope of @xmath255 for the stable cases at the same scale at which the @xmath95 solution becomes unstable ) .
we can see this analytically , by considering the approximated expressions of the scalar system inside the horizon .
we will obtain the scale @xmath256 times a numerical factor , where @xmath257 ( in the numerical examples that we have shown , @xmath258 ) .
let us discuss the approximation in more details .
first of all , using the slow roll conditions ( [ slow ] ) , we find @xmath259 in all cases that we have studied .
therefore , we can disregard @xmath260 in comparison to @xmath261 inside eqs .
( [ ks - dg0 ] ) and ( [ os - dg0 ] ) . from this , and from the slow roll approximations ( [ slow ] ) , we obtain , in the sub - horizon regime @xmath262 and @xmath263    we want to solve the evolution equations ( [ delta - eom ] ) until @xmath264 .
we can do so by rewriting them as @xmath265 by inserting only the terms written explicitly in ( [ k - s - sub ] ) and ( [ o - s - sub ] ) , and by treating these terms as constant .
all these approximations amount in considering only terms that contribute to the dynamics at @xmath266 .
we note for instance that @xmath267 in this regime , so that it is consistent to set @xmath268 in the approximated equation , while retaining @xmath269 and the dominant term of @xmath270 .
we also note that the fastest evolving coefficient in the matrices is @xmath271 that evolves on a @xmath272 timescale ; therefore time derivatives of the terms in ( [ k - s - sub ] ) and ( [ o - s - sub ] ) introduce at most terms with an additional factor @xmath4 , which are therefore suppressed in the @xmath273 and @xmath274 regimes .    performing these approximations
, the system ( [ app - sca - eq ] ) reduces to a set of linear second order equations with constant coefficients .
the equations split in three separate groups , one for the complex unknowns @xmath275 ( with @xmath101 ) , one for the the complex unknowns @xmath276 , and one for the three complex unknowns @xmath277 . in each group we need to solve three second order differential equations , and therefore we have @xmath227 possible solutions ; we note that the three groups have the identical set of equations , and they differ only in the initial conditions
. therefore the solutions are of the type @xmath278 where the matrices @xmath279 are integration constants ( so to match the initial conditions ) and @xmath280 are constant numbers .
the system is unstable if any of the @xmath280 has a real and positive part . by solving the system
, one can see that the only coefficient that can possibly be real is @xmath281^{1/2 } \label{xi}\end{aligned}\ ] ] and therefore @xmath282 where , using the slow roll condition , @xmath283    if the expression in parenthesis is negative the background solution is stable . otherwise , there is an instability at sufficiently large wavelengths .
therefore @xmath284 for any fixed @xmath72 , this condition is violated at sufficiently small @xmath285 .
this provides an upper bound to the amount of inflation in the model ( for any given choice of parameters ) .
one needs to verify that this upper bound is comparable with the required amount of inflation . for the evolutions shown in figure [ fig : bck ] , the inflaton field , at @xmath87 e - folds before the end of inflation , evaluates to @xmath286 for @xmath95 , to @xmath287 for @xmath288 , to @xmath289 for @xmath290 , and to @xmath291 for @xmath292 .
correspondingly , the rhs of the condition ( [ stability - param ] ) evaluates to @xmath293 for @xmath95 , to @xmath294 for @xmath288 , to @xmath295 for @xmath290 , and to @xmath296 for @xmath292 .
the criterion ( [ stability - param ] ) therefore indicates that the choice @xmath95 does not lead to a stable inflationary solution of @xmath87 e - folds ( recall that @xmath1 and @xmath297 are fixed to @xmath298 and to @xmath299 , respectively , in this example ) , while the other choices do .
this is in perfect agreement with figure [ fig : pkt ] .
( [ stability - k ] ) confirms that the instability , if it exists , takes place at @xmath300 times an order one factor .
for example , for the unstable @xmath95 choice , the criterion ( [ stability - k ] ) gives an instability starting at @xmath301 , corresponding to @xmath302 , in excellent agreement with the evolution seen in the figure .
finally , let us verify that the instability is extremely fast .
let us assume that we have a background solution with @xmath303 . at any given moment , modes with @xmath304 have @xmath305 , and are therefore stable .
moreover , the instability is also negligible at very large scales , given that @xmath306 for @xmath307 .
however , for @xmath308 , a quick study of eq .
( [ xi ] ) shows that @xmath309 for @xmath271 smaller than , but parametrically equal to @xmath310 .
( for @xmath311 , the maximum of @xmath312 is obtained for @xmath313 ) .
this corresponds to a very short instability time @xmath314 .
therefore , each mode experiences a strong instability while still inside the horizon , and we conclude that the background solution is unstable whenever the condition ( [ stability - param ] ) is violated .
we conclude this section by showing in figure [ fig : spe ] the power spectrum @xmath315 for the stable @xmath316 configurations , which we discuss in the next section .
we performed a complete study of the linear order quantum fluctuations for the chromo - natural inflation model .
we separated the metric , gauge and axion fluctuations into scalar , vector and tensor modes , verifying explicitly that they decouple at the linearized level .
we computed their equations of motion and worked out the quantization of the system .
the tensor sector consists of the two gravity wave polarizations plus two modes from the gauge field .
the gauge mode of one helicity becomes tachyonic for some finite interval of time next to horizon crossing , sourcing one gravity wave helicity .
the vector sector consists of two dynamical modes that originate from the vector field , and that remain perturbatively small .
the scalar sector contains three dynamical modes , one originating from the inflation and two from the gauge field .
we showed that , for some parameter choice , one of the eigenfrequencies of this system , that we denoted by @xmath312 , can become imaginary inside the hubble horizon , leading to a fast instability .
let us compare this with a standard result in inflation .
specifically , let us consider the gravity wave amplitude @xmath317 .
the canonical variable @xmath318 obeys @xmath319 with the dispersion relation @xmath320 . as a consequence @xmath321 outside the horizon , corresponding to a frozen amplitude of @xmath322
therefore , although the frequency of the canonical mode becomes tachyonic outside the horizon , its magnitude is not large enough to lead to a physical instability .
the situation is analogous for a test scalar field or for the scalar perturbations in the standard case .    for the case at hand ,
the eigenfrequency @xmath312 is given in eq .
( [ xi ] ) . in this discussion , for illustrative purposes , let us approximate the full expression obtained in eq .
( [ xi ] ) with @xmath323 , where @xmath324 is defined in ( [ lambd ] ) , and @xmath325 an order one factor .
no instability appears if @xmath325 is negative .
otherwise , an instability takes place for modes of physical momentum @xmath326 .
the timescale of the instability is @xmath327 , ( or , equivalently , the amplitude of the scalar modes grows proportionally to a large power of the scale factor ) which indeed corresponds to a fast instability .
this heuristic discussion reproduces the results obtained from the precise form of @xmath312 ( see section [ sec : scalars ] for the precise computation ) .
moreover , we also performed a fully numerical and exact study of the scalar perturbations that confirms these analytical results .
the parameters of the model can be chosen so that @xmath312 remains real inside the horizon ( equivalently , @xmath328 in the heuristic expression ) so that the instability is avoided .
this corresponds to choosing @xmath329 , where @xmath83 is the vev of the vector field , and @xmath36 the su(2 ) coupling .
quite interestingly , the quantity @xmath3 coincides with the mass of the fluctuations of the vector field in the @xmath219 regime @xcite , and it approximates well this mass for @xmath330 .
therefore , the inflationary solution is stable if and only if the vector field is sufficiently heavy .
we have obtained the boundary of the stable region @xmath219 analyzed in @xcite and our formalism can be readily employed to study the phenomenology of the boundary region .
this study is beyond the purposes of the present work . for illustrative purposes ,
we have however computed the power spectrum of @xmath331 for some sample choices of parameters .
specifically , figure [ fig : spe ] shows the power spectrum for an axion decay constant @xmath332 , for @xmath34 ( we stress that the main motivation for the model is to provide inflation for a sub - planckian axion decay constant , and that this can be obtained for sufficiently large @xmath297 ) , and for three choices of @xmath72 .
we note that , among those shown in figure [ fig : spe ] , only the power spectrum obtained for @xmath292 is sufficiently flat to meet the observational bounds @xcite .
we numerically found that the spectral tilt is a decreasing function of @xmath72 in the @xmath333 interval ( equivalently , it is an increasing function of @xmath35 in this interval ) .
this behavior can also be seen in the analytic relation given in @xcite in the large @xmath3 regime .
therefore , both the requirement of stability and of flatness of the power spectrum pose a lower bound on @xmath35 .    on the other hand , the discussion around eq .
( [ mgh ] ) leads to the conclusion that a large @xmath35 leads to a detectable or ruled out gravity wave signal .
the enhanced signal is parity violating , which results in nonvanishing tb and eb correlations .
a measure of the net handedness of the tensor modes is @xmath334 where @xmath335 is the power spectrum of the left / right - helicity gravity wave modes .
the corresponding observational bounds have been studied in @xcite . in figure
[ fig : r - dchi ] we compare the bounds presented in @xcite with the value of @xmath336 and @xmath337 obtained in this model , for our choice of @xmath332 and @xmath34 , and for different values of @xmath72 in the @xmath338 range ( the tensor power spectra are obtained from the @xmath235 element of ( [ power ] ) for the left and right helicity tensor sectors , following the procedure outlined in subsection [ subsec : quantization ] ) .
greater values of @xmath72 do not lead to a visible gravity wave signal even in a cosmic variance limited experiment , while lower values are ruled out by the current limit @xmath339 @xcite .
we point out that all values of @xmath72 considered in this plot are actually ruled out because they lead to a too small value of @xmath340 ( the largest value @xmath341 is obtained for @xmath342 , and @xmath340 then further decreases at larger @xmath72 ) .
we stress that our choice of @xmath332 and @xmath34 is only dictated by the requirement of @xmath343 , so that the coupling to the gauge field plays a relevant role for the inflation dynamics ( which , in turns , is the initial purpose of the model @xcite ) .
obviously , these values can be changed , and our phenomenology discussion has the only purpose of understanding what kind of limits can be imposed on the model .
we expect that a viable region will exist at larger @xmath1 ( and smaller @xmath297 ) , as the model becomes closer to a free inflaton model slowly rolling on a flat potential .    at the theoretical level
, it is interesting to note that the model can result in an observable gravity wave signal in the near future , even if the inflaton spans a range which is some order of magnitudes smaller than the planck scale .
this evades the lyth bound @xcite , that states that @xmath344 is possible only for an excursion of the inflaton of @xmath345 during the last @xmath346 e - folds of inflation .
we note that the model we have studied evades the bound because of the @xmath8 interaction .
quite interestingly , the only other examples that we are aware of where @xmath347 can be achieved even if the inflation evolution is orders of magnitude below the planck scale   have been constructed where the inflationary range is smaller than , but still of @xmath345 @xcite . ]
are those studied in refs .
@xcite , which are characterized by the same pseudo - scalar interaction . as we already mentioned , ref .
@xcite is also characterized by a @xmath348 evolution due to the damping from gauge field production . in this case
, an interesting parity - violating gravity wave signal can be generated @xcite from the produced gauge quanta , but one needs to evade the simultaneous generation of too many non - gaussian scalar density perturbations @xcite .
@xcite achieves this by considering a large number of gauge fields ( this reduces non - gaussianity by the central limit theorem ) or by introducing a curvaton field .
@xcite shows that @xmath349 can be obtained , and the non - gaussianity limit can be respected , if the rolling - scalar is not the inflaton .
an important difference is that , however , in the mechanism of @xcite the tensor modes are produced at the non - linear level by the vector fields produced by the rolling inflaton . for the model @xcite , the production occurs already at the linear level , due to its mixing with the vector mode induced by the vector vev and the non - abelian structure ( @xmath350 ) .
quite likely , also in this model the vector modes @xmath192 can source a significant amount of scalar density perturbations at the non - linear level . this may reduce @xmath336 from the level studied here , although a too large @xmath3 will still likely be ruled out by the significant gravity wave production .
it is possible that , for regimes resulting in acceptable @xmath336 , the sourced scalar modes will lead to interesting levels of non - gaussianity and primordial black holes as those obtained in @xcite .
it is also possible that , for some choice of parameters , the vector field production will `` self - regulate '' ( with a consequent decrease of @xmath336 ) due to its backreaction on the inflationary dynamics , as visible in the background evolutions studied in @xcite ( a too large production may slow down the inflaton , and this may in turns decrease the vector production ) .
all these interesting possibilities remain to be studied .    in the first version of this manuscript
, we pointed out that ( i ) the inflationary solution of @xcite is stable if and only if @xmath351 , and that ( ii ) the tensor - to - scalar ratio @xmath336 in this model is enhanced with respect to the case of a free inflaton , leading to violation of the lyth s bound @xcite for some choice of parameters .
both these claims are confirmed by the present analysis .
the stability study is unchanged with respect to the first version .
the study of the tensor modes , and the phenomenology considerations that follow from it , have instead been updated , to include the helicity violating terms in the tensor action that were erroneously missing in the first version .
such terms result in a further increase of @xmath336 and in a helicity violating gravity wave signal .
the relevance of these terms was pointed out in @xcite , that appeared on the archive between the first and the current version of this manuscript .
@xcite agrees with our limit @xmath351 for the stability of the inflationary solution , and our revised analysis agrees with the effect of the helicity violating effects found in @xcite .
mp acknowledges nicola bartolo and sabino matarrese for useful conversations . this work was supported in part by doe grant de - fg02 - 94er-40823 at the university of minnesota .
mp would like to thank the university of padova , and infn , sezione di padova , for their friendly hospitality and for partial support during his sabbatical leave .
in section [ sec : scalars ] we studied the system of linear perturbations for the model , disregarding the perturbations of the metric . here
we summarize the results for the full system , and we confirm the accuracy of the approximation made in the main text .
we start from the full set of @xmath227 scalar perturbations ( after gauge fixing ) , including the @xmath224 modes from the metric perturbations ( @xmath225 and @xmath226 ) that we have ( artificially ) set to zero in the main text . both these modes are non - dynamical , and we integrate them out , together with the other non - dynamical perturbation @xmath231 .
we are left with three dynamical modes , and we `` rotate '' them as in eq .
( [ canonical ] ) of the main text . in this way
, we obtain the full quadratic action for the perturbations . it is formally of the type @xmath352 \label{formal - s2sf}\end{aligned}\ ] ] namely , the modes ( [ canonical ] ) are not the exact canonical scalar variables . however , as we now show , they provide a very good approximation to the canonical modes in the slow roll regime .
more in general , the matrices of the full action ( [ formal - s2sf ] ) are extremely involved .
we studied them in slow roll approximation . specifically , using the slow roll approximation ( [ slow ] ) , we can cast all the elements of these matrices in the form @xmath353 where the coefficients @xmath212 and @xmath354 only depend on the parameters of the model , and on slowly evolving background quantities .
we computed the leading order expression for these coefficients in the slow roll approximation .
for example , we obtain @xmath355 } \nonumber\\ & = & 1 + \frac{{\rm o } \left ( \frac{f^2}{m_p^2 y^{3/2 } \lambda^3 } \right ) h^4}{2 h^2 p^2 + { \rm o } \left ( \frac{1}{y^{3/2 } } \right ) \left [ 1 + { \rm o } \left ( \frac{1}{\lambda } \right ) \right ] h^4 } \nonumber\\\end{aligned}\ ] ] this expression is extremely close to one , since , parametrically , the second term is @xmath356 outside the horizon , and even more suppressed inside the horizon .
in fact , the kinetic matrix differs from the identity only up to slow roll suppressed quantities : @xmath357 therefore , up to very small slow roll corrections , the modes ( [ canonical ] ) are also the canonical variable of the full scalar system .    performing the same procedure on @xmath358
, we obtain @xmath359}{p^2 + m_g^2 \left [ 1 + { \rm o } \left ( \frac{1}{\lambda } \right ) \right ] } \nonumber\\ & = & \left ( \frac{k_{\rm s}}{a } \right)_{12 }   \left [   1 + { \rm o } \left ( \frac{\lambda^{-1 } } { 1 + y^{3/2 } \ , p^2 / h^2 }   \right ) \right ] \nonumber\\ \left ( \frac{k_{\rm s , f}}{a } \right)_{13 } & = & - \frac{g^2 \lambda q^3}{\sqrt{2 } f \sqrt{p^2+m_g^2 } } \ , \frac{\left ( p^2 + m_g^2 \right ) \left [ 1 + { \rm o } \left (   \frac{1}{\lambda } \right ) \right ] } { p^2 + m_g^2 \left [ 1 + { \rm o } \left ( \frac{1}{\lambda } \right ) \right ] }   \nonumber\\ & = &   \left ( \frac{k_{\rm s}}{a } \right)_{13 } \left [ 1 + { \rm o } \left ( \lambda^{-1 } \right ) \right ] \nonumber\\ \left ( \frac{k_{\rm s , f}}{a } \right)_{23 }   & = & \frac{h^4}{\left ( p^2+m_g^2 \right)^{3/2 } } \left [ { \rm o } \left ( \frac{1}{y^{9/4 } \lambda } \right ) \frac{p^2}{h^2 } + { \rm o } \left ( \frac{1}{y^{15/4 } \lambda } \right )   \right ] \nonumber\\ & \sim & { \rm o } \left ( \frac{h}{\sqrt{p^2+h^2 } } \frac{f}{m_p } \frac{1}{\lambda^{3/2 } } \right ) \times   \left ( \frac{k_{\rm s , f}}{a } \right)_{12 } \end{aligned}\ ] ] where we recall that @xmath360 .
the @xmath97 and @xmath361 entries are in excellent agreement with those given in the main text .
we recall that the @xmath361 entry is much smaller than the @xmath97 entry inside the horizon , and it is negligible in the stability study .
we note that , for the full system , the @xmath96 element is nonvanishing , while @xmath362 . however , this element is strongly slow roll suppressed with respect to the other two , and completely negligible .    proceeding in the same way ( for brevity ,
we omit here the powers of @xmath72 ) , we obtain @xmath363 and we see that , for all these entries , the expressions of @xmath364 and of @xmath365 agree at all scales ( namely , for any value of @xmath234 ) up to subdominant @xmath366 corrections . for the @xmath235 entry , we obtain @xmath367 also on this entry , the two matrices are in perfect agreement during the full sub - horizon regime . however , while for the other entries the agreement continues also in the super - horizon regime , the @xmath235 entries differ from each other for @xmath368 . by evaluating the coefficients in ( [ om11-compa ] )
, we found that the disagreement starts only at @xmath369 for all the evolutions studied in the main text .    to conclude
, all the matrix elements , up to one exception , of the system of scalar perturbations studied in the main text ( in which we made the approximation of disregarding the scalar metric perturbations ) agree at all scales ( namely , for any value of @xmath234 ) with the corresponding entries of the full system up to subdominant terms in a slow - roll expansion .
the single exception is the @xmath235 entry of @xmath142 , for which the agreement persists during the entire sub - horizon regime , at horizon crossing , and also up to @xmath236 , but not further .
this proves that the stability study performed in the main text is valid also when the metric perturbations are included , given that the instability , when present , manifests itself deeply inside the horizon .
most likely , this guarantees that also the power spectra shown in the main text are accurate , since the disagreement manifests itself only after the powers have frozen ( see figure [ fig : pkt ] ) .
r.  kallosh , lect .
notes phys .
* 738 * , 119 ( 2008 ) [ hep - th/0702059 [ hep - th ] ] ; r.  kallosh , n.  sivanandam and m.  soroush , phys .
d * 77 * , 043501 ( 2008 ) [ arxiv:0710.3429 [ hep - th ] ] .
e.  silverstein and a.  westphal , phys
.  rev .
d * 78 * , 106003 ( 2008 ) [ arxiv:0803.3085 [ hep - th ] ] ; l.  mcallister , e.  silverstein and a.  westphal , phys .
d * 82 * , 046003 ( 2010 ) [ arxiv:0808.0706 [ hep - th ] ] ; n.  kaloper and l.  sorbo , phys .  rev
.  lett .
* 102 * , 121301 ( 2009 ) [ arxiv:0811.1989 [ hep - th ] ] . c.  germani and a.  kehagias , phys .
* 106 * , 161302 ( 2011 ) [ arxiv:1012.0853 [ hep - ph ] ] ; j.  ohashi and s.  tsujikawa , jcap * 1210 * , 035 ( 2012 ) [ arxiv:1207.4879 [ gr - qc ] ] ; d.  maity , arxiv:1209.6554 [ hep - th ] .
m.  m.  anber and l.  sorbo , phys .  rev .
d * 81 * , 043534 ( 2010 ) [ arxiv:0908.4089 [ hep - th ] ] .
s.  mohanty and a.  nautiyal , phys .
d * 78 * , 123515 ( 2008 ) [ arxiv:0807.0317 [ hep - ph ] ] ; l.  visinelli , jcap * 1109 * , 013 ( 2011 ) [ arxiv:1107.3523 [ astro-ph.co ] ] .
e.  dimastrogiovanni , n.  bartolo , s.  matarrese and a.  riotto , adv .
astron .
* 2010 * , 752670 ( 2010 ) [ arxiv:1001.4049 [ astro-ph.co ] ] ; k.  dimopoulos , int .
j.  mod .
d * 21 * , 1250023 ( 2012 ) [ erratum - ibid .
d * 21 * , 1292003 ( 2012 ) ] [ arxiv:1107.2779 [ hep - ph ] ] ; j.  soda , class .
* 29 * , 083001 ( 2012 ) [ arxiv:1201.6434 [ hep - th ] ] ; a.  maleknejad , m.  m.  sheikh - jabbari and j.  soda , arxiv:1212.2921 [ hep - th ] .
n.  bartolo , s.  matarrese , m.  peloso and a.  ricciardone , arxiv:1210.3257 [ astro-ph.co ] . m.  -a .
watanabe , s.  kanno and j.  soda , phys .
lett .   * 102 * , 191302 ( 2009 ) [ arxiv:0902.2833 [ hep - th ] ] ; s.  kanno , j.  soda and m.  -a .
watanabe , jcap * 1012 * , 024 ( 2010 ) [ arxiv:1010.5307 [ hep - th ] ] ; k.  murata and j.  soda , jcap * 1106 * , 037 ( 2011 ) [ arxiv:1103.6164 [ hep - th ] ] ; k.  yamamoto , phys .
d * 85 * , 123504 ( 2012 ) [ arxiv:1203.1071 [ astro-ph.co ] ] ; k.  -i .
maeda and k.  yamamoto , arxiv:1210.4054 [ astro-ph.co ] ; h.  funakoshi and k.  yamamoto , arxiv:1212.2615 [ astro-ph.co ] ; k.  dimopoulos , m.  karciauskas and j.  m.  wagstaff , phys .  rev .
d * 81 * , 023522 ( 2010 ) [ arxiv:0907.1838 [ hep - ph ] ] .
k.  dimopoulos , m.  karciauskas and j.  m.  wagstaff , phys .
b * 683 * , 298 ( 2010 ) [ arxiv:0909.0475 [ hep - ph ] ] ; j.  m.  wagstaff and k.  dimopoulos , phys .
d * 83 * , 023523 ( 2011 ) [ arxiv:1011.2517 [ hep - ph ] ] ; r.  namba , phys .
d * 86 * , 083518 ( 2012 ) [ arxiv:1207.5547 [ astro-ph.co ] ] .
j.  a.  r.  cembranos , c.  hallabrin , a.  l.  maroto and s.  j.  n.  jareno , phys .
d * 86 * , 021301 ( 2012 ) [ arxiv:1203.6221 [ astro-ph.co ] ] ; j.  a.  r.  cembranos , a.  l.  maroto and s.  j.  n.  jareno , arxiv:1212.3201 [ astro-ph.co ] .
b.  himmetoglu , c.  r.  contaldi and m.  peloso , phys .
lett .   * 102 * , 111301 ( 2009 ) [ arxiv:0809.2779 [ astro - ph ] ] .
b.  himmetoglu , c.  r.  contaldi and m.  peloso , phys .
d * 79 * , 063517 ( 2009 ) [ arxiv:0812.1231 [ astro - ph ] ]
. b.  himmetoglu , c.  r.  contaldi and m.  peloso , phys .
d * 80 * , 123530 ( 2009 ) [ arxiv:0909.3524 [ astro-ph.co ] ] .
d.  baumann _ et al .
_ [ cmbpol study team collaboration ] , aip conf .
* 1141 * , 10 ( 2009 ) [ arxiv:0811.3919 [ astro - ph ] ] .
d.  h.  lyth , phys .
* 78 * , 1861 ( 1997 ) [ hep - ph/9606387 ] .
a.  e.  gumrukcuoglu , b.  himmetoglu and m.  peloso , phys .
d * 81 * , 063528 ( 2010 ) [ arxiv:1001.4088 [ astro-ph.co ] ] .
s.  saito , k.  ichiki and a.  taruya , jcap * 0709 * , 002 ( 2007 ) [ arxiv:0705.3701 [ astro - ph ] ] .
v.  gluscevic and m.  kamionkowski , phys .
d * 81 * , 123529 ( 2010 ) [ arxiv:1002.1308 [ astro-ph.co ] ] .
r.  keisler , c.  l.  reichardt , k.  a.  aird , b.  a.  benson , l.  e.  bleem , j.  e.  carlstrom , c.  l.  chang and h.  m.  cho _ et al .
_ , astrophys .  j.   * 743 * , 28 ( 2011 ) [ arxiv:1105.3182 [ astro-ph.co ] ]
. i.  ben - dayan and r.  brustein , jcap * 1009 * , 007 ( 2010 ) [ arxiv:0907.2384 [ astro-ph.co ] ] .
s.  hotchkiss , a.  mazumdar and s.  nadathur , jcap * 1202 * , 008 ( 2012 ) [ arxiv:1110.5389 [ astro-ph.co ] ] .
p.  d.  meerburg and e.  pajer , arxiv:1203.6076 [ astro-ph.co ] .
a.  linde , s.  mooij and e.  pajer , arxiv:1212.1693 [ hep - th ] .
n.  barnaby , e.  pajer and m.  peloso , phys .
d * 85 * , 023525 ( 2012 ) [ arxiv:1110.3327 [ astro-ph.co ] ] .
p.  adshead , e.  martinec and m.  wyman , arxiv:1301.2598 [ hep - th ] .