NonGaussianity from violation of slowroll in multiple inflation
Abstract:
Multiple inflation is a model based on N=1 supergravity wherein there are sudden changes in the mass of the inflaton because it couples to ‘flat direction’ scalar fields which undergo symmetry breaking phase transitions as the universe cools. The resulting brief violations of slowroll evolution generate a nongaussian signal which we find to be oscillatory and yielding . This is potentially detectable by e.g. Planck but would require new bispectrum estimators to do so. We also derive a modelindependent result relating the period of oscillations of a phase transition during inflation to the period of oscillations in the primordial curvature perturbations generated by the inflaton.
1 Introduction
A major goal in modern observational cosmology is to detect any nongaussianity of the temperature anisotropies in the cosmic microwave background (CMB) [1, 2, 3]. All single field, slow roll, inflation models give rise to anisotropies that are gaussian [4]. Therefore detection of a significant nongaussian signal would falsify such models and provide new insights into the dynamics of inflation.
Quantifying such a signal is however not a straightforward task. Whilst gaussianity is a welldefined property, nongaussianity is not and the anisotropies can, in principle, deviate from gaussianity in many different ways [3]. Also, the CMB temperature anisotropies measured by WMAP are very close to gaussian [5]. Therefore any measure designed to detect nongaussianity needs to be sensitive to a very small signal.
One measure of nongaussianity that is particularly useful is the three point correlation function, or ‘bispectrum’ of the temperature anisotropies. Gaussian random variables (and functionals of gaussian random variables) have the property that all odd power correlation functions are zero. This makes the bispectrum the lowest order statistic for which any nonzero result would indicate a departure from gaussianity. The bispectrum contains much more information than the power spectrum as, in general, it depends on both scale and shape. Therefore, if a nonzero bispectrum is detected it will also be an extremely useful statistic for constraining models of the early universe.
However there are limitations to how much information can be inferred from the bispectrum:

Modern CMB experiments possess very large numbers of pixels of data. e.g. for WMAP and for Planck [6]. Considering that each higher order of correlation function will require an additional factor of calculations, exact determination of the bispectrum quickly becomes impossible. Therefore various estimators have been constructed [7], but each can look only for a specific type of bispectrum and so might miss a signal different from the type being searched for.

Higher order correlation functions suffer more from cosmic variance because more information is required at each scale to compute them. In this work we will deal with a scale dependent signal, therefore limitations due to cosmic variance are particularly relevant.
These problems are not necessarily insurmountable. Computing power will only continue to get stronger, making bispectrum estimators increasingly powerful with time. For each potential signal, estimators can also be constructed for that specific signal, ensuring processing time is used as efficiently as possible. The estimator can be optimised for detecting primordial nongaussianity while discriminating against secondary nongaussianities arising from e.g. unsubtracted point sources or residuals from component separation [8]. Finally, in models that have nontrivial scale dependence, a correlation will likely exist between the power spectrum and bispectrum as we illustrate in this paper.
One of the first methods for quantifying nongaussianity involved rewriting the Newtonian potential as [1],
(1) 
where is a constant, and is a gaussian variable. However, this parameterisation has a rather specific form and does not capture other possible deviations from gaussianity. A more general parameterisation can be obtained by allowing to depend explicitly on the wavevector, . In terms of the adiabatic curvature perturbation, , the parameterisation becomes [2, 9],
(2) 
Here, the product is used because has been allowed to depend on scale.^{1}^{1}1See Eq.(229) in Ref.[2] for how to define this product. The factor of 3/5 comes from the relationship between the adiabatic curvature and the Newtonian potential at matter domination.
For as defined in Eq.(1), the WMAP 5year constraint is [5]. This reinforces the point made already: the primordial temperature anisotropies are very close to gaussian. Given that the amplitude of these anisotropies is the nongaussian contribution to these anisotropies is at most (or of the overall anisotropy). In general the expectation from inflation is that the anisotropies should be close to gaussian, therefore this result can be considered as a success of the inflationary paradigm.
More precisely the expectation from inflation of nongaussianity as parameterised by Eq.(2) depends on the specific model. It is known that for single field, slowroll inflation with a canonical kinetic term and a vacuum initial state, the result is [4] where is the usual slowroll parameter defined later in Eq.(13). The best sensitivity expected from the Planck satellite is to of while the secondary contribution from postinflationary evolution is of [10]. Therefore the prediction of the simplest inflationary toy models ^{2}^{2}2By this we mean a generic finetuned potential such as the frequently used which has not yet been convincingly obtained from a physical theory, especially since inflation occurs in such models at . is that there should not be a detection of primordial nongaussianity in the near future.
Conversely a deviation from the simplest toy models can produce a larger value for . The study of these effects serves two purposes. Firstly, from the cosmological perspective, a detection of will immediately rule out single field slowroll inflation and focus attention on determining how inflation actually occured. Secondly, from the perspective of inflationary model building based on fundamental physics, as the bounds on grow tighter, some interesting models can be ruled out if there is no detection.
Much work has been done on nongaussianity generated due to multiple scalar fields (e.g. Ref.[11]), noncanonical kinetic terms (e.g. Ref.[12]) or nonvacuum initial states (e.g. Ref.[13]). However, surprisingly little attention has been paid to the possibility of nongaussianity generated by a violation of slowroll. Ref.[14] did consider this in a context when multiple fields are present; however, the nongaussianity itself is generated by the multiple fields, not the violation of slowroll (see also Ref.[15]). Refs.[16, 17] do consider the nongaussianity generated by violation of slowroll, but for a toy inflationary model. We follow their methods closely but consider a physical model of inflation. In both models the violation of slowroll generates sharp features in the power spectrum and also generates a ringing in the bispectrum with a characteristic period which we calculate below.
1.1 Features in the Power Spectrum of Primordial Fluctuations
It is an expectation in the simplest toy models of inflation that the power spectrum of temperature anisotropies in the CMB should be nearly scaleinvariant [18]. However this need not be the case for physical models. If the observed spectrum was indeed perfectly scaleinvariant, such models would be ruled out and naïvely this might seem to be the case. Assuming the “concordance” CDM cosmology, the primordial spectrum is well parameterised as a power law with spectral index [5]. This indicates that there cannot be any significant scale dependence that affects the spectrum over a large range of scales but it does not preclude the existence of e.g. localised oscillations in the spectrum, or other sharp features. Moreover since the observed anistrotopies arise from the convolution of the unknown primordial spectrum with the transfer function of the assumed cosmological model whose parameters are being determined, it is obvious that both unknowns cannot be determined simultaneously without further assumptions.
In fact, there are indications that the primordial spectrum might not be a scalefree power law, even assuming the CDM cosmology [19, 20, 21, 22, 23, 24]. At large angular scales () there is essentially no power and there are anomalous ‘glitches’, especially in the range of multipoles [25, 26]. The statistical significance of these anomalies is not sufficient to claim a definite detection, however it is strong enough to provide some tension with the fit to a scalefree powerlaw primordial spectrum which has only a 3% probability of being a good description of the WMAP 1year data [25], although this did improve to with the WMAP 3year data release [26]. Future measurements of the and mode polarisations by the Planck and (proposed) CMBPol satellites will throw light on whether these glitches are real or not [27].
“In the absence of an established theoretical framework in which to interpret these glitches (beyond the Gaussian, random phase paradigm), they will likely remain curiosities” [26]. Indeed if there were a model that had predicted, without ambiguity, the position and amplitude of the glitches, this would be seen as very strong evidence for the model. Although no such model exists, the general possibility of generating glitches over a range of scales had in fact been proposed prior to the WMAP observations in the context of ‘multiple inflation’ wherein the mass of the inflaton field undergoes sudden changes during inflation [28]. This generates characteristic localized oscillations in the spectrum, as was demonstrated numerically in a toy model of a inflationary potential with a ‘kink’ parameterised as [29]:
(3) 
By tuning the position, amplitude and gradient of the kink, the locations of the glitches can be varied to match the glitches seen in the power spectrum. One can then perform statistical likelihood tests to determine whether the fit is better with the glitches, but with the additional parameters, or with the simple scalefree spectrum. It was found by the WMAP team that the fit to the 1year data improves significantly (by ) for the model parameters , and , where GeV [30]. This analysis was repeated later using the WMAP 3year data, with similar results [31].
This seems encouraging, however in the toy model above is not the mass of the inflaton — in fact in all such monomial ‘chaotic’ inflation models with , inflation occurs at field values , hence the leading term in a Taylor expansion of the potential around is always linear in (since this is not a point of symmetry), rather than quadratic as for a mass term [32]. The effect of a change in the inflaton mass can be sensibly modelled only in ‘new’ inflation where inflation occurs at field values and an effective field theory description of the inflaton potential is possible. The ‘slowroll’ conditions are violated when the inflaton mass changes due to its (gravitational) coupling to ‘flat direction’ fields which undergo thermal phase transitions as the universe cools during inflation [28]. The resulting effect on the spectrum of the curvature perturbation was found by analytic solution of the governing equations to correspond to a ‘step’ followed by rapidly damped oscillations [33].^{3}^{3}3Although a similar phenomenon had been noted earlier for the case where the inflaton potential has a jump in its slope [34], such a discontinuity has no physical interpretation.
The next step should be to predict other observables, having used the power spectrum to constrain all the parameters in the model. In this paper we calculate the bispectrum of the multiple inflation model [28], using its predicted power spectrum [33] and the set of parameters which provide the best reduced in the CDM cosmology [35]. We also examine the effect on the bispectrum of varying these parameters over their full natural range.^{4}^{4}4We use the word ‘natural’ in this context to mean “stable towards radiative corrections”. The nongaussianity is found to be potentially detectable by the Planck satellite, or perhaps even a reanalysis of the WMAP data.
2 Multiple inflation
The biggest difficulty in inflationary model building is obtaining a potential that is both sufficiently flat to support inflation for efolds and stable towards radiative corrections. supergravity (SUGRA), the locally realised version of supersymmetry (SUSY), is a natural framework for achieving this (see Ref.[36] for a comprehensive review, especially of the cosmological issues discussed below). In SUSY there are usually many ‘flat directions’ in field space, i.e. scalar degrees of freedom which, while SUSY remains unbroken, have perfectly flat potentials. When SUSY is broken, the flat directions are lifted and usually acquire a masssquared related to the SUSY breaking scale. It would seem natural to identify one of these flat directions as the inflaton but to achieve a sufficient number of efolds of inflation, the masssquared of the inflaton needs to be much smaller than the Hubble expansion rate during inflation — if it is not, the evolution is too quick and inflation ends after only a few efolds. However, due to SUSY breaking by the vacuum energy driving inflation, the natural expectation is . This is commonly known as the problem in supergravity and superstring model building because it relates to the slowroll parameter (defined in Eq.14) being too large to support inflation.
There has been much work done to attempt to circumvent the problem (see Ref.[37] for a review and Ref.[38, 39] for some recent attempts in the framework of string/brane models). In this work we will assume that the inflaton mass has a symmetry protecting it against SUSYbreaking corrections. This can be done for example by having the inflaton in a ‘hidden sector’ where it interacts with the other sectors only gravitationally [40].
In the simplest inflationary models, scalar fields other than the inflaton itself are usually ignored. This is justified by the argument that, because of their larger masses, they will not contribute to the adiabatic density perturbations produced during inflation. If these fields decay before the end of inflation then their decay products will be diluted by the expansion of the universe and they will also not contribute any isocurvature perturbations. Such fields usually cannot couple to the inflaton either because any such coupling could endanger the flatness of the inflaton potential.
Multiple inflation [28] is based on SUGRA and assumes that while ‘flat direction’ scalar fields have their natural SUSYbreaking scale masses of order the Hubble parameter, there is just one scalar field (the inflaton) with a mass that is kept anomalously small through a symmetry. The amplitude of the potential during inflation will be the SUSY breaking scale (which need not be the electroweak scale as SUSY breaking can be different during and after inflation). The flat directions couple gravitationally to the inflaton and thus affect its evolution when their own masses change as the universe cools during inflation and they undergo symmetrybreaking phase transitions.
2.1 Evolution of the flat direction fields
It is assumed [28] that before inflation begins the universe is in a thermal state at temperature . At this point the flat directions, , are trapped at the origin by a thermal potential . Unless protected by a symmetry, the mass of the field will be of the order of the Hubble parameter, , which is itself determined by the scale of the vacuum energy driving inflation. Henceforth we work in units where the reduced Planck mass , therefore and .
Although the masssquared is negative, the flat directions will be lifted at large field values by nonrenormalisable operators (). For small field values, at nonzero temperatures, the full potential for the flat directions is thus [28],
(4) 
Before inflation, when the temperature is much larger than the mass of the field, it is trapped at the origin by a thermal barrier (quantum tunneling through which is negligible). When inflation starts, the temperature drops rapidly and after efolds of inflationary expansion the thermal barrier disappears. The field can then evolve to its minimum at and meanwhile there are efolds more of inflation [28].^{5}^{5}5The field evolves as so the occupation number of thermal states drops as . There is no nongaussianity generated due to the thermal fluctuations while the field is at the origin (c.f. [41]) since observable perturbations today leave the horizon much later. The field evolves (as a critically damped oscillator) according to the equation of motion:
(5) 
when it reaches the minimum the oscillations of the field are damped rapidly (within 1–2 cycles) by the term (see Fig.1 of [33]). Damping due to particle creation is far less effective since this is happening in an inflating background [28]. Consequently no isocurvature pertubations are generated by such oscillations.
The relationship of the flat directions to the inflaton
It was shown in Ref.[28] that a term in the Kähler potential with of is consistent with the underlying symmetries and corresponds to a term in the scalar potential where . The full potential, including the inflaton, is thus:
(6) 
There will be a running of the inflaton mass too, however the field value of the inflaton remains small throughout inflation so we can neglect the higherdimensional operators that would parameterise such running.
Hence as each field, , rolls to its minimum at , the effective inflaton masssquared changes by . The amplitude of the potential also changes as each field falls into its minimum, however, this change will be very small compared to the dominant term in the potential [28].
During periods of slowroll, the amplitude of the primordial perturbation spectrum can be expressed in terms of the potential and its derivatives as,
(7) 
Where the indicates that the quantities should be evaluated as the relevant scales ‘cross the horizon’, i.e. when . A transition in the mass of the inflaton, though not causing a significant change in , will cause a significant change in . If is positive, the mass will decrease causing a jump in the amplitude of the spectrum, while if it is negative then the amplitude will fall. As the flat direction oscillates in its minimum, the mass of the coupled inflaton will briefly oscillate as well. If the amplitude of the oscillations is large enough, one would expect a ‘ringing’ in the primordial power spectrum [35] for which there is tentative observational indication [21, 22, 23, 24].
The WMAP measurements indicate that, from through to the limit of set by signaltonoise considerations, the inferred primordial spectrum is nearly flat. However, as was noted earlier, there is tentative evidence of a step at very large scales of order the present Hubble radius [19, 21, 22]. It was in fact a prediction of the original paper on multiple inflation [28] that there should be of such feature in the efolds of inflationary history probed by the CMB.
The advantage of the multiple inflation model over the kink potential (3) considered in Ref.[16, 29] (or the toy model considered in Ref.[34]) is its grounding in a consistent, particle physics framework. While the latter model is able to improve the fit to the glitches in the data [30, 31] the constraints derived on the model parameters do not tell us anything about the relevant physical processes at these high energy scales. By contrast, constraints on multiple inflation are related directly to the masses and couplings of fields potentially present during inflation. For example, the scalar potential of the Minimal Supersymmetric Standard Model (MSSM) is flat along many directions in field space and a catalogue exists of all these directions [42]. These would be lifted during inflation due to SUSYbreaking by the large vacuum energy present and conceivably undergo symmetry breaking phase transitions.
3 Method
3.1 Single field justification
We follow closely the method of Ref.[16], who use the framework of Ref.[4] in calculating the bispectrum. This method uses the comoving gauge where, in single field inflation, the inflaton fluctuations are zero at all orders. When there are other, more massive, fields present this is no longer true because these extra fields will slightly perturb the instantaneous, adiabatic, direction in field space away from the direction of the main inflaton field. In multiple field scenarios a simpler gauge to describe the evolution of perturbations within the horizon is the uniform curvature gauge. In this gauge, neglecting gravitational waves, the slicing of spacetime into equal time hypersurfaces is made such that the spatial curvature on each of these surfaces is zero.
It was shown in Ref.[44] (and generalised beyond linear order in Ref.[45]) that there exists a gaugeindependent, nonperturbative quantity which, on sufficiently large scales, coincides with the adiabatic, scalar density perturbation at the perturbative level. Defining the number of efolds of expansion as , it was shown in Ref.[44, 45] that if we start in the uniform curvature gauge, this quantity can be expressed as . This allows us (on sufficiently large scales) to express the adiabatic curvature perturbation as [47],
(8) 
If we set the initial, uniform curvature, hypersurface a few efolds of expansion after the modes cross the horizon then due to the much smaller mass of the inflaton, we have , hence for all the flat direction fields. For all fields involved (including the inflaton), the higher order terms in Eq.(8) are also very small.^{6}^{6}6These terms become significant in models where nongaussianity is seeded by entropy perturbations, after the perturbations have left the horizon [46]. Under these conditions, the result holds to a very good approximation. The flat direction fields do not affect the adiabatic curvature perturbation and thus isocurvature perturbations do not arise gravitationally during inflation. Therefore the inflaton fluctuations will be zero in the comoving gauge during inflation and we are free to use the action derived in Ref.[4] using this gauge.
It should be stressed however that the flat directions do play a crucial role in the generation of the nongaussianities. The point of the preceding argument is that they do this only through their coupling to the inflaton and not directly through their coupling to gravity. What this means for the calculation is that when solving the zeroth order equations of motion for the inflaton we must include the effects of the flat direction fields. However, when computing the perturbative action to quadratic and cubic order in the perturbations, we need consider only the inflaton.
In principle it might have been less confusing to work initially in the uniform curvature gauge, where the contribution from all the fields is unambiguous, and to convert to the comoving gauge once the curvature perturbation has left the horizon. This would have the benefit of less ambiguity in the present case and be a necessity in genuine multifield models. In future work we intend to generalise the method of Ref.[16], as used in this paper, to inflationary models involving several, dynamically equally important, fields.
3.2 Calculating the bispectrum
We are interested in calculating the bispectrum (threepoint correlation function) of the quantum observable , corresponding to the scalar, adiabatic curvature perturbation. The model we are considering has the standard, purely gaussian vacuum state, and the quadratic action of the field, , is effectively just the free field action (i.e. a kinetic term and a mass term). A free gaussian field will remain gaussian, therefore in order to probe any nongaussianity in the fluctuations it is necessary to consider the cubic action. We know from observations of the power spectrum that i.e. higher terms in the perturbative expansion of the action will be less important in general.^{7}^{7}7This need not be true always, e.g. if a symmetry forces the cubic action to be anomalously small then the quartic action will produce the dominant contribution to any nongaussianity [48]. Therefore, we need only consider the action up to cubic order.
We proceed by rewriting the quantum operator in the interaction picture:
(9) 
where is the free field Hamiltonian derivable from the quadratic action and is the initial vacuum state value of the operator. The full quantum operator is then,
(10) 
where,
(11) 
is the time evolution operator for the interaction terms (and denotes the time ordering operator).
Putting these terms into the threepoint function gives,
(12)  
Where the last line of equality holds to first order in the expansion of the exponentials.
It is usually emphasised at this point that the vacuum with respect to which the expectation value is being taken is the fully interacting vacuum, not the free vacuum. This is important because the appearing in these equations will be expressed in terms of the ladder operators that annihilate the free vacuum. Fortunately, this is not an issue in this particular case because, similarly to what is done in Minkowski space, one can slightly deform the integral contour into Euclidean space. This projects out the free vacuum, as explicitly shown in Ref.[9].
Up to field redefinitions proportional to the free field equations of motion, the Hamiltonian in Eq.(12) can be obtained from the action derived in Ref.[4] (Eq.(3.9) in that paper and Eq.(3.4) in Ref.[16]). This action can be expressed as a sum of terms whose coefficients are functions of the slowroll parameters, defined as:
(13) 
(14) 
where the approximate equality holds when both and are small. Note that there are other common definitions for and in terms of derivatives of the potential: and ; the ones defined above are related through and . These equalities hold in the slowroll approximation when all four quantities are small.
In terms of these quantities, the cubic action for is:,
(15)  
where .
The term is proportional to the free field equations of motion, , and in Ref.[4] is removed by a field redefinition. This can have an important effect on the bispectrum, however, in the present case it does not, due to its relative size. There is just one term in (Eq.(3.10) in Ref.[4]) which does not include a derivative on one of the terms. As we can choose to apply the correction due to the redefinition after the relevant scales have left the horizon, only this term will have a nonzero effect — its coefficient is , which is much smaller than our leading term.^{8}^{8}8Note that although does become large during the phase transition in multiple inflation, we can choose to apply the correction after the phase transition is over, at which point is again . The redefined field will evolve outside the horizon while is significant so this can be done only if we also calculate the evolution of the redefined bispectrum out to this time, which we have indeed done.
With the field redefinition taken care of, the Hamiltonian in Eq.(12) just becomes . To use this to evaluate Eq.(12) we need to solve the commutator in this equation. We are free to use the standard commutation relations for ; to do so we first need to express in terms of the ladder operators of the free field vacuum
(16) 
with the usual commutation relations, and mode functions . The relationship between and is the usual
(17) 
When these definitions are substituted back into Eq.(12), the result is proportional to an integral over time of the free field mode functions , evaluated on the interaction Hamiltonian, (3.2.1).
The dominant term in multiple inflation
It has been stated already that multiple inflation violates slowroll, nevertheless remains throughout multiple inflation for all parameter values considered here.^{9}^{9}9If does become large, inflation will stop; this is potentially possible if there is a later inflationary epoch [28] but we do not consider this possibility here. There is only one term in the action (15) that does not have a factor of at least in it — this is the term with (note that is of order ) . However although remains small, due to the transition in the inflaton mass, and can temporarily become large. The size of the latter term is dictated by both the magnitude and the rate of the mass change in the phase transition. Multiple inflation is a small field model, hence is very small () and the term dominates over all others because it has one less factor of .
3.2.1 Numerical method for calculating bispectrum
Upon substituting this leading term into Eq.(12), using the definitions of and and working through the commutation, the result is (in conformal time ):
(18) 
This should be compared to Eq.(3.21) in Ref.[16].
To obtain the bispectrum for multiple inflation, we must solve this integral. To do this we first need to numerically solve for the free field mode functions , the scale factor and the slow roll parameters and . Once this has been done, the results are put into the above integral and evaluated out to a value of after the end of inflation. All that is needed after that is a rescaling of the full bispectrum to obtain a scale dependent generalisation of .
Fortunately, the slow roll parameters in the above integral are composite variables of other variables more fundamental to the model (e.g. ). Therefore, we need only solve for these more fundamental variables and then evaluate and from them. The full list of variables that need to be solved for, with their equations of motion, are [49]:
(19) 
with , , and:
(20)  
(21)  
(22) 
In all of the above, the indicates differentiation with respect to conformal time. These four equations are just the free field equation of motion for the observable , the two zeroth order equations of motion for and , and finally the Friedmann equation for the Hubble parameter .
We solve these equations numerically, using the Matlab function ode113 [50]. For initial conditions we take the BunchDavies vacuum state [49]. This amounts to:
(23) 
(24) 
The initial value of in the integrand was set such that integration of the mode functions, , begins when the following condition is first satisfied:
At earlier times, the mode functions are highly oscillatory and would cancel in the integral in Eq.(3.2.1). In fact there is only a small window in which this integral is nonzero — very early on, the mode functions are highly oscillatory and any contribution to this integral will be washed out, while very late on, the mode functions are frozen and thus the term will be zero. Finally, when the phase transition is not occurring, the term will be negligible. It is only those modes that are leaving the horizon during the phase transition that will give a a nonzero.contribution to the bispectrum.
A regularisation scheme is used in Ref.[16] to counter a sharp cutoff in the integral in Eq.(3.2.1) (and improved upon in Ref.[17]). We find that a regularisation scheme is useful for calculational efficiency but not necessary for accuracy. It was argued in Ref.[16] that although the integral is technically convergent as , this sharp cutoff would add a spurious contribution to the bispectrum, due to the highly oscillatory nature of the mode functions at early times. We find that this is a problem only when the cutoff also coincides with the violation of slowroll. If this occurs, there was indeed a spurious contribution due to the sharp nature of the cutoff; however, if we begin our integration earlier, when the slowroll contribution to the integral, is negligible, no spurious contribution is obtained.
In ordinary slowroll, to accurately calculate the contribution from the other, dominant terms in Eq.(15) by numerical means, a regularisation scheme would be necessary. For multiple inflation (and for the kink model considered in Ref.[16]) these other terms are small enough to be ignored. This is especially true in multiple inflation where is exceptionally small. What the regularisation scheme does do is allow the integration to begin later without sacrificing accuracy. Nevertheless, we did not use a regularisation in calculating our results and altered the integration start time when there was a risk of a spurious contribution.
3.2.2 Comparing a bispectrum to
We have defined a scaledependent generalisation of the parameter in Eq.(2). Now we discuss how to calculate from a given bispectrum. For a flat bispectrum a convention has been well established, however there is an ambiguity in the present case due to the scale dependence of both the power spectrum and the bispectrum. For a flat spectrum, following Refs.[4, 9], we start by reexpressing the bispectrum as,
(25) 
In this notation, we can rewrite in the form ^{10}^{10}10This follows directly from Eq.(2) if we rewrite it in space and then evaluate the bispectrum.
(26) 
However, as pointed out in Ref.[16], this definition entangles some of the ringing of the power spectrum in the definition of . The factors of in the bispectrum that are divided out to get are proportional to the power spectrum which is itself oscillatory. This can be overcome by factoring out a number that is equal to the square of the value of the power spectrum that is measured when assuming scale invariance [16]. This gives
(27) 
In the case of slowroll and a flat spectrum this definition is equivalent to our previous one in Eq.(25) (with ). It is then appropriate to define in the same way as in Eq.(26) and it is this that we discuss in the following section. Note that although Refs.[4, 9, 16] define differently, the definitions for coincide. We have followed the convention of Ref.[9] for .
4 Results
4.1 Reproduction of the kink model
To check the accuracy of our code and calculations we first set out to reproduce the bispectrum of the ‘kink’ model in Eq.(3) [16]. Note that dictates where the kink occurs, dictates the magnitude of the total shift in slope and dictates the rate of the change in slope. For models where the nongaussianity is generated as the modes cross the horizon, the bispectrum will be peaked over the equilateral shape (). The authors of Ref.[16] consider the statistic , where ; in the equilateral case this becomes . In Fig.1 we show the equilateral bispectrum for the kink model over the full range over which it is nonzero, to be compared with Fig.4 of Ref.[16] which shows this over a narrower range in . Note that although becomes as large as 10, it does not remain large over an extensive range of . This has important consequences for any attempt to detect the nongaussianity.
4.2 Bispectrum of multiple inflation
The most important, and difficult, part of the bispectrum calculation is obtaining the mode functions, . Therefore we first demonstrate that we can reproduce the power spectrum for multiple inflation accurately as shown in Fig.2 (to be compared with Fig.2 of Ref.[35]).
These spectra were calculated considering only one flat direction field, with a positive coupling () to the inflaton. The parameters used in the multiple inflation potential (6) are , , , and [35]. Note that is fixed through Eq.(7) by requiring the amplitude of perturbations to be the observed value () after the phase transition, and is different for each value of , the order of the nonrenormalisable operator which lifts the flat direction potential at large values of the field. The time at which the phase transition occurs also varies with . We emphasise that the value is particularly favoured as it is a well known flat direction in the MSSM scalar potential [42].
We can now confidently use our code to calculate the bispectrum of multiple inflation as shown in Fig.3. Note that as , the bispectrum begins to exhibit noise — this is a consequence of our not applying the regularisation scheme described in Ref.[16] and is only present at these scales for the reasons explained earlier. Secondly, the nongaussianity, while not as large as in the kink model, is nonzero over a wider range of . This can be traced to the fact that being a physical rather than a toy model, the features in the power spectrum generated by multiple inflation [35] are not as sharp as in the kink model [29]. Because the dominant term in the action is proportional to , if the transition is not sharp, then the contribution to the bispectrum will not be large either. In multiple inflation the field evolves to its minimum analogous to a critically damped oscillator and oscillates for a period determined by the Hubble damping, which means that the bispectrum is nonzero for longer.
Due to cosmic variance (the transition occurs at large scales) and secondary nongaussianities, it would be difficult to detect this signal even with a perfect estimator. We have a similar concern about the bispectrum of the kink model [16] because although it has a larger amplitude, it is nonzero over only a decade in making it very hard to acquire sufficient statistical significance for a convincing detection.
4.2.1 The effect of varying parameters
In Ref.[35], when finding the bestfit for multiple inflation to the WMAP and data, the only parameters varied were the amplitude of the potential, the location of the phase transition and the power of the nonrenormalisable term lifting the flat direction. All other parameters were held fixed at their most natural values, in order to maximise the predictive power of the model. Unlike the toy kink model [16] or other empirical models [43], the parameters of multiple inflation are not free to range over any set of values because of its grounding in the effective field theory framework. If one simply allowed all parameters to vary and then performed a naïve test, the physically constrained (and thus predictive) nature of multiple inflation would not be factored in.
We note in passing that there is a tension between the two methods of determining the likelihood of a model being correct. From the fundamental physics perspective one can calculate the parameter values that give the best fit and then estimate how natural they are. From the cosmology model fitting perspective one can determine the goodness of fit, given the number of parameters needed to achieve it. Neither test is fully satisfactory because natural valued parameters could give a bad fit and a good fit could come from an unnatural model. A Bayesian treatment that can quantify the naturalness of a model as a prior probability distribution would solve this problem. The resulting Bayesian likelihood would incorporate both the naturalness of the underlying model and the goodness of the fit to the data.
Irrespective of this issue, if multiple inflation is true, we do not know for certain what the values of the parameters will be. Therefore we have explored the effect of varying the parameters , and in turn and show the results in Figs. 4, 5 and 6. In each plot, we change as we vary the corresponding parameter to ensure that the power spectrum is always normalised to its observed value of . For ease of comparison we also alter the starting point of the phase transition so that the oscillations in each figure begin at the same point (this does not alter the shape of the bispectrum in any way). Each plot has been approximately equated to spatial scales probed by the CMB observations, although we are not here comparing with actual data.
The simplest parameter to consider is which does not affect the field’s dynamics but changes the effective mass of the inflaton through , so increasing increases the magnitude of the oscillations. Therefore we expect that increasing has no effect on the bispectrum except to increase its amplitude and this is just what we see in Fig.4. It is clear that a large bispectrum can be generated from natural values of the parameter (although values above can probably be ruled out already from the power spectrum).
4.2.2 The relationship between and
To consider the effects of and on the bispectrum it is necessary first to derive an identity relating the period of the oscillations over time during the phase transition to the period of the oscillations over in the power spectrum and bispectrum. Eq.(19) describes the time evolution of the mode functions of the operator . In the slowroll approximation, . When this quantity is small, oscillates as , and when this quantity is large we have , with constant. To a good approximation, the value of the constant is determined by equating the two solutions when . This approximation will continue to hold even if so long as continues to grow exponentially.
The quantity can be rewritten as . Since is almost constant we can write and . If there are oscillations induced on they will be seen more sensitively by higher derivatives of this quantity, hence will have a relative amplitude of oscillation much greater than . This means that will oscillate with the same frequency as , which in turn will oscillate with the same frequency as in its minimum.
In the slowroll limit, grows as and this holds even during the phase transition since the vacuum energy remains sensibly constant.^{11}^{11}11This will of course not hold precisely, however the relative deviation will be much smaller than the deviation of from its slowroll value, due to the term in . Therefore, if the oscillation repeats itself every units of time, then it will repeat itself every units of . It was established earlier that this means will also repeat itself every units of . We can thus parameterise the deviation of from by the equation:
(28) 
where is a function with period . When we do the matching of the solutions inside and outside the horizon and equate , we get,
(29) 
If the period of is then the period of will be because it will repeat itself half as often. If we ask what the difference in will be in the matching solution when undergoes one period, from to , with , we get,
(30) 
Unfortunately this result is not independent of the values of and and will vary as the oscillations occur. However if we consider what will happen to during this same time period we get, with :
(31) 
Now the only time dependence comes in the difference between the function at the two repeating points. These points were chosen because they were at the same point of the oscillation in , therefore to a good approximation will be equal at each point. This approximation will only hold exactly if the oscillations in are exactly repeating.
It is thus expected that the mode functions will oscillate with a period over equal to , with a small drift due to the change of the function during each oscillation. This result is general to any oscillatory phase transition during inflation and is not specific to multiple inflation.
It is a simple matter now to calculate the expected period, , of the oscillations in multiple inflation. The minimum of the potential for occurs at the point, . If we expand the potential around this point, with respect to the variable , then the potential becomes,
(32) 
It is useful to note that drops out of the quadratic term and only appears in the higher order terms. Seen from this perspective, the equation of motion for the field , when in its minimum, is to first approximation a damped harmonic oscillator. The period, , of a damped harmonic oscillator with damping term and potential is simply:
From the results earlier we expect then that the period of oscillations in both the power spectrum and bispectrum should depend only on the ratio and the power, , of the nonrenormalisable term lifting the potential of the flat direction, . A brief check of the period of confirms the expectation that for the mode functions.^{12}^{12}12 should have period because it involves the square of the mode functions . To a good approximation, the bispectrum will have a period ; however the bispectrum is more complicated because the are the mode functions of the , which are the free field operators. The bispectrum is a correlation function between the full operators and there will be additional oscillations generated by the terms. This can be read off Fig.2 for the parameters listed in the figure.
Therefore we expect that the result of increasing will be to decrease the period of oscillations in the bispectrum. Since this means more rapid oscillations, we also expect a magnification in during the phase transition, hence a magnification in the amplitude of the bispectrum. This is just what we see in Fig.5 (there is also a small effect on the shape of the bispectrum).
Based on this calculation we expect that changing will have only a secondary effect on the bispectrum. Indeed in Fig.6 the amplitude changes from to when is changed by over an order of magnitude. By comparison, when is similarly altered, the change in the bispectrum amplitude was from to . The increase in the amplitude is because changing changes the minimum of the field, which alters the magnitude of the change in the mass of the inflaton.
4.3 The bump model of [35]
Finally, we point out that in Ref.[35] a second example of multiple inflation was considered with two flat directions which have opposite sign couplings to the inflaton and produce a ‘bump’ in the primordial power spectrum. Using this model it was found that the WMAP 3year data can be fitted by an Einstein deSitter cosmology without a cosmological constant. We have calculated the bispectrum for this model which has , , , , , and find the amplitude of the oscillations induced by the term to be only . Such a small value will be swamped by secondary non gaussianities and could never be observed.
5 Conclusions
There is tentative evidence that the primordial power spectrum of scalar perturbations is imprinted with sharp features on large scales [25, 26]. These could have resulted from one or more phase transitions early on in the inflationary epoch as happens in the multiple inflation model [28]. This model is well motivated by fundamental physics ( supergravity) and such spectral features were predicted before WMAP provided observational indications for them.
It was shown [16] that any such departure from slowroll during inflation should also generate nongaussianity. We have calculated the nongaussianity in mutliple inflation for the parameter values that best fit the features in the power spectrum [35]. This is on the edge of being observable but a clear detection would require a new type of bispectrum estimator due to its scale dependence and non factorisability [7]. Significantly larger nongaussianities can be generated during multiple inflation if the parameters are allowed to range over values which are technically natural (i.e. stable towards radiative corrections).
The form of this nongaussianity is tightly correlated with the power spectrum — the oscillations in the bispectrum should begin at the same multipole as the oscillations in the power spectrum and have two thirds of the period. There have been a number of attempts, to deconvolve the primordial power spectrum directly from the WMAP data [20, 21, 22] (see also Ref.[51]). It is generally found that there is a supression of power on the scale of the present Hubble radius, followed by a ‘ringing’ at medium scales. By measuring the period of the oscillations in the power spectrum one can predict the period (and phase) of oscillations in the bispectrum, in a model independent manner.
Forthcoming measurements of CMB polarisation by Planck ought to shed light on whether these features in the power spectrum are systematic errors or genuine evidence of nontrivial dynamics during the inflationary era. The associated nongaussian signal should also be detectable according to our model calculation and provide insight into the dynamics.
5.1 Acknowledgements
SH acknowledges a Clarendon Fellowship and support from Balliol College, Oxford and the EU Marie Curie Network “UniverseNet” (HPRNCT2006035863). We thank Paul Hunt, David Lyth, Graham Ross, Misao Sasaki and David Wands for helpful comments and discussions.
References
 [1] E. Komatsu and D. N. Spergel, Phys. Rev. D 63 (2001) 063002.
 [2] N. Bartolo, E. Komatsu, S. Matarrese and A. Riotto, Phys. Rept. 402 (2004) 103.
 [3] X. Chen, M. x. Huang, S. Kachru and G. Shiu, J. Cosmo. Astropart. Phys. 01 (2007) 002.
 [4] J. M. Maldacena, J. High Energy Phys. 05 (2003) 013.
 [5] E. Komatsu et al. [WMAP Collaboration], Astrophys. J. Suppl. 180 (2009) 330.
 [6] M. A. J. Ashdown et al. [Planck CTP Collaboration], Astron. Astrophys. 467 (2007) 761.
 [7] J. R. Fergusson and E. P. S. Shellard, Phys. Rev. D 80 (2009) 043510.
 [8] D. Munshi and A. Heavens, \arXivid0904.4478 [astroph.CO].
 [9] D. Seery and J. E. Lidsey, J. Cosmo. Astropart. Phys. 06 (2005) 003.
 [10] E. Komatsu et al., \arXivid0902.4759 [astroph.CO].
 [11] M. Sasaki, J. Valiviita and D. Wands, Phys. Rev. D 74 (2006) 103003.
 [12] D. Langlois, S. RenauxPetel, D. A. Steer and T. Tanaka, Phys. Rev. D 78 (2008) 063523.
 [13] J. Martin, A. Riazuelo and M. Sakellariadou, Phys. Rev. D 61 (2000) 083518.
 [14] C. T. Byrnes and G. Tasinato, J. Cosmo. Astropart. Phys. 08 (2009) 016.
 [15] Y. F. Cai and H. Y. Xia, Phys. Lett. B 677 (2009) 226.
 [16] X. Chen, R. Easther and E. A. Lim, J. Cosmo. Astropart. Phys. 06 (2007) 023.
 [17] X. Chen, R. Easther and E. A. Lim, J. Cosmo. Astropart. Phys. 04 (2008) 010.
 [18] D. H. Lyth and A. Riotto, Phys. Rept. 314 (1999) 1.
 [19] J. Martin and C. Ringeval, Phys. Rev. D 69 (2004) 083515; Phys. Rev. D 69 (2004) 127303.
 [20] N. Kogo, M. Sasaki and J. Yokoyama, Phys. Rev. D 70 (2004) 103001; Prog. Theor. Phys. 114 (2005) 555.
 [21] A. Shafieloo and T. Souradeep, Phys. Rev. D 70 (2004) 043523; A. Shafieloo, T. Souradeep, P. Manimaran, P. K. Panigrahi and R. Rangarajan, Phys. Rev. D 75 (2007) 123502.
 [22] D. TocchiniValentini, M. Douspis and J. Silk, Mon. Not. R. Astron. Soc. 359 (2005) 31; D. TocchiniValentini, Y. Hoffman and J. Silk, Mon. Not. R. Astron. Soc. 367 (2006) 1095.
 [23] G. Nicholson and C. R. Contaldi, J. Cosmo. Astropart. Phys. 07 (2009) 011.
 [24] K. Ichiki and R. Nagata, Phys. Rev. D 80 (2009) 083002.
 [25] D. N. Spergel et al., Astrophys. J. Suppl. 148 (2003) 175.
 [26] G. Hinshaw et al. [WMAP Collaboration], Astrophys. J. Suppl. 170 (2007) 288.
 [27] M. J. Mortonson, C. Dvorkin, H. V. Peiris and W. Hu, Phys. Rev. D 79 (2009) 103519.
 [28] J. A. Adams, G. G. Ross and S. Sarkar, Nucl. Phys. B 503 (1997) 405.
 [29] J. Adams, B. Cresswell and R. Easther, Phys. Rev. D 64 (2001) 123514.
 [30] H. V. Peiris et al., Astrophys. J. Suppl. 148 (2003) 213.
 [31] L. Covi, J. Hamann, A. Melchiorri, A. Slosar and I. Sorbera, Phys. Rev. D 74 (2006) 083509; J. Hamann, L. Covi, A. Melchiorri and A. Slosar, Phys. Rev. D 76 (2007) 023503.
 [32] G. German, G. G. Ross and S. Sarkar, Phys. Lett. B 469 (1999) 46; Nucl. Phys. B 608 (2001) 423.
 [33] P. Hunt and S. Sarkar, Phys. Rev. D 70 (2004) 103518.
 [34] A. A. Starobinsky, JETP Lett. 55 (1992) 489.
 [35] P. Hunt and S. Sarkar, Phys. Rev. D 76 (2007) 123504.
 [36] D. J. H. Chung, L. L. Everett, G. L. Kane, S. F. King, J. D. Lykken and L. T. Wang, Phys. Rept. 407 (2005) 1.
 [37] L. Randall, \arXividhepph/9711471.
 [38] S. Kachru, R. Kallosh, A. Linde, J. M. Maldacena, L. P. McAllister and S. P. Trivedi, J. Cosmo. Astropart. Phys. 0310 (2003) 013.
 [39] D. A. Easson and R. Gregory, Phys. Rev. D 80 (2009) 083518.
 [40] G. G. Ross and S. Sarkar, Nucl. Phys. B 461 (1996) 597; J. A. Adams, G. G. Ross and S. Sarkar, Phys. Lett. B 391 (1997) 271.
 [41] S. Das and S. Mohanty, Phys. Rev. D 80 (2009) 123537.
 [42] T. Gherghetta, C. F. Kolda and S. P. Martin, Nucl. Phys. B 468 (1996) 37.
 [43] M. Joy, V. Sahni and A. A. Starobinsky, Phys. Rev. D 77 (2008) 023514; M. Joy, A. Shafieloo, V. Sahni and A. A. Starobinsky, J. Cosmo. Astropart. Phys. 0906 (2009) 028.
 [44] M. Sasaki and E. D. Stewart, Prog. Theor. Phys. 95 (1996) 71.
 [45] D. H. Lyth, K. A. Malik and M. Sasaki, J. Cosmo. Astropart. Phys. 05 (2005) 004.
 [46] D. H. Lyth, C. Ungarelli and D. Wands, Phys. Rev. D 67 (2003) 023503.
 [47] D. H. Lyth and Y. Rodriguez, Phys. Rev. Lett. 95 (2005) 121302; D. Seery and J. E. Lidsey, J. Cosmo. Astropart. Phys. 0509 (2005) 011.
 [48] K. T. Engel, K. S. M. Lee and M. B. Wise, Phys. Rev. D 79 (2009) 103530.
 [49] V. F. Mukhanov, H. A. Feldman and R. H. Brandenberger, Phys. Rept. 215 (1992) 203.
 [50] L. F. Shampine and M. W. Reichelt, SIAM Journal on Scientific Computing 18 (1997) 1.
 [51] M. Bridges, A. N. Lasenby and M. P. Hobson, Mon. Not. R. Astron. Soc. 369 (2006) 1123; Mon. Not. R. Astron. Soc. 381 (2007) 68.