# Review of Local Non-Gaussianity from Multifield Inflation

Authors: Christian T. Byrnes; Ki-Young Choi
Journal: Advances in Astronomy
 (2010)

---

## Abstract

We review models which generate a large non-Gaussianity of the local form. We first briefly
consider three models which generate the non-Gaussianity either at or after the end of inflation; the
curvaton scenario, modulated (p)reheating, and an inhomogeneous end of inflation. We then focus
on ways of generating the non-Gaussianity during inflation. We derive general conditions which
a product or sum separable potential must satisfy in order to generate a large local bispectrum
during slow-roll inflation. As an application, we consider two-field hybrid inflation. We then derive
a formalism not based on slow roll which can be applied to models in which the slow-roll parameters
become large before inflation ends. An exactly soluble two-field model is given in which this
happens. Finally, we also consider further non-Gaussian observables, a scale dependence offNL and
the trispectrum.


---

## Body

## 1. Introduction

There are many models of the universe which can predict a large non-Gaussianity. However the predicted amplitude and the shape of the non-Gaussianity are different among different classes of models. One category is those which generate the non-Gaussianity due to nontrivial classical dynamics on superhorizon scales. These models predict the shape of the bispectrum to be of the so-called “local type”, which can be expressed as an expansion of the Bardeen potential [1] (1)Φ(x)=ΦL(x)+fNL(ΦL2(x)-〈ΦL2(x)〉),
where Φ is the curvature perturbation on a Newtonian slice and ΦL is its linear and Gaussian part. 〈ΦL2(x)〉 denotes the ensemble average in a statistically homogeneous distribution. The current limit on the local type of the nonlinearity parameter fNL from seven years of WMAP data [2] is -10<fNL<74 at the 95% confidence level. Constraints are expected to improve rapidly and significantly, first with Planck data and later using large scale structure data, see the recent reviews [3–5]. The Bardeen potential is related to the primordial curvature perturbation of ζ on large scales and in the matter dominated era by Φ=(3/5)ζ.The curvature perturbation at horizon exit is determined by the classical perturbations of the scalar fields,δϕi(x). The subsequent evolution of ζ can be conveniently described by the δN formalism [6–11]. The curvature perturbation is given by up to quadratic terms [11] (2)ζ=δN=∑IN,IδφI*+12∑IJN,IJδφI*δφJ*+⋯,
where N(x,t) is the e-folding number evaluated in an unperturbed Universe, from the epoch of horizon exit to later epoch of uniform energy density hypersurface (for an extension to include gradient terms, see [12]). The power spectrum 𝒫ζ and the bispectrum Bζ are defined by
(3)〈ζk1ζk2〉≡(2π)3δ3(k1+k2)2π2k13𝒫ζ(k1),〈ζk1ζk2ζk3〉≡(2π)3δ3(k1+k2+k3)Bζ(k1,k2,k3).
From this, we can define the observable quantities, the spectral index, the tensor-to-scalar ratio, and the nonlinearity parameter (4)nζ-1≡∂log𝒫ζ∂logk,(5)r=𝒫T𝒫ζ=8𝒫*MP2𝒫ζ,(6)fNL=56k13k23k33k13+k23+k33Bζ(k1,k2,k3)4π4𝒫ζ2,
where 𝒫T=8𝒫*/MP2=8H*2/(4π2MP2) is the power spectrum of the tensor metric fluctuations. It is well known that single-field inflation does not lead to a detectably large non-Gaussianity, in fact fNL is suppressed by slow-roll parameters [13]. Observably large non-Gaussianity can be obtained by breaking the slow-roll conditions during inflation [14], using extended kinetic terms [15–25], see also the reviews [26, 27], or going beyond models of single-field inflation [11, 28–34]It is natural to consider multiple scalar field since they are ubiquitous in many beyond the standard model of particle physics, such as supersymmetry and string theory. These scalar fields generate nonadiabatic perturbations during inflation and change the evolution of the curvature perturbation after horizon exit. The residual isocurvature perturbation may be present in the primordial density fluctuation and can be correlated with the curvature perturbation or may be responsible for an observably large non-Gaussianity in the cosmic microwave background and large scale structure, for observational limits on isocurvature perturbations see [1, 35, 36]. In this paper, we will only consider models with adiabatic primordial perturbations, in which the isocurvature perturbation present during inflation is converted into an adiabatic perturbation. We also neglect the secondary non-Gaussianities generated at later times, for example see [37, 38].There are popular multifield models [39] which may generate observably large non-Gaussianity. These include the curvaton scenario, modulated (p)reheating, and an inhomogeneous end of inflation, see Section 2. In these scenarios, large non-Gaussianity is generated either by the means of ending inflation, or after inflation. It was shown recently that it also possible to generate large non-Gaussianity during the evolution of slow-roll multifield inflation, see Section 3.All of these models generate the large non-Gaussianity after horizon exit, such as after reheating, at the end of inflation, at the phase transition, or during inflation after horizon exit and involves the perturbation of the nonadiabatic mode. Therefore, the non-Gaussianity of these models is of the local type which is distinguishable from other shapes of non-Gaussianity (for a list of possibilities see, e.g., [40–42]), in which the non-Gaussianity is generated intrinsically from the quantum fluctuations, or during horizon exit.In Section2, we summarise the aforementioned three models, which are popular methods of generating a large non-Gaussianity. In Section 3, we review the possibility of generating a large non-Gaussianity from multifield slow-roll inflation, and in Section 4, we consider hybrid inflation with two inflaton fields as an application. Then in Section 5, we discuss multifield models of inflation without assuming the slow-roll conditions and present an exact solution. Non-Gaussian observables beyond fNL, such as its scale dependence and the trispectrum are introduced in Section 6. Finally, we conclude in Section 7.


## 2. Summary of Popular Models Generating a Large Local Non-Gaussianity

Multifield inflationary models [39] can lead to an observable signature of non-Gaussianity. These include models in which the large non-Gaussianity is generated either by the means of ending inflation, or after inflation. We review the three most popular of these models in this section, for a discussion on how these models are related see [43]. In the next section, we review the possibility of generating large non-Gaussianity during multifield slow-roll inflation. Later in Section 5, we discuss multifield models of inflation without assuming the slow-roll conditions. This does not exhaust all of the possible ways of generating a large local non-Gaussianity from inflation, see also the reviews [44–47] and the papers [48, 49]. Non-Guassianity of the local form can also be generated in the ekpyrotic scenario [50, 51], although in the simplest case fNL is large and negative [52], which is observationally ruled out.


### 2.1. Curvaton Scenario

In this scenario, there is a light weakly-interacting inhomogeneous scalar field, the curvatonχ, in addition to the inflaton field, ϕ. The additional scalar field is completely subdominant during inflation however, it can dominate the energy density of the Universe later since the time-averaged equation of state becomes that of the pressureless matter (P=0), whose energy density decreases more slowly than that of the radiation which the inflaton field has decayed into. After decay, the inhomogeneity of the scalar field leads to the density perturbation of radiation [53, 54]. This is the basis of the curvaton scenario [55–59] which liberates the inflaton field from being required to generate the observed spectrum of perturbations (although there is instead an upper bound on their allowed magnitude).The significant non-Gaussianity of the local type can be generated in the curvaton model [60]. A study of non-Gaussianity in the curvaton scenario using second-order cosmological perturbation theory was done by Bartolo et al. [31] with a sudden decay approximation and also by Malik and lyth [32] using a fully numeric approach, which accounted for both the sudden decay and nonsudden decay approximations. Lyth and Rodríguez [11] used the δN formalism including the quadratic terms (second order terms in (2)) to calculate the non-Gaussianity. The nonlinear generalisation of the curvature perturbation and a numerical study was done to give a full probability distribution function in [61]. This simple curvaton model was extended into the mixed curvaton-inflaton scenario in which fluctuations of both the inflaton and a curvaton contribute to the primordial density perturbation studying the isocurvature perturbation in [62]. The multiple curvaton fields which contribute to the primordial density perturbation was studied in [63–65]. Observational constraints on curvaton models with isocurvature perturbation and non-Gaussianity can be found in [66–68].The generation of large non-Gaussianity in the curvaton scenario can be easily understood using the result of the nonlinearδN formalism [61, 62]. Assuming the sudden decay of the curvaton field on the H=Γ uniform energy-density hypersurface, that is, when the local Hubble rate equals the decay rate for the curvaton, leads to a nonlinear relation between the local curvaton density and the radiation density before and after the decay (7)Ωre4(ζinf-ζ)+Ωχe3(ζχ-ζ)=1,
where Ωr+Ωχ=1. Here, ζ is the primordial curvature perturbation which remains constant on large scales in the radiation-dominated era after the curvaton decays and ζχ is the curvaton perturbation. Expanding this equation order by order yields up to second order [61] (8)ζ=(1-rχ)ζinf+rχζχ+rχ(1-rχ)(3+rχ)2(ζχ-ζinf)2+higher-orderterms,
where rχ≡3Ωχ/(4-Ωχ) at the time of the curvaton decay. The curvature perturbation of curvaton field [61] (9)ζχ=23δ1χχ0-13(δ1χχ0)2+29(δ1χχ0)3,
is almost constant before the curvaton decays from the start of curvaton oscillation. Note that the intrinsic non-Gaussianity in the curvaton field during oscillation is fNL,χ=-5/4, which comes from the nonlinear relation between δ1χ and ζχ. After the curvaton decays into radiation, the curvature perturbation of radiation has the following non-Gaussianity parameters: [11, 31, 61]
(10)fNL=r̃2(1+r̃)2[54rχ-53-5rχ6],
where rχ is evaluated when the curvaton decays based on the sudden decay approximation but gives good agreement with a full numerical study [32, 61], and r̃ is defined as (11)r̃≡|rχζχ|2|(1-rχ)ζinf|2.
In the limit of the pure curvaton scenario, ζinf→0, it recovers the usual normalization for fNL, that is, r̃2/(1+r̃)2→1. Therefore, the large fNL is obtained for small rχ. This large fNL, in spite of the order of unity non-Gaussianity in ζχ, has the origin in the nonlinear relation between ζχ (or δ1χ) and ζ in (8).In fact (7), and thus (8), can be used “at any time” before the curvaton decay with rχ evaluated at that time to find the total curvature perturbation at that moment, as well as at the epoch of curvaton decay. This enables us to understand the evolution of the curvature perturbation ζ on the uniform energy density hypersurface before the curvaton decay. Initially when the primordial radiation was generated, for example, from the inflaton decay, rχ is negligibly small so we cannot ignore the perturbations of the inflaton field. At this early time, ζ=ζinf and fNL=0 since r̃∝rχ2 which makes fNL negligible in (10). As time goes on, rχ increases while ζinf and ζχ are almost constant since there is no energy flow between two fluids satisfying adiabatic condition except around the decay time, thus the curvature perturbation evolves gradually. If the curvaton survives long enough, rχ reaches a value with which the rχζχ dominates the first term, ζinf, and the curvaton perturbation is then responsible for the primordial inhomogeneity of the Universe. The large non-Gaussianity is possible only when rχ at the decay time is small that is, the energy density of curvaton is subdominant compared to that of radiation. This is because when rχ=1, the curvaton energy density dominates over the radiation component before decay, then ζ is effectively the curvature perturbation of a single field (curvaton), ζχ, and therefore non-Gaussianity fNL=fNL,χ=-5/4. We will also see this phenomena in the multifield slow-roll inflation later, in Section 3.


### 2.2. Modulated (P)reheating

After inflation, the energy density in the inflaton field must be transferred into radiation. In the simplest case of adiabatic perturbations, this process does not affect the primordial curvature perturbation on scales which are observable today, because these scale were much larger than the horizon at the time of reheating. However, in a spirit similar to the curvaton model, there may be a subdominant light scalar field present during inflation which modulates the efficiency of reheating. This makes the efficiency of reheating a spatially dependent process. The quasiscale invariant perturbations in this field, which during inflation are an isocurvature perturbation, may be converted into the primordial curvature perturbation during this process. For a review of reheating after inflation, see for example [69, 70].As a simple illustration of this, in the “old” models of perturbative reheating, the decay rate of the inflaton is given byΓ~λ2m, where λ is the strength of the coupling between the inflaton and reheating fields and m is the mass of the inflaton. Provided that Γ is much less than the Hubble rate at the end of inflation, so that reheating takes place slowly, the reheating temperature is given by
(12)Treh~ΓMP~λmMP.
If the coupling strength is a function of the local value of a second scalar field χ then this will give rise to a perturbation in the decay rate of the inflaton field, and thus in the reheating temperature which is responsible for the density perturbation after reheating. Provided that the χ field is perturbed, the curvature perturbation can be written as (13)ζ=αδΓΓ,
where α is a parameter depending on the ratio of Γ to the Hubble expansion at the end of inflation, whose value increases with decreasing Γ/H and leads to 1/6 in the limit Γ≪Hinf [33, 71]. Hence, perturbations of the light scalar field χ can be imprinted into the radiation temperature through the reheating process [71, 72]. In a similar way, the mass of the decaying particle may be modulated [73]. The above equation can be extended to second order, which is required to calculate fNL. In the simple quadratic case that Γ=Γ0+Γ1(χ/χ¯)2=Γ0+Γ1(1+δχ/χ¯)2, where χ¯ and δχ are the homogeneous background part and the perturbation of the field χ, respectively, the result is [33, 74–76]
(14)fNL=512Γ0+Γ1Γ1α.
Hence, we can see that the non-Gaussianity can easily be larger than unity in this model by two sources. One is due to the small value of α which occurs when the decay rate of the inflation is not much smaller than the Hubble constant at the end of inflation. The other is when Γ1 is small in which χ field only controls part of the decay channels of inflaton field.It is now generally considered that a period of rapid and highly nonlinear preheating preceded the perturbative reheating phase [69]. In this case, one can instead consider modulated preheating, see, for example, [77–79]. A particularly rapid form of preheating, known as instant preheating [80], has attracted more attention in this area, as it is possible to find analytic estimates [81–85]. It is possible to generate the primordial curvature perturbation in these models by having either a weak [81] or a strong [83] symmetry breaking of the two-field potential during inflation which affects the efficiency of instant preheating, or by having two inflaton fields which are coupled with different strengths to the preheat field [82]. However, all of the above models of modulated instant preheating are only allowed, even at the linear level in perturbations, for a small corner of parameter space and even then only for a special choice of initial conditions [85].


### 2.3. Inhomogeneous End of Inflation

An inhomogeneous phase transition can occur also at the end of inflation between inflationary and non-inflationary expansion due to the condition of how inflation ends [34, 86]. In the case of multicomponent inflation, there is a family of inflationary trajectories. If the family of trajectories is not perpendicular to the surface on which inflation ends then there is a possibility that inflation ends on a slice which is not of uniform density and the generation of the curvature perturbation at the end of inflation is expected [34] (This was also discussed in [87–89] using the linear relation between two fields at the end of inflation and the non-Gaussianity is transferred from that of the light field which is generated from its nonlinear self-coupling). In single field inflation, the energy density is determined solely by the single field ϕ, and the end of inflation is parameterised by the unique value ϕe. If there is another field σ, then ϕe(σ) will depend on position through the perturbation δσ(x). As a result, the change of e-folding number from a spacetime slice of uniform energy density just before the end of inflation to a spacetime slice of uniform density just after the end of inflation will have a perturbation δNe=ζe. This ζe is the curvature perturbation generated by the end of inflation and can dominate over the contribution from inflation [34].Using the perturbation ofϕe(σ)(15)δϕe=ϕe′δσ+12ϕe′′(δσ)2,
the curvature perturbation generated at the end of inflation is [34]
(16)ζe=Ne′δϕe+12Ne′′(δϕe)2=Ne′ϕe′δσ+12[2Ne′′ϕe′2+ne′ϕe′′](δσ)2.Whenζe dominates the perturbation generated during the inflation, the power spectrum is given by
(17)𝒫e=ϕe′22ϵe(H*2π)2
The non-Gaussianity can also be generated at the end of inflation since the perturbation at the end of inflation is not required to be Gaussian. Thus, we expect that the observable large non-Gaussianity can be generated in this process. From (16) the non-Gaussianity is [34]
(18)fNL=-53ϵe2ϕe′′(ϕe′)2.
This was applied to the simple example of extended hybrid inflation using sudden end approximation [90, 91].Recently Sasaki [92], Naruko and Sasaki [93]  calculated analytically the curvature perturbation and the non-Gaussianity in the multi-brid inflation models from horizon exit to the end of inflation considering the general couplings to the waterfall field. In general, we cannot just ignore the curvature perturbation generated during inflation. Furthermore, that is also affected by the end of inflation. A deeper understanding of non-Gaussianity both from the evolution during multifield hybrid inflation and effects from the end of inflation is considered by Byrnes et al. [94, 95] which is summarised in Section 3.


## 2.1. Curvaton Scenario

In this scenario, there is a light weakly-interacting inhomogeneous scalar field, the curvatonχ, in addition to the inflaton field, ϕ. The additional scalar field is completely subdominant during inflation however, it can dominate the energy density of the Universe later since the time-averaged equation of state becomes that of the pressureless matter (P=0), whose energy density decreases more slowly than that of the radiation which the inflaton field has decayed into. After decay, the inhomogeneity of the scalar field leads to the density perturbation of radiation [53, 54]. This is the basis of the curvaton scenario [55–59] which liberates the inflaton field from being required to generate the observed spectrum of perturbations (although there is instead an upper bound on their allowed magnitude).The significant non-Gaussianity of the local type can be generated in the curvaton model [60]. A study of non-Gaussianity in the curvaton scenario using second-order cosmological perturbation theory was done by Bartolo et al. [31] with a sudden decay approximation and also by Malik and lyth [32] using a fully numeric approach, which accounted for both the sudden decay and nonsudden decay approximations. Lyth and Rodríguez [11] used the δN formalism including the quadratic terms (second order terms in (2)) to calculate the non-Gaussianity. The nonlinear generalisation of the curvature perturbation and a numerical study was done to give a full probability distribution function in [61]. This simple curvaton model was extended into the mixed curvaton-inflaton scenario in which fluctuations of both the inflaton and a curvaton contribute to the primordial density perturbation studying the isocurvature perturbation in [62]. The multiple curvaton fields which contribute to the primordial density perturbation was studied in [63–65]. Observational constraints on curvaton models with isocurvature perturbation and non-Gaussianity can be found in [66–68].The generation of large non-Gaussianity in the curvaton scenario can be easily understood using the result of the nonlinearδN formalism [61, 62]. Assuming the sudden decay of the curvaton field on the H=Γ uniform energy-density hypersurface, that is, when the local Hubble rate equals the decay rate for the curvaton, leads to a nonlinear relation between the local curvaton density and the radiation density before and after the decay (7)Ωre4(ζinf-ζ)+Ωχe3(ζχ-ζ)=1,
where Ωr+Ωχ=1. Here, ζ is the primordial curvature perturbation which remains constant on large scales in the radiation-dominated era after the curvaton decays and ζχ is the curvaton perturbation. Expanding this equation order by order yields up to second order [61] (8)ζ=(1-rχ)ζinf+rχζχ+rχ(1-rχ)(3+rχ)2(ζχ-ζinf)2+higher-orderterms,
where rχ≡3Ωχ/(4-Ωχ) at the time of the curvaton decay. The curvature perturbation of curvaton field [61] (9)ζχ=23δ1χχ0-13(δ1χχ0)2+29(δ1χχ0)3,
is almost constant before the curvaton decays from the start of curvaton oscillation. Note that the intrinsic non-Gaussianity in the curvaton field during oscillation is fNL,χ=-5/4, which comes from the nonlinear relation between δ1χ and ζχ. After the curvaton decays into radiation, the curvature perturbation of radiation has the following non-Gaussianity parameters: [11, 31, 61]
(10)fNL=r̃2(1+r̃)2[54rχ-53-5rχ6],
where rχ is evaluated when the curvaton decays based on the sudden decay approximation but gives good agreement with a full numerical study [32, 61], and r̃ is defined as (11)r̃≡|rχζχ|2|(1-rχ)ζinf|2.
In the limit of the pure curvaton scenario, ζinf→0, it recovers the usual normalization for fNL, that is, r̃2/(1+r̃)2→1. Therefore, the large fNL is obtained for small rχ. This large fNL, in spite of the order of unity non-Gaussianity in ζχ, has the origin in the nonlinear relation between ζχ (or δ1χ) and ζ in (8).In fact (7), and thus (8), can be used “at any time” before the curvaton decay with rχ evaluated at that time to find the total curvature perturbation at that moment, as well as at the epoch of curvaton decay. This enables us to understand the evolution of the curvature perturbation ζ on the uniform energy density hypersurface before the curvaton decay. Initially when the primordial radiation was generated, for example, from the inflaton decay, rχ is negligibly small so we cannot ignore the perturbations of the inflaton field. At this early time, ζ=ζinf and fNL=0 since r̃∝rχ2 which makes fNL negligible in (10). As time goes on, rχ increases while ζinf and ζχ are almost constant since there is no energy flow between two fluids satisfying adiabatic condition except around the decay time, thus the curvature perturbation evolves gradually. If the curvaton survives long enough, rχ reaches a value with which the rχζχ dominates the first term, ζinf, and the curvaton perturbation is then responsible for the primordial inhomogeneity of the Universe. The large non-Gaussianity is possible only when rχ at the decay time is small that is, the energy density of curvaton is subdominant compared to that of radiation. This is because when rχ=1, the curvaton energy density dominates over the radiation component before decay, then ζ is effectively the curvature perturbation of a single field (curvaton), ζχ, and therefore non-Gaussianity fNL=fNL,χ=-5/4. We will also see this phenomena in the multifield slow-roll inflation later, in Section 3.


## 2.2. Modulated (P)reheating

After inflation, the energy density in the inflaton field must be transferred into radiation. In the simplest case of adiabatic perturbations, this process does not affect the primordial curvature perturbation on scales which are observable today, because these scale were much larger than the horizon at the time of reheating. However, in a spirit similar to the curvaton model, there may be a subdominant light scalar field present during inflation which modulates the efficiency of reheating. This makes the efficiency of reheating a spatially dependent process. The quasiscale invariant perturbations in this field, which during inflation are an isocurvature perturbation, may be converted into the primordial curvature perturbation during this process. For a review of reheating after inflation, see for example [69, 70].As a simple illustration of this, in the “old” models of perturbative reheating, the decay rate of the inflaton is given byΓ~λ2m, where λ is the strength of the coupling between the inflaton and reheating fields and m is the mass of the inflaton. Provided that Γ is much less than the Hubble rate at the end of inflation, so that reheating takes place slowly, the reheating temperature is given by
(12)Treh~ΓMP~λmMP.
If the coupling strength is a function of the local value of a second scalar field χ then this will give rise to a perturbation in the decay rate of the inflaton field, and thus in the reheating temperature which is responsible for the density perturbation after reheating. Provided that the χ field is perturbed, the curvature perturbation can be written as (13)ζ=αδΓΓ,
where α is a parameter depending on the ratio of Γ to the Hubble expansion at the end of inflation, whose value increases with decreasing Γ/H and leads to 1/6 in the limit Γ≪Hinf [33, 71]. Hence, perturbations of the light scalar field χ can be imprinted into the radiation temperature through the reheating process [71, 72]. In a similar way, the mass of the decaying particle may be modulated [73]. The above equation can be extended to second order, which is required to calculate fNL. In the simple quadratic case that Γ=Γ0+Γ1(χ/χ¯)2=Γ0+Γ1(1+δχ/χ¯)2, where χ¯ and δχ are the homogeneous background part and the perturbation of the field χ, respectively, the result is [33, 74–76]
(14)fNL=512Γ0+Γ1Γ1α.
Hence, we can see that the non-Gaussianity can easily be larger than unity in this model by two sources. One is due to the small value of α which occurs when the decay rate of the inflation is not much smaller than the Hubble constant at the end of inflation. The other is when Γ1 is small in which χ field only controls part of the decay channels of inflaton field.It is now generally considered that a period of rapid and highly nonlinear preheating preceded the perturbative reheating phase [69]. In this case, one can instead consider modulated preheating, see, for example, [77–79]. A particularly rapid form of preheating, known as instant preheating [80], has attracted more attention in this area, as it is possible to find analytic estimates [81–85]. It is possible to generate the primordial curvature perturbation in these models by having either a weak [81] or a strong [83] symmetry breaking of the two-field potential during inflation which affects the efficiency of instant preheating, or by having two inflaton fields which are coupled with different strengths to the preheat field [82]. However, all of the above models of modulated instant preheating are only allowed, even at the linear level in perturbations, for a small corner of parameter space and even then only for a special choice of initial conditions [85].


## 2.3. Inhomogeneous End of Inflation

An inhomogeneous phase transition can occur also at the end of inflation between inflationary and non-inflationary expansion due to the condition of how inflation ends [34, 86]. In the case of multicomponent inflation, there is a family of inflationary trajectories. If the family of trajectories is not perpendicular to the surface on which inflation ends then there is a possibility that inflation ends on a slice which is not of uniform density and the generation of the curvature perturbation at the end of inflation is expected [34] (This was also discussed in [87–89] using the linear relation between two fields at the end of inflation and the non-Gaussianity is transferred from that of the light field which is generated from its nonlinear self-coupling). In single field inflation, the energy density is determined solely by the single field ϕ, and the end of inflation is parameterised by the unique value ϕe. If there is another field σ, then ϕe(σ) will depend on position through the perturbation δσ(x). As a result, the change of e-folding number from a spacetime slice of uniform energy density just before the end of inflation to a spacetime slice of uniform density just after the end of inflation will have a perturbation δNe=ζe. This ζe is the curvature perturbation generated by the end of inflation and can dominate over the contribution from inflation [34].Using the perturbation ofϕe(σ)(15)δϕe=ϕe′δσ+12ϕe′′(δσ)2,
the curvature perturbation generated at the end of inflation is [34]
(16)ζe=Ne′δϕe+12Ne′′(δϕe)2=Ne′ϕe′δσ+12[2Ne′′ϕe′2+ne′ϕe′′](δσ)2.Whenζe dominates the perturbation generated during the inflation, the power spectrum is given by
(17)𝒫e=ϕe′22ϵe(H*2π)2
The non-Gaussianity can also be generated at the end of inflation since the perturbation at the end of inflation is not required to be Gaussian. Thus, we expect that the observable large non-Gaussianity can be generated in this process. From (16) the non-Gaussianity is [34]
(18)fNL=-53ϵe2ϕe′′(ϕe′)2.
This was applied to the simple example of extended hybrid inflation using sudden end approximation [90, 91].Recently Sasaki [92], Naruko and Sasaki [93]  calculated analytically the curvature perturbation and the non-Gaussianity in the multi-brid inflation models from horizon exit to the end of inflation considering the general couplings to the waterfall field. In general, we cannot just ignore the curvature perturbation generated during inflation. Furthermore, that is also affected by the end of inflation. A deeper understanding of non-Gaussianity both from the evolution during multifield hybrid inflation and effects from the end of inflation is considered by Byrnes et al. [94, 95] which is summarised in Section 3.


## 3. Non-Gaussianity during Slow-Roll Inflation

It is well known that the single field slow-roll inflation with canonical kinetic terms generates a nonlinear parameter of the order of the spectral tilt of the spectrum and thus too small to be observed [13]. On the other hand, in multifield inflationary models the nonadiabatic perturbation may change this conclusion. Furthermore, even within slow-roll inflation an observably large non-Gaussianity can be generated if certain conditions are satisfied.Rigopoulos et al. introduced a formalism [96, 97] to deal with nonlinearity based on a system of fully nonlinear equations for long wavelengths and applied it to a two-field inflation model. In the case of two-field quadratic inflation with two massive fields, an analytic calculation is possible with slow-roll conditions, and this shows that the nonlinear parameter is small [98, 99]. The numerical calculation confirms the analytical prediction that it is possible to generate a narrow spike of large non-Gaussianity while the inflaton trajectory turns a corner, but the non-Gaussianity decays quickly after the corner [99, 100]. This spike is due to the temporary jump in the slow-roll parameters [101, 102].Recently, an analysis based on the general analytic formula was done by Byrnes et al. [94, 95], and they showed that it is possible to generate observable large non-Gaussianity even during slow-roll multifield inflation models and explicitly wrote the conditions as well as concrete examples (There is a claim [103] that the sizeable value of fNL can arise from a loop corrections). In this section, we summarise the general conditions for two-field slow-roll inflation to generate a large non-Gaussianity and briefly give a specific example. In Section 4 we consider the multifield hybrid inflation model which has been studied in depth, and we include the effects from the end of inflation.


### 3.1. General Formulas

In the case of a potential which is either sum or product separable, it is possible to compute the curvature perturbation using theδN-formalism and slow-roll conditions, and an analytic formula for the nonlinearity parameter can be obtained. The general formula for the nonlinearity parameter fNL was calculated in the case of a separable potential by sum [99], and this was later extended to the non-canonical kinetic terms as well as to the separable potential by product [104]. This has been further generalized to the arbitrary number of fields [105] and to the trispectrum [106].For a product separable potential,W(φ,χ)=U(φ)V(χ), we can find the number of e-foldings analytically in the slow-roll limit [107], (19)N(φ*,χ*)=-1MP2∫*eUUφdφ=-1MP2∫*eVVχdχ.
Then, the power spectrum and the nonlinear parameters can be calculated straightforwardly using the derivatives of e-folding number [99, 104, 107]. In the case of two inflatons with canonical kinetic terms, the power spectrum and spectral index are found to be [104] (20)𝒫ζ=W*24π2MP4(u2ϵφ*+v2ϵχ*),nζ-1=-2ϵ*-4u2(1-ηφφ*/2ϵφ*)+v2(1-ηχχ*/2ϵχ*)u2/ϵφ*+v2/ϵχ*,r=16(u2ϵφ*+v2ϵχ*)-1,
where the slow-roll parameters are (21)ϵφ=MP22(UφU)2=ϵcos2θ,ϵχ=MP22(VχV)2=ϵsin2θ,
where θ is the angle between the adiabatic perturbation and one of the fields and (22)ηφφ=MP2VφφV,ηφχ=MP2UφVχW,ηχχ=MP2VχχV.
Note that the superscript “*” and “e” denotes that the values are evaluated at horizon exit and at end of inflation, respectively. In the above, u and v are evaluated at the end of inflation as (23)u≡ϵφeϵe=cos2θe,v≡ϵχeϵe=sin2θe.The nonlinearity parameterfNL becomes [104](24)fNL=562(u2/ϵφ*+v2/ϵχ*)2×[u3ϵφ*(1-ηφφ*2ϵφ*)+v3ϵχ*(1-ηχχ*2ϵχ*)-(uϵφ*-vϵχ*)2𝒜P]
where
(25)𝒜P≡-ϵφeϵχe(ϵe)2[η̂e-4ϵφeϵχeϵe],η̂≡ϵχηφφ+ϵφηχχϵ.
Similar formulas are found in the case of a sum separable potential [99, 104].


### 3.2. Conditions for Generating a LargefNL

Using an analytic formula for the nonlinear parameterfNL in the case of a sum or product separable potential, it is possible to generate significant non-Gaussianity even during slow-roll inflation with Gaussian perturbations at Hubble exit and the general conditions for it can be written explicitly [94].For the separable potential by product large nonlinear parameter|fNL|≳1 is possible when (there is a symmetric region of large non-Gaussianity under the exchange of the two fields).(26)sin2θ*≲sin4θe(1sin2θe𝒢p-1),𝒢p=65|-ηχχ*+2ηχχe|-1.
with this condition, the large fNL is given by (27)fNL≃56sin6θe(sin2θ*+sin4θe)2[-ηχχ*+2ηχχe].
For large fNL, the substantial increase in sinθ is necessary between horizon exit and the end of inflation, which corresponds to the curve of the trajectory in the field space during the slow-roll inflation. However, the background trajectory must be almost entirely in the one-field direction.One specific toy model is the quadratic times exponential potential [94], where the potential is given by (28)W(φ,χ)=12e-λφ2/MP2m2χ2.
With expansion of exponential, this potential can be understood as the quadratic chaotic inflation model which has a coupling to another light scalar field. For this potential, considering the condition for large non-Gaussianity, (26), large fNL is generated when the φ at horizon exit is small. One example is found to give |fNL|~30 with φ*~10-3MP for λ=0.05 and χ*=16MP.For a sum separable potential, we can find similar conditions for generating a largefNL [94].


## 3.1. General Formulas

In the case of a potential which is either sum or product separable, it is possible to compute the curvature perturbation using theδN-formalism and slow-roll conditions, and an analytic formula for the nonlinearity parameter can be obtained. The general formula for the nonlinearity parameter fNL was calculated in the case of a separable potential by sum [99], and this was later extended to the non-canonical kinetic terms as well as to the separable potential by product [104]. This has been further generalized to the arbitrary number of fields [105] and to the trispectrum [106].For a product separable potential,W(φ,χ)=U(φ)V(χ), we can find the number of e-foldings analytically in the slow-roll limit [107], (19)N(φ*,χ*)=-1MP2∫*eUUφdφ=-1MP2∫*eVVχdχ.
Then, the power spectrum and the nonlinear parameters can be calculated straightforwardly using the derivatives of e-folding number [99, 104, 107]. In the case of two inflatons with canonical kinetic terms, the power spectrum and spectral index are found to be [104] (20)𝒫ζ=W*24π2MP4(u2ϵφ*+v2ϵχ*),nζ-1=-2ϵ*-4u2(1-ηφφ*/2ϵφ*)+v2(1-ηχχ*/2ϵχ*)u2/ϵφ*+v2/ϵχ*,r=16(u2ϵφ*+v2ϵχ*)-1,
where the slow-roll parameters are (21)ϵφ=MP22(UφU)2=ϵcos2θ,ϵχ=MP22(VχV)2=ϵsin2θ,
where θ is the angle between the adiabatic perturbation and one of the fields and (22)ηφφ=MP2VφφV,ηφχ=MP2UφVχW,ηχχ=MP2VχχV.
Note that the superscript “*” and “e” denotes that the values are evaluated at horizon exit and at end of inflation, respectively. In the above, u and v are evaluated at the end of inflation as (23)u≡ϵφeϵe=cos2θe,v≡ϵχeϵe=sin2θe.The nonlinearity parameterfNL becomes [104](24)fNL=562(u2/ϵφ*+v2/ϵχ*)2×[u3ϵφ*(1-ηφφ*2ϵφ*)+v3ϵχ*(1-ηχχ*2ϵχ*)-(uϵφ*-vϵχ*)2𝒜P]
where
(25)𝒜P≡-ϵφeϵχe(ϵe)2[η̂e-4ϵφeϵχeϵe],η̂≡ϵχηφφ+ϵφηχχϵ.
Similar formulas are found in the case of a sum separable potential [99, 104].


## 3.2. Conditions for Generating a LargefNL

Using an analytic formula for the nonlinear parameterfNL in the case of a sum or product separable potential, it is possible to generate significant non-Gaussianity even during slow-roll inflation with Gaussian perturbations at Hubble exit and the general conditions for it can be written explicitly [94].For the separable potential by product large nonlinear parameter|fNL|≳1 is possible when (there is a symmetric region of large non-Gaussianity under the exchange of the two fields).(26)sin2θ*≲sin4θe(1sin2θe𝒢p-1),𝒢p=65|-ηχχ*+2ηχχe|-1.
with this condition, the large fNL is given by (27)fNL≃56sin6θe(sin2θ*+sin4θe)2[-ηχχ*+2ηχχe].
For large fNL, the substantial increase in sinθ is necessary between horizon exit and the end of inflation, which corresponds to the curve of the trajectory in the field space during the slow-roll inflation. However, the background trajectory must be almost entirely in the one-field direction.One specific toy model is the quadratic times exponential potential [94], where the potential is given by (28)W(φ,χ)=12e-λφ2/MP2m2χ2.
With expansion of exponential, this potential can be understood as the quadratic chaotic inflation model which has a coupling to another light scalar field. For this potential, considering the condition for large non-Gaussianity, (26), large fNL is generated when the φ at horizon exit is small. One example is found to give |fNL|~30 with φ*~10-3MP for λ=0.05 and χ*=16MP.For a sum separable potential, we can find similar conditions for generating a largefNL [94].


## 4. Hybrid Inflation with Two Inflaton Fields

We consider a model of two field hybrid inflation, whose potential is given by(29)W(φ,χ)=W0exp(12ηφφφ2MP2)exp(12ηχχχ2MP2)≃W0(1+12ηφφφ2MP2+12ηχχχ2MP2),
which is vacuum dominated, that is, which satisfies |ηφφφ2|≪MP2 and |ηχχχ2|≪MP2. It is in this regime that the two ways of writing the potential given above are equivalent (at leading order in slow roll), so we can treat the potential as being both sum and product separable. We assume that inflation ends abruptly by a waterfall field which is heavy during inflation and hence does not affect the dynamics during inflation. First, we calculate observables during slow-roll inflation. We will consider the full potential including the waterfall field in Section 4.2, which incorporate the effects from the end of inflation considered in Section 2.3. We will see that this can lead to a change in observables on the surface where the waterfall field is destabilised.In the vacuum-dominated regime, the slow-roll solutions are(30)φ(N)=φ*e-ηφφN,χ(N)=χ*e-ηχχN,
where “*” denotes the value at the horizon exit. Throughout this section, whenever we write a quantity without making it explicit at which time it should be evaluated, we mean the equation to be valid at any time N  e-foldings after Hubble exit and while slow roll is valid. Generally, we will be interested in quantities at the end of inflation, in which case we take N=60.The slow-roll parameters are(31)ϵφ=12ηφφ2φ2MP2,ϵχ=12ηχχ2χ2MP2,ϵ=ϵφ+ϵχ.
We note that the dominant slow-roll parameters ηφφ and ηχχ are constants during inflation in the vacuum dominated regime and that they are much larger than the slow-roll parameters ϵφ and ϵχ throughout inflation.From the previous section, large non-Gaussianity can be realised in either of two regions(32)cos2θ≡φ̇2φ̇2+χ̇2≃ϵφϵφ+ϵχ≪1,orsin2θ≡χ̇2φ̇2+χ̇2≃ϵχϵφ+ϵχ≪1.
Since the two regions are symmetrical [94] (before specifying the values of ηφφ and ηχχ), in the rest of this section, we will focus on the second region. In this region where ϵφ≫ϵχ, |fNL|>1 is fulfilled by the condition (33)sin2θ*≲sin4θ(5|ηχχ|6sin2θ-1),
in other words, (34)|ηχχ|-1e-4(ηφφ-ηχχ)N≲sin2θ≃ϵχϵφ≲|ηχχ|.
This implies three inequalities on θ(35)sin2θ*<13(56)2(34)4|ηχχ|2,sin2θ<56|ηχχ|,sin2θsin2θ*>2451|ηχχ|.
Note that in this region sinθ≃ηχχχ/(ηφφφ) from (30), we require N(ηφφ-ηχχ)>1 so that sin2θ grows significantly during inflation.


### 4.1. Simplified Formula for the Observables WhenfNL Is Large

We can substantially simplify all of the above formula in the case wherefNL is large. We define the quantity (36)r̃≡(∂N/∂χ*)2(∂N/∂φ*)2=ϵχϵφe2(ηφφ-ηχχ)N.
In the region, we are considering where fNL is large, this is approximately given by the initial and final angles of the background trajectory with different exponents (37)r̃≃sin4θsin2θ*.In the case of large non-Gaussianity it follows that(38)𝒫ζ≃W*24π2MP4ϵ*(1+ϵχϵφe2(ηφφ-ηχχ)N)=8r(H*2π)2,(39)nζ-1≃2ηφφ+r̃ηχχ1+r̃,(40)r≃16ϵ*1+r̃,(41)fNL≃56sin6θe(sin2θ*+sin4θe)2ηχχ=56r̃(1+r̃)2ηχχe2(ηφφ-ηχχ)N.The first condition in (35) implies that (42)χ*≪φ*.
We, therefore, require a very small value of χ* in order to have a large non-Gaussianity. While this may be a fine tuning, we comment that requiring a large ratio of the initial field values is common in models where the isocurvature perturbation generates a significant part of the primordial adiabatic perturbation, see the discussion in the conclusion of [85]. If the inflaton perturbations are neglected a priori, such as in the pure (i.e., not mixed) curvaton scenario, this remains implicitly true. For an explicit discussion of this, in the simple case that both the inflaton and curvaton fields have a quadratic potential, see [108].The sign offNL is determined by the sign of ηχχ. The amplitude of fNL depends exponentially on the difference of the slow-roll parameters, ηφφ-ηχχ, which we require to be positive to be in the branch of large non-Gaussianity where sin2θ≪1, while the spectral index depends on a weighted sum of the slow-roll parameters, so it is possible to have a large non-Gaussianity and a scale invariant spectrum. However, it is not possible to have a large and positive fNL and a red spectrum of perturbations. We will see in Section 4.2.2 that by including the effect of the waterfall field this conclusion may change, depending on the values of the coupling constants between the two inflaton fields and the waterfall field.In Table1, we give some explicit examples of values of ηφφ,ηχχ,φ* and χ* which lead to a large non-Gaussianity. Using (4), we also calculate the spectral index. The first example in Table 1 shows that it is possible to have |fNL|≃100 and a scale invariant spectrum. We also see that it is possible to generate a large non-Gaussianity during slow roll with ηφφ and ηχχ both positive or both negative, or when one is positive and the other negative corresponding to a saddle point. The results for this model were verified using a novel calculational method of momentum transport equations in [109].Table 1
Table showing some initial conditions for the hybrid inflation model that lead to large levels of non-Gaussianity. The table shows the parameterr̃, the bispectrum and the trispectrum nonlinearity parameters, the spectral index and tensor-to-scalar ratio. They are evaluated when the number of e-foldings from the end of inflation is Nk=60. The trispectrum parameters τNL and gNL are defined and calculated in Section 6.1.
ηφφηχχφ*χ*r̃fNLτNLgNLnζ-1r0.04-0.0416.8×10-51-1234.4×104-3300.0060.04-0.0411.5×10-45-688×103-24-0.050.0020.080.0110.001819.272470.770.090.0260.02-0.0410.000371-11.1357-2.6-0.020.002-0.01-0.0913×10-60.16-1321.8×105-44-0.040.00070.06-0.0114.3×10-40.1-3148-0.20.110.0260.01-0.0617.5×10-60.04-82.5×103-20.010.0008


### 4.2. Effect of the Waterfall Field and Further Evolution after Inflation

In this section, we include the effects of the waterfall fieldρ which is required to end hybrid inflation. Inflation ends when the waterfall field is destabilised that is, when its effective mass becomes negative. During inflation, the waterfall field is heavy, and it is trapped with a vacuum expectation value of zero, so we can neglect it during inflation. The end of inflation occurs when the effective mass of the waterfall field is zero, which occurs on a hypersurface defined in general by [92, 93]
(43)σ2=G(φ,χ)≡g12φ2+g22χ2,
which is realised by the potential W(φ,χ), defined by (29), where W0 is given by
(44)W0=12G(φ,χ)ρ2+λ4(ρ2-σ2λ)2.
Here, g1 (g2) is the coupling between the φ (χ) field and ρ is the waterfall field with self-coupling λ. In general, the hypersurface defined by this end condition is not a surface of uniform energy density.As discussed earlier, this is an example of a model with an inhomogeneous end of inflation, that is, where inflation ends at slightly different times in different places. It has also been shown for the hybrid potential we are considering that this can be used to generate a large amount of non-Gaussianity, for certain parameters values and fine tuning of the parameters [90, 93]. However, these papers concern the large non-Gaussianity generated at the end of inflation rather than during slow-roll inflation, by having a very large ratio of couplings g1/g2≪1. Here, we consider the case where g1 and g2 have the same order of magnitude with g12/g22=ηφφ/ηχχ in Section 4.2.1 and with g12=g22 in Section 4.2.2.


#### 4.2.1.g12/g22=ηφφ/ηχχ

In this case, we have chosen the coupling constants (which can satisfyg12<0 and/or g22<0) such that the surface where the waterfall field is destabilised corresponds to a surface of uniform energy density. This is because the end condition can be rewritten as σ2=2g22/ηχχ((1/2)ηφφφ2+(1/2)ηχχχ2)∝W. In this case, the value of all observable quantities such as the power spectrum and non-Gaussianity are the same as those we calculated previously which were valid at the final hypersurface of uniform energy density during inflation.


#### 4.2.2.g12=g22

In this case, the end of inflation given by the condition in (43) does not occur on a uniform energy density hypersurface [93]. We will show how the non-Gaussianity is modified by the condition at the end of inflation in this example. In general, we expect there to be some modification to non-Gaussianity from the end of inflation, except in the special case we considered in Section 4.2.1. We have checked in [95] that the correction from the extra expansion which occurs from the surface on which inflation ends up to a surface of uniform energy density is small, and therefore one can still use the δN formalism in this case.In this case, the power spectrum andfNL in the limit of large non-Gaussianity are [93] (45)𝒫ζ≃W*24π2MP4ϵ*(1+ηφφ2ηχχ2ϵχϵφe2(ηφφ-ηχχ)N)=8r(H*2π)2,(46)nζ-1≃2ηφφ+(ηφφ2/ηχχ2)r̃ηχχ1+(ηφφ2/ηχχ2)r̃,(47)r≃16ϵ*(1+ηφφ2ηχχ2r̃)-1,(48)fNL≃56(ηφφ2/ηχχ2)r̃(1+(ηφφ2/ηχχ2)r̃)2ηφφe2(ηφφ-ηχχ)N.
We can see that the observables at the end of inflation are changed by the additional ratio of ηφφ2/ηχχ2 before r̃. Furthermore, fNL has an additional factor of ηφφ/ηχχ.In Table2, we show the values of fNL, nζ-1 and r with end condition g12=g22 for the same parameter values that we used in Table 1, where the inflation ends on the uniform energy density hypersurface. The first two examples show that if ηφφ=-ηχχ, then the observables are unchanged for two different end conditions except that the sign of fNL is switched. The second example in the table shows that in this case it is possible to have a red spectral index and a positive value of fNL. For many values of the initial parameters when |ηφφ/ηχχ|≠1, the magnitude of fNL decreases compared to Table 1, but the reverse can also happen, an example of this where |fNL| grows by more than an order of magnitude is shown in the final row of Tables 1 and 2.Table 2
Same as Table1, but with different end condition, g12=g22 as used in Section 4.2.2.
ηφφηχχφ*χ*r̃fNLnζ-1r0.04-0.0416.8×10-5112300.0060.04-0.0411.5×10-4568-0.050.0020.080.0110.001814.590.020.00080.02-0.0410.0003713.50.020.026-0.01-0.0913×10-60.16-0.2-0.020.00080.06-0.0114.3×10-40.1380.010.006


#### 4.2.3. Further Evolution after Inflation

So far, we have assumed a quick transition to the radiation epoch at the end of inflation, thereby neglecting the dynamics of the waterfall field. However, if we consider the role of the waterfall field, then after the waterfall field is destabilised, there may be a further evolution of the primordial curvature perturbation, which will lead to a change of the observable parameters. This applies to any model with an inhomogeneous end of inflation, since there are isocurvature perturbations still present after the waterfall field is destabilised and inflation has ended. Further evolution will depend on the details of reheating in a model dependent way. To the best of our knowledge, this issue has not been considered in depth in any paper. If we assume an instantaneous transition to radiation domination (so a completely efficient and immediate decay of the waterfall and inflaton fields), then there will be no further change to the observables as we have argued in the previous section. However, this is clearly an idealised case.In the special case where the waterfall field is also light during inflation, Barnaby and Cline [110] have shown there is the possibility of generating a large non-Gaussianity during preheating for certain parameter values. This is possible even if there is only one inflaton field and the waterfall field present. However, in this case, inflation does not end abruptly when the waterfall field is destabilised so this is not the scenario we have considered in this paper. A recent claim that even a heavy waterfall field may lead to a scale invariant local non-Gaussianity has been withdrawn [111].


## 4.1. Simplified Formula for the Observables WhenfNL Is Large

We can substantially simplify all of the above formula in the case wherefNL is large. We define the quantity (36)r̃≡(∂N/∂χ*)2(∂N/∂φ*)2=ϵχϵφe2(ηφφ-ηχχ)N.
In the region, we are considering where fNL is large, this is approximately given by the initial and final angles of the background trajectory with different exponents (37)r̃≃sin4θsin2θ*.In the case of large non-Gaussianity it follows that(38)𝒫ζ≃W*24π2MP4ϵ*(1+ϵχϵφe2(ηφφ-ηχχ)N)=8r(H*2π)2,(39)nζ-1≃2ηφφ+r̃ηχχ1+r̃,(40)r≃16ϵ*1+r̃,(41)fNL≃56sin6θe(sin2θ*+sin4θe)2ηχχ=56r̃(1+r̃)2ηχχe2(ηφφ-ηχχ)N.The first condition in (35) implies that (42)χ*≪φ*.
We, therefore, require a very small value of χ* in order to have a large non-Gaussianity. While this may be a fine tuning, we comment that requiring a large ratio of the initial field values is common in models where the isocurvature perturbation generates a significant part of the primordial adiabatic perturbation, see the discussion in the conclusion of [85]. If the inflaton perturbations are neglected a priori, such as in the pure (i.e., not mixed) curvaton scenario, this remains implicitly true. For an explicit discussion of this, in the simple case that both the inflaton and curvaton fields have a quadratic potential, see [108].The sign offNL is determined by the sign of ηχχ. The amplitude of fNL depends exponentially on the difference of the slow-roll parameters, ηφφ-ηχχ, which we require to be positive to be in the branch of large non-Gaussianity where sin2θ≪1, while the spectral index depends on a weighted sum of the slow-roll parameters, so it is possible to have a large non-Gaussianity and a scale invariant spectrum. However, it is not possible to have a large and positive fNL and a red spectrum of perturbations. We will see in Section 4.2.2 that by including the effect of the waterfall field this conclusion may change, depending on the values of the coupling constants between the two inflaton fields and the waterfall field.In Table1, we give some explicit examples of values of ηφφ,ηχχ,φ* and χ* which lead to a large non-Gaussianity. Using (4), we also calculate the spectral index. The first example in Table 1 shows that it is possible to have |fNL|≃100 and a scale invariant spectrum. We also see that it is possible to generate a large non-Gaussianity during slow roll with ηφφ and ηχχ both positive or both negative, or when one is positive and the other negative corresponding to a saddle point. The results for this model were verified using a novel calculational method of momentum transport equations in [109].Table 1
Table showing some initial conditions for the hybrid inflation model that lead to large levels of non-Gaussianity. The table shows the parameterr̃, the bispectrum and the trispectrum nonlinearity parameters, the spectral index and tensor-to-scalar ratio. They are evaluated when the number of e-foldings from the end of inflation is Nk=60. The trispectrum parameters τNL and gNL are defined and calculated in Section 6.1.
ηφφηχχφ*χ*r̃fNLτNLgNLnζ-1r0.04-0.0416.8×10-51-1234.4×104-3300.0060.04-0.0411.5×10-45-688×103-24-0.050.0020.080.0110.001819.272470.770.090.0260.02-0.0410.000371-11.1357-2.6-0.020.002-0.01-0.0913×10-60.16-1321.8×105-44-0.040.00070.06-0.0114.3×10-40.1-3148-0.20.110.0260.01-0.0617.5×10-60.04-82.5×103-20.010.0008


## 4.2. Effect of the Waterfall Field and Further Evolution after Inflation

In this section, we include the effects of the waterfall fieldρ which is required to end hybrid inflation. Inflation ends when the waterfall field is destabilised that is, when its effective mass becomes negative. During inflation, the waterfall field is heavy, and it is trapped with a vacuum expectation value of zero, so we can neglect it during inflation. The end of inflation occurs when the effective mass of the waterfall field is zero, which occurs on a hypersurface defined in general by [92, 93]
(43)σ2=G(φ,χ)≡g12φ2+g22χ2,
which is realised by the potential W(φ,χ), defined by (29), where W0 is given by
(44)W0=12G(φ,χ)ρ2+λ4(ρ2-σ2λ)2.
Here, g1 (g2) is the coupling between the φ (χ) field and ρ is the waterfall field with self-coupling λ. In general, the hypersurface defined by this end condition is not a surface of uniform energy density.As discussed earlier, this is an example of a model with an inhomogeneous end of inflation, that is, where inflation ends at slightly different times in different places. It has also been shown for the hybrid potential we are considering that this can be used to generate a large amount of non-Gaussianity, for certain parameters values and fine tuning of the parameters [90, 93]. However, these papers concern the large non-Gaussianity generated at the end of inflation rather than during slow-roll inflation, by having a very large ratio of couplings g1/g2≪1. Here, we consider the case where g1 and g2 have the same order of magnitude with g12/g22=ηφφ/ηχχ in Section 4.2.1 and with g12=g22 in Section 4.2.2.


### 4.2.1.g12/g22=ηφφ/ηχχ

In this case, we have chosen the coupling constants (which can satisfyg12<0 and/or g22<0) such that the surface where the waterfall field is destabilised corresponds to a surface of uniform energy density. This is because the end condition can be rewritten as σ2=2g22/ηχχ((1/2)ηφφφ2+(1/2)ηχχχ2)∝W. In this case, the value of all observable quantities such as the power spectrum and non-Gaussianity are the same as those we calculated previously which were valid at the final hypersurface of uniform energy density during inflation.


### 4.2.2.g12=g22

In this case, the end of inflation given by the condition in (43) does not occur on a uniform energy density hypersurface [93]. We will show how the non-Gaussianity is modified by the condition at the end of inflation in this example. In general, we expect there to be some modification to non-Gaussianity from the end of inflation, except in the special case we considered in Section 4.2.1. We have checked in [95] that the correction from the extra expansion which occurs from the surface on which inflation ends up to a surface of uniform energy density is small, and therefore one can still use the δN formalism in this case.In this case, the power spectrum andfNL in the limit of large non-Gaussianity are [93] (45)𝒫ζ≃W*24π2MP4ϵ*(1+ηφφ2ηχχ2ϵχϵφe2(ηφφ-ηχχ)N)=8r(H*2π)2,(46)nζ-1≃2ηφφ+(ηφφ2/ηχχ2)r̃ηχχ1+(ηφφ2/ηχχ2)r̃,(47)r≃16ϵ*(1+ηφφ2ηχχ2r̃)-1,(48)fNL≃56(ηφφ2/ηχχ2)r̃(1+(ηφφ2/ηχχ2)r̃)2ηφφe2(ηφφ-ηχχ)N.
We can see that the observables at the end of inflation are changed by the additional ratio of ηφφ2/ηχχ2 before r̃. Furthermore, fNL has an additional factor of ηφφ/ηχχ.In Table2, we show the values of fNL, nζ-1 and r with end condition g12=g22 for the same parameter values that we used in Table 1, where the inflation ends on the uniform energy density hypersurface. The first two examples show that if ηφφ=-ηχχ, then the observables are unchanged for two different end conditions except that the sign of fNL is switched. The second example in the table shows that in this case it is possible to have a red spectral index and a positive value of fNL. For many values of the initial parameters when |ηφφ/ηχχ|≠1, the magnitude of fNL decreases compared to Table 1, but the reverse can also happen, an example of this where |fNL| grows by more than an order of magnitude is shown in the final row of Tables 1 and 2.Table 2
Same as Table1, but with different end condition, g12=g22 as used in Section 4.2.2.
ηφφηχχφ*χ*r̃fNLnζ-1r0.04-0.0416.8×10-5112300.0060.04-0.0411.5×10-4568-0.050.0020.080.0110.001814.590.020.00080.02-0.0410.0003713.50.020.026-0.01-0.0913×10-60.16-0.2-0.020.00080.06-0.0114.3×10-40.1380.010.006


### 4.2.3. Further Evolution after Inflation

So far, we have assumed a quick transition to the radiation epoch at the end of inflation, thereby neglecting the dynamics of the waterfall field. However, if we consider the role of the waterfall field, then after the waterfall field is destabilised, there may be a further evolution of the primordial curvature perturbation, which will lead to a change of the observable parameters. This applies to any model with an inhomogeneous end of inflation, since there are isocurvature perturbations still present after the waterfall field is destabilised and inflation has ended. Further evolution will depend on the details of reheating in a model dependent way. To the best of our knowledge, this issue has not been considered in depth in any paper. If we assume an instantaneous transition to radiation domination (so a completely efficient and immediate decay of the waterfall and inflaton fields), then there will be no further change to the observables as we have argued in the previous section. However, this is clearly an idealised case.In the special case where the waterfall field is also light during inflation, Barnaby and Cline [110] have shown there is the possibility of generating a large non-Gaussianity during preheating for certain parameter values. This is possible even if there is only one inflaton field and the waterfall field present. However, in this case, inflation does not end abruptly when the waterfall field is destabilised so this is not the scenario we have considered in this paper. A recent claim that even a heavy waterfall field may lead to a scale invariant local non-Gaussianity has been withdrawn [111].


## 4.2.1.g12/g22=ηφφ/ηχχ

In this case, we have chosen the coupling constants (which can satisfyg12<0 and/or g22<0) such that the surface where the waterfall field is destabilised corresponds to a surface of uniform energy density. This is because the end condition can be rewritten as σ2=2g22/ηχχ((1/2)ηφφφ2+(1/2)ηχχχ2)∝W. In this case, the value of all observable quantities such as the power spectrum and non-Gaussianity are the same as those we calculated previously which were valid at the final hypersurface of uniform energy density during inflation.


## 4.2.2.g12=g22

In this case, the end of inflation given by the condition in (43) does not occur on a uniform energy density hypersurface [93]. We will show how the non-Gaussianity is modified by the condition at the end of inflation in this example. In general, we expect there to be some modification to non-Gaussianity from the end of inflation, except in the special case we considered in Section 4.2.1. We have checked in [95] that the correction from the extra expansion which occurs from the surface on which inflation ends up to a surface of uniform energy density is small, and therefore one can still use the δN formalism in this case.In this case, the power spectrum andfNL in the limit of large non-Gaussianity are [93] (45)𝒫ζ≃W*24π2MP4ϵ*(1+ηφφ2ηχχ2ϵχϵφe2(ηφφ-ηχχ)N)=8r(H*2π)2,(46)nζ-1≃2ηφφ+(ηφφ2/ηχχ2)r̃ηχχ1+(ηφφ2/ηχχ2)r̃,(47)r≃16ϵ*(1+ηφφ2ηχχ2r̃)-1,(48)fNL≃56(ηφφ2/ηχχ2)r̃(1+(ηφφ2/ηχχ2)r̃)2ηφφe2(ηφφ-ηχχ)N.
We can see that the observables at the end of inflation are changed by the additional ratio of ηφφ2/ηχχ2 before r̃. Furthermore, fNL has an additional factor of ηφφ/ηχχ.In Table2, we show the values of fNL, nζ-1 and r with end condition g12=g22 for the same parameter values that we used in Table 1, where the inflation ends on the uniform energy density hypersurface. The first two examples show that if ηφφ=-ηχχ, then the observables are unchanged for two different end conditions except that the sign of fNL is switched. The second example in the table shows that in this case it is possible to have a red spectral index and a positive value of fNL. For many values of the initial parameters when |ηφφ/ηχχ|≠1, the magnitude of fNL decreases compared to Table 1, but the reverse can also happen, an example of this where |fNL| grows by more than an order of magnitude is shown in the final row of Tables 1 and 2.Table 2
Same as Table1, but with different end condition, g12=g22 as used in Section 4.2.2.
ηφφηχχφ*χ*r̃fNLnζ-1r0.04-0.0416.8×10-5112300.0060.04-0.0411.5×10-4568-0.050.0020.080.0110.001814.590.020.00080.02-0.0410.0003713.50.020.026-0.01-0.0913×10-60.16-0.2-0.020.00080.06-0.0114.3×10-40.1380.010.006


## 4.2.3. Further Evolution after Inflation

So far, we have assumed a quick transition to the radiation epoch at the end of inflation, thereby neglecting the dynamics of the waterfall field. However, if we consider the role of the waterfall field, then after the waterfall field is destabilised, there may be a further evolution of the primordial curvature perturbation, which will lead to a change of the observable parameters. This applies to any model with an inhomogeneous end of inflation, since there are isocurvature perturbations still present after the waterfall field is destabilised and inflation has ended. Further evolution will depend on the details of reheating in a model dependent way. To the best of our knowledge, this issue has not been considered in depth in any paper. If we assume an instantaneous transition to radiation domination (so a completely efficient and immediate decay of the waterfall and inflaton fields), then there will be no further change to the observables as we have argued in the previous section. However, this is clearly an idealised case.In the special case where the waterfall field is also light during inflation, Barnaby and Cline [110] have shown there is the possibility of generating a large non-Gaussianity during preheating for certain parameter values. This is possible even if there is only one inflaton field and the waterfall field present. However, in this case, inflation does not end abruptly when the waterfall field is destabilised so this is not the scenario we have considered in this paper. A recent claim that even a heavy waterfall field may lead to a scale invariant local non-Gaussianity has been withdrawn [111].


## 5. Multiple-Field Inflation without Slow Roll

In this section, we demonstrate a new method that provides, for certain classes of models, analytical expressions forfNL valid in regimes beyond a slow-roll approximation [112]. This method is based on the first-order Hamilton-Jacobi formalism developed by Salopek and Bond [113], which allows us to express inflationary observables in the multifield case, without having to focus on a slow-roll regime (see [114] for a similar application of this formalism to the single-field case).The main idea on which our arguments are based is to focus on inflationary multifield trajectories in which theHubble rate, and not the potential, is separable. This approach allows us to analytically study cases where the non-Gaussian parameter fNL becomes large, in regimes in which the slow-roll parameters are enhanced after Hubble exit, but before inflation ends. In particular, we provide analytical equations that express the nonlinearity parameter fNL, in terms of quantities that generalise the slow-roll parameters during inflation. Using this method, we demonstrate in the next subsection an exact solution of multifield inflation which can give rise to a large non-Gaussianity due to a strong break down in slow roll shortly before inflation ends.The first-order Hamilton-Jacobi equations of motion are(49)H2=13W(ϕ,χ)+23[(∂H∂ϕ)2+(∂H∂χ)2],ϕ̇=-2∂H∂ϕ,χ̇=-2∂H∂χ.We concentrate on solutions satisfying the following Ansatz for the Hubble parameterH(ϕ,χ)(50)H(ϕ,χ)=H(1)(ϕ)+H(2)(χ),
that is, we demand that it can be split as a sum (hence, we follow an analysis related to [99]. In an analogous way, one could also consider a situation in which the Hubble parameter splits as a product of single-field pieces [104]. The formula's in this section were extended to an arbitrary number of fields with a sum separable Hubble parameter in [115]) of two pieces, each one depending on a single field. Notice that this is different with respect to the requirement of separable potentials usually done in the literature and earlier in this paper.It is convenient to define the following quantities:(51)δϕ=(H,ϕ(1)H)2,δχ=(H,χ(2)H)2,γϕ=H,ϕϕ(1)H,γχ=H,χχ(2)H,
and δ=δϕ+δχ. Although their definition resembles the corresponding one for the usual slow-roll parameters ϵ and η, they do not coincide with them when taking a slow roll limit (it is nevertheless simple to work out the relation with the slow roll parameters, defined by (21) and (22). In a slow-roll regime, one finds δϕ≃ϵϕ/2,γϕ≃(ηϕϕ-ϵϕ)/2). We assume that the quantities δ and γ are much smaller than unity at Hubble exit t=t*, in order that we can use the δN formalism. During inflation, by definition, we have to ensure that the quantity ϵH≡-Ḣ/H2=2δ<1. The quantities γ can, however, become much larger than unity during inflation.Using similar techniques to those applied to a separable potential, one can calculate the spectral index andfNL, the details were given in [112] and the results are (52)nζ-1=-4δ*-4uH(1-(γ*ϕ/δ*ϕ)uH)uH2/δ*ϕ+vH2/δ*χ+(-4)vH(1-(γ*χ/δ*χ)vH)uH2/δ*ϕ+vH2/δ*χ(53)65fNL=2(uH2/δ*ϕ)(1-(γ*ϕ/δ*ϕ)uH)(uH2/δ*ϕ+vH2/δ*χ)2+2(v2/δ*χ)(1-(γ*χ/δ*χ)vH)(uH2/δ*ϕ+vH2/δ*χ)2+22(uH/δ*ϕ-vH/δ*χ)2𝒜H(uH2/δ*ϕ+vH2/δ*χ)2,
where the new symbols are defined by (54)uH≡H*(1)+ZeH*,vH≡H*(2)-ZeH*,Ze=(He(2)δeϕ-He(1)δeχ)δe,𝒜H=-He2H*2δeϕδeχδe(12-γessδe),γss=(δχγϕ+δϕγχ)δ.We reiterate that the formula forfNL is exact (apart from neglecting a slow-roll suppressed contribution due to the non-Gaussianity of the fields at Hubble exit) and is not based on a slow-roll expansion. Although δ* and γ* are necessarily small, the quantity 𝒜 is not suppressed by these parameters and can assume large values during inflation, enhancing fNL. We note in agreement with [99] that if one of the fields has reached a minimum so ϕ̇=0 or χ̇=0 at the end of inflation then 𝒜H=0 and the non-Gaussianity at the end of inflation will be very small. Therefore, any model of inflation with a separable potential or separable Hubble factor with a large non-Gaussianity present at the end of inflation must have both fields still evolving, and therefore the presence of isocurvature modes. It would, therefore, also be interesting to study the evolution of the perturbations after inflation, to see if this leads to an important change in the observables [116].


### 5.1. Exact Solution

Very few exact solutions in multiple-field inflation are known [117]. Here, we present an exact solution which can give rise to a very large non-Gaussianity at the end of inflation. A potential with a similar form may be motivated in string theory, in the context of Kähler moduli inflation [118–125].The potential we consider is [112] (55)W(ϕ,χ)=U0(A1B12-A1e-αϕ+A2e-2αϕ-B1e-βχ+B2e-2βχ+A1B12e-αϕ-βχ),
where the parameters α and β satisfy (56)α2=32-6A2A12,β2=32-6B2B12.The Hamilton-Jacobi equations have the following exact solution forH and the fields:(57)H=H0(1-A12e-αϕ-B12e-βχ),ϕ=1αln[eαϕ*-A1α2H0t],χ=1βln[eβχ*-B1β2H0t],
where we have defined H0=U0/3, t is cosmic time, and we have set t*=0. The overall factor to the potential U0 can be freely chosen so that the amplitude of the scalar power spectrum matches the observed amplitude of perturbations in the CMB. The scale factor results (58)a(t)=a0(eαϕ*-α2A1H0t)1/2α2(eβχ*-β2B1H0t)1/2β2eH0t.
Notice that the solution becomes singular at late times, when the scale factor vanishes and the field values diverge. This singularity occurs well after inflation ends, and we will not need to discuss it in our analysis. From the relations (59)γϕ=-αδϕ,γχ=-βδχ,
we notice that in this example γϕ and γχ can become much larger than unity at the end of inflation if |α| and/or |β| are much greater than one. This corresponds to a break down in slow roll, although ϵH remains smaller than unity during the inflationary era by definition.As a simple, concrete example of parameter choice which gives rise to a large non-Gaussianity, we choose a regime in which the parameters|α| and |β| are both large, let us say larger than some quantity R≫1. We note from (56) that we are, therefore, required to have A2<0 and B2<0. In this case, the potential does not have a minimum, but inflation still ends through ϵH growing larger than unity, and we can trust our results in this regime. The potential would need modifications which apply after inflation in order for reheating to take place and these modifications may provide a minimum for the potential. We parameterise the values of the quantities δeϕ and δeχ at the end of inflation as
(60)δeϕ=12m,δeχ=12nsuchthat1m+1n=1,
where the final equality follows from ϵH=1. We assume that the quantities m and n are not too big. Namely, they satisfy the inequality m,n≪R2. This implies that H0≃H*≃He. We choose |γ*ϕ|=1/40p, |γ*χ|=1/40q, such that we can write (61)Ntot=60≃10(p+q).
This fixes the initial values for the fields ϕ* and χ*, to the values (62)eαϕ*=20pα2A1,eβχ*=20qβ2B1.The conditions (60) imply (63)H0te=20p-m2α2,H0te=20q-n2β2.
Since m and n are much smaller than R2, we expect p≃q. We hence have the simple relations (64)δ*ϕ≃14Ntot2α2,δ*χ≃14Ntot2β2,γ*ϕ≃γ*χ≃-12Ntot.
From (52), we find (65)nζ-1≃-2Ntot≃-0.04.
So, we have a red spectral index in agreement with present day observations [126] and in [112] it was also shown that the tensor-to-scalar ratio for this model is negligible.When evaluated at the end of inflation, we find(66)65fNL≃-22mn[α2n-β2mα2n2+β2m2]2(αn+βm).
We stress that this formula provide only the dominant contributions to fNL. It is valid in the case |α|,|β|>R≫1, and m,n<R2.As an explicit case, we takeα=-100,β=20,m=6,n=6/5 and the initial conditions to satisfy (62) with Ntot=60. Then, independently of the values of A1 and B1, we find fNL≃58 from the simplified formula (66), and fNL≃53 from the full formula (53). See Figures 1 and 2. Notice that the results for fNL differ from each other by quantities of order 1/R=1/β as expected by the approximations we have made. At the end of inflation, the parameters are γϕ=29, γχ=-13 and this shows that the slow-roll approximation has been significantly broken, (it is nevertheless simple to work out the relation with the slow roll parameters, defined by (21) and (22). In a slow-roll regime, one finds δϕ≃ϵϕ/2,γϕ≃(ηϕϕ-ϵϕ)/2). From Figure 2, one can see that the trajectory is straight for much of inflation but it turn near the end of inflation (during the last e-folding before inflation ends), and it is during this time that fNL grows larger than unity.Figure 1
Plot showingfNL as a function of ϵH towards the end of inflation, for the values of the parameters given in the text. Inflation ends when ϵH=1; for this example fNL≃59 at that time.(a) shows the trajectory considered for the parameters given after (66) superimposed on a contour plot of the potential. The trajectory is shown from the time of horizon crossing at t* and runs until the end of inflation at N=60. The square on the trajectory, labelled with N=59, indicates a point along the trajectory one e-folding before inflation ends as ϕ and χ roll towards zero. This shows that the fields roll much more quickly during the final stage of inflation, and the trajectory curves near the end. (b) shows the potential for the same parameter values. Notice that inflation ends on the plateau long before the potential becomes negative.
(a)(b)Notice that our requirements of final values for the quantitiesδeϕ and δeχ (both much bigger than 1/R) imposes fine-tuning constraints on the initial conditions, since at leading order in 1/R the values for γ*ϕ and γ*χ must coincide (see (62) and recall that p≃q).We stress that the enhancement of non-Gaussianity occurs towards the end of inflation. What happens justafter inflation is a model-dependent issue; we cannot address this question within the approximations used in this concrete model. For our choice of parameters, we notice that fNL is still increasing at the end of inflation, see Figure 1, but for other choices it may start to decrease before inflation ends, see [115]. We stress that in the regime where |α|,|β|≫1 the slow-roll parameters will necessarily become much greater than unity by the end of inflation, which may correspond to |γss|≫1: this is one reason by which fNL can become large. Our formalism is at least in principle suitable to study the evolution of non-Gaussianity after inflation, when the parameter ϵH becomes larger than unity (although we restate that the specific potential being considered here anyway needs to be modified after inflation ends).


## 5.1. Exact Solution

Very few exact solutions in multiple-field inflation are known [117]. Here, we present an exact solution which can give rise to a very large non-Gaussianity at the end of inflation. A potential with a similar form may be motivated in string theory, in the context of Kähler moduli inflation [118–125].The potential we consider is [112] (55)W(ϕ,χ)=U0(A1B12-A1e-αϕ+A2e-2αϕ-B1e-βχ+B2e-2βχ+A1B12e-αϕ-βχ),
where the parameters α and β satisfy (56)α2=32-6A2A12,β2=32-6B2B12.The Hamilton-Jacobi equations have the following exact solution forH and the fields:(57)H=H0(1-A12e-αϕ-B12e-βχ),ϕ=1αln[eαϕ*-A1α2H0t],χ=1βln[eβχ*-B1β2H0t],
where we have defined H0=U0/3, t is cosmic time, and we have set t*=0. The overall factor to the potential U0 can be freely chosen so that the amplitude of the scalar power spectrum matches the observed amplitude of perturbations in the CMB. The scale factor results (58)a(t)=a0(eαϕ*-α2A1H0t)1/2α2(eβχ*-β2B1H0t)1/2β2eH0t.
Notice that the solution becomes singular at late times, when the scale factor vanishes and the field values diverge. This singularity occurs well after inflation ends, and we will not need to discuss it in our analysis. From the relations (59)γϕ=-αδϕ,γχ=-βδχ,
we notice that in this example γϕ and γχ can become much larger than unity at the end of inflation if |α| and/or |β| are much greater than one. This corresponds to a break down in slow roll, although ϵH remains smaller than unity during the inflationary era by definition.As a simple, concrete example of parameter choice which gives rise to a large non-Gaussianity, we choose a regime in which the parameters|α| and |β| are both large, let us say larger than some quantity R≫1. We note from (56) that we are, therefore, required to have A2<0 and B2<0. In this case, the potential does not have a minimum, but inflation still ends through ϵH growing larger than unity, and we can trust our results in this regime. The potential would need modifications which apply after inflation in order for reheating to take place and these modifications may provide a minimum for the potential. We parameterise the values of the quantities δeϕ and δeχ at the end of inflation as
(60)δeϕ=12m,δeχ=12nsuchthat1m+1n=1,
where the final equality follows from ϵH=1. We assume that the quantities m and n are not too big. Namely, they satisfy the inequality m,n≪R2. This implies that H0≃H*≃He. We choose |γ*ϕ|=1/40p, |γ*χ|=1/40q, such that we can write (61)Ntot=60≃10(p+q).
This fixes the initial values for the fields ϕ* and χ*, to the values (62)eαϕ*=20pα2A1,eβχ*=20qβ2B1.The conditions (60) imply (63)H0te=20p-m2α2,H0te=20q-n2β2.
Since m and n are much smaller than R2, we expect p≃q. We hence have the simple relations (64)δ*ϕ≃14Ntot2α2,δ*χ≃14Ntot2β2,γ*ϕ≃γ*χ≃-12Ntot.
From (52), we find (65)nζ-1≃-2Ntot≃-0.04.
So, we have a red spectral index in agreement with present day observations [126] and in [112] it was also shown that the tensor-to-scalar ratio for this model is negligible.When evaluated at the end of inflation, we find(66)65fNL≃-22mn[α2n-β2mα2n2+β2m2]2(αn+βm).
We stress that this formula provide only the dominant contributions to fNL. It is valid in the case |α|,|β|>R≫1, and m,n<R2.As an explicit case, we takeα=-100,β=20,m=6,n=6/5 and the initial conditions to satisfy (62) with Ntot=60. Then, independently of the values of A1 and B1, we find fNL≃58 from the simplified formula (66), and fNL≃53 from the full formula (53). See Figures 1 and 2. Notice that the results for fNL differ from each other by quantities of order 1/R=1/β as expected by the approximations we have made. At the end of inflation, the parameters are γϕ=29, γχ=-13 and this shows that the slow-roll approximation has been significantly broken, (it is nevertheless simple to work out the relation with the slow roll parameters, defined by (21) and (22). In a slow-roll regime, one finds δϕ≃ϵϕ/2,γϕ≃(ηϕϕ-ϵϕ)/2). From Figure 2, one can see that the trajectory is straight for much of inflation but it turn near the end of inflation (during the last e-folding before inflation ends), and it is during this time that fNL grows larger than unity.Figure 1
Plot showingfNL as a function of ϵH towards the end of inflation, for the values of the parameters given in the text. Inflation ends when ϵH=1; for this example fNL≃59 at that time.(a) shows the trajectory considered for the parameters given after (66) superimposed on a contour plot of the potential. The trajectory is shown from the time of horizon crossing at t* and runs until the end of inflation at N=60. The square on the trajectory, labelled with N=59, indicates a point along the trajectory one e-folding before inflation ends as ϕ and χ roll towards zero. This shows that the fields roll much more quickly during the final stage of inflation, and the trajectory curves near the end. (b) shows the potential for the same parameter values. Notice that inflation ends on the plateau long before the potential becomes negative.
(a)(b)Notice that our requirements of final values for the quantitiesδeϕ and δeχ (both much bigger than 1/R) imposes fine-tuning constraints on the initial conditions, since at leading order in 1/R the values for γ*ϕ and γ*χ must coincide (see (62) and recall that p≃q).We stress that the enhancement of non-Gaussianity occurs towards the end of inflation. What happens justafter inflation is a model-dependent issue; we cannot address this question within the approximations used in this concrete model. For our choice of parameters, we notice that fNL is still increasing at the end of inflation, see Figure 1, but for other choices it may start to decrease before inflation ends, see [115]. We stress that in the regime where |α|,|β|≫1 the slow-roll parameters will necessarily become much greater than unity by the end of inflation, which may correspond to |γss|≫1: this is one reason by which fNL can become large. Our formalism is at least in principle suitable to study the evolution of non-Gaussianity after inflation, when the parameter ϵH becomes larger than unity (although we restate that the specific potential being considered here anyway needs to be modified after inflation ends).


## 6. Higher Order Non-Gaussian Observables

Non-Gaussianity is most popularly parameterised in terms of one nonlinearity parameter,fNL. There are several different fNL's used in the literature which parameterise different shape dependences of the bispectrum [41, 42], the most popular being the local model discussed in this paper and the so-called equilateral model which can arise from models of inflation with noncanonical kinetic terms (a popular example is DBI inflation [15–25]). But can we learn more than one number from an observation of non-Gaussianity? And if fNLlocal is detected how can we hope to distinguish between the models which can generate this? Fortunately, observations of non-Gaussianity have the possibility of teaching us much more than the value of one parameter. Two examples of further non-Gaussian observables are a scale dependence of fNL and the trispectrum.


### 6.1. Trispectrum

The four point function of the primordial curvature perturbation is defined by(67)〈ζk1ζk2ζk3ζk4〉e≡Tζ(k1,k2,k3,k4)(2π)3δ3(k1+k2+k3+k4),
which using the δN formalism, and assuming that the fields have a Gaussian distribution at Hubble exit, is given by(68)Tζ(k1,k2,k3,k4)=τNL[Pζ(k13)Pζ(k3)Pζ(k4)+(11perms)]+5425gNL[Pζ(k2)Pζ(k3)Pζ(k4)+(3perms)],
where k13≡|k1+k3|, and the trispectrum nonlinearity parameters can be calculated using [98, 106, 127]
(69)τNL=NABNACNBNC(NDND)3,gNL=2554NABCNANBNC(NDND)3.
Hence, we see that the trispectrum depends on two nonlinearity parameters (as opposed to one, fNL, for the bispectrum), and they may be observationally distinguishable since they are prefactors of terms with different shape dependences (68). The current observational bound on the local type of the bispectrum from seven years of WMAP data is -10<fNL<74 at the 2σ level [2]. Recently, there has been a first constraint on both τNL and gNL, which came from using WMAP5 data, with the bounds -7.4<gNL/105<8.2 and -0.6<τNL/104<3.3 at 95% confidence [128]. There have been two other observational constraints on the trispectrum through gNL, setting τNL=0. The bounds, which in both cases are roughly |gNL|≲  105–106, come from large scale structure [129], and the CMB [130]. If there is no detection of non-Gaussianity it is expected that with Planck data the bounds will be reduced to about |fNL|≲10, τNL≲103 and gNL≲105 at the 2σ level and future observations may become even tighter [131–135]. Note that the expected future constraint on gNL is about two orders of magnitude weaker than that on τNL [135]. We will see that it is also possible for some models and parameter ranges that the first observational signature will come through the trispectrum rather than the bispectrum.


#### 6.1.1. Trispectrum in Two-Field Hybrid Inflation

In the regime where|gNL| and τNL are greater than unity, they are given by(70)gNL=103r̃(ηφφ-2ηχχ)-ηχχ1+r̃fNL,(71)τNL=r̃(1+r̃)3ηχχ2e4N(ηφφ-ηχχ)=1+r̃r̃(65fNL)2.
We see that gNL is subdominant to fNL and hence will not provide a competitive observational signature. It follows from (71) that τNL>(6fNL/5)2, so τNL may be large and provide an extra observable parameter for this model. This inequality between τNL and fNL is true in general [75], and equality is reached whenever a single field direction during inflation generates the primordial curvature perturbation. However, it is usually assumed that τNL~fNL2 since both arise from second derivatives in the δN formalism. In fact, for our model it is possible to have a small fNL (and hence also a small gNL) but a large and potentially observable τNL. For this, we require that r̃≪1, although in practice if we make it too small it may no longer be possible to satisfy a constraint on a minimum possible field velocity consistent with a classical slow-roll trajectory, as discussed in [95] (see also [136] and a discussion which reaches a different conclusion is given in [137]). In the final example in Table 1, we give an explicit example of parameter values which give rise to an fNL which is probably too small to be detected with Planck but with a very large trispectrum through τNL>103 that should be detectable at a high significance. For another example with fNL,gNL≲O(1) but τNL≫1, see [76]. In contrast, it has been shown in several papers [127, 138–142] that in the curvaton scenario where the curvaton has a nonquadratic potential it is possible to realise |gNL|≫1 while τNL=(6fNL/5)2 is small with some tuning of parameters. This is also possible in the exact solution which we presented in Section 5 and we consider this next.


#### 6.1.2. Trispectrum in the Exact Solution

In this regime, where|α|≫|β|>R and m,n>1/R2, it is also possible to give compact expressions for the trispectrum (4-point function) nonlinearity parameters, in terms of fNL2 as (72)τNL≃(65fNL)2,5425gNL≃-32n-mm(65fNL)2.
The complete result was given in [115]. Hence, both of the trispectrum nonlinearity parameters are generally large whenever fNL is. We note that if m≃1 then from (60) n≫1, so the trispectrum through gNL will give the dominant signal of non-Gaussianity through a large, negative gNL.


### 6.2. Scale Dependence offNL

In its simplest form, the local form offNL defined by (1) is a constant parameter, independent of both position and scale. However, in realistic models fNL is likely to be mildly scale-dependent. This may happen in two ways. In general, single-field models (such as the curvaton scenario) where the scalar field which generates the primordial curvature perturbation has a nonquadratic potential the non-linearities this generates will give rise to scale dependence of fNL. Alternatively, even in models where all of the fields have a quadratic potential; if the primordial curvature perturbation has contributions from more than one field, and the fields do not all have the same mass, then fNL will again have a scale dependence. This is because the correlation between the first and second order terms of ζ, which the bispectrum depends on, will become scale-dependent. This is indeed what happens in the two-field hybrid model presented earlier. For details about a scale dependence of local fNL, see Byrnes et al. [143].We define a second observable parameter derived from the bispectrum(73)nfNL=dln|fNL|dlnk,
which is analogous to the scale dependence of the power spectrum. There is a subtlety here, which is that in general fNL may depend on three independent parameters, k1,k2, and k3. However, it was shown in [143] that provided one takes the derivative while keeping the ratio of the three k vectors fixed then nfNL is independent of the shape of the triangle described by the three k vectors, which makes this a well defined quantity. Observational prospects for this quantity were considered in [144], who showed that the Planck satellite is sensitive to nfNL≃0.1, assuming a fiducial value of fNL=50. This observational sensitivity to nfNL is about a factor of two larger than the current preferred value of the power spectrums spectral index and may therefore provide an interesting extra constraint on non-Gaussian models. The scale dependence of an equilateral form of fNL has also been considered from both a theoretical and observational perspective, see for example [144–147].In the example of two-field hybrid inflation,fNL has a scale dependence both because of the exponential term in fNL, (41), and because r̃ will vary through the change of the initial value of sin2θ*. We find (74)∂lnr̃∂lnk=∂lne2N(ηφφ-ηχχ)∂lnk=-2(ηφφ-ηχχ).
Using this we find from (41) that (75)nfNL≡dlogfNLdlogk=-4ηφφ-ηχχ1+r̃.
In the case that we include the effect from the surface where the waterfall field is destabilised and g12=g22, we find from (48) that (76)nfNL=-4ηφφ-ηχχ1+(ηφφ/ηχχ)2r̃.
For both cases, the spectral index of fNL satisfies (77)-4(ηφφ-ηχχ)<nfNL<0,
for any value of r̃ and hence fNL will be smaller on small scales.Because we require a relatively large value ofηφφ-ηχχ>1/N for our model to generate a large non-Gaussianity it is quite possible for our model to generate a relatively significant scale dependence of fNL. However, the amount also depends on r̃, and when this is large, then the χ field is almost solely responsible for generating ζ at both first and second order and nfNL is suppressed. We can also see in agreement with the statement at the beginning of this subsection that in this case, if the mass of the two fields are equal, then the two fields have the same scale dependence and nfNL=0.We note that this is in contrast to the large non-Gaussianity from an inhomogeneous end of inflation found in [93]. In the specific cases, they considered to generate a large non-Gaussianity the non-Gaussianity was generated purely at the end of inflation and fNL is scale-independent. In detail, we see from (4.4) and (4.24) in [93] that their formulas for fNL does not depend on N or on any quantities evaluated at Hubble exit. This is also in contrast to the exact solution considered in Section 5.1. Our expression for fNL in (66) depends on the initial values through p and q, and this dependence drops out at leading order in α. Hence, fNL is independent of the number of e-foldings and scale-independent.


## 6.1. Trispectrum

The four point function of the primordial curvature perturbation is defined by(67)〈ζk1ζk2ζk3ζk4〉e≡Tζ(k1,k2,k3,k4)(2π)3δ3(k1+k2+k3+k4),
which using the δN formalism, and assuming that the fields have a Gaussian distribution at Hubble exit, is given by(68)Tζ(k1,k2,k3,k4)=τNL[Pζ(k13)Pζ(k3)Pζ(k4)+(11perms)]+5425gNL[Pζ(k2)Pζ(k3)Pζ(k4)+(3perms)],
where k13≡|k1+k3|, and the trispectrum nonlinearity parameters can be calculated using [98, 106, 127]
(69)τNL=NABNACNBNC(NDND)3,gNL=2554NABCNANBNC(NDND)3.
Hence, we see that the trispectrum depends on two nonlinearity parameters (as opposed to one, fNL, for the bispectrum), and they may be observationally distinguishable since they are prefactors of terms with different shape dependences (68). The current observational bound on the local type of the bispectrum from seven years of WMAP data is -10<fNL<74 at the 2σ level [2]. Recently, there has been a first constraint on both τNL and gNL, which came from using WMAP5 data, with the bounds -7.4<gNL/105<8.2 and -0.6<τNL/104<3.3 at 95% confidence [128]. There have been two other observational constraints on the trispectrum through gNL, setting τNL=0. The bounds, which in both cases are roughly |gNL|≲  105–106, come from large scale structure [129], and the CMB [130]. If there is no detection of non-Gaussianity it is expected that with Planck data the bounds will be reduced to about |fNL|≲10, τNL≲103 and gNL≲105 at the 2σ level and future observations may become even tighter [131–135]. Note that the expected future constraint on gNL is about two orders of magnitude weaker than that on τNL [135]. We will see that it is also possible for some models and parameter ranges that the first observational signature will come through the trispectrum rather than the bispectrum.


### 6.1.1. Trispectrum in Two-Field Hybrid Inflation

In the regime where|gNL| and τNL are greater than unity, they are given by(70)gNL=103r̃(ηφφ-2ηχχ)-ηχχ1+r̃fNL,(71)τNL=r̃(1+r̃)3ηχχ2e4N(ηφφ-ηχχ)=1+r̃r̃(65fNL)2.
We see that gNL is subdominant to fNL and hence will not provide a competitive observational signature. It follows from (71) that τNL>(6fNL/5)2, so τNL may be large and provide an extra observable parameter for this model. This inequality between τNL and fNL is true in general [75], and equality is reached whenever a single field direction during inflation generates the primordial curvature perturbation. However, it is usually assumed that τNL~fNL2 since both arise from second derivatives in the δN formalism. In fact, for our model it is possible to have a small fNL (and hence also a small gNL) but a large and potentially observable τNL. For this, we require that r̃≪1, although in practice if we make it too small it may no longer be possible to satisfy a constraint on a minimum possible field velocity consistent with a classical slow-roll trajectory, as discussed in [95] (see also [136] and a discussion which reaches a different conclusion is given in [137]). In the final example in Table 1, we give an explicit example of parameter values which give rise to an fNL which is probably too small to be detected with Planck but with a very large trispectrum through τNL>103 that should be detectable at a high significance. For another example with fNL,gNL≲O(1) but τNL≫1, see [76]. In contrast, it has been shown in several papers [127, 138–142] that in the curvaton scenario where the curvaton has a nonquadratic potential it is possible to realise |gNL|≫1 while τNL=(6fNL/5)2 is small with some tuning of parameters. This is also possible in the exact solution which we presented in Section 5 and we consider this next.


### 6.1.2. Trispectrum in the Exact Solution

In this regime, where|α|≫|β|>R and m,n>1/R2, it is also possible to give compact expressions for the trispectrum (4-point function) nonlinearity parameters, in terms of fNL2 as (72)τNL≃(65fNL)2,5425gNL≃-32n-mm(65fNL)2.
The complete result was given in [115]. Hence, both of the trispectrum nonlinearity parameters are generally large whenever fNL is. We note that if m≃1 then from (60) n≫1, so the trispectrum through gNL will give the dominant signal of non-Gaussianity through a large, negative gNL.


## 6.1.1. Trispectrum in Two-Field Hybrid Inflation

In the regime where|gNL| and τNL are greater than unity, they are given by(70)gNL=103r̃(ηφφ-2ηχχ)-ηχχ1+r̃fNL,(71)τNL=r̃(1+r̃)3ηχχ2e4N(ηφφ-ηχχ)=1+r̃r̃(65fNL)2.
We see that gNL is subdominant to fNL and hence will not provide a competitive observational signature. It follows from (71) that τNL>(6fNL/5)2, so τNL may be large and provide an extra observable parameter for this model. This inequality between τNL and fNL is true in general [75], and equality is reached whenever a single field direction during inflation generates the primordial curvature perturbation. However, it is usually assumed that τNL~fNL2 since both arise from second derivatives in the δN formalism. In fact, for our model it is possible to have a small fNL (and hence also a small gNL) but a large and potentially observable τNL. For this, we require that r̃≪1, although in practice if we make it too small it may no longer be possible to satisfy a constraint on a minimum possible field velocity consistent with a classical slow-roll trajectory, as discussed in [95] (see also [136] and a discussion which reaches a different conclusion is given in [137]). In the final example in Table 1, we give an explicit example of parameter values which give rise to an fNL which is probably too small to be detected with Planck but with a very large trispectrum through τNL>103 that should be detectable at a high significance. For another example with fNL,gNL≲O(1) but τNL≫1, see [76]. In contrast, it has been shown in several papers [127, 138–142] that in the curvaton scenario where the curvaton has a nonquadratic potential it is possible to realise |gNL|≫1 while τNL=(6fNL/5)2 is small with some tuning of parameters. This is also possible in the exact solution which we presented in Section 5 and we consider this next.


## 6.1.2. Trispectrum in the Exact Solution

In this regime, where|α|≫|β|>R and m,n>1/R2, it is also possible to give compact expressions for the trispectrum (4-point function) nonlinearity parameters, in terms of fNL2 as (72)τNL≃(65fNL)2,5425gNL≃-32n-mm(65fNL)2.
The complete result was given in [115]. Hence, both of the trispectrum nonlinearity parameters are generally large whenever fNL is. We note that if m≃1 then from (60) n≫1, so the trispectrum through gNL will give the dominant signal of non-Gaussianity through a large, negative gNL.


## 6.2. Scale Dependence offNL

In its simplest form, the local form offNL defined by (1) is a constant parameter, independent of both position and scale. However, in realistic models fNL is likely to be mildly scale-dependent. This may happen in two ways. In general, single-field models (such as the curvaton scenario) where the scalar field which generates the primordial curvature perturbation has a nonquadratic potential the non-linearities this generates will give rise to scale dependence of fNL. Alternatively, even in models where all of the fields have a quadratic potential; if the primordial curvature perturbation has contributions from more than one field, and the fields do not all have the same mass, then fNL will again have a scale dependence. This is because the correlation between the first and second order terms of ζ, which the bispectrum depends on, will become scale-dependent. This is indeed what happens in the two-field hybrid model presented earlier. For details about a scale dependence of local fNL, see Byrnes et al. [143].We define a second observable parameter derived from the bispectrum(73)nfNL=dln|fNL|dlnk,
which is analogous to the scale dependence of the power spectrum. There is a subtlety here, which is that in general fNL may depend on three independent parameters, k1,k2, and k3. However, it was shown in [143] that provided one takes the derivative while keeping the ratio of the three k vectors fixed then nfNL is independent of the shape of the triangle described by the three k vectors, which makes this a well defined quantity. Observational prospects for this quantity were considered in [144], who showed that the Planck satellite is sensitive to nfNL≃0.1, assuming a fiducial value of fNL=50. This observational sensitivity to nfNL is about a factor of two larger than the current preferred value of the power spectrums spectral index and may therefore provide an interesting extra constraint on non-Gaussian models. The scale dependence of an equilateral form of fNL has also been considered from both a theoretical and observational perspective, see for example [144–147].In the example of two-field hybrid inflation,fNL has a scale dependence both because of the exponential term in fNL, (41), and because r̃ will vary through the change of the initial value of sin2θ*. We find (74)∂lnr̃∂lnk=∂lne2N(ηφφ-ηχχ)∂lnk=-2(ηφφ-ηχχ).
Using this we find from (41) that (75)nfNL≡dlogfNLdlogk=-4ηφφ-ηχχ1+r̃.
In the case that we include the effect from the surface where the waterfall field is destabilised and g12=g22, we find from (48) that (76)nfNL=-4ηφφ-ηχχ1+(ηφφ/ηχχ)2r̃.
For both cases, the spectral index of fNL satisfies (77)-4(ηφφ-ηχχ)<nfNL<0,
for any value of r̃ and hence fNL will be smaller on small scales.Because we require a relatively large value ofηφφ-ηχχ>1/N for our model to generate a large non-Gaussianity it is quite possible for our model to generate a relatively significant scale dependence of fNL. However, the amount also depends on r̃, and when this is large, then the χ field is almost solely responsible for generating ζ at both first and second order and nfNL is suppressed. We can also see in agreement with the statement at the beginning of this subsection that in this case, if the mass of the two fields are equal, then the two fields have the same scale dependence and nfNL=0.We note that this is in contrast to the large non-Gaussianity from an inhomogeneous end of inflation found in [93]. In the specific cases, they considered to generate a large non-Gaussianity the non-Gaussianity was generated purely at the end of inflation and fNL is scale-independent. In detail, we see from (4.4) and (4.24) in [93] that their formulas for fNL does not depend on N or on any quantities evaluated at Hubble exit. This is also in contrast to the exact solution considered in Section 5.1. Our expression for fNL in (66) depends on the initial values through p and q, and this dependence drops out at leading order in α. Hence, fNL is independent of the number of e-foldings and scale-independent.


## 7. Conclusions

We have reviewed various models which can generate a large local non-Gaussianity. A feature shared by all of these models is that they have more than one light scalar field present during inflation. This extra degree of freedom generates an isocurvature perturbation which is at least partially converted into the primordial curvature perturbation after horizon exit of the modes which are observable today. In the curvaton and modulated reheating scenarios, this conversion occurs after the end of inflation, while in the inhomogeneous end of inflation scenario, this conversion occurs on the non-uniform energy density hypersurface on which inflation ends. For these three scenarios, the light field which generates the primordial curvature perturbation after or at the end of inflation can be treated as a test field which does not affect the inflationary dynamics.Our main focus has been on models in which a large non-Gaussianity is generated during inflation. This can occur even within slow-roll inflation for certain potentials and certain trajectories. We have shown, at least in the case of a separable potential, that the trajectory is required to be almost entirely along the direction of one field but that the orthogonal field must become more important towards the end of inflation, and hence the inflationary trajectory must curve. In absolute terms, the change to the angle of the background trajectory is small (compared to a trajectory which turns by a right angle during inflation), but in relative terms it must grow by at least an order of magnitude. This is in contrast to the previous three scenarios.We have reviewed two-field hybrid inflation as an explicit model of a separable potential where the conditions required to generate a large non-Gaussianity can be satisfied. The conditions can be satisfied for any possible combination of positive and negativeη parameters, so the potential can be bowl-shaped, a hill-top, or have a saddle point. The main conditions which must be satisfied is that the difference of the two η parameters must not be too small, ηφφ-ηχχ~0.1, and the value of the χ field must be very subdominant to that of the φ field initially (or vice versa). In general, but depending on the coupling constants between the two inflaton fields and the waterfall field, there is a change to observables at the end of inflation, due to the fact that the surface on which the waterfall field is destabilised and inflation ends might not be a surface of uniform energy density. This effect is responsible for the inhomogeneous end of inflation scenario. It is then a model, dependent question whether there will be further evolution to the observables during reheating in this model, this deserves further attention.One similarity that this hybrid inflation model has together with the quadratic curvaton scenario is that in both cases the initially subdominant (approximately isocurvature) fieldχ has the ratio δχ/χ~ζχ approximately constant and the field fluctuations do not become more non-Gaussian with time. However, the effect of this light field on the primordial curvature perturbation grows, during inflation in the hybrid scenario and before the curvaton decay in the curvaton scenario. It is this nonlinear transfer between the field fluctuation and ζ, described by the δN formalism, which can generate a large non-Gaussianity. Therefore, the non-Gaussianity in the hybrid scenario which we have studied is generated on super horizon scales during slow-roll inflation, in a similar way to which non-Gaussianity is generated over time in the curvaton scenario before the decay of the curvaton. The evolution of fNL during inflation is explicitly calculated and plotted in [95]. This conclusion is somewhat different from that in [47], and we plan to elaborate on this point in a future work. For more discussion on the distinction between non-Gaussianity generated by a non-Gaussian field perturbation, and non-Gaussianity generated by a nonlinear transfer between a Gaussian field perturbation and ζ, see [109]. An example where the subdominant fields fluctuations can become non-Gaussian due to a large self interaction was discussed by Bernardeau [45].In order to study models where slow-roll breaks down before the end of inflation, it is clearly necessary to go beyond a formalism based on the slow-roll approximation. We have shown how this can be done in the context a separable Hubble parameter instead of a separable potential and this leads to an exact expression for (the local part of)fNL in these models. As an explicit example, an exact two-field solution with an exponential potential was given. For some parameter choices, this leads to a strong break down of slow roll before the end of inflation, which may give rise to a large non-Gaussianity. Further work is also required for this model to understand how the potential may be modified after the end of inflation in order that reheating occurs.Non-Gaussianity is a topical field, in which observations have improved greatly over the last decade through both studies of the CMB and large-scale structure. Observations so far have heavily focused on constraining the bispectrum nonlinearity parameterfNL. Currently, the tightest constraint comes from the WMAP satellite, assuming the local model of non-Gaussianity this constrains the amplitude of the non-Gaussian part of the primordial curvature perturbation to be less than about one thousandth the amplitude of the Gaussian perturbation. This constraint is likely to be tightened considerably by the Planck satellite, which is currently taking data, or instead there might be a detection. A detection of fNL at this level would rule out simplest models of inflation, which are single field with a canonical kinetic term. Clearly, this would be an extremely exciting result.However, even if we are in the fortunate position of having a detection offNL as well as improved constraints/detection of the scalar-to-tensor ratio and the spectral index, there will probably still be several viable scenarios, as detailed in this paper, which for suitable parameter choices and initial conditions can match the observations. Fortunately, non-Gaussianity is about much more than one number. The trispectrum (four-point function) depends on two nonlinearity parameters. In general, τNL≥(6fNL/5)2. If the current observational hints (which are not statistically significant) that fNL~40 turn out to be true, then both the bispectrum and the trispectrum should be large enough for Planck to detect. Even if the bispectrum turns out to be much smaller, although for many models τNL is close to the lower bound, we have seen that in the model of hybrid inflation it is possible to have τNL≫fNL2, so the trispectrum might even be the first observational signature of non-Gaussianity. Alternatively, the trispectrum through a large gNL might give the first observational signature, as is possible in self-interacting curvaton models or the exact solution with an exponential potential. If fNL is detected, it will also be possible to either constrain or detect a scale dependence of this parameter. Although it is often assumed to be constant, this is only true for certain simple models, and for example in the two-field hybrid inflation model it generally has a significant scale dependence. We have, therefore, seen that non-Gaussianity is an important and powerful method of constraining and distinguishing between the many models of inflation.



---

*Source: 724525-2010-09-16.xml*