Revisiting non-Gaussianity from non-attractor inflation models

Non-attractor inflation is known as the only single field inflationary scenario that can violate non-Gaussianity consistency relation with the Bunch-Davies vacuum state and generate large local non-Gaussianity. However, it is also known that the non-attractor inflation by itself is incomplete and should be followed by a phase of slow-roll attractor. Moreover, there is a transition process between these two phases. In the past literature, this transition was approximated as instant and the evolution of non-Gaussianity in this phase was not fully studied. In this paper, we follow the detailed evolution of the non-Gaussianity through the transition phase into the slow-roll attractor phase, considering different types of transition. We find that the transition process has important effect on the size of the local non-Gaussianity. We first compute the net contribution of the non-Gaussianities at the end of inflation in canonical non-attractor models. If the curvature perturbations keep evolving during the transition—such as in the case of smooth transition or some sharp transition scenarios—the 𝒪(1) local non-Gaussianity generated in the non-attractor phase can be completely erased by the subsequent evolution, although the consistency relation remains violated. In extremal cases of sharp transition where the super-horizon modes freeze immediately right after the end of the non-attractor phase, the original non-attractor result can be recovered. We also study models with non-canonical kinetic terms, and find that the transition can typically contribute a suppression factor in the squeezed bispectrum, but the final local non-Gaussianity can still be made parametrically large.


Introduction
Inflationary cosmology is the leading paradigm of the very early universe [1][2][3][4][5][6][7], in which the universe has experienced a primordial phase of quasi-de Sitter expansion.The simplest inflation model is realized by a canonical scalar field slowly rolling along a sufficiently flat potential.The associated perturbation theory successfully predicted a nearly scale-invariant power spectrum of primordial curvature perturbation, which is favoured by the latest cosmic microwave background (CMB) observations [8,9].Moreover, it is widely acknowledged that the primordial non-Gaussianity, which encodes information about the very early universe, could be a powerful tool to discriminate different inflation models or alternative scenarios [10][11][12][13].Remarkably, there is a consistency relation for non-Gaussianity in single-field slow-roll inflation models pointed out by Maldacena [14,15].The consistency relation states that the amplitude of the primordial non-Gaussianity in squeezed configuration -where the wavelength of one mode is much larger than the other two in the three point correlation fucntion -is proportional to the spectral index of the power spectrum of scalar perturbations, i.e. f NL = 5(1 − n s )/12.Accordingly, the observation of the almost scale invariant power spectrum of linear perturbation indicates extremely small amount of nonlinear correlations in squeezed limit.As a result, one expects that the simplest inflation model in terms of single slow-roll scalar field would be ruled out if any squeezed limit non-Gaussianity could be detected.
It is, however, interesting to notice that there exists a nontrivial inflationary scenario, dubbed as non-attractor inflation [16][17][18][19][20][21], that can violate Maldecena's consistency relation even in the framework of single scalar field with Bunch-Davies initial states.This is due to the fact that curvature perturbations generated from quantum fluctuations during the non-attractor phase are dominated by the growing modes at super-Hubble scales, of which the behaviour is much similar to the matter bounce cosmology [22][23][24] rather than the cosmology of slow-roll inflation.Accordingly, similar to the matter bounce cosmology [25,26], large amount of local non-Gaussianity -which contributes dominantly to the squeezed limit bispectrum -can be achieved in non-attractor inflation models.Ref. [17] considers a simple model with canonical kinetic term which predicts f NL 5/2.The idea is then further generalized to the models with non-canonical kinetic terms in [19][20][21] where it has been shown that the non-Gaussianity can be arbitrarily large.Inspired by this unconventional behaviour of primordial perturbations, many studies have been devoted to understand the possible violation of the consistency relation during the non-attractor phase from a variety of theoretical perspectives [27][28][29][30].
Furthermore, it is important to notice that the non-attractor inflation alone is not phenomenologically viable [17,18,31].Namely, without a conventional attractor phase, the non-attractor inflation does not provide enough e-folds or cannot fit the COBE normalization of the density perturbations.For a more realistic consideration, the phase of non-attractor inflation shall be regarded as some initial stage of the whole inflationary era, and a phase transition from non-attractor to the slow-roll attractor evolution becomes essential for this class of models.Therefore, the non-attractor inflation model consists of at least three different kinds of phases: the non-attractor phase, the transition phase, and the slow-roll phase.We shall define these phases more explicitly in models we study.During the transition phase, modes that exited the horizon may not freeze, the main focus of this paper is to understand how the transition process would influence primordial non-Gaussianities generated in the non-attractor phase.
In this work, we revisit primordial non-Gaussianities from non-attractor inflation by focusing on the impact of the non-attractor to attractor transition.We begin with a detailed analysis of the non-attractor inflation model with a canonical scalar field, which was previously studied in Ref. [17,31].Here the transition processes are classified into two different cases, depending on whether the background evolution around the transition is smooth or sharp.We first apply the in-in formalism to study the bispectrum in these two cases separately.For the smooth transition, our calculation shows that the non-Gaussianity generated in the nonattractor phase cannot survive through the transition to the slow-roll attractor phase.So the value f NL = 5/2 generated during the non-attractor phase returns to ∼ 0 (slow-roll-suppressed) in the slow-roll phase, and the net contribution to the local f NL is negligible as in the slow-roll attractor case.The situation is more complicated in the sharp transition.After a detailed analysis on the background and perturbations, we find that, in general the non-Gaussianity generated in the non-attractor phase is also suppressed after the transition.But in extremal cases where the curvature perturbations freeze out immediately at the transition time, the original result f NL 5/2 can be recovered.We confirm all these results by employing the simple and intuitive calculation of the δN formalism.Note that despite the non-trivial evolution of non-Gaussianity during the transition phase, the consistency relation is still violated even though the amplitude of non-Gaussianity might be slow-roll suppressed.This is a consequence of the fact that the curvature perturbation modes keep evolving after they crossed the Hubble horizon; in contrast with the conventional, slow-roll models where curvature perturbation is conserved on super-horizon scales.
We further study the transition process in the non-attractor inflation model driven by a non-canonical scalar field, as constructed in [19,20].The background evolution shows that, the inflaton field first becomes canonical before the cosmological system enters into the phase of slow-roll attractor through a smooth transition phase.The difference between these models and the above canonical model is that now we have two types of terms in the non-canonical models.The first type behaves very similarly to the interaction term in the canonical model, and it does not contribute to large local non-Gaussianity either when a smooth transition is taken into account.However, the non-canonical models have another set of qualitatively different terms.These second type of terms are unique due to the presence of the non-canonical kinetic terms.The contribution to large local non-Gaussianity from these terms do not get exactly erased by the smooth transition period, but instead gets an additional suppression factor.Since the suppression factor and the amplitude of primordial non-Gaussianity generated in the non-attractor phase are independent of each other, the large local non-Gaussianity is still possible for certain model parameters.So the main conclusions of [19,20] remain unchanged.
The paper is organized as follows.In Section 2 we study the canonical model of nonattractor inflation.After reviewing previous works, we focus on the detailed transition process from the initial non-attractor phase to the subsequent phase of slow-roll attractor.Then we elaborate on the behaviour of local non-Gaussianity in two different cases -smooth transition and sharp transition, via both in-in formalism and δN formalism.In Section 3 we generalize the study of the non-attractor inflation to models with non-canonical kinetic terms, where we only consider smooth transition case.The detailed transition process in these models is shown by full analysis of the background dynamics.After that, we estimate the size of the non-Gaussianity and find a suppression effect caused by the background evolution of the transition process.We summarize our conclusions with a discussion in Section 4. Throughout the paper we take the convention of the reduced Planck mass to be M 2 pl = 1/8πG = 1.

The canonical model
In this section we revisit the calculation of primordial non-Gaussianities in the model of canonical non-attractor inflation, and show how the different transition processes may change the non-Gaussianity generated in the non-attractor phase.

The non-attractor phase and local non-Gaussianity
The canonical non-attractor model is constructed by assuming that the inflaton's potential is almost a constant, i.e. for sufficiently large regime one has V (φ) V 0 [16,17].Accordingly, the background equations in this model are given by where a dot denotes the derivative with respect to cosmic time t, and H ≡ ȧ/a is the Hubble parameter.This leads to the following behaviour for the slow-roll parameters As shown in the above equation, the slow-roll parameter decays very quickly during the non-attractor phase, and thus, one can take the limit → 0 as a good approximation here.As a result, the Hubble parameter H is nearly constant during the non-attractor phase and in terms of conformal time τ the scale factor takes a −1/(Hτ ).In addition, the second slow-roll parameter η is of order O(1).
For the primordial curvature perturbation R, we define z ≡ a √ 2 and u k ≡ zR k .Then at the linear level, the perturbation variable u k is governed by the Mukhanov-Sasaki equation where the prime denotes the derivative to conformal time τ .Following the standard treatment, the effective mass can be written as z /z (ν 2 − 1/4)/τ 2 , where for 1, ν is given by In the non-attractor stage, η = −6, and thus, ν = 3/2.Consequently, Eq. (2.3) yields the mode function of curvature perturbation as follows, which looks the same as the one in the lowest order slow-roll approximation.But notice that is rapidly evolving here in contrary to the slow-roll case.After Hubble-exit, one can get a scale-invariant power spectrum of primordial curvature perturbation, of which the form takes P R (k) ≡ H 2 8π 2 .However, since ∝ a −6 , the amplitude of curvature perturbation grows as R k ∝ a 3 at super-Hubble scales.As a result, the final form of the power spectrum ought to be evaluated after the end of the non-attractor phase.
In order to calculate the non-Gaussianity, one needs to study the three-point correlation function of primordial curvature perturbation (2.6) At the squeezed limit k 1 k 2 k 3 , the bispectrum B R can be expressed as where f NL is the amplitude of non-Gaussianity in squeezed limit.The consistency relation, predicts f NL 5 12 (1 − n s ) which we will see is violated in non-attractor models.Notice that the local shape has the same scaling behaviour in squeezed limit, although it is well defined in any configuration [12].The non-Gaussianity that is generated during the non-attractor phase is indeed in the local shape but we are only interested in the squeezed limit (which tells us whether the consistency relation is violated or not); therefore we will not discuss non-Gaussianities in general configurations.
Ref. [17] uses two methods to compute the size of local non-Gaussianity.The first method focuses on the non-attractor phase alone.Because the contributions from the terms in cubic Lagrangian are slow-roll suppressed in this phase, Ref. [17] focuses on the contribution from a field-redefinition term in which yields at the end of the non-attractor phase τ e .If these perturbations got frozen immediately at the end of this phase and were carried along to the attractor slow-roll phase, we would end up with this order-one non-Gaussianity.However, the transition from the non-attractor phase to the slow-roll phase may not be an instant process and the process is not generically an attractor solution either.It turns out that the evolution of modes at the super-horizon scales can be non-negligible during this transition period.
The second method used in Ref. [17] indeed considers this transition, but treating it as an instant process.In this method, the field redefinition term no longer contributes because the parameter η should now be evaluated at the end of inflation instead of at the end of the non-attractor phase.This value of η is negligible.The corresponding contribution should now, equivalently, come from an interaction term in the cubic Lagrangian, (2.10) as correctly considered in Ref. [17].The bispectrum coming from this interaction term is ) where τ 0 is the conformal time after which the super-horizon curvature perturbation as well as the corresponding bispectrum cease evolving.We also remind that τ e denotes the end of the non-attractor phase.The η parameter goes from −6 to nearly zero and then the coefficient η can be comparably large.If the transition is approximated as an instant process that takes place suddenly at the time τ e when the non-attractor phase ends [17], then one may expect τ 0 = τ e and the behaviour of η during the transition period can be approximated by a step function (2.12) As a result, the interaction term (2.10) leads to and the value f NL = 5/2 will be recovered.However, one may still wonder whether this conclusion holds true if we consider a complete transition process.In the next subsections, we will study various transition cases in details, and show that, for a smooth transition the actual contribution from (2.10) is negligible; while the O(1) local non-Gaussianity can be recovered from a sharp transition.

The non-attractor to slow-roll transitions
The reason that the ultra-slow-roll inflation with a constant potential cannot be a complete model (even if we impose an abrupt cutoff and start the reheating instantly in the non-attractor phase) is that, after 40 ∼ 60 efolds, the density perturbation cannot produce the observed value. 1 So a transition to a slow-roll phase is needed.In the following, we construct a model that describes such a transition.The advantage of our model is that the exact analytical solutions can be obtained, in which the inflaton field begins the evolution in the non-attractor phase and then joins the slow-roll phase gradually.Suppose that the non-attractor phase ends at φ e = φ(τ e ), and after that, a slow-roll potential V (φ) is attached to the constant one.Since the transition process is very short and the inflaton field excursion is very tiny, during this period, the attached slow-roll potential can be expanded as follows,

.14)
Figure 1.A sketch plot of the potentials of non-attractor inflation with smooth and sharp transitions.Note that the inflaton rolls from right to left, i.e. φ is decreasing during the evolution.
Here we have introduced the potential slow-roll parameters V ≡ 1/2(V (φ e )/V (φ e )) 2 and η V ≡ V (φ e )/V (φ e ), which are expected to be small constants such that the slow-roll dynamics can be triggered after the transition.Accordingly, one can sketch the possible form of potentials depending on different values of parameters which correspond to different types of transition, as shown in Figure 1.We may distinguish two extreme possibilities: if we require the derivative of the potential to be continuous, then V = 0 and thus the transition is smooth; whereas for other cases, such as √ 2 V |η V |, we get sharp transition.Note that by considering the above potential we restricted ourselves to the case with continuous potential and the positivity of the second term also implies that the inflaton rolls-down instead of jumping up.By the end of this section, however, we will discuss how the results may change by considering non-standard cases of discontinuous potential or negative slope.Finally, notice that the above additional potential in a single field model of inflation, breaks the internal shift symmetry explicitly; therefore even the generalized consistency relations [29,30] are not applicable, unless if the bispectrum does not evolve when the potential (2.14) switches on.
In this type of inflation model, initially the inflaton field rolls along the constant potential V = V 0 for φ > φ e , which we define as the non-attractor phase.After the inflaton field reaches φ e , the potential becomes (2.14), on which inflation transits to the slow-roll attractor.We define this period as the transition phase, as shown by the light green region in Figure 1.Using e-folding number N as variable (with the convention dN = Hdt), the background equations become where we have assumed that the Hubble parameter is a constant.Without losing generality, we can set N = 0 at φ e and the field velocity at the same moment is introduced to be π e , then we have the following analytical solution ) with the parameters being introduced.Notice that in our convention π e < 0 (because φ is decreasing throughout the evolution) and hence h < 0. After some simple algebra, the slow-roll parameters defined in (2.2) during the transition are given by ) We can see from the above that, as N increases, η goes from −6 − h to −2η V during the transition.

Non-attractor initial condition
Slow-roll attractor Note that, the background evolutions in smooth and sharp transitions behave manifestly different, and h is a crucial parameter to characterize their difference.For the smooth transition, h → 0, and thus at the beginning of the transition phase η = −6, which continuously follows the non-attractor phase and then smoothly evolves to the slow-roll attractor.Figure 2 shows this behaviour via the phase space diagram and the evolution of η, where the smooth transition is depicted by the numerical solution of the non-attractor initial condition on a plateau-like potential 2 .
For the sharp transition, h is a negative constant determined by the field velocity π e at the end of the non-attractor phase.From (2.19) we see that, when the attractor is reached after the sharp transition, we have 0 V with 0 = π 2 0 /2, where π 0 is the field velocity dφ dN during the slow-roll phase.Therefore, the parameter h can be described also by the ratio between π 0 and π

.21)
From the relative magnitudes of π e and π 0 , it is straightforward to see that, the value of |h| can be of order unity or even bigger, and there are three possible cases in sharp transition: h < −6, h = −6 and −6 < h < 0, as shown in the phase space diagram in Figure 3. Consequently at the beginning of the transition η = −6 − h can be quite large, which differs from its value during the non-attractor phase (where it is η = −6).Thus there is a sudden change of η at the transition time, from −6 to −6 − h, as shown by the numerical examples in the right panel of Figure 3.For the later convenience, we formulate the evolution of η around τ e as Therefore, typically a sharp transition process consists of an instant transition at the beginning and a following period of relaxation described by (2.16) -(2.20).One special case is h = −6 + 2η V −6, where inflaton joins the slow-roll attractor immediately after the instant transition and there is no relaxation process.However, it still differs from the oversimplified case in (2.12).As we shall show in Section 2.4, this realistic instant transition does not imply immediate freezing of the curvature perturbation (i.e.τ 0 = τ e ), and the evolving super-horizon mode after the instant transition can still modify the non-Gaussianity generated during the non-attractor phase.With these background solutions of transitions, in the following we shall perform a detailed study of non-Gaussianities.The in-in formalism is applied in Section 2.3 and 2.4, for smooth and sharp transitions respectively.In Section 2.5, we further confirm the in-in results in both cases via δN formalism.

Non-Gaussianity in a smooth transition
In this subsestion, we focus on in-in calculation of the smooth transition case, which corresponds to the limit V → 0 in the potential (2.14), and demonstrate that there is a cancellation for the local non-Gaussianity generated during the non-attractor stage.Then in Section 2.3.1, we confirm this conclusion by the numerical study of a realistic model.At last, in Section 2.3.2,we perform an extended analysis to show that this conclusion holds true for smooth transition in general.
Before the in-in calculation, we should first check the behaviour of the mode function during the transition, which is governed by the Mukhanov-Sasaki equation (2.3) and the index ν in (2.4).Even though η and η varies dramatically during the transition, surprisingly the exact solution (2. 19) and (2.20) gives us ν 2 = 9/4 − 3η V , which is constant and the same as the result in slow-roll attractors 3 .Therefore, the mode function in (2.5) still applies here as the leading order approximation, and the resulting power spectrum in this period is still nearly scale-invariant.We should further remark that, the curvature perturbation still evolves during the transition, and should be fixed after the slow-roll attractor is reached.That is to say the final amplitude of the power spectrum is P R (k) ≡ H 2 8π 2 0 , where 0 is the in the slow-roll stage.
With this analytical description of the smooth transition, now let us look at the bispectrum caused by the cubic interaction term (2.10).We can substitute the mode function (2.5) into the in-in integral in (2.11).Notice that, even though is small, it varies fast during the transition, thus where the second term is due to the 's evolution4 .Taking the squeezed limit (2.24) Since η is negligible during the non-attractor and slow-roll phase, we only need to compute this integral during the transition process (from τ e to τ 0 ).As mentioned earlier, the evolution of the bispectrum after τ 0 is suppressed, as is well-known in the attractor case where the super-horizon curvature perturbation freezes out.Since we are mainly interested in the perturbation modes which exit the Hubble radius during the non-attractor phase, we can use Thus the leading order contribution of the above bispectrum becomes Plugging in the analytical expressions for and η in (2.19) and (2.20), we find after the transition Here √ 2 0 can be expressed as the field velocity dφ dN at the beginning of the slow-roll phase τ 0 , thus √ 2 0 |π e |.As a result, the local non-Gaussianity becomes negligible after the transition.
If we compare this calculation with the result (2.13) in the instant transition approximation, we can just identify τ = τ 0 and = 0 in (2.25) using the step function (2.12) for η.However, here when we compute the smooth transition explicitly, the third term in the bracket becomes negligible, since τ 0 /τ < 1 during the transition.And we have seen that there is a cancellation between the first two terms, in contrast with the instant transition approximation which gives order one result.This cancellation is also demonstrated numerically as follows.

Numerical study on a plateau-like potential
In a realistic case, non-attractor inflation can be seen as imposing the ultra-slow-roll initial condition on a plateau-like inflaton potential (such as Starobinsky inflation [3] and α-attractors [32,33]).In such a situation, the smooth transition to slow-roll attractor occurs automatically.
In addition, due to the scale-invariant power spectrum generated in the initial non-attractor phase, the primordial perturbations can be suppressed on large scales, which is favored by current CMB observations [31].Now we study the background evolution of this realistic model numerically, and then further check the analytical results above.Consider the following potential of Starobinsky inflation [3] which is very flat for large φ.In the slow-roll attractor, the field velocity satisfies φsr = −V /(3H).However, if inflation starts with a much larger velocity | φ| | φsr | on this very flat potential, initially it would be in the non-attractor phase.Solving the background equations numerically, we get the results shown in Figure 2. As we can see from the phase space diagram and the evolution of η, this realistic model indeed has a non-attractor initial phase, and then it will join the slow-roll attractor very quickly.
With this numerical solution, we can go back to do the full computation for the integral in (2.24), not only for the perturbation which exit the Hubble radius before the transition (non-attractor modes), but also for those small scale modes (slow-roll modes).The final bispectrum receives contributions from both terms in (2.24).The numerical result of local f NL as a function of k is shown in Figure 4.As we see, if we only consider one contribution, the local non-Gaussianity is O(1) for the non-attractor modes, and then it vanishes for the slow-roll modes.However, when we combine these two contributions together, they cancel each other and yield vanishing f NL even for non-attractor modes.This result confirms the analytical calculation above.In summary, both analytical and numerical calculation of the smooth transition process show that, there is a mysterious cancellation happening during this transition period.In the following subsection, we shall understand this cancellation in a more general way.

A more general analysis
To understand what is going on during a smooth transition, we present a more general analysis as follows.First of all, let us remind of the background equations Here the potential is required to have a slow-roll attractor, but for now we do not assume any slow-roll conditions.Then using the background equations and the Hubble "slow-roll" parameters defined in (2.2), the second and third order derivatives of the slow-roll potential can be exactly expressed as which respectively correspond to the inflaton mass and self-coupling.Note that these derivatives of the potential should be suppressed so that the slow-roll attractor is possible.Due to this requirement, some useful combinations of η and η, that we will soon encounter, should be much smaller than unity, even though η and η can be individually large during the non-attractor and transition stages.One consequence of this observation is the behaviour of the effective mass in the Mukhanov-Sasaki equation (2.3).As we mentioned in the last subsection, the coefficient ν 2 − 9/4 there is always small, even during the transition where η and η are big.Now we see this parameter is directly related to inflaton mass which does not care if inflation is in the attractor or not.With this knowledge, let us look at the cubic interaction term (2.10) again.In our in-in calculation above, one subtlety is caused by the evolution behaviour of Ṙ.We can remove it via integration by part, and express (2.10) as Since there is no more time derivative on R, the only important effect lies in the cubic coupling.
Here we are encouraged to introduce the effective coupling as (2.32) Again it looks like due to the drastic variation of η, these terms could be large during the transition.However, if we plug in our analytical and numerical solutions in the last section, this combination is shown to be negligible.Interestingly, they are also present in V , and can be written as Therefore the contribution from this term is always small, no matter how big η and η are during the transition.The presence of V is not a coincidence here.In the flat gauge, the operator which contributes to the cubic Lagrangian (2.10) comes from the self-interaction of field fluctuations (2.34) Taking the decoupling limit, we have (2.35) And after integration by parts, the self-interaction term exactly gives us the cubic term (2.10).In summary, for a smooth non-attractor to slow-roll transition, as long as we have a slow-roll potential (V is small), the non-Gaussianities would be always small.The cubic interaction term (2.10), which was previously thought to contribute sizable f NL in the instant transition approximation, actually never contributes in the realistic smooth transition case.However, this argument may not work in sharp transition cases.There the potential is unsmooth around the transition, which may yield large V .Furthermore, the unconventional behaviour of the mode function will add extra complications.These issues of the sharp transition will be addressed in the next subsection.

Non-Gaussianity in a sharp transition
As we discussed previously, the background of sharp transition differs from the smooth transition case.Now we come to study the effect of a sharp transition on the evolution of perturbations, especially on the local non-Gaussianity.The sharp transition corresponds to the case where the second term in the potential (2.14) is also important.In this subsection, we shall study the case The form of the potential (2.14) can be invalid after the inflaton field evolves to sufficiently large distances from the transition point.But we can assume that the slow-roll limit is already reached before that happens so that we do not need to keep track of perturbations any more.
First of all, unlike the smooth transition case, the behaviour of the mode function in the sharp transition is more complicated.If we look at the Mukhanov-Sasaki equation and the index ν in (2.4), the analytical solution of the sharp transition (2.19) and (2.20) still gives us ν 2 = 9/4 − 3η V .However, due to the sudden change of η at the transition time τ e , one cannot simply continue using the initial mode function (2.5) after τ e .Here when the transition happens, the matching condition requires the mode function and its first derivative to be continuous, i.e.R(τ e − ) = R(τ e + ) and R (τ e − ) = R (τ e + ).This gives us the following behaviour of curvature perturbation after where

.38)
We can easily check the long wavelength behaviour of the mode function after This solution satisfies the super-horizon EoM: R + (3 + η)H Ṙ = 0.At the time τ 0 of the slow-roll stage, we get the freezed amplitude For large values of |h| the above relation reduces to R k (τ 0 ) H √ 4 ek 3 which is similar to the mode function at the transition time τ e .Thus, for |h| 1, the final power spectrum does not change much by the transition and we have P R H 2 8π 2 e .Therefore, we expect to recover the previously calculated non-Gaussianity f NL = 5/2 in the |h| 1 limit where the mode function is assumed to freeze instantly after transition.We will confirm this expectation explicitly below.Note also that, in the h = −6 case with only instant transition, the super-horizon modes still evolve from τ e to τ 0 , as can be seen from (2.39).This shows that a realistic instant transition to the slow-roll evolution (which corresponds to h = −6) does not imply an instant freezing of the mode function, thus we do not expect to recover f NL = 5/2 after this transition.On the other hand, for |h| 1, the adiabatic limit is reached instantly and the mode function freezes out immediately whereas the background evolution experiences a transition period before it relaxes to the slow-roll dynamics.
For the sharp transition, the in-in integral in the bispectrum (2.11) can be divided into two nontrivial pieces : one is the contribution from the instant transition at τ e , where η can be approximated by the step function as in (2.22); and the other one is the relaxation period from τ e to τ 0 , which is described by the analytical solution in Section 2.2.
For the first piece, the integral goes from τ e − to τ e + .At τ e − , the mode function is described by (2.5), and η = −6.At τ e + , the mode function is given by (2.36), and η = −6 − h.Thus taking the squeezed limit and focusing on perturbation modes which exit the Hubble radius during the non-attractor phase, we can write this contribution to the bispectrum as Then via (2.22), the integral above yields where in the last step we took an integration by parts to reduce the integral to boundary terms.
The second part of the integral corresponds to the relaxation process after τ e .Substituting the mode function (2.36) and (2.37) into (2.11), its contribution to the squeezed bispectrum is given by Using the analytical solution during the relaxation (2.16) and (2.17), the above integral becomes (2.44) where we used 0 V which holds in the sharp transition with Adding these two contributions together, we get (2.45) As we see, the amplitude of local non-Gaussianity is mainly determined by the h parameter in sharp transition case.For |h| 1, it yields the maximum value f NL 5/2, which recovers the result in the initial non-attractor phase.For the instant transition (h = −6), we get a reduced value f NL = 5/8.In general, the sharp transition suppresses the amount of local non-Gaussianity generated during the non-attractor phase.The extremal case is h → 0, where we have negligible contribution 3  5 f NL = −hη V /12, similar to the smooth transition result.Concluding the subsection, we remark that the sharp transition of non-attractor inflation is different from the inflationary feature models, where due to the kink or step on the potential, one may have a short non-slow-roll period which connects two slow-roll stages before and after the local feature.Since initially inflation is on the slow-roll attractor, long wavelength modes will remain constant during the non-slow-roll period.Therefore these feature models cannot result in nontrivial local non-Gaussianity for large scale perturbations, and the consistency relation is still valid.The reason is that once the mode is frozen in the adiabatic limit it remains so regardless of what may happen after, because a constant is a solution of the EoM for the super-horizon mode function.However, in the sharp transition here, because of the initial non-attractor phase, long wavelength modes may continue to evolve on super-horizon scales.As a consequence, local non-Gaussianity can be modified on large scales due to the transition.
Related to this issue, it is also known that the presence of sharp feature on potential will generate scale-dependent oscillatory signals in power spectrum and non-Gaussianities (See e.g.[12] for a review).The argument is very general and should apply here as well.However, this sinusoidal oscillation starts to appear around the scale k ∼ 1/τ e and has a wavelength ∆k ∼ 1/τ e .So they appear at much shorter scales than what we are interested in in this paper.

δN calculation
The δN formalism [34][35][36][37][38][39][40] is a simple and intuitive approach to the non-linear behaviour of curvature perturbations.Based on the separate universe assumption, it mainly captures the super-horizon effects of the perturbation modes, thus it just provides what we need for the calculation of local non-Gaussianity.For non-attractor inflation, one extra subtlety one should take care of is that the number of e-folds N does not only depend on the initial field value φ, but also on the initial field velocity π [17].In the following, via δN formalism we give a unified calculation of local non-Gaussianity that captures both smooth and sharp transition cases, and recovers the in-in results in Section 2.3 and 2.4 in two extreme limits.
For the non-attractor phase, the number of e-folds N can be easily worked out.As in Section 2.2, we set N = 0, φ = φ e and dφ/dN = π e at the end of the non-attractor phase, then the background equations ( 2 Next we can invert this solution and obtain the e-folds of the non-attractor phase in terms of the initial φ and π .47)where in the second equality we used the following relation of the non-attractor phase (2.48) For the subsequent transition and slow-roll stages, the analytical solutions are already worked out in (2.16) and (2.17).Here we need to study the evolution until the end of the transition, where the slow-roll attractor is reached.Let us set N = N f and φ = φ f at that time.Then N f is big, and (2.16) yields the following approximation which gives us (2.50) In the second equality, we separate out the parts unrelated with initial condition (φ, π) as a constant.Note here, due to the relation (2.48), π e and also h are determined by the initial φ and π in the non-attractor phase.
Finally, the total e-folding number from the non-attractor phase to the slow-roll stage counted backward in time is given by (2.51) The δN formula is simply given by Since δφ is approximately constant on super-horizon scales, δπ is exponentially suppressed and thus can be neglected.As a result, from (2.51) we get where again we used the initial condition dependence of π e (φ, π) from (2.48).And the local non-Gaussianity directly follows This calculation is valid for both smooth transition (h → 0) and sharp transition (h = 0).As we discussed previously, η V is always small, but |h| can be large for the sharp transition.
Thus similar with the in-in result (2.45), when |h| 1, we recover f NL = 5/2.For the smooth transition or sharp transition with small h, we get f NL −5η V /6 = 5η 0 /12, where η 0 is the second Hubble slow-roll parameter in the slow-roll stage.Note that this also agrees with the full in-in calculation.In such cases, the in-in result from cubic interaction term (2.10) is sub-dominant, and thus the leading contribution comes from the field redefinition (2.8), which yields the same result as above.
We close this section by some concluding remarks.It is interesting to discuss the implications of our results on the consistency relation violation in canonical non-attractor inflation.As we know, the power spectrum generated in the non-attractor phase is scaleinvariant with n s − 1 = 0.However, the final result (2.55) yields nonzero value for f NL after the transition.Even in the smooth transition case where f NL is slow-roll suppressed, we do not have f NL = 5  12 (1 − n s ).Therefore, the consistency relation is still violated in the non-attractor inflation with full consideration of the transition process.
It is also interesting to notice that f NL = 5/2 is the maximum non-Gaussianity that one can obtain from such a model irrespective of the details of the transition period.This upper bound holds true even if one considers either a bump potential (where the slope at the transition point is negative) or a step potential (where the potential is discontinuous).That is, although the final f NL as a function of parameters is clearly different for these cases, its value cannot exceed the f NL that is generated purely during the non-attractor phase.In terms of δN formalism, there can be two contributions to the final non-Gaussianity: one from the non-attractor e-folds N i in (2.51), another one from N f .When N i terms are the dominant contribution in the δN expansion (2.53), we recover the O(1) non-Gaussianity of the non-attractor phase.In the opposite limit, where N f terms are dominating, it turns out that the non-Gaussianity is small.This is an interesting observation without rigorous proof.But we remark that the N f part of the evolution is basically the case with non-slow-roll initial condition on a slow-roll potential, which is generically expected to yield small non-Gaussianity, as we argued in Section 2.3.2.Thus if N f terms dominate in δN expansion (2.53), we expect a slow-roll suppressed f NL .As a consequence, the upper bound is given by the non-attractor result f NL = 5/2 when N i terms contribute.
Having discussed the difficulties in obtaining larger than 5/2 non-Gaussianity, it is worth mentioning that it is not impossible in canonical non-attractor models.As a concrete example, consider an upward step in an otherwise constant potential.If the initial velocity of the inflaton is sufficiently high, it climbs up the step during which the non-Gaussianity blows up momentarily but rapidly relaxes to the f NL = 5/2 afterwards.Thus, if we cut the potential right after the step and attach it to the slow-roll potential with the condition that the adiabatic limit is reached quickly (i.e. the case with |h| 1 as discussed above), we can obtain arbitrarily large non-Gaussianity.Although this is a very restrictive and fine-tuned scenario, it nevertheless shows that there is no physical reason to believe that large local non-Gaussianity cannot be obtained from a canonical, non-attractor, single field model of inflation.
inflation due to the fine tuning of its initial conditions, we study the perturbations only in a specific limit where the analytic calculation is still tractable.
In this model the non-attractor inflation is realized by a k-essence field with the following Lagrangian where X = − 1 2 (∂φ) 2 , and α, M , V 0 , v, β are free parameters.In this model, the sound speed c s is given by The following variables are also defined here for future reference: To the best of our knowledge, so far this is the only model which can give us f local NL 1 in single-field inflation with Bunch-Davies initial state.In this section, we will give a detailed analysis for the transition process in this model, and perform the full calculation to test whether large non-Gaussianity remains or not.

Background evolution of k-essence non-attractor model
First of all, let us focus on the background dynamics of this model.The equation of motion for inflaton can be written as From the above equation and (3.2) we can see that one important parameter here for the evolution is the ratio X/M 4 .For X M 4 , this model is non-canonical with c 2 s 1/(2α − 1); but for X M 4 , it returns to the canonical case.In this model initially the inflaton field climbs up the hilltop potential, with the kinetic energy dominated by the non-canonical term.Later on, as X decreases dramatically in the non-attractor phase, the system would go from the non-canonical regime to the canonical regime.
For k-essence field, the slow-roll parameters are expressed as As we know, a non-attractor phase happens when ∝ a −6 and η −6.In the original papers [19,20], an ansatz φ(t) ∝ a κ was used to get the initial non-attractor stage.This was achieved by letting the V φ term compete with the φ and φ terms in the equation of motion (3.5).And the following conditions for parameter choices are required However, it is still not clear how the system transits to the attractor phase in details.In the following we perform a full numerical study of the "non-attractor to slow-roll" transition in the k-essence model.Before that, we summarize the generic behaviour for the evolution first: The main results of the numerical solution are shown in Figure 5.At the beginning, since the potential is tuned to accommodate with the ansatz as shown in (3.8), inflation occurs in the phase with η = −6, while X M 4 gives a small sound speed.We call this initial stage the non-attractor I. Then as the inflaton approaches the hilltop, the V φ term in (3.5) becomes subdominant, and thus the equation of motion becomes φ + 3Hc 2 s φ 0, which according to (3.7) yields η = −3(c 2 s + 1) . (3.9) Since the inflaton field is still non-canonical (c s 1), we have η −3.We dub this period as the non-attractor II phase.Next, X continues decreasing and becomes smaller than M 4 , then the canonical term in P (X, φ) begins to dominate the kinetic energy of inflaton.After that, the scalar field becomes canonical, and we call this moment the canonical transition.And from (3.9), we see the system goes to the canonical non-attractor regime with η = −6.This stage has the same behaviour with the canonical non-attractor model, and is called non-attractor III here.Finally, the following transition to the slow-roll attractor is the same as what we discussed in Section 2. As we see, the "non-attractor to slow-roll" transition is much more complicated in the non-canonical model.One important feature is that there is also a canonical transition prior to the slow-roll attractor phase.This qualitative description is confirmed by the numerical analysis below.
Numerical Study.Following the choice of parameter values in [19,20], here we take α = 10, M = 5 × 10 −5 , V 0 = 6.25 × 10 −4 , while v and β are given by the relation in (3.8).Initially inflaton field is set to roll up the hilltop potential from φ i = 2 × 10 −6 .Then via varying the initial field velocity, we find different transition behaviours.The numerical solutions of background dynamics are shown here.Figure 6 gives us the phase space diagram.In Figure 7, we focus on the evolution of two parameters: the second slow-roll parameter η, which is important for the non-attractor behaviour, and the sound speed c s , which tells if inflaton field is canonical or not.
From these figures, we can see a generic pattern for the transition process: after the non-attractor I stage (η −6), inflation first enters the non-attractor II phase (η −3), and later as shown by the evolution of c s , the canonical transition happens.Here we introduce a critical field velocity φc , for which inflaton can just reach the top of the potential and will stay there forever.Then accordingly the numerical analysis can be classified into three representative cases: • Undershoot (blue curves).This corresponds to the case where the initial field velocity is smaller than | φc |.After a very short non-attractor II phase, in the canonical nonattractor regime, inflaton stops somewhere before reaching the top of the potential and then rolls backward to initiate the slow-roll phase.At the turning point, since φ = 0, we have = 0 and η = ∞.
• Critical case (red curves).The initial field velocity is set to be the critical value.In this case, after the canonical transition, the system reaches an eternal non-attractor stage with η = −6 and c s = 1.
• Overshoot (orange curves).This is the case where the initial field velocity is larger than | φc |.As we see, here the non-attractor I stage is very short, while the non-attractor II phase lasts for a longer time, during which the inflaton field rolls over the top of the potential.After this, inflation goes into the canonical non-attractor regime and then transits to the slow-roll stage as we discussed before.
In these three cases, only undershoot and overshoot can give us successful "non-attractor to slow-roll" transition.Although the details can be very different, both these two numerical results verify the evolution in Figure 5 and the qualitative description there, i.e. in these non-canonical models the canonical transition always occurs before the relaxation to slow-roll.This holds true at least for our choice of parameters which are consistent with [19].It would be interesting to see whether it is also true for other values of the parameters; however, we do not go further in this direction here.In the following rough calculation of non-Gaussianity, we shall use this general transition behaviour as the basic setup and refer to these two cases (overshoot and undershoot) for details.

Non-Gaussianities
With the above background analysis, we are ready to study the primordial perturbations.At first glance, a full calculation could be very difficult, since the transition behaviour is quite complicated.Numerical calculation also faces a technical UV-convergence problem because the non-attractor phase is rather short.
However, the problem can be simplified if we focus on the generic pattern of the transition.As we see, the main difficulty comes from the occurrence of the non-attractor II phase, during which we have η = −3.If this period lasts for a long time (as in the overshoot case), we cannot get a scale-invariant power spectrum for curvature modes that are leaving the horizon during this period.This can be interesting for the research of features in the primordial perturbations, but in this paper we keep focusing on the behaviour of non-Gaussianity during the transition.And for the analysis of the bispectrum, it is the canonical transition that plays a crucial role here.
Therefore we propose the following limit case for analytical study: The non-attractor II phase is so short such that its effect can be neglected.In this approximation, before the slow-roll attractor, η can always be seen as −6 and the canonical transition occurs at some time in this stage.In principle this does not agree with the numerical results since it breaks the relation (3.9), but it can be seen as an approximated description of the undershoot case.
Based on the qualitative analysis above, next we focus on the modes which exit the Hubble radius during the non-attractor I stage, and do the back-of-the-envelope estimates for the non-Gaussianities.The starting point is the cubic action for a general k-essence field in the comoving gauge [41,42].Since 1 always holds true during the whole transition process, again we can focus on the decoupling limit with only three operators left The coefficient of the Ṙ3 term is given by where (3.2), (3.3), (3.4), (3.6) and (3.7) are used for the second equality.Before the canonical transition we have Ξ = 2(c 2 s − 1)/3c 2 s .This coefficient and the one for the second term in (3.10) both vanish after the canonical transition.At the same time, the following field redefinition is also considered in [20] which can give large non-Gaussianity in the non-attractor I phase.However, since this term should be evaluated in the slow-roll stage where η 0 and Ṙ 0, its contribution can be neglected.For the last operator in (3.10), again we re-express it via integration by part Plugging in the numerical solution, we confirm that the effective coupling here is also negligible during the transitions as in the canonical case.
The difference from the canonical case arises due to a couple of interaction terms in the Lagrangian that are unique for the non-canonical models.Let us estimate their contributions.The first term in (3.10) gives the bispectrum of where τ c is the conformal time at the canonical transition, and τ 0 is the one at the beginning of the slow-roll phase.Since Ξ vanishes after the canonical transition, the in-in integral stops at τ c .Another subtlety here is the mode function R k .Since there is a sudden change of c s around the transition, in principle one has to use the general slow-roll formalism to solve its behaviour, taking into account the discontinuity around the canonical transition.However, since the integral above vanishes right after the canonical transition, the mode function after transition becomes irrelevant for that integral; and it is easy to check that it does not affect the prefactors R k i (τ 0 ) in (3.14) at leading order either.Therefore, for a rough estimate, here we take the following zeroth order approximation Since we mainly care about the modes crossing the Hubble radius during the initial nonattractor phase, we have −kτ 0 −kτ c 1. Meanwhile in this limited case we assume η = −6 before the time τ 0 , which means for this whole period (τ ) = 0 τ 6 /τ 6 0 .As a result, the bispectrum becomes which is in the local shape.As we see, when τ 0 = τ c it returns to the previous result in [20].However, if the canonical transition occurs ∆N e-folds before the slow-roll phase, (τ 0 /τ c ) 6 would give a suppression factor ∼ e −6∆N .Correspondingly in the squeezed limit, we get the following amplitude of non-Gaussianity This suppression is caused by the super-Hubble evolution of the curvature perturbation after the canonical transition.Since R keeps growing until the end of the non-attractor phase, the difference between R(τ c ) and R(τ 0 ) yields the suppression factor above.With a similar procedure, the second term in (3.10) gives where in the last step we ignored the e −6∆N suppression term.Again, the above result agrees with the one in [20] when τ 0 = τ c .In general, the duration of the non-attractor stage after the canonical transition can be ∆N ∼ O(1), thus the large non-Gaussianity generated in the non-attractor stage can be suppressed a lot.Summing up the leading terms of two contributions above, we get the following overall amplitude This estimate shows us how the non-Gaussianity generated in the initial non-attractor stage is suppressed after the canonical transition.Notice that the sound speed is determined by the model parameters, while the duration of the non-attractor III stage is related to the choice of initial conditions, thus c s and ∆N are two independent parameters.Thus, we conclude that, it is still possible to have large non-Gaussianity in single field inflation.

Conclusion and discussion
In this paper, we investigated the production of primordial non-Gaussianities from models of non-attractor inflation.We revisited various non-attractor models constructed in the literature in order to understand the evolution of large local non-Gaussianity when the models undergo the transition from the non-attractor phase to slow-roll phase.The purpose of this study is less of trying to present these fine-tuned toy-models as phenomenological candidates for data fitting, rather trying to understand more precisely the physical implications of Maldacena's single field consistency relation and various counter-examples that have been constructed.
Comparing with previous studies, we pay special attention to the transition period from the non-attractor phase to the conventional slow-roll phase.Such a transition is necessary for these models to have sufficient efolds or have the correct amplitude of density perturbations.We considered two types of non-attractor inflation, which are driven by a canonical scalar field and a non-canonical k-essence field, respectively.
For models with canonical kinetic terms, we consider two different evolutionary processes after the non-attractor phase: smooth transition and sharp transition.Through the calculation of both in-in and δN formalism, we find that a full consideration of the transition process generically suppresses the local non-Gaussianity generated in the non-attractor phase, but Maldacena's consistency condition is still violated.In the smooth transition, the super-horizon modes continue evolving after the non-attractor phase, and the O(1) non-Gaussian signals are completely erased during the transition period and the final f NL at the end of inflation is slow-roll suppressed.Meanwhile for sharp transition, the final amplitude of the local non-Gaussianity generated in the non-attractor phase depends on the details of the transition process.In the extremal cases where the curvature perturbation freezes immediately right after the non-attractor phase, we get the maximum possibility of local non-Gaussianity, which recovers the original result in the non-attractor phase f NL = 5/2.
For models with non-canonical kinetic terms, although similar situation applies to one of the terms in the Lagrangian, the non-Gaussianities coming from two other terms, which are unique to non-canonical models, survive.Nonetheless, our rough estimations of this case show that the effect of smooth transition is still non-negligible.In addition to the contribution ∼ 1/c 2 s obtained in the previous studies, the transition period contributes to an extra suppression factor due to mode evolution outside the horizon during the transition phase.Since these two contributions are independent of each other, the conclusion, that the large local non-Gaussianity can be obtained in such single field models, remains the same; but the expression of f NL should be revised.
As a final remark, we note that, recently, Ref. [43] argued that the O(1) local bispectrum generated from the canonical non-attractor inflation model, as calculated in Ref. [17], is not locally observable.The study of Ref. [43] focuses on the non-attractor phase.Here we will not analyze their argument in detail which is beyond the scope of this paper.For our purpose, we simply point out that one of the main differences between their work and ours is that we have analyzed in details the subsequent transition process from the non-attractor phase to the standard single field slow-roll inflation, in order to be able to discuss the observability at all.As we have concluded, the final f NL can range anywhere between zero and a value much larger than 1.If the value of f NL is much larger than 1 − n s , these local bispectra should be in principle observable.At the reheating surface, these local bispectra are indistinguishable from those arising from models in which we replace the single field non-attractor phase with a multifield phase and use the multifield phase to generate the same amount of primordial local bispectra.

Figure 2 .
Figure 2. Smooth transition.Left Panel: the phase space diagram of non-attractor to slow-roll transition on a plateau-like potential.Right Panel: the evolution of η parameter during the transition.

Figure 3 .
Figure 3. Sharp transition.Left Panel: the phase space diagram of sharp transition for three different cases.Right Panel: the evolution of η parameter during the sharp transition.

Figure 4 .
Figure 4.The cancellation in the in-in integral (2.24).

3 ,
(3.18)   Therefore, it is still in the local form and the final amplitude of non-Gaussianity is given by