We study the impact of different bias and redshift-space models on the halo power spectrum, quantifying their effect by comparing the fit to a subset of realizations taken from the WizCOLA suite. These provide simulated power spectrum measurements between $k_{\rm min}$ = 0.03 h/Mpc and $k_{\rm max}$ = 0.29 h/Mpc, constructed using the comoving Lagrangian acceleration method. For the bias prescription we include (i) simple linear bias; (ii) the McDonald & Roy model and (iii) its coevolution variant introduced by Saito et al.; and (iv) a very general model including all terms up to one-loop and corrections from advection. For the redshift-space modelling we include the Kaiser formula with exponential damping and the power spectrum provided by (i) tree-level perturbation theory and (ii) the Halofit prescription; (iii) one-loop perturbation theory, also with exponential damping; and (iv) an effective field theory description, also at one-loop, with damping represented by the EFT subtractions. We quantify the improvement from each layer of modelling by measuring the typical improvement in chi-square when fitting to a member of the simulation suite. We attempt to detect overfitting by testing for compatibility between the best-fit power spectrum per realization and the best-fit over the entire WizCOLA suite. For both bias and the redshift-space map we find that increasingly permissive models yield improvements in chi-square but with diminishing returns. The most permissive models show modest evidence for overfitting. Accounting for model complexity using the Bayesian Information Criterion, we argue that standard perturbation theory up to one-loop, or a related model such as that of Taruya, Nishimichi & Saito, coupled to the coevolution bias model, is likely to provide a good compromise for near-future galaxy surveys operating with comparable $k_{\rm max}$.

We extend the public CppTransport code to calculate the statistical properties of fluctuations in multiple-field inflationary models with curved field space. Our implementation accounts for all physical effects at tree-level in the 'in-in' diagrammatic expansion. This includes particle production due to time-varying masses, but excludes scenarios where the curvature perturbation is generated by averaging over the decay of more than one particle. We test our implementation by comparing results in Cartesian and polar field-space coordinates, showing excellent numerical agreement and only minor degradation in compute time. We compare our results with the PyTransport 2.0 code, which uses the same computational approach but a different numerical implementation, finding good agreement. Finally, we use our tools to study a class of gelaton-like models which could produce an enhanced non-Gaussian signal on equilateral configurations of the Fourier bispectrum. We show this is difficult to achieve using hyperbolic field-space manifolds and simple inflationary potentials.

The use of Eulerian 'standard perturbation theory' to describe mass assembly in the early universe has traditionally been limited to modes with k $\lesssim$ 0.1 h/Mpc at z = 0. At larger k the SPT power spectrum deviates from measurements made using N-body simulations. Recently, there has been progress in extending the reach of perturbation theory to larger k using ideas borrowed from effective field theory. We revisit the computation of the redshift-space matter power spectrum within this framework, including for the first time the full one-loop time dependence. We use a resummation scheme proposed by Vlah et al. to account for damping of baryonic acoustic oscillations due to large-scale random motions and show that this has a significant effect on the multipole power spectra. We renormalize by comparison to a suite of custom N-body simulations matching the MultiDark MDR1 cosmology. At z = 0 and for scales k $\lesssim$ 0.4 h/Mpc we find that the EFT furnishes a description of the real-space power spectrum up to $\sim$2%, for the $\ell$ = 0 mode up to $\sim$5%, and for the $\ell$ = 2, 4 modes up to $\sim$25%. We argue that, in the MDR1 cosmology, positivity of the $\ell$ = 0 mode gives a firm upper limit of k $\approx$ 0.75 h/Mpc for the validity of the one-loop EFT prediction in redshift space using only the lowest-order counterterm. We show that replacing the one-loop growth factors by their Einstein-de Sitter counterparts is a good approximation for the $\ell$ = 0 mode, but can induce deviations as large as 2% for the $\ell$ = 2, 4 modes. An accompanying software bundle, distributed under open source licenses, includes Mathematica notebooks describing the calculation, together with parallel pipelines capable of computing both the necessary one-loop SPT integrals and the effective field theory counterterms.

Over the next decade, improvements in cosmological parameter constraints will be driven by surveys of large-scale structure. Its inherent non-linearity suggests that significant information will be embedded in higher correlations beyond the two-point function. Extracting this information is extremely challenging: it requires accurate theoretical modelling and significant computational resources to estimate the covariance matrix describing correlations between different Fourier configurations. We investigate whether it is possible to reduce the covariance matrix without significant loss of information by using a proxy that aggregates the bispectrum over a subset of Fourier configurations. Specifically, we study the constraints on $\Lambda$CDM parameters from combining the power spectrum with (a) the modal bispectrum decomposition, (b) the line correlation function and (c) the integrated bispectrum. We forecast the error bars achievable on $\Lambda$CDM parameters using these proxies in a future galaxy survey and compare them to those obtained from measurements of the Fourier bispectrum, including simple estimates of their degradation in the presence of shot noise. Our results demonstrate that the modal bispectrum performs as well as the Fourier bispectrum, even with considerably fewer modes than Fourier configurations. The line correlation function has good performance but does not match the modal bispectrum. The integrated bispectrum is comparatively insensitive to changes in the background cosmology. We find that adding bispectrum data can improve constraints on bias parameters and the normalization $\sigma_8$ by up to 5 compared to power spectrum measurements alone. For other parameters, improvements of up to $\sim$ 20% are possible. Finally, we use a range of theoretical models to explore how the sophistication required for realistic predictions varies with each proxy. (abridged)

In 'modified' gravity the observed acceleration of the universe is explained by changing the gravitational force law or the number of degrees of freedom in the gravitational sector. Both possibilities can be tested by measurements of cosmological structure formation. In this paper we elaborate the details of such tests using the Galileon model as a case study. We pay attention to the possibility that each new degree of freedom may have stochastically independent initial conditions, generating different types of potential well in the early universe and breaking complete correlation between density and velocity power spectra. This 'stochastic bias' can confuse schemes to parametrize the predictions of modified gravity models, such as the use of the growth parameter f alone. Using data from the WiggleZ Dark Energy Survey we show that it will be possible to obtain constraints using information about the cosmological-scale force law embedded in the multipole power spectra of redshift-space distortions. As an example, we obtain an upper limit on the strength of the conformal coupling to matter in the cubic Galileon model, giving |1/M| < 200 / Mp. This allows the fifth-force to be stronger than gravity, but is consistent with zero coupling.

We present a complete framework for numerical calculation of the power spectrum and bispectrum in canonical inflation with an arbitrary number of light or heavy fields. Our method includes all relevant effects at tree-level in the loop expansion, including (i) interference between growing and decaying modes near horizon exit; (ii) correlation and coupling between species near horizon exit and on superhorizon scales; (iii) contributions from mass terms; and (iv) all contributions from coupling to gravity. We track the evolution of each correlation function from the vacuum state through horizon exit and the superhorizon regime, with no need to match quantum and classical parts of the calculation; when integrated, our approach corresponds exactly with the tree-level Schwinger or 'in-in' formulation of quantum field theory. In this paper we give the equations necessary to evolve all two- and three-point correlation functions together with suitable initial conditions. The final formalism is suitable to compute the amplitude, shape, and scale dependence of the bispectrum in models with |fNL| of order unity or less, which are a target for future galaxy surveys such as Euclid, DESI and LSST. As an illustration we apply our framework to a number of examples, obtaining quantitatively accurate predictions for their bispectra for the first time. Two accompanying reports describe publicly-available software packages that implement the method.

CppTransport is a numerical platform that can automatically generate and solve the evolution equations for the 2- and 3-point correlation functions (in field space and for the curvature perturbation) for any inflationary model with canonical kinetic terms. It makes no approximations beyond the applicability of tree-level perturbation theory. Given an input Lagrangian, CppTransport performs symbolic calculations to determine the 'Feynman rules' of the model and generates efficient C++ to integrate the correlation functions of interest. It includes a visualization suite that automates extraction of observable quantities from the raw n-point functions and generates high quality plots with minimal manual intervention. It is intended to be used as a collaborative platform, promoting the rapid investigation of models and systematizing their comparison with observation. This guide describes how to install and use the system, and illustrates its use through some simple examples.

Observations of the microwave background fluctuations suggest a scale-dependent amplitude asymmetry of roughly 2.5 sigma significance. Inflationary explanations for this 'anomaly' require non-Gaussian fluctuations which couple observable modes to those on much larger scales. In this Letter we describe an analysis of such scenarios which significantly extends previous treatments. We identify the non-Gaussian 'response function' which characterizes the asymmetry, and show that it is non-trivial to construct a model which yields a sufficient amplitude: many independent fine tunings are required, often making such models appear less likely than the anomaly they seek to explain. We present an explicit model satisfying observational constraints and determine for the first time how large its bispectrum would appear to a Planck-like experiment. Although this model is merely illustrative, we expect it is a good proxy for the bispectrum in a sizeable class of models which generate a scale-dependent response using a large eta parameter.

If the primordial bispectrum is sufficiently large then the CMB hemispherical asymmetry may be explained by a large-scale mode of exceptional amplitude which perturbs the zeta two-point function. We extend previous calculations, which were restricted to one- or two-source scenarios, by providing a method to compute the response of the two-point function in any model yielding a 'local-like' bispectrum. In general, this shows that it is not the reduced bispectrum fNL which sources the amplitude and scale-dependence of the mode coupling but rather a combination of 'response functions'. We discuss why it is difficult to construct successful scenarios and enumerate the fine-tunings which seem to be required. Finally, we exhibit a concrete model which can be contrived to match the observational constraints and show that to a Planck-like experiment it would appear to have |fNL-local| ~ |fNL-equi| ~ |fNL-ortho| ~ 1. Therefore, contrary to previous analyses, we conclude that it is possible to generate the asymmetry while respecting observational constraints on the bispectrum and low-ell multipoles even without tuning our location on the long-wavelength mode.

We give an explicit relation, up to second-order terms, between scalar-field fluctuations defined on spatially-flat slices and the curvature perturbation on uniform-density slices. This expression is a necessary ingredient for calculating observable quantities at second-order and beyond in multiple-field inflation. We show that traditional cosmological perturbation theory and the `separate universe' approach yield equivalent expressions for superhorizon wavenumbers, and in particular that all nonlocal terms can be eliminated from the perturbation-theory expressions.

We describe how to apply the transport method to compute inflationary observables in a broad range of multiple-field models. The method is efficient and encompasses scenarios with curved field-space metrics, violations of slow-roll conditions and turns of the trajectory in field space. It can be used for an arbitrary mass spectrum, including massive modes and models with quasi-single-field dynamics. In this note we focus on practical issues. It is accompanied by a Mathematica code which can be used to explore suitable models, or as a basis for further development.

In this short paper, we present constraints on the Galileon inflationary model from the CMB bispectrum. We employ a principal-component analysis of the independent degrees of freedom constrained by data and apply this to the WMAP 9-year data to constrain the free parameters of the model. A simple Bayesian comparison establishes that support for the Galileon model from bispectrum data is at best weak.

We develop a needlet approach to estimate the amplitude of general (including non-separable) bispectra and trispectra in the cosmic microwave background, and apply this to the WMAP 9-year data. We obtain estimates for the `orthogonal' bispectrum mode, yielding results which are consistent with the WMAP 7-year data. We do not observe the frequency-dependence suggested by the WMAP team's analysis of the 9-year data. We present 1-$\sigma$ constraints on the `local' trispectrum shape $\gnl/10^5= -4.1\pm 2.3$, the `$c1$' equilateral model $\gnl^{c_1}/10^6= -0.8\pm 2.9$, and the constant model $\gnl^{\rm{const}}/10^6= -0.2\pm 1.8$, together with a $95\%$ confidence-level upper bound on the multifield local parameter $\taunl<22000$. We estimate the bias on these parameters produced by point sources. The techniques developed in this paper should prove useful for other datasets such as Planck.

We use WMAP 9-year bispectrum data to constrain the free parameters of an 'effective field theory' describing fluctuations in single-field inflation. The Lagrangian of the theory contains a finite number of operators associated with unknown mass scales. Each operator produces a fixed bispectrum shape, which we decompose into partial waves in order to construct a likelihood function. Based on this likelihood we are able to constrain four linearly independent combinations of the mass scales. As an example of our framework we specialize our results to the case of 'Dirac-Born-Infeld' and 'ghost' inflation and obtain the posterior probability for each model, which in Bayesian schemes is a useful tool for model comparison. Our results suggest that DBI-like models with two or more free parameters are disfavoured by the data by comparison with single parameter models in the same class.

We derive the 'separate universe' method for the inflationary bispectrum, beginning directly from a field-theory calculation. We work to tree-level in quantum effects but to all orders in the slow-roll expansion, with masses accommodated perturbatively. Our method provides a systematic basis to account for novel sources of time-dependence in inflationary correlation functions, and has immediate applications. First, we use our result to obtain the correct matching prescription between the 'quantum' and 'classical' parts of the separate universe computation. Second, we elaborate on the application of this method in situations where its validity is not clear. As a by-product of our calculation we give the leading slow-roll corrections to the three-point function of field fluctuations on spatially flat hypersurfaces in a canonical, multiple-field model.

We present the first 3+1-dimensional numerical simulations of scalar fields with non-minimal kinetic terms. As an example, we examine the existence and stability of preheating in the presence of a Dirac-Born-Infeld (DBI) inflaton coupled to a canonical matter field. The simulations represent the full non-linear theory in the presence of an expanding Universe. We show that parametric resonance in the matter field, along with self-resonance in the inflaton, repopulate the Universe with matter particles as efficiently as in traditional preheating.

We provide a formula for the scaling behaviour of the inflationary bispectrum in the 'squeezed' limit where one momentum becomes much smaller than the other two. This determines the scaling of the halo bias at low wavenumber and will be an important observable for the next generation of galaxy surveys. Our formula allows it to be predicted for the first time for a generic inflationary model with multiple light, canonically-normalized scalar fields.

In this paper we combine partial-wave (`modal') methods with a wavelet analysis of the CMB bispectrum. Our implementation exploits the advantages of both approaches to produce robust, reliable and efficient estimators which can constrain the amplitude of arbitrary primordial bispectra. This will be particularly important for upcoming surveys such as \emph{Planck}. A key advantage is the computational efficiency of calculating the inverse covariance matrix in wavelet space, producing an error bar which is close to optimal. We verify the efficacy and robustness of the method by applying it to WMAP7 data, finding $\fnllocal=38.4\pm 23.6$ and $\fnlequil=-119.2\pm123.6$.

We use transport techniques to calculate the trispectrum produced in multiple-field inflationary models with canonical kinetic terms. Our method allows the time evolution of the local trispectrum parameters, tauNL and gNL, to be tracked throughout the inflationary phase. We illustrate our approach using examples. We give a simplified method to calculate the superhorizon part of the relation between field fluctuations on spatially flat hypersurfaces and the curvature perturbation on uniform density slices, and obtain its third-order part for the first time. We clarify how the 'backwards' formalism of Yokoyama et al. relates to our analysis and other recent work. We supply explicit formulae which enable each inflationary observable to be computed in any canonical model of interest, using a suitable first-order ODE solver.

A pressing problem in comparing inflationary models with observation is the accurate calculation of correlation functions. One approach is to evolve them using ordinary differential equations ("transport equations"), analogous to the Schwinger-Dyson hierarchy of in-out quantum field theory. We extend this approach to the complete set of momentum space correlation functions. A formal solution can be obtained using raytracing techniques adapted from geometrical optics. We reformulate inflationary perturbation theory in this language, and show that raytracing reproduces the familiar "delta N" Taylor expansion. Our method produces ordinary differential equations which allow the Taylor coefficients to be computed efficiently. We use raytracing methods to express the gauge transformation between field fluctuations and the curvature perturbation, zeta, in geometrical terms. Using these results we give a compact expression for the nonlinear gauge-transform part of fNL in terms of the principal curvatures of uniform energy-density hypersurfaces in field space.

We study inflationary perturbations in multiple-field models, for which zeta typically evolves until all isocurvature modes decay--the "adiabatic limit". We use numerical methods to explore the sensitivity of the nonlinear parameter fNL to the process by which this limit is achieved, finding an appreciable dependence on model-specific data such as the time at which slow-roll breaks down or the timescale of reheating. In models with a sum-separable potential where the isocurvature modes decay before the end of the slow-roll phase we give an analytic criterion for the asymptotic value of fNL to be large. Other examples can be constructed using a waterfall field to terminate inflation while fNL is transiently large, caused by descent from a ridge or convergence into a valley. We show that these two types of evolution are distinguished by the sign of the bispectrum, and give approximate expressions for the peak fNL.

We compute the covariant three-point function near horizon-crossing for a system of slowly-rolling scalar fields during an inflationary epoch, allowing for an arbitrary field-space metric. We show explicitly how to compute its subsequent evolution using a covariantized version of the separate universe or "delta-N" expansion, which must be augmented by terms measuring curvature of the field-space manifold, and give the nonlinear gauge transformation to the comoving curvature perturbation. Nonlinearities induced by the field-space curvature terms are a new and potentially significant source of non-Gaussianity. We show how inflationary models with non-minimal coupling to the spacetime Ricci scalar can be accommodated within this framework. This yields a simple toolkit allowing the bispectrum to be computed in models with non-negligible field-space curvature.

We present a simple and efficient method to compute the superhorizon evolution of the spectral index in multi-field inflationary models, using transport equation techniques. We illustrate the evolution of n_s with time for various interesting potentials.

Galileon fields arise naturally from the decoupling limit of massive gravities, and possess special self-interactions which are protected by a spacetime generalization of Galilean symmetry. We briefly revisit the inflationary phenomenology of Galileon theories. Working from recent computations of the fluctuation Lagrangian to cubic order in the most general model with second-order equations of motion, we show that a distinct shape is present but with suppressed amplitude. A similar shape has been found in other higher-derivative models. It may be visible in a theory tuned to suppress the leading-order shapes, or if the overall bispectrum has large amplitude. Using a partial-wave expansion of the bispectrum, we suggest a possible origin for the frequent appearance of this shape. It follows that models with very disparate microphysics can produce very similar bispectra. We argue that it may be more profitable to distinguish these models by searching for relations between the amplitudes of these common shapes. We illustrate this method using the example of DBI and k-inflation.

We have recently shown that multi-field axion N-flation can lead to observable non-gaussianity in much of its parameter range, with the assisted inflation mechanism ensuring that the density perturbations are sufficiently close to scale invariance. In this paper we extend our analysis in several directions. In the case of equal-mass axions, we compute the probability distributions of observables and their correlations across the parameter space. We examine the case of unequal masses, and show that the mass spectrum must be very densely packed if the model is to remain in agreement with observations. The model makes specific testable predictions for all major perturbative observables, namely the spectral index, tensor-to-scalar ratio, bispectrum, and trispectrum.

Nongaussian statistics are a powerful discriminant between inflationary models, particularly those with noncanonical kinetic terms. Focusing on theories where the Lagrangian is an arbitrary Lorentz-invariant function of a scalar field and its first derivatives, we review and extend the calculation of the observable three-point function. We compute the "next-order" slow-roll corrections to the bispectrum in closed form, and obtain quantitative estimates of their magnitude in DBI and power-law k-inflation. In the DBI case our results enable us to estimate corrections from the shape of the potential and the warp factor: these can be of order several tens of percent. We track the possible sources of large logarithms which can spoil ordinary perturbation theory, and use them to obtain a general formula for the scale dependence of the bispectrum. Our result satisfies the next-order version of Maldacena's consistency condition and an equivalent consistency condition for the scale dependence. We identify a new bispectrum shape available at next-order, which is similar to a shape encountered in Galileon models. If fNL is sufficiently large this shape may be independently detectable.

We study the evolution of non-Gaussianity in multiple-field inflationary models, focusing on three fundamental questions: (a) How is the sign and peak magnitude of the non-linearity parameter fNL related to generic features in the inflationary potential? (b) How sensitive is fNL to the process by which an adiabatic limit is reached, where the curvature perturbation becomes conserved? (c) For a given model, what is the appropriate tool -- analytic or numerical -- to calculate fNL at the adiabatic limit? We summarise recent results obtained by the authors and further elucidate them by considering an inflection point model.

Galileon inflation is a radiatively stable higher derivative model of inflation. The model is determined by a finite number of relevant operators which are protected by a covariant generalization of the Galileon shift symmetry. We show that the nongaussianity of the primordial density perturbation generated during an epoch of Galileon inflation is a particularly powerful observational probe of these models and that, when the speed of sound is small, fNL can be larger than the usual result fNL ~ 1/c_s^2.

We study the transformation properties of a scalar-tensor theory, coupled to fermions, under the Weyl rescaling associated with a transition from the Jordan to the Einstein frame. We give a simple derivation of the corresponding modification to the gauge couplings. After changing frame, this gives rise to a direct coupling between the scalar and the gauge fields.

In a recent publication, we proposed that inflationary perturbation theory can be reformulated in terms of a probability transport equation, whose moments determine the correlation properties of the primordial curvature perturbation. In this paper we generalize this formulation to an arbitrary number of fields. We deduce ordinary differential equations for the evolution of the moments of zeta on superhorizon scales, which can be used to obtain an evolution equation for the dimensionless bispectrum, fNL. Our equations are covariant in field space and allow identification of the source terms responsible for evolution of fNL. In a model with M scalar fields, the number of numerical integrations required to obtain solutions of these equations scales like O(M^3). The performance of the moment transport algorithm means that numerical calculations with M >> 1 fields are straightforward. We illustrate this performance with a numerical calculation of fNL in Nflation models containing M ~ 10^2 fields, finding agreement with existing analytic calculations. We comment briefly on extensions of the method beyond the slow-roll approximation, or to calculate higher order parameters such as gNL.

We study perturbations in the multi-field axion Nflation model, taking account of the full cosine potential. We find significant differences with previous analyses which made a quadratic approximation to the potential. The tensor-to-scalar ratio and the scalar spectral index move to lower values, which nevertheless provide an acceptable fit to observation. Most significantly, we find that the bispectrum non-gaussianity parameter f_NL may be large, typically of order 10 for moderate values of the axion decay constant, increasing to of order 100 for decay constants slightly smaller than the Planck scale. Such a non-gaussian fraction is detectable. We argue that this property is generic in multi-field models of hilltop inflation.

A Galileon field is one which obeys a spacetime generalization of the non-relativistic Galilean invariance. Such a field may possess non-canonical kinetic terms, but ghost-free theories with a well-defined Cauchy problem exist, constructed using a finite number of relevant operators. The interactions of this scalar with matter are hidden by the Vainshtein effect, causing the Galileon to become weakly coupled near heavy sources. We revisit estimates of the fifth force mediated by a Galileon field, and show that the parameters of the model are less constrained by experiment than previously supposed.

We study how the universe reheats following an era of chaotic Dirac-Born-Infeld inflation, and compare the rate of particle production with that in models based on canonical kinetic terms. Particle production occurs through non-perturbative resonances whose structure is modified by the nonlinearities of the Dirac-Born-Infeld action. We investigate these modifications and show that the reheating process may be efficient. We estimate the initial temperature of the subsequent hot, radiation-dominated phase.

In this article, I briefly review the status of infrared effects which occur when using inflationary models to calculate initial conditions for a subsequent hot, dense plasma phase. Three types of divergence have been identified in the literature: secular, "time-dependent" logarithms, which grow with time spent outside the horizon; "box-cutoff" logarithms, which encode a dependence on the infrared cutoff when calculating in a finite-sized box; and "quantum" logarithms, which depend on the ratio of a scale characterizing new physics to the scale of whatever process is under consideration, and whose interpretation is the same as conventional field theory. I review the calculations in which these divergences appear, and discuss the methods which have been developed to deal with them.

In this paper we study various particle physics effects of a light, scalar dark energy field with chameleon-like couplings to matter. We show that a chameleon model with only matter couplings will induce a coupling to photons. In doing so, we derive the first microphysical realization of a chameleonic dark energy model coupled to the electromagnetic field strength. This analysis provides additional motivation for current and near-future tests of axion-like and chameleon particles. We find a new bound on the coupling strength of chameleons in uniformly coupled models. We also study the effect of chameleon fields on Higgs production, which is relevant for hadron colliders. These are expected to manufacture Higgs particles through weak boson fusion, or associated production with a Z or W. We show that, like the Tevatron, the LHC will not be able to rule out or observe chameleons through this mechanism, because gauge invariance of the low energy Lagrangian suppresses the corrections that may arise.

In hybrid inflationary models, inflation ends by a sudden instability associated with a steep ridge in the potential. Here we argue that this feature can generate a large contribution to the curvature perturbation on observable scales. This contribution is almost scale-invariant but highly non-Gaussian. The degree of non-Gaussianity can exceed current observational bounds, unless the inflationary scale is extremely low or the hybrid potential contains very large coupling constants. Non-linear effects on small scales may quench the non-Gaussian signal, and while we find no compelling evidence that this occurs, full lattice simulations are required to definitively address this issue. Note added: We now believe that nonlinear effects will invalidate the original computation in this paper essentially instantaneously after the short-wavelength modes reach the minimum of their potential. This means that the mechanism described in this paper will not lead to appreciatable curvature perturbations on long wavelengths, and no useful constraints on hybrid inflation will result. We have inserted a brief calculation on p2 of this manuscript to explain this fact, but have otherwise left the manuscript unchanged.

An estimate of the one-loop correction to the power spectrum of the primordial curvature perturbation is given, assuming it is generated during a phase of single-field, slow-roll inflation. The loop correction splits into two parts, which can be calculated separately: a purely quantum-mechanical contribution which is generated from the interference among quantized field modes around the time when they cross the horizon, and a classical contribution which comes from integrating the effect of field modes which have already passed far beyond the horizon. The loop correction contains logarithms which may invalidate the use of naive perturbation theory for cosmic microwave background (CMB) predictions when the scale associated with the CMB is exponentially different from the scale at which the fundamental theory which governs inflation is formulated.

The leading quantum correction to the power spectrum of a gravitationally-coupled light scalar field is calculated, assuming that it is generated during a phase of single-field, slow-roll inflation.

We present a novel method for calculating the primordial non-Gaussianity produced by super-horizon evolution during inflation. Our method evolves the distribution of coarse-grained inflationary field values using a transport equation. We present simple evolution equations for the moments of this distribution, such as the variance and skewness. This method possesses some advantages over existing techniques. Among them, it cleanly separates multiple sources of primordial non-Gaussianity, and is computationally efficient when compared with popular alternatives, such as the "delta N" framework. We adduce numerical calculations demonstrating that our new method offers good agreement with those already in the literature. We focus on two fields and the fNL parameter, but we expect our method will generalize to multiple scalar fields and to moments of arbitrarily high order. We present our expressions in a field-space covariant form which we postulate to be valid for any number of fields.

We study models in which a light scalar dark energy particle couples to the gauge fields of the electroweak force, the photon, Z and W bosons. Our analysis applies to a large class of interacting dark energy models, including those in which the dark energy mass can be adjusted to evade fifth-force bounds by the so-called chameleon mechanism. We conclude that--with the usual choice of Higgs sector--electroweak precision observables are screened from the indirect effects of dark energy, making such corrections effectively unobservable at present-day colliders, and limiting the dark energy discovery potential of any future International Linear Collider. We show that a similar screening effect applies to processes mediated by flavour-changing neutral currents, which can be traced to the Glashow-Iliopoulos-Maiani mechanism. However, Higgs boson production at the Large Hadron Collider via weak boson fusion may receive observable corrections.

The primordial density fluctuation inevitably couples to all forms of matter via loop corrections and depends on the ambient conditions while inflation was ongoing. This gives us an opportunity to observe processes which were in progress while the universe was inflating, provided they were sufficiently dramatic to overcome suppression by powers of (H/MP)^2 ~ 10^(-9), where H is the Hubble scale during inflation and MP is the Planck mass. As an example, if a primordial magnetic field was synthesized during inflation, as suggested by some interpretations of the apparently universal 10^(-6) gauss field observed on galactic scales, then this could leave traces in inflationary observables. In this paper, I compute corrections to the spectrum and bispectrum generated by a varying electromagnetic coupling during inflation, assuming that the variation in this coupling is mediated by interaction with a collection of light scalar fields. If the mass scale associated with this interaction is too far below the Planck scale then the stability of perturbation theory can be upset. For the mass-scale which is relevant in the standard magnetogenesis scenario, however, the theory is stable and the model is apparently consistent with observational constraints.

We study signatures in the Cosmic Microwave Background (CMB) induced by the presence of strong spatial curvature prior to the epoch of inflation which generated our present universe. If inflation does not last sufficiently long to drive the large-scale spatial curvature to zero, then presently observable scales may have left the horizon while spatial slices could not be approximated by a flat, Euclidean geometry. We compute corrections to the power spectrum and non-gaussianity of the CMB temperature anisotropy in this scenario. The power spectrum does not receive significant corrections and is a weak diagnostic of the presence of curvature in the initial conditions, unless its running can be determined with high accuracy. However, the bispectral non-gaussianity parameter f_NL receives modifications on the largest observable scales. We estimate that the maximum signal would correspond to f_NL ~ 0.3, which is out of reach for present-day microwave background experiments.

It is well-known that expectation values in de Sitter space are afflicted by infra-red divergences. Long ago, Starobinsky proposed that infra-red effects in de Sitter space could be accommodated by evolving the long-wavelength part of the field according to the classical field equations plus a stochastic source term. I argue that--when quantum-mechanical loop corrections are taken into account--the separate-universe picture of superhorizon evolution in de Sitter space is equivalent, in a certain leading-logarithm approximation, to Starobinsky's stochastic approach. In particular, the time evolution of a box of de Sitter space can be understood in exact analogy with the DGLAP evolution of partons within a hadron, which describes a slow logarithmic evolution in the distribution of the hadron's constituent partons with the energy scale at which they are probed.

We compute the connected four-point correlation function of the primordial curvature perturbation generated during inflation with standard kinetic terms, where the correlation is established via exchange of a graviton between two pairs of scalar fluctuations. Any such correlation yields a contribution to the scalar trispectrum of the order of the tensor to scalar ratio r. This contribution is numerically one order of magnitude larger than the one previously calculated on the basis of scalar perturbations interacting at a point and satisfies a simple relation in the limit where the momentum of the graviton which is exchanged becomes much smaller than the external momenta. We conclude that the total non-linearity parameter generated by single-field models of slow-roll inflation is at maximum tauNL ~ r.

A new and powerful probe of the origin and evolution of structures in the Universe has emerged and been actively developed over the last decade. In the coming decade, non-Gaussianity, i.e., the study of non-Gaussian contributions to the correlations of cosmological fluctuations, will become an important probe of both the early and the late Universe. Specifically, it will play a leading role in furthering our understanding of two fundamental aspects of cosmology and astrophysics: (i) the physics of the very early universe that created the primordial seeds for large-scale structures, and (ii) the subsequent growth of structures via gravitational instability and gas physics at later times. To date, observations of fluctuations in the Cosmic Microwave Background (CMB) and the Large-Scale Structure of the Universe (LSS) have focused largely on the Gaussian contribution as measured by the two-point correlations (or the power spectrum) of density fluctuations. However, an even greater amount of information is contained in non-Gaussianity and a large discovery space therefore still remains to be explored. Many observational probes can be used to measure non-Gaussianity, including CMB, LSS, gravitational lensing, Lyman-alpha forest, 21-cm fluctuations, and the abundance of rare objects such as clusters of galaxies and high-redshift galaxies. Not only does the study of non-Gaussianity maximize the science return from a plethora of present and future cosmological experiments and observations, but it also carries great potential for important discoveries in the coming decade.

We calculate the tree-level bispectrum of the inflaton field perturbation directly from the field equations, and construct the corresponding f_NL parameter. Our results agree with previous ones derived from the Lagrangian. We argue that quantum theory should only be used to calculate the correlators when they first become classical a few Hubble times after horizon exit, the classical evolution taking over thereafter.

We show that during inflation, a quantum fluctuation becomes classical at all orders if it becomes classical at first order. Implications are discussed.

We derive the Klein--Gordon equation for a single scalar field coupled to gravity at second order in perturbation theory and leading order in slow-roll. This is done in two ways: we derive the Klein--Gordon equation first using the Einstein field equations, and then directly from the action after integrating out the constraint equations. We also point out an unexpected result regarding the treatment of the field equations.

Inflationary observables, like the power spectrum, computed at one- and higher-order loop level seem to be plagued by large infra-red corrections. In this short note, we point out that these large infra-red corrections appear only in quantities which are not directly observable. This is in agreement with general expectations concerning infra-red effects.

We study brane inflation scenarios in a warped throat geometry and show that there exists a consistency condition between the non-Gaussianity of the curvature perturbation and the amplitude and scale-dependence of the primordial gravitational waves. This condition is independent of the warping of the throat and the form of the inflaton potential. We find that such a relation could be tested by a future CMB polarization experiment if the Planck satellite is able to detect both a gravitational wave background and a non-Gaussian statistic. In models where the observable stage of inflation occurs when the brane is in the tip region of the throat, we derive a further consistency condition involving the scalar spectral index, the tensor-scalar ratio and the curvature perturbation bispectrum. We show that when such a relation is combined with the WMAP3 results, it leads to a model-independent bound on the gravitational wave amplitude given by 0.001 < r < 0.01. This corresponds to the range of sensitivity of the next generation of CMB polarization experiments.

We calculate the trispectrum of the primordial curvature perturbation generated by an epoch of slow-roll inflation in the early universe, and demonstrate that the non-gaussian signature imprinted at horizon crossing is unobservably small, of order tau_NL < r/50, where r < 1 is the tensor-to-scalar ratio. Therefore any primordial non-gaussianity observed in future microwave background experiments is likely to have been synthesized by gravitational effects on superhorizon scales. We discuss the application of Maldacena's consistency condition to the trispectrum.

We present an estimate for the non-linear parameter \tau_NL, which measures the non-gaussianity imprinted in the trispectrum of the comoving curvature perturbation, \zeta. Our estimate is valid throughout the inflationary era, until the slow-roll approximation breaks down, and takes into account the evolution of perturbations on superhorizon scales. We find that the non-gaussianity is always small if the field values at the end of inflation are negligible when compared to their values at horizon crossing. Under the same assumption, we show that in Nflation-type scenarios, where the potential is a sum of monomials, the non-gaussianity measured by \tau_NL is independent of the couplings and initial conditions.

We show how to obtain the probability density function for the amplitude of the curvature perturbation, R, produced during an epoch of slow-roll, single-field inflation, working directly from n-point correlation functions of R. These n-point functions are the usual output of quantum field theory calculations, and as a result we bypass approximate statistical arguments based on the central limit theorem. Our method can be extended to deal with arbitrary forms of non-Gaussianity, appearing at any order in the n-point hierarchy. We compute the probability density for the total smoothed perturbation within a Hubble volume, \epsilon, and for the spectrum of \epsilon. When only the two-point function is retained, exact Gaussian statistics are recovered. When the three-point function is taken into account, we compute explicitly the leading slow-roll correction to the Gaussian result.

We use the dS/CFT correspondence and bulk gravity to predict the form of the renormalized holographic three-point correlation function of the operator which is dual to the inflaton field perturbation during single-field, slow-roll inflation. Using Maldcaena's formulation of the correspondence, this correlator can be related to the three-point function of the curvature perturbation generated during single-field inflation, and we find exact agreement with previous bulk QFT calculations. This provides a consistency check on existing derivations of the non-Gaussianity from single-field inflation and also yields insight into the nature of the dS/CFT correspondence. As a result of our calculation, we obtain the properly renormalized dS/CFT one-point function, including boundary contributions where derivative interactions are present in the bulk. In principle, our method may be employed to derive the n-point correlators of the inflationary curvature perturbation within the context of (n-1)th-order perturbation theory, rather than nth-order theory as in conventional approaches.

Observational consequences of inflationary cosmology in the holographic dual of the Randall-Sundrum type II braneworld scenario, as motivated by the AdS/CFT correspondence, are investigated. High energy corrections to the standard four-dimensional Friedmann equation induce a corresponding modification to the form of the single-field inflationary consistency equation based on Einstein gravity. The degree of departure from the standard expression is determined by the ratio, r, of the primordial tensor and scalar perturbation amplitudes and the coefficient, c, of the conformal anomaly in the dual gauge theory. It is found that a necessary condition for detecting such a correction with the next generation of cosmic microwave background (CMB) polarization experiments is that r >= 0.06. The bound is tightened to r > 0.3 for values of the central charge that are compatible with known compactifications of type IIB string theory as parametrized in terms of F-theory compactification on Calabi-Yau four-folds. This is close to the present upper bounds inferred from combined observations of the CMB anisotropy power spectrum and high redshift surveys. We conclude that if such modifications to the inflationary consistency equation are to be observable, the gravitational wave background should be detected in the near future. A further consequence of the non-standard dynamics at high energies is that the initial state of the universe is a quiescent singularity with a finite density and pressure.

Braneworld cosmology supported by a bulk scalar field with an exponential potential is developed. A general class of separable backgrounds for both single and two-brane systems is derived, where the bulk metric components are given by products of world-volume and bulk coordinates and the world-volumes represent any anisotropic and inhomogeneous solution to an effective four-dimensional Brans-Dicke theory of gravity. We deduce a cosmic no hair theorem for all ever expanding, spatially homogeneous Bianchi world-volumes and find that the spatially flat and isotropic inflationary scaling solution represents a late-time attractor when the bulk potential is sufficiently flat. The dependence of this result on the separable nature of the bulk metric is investigated by applying the techniques of Hamilton-Jacobi theory to five-dimensional Einstein gravity. We employ the spatial gradient expansion method to determine the asymptotic form of the bulk metric up to third-order in spatial gradients. It is found that the condition for the separable form of the metric to represent the attractor of the system is precisely the same as that for the four-dimensional world-volume to isotropize. We also derive the fourth-order contribution to the Hamilton-Jacobi generating functional. Finally, we conclude by placing our results within the context of the holographic approach to braneworld cosmology.

We calculate the three-point correlation function evaluated at horizon crossing for a set of interacting scalar fields coupled to gravity during inflation. This provides the initial condition for the three-point function of the curvature perturbation in the Sasaki--Stewart \delta N formulation. We find that the effect is small, of the order of a slow-roll parameter, and that the non-gaussianity can be determined on large scales once the unperturbed background evolution is known. As an example of the use of our formalism, we calculate the primordial non-gaussianity arising in a model of assisted inflation.

We calculate the three-point function for primordial scalar fluctuations in a single field inflationary scenario where the scalar field Lagrangian is a completely general function of the field and its first derivative. We obtain an explicit expression for the three-point correlation function in a self-consistent approximation scheme where the expansion rate varies slowly, analogous to the slow-roll limit in standard, single-field inflation. The three-point function can be written in terms of the familiar slow-roll parameters and three new parameters which measure the non-trivial kinetic structure of the scalar field, the departure of the sound speed from the speed of light, and the rate of change of the sound speed.

We investigate the constraints on quintessence arising from both renormalisable and non-renormalisable couplings where the 5d Planck mass is around the TeV scale. The quintessence field vacuum expectation value is typically of order the 4d Planck mass while non-renormalisable operators are expected to be suppressed by the 5d Planck mass. Non-renormalisable operators are therefore important in computing the 4d effective quintessence potential. We then study the quantum corrections to the quintessence potential due to fermion and graviton loops. The tower of Kaluza-Klein modes competes with the TeV-scale cut-off, altering the graviton contribution to the vacuum polarization of quintessence. Nevertheless we show that, as in four dimensions, the classical potential is stable to such radiative corrections.

The braneworld cosmology, in which our universe is imbedded as a hypersurface in a higher dimensional bulk, has the peculiar property that the inflationary consistency relation derived in a four-dimensional cosmology persists. This consistency condition relates the ratio of tensor and scalar perturbation amplitudes to the tensor spectral index produced during an epoch of slow-roll scalar field inflation. We attempt to clarify this surprising degeneracy. Our argument involves calculating the power spectrum of scalar field fluctuations around geometries perturbed away from the exact de Sitter case. This calculation is expected to be valid for perturbations which would not cause a late-time acceleration of the universe. We use these results to argue that the emergence of the same consistency relation in the braneworld can be connected with a specific property, that five-dimensional observables smoothly approach their four-dimensional counterparts as one takes the brane to infinite tension. We exhibit an explicit example where this does not occur, and in which a consistency relation does not persist.