over the last half - a - dozen years , it has been increasingly realized that the detection of non - gaussianities in the primordial perturbations can considerably help in constraining the inflationary models ( see refs .
@xcite ; for early efforts in this direction , see refs .
@xcite ) . in particular ,
the detection of a high value for the dimensionless non - gaussianity parameter @xmath0 that is often used to describe the amplitude of the scalar bi - spectrum can rule out a wide class of models . for instance , if the extent of non - gaussianity actually proves to be as large as the mean values of @xmath0 arrived at from the wmap data ( see refs .
@xcite ; in this context , also see refs .
@xcite ) , then canonical scalar field models that lead to slow roll inflation and nearly scale invariant primordial spectra will cease to be consistent with the data
.    however , the evaluation of the scalar bi - spectrum in a generic inflationary model proves to be a non - trivial task  @xcite .
often , it is the slow roll approximation that is resorted to in order to arrive at analytical expressions for the scalar bi - spectrum and the non - gaussianity parameter  @xmath0  @xcite .
while the slow roll approximation can actually encompass a relatively wide class of models , evidently , the approach will cease to be valid when departures from slow roll occur . in such situations , one
is often forced to approach the problem numerically . in this work
, we shall present a new and accurate fortran code , the bi - spectra and non - gaussianity operator ( bingo ) , to numerically evaluate the scalar bi - spectrum and the non - gaussianity parameter  @xmath0 for single field inflationary models involving the canonical scalar field .
although , some partial numerical results have already been published in the literature , we believe that it is for the first time that a general and efficient code has been put together to evaluate the complete scalar bi - spectrum .
though , we shall largely restrict our discussion in this paper to the equilateral case , as we shall illustrate , the code can compute the scalar bi - spectrum and the parameter  @xmath0 for any triangular configuration of the wavevectors .
also , as we shall demonstrate , bingo can also compute _ all _ the different contributions to the bi - spectrum .
moreover , it is worth noting that , under certain conditions , the code can arrive at the results within a matter of a few minutes .
we should mention here that we have made a limited version of bingo , viz .
one that focuses on the equilateral limit , available online at https://www.physics.iitm.ac.in/~sriram/bingo/bingo.html .    as an immediate application
, we shall utilize the code to examine of the power of the non - gaussianity parameter @xmath0 to discriminate between different inflationary models that admit deviations from slow roll and generate similar features in the scalar power spectrum .
recall that , most single field inflationary models naturally lead to an extended period of slow roll and hence to nearly scale invariant primordial power spectra ( see , for example , any of the following texts  @xcite or reviews  @xcite ) , which seem to be fairly consistent with the recent data on the anisotropies in the cosmic microwave background ( cmb ) ( in this context , see refs .
@xcite ; for a comparison of a class of inflationary models with the data , see , for example , ref .
@xcite and references therein ) as well as other observational constraints .
however , it has been repeatedly noticed that certain features in the inflationary scalar power spectrum can improve the fit to the cmb data at the cost of some additional parameters ( see , for instance , refs .
though the statistical significance of these features remain to be understood satisfactorily  @xcite , they gain importance from the phenomenological perspective of comparing the models with the data , because only a smaller class of single field inflationary models , which allow for departures from slow roll , can generate them .
interestingly , demanding the presence of features in the scalar power spectrum seems to generically lead to larger non - gaussianities ( see , for example , refs .
therefore , features may offer the only route ( unless one works with non - vacuum initial states  @xcite ) for the canonical scalar fields to remain viable if @xmath0 turns out to be significant .
if indeed the presence of features turns out to be the correct reason behind possibly large non - gaussianities , can we observationally identify the correct underlying inflationary scenario , in particular , given the fact that different models can lead to similar features in the scalar power spectrum ?
in other words , to what extent can the non - gaussianity parameter @xmath0 help us discriminate between the inflationary models that permit features ?
to address this question , we shall consider a few typical inflationary models leading to features , assuming that they can be viewed as representatives of such a class of scenarios . concretely
, we shall consider the starobinsky model  @xcite and the punctuated inflationary scenario  @xcite , both of which result in a sharp drop in power at large scales that is followed by oscillations .
we shall also study large and small field models with an additional step introduced in the inflaton potential  @xcite .
the step leads to a burst of oscillations in the scalar power spectrum which improve the fit to the outliers near the multipole moments of @xmath1 and @xmath2 in the cmb anisotropies .
we shall also consider oscillating inflaton potentials such as the one that arises in the axion monodromy model which lead to modulations in the power spectrum over a wide range of scales  @xcite .
the plan of this paper is as follows . in the following section
, we shall quickly describe a few essential details pertaining to the power spectrum and the bi - spectrum as well as the non - gaussianity parameter @xmath0 in inflationary models involving a single , canonical , scalar field . in sec .
[ sec : nc - sbs ] , after demonstrating that the super hubble contributions to the bi - spectrum during inflation prove to be negligible , we shall describe the method that bingo adopts to numerically compute the bi - spectrum and the non - gaussianity parameter @xmath0 .
we shall also illustrate the extent of accuracy of bingo by comparing the numerical results with the analytical results available in three situations .
firstly , restricting ourselves to the equilateral limit , we shall compare the results from the numerical computations with the spectral dependence of the bi - spectrum expected in power law inflation .
secondly , focusing on the case of the archetypical quadratic potential , we shall contrast the numerical results with the analytical expressions obtained using the slow roll approximation for an arbitrary triangular configuration of the wavevectors .
thirdly , considering a more non - trivial situation involving a brief period of fast roll , we shall compare the numerical results with the analytical results that have recently been obtained in the case of the starobinsky model in the equilateral limit  @xcite . in the succeeding two sections
, we shall utilize bingo to examine of the power of the non - gaussianity parameter @xmath0 to discriminate between different inflationary models that permit deviations from slow roll and generate similar features in the scalar power spectrum .
after a quick outline of the inflationary models of our interest , we shall discuss the scalar power spectra that arise in these models .
we shall then present the main results , and compare the @xmath0 that arise in the various models that we consider .
we shall finally conclude in sec .
[ sec : d ] with a brief summary and outlook .
a few words on our conventions and notations are in order before we proceed .
we shall work with units such that @xmath3 , and we shall set @xmath4 . as is often done in the context of inflation , we shall assume the background to be described by the spatially flat , friedmann - lemaitre - robertson - walker line - element .
we shall denote the conformal time as @xmath5 , while @xmath6 shall represent the number of e - folds .
moreover , an overprime shall denote differentiation with respect to the conformal time coordinate @xmath5 .
lastly , @xmath7 shall denote the scale factor , with the corresponding hubble parameter being defined as @xmath8 .
in this section , we shall first rapidly summarize the essential definitions of the power spectrum and the bi - spectrum as well as the non - gaussianity parameter @xmath0 .
we shall then list the various contributions to the bi - spectrum that arise in the maldacena formalism ( for the original discussions , see refs .
@xcite ; for more recent discussions , see , for example , refs .
@xcite ) .
let us begin by recalling some essential points concerning the power spectrum .
on quantization , the operator corresponding to the curvature perturbation @xmath9 can be expressed as @xmath10,\label{eq : cr - d}\end{aligned}\ ] ] where @xmath11 and @xmath12 are the usual creation and annihilation operators that satisfy the standard commutation relations .
the modes @xmath13 are governed by the differential equation  @xcite @xmath14 where @xmath15 , with @xmath16 being the first slow roll parameter .
the dimensionless scalar power spectrum @xmath17 is defined in terms of the correlation function of the fourier modes of the curvature perturbation @xmath18 as follows : @xmath19 where @xmath20 denotes the bunch - davies vacuum which is defined as @xmath21 @xmath22 @xmath23  @xcite . in terms of the modes
@xmath13 , the scalar power spectrum is given by @xmath24 as is well known , numerically , the initial conditions are imposed on the modes @xmath13 when they are well inside the hubble radius , and the power spectrum is evaluated at suitably late times when the modes are sufficiently outside the hubble radius ( see , for instance , refs .
we shall discuss the details concerning the numerical evolution of the background as well as the perturbations and the computation of the corresponding power spectrum in the next section .
the scalar bi - spectrum @xmath25 is related to the three point correlation function of the fourier modes of the curvature perturbation @xmath26 , evaluated at the end of inflation , say , at the conformal time @xmath27 , as follows  @xcite : @xmath28 in our discussion below , for the sake of convenience , we shall set @xmath29 the observationally relevant non - gaussianity parameter  @xmath0 is introduced through the relation @xmath30 , \label{eq : fnl - i}\ ] ] where @xmath31 denotes the gaussian quantity , and the factor of @xmath32 arises due to the relation between the bardeen potential and the curvature perturbation during the matter dominated epoch . upon making use of the corresponding relation between @xmath33 and @xmath31 in fourier space and the wick s theorem
, one obtains that  @xcite @xmath34.\end{aligned}\ ] ] this expression can then be utilized to arrive at the following relation between the non - gaussianity parameter @xmath0 and the scalar bi - spectrum @xmath25 or , equivalently , the quantity @xmath35 that we have introduced  @xcite : @xmath36^{-1}.\label{eq : fnl - d}\end{aligned}\ ] ] in particular , in the equilateral limit , i.e. when @xmath37 , this expression simplifies to @xmath38 two points need to be emphasized here regarding the above expressions for the quantities @xmath0 and @xmath39 .
firstly , it ought to be noted that the non - gaussianity @xmath0 parameter introduced in eq .
( [ eq : fnl - i ] ) [ as well as the expression arrived at in eq .
( [ eq : fnl - d ] ) ] corresponds to the so - called local form of the parameter .
secondly , the quantity @xmath39 corresponds to the equilateral limit of the _ local _
@xmath0 parameter , and is intrinsically different from a similar parameter introduced when the bi - spectrum has the equilateral form ( in this context , see , for instance , ref .
@xcite ) .      in the maldacena formalism
@xcite , the bi - spectrum is evaluated using the standard rules of perturbative quantum field theory , based on the interaction hamiltonian that depends cubically on the curvature perturbation .
it can be shown that the bi - spectrum that results from the interaction hamiltonian can be expressed as  @xcite @xmath40\ ;   \cg_{_{c}}(\vka,\vkb,\vkc)\nn\\   & & + \l[f_{\vka}^{\ast}(\ee)\ , f_{\vkb}^{\ast}(\ee)\,f_{\vkc}^{\ast}(\ee)\r]\ ; \cg_{_{c}}^{\ast}(\vka,\vkb,\vkc)\biggr\ } \nonumber \\ & & + g_{7}(\vka,\vkb,\vkc).\label{eq : g}\end{aligned}\ ] ] the quantities @xmath41 with @xmath42 correspond to the six terms in the interaction hamiltonian , and are described by the integrals @xmath43,\label{eq : cg3}\\
\cg_{4}(\vka,\vkb,\vkc ) & = & i\,\int_{\ei}^{\ee } \d\eta\ , a^2\,\epsilon_{1}\ , \epsilon_{2}'\ , \l(f_{\vka}^{\ast}\,f_{\vkb}^{\ast}\ , f_{\vkc}'^{\ast}+{\rm two~permutations}\r),\label{eq : cg4}\\ \cg_{5}(\vka,\vkb,\vkc ) & = & \frac{i}{2}\,\int_{\ei}^{\ee } \d\eta\ ,   a^2\ , \epsilon_{1}^{3}\ , \biggl[\l(\f{\vka\cdot\vkb}{\kb^{2}}\r)\ , f_{\vka}^{\ast}\,f_{\vkb}'^{\ast}\ , f_{\vkc}'^{\ast } \nonumber \\ & &   + { \rm five~permutations}\biggr],\label{eq : cg5}\\ \cg_{6}(\vka,\vkb,\vkc )   & = & \frac{i}{2}\,\int_{\ei}^{\ee } \d\eta\ , a^2\ ,   \epsilon_{1}^{3}\ , \biggl\{\l[\f{\ka^{2}\,\l(\vkb\cdot\vkc\r)}{\kb^{2}\,\kc^{2}}\r]\ ,   f_{\vka}^{\ast}\ , f_{\vkb}'^{\ast}\ , f_{\vkc}'^{\ast } \nonumber \\ & & + { \rm two~permutations}\biggr\},\label{eq : cg6}\end{aligned}\ ] ] where @xmath44 is the second slow roll parameter that is defined with respect to the first as follows : @xmath45 . the lower limit of the above integrals , viz .
@xmath46 , denotes a sufficiently early time when the initial conditions are imposed on the modes .
the additional , seventh term @xmath47 arises due to the field redefinition , and its contribution to @xmath35 is given by @xmath48
in this section , after illustrating that the super - hubble contributions to the complete bi - spectrum during inflation proves to be negligible , we shall outline the methods that bingo adopts to numerically evolve the equations governing the background and the perturbations , and eventually evaluate the inflationary scalar power and bi - spectra . also , we shall illustrate the extent of accuracy of bingo by comparing the results from the code with the expected form of the bi - spectrum in the equilateral limit in power law inflation and the analytical results that are available in the case of the starobinsky model  @xcite .
we shall also contrast the numerical results with the analytical results obtained under the slow roll approximation in the case of the popular quadratic potential for an arbitrary triangular configuration of the wavevectors .
it is clear from the expressions in the previous section that the evaluation of the bi - spectrum involves integrals over the mode @xmath13 and its derivative @xmath49 as well as the slow roll parameters @xmath50 , @xmath44 and the derivative @xmath51 .
while evaluating the power spectra , it is well known that , in single field inflationary models , it suffices to evolve the curvature perturbation from an initial time when the modes are sufficiently inside the hubble radius to a suitably late time when the amplitude of the curvature perturbation settles down to a constant value on super - hubble scales ( see , for example , refs .
we shall illustrate that many of the contributions to the bi - spectrum prove to be negligible when the modes evolve on super - hubble scales .
interestingly , we shall also show that , those contributions to the bi - spectrum which turn out to be significant at late times when the modes are well outside the hubble radius are canceled by certain other contributions that arise . as a consequence , we shall argue that , numerically , it suffices to evaluate the integrals over the period of time during which the curvature perturbations have been conventionally evolved to arrive at the power spectra , viz . from the sub - hubble to the super - hubble scales .
in fact , we had recently shown that , for cosmological modes that are on super - hubble scales during preheating , the contributions to the bi - spectrum and , equivalently , to the non - gaussianity parameter @xmath0 , prove to be negligible  @xcite .
the calculations we had presented in the work can be easily extended to the case of the modes on super - hubble scales during inflation .
therefore , we shall rapidly sketch the essentials arguments below , referring the readers to our recent work for further details  @xcite .      during inflation ,
when the modes are on super - hubble scales , it is well known that the solution to @xmath13 can be written as  @xcite @xmath52 where @xmath53 and @xmath54 are @xmath23-dependent constants which are determined by the initial conditions imposed on the modes in the sub - hubble limit .
the first term involving @xmath53 is the growing mode , which is actually a constant , while the term containing @xmath54 represents the decaying mode .
therefore , on super - hubble scales , the mode @xmath13 simplifies to @xmath55 moreover , the leading non - zero contribution to its derivative is determined by the decaying mode , and is given by @xmath56 where we have set @xmath57 .
it is now a matter of making use of the above solutions for @xmath13 and @xmath49 to determine the super - hubble contributions to the bi - spectrum during inflation .      to begin with ,
note that , each of the integrals @xmath58 , where @xmath59 , can be divided into two parts as follows : @xmath60 the integrals in the first term @xmath61 run from the earliest time ( i.e. @xmath46 ) when the smallest of the three wavenumbers @xmath62 , @xmath63 and @xmath64 is sufficiently inside the hubble radius [ typically corresponding to @xmath65 to the time ( say , @xmath66 ) when the largest of the three wavenumbers is well outside the hubble radius [ say , when @xmath67 .
then , evidently , the second term @xmath68 will involve integrals which run from the latter time @xmath66 to the end of inflation at  @xmath69 ( in this context , see fig .
[ fig : eta ] ) .
-10pt    in what follows , we shall discuss the various contributions to the bi - spectrum due to the terms @xmath68 .
we shall show that the corresponding contribution either remains small or , when it proves to be large , it is exactly canceled by another contribution to the bi - spectrum .
let us first focus on the fourth term @xmath70 since it has often been found to lead to the largest contribution to the bi - spectrum when deviations from slow roll occur  @xcite .
as the slow roll parameters turn large towards the end of inflation , we can expect this term to contribute significantly at late times .
however , as we shall quickly illustrate , such a late time contribution is exactly canceled by the contribution from @xmath47 which arises due to the field redefinition  @xcite . upon using the form  ( [ eq : fk - shs ] ) of the mode @xmath71 and its derivative  ( [ eq : fkp - shs ] ) on super - hubble scales in the expression  ( [ eq : cg4 ] ) ,
the integral involved can be trivially carried out with the result that the corresponding contribution to the bi - spectrum can be expressed as @xmath72\ , \biggl[\vert a_{\vka}\vert^2\ , \vert a_{\vkb}\vert^2\ , \l(a_{\vkc}\ , { \bar b}_{\vkc}^{\ast}-a_{\vkc}^{\ast}\ ,   { \bar b}_{\vkc}\r)\nn\\ & & + \,{\rm two~permutations}\biggr].\label{eq : g4-d - prh}\end{aligned}\ ] ] the wronskian corresponding to the equation of motion  ( [ eq : de - fk ] ) and the standard bunch - davies initial condition leads to the relation @xmath73 , which can then be utilized to arrive at the following simpler expression  @xcite : @xmath74\nn\\ & & \times\ , \biggl[\vert a_{\vka}\vert^2\ , \vert a_{\vkb}\vert^2 + \,~{\rm two~permutations}\biggr].\label{eq : g4-at - es}\end{aligned}\ ] ] the first of these terms involving the value of @xmath44 at the end of inflation _ exactly _ cancels the contribution @xmath47 [ with @xmath13 set to @xmath53 in eq .
( [ eq : g7 ] ) ] that arises due to the field redefinition .
but , the remaining contribution can not be ignored and needs to be taken into account .
it is useful to note that this term is essentially the same as the one due to the field redefinition , but which is now evaluated on super - hubble scales ( i.e. at @xmath66 ) rather than at the end of inflation . in other words , if we consider the fourth and the seventh terms together , it is equivalent to evaluating the contribution to the bi - spectrum corresponding to @xmath75 , and adding to it the contribution due to the seventh term @xmath47 evaluated at @xmath66 rather than at the end of inflation .
let us now turn to the contribution due to the second term , which can occasionally prove to be comparable to the contribution due to the fourth term  @xcite . upon making use of the behavior of the mode @xmath13 and its derivative on super - hubble scales in the integral  ( [ eq : cg2 ] )
, one obtains the contribution to the bi - spectrum due to @xmath76 to be  @xcite @xmath77,\end{aligned}\ ] ] where the quantity @xmath78 is described by the integral @xmath79 note that , due to the quadratic dependence on the scale factor , actually , @xmath78 is a rapidly growing quantity at late times .
however , the complete super - hubble contribution to the bi - spectrum vanishes identically since the integral @xmath78 is a purely real quantity  @xcite .
hence , in the case of the second term , it is sufficient to evaluate the contribution to the bi - spectrum due to @xmath80 .
due to their structure , one finds that the first and the third terms and the fifth and the sixth terms can be evaluated together . on super - hubble scales , one can easily show that the contributions due to the first and the third terms can be written as @xmath81\;i_{13}(\ee,\es),\end{aligned}\ ] ] while the corresponding contributions due to the fifth and the sixth terms can be obtained to be @xmath82\;i_{56}(\ee,\es ) , \label{eq : g56-shc}\end{aligned}\ ] ] where the quantities @xmath83 and @xmath84 are described by the integrals @xmath85 and @xmath86 hence , the non - zero , super - hubble contribution to the bi - spectrum is determined by the complete contribution due to the first , the third , the fifth and the sixth terms arrived at above . in order to illustrate that this contribution is insignificant , we shall now turn to estimating the amplitude of the corresponding contribution to the non - gaussianity parameter  @xmath0 .
let us restrict ourselves to the equilateral limit for simplicity .
in such a case , the super - hubble contributions due to the first , the third , the fifth and the sixth terms to the non - gaussianity parameter  @xmath0 can be obtained by substituting the quantities  ( [ eq : g13-shc ] ) and  ( [ eq : g56-shc ] ) above in the expression  ( [ eq : fnl - eq ] ) . it is straightforward to show that the corresponding @xmath0 is given by @xmath87,\end{aligned}\ ] ] where we have made use of the fact that @xmath88 at late times in order to arrive at the power spectrum .    to estimate the above super - hubble contribution to the non - gaussianity parameter @xmath89 ,
let us focus on inflation of the power law form . in power law inflation ,
the scale factor is given by @xmath90 with @xmath91 and @xmath92 being constants , while @xmath93 is a free index .
it is useful to note that , in such a case , the first slow roll parameter is a constant , and is given by @xmath94 . under these conditions , the quantity within the brackets involving @xmath53 and @xmath95 in the expression above can be easily evaluated ( in this context , see ref .
@xcite ) to arrive at @xmath96\ , \gamma^2\l(\gamma + \f{1}{2}\r)\ , 2^{2\,\gamma+1}\ , \l(2\,\gamma+1\r)\ , ( \gamma+2)\nn\\   & & \times\,\l(\gamma+1\r)^{-2\,(\gamma+1)}\ , \sin\,(2\,\pi\,\gamma)\ , \l[1-\f{h_{\rm s}}{h_{\rm e}}\,{\rm e}^{-3\,(\ne-\ns)}\r]\nn\\ & & \times\,\l(\frac{k}{a_{\rm s}\ , h_{\rm s}}\r)^{-(2\,\gamma+1)}.\end{aligned}\ ] ] it should be mentioned that , in arriving at this expression , for convenience , we have set @xmath92 to be @xmath97 , which corresponds to @xmath91 being @xmath98 , viz .
the scale factor at @xmath66 .
moreover , while @xmath99 and @xmath100 denote the e - folds corresponding to @xmath66 and @xmath69 , @xmath101 and @xmath102 represent the hubble scales at these times , respectively .
recall that , @xmath66 denotes the conformal time when the largest wavenumber of interest , say , @xmath103 , is well outside the hubble radius , i.e. when @xmath104 .
since @xmath105 is expected to be at least @xmath2 for the smallest cosmological scale , it is clear that the factor involving @xmath106 $ ] can be completely neglected .
observations point to the fact that @xmath107 . therefore , if we further assume that @xmath108 , where @xmath109 , we find that the above estimate for the non - gaussianity parameter reduces to @xmath110 where , in obtaining the final value , we have set @xmath111 .
the inequality above arises due to the fact that , for larger scales , i.e. when @xmath112 , @xmath113 at @xmath66 . in models involving the canonical scalar field , the smallest values of @xmath0 are typically generated in slow roll inflationary scenarios , wherein the non - gaussianity parameter has been calculated to be of the order of the first slow roll parameter  @xcite .
the above estimate clearly points to fact that the super - hubble contributions to the complete bi - spectrum and the non - gaussianity parameter @xmath0 can be entirely ignored .    in summary ,
to determine the scalar bi - spectrum , it suffices to evaluate the contributions to the bi - spectrum due to the quantities @xmath114 , with @xmath59 , which involve integrals running from the initial time @xmath46 to the time @xmath66 when the smallest of the three modes reaches super - hubble scales ( cf .
[ fig : eta ] ) .
further , the addition of the contribution due to the field redefinition evaluated at @xmath66 ensures that no non - trivial super - hubble contributions are ignored . in the following sub - section , with the help of a specific example
, we shall also corroborate these conclusions numerically .
the scalar bi - spectrum and the parameter @xmath0 can be easily evaluated analytically in the slow roll inflationary scenario  @xcite . however , barring some exceptional cases  @xcite , it often proves to be difficult to evaluate the bi - spectrum analytically when departures from slow roll occur .
hence , one has to resort to numerical computations in such cases .
bingo solves the background as well as the perturbation equations using rksuite  @xcite , which is a publicly available routine to solve ordinary differential equations .
we shall treat the number of e - folds as the independent variable , which allows for an efficient and accurate computation . to obtain the power spectrum ,
we impose the standard bunch - davies initial conditions on the perturbations when the modes are well inside the hubble radius , and evolve them until suitably late times .
typically , in the case of smooth inflaton potentials , it suffices to evolve the modes @xmath13 from an initial time when @xmath115 .
however , in certain models wherein the potentials contain oscillatory terms , the modes may have to be evolved from deeper inside the hubble radius .
for example , in the case of the axion monodromy model that we shall discuss in due course ( see sub - sec .
[ subsec : oip ] ) , for values of the model parameters of our interest , we find that we need to evolve the modes from an initial time when @xmath116 , so that the resonance that occurs in the model due to the oscillations in the potential is captured  @xcite .
the scalar power spectra that arise in the various models of our interest , as displayed in fig .
[ fig : sps - all ] , have all been evaluated at super - hubble scales , say , when @xmath117 , which is typically when the amplitude of the curvature perturbations freeze in .    having obtained the behavior of the background and the modes , bingo carries out the integrals involved in arriving at the bi - spectrum using the method of adaptive quadrature  @xcite .
it is useful to note that , in the equilateral limit of the bi - spectrum , which we shall largely concentrate on in this work , we can evolve each of the modes of interest independently and calculate the integrals for the modes separately .
the integrals @xmath118 actually contain a cut off in the sub - hubble limit , which is essential for singling out the perturbative vacuum  @xcite .
numerically , the presence of the cut off is fortunate since it controls the contributions due to the continuing oscillations that would otherwise occur .
we should highlight here that such a cut off procedure was originally introduced in the earliest numerical efforts to evaluate the scalar bi - spectrum in situations involving deviations from slow roll ( in this context , see refs .
@xcite ) . generalizing the cut off that is often introduced analytically in the slow roll case
, we shall impose a cut off of the form @xmath119 $ ] , where @xmath120 is a small parameter .
in the previous two sub - sections , we had discussed as to how the integrals need to be carried out from the early time @xmath46 when the largest scale is well inside the hubble radius to the late time @xmath66 when the smallest scale is sufficiently outside .
if one is focusing on the equilateral configuration , rather than integrate from @xmath46 to @xmath66 , it suffices to compute the integrals for the modes from the time when each of them satisfy the sub - hubble condition , say , @xmath121 , to the time when they are well outside the hubble radius , say , when @xmath122 . in other words ,
one carries out the integrals exactly over the period the modes are evolved to obtain the power spectrum .
the presence of the cut - off ensures that the contributions at early times , i.e. near @xmath46 , are negligible .
furthermore , it should be noted that , in such a case , the corresponding super - hubble contribution to @xmath89 will saturate the bound  ( [ eq : fnl - shs - ub ] ) in power law inflation for all the modes .    with the help of a specific example ,
let us now illustrate that , for a judicious choice of  @xmath120 , the results from bingo are largely independent of the upper and the lower limits of the integrals .
in fact , working in the equilateral limit , we shall demonstrate these points in two steps for the case of the standard quadratic potential [ see eq .
( [ eq : qp ] ) ] .
firstly , focusing on a specific mode , we shall fix the upper limit of the integral to be the time when @xmath123 . evolving the mode from different initial times
, we shall evaluate the integrals involved from these initial times to the fixed final time for different values of @xmath120 .
this exercise helps us to identify an optimal value for @xmath120 when we shall eventually carry out the integrals from @xmath124 .
secondly , upon choosing the optimal value for @xmath120 and integrating from @xmath124 , we shall calculate the integrals for different upper limits . for reasons
outlined in the previous two sub - sections , it proves to be necessary to consider the contributions to the bi - spectrum due to the fourth and the seventh terms together .
moreover , since the first and the third , and the fifth and the sixth , have similar structure in the equilateral limit , it turns out to be convenient to club these terms as we have discussed before . in fig .
[ fig : g - delta ] , we have plotted the value of @xmath125 times the different contributions to the bi - spectrum , viz .
@xmath126 , @xmath127 , @xmath128 and @xmath129 , as a function of @xmath120 when the integrals have been carried out from @xmath130 of @xmath131 , @xmath132 and @xmath133 for a mode which leaves the hubble radius around @xmath134 e - folds before the end of inflation in the case of the quadratic potential .
the figure suggests that the term @xmath128 is fairly independent of the cut - off parameter @xmath120 .
this can be attributed to the fact that the integral @xmath135 depends on the quantity @xmath51 , which effectively behaves as a cut off .
for this reason , actually , we shall _ not _ introduce the cut off while evaluating the term @xmath136 .
moreover , the figure clearly indicates @xmath137 to be a highly suitable value for the other terms . a larger @xmath120 leads to a sharper cut off reducing the value of the integrals .
one could work with a smaller @xmath120 , in which case , the figure suggests that , one would also need to necessarily integrate from deeper inside the hubble radius .
[ 0][1][1.0]@xmath138 [ 0][1][1.0]@xmath120    -10pt    in fig .
[ fig : g - nshs - all ] , after fixing @xmath120 to be @xmath139 and , with the initial conditions imposed at @xmath140 , we have plotted the four contributions to the bi - spectrum for a mode that leaves the hubble radius at @xmath141 e - folds before the end of inflation as a function of the upper limit of the integrals .
[ 0][1][1.0]@xmath138 [ 0][1][0.9]@xmath142 [ 0][1][0.9]@xmath130    -10pt    it is evident from the figure that the values of the integrals converge quickly once the mode leave the hubble radius . for efficient numerical integration , as in the case of the power spectrum , we have chosen the super - hubble limit to correspond to @xmath143 .
we have repeated similar tests for the other models of our interest too .
these tests confirm the conclusions that we have arrived at above , indicating the robustness of the numerical methods and procedures that we have adopted .
it is important that we stress three points further at this stage of our discussion .
firstly , in situations involving departures from slow roll , the dominant contribution , viz .
the @xmath128 term , does not require the introduction of the cut off at all .
secondly , in the absence of analytical expressions to compare the numerical results with , it is imperative that a suitable value for the cut off parameter @xmath120 is arrived at by repeating the exercise that we have carried out in fig .  [ fig : g - delta ] for the specific model of interest . in this context
, it is further important to appreciate the fact that the integrals involved will depend on both the cut off parameter @xmath120 as well as the initial value of @xmath144 .
a complete check seems to be mandatory in order to ensure that a unique result is arrived at .
( as we mentioned above , we have indeed carried out this exercise for the models that we consider here . )
lastly , we also need to emphasize that the above comments and conclusions would not apply to certain special cases wherein the amplitude of the curvature perturbation can evolve on super - hubble scales such as , say , in ultra slow roll inflation ( in this context , see , for instance , ref .
@xcite ) .
before we go on to consider the bi - spectra generated in the inflationary models of our interest , in order to demonstrate the accuracy of bingo , we shall compare the numerical results from the code with the analytical results that can be arrived at in three situations . in the following sub - sections
, we shall compare the results from bingo with : ( i )  the spectral dependence that can be arrived at in the context of power law inflation in the equilateral limit , ( ii )  the slow roll results as applied to the case of the conventional quadratic potential for a generic triangular configuration of the wavevectors , and ( iii )  the analytic results that can be obtained for the starobinsky model in the equilateral limit .
the first case that we shall discuss is power law inflation wherein , as we shall soon outline , the spectral shape of the non - zero contributions to the bi - spectrum can be easily arrived at in the equilateral limit .
consider the case of power law inflation described by the scale factor  ( [ eq : pli ] ) with @xmath145 .
in such a case , as we have seen , @xmath50 is a constant and , hence , @xmath44 and @xmath51 , which involve derivatives of @xmath50 , reduce to zero .
since the contributions due to the fourth and the seventh terms , viz . @xmath146 and @xmath147 ,
depend on @xmath51 and @xmath44 , respectively [ cf .
( [ eq : cg4 ] ) and  ( [ eq : g7 ] ) ] , these terms vanish identically in power law inflation . as is well known , one can express the modes @xmath13 as @xmath148 , where @xmath149 is the mukhanov - sasaki variable . in power law inflation
, one finds that , the variable @xmath149 depends only on the combination @xmath150 ( see , for example , ref .
moreover , since @xmath50 is a constant in power law expansion , we have @xmath151 . under these conditions , in the equilateral limit , with a simple rescaling of the variable of integration in the expressions  ( [ eq : cg1 ] ) , ( [ eq : cg2 ] ) , ( [ eq : cg3 ] ) , ( [ eq : cg5 ] ) and  ( [ eq : cg6 ] ) , it is straightforward to show that the quantities @xmath152 , @xmath153 , @xmath154 , @xmath155 and @xmath156 , all depend on the wavenumber as @xmath157 .
then , upon making use of the asymptotic form of the modes @xmath13 , it is easy to illustrate that the corresponding contributions to the bi - spectrum , viz .
@xmath158 , @xmath159 and @xmath129 , all behave as @xmath160 . since the power spectrum in power law inflation is known to have the form @xmath161 ( see , for example , refs .
@xcite ) , the expression  ( [ eq : fnl - eq ] ) for @xmath39 then immediately suggests that the quantity will be strictly scale invariant for all @xmath93 .
in fact , apart from these results , it is also simple to establish the following relation between the different contributions : @xmath162 , a result , which , in fact , also holds in slow roll inflation  @xcite . in other words , in power law inflation
, it is possible to arrive at the spectral dependence of the non - zero contributions to the bi - spectrum without having to explicitly calculate the integrals involved .
further , one can establish that the non - gaussianity parameter @xmath39 is exactly scale independent for any value of @xmath93 . while these arguments do not help us in determining the amplitude of the various contributions to the bi - spectrum or the non - gaussianity parameter , their spectral shape and the relative magnitude of the above - mentioned terms provide crucial analytical results to crosscheck our numerical code . in fig .
[ fig : pl ] , we have plotted the different non - zero contributions to the bi - spectrum computed using our numerical code and the spectral dependence we have arrived at above analytically for two different values of  @xmath93 in the case of power law inflation .
we have also indicated the relative magnitude of the first and the third and the fifth and the sixth terms arrived at numerically .
lastly , we have also illustrated the scale independent behavior of the non - gaussianity parameter @xmath39 for both the values of @xmath93 .
it is clear from the figure that the numerical results agree well with the results and conclusions that we arrived at above analytically .
[ 0][1][1.5]@xmath138 [ 0][1][1.5]@xmath163 [ 0][1][1.5]@xmath39    -10pt      let us now turn to the example of the conventional quadratic potential .
the bi - spectrum in such a case can be evaluated analytically in the slow roll approximation for an arbitrary triangular configuration of the wavevectors  @xcite . in the slow roll limit
, one finds that , the contributions due to the first , the second , the third and the seventh terms are of the same order , while the remaining terms turn out to be comparatively insignificant . for convenience ,
let us explicitly write down the contributions to the bi - spectrum due to the first , the second , the third and the seventh terms arrived at in the slow roll approximation .
they are given by @xmath164,\label{eq : g1-sr}\\ g_2(\vka,\vkb,\vkc ) & = & \f{h_{_{\rm i}}^4}{16\ , \mpl^4\ , \epsilon_1}\ , \l(\f{1}{k_1\ , k_2\,k_3}\r)^3\,\l(\vka\cdot\vkb+\rm two~permutations\r)\nn\\ & & \times\ , \l[-k_{_{\rm t}}+\f{1}{k_{_{\rm t}}}\ ,   \l(k_1\ , k_2+k_1\ , k_3+k_2\ , k_3\r ) + \f{k_1\,k_2\ ,
k_3}{k_{_{\rm t}}^{2}}\r],\\ g_3(\vka,\vkb,\vkc ) & = & -\f{h_{_{\rm i}}^4}{16\ , \mpl^4\ , \epsilon_1}\ , \l(\f{1}{k_1\ , k_2\,k_3}\r)^3\ , \biggl[\l(\vka\cdot\vkb\r)\ , \f{k_3 ^ 2}{k_{_{\rm t}}}\ , \l(2+\f{k_1+k_2}{k_{_{\rm t}}}\r)\nn\\ & & + \,{\rm two~permutations}\biggr],\\ g_7(\vka,\vkb,\vkc ) & = & 2\ , \pi^{4}\ , \epsilon_2\,\l(\f{1}{k_1\ , k_2\,k_3}\r)^3\nn\\ & & \times\ , \l[k_1^{3}\ , { \cal p}_{_{\rm s}}(k_2)\,{\cal p}_{_{\rm s}}(k_3 ) + \rm two~permutations\r],\label{eq : g7-sr}\end{aligned}\ ] ] respectively , with @xmath165 .
the quantities @xmath50 and @xmath44 are the first two slow roll parameters which are assumed to be largely constant , while @xmath166 denotes the hubble scale during slow roll inflation . the scalar power spectrum during slow roll can be written as  @xcite @xmath167 where @xmath168 and @xmath169 denotes a suitable , so - called , pivot scale .    -100pt    -50pt
-50pt    from the above expressions and , upon using the fact that @xmath170 has to be zero , one can arrive at the corresponding non - gaussianity parameter @xmath0 for an arbitrary triangular configuation of the wavevectors .
we find that it can be written as @xmath171^{-1}\nn\\ & & \times\,\biggl[\epsilon_1\,{\cal f}_{1}(\vka,\vkb,\vkc ) + \f{\epsilon_1}{2}\,{\cal f}_{2}(\vka,\vkb,\vkc)\nn\\ & & -\,\f{\epsilon_1}{2}\,{\cal f}_{3}(\vka,\vkb,\vkc ) + \frac{\epsilon_2}{2}\,{\cal f}_{7}(\vka,\vkb,\vkc)\biggr ] , \end{aligned}\ ] ] with @xmath172\biggr\},\\ { \cal f}_{2}(\vka,\vkb,\vkc )    & = & \l[1+\l(\frac{k_2}{k_1}\r)^2+\l(\f{k_3}{k_1}\r)^2\r]\nn\\   & & \times\,\biggl[1+\f{k_2}{k_1}+\f{k_3}{k_1}\nn\\ & & -\,\l(1+\f{k_2}{k_1}+\f{k_3}{k_1}\r)^{-1}\ , \l(\f{k_2}{k_1}+\f{k_3}{k_1}+\f{k_2}{k_1}\,\f{k_3}{k_1}\r)\nn\\   & & -\l(1+\f{k_2}{k_1}+\frac{k_3}{k_1}\r)^{-2}\ , \f{k_2}{k_1}\,\f{k_3}{k_1}\biggr],\\ { \cal f}_{3}(\vka,\vkb,\vkc )   & = & \l(1+\frac{k_2}{k_1}+\f{k_3}{k_1}\r)^{-1}\nn\\ & & \times\,\biggl\{\l[2+\l(1+\f{k_2}{k_1}+\f{k_3}{k_1}\r)^{-1}\ , \l(\f{k_2}{k_1}+\f{k_3}{k_1}\r)\r]\nn\\   & & \times\,\l[1-\l(\f{k_2}{k_1}\r)^2 -\l(\frac{k_3}{k_1}\r)^2\r]\nn\\ & & -\,\l[2+\l(1+\f{k_2}{k_1}+\f{k_3}{k_1}\r)^{-1}\ , \l(1+\f{k_2}{k_1}\r)\r]\nn\\   & & \times\,\l[1+\l(\f{k_2}{k_1}\r)^2 -\l(\frac{k_3}{k_1}\r)^2\r]\,\l(\f{k_3}{k_1}\r)^2\nn\\ & & -\,\l[2+\l(1+\f{k_2}{k_1}+\f{k_3}{k_1}\r)^{-1}\ , \l(1+\f{k_3}{k_1}\r)\r]\nn\\   & & \times\,\l[1-\l(\f{k_2}{k_1}\r)^2 + \l(\frac{k_3}{k_1}\r)^2\r]\,\l(\f{k_2}{k_1}\r)^2\biggr\},\\   { \cal f}_{7}(\vka,\vkb,\vkc )   & = & \l(\f{k_1}{k_\ast}\r)^{2\ , \l(n_{_{\rm s}}-1\r)}\ , \biggl[\l(\frac{k_2}{k_1}\r)^{n_{_{\rm s}}-1}\ , \l(\frac{k_3}{k_1}\r)^{n_{_{\rm s}}-1}\nn\\ & & + \l(\f{k_2}{k_1}\r)^3\,\l(\f{k_3}{k_1}\r)^{n_{_{\rm s}}-1 } + \l(\f{k_3}{k_1}\right)^3\ , \l(\f{k_2}{k_1}\r)^{n_{_{\rm s}}-1}\biggr].\end{aligned}\ ] ] in fig .
[ fig : qp - ac ] , we have plotted the above analytical expression for @xmath0 in the case of the quadratic potential as well as the corresponding numerical result from bingo for an arbitrary triangular configuration of the wavevectors  @xcite .
the two plots look remarkably alike . also , in fig .
[ fig : qp - ac - c ] , we have plotted the percentage difference between the analytical and the numerical results .
the figure suggests that the analytical and the numerical results match within @xmath173 .
it is worth stressing that part of the difference can be attributed to the slow roll approximation utilized in arriving at the analytical result .
the extent of the agreement reflects the fact that bingo performs rather well .
however , we should stress here that the level of accuracy of bingo can be improved by essentially carrying out the integrals from a larger initial value of @xmath144 ( i.e. from an earlier point in time ) and simultaneously working with a smaller value of the cut off parameter  @xmath120 .
the last of the examples that we shall discuss is the starobinsky model which involves a linear potential with a sharp change in its slope [ cf .
( [ eq : p - sm ] ) ] . in the starobinsky model , under certain conditions , the complete scalar bi - spectrum can be evaluated analytically in the equilateral limit  @xcite .    as we shall discuss in the next section , in the starobinsky model
, the change in the slope causes a brief period of fast roll , which leads to sharp features in the scalar power spectrum ( see fig .  [
fig : sps - all ] ) .
it was known that , for certain range of parameters , one could evaluate the scalar power spectrum analytically in the starobinsky model , which matches the actual , numerically computed spectrum exceptionally well  @xcite .
interestingly , it has been recently shown that , in the equilateral limit , the model allows the analytic evaluation of the complete scalar bi - spectrum too  @xcite . in fig .
[ fig : sm - all ] , we have plotted the numerical as well as the analytical results for the functions @xmath158 , @xmath159 , @xmath128 , and @xmath129 for the starobinsky model .
we have plotted for parameters of the model for which the analytical results are considered to be a good approximation  @xcite .
it is evident from the figure that the numerical results match the analytical ones very well .
importantly , the agreement proves to be excellent in the case of the dominant contribution  @xmath128 .
[ 0][1][1.0]@xmath138 [ 0][1][1.0]@xmath174    a couple of points concerning concerning the numerical results in the case of the starobinsky model ( both in fig .
[ fig : sps - all ] wherein we have plotted the power spectrum as well as in fig .
[ fig : sm - all ] which contains the bi - spectrum ) require some clarification .
the second ( and higher ) derivatives of the potential  ( [ eq : p - sm ] ) evidently contain a discontinuity .
these discontinuities needs to be smoothened in order for the problem to be numerically tractable .
the spectra and the bi - spectra in the starobinsky model we have illustrated have been computed with a suitable smoothing of the discontinuity , while at the same time retaining a sufficient level of sharpness so that they closely correspond to the analytical results that have been arrived at  @xcite .
in the remainder of this paper , as an immediate application , we shall utilize bingo to examine the power of the non - gaussianity parameter @xmath0 to help us discriminate between inflationary models that lead to similar features at the level of the scalar power spectrum .
we shall list the various inflationary models of our interest in this section , and describe the power spectra that arise in these models . in the next section
, we shall present the essential results and compare the @xmath0 generated in the different models .
broadly , the models that we shall consider can be categorized into three classes .
the first class shall involve potentials which admit a relatively mild and brief departure from slow roll .
the second class shall contain small but repeated deviations from slow roll , while the third and the last class shall involve a short but rather sharp departure from slow roll . let us now briefly outline the different inflationary models that we shall consider under these classes and discuss the scalar power spectra that are generated by them .      under the first class
, we shall consider models wherein a step has been introduced in potentials that otherwise admit only slow roll . given a potential , say , @xmath175 , we shall introduce the step by _ multiplying _ the potential by the following function : h _ ( ) = 1 +  ( ) , [ eq : step ] as is often done in the literature  @xcite .
clearly , the quantities @xmath176 , @xmath177 and @xmath178 denote the height , the location and the width of the step , respectively .
we shall consider the effects of the introduction of the step in the archetypical quadratic large field model , viz .
v ( ) = m^2 ^ 2,[eq : qp ] where @xmath179 represents the mass of the inflaton , and a small field model governed by the potential v ( ) = v_0 [1-()^p ] .
[ eq : sfm ] we shall specifically focus on the case wherein @xmath180 and @xmath181 , as it leads to a tilt that is consistent with the observations , and a smaller tensor - to - scalar ratio ( of @xmath182 ) than the above quadratic potential ( in planck units ) can be considered to be a drawback of the small field model . ] .
the introduction of the step has been shown to lead to a short period of deviation from slow roll inflation , which , in turn , produces a burst of oscillations in the scalar power spectrum , resulting in an improved fit to the data , in both these models  @xcite .
the best fit parameters in the case of the quadratic potential  ( [ eq : qp ] ) prove to be @xmath183 , @xmath184 , @xmath185 and @xmath186 . the field is assumed to start on the inflationary attractor at an initial value of @xmath187 , so that at least @xmath188 e - folds of inflation takes place .
the small field model  ( [ eq : sfm ] ) too is found to lead to a very similar spectrum and an almost identical extent of improvement in the fit to the data .
the best fit values of the parameters in this case are found to be @xmath189 , @xmath190 , @xmath191 and @xmath192 .
further , as in the quadratic case , the field is set to start on the attractor .
the starting value of the field is chosen to be @xmath193 , and it rolls through the step towards @xmath194 .
[ fig : p ] contains an illustration of the potentials corresponding to the quadratic case  ( [ eq : qp ] ) and the small field model  ( [ eq : sfm ] ) in the presence of the step  ( [ eq : step ] ) for the values of the parameters discussed above .
-10pt      the second class of models that we shall consider are those which lead to small but repeated deviations from slow roll as in the case of models which contain oscillatory terms in the inflaton potential .
the first such model that we shall consider is the conventional quadratic potential  ( [ eq : qp ] ) that is modulated by sinusoidal oscillations as follows  @xcite : v ( ) = m^2 ^ 2 [1 + (+)].[eq : qpsin ] the second model that we shall consider in this context is the axion monodromy model which is motivated by string theory  @xcite .
the inflaton potential in such a case is given by v ( ) = [+(+ ) ] .
[ eq : amm ] evidently , in these potentials , while the parameters  @xmath176 and @xmath195 describe the amplitude and the frequency of the superimposed oscillations , the parameter  @xmath196 shifts the oscillations within one period .
note that , whereas the amplitude of the oscillation is fixed in the axion monodromy model , the amplitude depends quadratically on the field in the chaotic model described by the potential  ( [ eq : qpsin ] ) .
the inflaton oscillates as it rolls down these potentials , and these oscillations persist all the way until the end of inflation .
the oscillating inflaton leads to small oscillations in the slow roll parameters ( in this context , see , for instance , ref .
@xcite ) , which , as we shall soon illustrate , results in continuing oscillations in the primordial scalar power spectrum .
it is found that , in the case of the potential  ( [ eq : qpsin ] ) , the power spectrum corresponding to @xmath197 , @xmath198 , @xmath199 and @xmath200 leads to a modest improvement in the fit to the cmb data  @xcite . whereas , the axion monodromy model  ( [ eq : amm ] ) leads to a considerable improvement in the fit for the following values of the parameters @xmath201 , @xmath202 , @xmath203 and @xmath204 .
while the field is assumed to start on the inflationary attractor at the initial value of about @xmath205 in the monodromy model , the field is evolved from @xmath206 in the case of the quadratic potential with sinusoidal oscillations  @xcite .
we have plotted the two potentials  ( [ eq : qpsin ] ) and  ( [ eq : amm ] ) in fig .
[ fig : p ] for the above values of the parameters .
we shall consider two models under the last class , both of which are known to lead to brief but sharp departures from slow roll .
the first of the inflationary models that we shall consider in this class is described by the following potential containing two parameters @xmath179 and @xmath207 : v ( ) = ^2 - ^n + ^2(n-1 ) .
[ eq : mssm - p ] the third quantity @xmath208 that appears in the potential is an integer which takes values greater than two .
such potentials are known to arise in certain minimal supersymmetric extensions of the standard model  @xcite .
it is worthwhile noting here that the case of @xmath209 has been considered much earlier for reasons similar to what we shall consider here , viz . towards producing certain features in the scalar power spectrum  @xcite . in the above potential ,
the coefficient of the @xmath210 term is chosen such that the potential contains a point of inflection at , say , @xmath211 ( i.e. the location where both @xmath212 and @xmath213 vanish ) , so that @xmath177 given by _
0 = ^.    if one starts at a suitable value of the field beyond the point the inflection in the above potential , it is found that one can achieve two epochs of slow roll inflation sandwiching a brief period of departure from inflation ( lasting for a little less than a e - fold ) , a scenario which has been dubbed as punctuated inflation  @xcite .
in fact , it is the point of inflection , around which the potential exhibits a plateau with an extremely small curvature , which permits such an evolution to be possible .
it is found that the following values for the potential parameters results in a power spectrum that leads to an improved fit to the cmb data in the @xmath209 case : @xmath214 and @xmath215 .
[ fig : p ] contains a plot of the potential  ( [ eq : mssm - p ] ) corresponding to these values of the parameters .
it should be added that the field is assumed to start from rest at an initial value of @xmath216 on the potential to arrive at the required behavior .
the second model that we shall consider in the last class is the starobinsky model  @xcite , which we had briefly discussed in the previous section .
as we shall soon illustrate , the model leads to a scalar power spectrum that has a certain resemblance to the spectrum generated by punctuated inflation .
the model consists of a linear potential with a sharp change in its slope at a given point , and can be described as follows : @xmath217 we have illustrated this potential in fig .
[ fig : p ] . evidently , while the value of the scalar field where the slope changes abruptly is @xmath218 , the slope of the potential above and below @xmath177 are given by @xmath219 and @xmath220 , respectively
moreover , the quantity @xmath221 denotes the value of the potential at @xmath222 .
a crucial assumption of the starobinsky model is that the value of @xmath223 is sufficiently large so that the behavior of the scale factor always remains close to that of de sitter .
the change in the slope causes a short period of deviation from slow roll as the field crosses  @xmath218 .
however , in contrast to the case of the punctuated inflationary scenario , where one encounters a brief departure from inflation , inflation continues uninterrupted in the starobinsky model .
we have not compared the starobinsky model with the data , and we shall work with two different sets of values for the parameters of the model .
one set corresponds to values of the parameters that allow for the comparison of the analytical results for the bi - spectrum that have been obtained in this case  @xcite with the corresponding numerical ones .
the other set shall be chosen to lead to a spectrum that closely mimics the power spectrum encountered in punctuated inflation . in the case of the former , as we had already mentioned in the caption of fig .
[ fig : sm - all ] , we shall choose the following values of the parameters : @xmath224 , @xmath225 , @xmath226 and @xmath227 , while , in the case of the latter , we shall work with the same values of @xmath228 and @xmath218 , but shall set @xmath229 , and @xmath230 . also , we shall work with an initial value of @xmath231 in the first instance and with @xmath232 in the second .
further , we shall start with field velocities that are determined by the slow roll conditions in both the cases .
we shall now discuss the scalar power spectra generated by the different inflationary models that we have listed above .
we shall highlight here certain aspects of the power spectra that arise in the different class of models .    in fig .
[ fig : sps - all ] , we have illustrated the scalar power spectrum that arises in the different models that we consider . as we had discussed in sub - sec .
[ subsec : dnm ] , the standard bunch - davies initial conditions are imposed on the perturbations when the modes are well inside the hubble radius .
the modes are evolved on to super hubble scales , and the spectrum is evaluated when the amplitude of the curvature perturbation has frozen when they are sufficiently outside the hubble radius .    [ 0][1][1.0]@xmath233 [ 0][1][1.0]@xmath163    let us first focus on the power spectra that arise in the starobinsky model  @xcite and the punctuated inflationary scenario  @xcite .
note that , both these models lead to a step like feature as well as a spike in the power spectrum .
the spikes arise due to the sharp departure from slow roll that occurs in these models .
while the first slow roll parameter @xmath50 remains small in the starobinsky model as the field crosses the transition , the second slow parameter @xmath44 turns very large briefly  @xcite . in the case of punctuated inflation
, @xmath50 itself grows to a large value thereby actually interrupting inflation for about a e - fold .
it is this property that results in a sharper spike in the power spectrum in the case of punctuated inflation than the starobinsky model .
the overall step in these models is easier to understand , and it simply arises due to the difference in the hubble scales associated with the slow roll epochs before and after the period of fast roll .
both these models also lead to oscillations before the spectra turn nearly scale invariant on small scales .
the spectra that arises in punctuated inflation , in addition to leading to a better fit to the outliers at very small multipoles ( because of the drop in power on these scales ) , also provides an improvement in the fit to the outlier at @xmath234  @xcite .
let us now turn to the cases of the models with a step in the potential and the potentials which contain oscillatory terms .
in contrast to the starobinsky model and the punctuated inflationary scenario , these models lead to much milder deviations from slow roll . as a result
, they also lead to smaller deviations from a nearly scale invariant power spectrum . in the context of models with a step in the potential
, a step at the correct location and of a suitable amplitude leads to a burst of oscillations in the power spectrum which in turn provides an improved fit to the outliers at @xmath234 and @xmath2  @xcite .
it is interesting to note in fig .
[ fig : sps - all ] that the spectra from punctuated inflation and the model with a step in the potential match briefly as they oscillate near scales corresponding to @xmath234 . clearly , it is this behavior that leads to a better fit to the outlier in the data at these scales .
while all the models that we have discussed on until now result in localized features , the models with oscillating inflation potentials are somewhat special as they contain non - local features , i.e. they lead to continuing oscillations that extend over a wide range of scales in the power spectrum .
needless to add , they arise due to the persisting oscillations of the scalar field as it rolls up and down the modulations in the potential .
interestingly , these features provide an improved fit to the data over a large range of all scales  @xcite .
we shall now discuss the bi - spectra and the non - gaussianity parameter @xmath0 arrived at numerically in the various models of our interest .    in fig .
[ fig : g - all.ps ] , we have plotted the various contributions , viz .
@xmath158 , @xmath159 , @xmath128 and @xmath129 ( in the equilateral limit ) for the punctuated inflationary scenario driven by the potential  ( [ eq : mssm - p ] ) , the quadratic potential  ( [ eq : qp ] ) with the step  ( [ eq : step ] ) , and the axion monodromy model  ( [ eq : amm ] ) which contains oscillations in the inflaton potential .    [ 0][1][1.0]@xmath138 [ 0][1][1.0]@xmath163 -15pt -15pt    these plots ( and also the ones in fig .
[ fig : sm - all ] ) clearly point to the fact that it is the combination @xmath128 that contributes the most to the scalar bi - spectrum in these cases  @xcite .    in fig .
[ fig : fnl - all ] , we have plotted the quantity @xmath39 due to the dominant contribution that arises in the various models that we have considered .
[ 0][1][1.0]@xmath39 [ 0][1][1.0]@xmath235 [ 0][1][1.0]@xmath163    -12pt    it is clear from this figure that , while in certain cases @xmath89 can prove to be a good discriminator , it can not help in others , and its ability to discriminate depends strongly on the differences in the background dynamics .
for instance , the evolution of the first two slow roll parameters are very similar when a step is introduced in either the quadratic potential or a small field model  @xcite .
hence , it is not surprising that the @xmath39 behaves in a similar fashion in both these models .
whereas , @xmath39 proves to be substantially different in punctuated inflation and the starobinsky model .
recall that , in the starobinsky model , the first slow roll parameter remains small throughout the evolution .
in contrast , it grows above unity for a very short period ( leading to a brief interruption of the accelerated expansion ) in the punctuated inflationary scenario .
it is this departure from inflation that leads to a sharp drop in the power spectrum and a correspondingly sharp rise in the parameter @xmath39 in punctuated inflation .
in fact , this occurs for modes that leave the hubble radius just before inflation is interrupted  @xcite .
however , note that , @xmath39 grows with @xmath163 at large wavenumbers in the starobinsky model .
this can be attributed to the fact that @xmath51 , which determines the contribution due to the fourth term , diverges due to the discontinuity in the second derivative of the potential  @xcite .
but , we should hasten to clarify that this can be shown to be an unphysical artifact arising due to the presence of a heaviside function in the description of the potential .
similarly , we find that @xmath39 is rather large in the axion monodromy model in contrast to the case wherein the conventional quadratic potential is modulated by an oscillatory term . the large value of @xmath89 that arises in the monodromy model can be attributed to the resonant behavior encountered in the model  @xcite .
in fact , we have also evaluated the @xmath39 for the case of quadratic potential modulated by sinusoidal oscillations , which too leads to continuing , periodic features in the scalar power spectrum  @xcite .
however , we find that the @xmath39 in such a case proves to be rather small ( of the order @xmath236 or so ) .
in this work , we have presented bingo , a code for the efficient computation of the bi - spectrum and the non - gaussianity parameter @xmath0 in single field inflationary models involving the canonical scalar field . while there has been previous efforts in the literature towards numerically calculating the bi - spectrum and the non - gaussianity parameter @xmath0 , our effort in arriving at bingo
can be said to be more complete for the following four reasons .
first and foremost , we have explicitly illustrated that the procedures that bingo adopts to evaluate the integrals involved are robust .
on the one hand , we have converged on a suitable form for the cut off required in the sub - hubble domain ( in particular , the value of the cut off parameter @xmath120 ) and the corresponding lower limit of the integrals , after a careful investigation of their effects on the amplitude of the bi - spectrum . on the other , we have shown that the super - hubble contributions to the complete bi - spectrum are negligible , which allows us to arrive at a convenient upper limit for the integrals . secondly , bingo can calculate all the contributions to the bi - spectrum .
thirdly , it can evaluate the bi - spectrum for an arbitrary triangular configuation of the wavevectors .
lastly , we have compared the numerical results from bingo with the analytical expressions available in certain specific cases , an exercise which illustrates that bingo can be accurate to better than @xmath173 . as we had mentioned in the introductory section ,
we have made a version of bingo , one that focuses on the equilateral limit , available online at https://www.physics.iitm.ac.in/~sriram/bingo/bingo.html . moreover ,
if needed , as we have pointed out earlier , the level of accuracy of bingo can be improved by carrying out the integrals from a larger initial value of @xmath144 ( i.e. from an earlier point in time ) and simultaneously working with a smaller value of the cut off parameter  @xmath120 .
in fact , in the code that is publicly available , these parameters have been allowed to be set by the user , as is desired .    after presenting bingo and illustrating the extent of its accuracy
, we had made use of the code to examine the power of the non - gaussianity parameter @xmath0 to lift degeneracies between various single field inflationary models involving the canonical scalar field at the level of the power spectrum . with this goal in mind ,
using bingo , we have evaluated the quantity @xmath39 in a slew of models that generate features in the scalar perturbation spectrum .
we find that the amplitude of @xmath39 proves to be rather different when the dynamics of the background turns out reasonably different , which , in retrospect , need not be surprising at all . for instance , models such as the punctuated inflationary scenario and the starobinsky model which lead to very sharp features in the power spectrum also lead to substantially large @xmath0 .
such possibilities can aid us discriminate between the models to some extent .
we had focused on evaluating the quantity @xmath0 in the equilateral limit .
it will be worthwhile to compute the corresponding values in the other limits as well as for the other forms ( such as the squeezed and the orthogonal ones ) of the non - gaussianity parameter  @xcite .
in particular , it will be interesting to examine if the so - called consistency relation between the local non - gaussianity parameter @xmath0 and the scalar spectral index in the squeezed limit is valid even in situations wherein extreme deviations from slow roll occur ( in this context , see refs .
@xcite ) . actually , since many of the models leading to features have interesting bi - spectral shapes , it is important that we arrive at the spectral shapes for the various models for an arbitrary triangular configurations of the wavenumbers .
it is also important that , apart from computing the local non - gaussianity parameter as we have done here , one evaluates the corresponding equilateral and the orthogonal parameters as well .
we are currently working on extending the scope of bingo to these cases .
we would like to conclude by highlighting one important point .
having computed the primordial bi - spectrum , the next logical step would be to compute the corresponding cmb bi - spectrum , an issue which we have not touched upon as it is beyond the scope of the current work . while tools seem to be available to evaluate the cmb bi - spectrum based on the first order brightness function , the contribution due to the brightness function at the second order remains to be understood satisfactorily ( in this context , see ref .
@xcite and the last reference in ref .
this seems to be an important aspect that is worth investigating closer .
we would like to acknowledge the use of the high performance computing facilities at the harish - chandra research institute , allahabad , india ( http://cluster.hri.res.in ) .
dkh also wishes to acknowledge support from the korea ministry of education , science and technology , gyeongsangbuk - do and pohang city for independent junior research groups at the asia pacific center for theoretical physics .
99 j.  maldacena , jhep  * 0305 * , 013 ( 2003 ) .
d.  seery and j.  e.  lidsey , jcap * 0506 * , 003 ( 2005 ) ; x.  chen , phys .
d * 72 * , 123518 ( 2005 ) ; x.  chen , m .- x .
huang , s.  kachru and g.  shiu , jcap * 0701 * , 002 ( 2007 ) ; d.  langlois , s.  renaux - petel , d.  a.  steer and t.  tanaka , phys .  rev .
* 101 * , 061301 ( 2008 ) ; phys .
d  * 78 * , 063523 ( 2008 ) .
x.  chen , adv .
astron .  *
2010 * , 638979 ( 2010 ) .
a.  gangui , f.  lucchin , s.  matarese and s.  mollerach , astrophys .
j.   * 430 * , 447 ( 1994 ) ; a.  gangui , phys .  rev .
d  * 50 * , 3684 ( 1994 ) ; a.  gangui and j.  martin , mon .  not .
* 313 * , 323 ( 2000 ) ; l.  wang and m.  kamionkowski , phys .
d , 063504 ( 2010 ) .
j.  dunkley _ et al .
_ , astrophys
.  j.  suppl .
* 180 * , 306 ( 2009 ) ; e.  komatsu _ et al .
_ , astrophys .
j.  suppl .
* 180 * , 330 ( 2009 ) .
d.  larson _
et al . _ ,
astrophys .
j.  suppl .
* 192 * , 16 ( 2011 ) ; e.  komatsu _ et al . _ ,
astrophys .
j.  suppl .
* 192 * , 18 ( 2011 ) . c.  l.  bennett _ et al . _ , arxiv:1212.5225v1 [ astro-ph.co ] ; g.  hinshaw _ et al .
_ , arxiv:1212.5226v1 [ astro-ph.co ] . e.  komatsu and d.  n.  spergel , phys .
d  * 63 * , 063002 ( 2001 ) ; e.  komatsu , d.  n.  spergel and b.  d.  wandelt , astrophys .
j.  * 634 * , 14 ( 2005 ) ; d.  babich and m.  zaldarriaga , phys .
d  * 70 * , 083005 ( 2004 ) ; m.  liguori , f.  k.  hansen , e.  komatsu , s.  matarrese and a.  riotto , phys .
d * 73 * , 043505 ( 2006 ) ; c.  hikage , e.  komatsu and t.  matsubara , astrophys .  j.  * 653 * ( 2006 ) 11 ( 2006 ) ; j.  r.  fergusson and e.  p.  s. shellard , phys .
d  * 76 * , 083523 ( 2007 ) ; a.  p.  s.  yadav , e.  komatsu and b.  d.  wandelt , astrophys .
j.  * 664 * , 680 ( 2007 ) ; p.  creminelli , l.  senatore and m.  zaldarriaga , jcap * 0703 * , 019 ( 2007 ) ; a.  p.  s.  yadav and b.  d.  wandelt , phys .  rev .  lett .  *
100 * , 181301 ( 2008 ) ; c.  hikage , t.  matsubara , p.  coles , m.  liguori , f.  k.  hansen and s.  matarrese , mon .  not .
soc .  * 389 * , 1439 ( 2008 ) ; o.  rudjord , f.  k.  hansen , x.  lan , m.  liguori , d.  marinucci and s.  matarrese , astrophys
.  j.  * 701 * , 369 ( 2009 ) ; k.  m.  smith , l.  senatore and m.  zaldarriaga , jcap * 0909 * , 006 ( 2009 ) ; j.  smidt , a.  amblard , c.  t.  byrnes , a.  cooray , a.  heavens and d.  munshi , phys .
d  * 81 * , 123007 ( 2010 ) ; j.  r.  fergusson , m.  liguori and e.  p.  s.  shellard , arxiv:1006.1642v2 [ astro-ph.co ] .
m.  liguori , e.  sefusatti , j.  r.  fergusson and e.  p.  s.  shellard , adv .
astron .  *
2010 * , 980523 ( 2010 ) ; a.  p.
s.  yadav and b.  d.  wandelt , arxiv:1006.0275v3 [ astro-ph.co ] ; e.  komatsu , class .
quantum grav . , 124010 ( 2010 ) . h.  funakoshi and s.  renaux - petel , arxiv:1211.3086v1 [ astro-ph.co ] .
e.  w.  kolb and m.  s.  turner , _ the early universe _ ( addison - wesley , redwood city , california , 1990 ) ; s.  dodelson , _ modern cosmology _ ( academic press , san diego , u.s.a . , 2003 ) ; v.  f.  mukhanov , _ physical foundations of cosmology _
( cambridge university press , cambridge , england , 2005 ) ; s.  weinberg , _ cosmology _ ( oxford university press , oxford , england , 2008 ) ; r.  durrer , _ the cosmic microwave background _ ( cambridge university press , cambridge , england , 2008 ) ; d.  h.  lyth and a.  r.  liddle , _ the primordial density perturbation _ ( cambridge university press , cambridge , england , 2009 ) ; p.  peter and j - p .
uzan , _ primordial cosmology _
( oxford university press , oxford , england , 2009 ) ; h.  mo , f.  v.  d.  bosch and s.  white , _ galaxy formation and evolution _
( cambridge university press , cambridge , england , 2010 ) .
h.  kodama and m.  sasaki , prog .
* 78 * , 1 ( 1984 ) ; v.  f.  mukhanov , h.  a.  feldman and r.  h.  brandenberger , phys
.  rep .
* 215 * , 203 ( 1992 ) ; j.  e.  lidsey , a.  liddle , e.  w.  kolb , e.  j.  copeland , t.  barreiro and m.  abney , rev .  mod .   phys .
* 69 * , 373 ( 1997 ) ; d.  h.  lyth and a.  riotto , phys .  rep .   *
314 * , 1 ( 1999 ) ; a.  riotto , arxiv : hep - ph/0210162 ; j.  martin , braz .
j.  phys .
* 34 * , 1307 ( 2004 ) ; lect .
notes phys .   *
669 * , 199 ( 2005 ) ; lect .
notes phys .
* 738 * , 193 ( 2008 ) ; b.  bassett , s.  tsujikawa and d.  wands , rev .
* 78 * , 537 ( 2006 ) ; w.  h. kinney , arxiv:0902.1529 [ astro-ph.co ] ; l.  sriramkumar , curr .
* 97 * , 868 ( 2009 ) ; d.  baumann , arxiv:0907.5424v1 [ hep - th ] .
j.  martin , c.  ringeval and r.  trotta , phys .
d  * 83 * , 063514 ( 2011 ) . h.  v.  peiris _ et al .
_ , astrophys .
j.  suppl .
* 148 * , 213 ( 2003 ) .
s.  hannestad , phys .
d  * 63 * , 043009 ( 2001 ) ; s.  l.  bridle , a.  m.  lewis , j.  weller and g.  efstathiou , mon .  not .
astron .
soc .   * 342 * , l72 ( 2003 ) ; p.  mukherjee and y.  wang , astrophys .
j.  * 599 * , 1 ( 2003 ) ; s.  hannestad , jcap  * 0404 * , 002 ( 2004 ) ; a.  shafieloo and t.  souradeep , phys .
d  * 70 * , 043523 ( 2004 ) ; d.  tocchini - valentini , y.  hoffman and j.  silk , mon .  not .
soc .  * 367 * , 1095 ( 2006 ) ; a.  shafieloo , t.  souradeep , p.  manimaran , p.  k.  panigrahi and r.  rangarajan , phys .
d  * 75 * , 123502 ( 2007 ) ; a.  shafieloo and t.  souradeep , phys .  rev .
d  * 78 * , 023511 ( 2008 ) ; r.  nagata and j.  yokoyama , phys .
d  * 79 * , 043010 ( 2009 ) ; g.  nicholson and c.  r.  contaldi , jcap  * 0907 * , 011 ( 2009 ) . b.  feng and x.  zhang , phys .  lett .
b  * 570 * , 145 ( 2003 ) ; m.  kawasaki and f.  takahashi , phys .
b  * 570 * , 151 ( 2003 ) ; r.  sinha and t.  souradeep , phys .
d  * 74 * , 043518 ( 2006 ) ; m.  j.  mortonson and w.  hu , phys .  rev .
d  * 80 * , 027301 ( 2009 ) .
r.  k.  jain , p.  chingangbam , j .- o .
gong , l.  sriramkumar and t.  souradeep , jcap * 0901 * , 009 ( 2009 ) ; r.  k.  jain , p.  chingangbam , l.  sriramkumar and t.  souradeep , phys .
d  * 82 * , 023509 ( 2010 ) .
j.  a.  adams , b.  cresswell and r.  easther , phys .
d  * 64 * , 123514 ( 2001 ) ; l.  covi , j.  hamann , a.  melchiorri , a.  slosar and i.  sorbera , phys .
d  * 74 * , 083509 ( 2006 ) ; j.  hamann , l.  covi , a.  melchiorri and a.  slosar , phys .
d  * 76 * , 023503 ( 2007 ) ; m.  j.  mortonson , c.  dvorkin , h.  v.  peiris and w.  hu , phys .
d  * 79 * , 103519 ( 2009 ) .
d.  k.  hazra , m.  aich , r.  k.  jain , l.  sriramkumar and t.  souradeep , jcap * 1010 * , 008 ( 2010 ) .
m.  benetti , m.  lattanzi , e.  calabrese and a.  melchiorri , phys .
d   * 84 * , 063509 ( 2011 ) .
m.  joy , v.  sahni , a.  a.  starobinsky , phys .
d  * 77 * , 023514 ( 2008 ) ; m.  joy , a.  shafieloo , v.  sahni , a.  a.  starobinsky , jcap * 0906 * , 028 ( 2009 ) .
hunt and s.  sarkar , phys .
d  * 70 * , 103518 ( 2004 ) ; phys .
d  * 76 * , 123504 ( 2007 ) ; m.  kawasaki , f.  takahashi and t.  takahashi , phys .
b  * 605 * , 223 ( 2005 ) ; j .- o .
gong , jcap  * 0507 * , 015 ( 2005 ) ; a.  ashoorioon and a.  krause , arxiv : hep - th/0607001v1 ; a.  ashoorioon , a.  krause and k.  turzynski , jcap * 0902 * , 014 ( 2009 ) ; m.  kawasaki and k.  miyamoto , arxiv:1010.3095v2 [ astro-ph.co ] ; a.  achucarro , j .- o .  gong , s.  hardeman , g.  a.  palma and s.  p.  patil , jcap * 1101 * , 030 ( 2011 ) ; k.  kumazaki , s.  yokoyama and n.  sugiyama , jcap * 1112 * , 008 ( 2011 )
. j.  m.  cline , p.  crotty and j.  lesgourgues , jcap  * 0309 * , 010 ( 2003 ) ; c.  r.  contaldi , m.  peloso , l.  kofman and a.  linde , jcap  * 0307 * , 002 ( 2003 ) ; l.  sriramkumar and t.  padmanabhan , phys .  rev .
d  * 71 * , 103512 ( 2005 ) ; d.  boyanovsky , h.  j.  de vega and n. g. sanchez , phys .
d  * 74 * , 123006 ( 2006 ) ; phys .
d  * 74 * , 123007 ( 2006 ) ; b.  a.  powell and w.  h.  kinney , phys .
d  * 76 * , 063512 ( 2007 ) ; c.  destri , h.  j.  de  vega and n.  g.  sanchez , phys .
d  * 78 * , 023013 ( 2008 ) .
j.  martin and c.  ringeval , phys .
d  * 69 * , 083515 ( 2004 ) ; phys .
d  * 69 * , 127303 ( 2004 ) ; jcap * 0501 * , 007 ( 2005 ) ; m.  zarei , phys .  rev
.  d  * 78 * , 123502 ( 2008 ) . c.  pahud , m.  kamionkowski and a.  r.  liddle , phys .
d  * 79 * , 083503 ( 2009 ) .
r.  flauger , l.  mcallister , e.  pajer , a.  westphal and g.  xu , jcap * 1006 * , 009 ( 2010 ) ; r.  flauger and e.  pajer , jcap * 1101 * , 017 ( 2011 ) .
t.  kobayashi and f.  takahashi , jcap * 1101 * , 026 ( 2011 ) .
m.  aich , d.  k.  hazra , l.  sriramkumar and t.  souradeep , arxiv:1106.2798v2 [ astro-ph.co ] m.  bridges , f.  feroz , m.  p.  hobson and a.  n.  lasenby , mon .  not .
soc .  * 400 * , 1075 ( 2009 ) ; h.  v.  peiris and l.  verde , phys .
d  * 81 * , 021302 ( 2010 ) ; z - k .
guo , d.  j.  schwarz and y .- z .
zhang , jcap  * 1108 * , 031 ( 2011 ) .
x.  chen , r.  easther and e.  a.  lim , jcap * 0706 * , 023 ( 2007 ) ; jcap * 0804 * , 010 ( 2008 ) . s.  hotchkiss and s.  sarkar , jcap * 1005 * , 024 ( 2010 ) ; s.  hannestad , t.  haugbolle , p.  r.  jarnhus and m.  s.  sloth , jcap * 1006 * , 001 ( 2010 ) ; r.  flauger and e.  pajer , jcap * 1101 * , 017 ( 2011 ) ; p.
adshead , w.  hu , c.  dvorkin and h.  v.  peiris , phys .
d  * 84 * , 043519 ( 2011 ) ; p.  adshead , c.  dvorkin , w.  hu and e.  a.  lim , phys .  rev .
d  * 85 * , 023531 ( 2012 ) .
a.  gangui , j.  martin and m.  sakellariadou , phys .
d  * 66 * , 083502 ( 2002 ) ; r.  holman and a.  j.  tolley , jcap * 0805 * , 001 ( 2008 ) ; w.  xue and b.  chen , phys .
d  * 79 * , 043518 ( 2009 ) ; p.  d.  meerburg ,
j.  p.  van der schaar and p.  s.  corasaniti , jcap * 0905 * , 018 ( 2009 ) ; x.  chen , jcap * 1012 * , 003 ( 2010 ) ; a.  ashoorioon and g.  shiu , jcap * 1103 * , 025 ( 2011 ) ; s.  kundu , jcap * 1202 * , 005 ( 2012 ) ; p.
d.  meerburg , r.  wijers and j.  p.  van der schaar , mon .  not .
soc .  * 421 * , 369 ( 2012 ) .
a.  a.  starobinsky , sov .
jetp  lett .
* 55 * , 489 ( 1992 ) .
j.  martin and l.  sriramkumar , jcap * 1201 * , 008 ( 2012 ) .
f.  arroja , a.  e.  romano and m.  sasaki , phys .
d  * 84 * , 123503 ( 2011 ) ; f.  arroja and m.  sasaki , jcap * 1208 * , 012 ( 2012 ) .
d.  k.  hazra , j.  martin and l.  sriramkumar , phys .
d  * 86 * , 063523 ( 2012 ) .
t.  bunch and p.  c.  w.  davies , proc .
lond .  a  * 360 * , 117 ( 1978 ) .
d.  s.  salopek , j.  r.  bond and j.  m.  bardeen , phys .
d  * 40 * , 1753 ( 1989 ) ; c.  ringeval , lect .
notes phys .
* 738 * , 243 ( 2008 ) . c.  dvorkin and w.  hu , phys .
d  * 81 * , 023518 ( 2010 ) ; w.  hu , phys .  rev .
d  * 84 * , 027303 ( 2011 ) .
r.  w.  brankin , i.  gladwell and l.  f.  shampine , _ rksuite :  a suite of runge - kutta codes for the initial value problem for odes _ , softreport 92-s1 , department of mathematics , southern methodist university , dallas , texas , u.s.a , 1992 .
see , for instance , http://www.netlib.org/ and http://www.nag.com/. j.  martin , h.  motohashi and t.  suyama , arxiv:1211.0083v1 [ astro-ph.co ] .
l.  f.  abbott and m.  b.  wise , nucl .
b  * 244 * , 541 ( 1984 ) ; d.  h.  lyth and e.  d.  stewart , phys .
b  * 274 * , 168 ( 1992 ) ; j.  martin and d.  j.  schwarz , phys .
d * 57 * , 3302 ( 1998 ) .
see , http://www.wolfram.com/. d.  k.  hazra , arxiv:1210.7170v1 [ astro-ph.co ] .
r.  allahverdi , j.  garcia - bellido , k.  enqvist and a.  mazumdar phys .
lett .  * 97 * , 191304 ( 2006 ) ; r.  allahverdi , k.  enqvist , j.  garcia - bellido , a.  jokinen and a.  mazumdar , jcap  * 0706 * , 019 ( 2007 ) ; j.  c.  b.  sanchez , k.  dimopoulos and d.  h.  lyth , jcap * 0701 * , 015 ( 2007 ) .
h.  m.  hodges , g.  r.  blumenthal , l.  a.  kofman and j.  r.  primack , nucl .
b  * 335 * , 197 ( 1990 ) .
r.  k.  jain , p.  chingangbam and l.  sriramkumar , jcap  * 0710 * , 003 ( 2007 ) .
p.  creminelli and m.  zaldarriaga , jcap 0410 , 006 ( 2004 ) .
s.  renaux - petel , jcap 1010 , 020 ( 2010 ) .
x.  chen , h.  firouzjahi , m.  h.  namjoo and m.  sasaki , arxiv:1301.5699v1 [ hep - th ] . c.  pitrou , j .-
uzan and f.  bernardeau , jcap * 1007 * , 003 ( 2010 ) .