Search

In inflationary models that produce a spike of power on short scales, back-reaction of small-scale substructure onto large-scale modes is enhanced. We argue that the separate universe framework provides a highly convenient tool to compute loop corrections that quantify this back-reaction. Each loop of interest is characterized by large hierarchies in wavenumber and horizon exit time. The separate universe framework highlights important factorizations involving these hierarchies. We interpret each loop correction in terms of a simple, classical, back-reaction model, and clarify the meaning of the different volume scalings that have been reported in the literature. We argue that significant back-reaction requires both short-scale nonlinearities and long-short couplings that modulate the short-scale power spectrum. In the absence of long-short couplings, only incoherent shot noise-like effects are present, which are volume-suppressed. Dropping the shot noise, back-reaction from a particular scale is controlled by a product of $f_{NL}$-like parameters: an equilateral configuration measuring the nonlinearity of the short-scale modes, and a squeezed configuration measuring the long-short coupling. These may carry important scale dependence controlling the behaviour of the loop in the decoupling limit where the hierarchy of scales becomes large. In single-field models the long-short coupling may be suppressed by this hierarchy, in which case the net back-reaction would be safely suppressed. We illustrate our framework using explicit computations in a 3-phase ultra-slow-roll scenario. Finally, we discuss different choices for the smoothing scale used in the separate universe framework and argue the effect can be absorbed into a renormalization of local operators., Comment: Missing terms in the 13-type loop have been added. The conclusions are strengthened, and the general message is unchanged

We update predictions for observables in the "delicate" D3/anti-D3 inflationary model on the conifold. We use a full CMB likelihood calculation to assess goodness-of-fit, which is necessary because in this model the zeta power spectrum often cannot be modelled as a power-law over observable scales. For the first time we are able to provide accurate forecasts for the amplitude of three-point correlations. In a significant portion of its parameter space the model follows Maldacena's single-field prediction fNL ~ -(5/12)(ns-1) if nt << 1. Therefore |fNL| is usually small when the power spectrum satisfies observational constraints. In a small number of cases the bispectrum is instead dominated by effects from rapid switching between angular minima. The resulting amplitudes are larger, but mostly with unacceptable spectral behaviour. In the most extreme case we obtain |fNLeq| ~ 75 at kt/3 = 0.002/Mpc. It has been suggested that the quasi-single field inflation ("QSFI") mechanism could produce significant 3-point correlations in this model. We do observe rare shifts in amplitude between equilateral and squeezed configurations that could possibly be associated with QSFI effects, but more investigation is needed to establish the full bispectrum shape. There is evidence of "shape" running between equilateral and squeezed configurations that may be inherited from the scale dependence of the spectrum. We explore the dependence of observables on discrete choices such as the truncation point of the potential. Our analysis illustrates the advantages of a standard format for information exchange within the inflationary model-building and testing community., Comment: 65 pages plus appendices. See Code & Data tab for data deposit. v2: Minor changes. This is an author-created, un-copyedited version of an article published in JCAP. IOP Publishing Ltd is not responsible for any errors or omissions in this version of the manuscript or any version derived from it. The Version of Record is available online at 10.1088/1475-7516/2022/02/013

We consider a disformal coupling between Standard Model matter and a cubic Galileon scalar sector, assumed to be a relict of some other physics that solves the cosmological constant problem rather than a solution in its own right. This allows the energy density carried by the Galileon scalar to be sufficiently small that it evades stringent constraints from the integrated Sachs-Wolfe effect, which otherwise rules out the cubic Galileon theory. Although the model with disformal coupling does not exhibit screening, we show there is a `screening-like' phenomenon in which the energy density carried by the Galileon scalar is suppressed during matter domination when the quadratic and cubic Galileon operators are both relevant and the quadratic sector has a stable kinetic term. We obtain the explicit 3+1 form of Maxwell's equations in the presence of the disformal coupling, and the wave equations that govern electromagnetic waves. The disformal coupling is known to generate a small mass that modifies their velocity of propagation. We use the WKB approximation to study electromagnetic waves in this theory and show that, despite remarkable recent constraints from the LIGO/Virgo observatories that restrict the difference in propagation velocity between electromagnetic and gravitational radiation to roughly 1 part in $10^{15}$, the disformal coupling is too weak to be constrained by events such as GW170817 or by the dispersion of electromagnetic radiation at different wavelengths., Comment: 18 pages, 1 figure. v2: Stability analysis added. This is an author-created, un-copyedited version of an article published in JCAP. IOP Publishing Ltd is not responsible for any errors or omissions in this version of the manuscript or any version derived from it. The Version of Record is available online at 10.1088/1475-7516/2021/10/085

We revisit the "sequester" proposal of Kaloper, Padilla and collaborators, in which the amplitude of the cosmological constant is decoupled from large contributions due to loops containing Standard Model particles. We review the different formulations of the model that have appeared in the literature, and estimate the importance of a particular class of quantum corrections-those that dress the interaction between the "rigid" scalars and infrared properties of the spacetime such as its 4-volume and integrated curvature. In formulations that do not adequately sequester graviton loops we argue that dressing of these interactions causes further failures of complete sequestration. We estimate the size of the effect and find that it is typically smaller than the cosmological term directly induced by loops containing a single virtual graviton. Meanwhile, in the most developed formulation of the scenario (where a rigid scalar couples to the Gauss-Bonnet density), this dressing can be absorbed into a rescaling of the rigid fields and is therefore harmless., Comment: 11 pages, 2 figures. v2: Minor changes to notation, accepted in PRD

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}$., Comment: 34 pages, 8 figures. Appendices give further details of the one-loop SPT calculations

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., Comment: 25 pages, 11 figures + 17 page appendix and references

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), Comment: 38 pages, 16 figures. Supporting dataset made available at https://doi.org/10.5281/ZENODO.438187

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., Comment: 55 pages excluding appendices and references, 13 figures. v2: minor clarifications. This is an author-created, un-copyedited version of an article published in JCAP. IOP Publishing Ltd is not responsible for any errors or omissions in this version of the manuscript or any version derived from it. The Version of Record is available online at https://doi.org/10.1088/1475-7516/2017/11/039

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., Comment: 128 pages. Matches 2013.3 release version deposited at zenodo.org (see linked DOI) up to changes in typesetting and correction of some minor typos

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., Comment: 59 pages, 16 figures. v2: Typos corrected in Eqs. (6.7c,d) but otherwise no changes. This is an author-created, un-copyedited version of an article published in JCAP. IOP Publishing Ltd is not responsible for any errors or omissions in this version of the manuscript or any version derived from it. The Version of Record is available online at http://dx.doi.org/10.1088/1475-7516/2016/12/033

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., Comment: 5 pages. v2: Minor changes to match version published in Phys. Rev. D

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., Comment: 38 pages plus references, 7 figures. v2: Minor changes to match version published in JCAP. This is an author-created, un-copyedited version of an article published in JCAP. IOP Publishing Ltd is not responsible for any errors or omissions in this version of the manuscript or any version derived from it. The Version of Record is available online at http://dx.doi.org/10.1088/1475-7516/2016/06/025

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., Comment: 17 pages plus references; v2: conclusions unchanged; mainly involves new numerics avoiding obscure SciPy bug that occurred in v1, and updates Alcock-Paczynski methodology

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., Comment: 22 pages; accompanying code available from transportmethod.com. v3: improved treatment of tensor power spectrum and tensor spectral index

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., Comment: 11 pages. Matches published version

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., Comment: 16 pages, no figures. This is an author-created, un-copyedited version of an article published in JCAP. IOP Publishing Ltd is not responsible for any errors or omissions in this version of the manuscript or any version derived from it. The Version of Record is available online at http://dx.doi.org/10.1088/1475-7516/2015/02/040. v3: Updated to match version published in JCAP

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., Comment: 25 pages, plus references and appendix

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., Comment: 21 pages - matches published version

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., Comment: 5 pages, 3 figures; note correction to parameter value in first paragraph of Numerical Simulations section (typographical error in published version)

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., Comment: 5 pages, 1 figure

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$., Comment: 19 pages, 6 figures

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., Comment: v1: 33 pages, plus appendix and references; 5 figures. v2: typographical typos fixed, minor changes to the main text and abstract, reference added; matches version published in JCAP

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., Comment: 22 pages, plus appendix and references

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., Comment: 24 pages, plus references and appendix. v2: matches version published in JCAP; typo fixed in Eq. (54)

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., Comment: 22 pages, plus bibliography and appendix. v2: minor changes, matches version published in JCAP

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., Comment: 6 pages; minor changes to match published version

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., Comment: v1: 25 pages, including tables, an appendix and references. v2: minor clarifications about the lowest-order consistency relations; matches version published in JCAP

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., Comment: 9 pages, 12 fgures

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., Comment: 12 pages, 2 figures, to appear in International Journal of Modern Physics A

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., Comment: v1: 25 pages, plus Appendix and bibliography, 6 figures. v2: minor edits to match published version in JCAP

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., Comment: v1: 37 pages, plus tables, figures and appendices. v2: supersedes version published in JCAP; some clarifications and more detailed comparison with earlier literature. All results unchanged. v3:improvements to some plots; text unchanged

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., Comment: 5 pages. v2: minor changes, matches version published in PLB

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., Comment: 13 pages, plus appendices and references. v2: some clarifications in introduction; missing references added. v3: discussion of Sec. IV.B clarified. Matches version published in JCAP. v4: typos corrected in (44)-(47). v5: some improvements to IV.B; sign error fixed in h4 term of (B25) and related numerical error in (47)

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., Comment: 23 pages, plus appendices and references; 4 figures. v2: incorrect statements regarding numerical delta N removed from Sec. 4.3. Minor modifications elsewhere

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., Comment: 5 pages with 3 figures incorporated. v2: minor updates, new analytic formulae for spectral index and trispectrum, references added, typo corrections. v3: minor updates to match version accepted by Physical Review Letters

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., Comment: 23 pages

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., Comment: 14 pages

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., Comment: Invited review for focus section of Classical & Quantum Gravity on nonlinear and nongaussian perturbation theory. Some improvements compared to version which will appear in CQG, especially in Sec. 2.3. 30 pages + references.

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., Comment: 26pp, 2 figs v2: note added (see abstract)

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., Comment: 24 pages, 4 figures. Uses iopart.cls and feynmf.sty. v2: Text revised, but conclusions unchanged. v3: Replaced with version accepted by PRD.

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., Comment: 24 pages, 4 colour figures; uses iopart.cls. v2: Erroneous statements about delta N method in Sec. 2 removed. Correction to gauge transformation in Eq. (12) brings numerical results of Sec. 4 into better agreement with the delta N formula. Conclusions remain the same. v3: minor changes

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., Comment: 23 pages, uses ioplatex.sty. v2: only bibliographic changes

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., Comment: 39 pages, 8 figures, uses ioplatex.sty and feynmp. v2: Matches version accepted by JHEP; extended discussion of straight corrections, summarized in the text and described in detail in a new Appendix A. Conclusions unchanged

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., Comment: 36 pages; uses iopart.cls and feynmp.sty. v2: Minor typos corrected. Matches version published in JCAP

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., Comment: 23 pages, 2 figures. Uses iopart.sty and feynmp.sty. v2: Bibliography updated, no change to text. v3: Updated to match version published in JCAP, with one extra sign error corrected in the Appendix

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., Comment: 37 pages, uses feynmp.sty and iopart.cls. v2: minor improvements in comparison with version submitted to (and accepted by) JCAP. Improves v1 with a more refined discussion of cutoffs and Lorentz invariance in Sections 4-5, but calculations are unchanged. Minor textual improvements throughout

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., Comment: 16 pages, uses iopart.sty. v2: replaced with version accepted by JCAP; minor changes of wording only. v3: supersedes version published by journal; typo fixed in Eq. (20) and updated references. v4: sign errors in Eqs. (32) and (38) corrected

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., Comment: 9 pages, revtex4

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., Comment: 33 pages, uses feynmp.sty and ioplatex journal style. v2: matches version published in JCAP. v3: corrects sign error in Eq. (58). Corrects final coefficient of the logarithm in Eq. (105). Small corrections to discussion of divergences in 1-point function. Minor improvements to discussion of UV behaviour in Sec. 4.2

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., Comment: 28 pages, uses feynmp.sty and ioplatex journal style. v2: supersedes version published in JCAP. Some corrections and refinements to the discussion and conclusions. v3: Corrects misidentification of quantum correction with an IR effect. Improvements to the discussion