We present a method to introduce relativistic corrections including linear dark energy perturbations in Horndeski theory into Newtonian simulations based on the N-body gauge approach. We assume that standard matter species (cold dark matter, baryons, photons and neutrinos) are only gravitationally-coupled with the scalar field and we then use the fact that one can include modified gravity effects as an effective dark energy fluid in the total energy-momentum tensor. In order to compute the scalar field perturbations, as well as the cosmological background and metric perturbations, we use the Einstein-Boltzmann code \hiclass. As an example, we study the impact of relativistic corrections on the matter power spectrum in k-essence, a subclass of Horndeski theory, including the effects of massless and massive neutrinos. For massive neutrinos with $\sum m_{\nu} = 0.1$ eV, the corrections due to relativistic species (photons, neutrinos and dark energy) can introduce a maximum deviation of approximately $8\%$ to the power spectrum at $k \sim 10^{-3} \ \textrm{Mpc}^{-1}$ at $z=0$, for a scalar field with sound speed $c_{s}^{2}\sim0.013$ during matter domination epoch. Our formalism makes it possible to test beyond \lcdm \ models probed by upcoming large-scale structure surveys on very large scales.

We examine the growth of structure in three different cosmological models with interacting dark matter and vacuum energy. We consider the case of geodesic dark matter with zero sound speed, where the relativistic growing mode in comoving-synchronous gauge coincides with the Newtonian growing mode at first order in $\Lambda$CDM. We study corrections to the linearly growing mode in the presence of interactions and the linear matter growth rate, $f_1$, contrasting this with the velocity divergence, $f_{rsd}\sigma_8$, observed through redshift-space distortions. We then derive second-order density perturbations in these interacting models. We identify the reduced bispectrum that corresponds to the non-linear growth of structure and show how the shape of the bispectrum is altered by energy transfer to or from the vacuum. Thus the bispectrum, or higher-order correlators, might in future be used to identify dark matter interactions.

We investigate cosmological models described by a scalar field with an exponential potential, and apply the stochastic formalism, which allows us to study how quantum field fluctuations give rise to stochastic noise. This modifies the classical dynamics of the scalar field at large scales, above a coarse-graining scale. In particular we explore how quantum field fluctuations perturb the equation of state on large scales which can lead to a quantum instability of the classical collapse solution below the Planck scale in the case of a pressureless collapse.

We study the stochastic formalism of inflation beyond the usual slow-roll approximation. We verify that the assumptions on which the stochastic formalism relies still hold even far from the slow-roll attractor. This includes demonstrating the validity of the separate universe approach to evolving long-wavelength scalar field perturbations beyond slow roll. We also explain that, in general, there is a gauge correction to the amplitude of the stochastic noise. This is because the amplitude is usually calculated in the spatially-flat gauge, while the number of e-folds is used as the time variable (hence one works in the uniform-$N$ gauge) in the Langevin equations. We show that these corrections vanish in the slow-roll limit, but we also explain how to calculate them in general. We compute them in difference cases, including ultra-slow roll and the Starobinsky model that interpolates between slow roll and ultra-slow roll, and find the corrections to be negligible in practice. This confirms the validity of the stochastic formalism for studying quantum backreaction effects in the very early universe beyond slow roll.

With an energy scale that can be as high as $10^{14}\,{\rm GeV}$, inflation may provide a unique probe of high-energy physics. Both scalar and tensor fluctuations generated during this early accelerated expansion contain crucial information about the particle content of the primordial universe. The advent of ground- and space-based interferometers enables us to probe primordial physics at length-scales much smaller than those corresponding to current CMB constraints. One key prediction of single-field slow-roll inflation is a red-tilted gravitational wave spectrum, currently inaccessible at interferometer scales. Interferometers probe directly inflationary physics that deviates from the minimal scenario and in particular additional particle content with sizeable couplings to the inflaton field. We adopt here an effective description for such fields and focus on the case of extra spin-2 fields. We find that a time-dependent sound speed for the helicity-2 modes can generate primordial gravitational waves with a blue-tilted spectrum, potentially detectable at interferometer scales.

We investigate an interacting dark sector scenario in which the vacuum energy is free to interact with cold dark matter (CDM), which itself is assumed to cluster under the sole action of gravity, i.e. it is in free fall (geodesic), as in $\Lambda$CDM. The interaction is characterised by a dimensionless coupling $q_{\rm V}$ that we constrain using cosmic microwave background data from the Planck 2015 data release, along with baryon acoustic oscillation, redshift space distortion and Type Ia supernova measurements. We present the full linear perturbation theory of this interacting scenario and use MCMC sampling to study five different cases: two cases in which we have $\Lambda$CDM evolution in the distant past, until a set redshift $z_{\rm trans}$, below which the interaction switches on and $q_{\rm V}$ is the single sampled parameter, with $z_{\rm trans}$ fixed at $z_{\rm trans}=3000$ and $z_{\rm trans}=0.9$ respectively; a case where we allow this transition redshift to vary along with $q_{\rm V}$; a case in which the vacuum energy is zero for $z>z_{\rm trans}$ and then begins to grow once the interaction switches on; and the final case in which we bin $q_{\rm V}(z)$ in four redshift bins to investigate the possibility of a dynamical interaction, reconstructing the redshift evolution of the function using Gaussian processes. We find that, in all cases where the high redshift evolution is not modified, the results are compatible with a vanishing coupling, thus finding no significant deviation from $\Lambda$CDM.

Our current understanding of the Universe is established through the pristine measurements of structure in the cosmic microwave background (CMB) and the distribution and shapes of galaxies tracing the large scale structure (LSS) of the Universe. One key ingredient that underlies cosmological observables is that the field that sources the observed structure is assumed to be initially Gaussian with high precision. Nevertheless, a minimal deviation from Gaussianityis perhaps the most robust theoretical prediction of models that explain the observed Universe; itis necessarily present even in the simplest scenarios. In addition, most inflationary models produce far higher levels of non-Gaussianity. Since non-Gaussianity directly probes the dynamics in the early Universe, a detection would present a monumental discovery in cosmology, providing clues about physics at energy scales as high as the GUT scale.

Tensor non-Gaussianity represents an important future probe of the physics of inflation. Inspired by recent works, we elaborate further on the possibility of significant primordial tensor non-Gaussianities sourced by extra fields during inflation. Unitarity constraints limit the impact of extra (spinning) particle content by means of a lower bound on the corresponding mass spectrum. For spin-2 particles, this takes the form of the well-known Higuchi bound. Massive ($m\gtrsim H$) particles will typically decay during inflation unless they are non-minimally coupled to the inflaton sector: the inflating field "lifts" the dynamics of the extra field(s), effectively getting around the limits imposed by unitarity. There exist several models that realize such a mechanism, but we focus here on the set-up of [1] where, through an EFT approach, one is able to capture the essential features common to an entire class of theories. In the presence of an extra massive spin-2 particle, the interactions in the tensor sector mimic very closely those in the scalar sector of quasi-single-field inflationary models. We calculate the tensor bispectrum in different configurations and extract its dependence on the extra tensor sound speed. We show in detail how one may obtain significant tensor non-Gaussianities whose shape-function interpolates between local and equilateral, depending on the mass of the extra field. We also estimate the LISA response functions to a tensor bispectrum supporting the intermediate-type shapes we find.

The relative probability to decay towards different vacua during inflation is studied. The calculation is performed in single-field slow-roll potentials using the stochastic inflation formalism. Various situations are investigated, including falling from a local maximum of the potential and escaping from a local minimum. In the latter case, our result is consistent with that of Hawking and Moss, but is applicable to any potential. The decay rates are also computed, and the case of a generic potential with multiple minima and maxima is discussed.

Simple models of single-field inflation in the very early universe can generate the observed amplitude and scale dependence of the primordial density perturbation, but models with multiple fields can provide an equally good fit to current data. We show how future observations will be able to distinguish between currently favoured models. If a curvaton field generates the primordial perturbations after inflation, we show how the total duration of inflation can be measured.

It is often claimed that the ultra-slow-roll regime of inflation, where the dynamics of the inflaton field are friction dominated, is a non-attractor and/or transient. In this work we carry out a phase-space analysis of ultra-slow roll in an arbitrary potential, $V(\phi)$. We show that while standard slow roll is always a dynamical attractor whenever it is a self-consistent approximation, ultra-slow roll is stable for an inflaton field rolling down a convex potential with $M_{\scriptscriptstyle{\mathrm{Pl}}} V''>|V'|$ (or for a field rolling up a concave potential with $M_{\scriptscriptstyle{\mathrm{Pl}}} V''<-|V'|$). In particular, when approaching a flat inflection point, ultra-slow roll is always stable and a large number of $e$-folds may be realised in this regime. However, in ultra-slow roll, $\dot{\phi}$ is not a unique function of $\phi$ as it is in slow roll and dependence on initial conditions is retained. We confirm our analytical results with numerical examples.

We study the scalar-tensor-tensor non-Gaussian signal in an inflationary model comprising also an axion coupled with SU(2) gauge fields. In this set-up, metric fluctuations are sourced by the gauge fields already at the linear level providing an enhanced chiral gravitational waves spectrum. The same mechanism is at work in generating an amplitude for the three-point function that is parametrically larger than in standard single-field inflation.

How much more will we learn about single-field inflationary models in the future? We address this question in the context of Bayesian design and information theory. We develop a novel method to compute the expected utility of deciding between models and apply it to a set of futuristic measurements. This necessarily requires one to evaluate the Bayesian evidence many thousands of times over, which is numerically challenging. We show how this can be done using a number of simplifying assumptions and discuss their validity. We also modify the form of the expected utility, as previously introduced in the literature in different contexts, in order to partition each possible future into either the rejection of models at the level of the maximum likelihood or the decision between models using Bayesian model comparison. We then quantify the ability of future experiments to constrain the reheating temperature and the scalar running. Our approach allows us to discuss possible strategies for maximising information from future cosmological surveys. In particular, our conclusions suggest that, in the context of inflationary model selection, a decrease in the measurement uncertainty of the scalar spectral index would be more decisive than a decrease in the uncertainty in the tensor-to-scalar ratio. We have incorporated our approach into a publicly available python class, foxi (https://sites.google.com/view/foxicode), that can be readily applied to any survey optimisation problem.

We show that in the Feebly Interacting Massive Particle (FIMP) model of Dark Matter (DM), one may express the inflationary energy scale $H_*$ as a function of three otherwise unrelated quantities, the DM isocurvature perturbation amplitude, its mass and its self-coupling constant, independently of the tensor-to-scalar ratio. The FIMP model assumes that there exists a real scalar particle that alone constitutes the DM content of the Universe and couples to the Standard Model via a Higgs portal. We consider carefully the various astrophysical, cosmological and model constraints, accounting also for variations in inflationary dynamics and the reheating history, to derive a robust estimate for $H_*$ that is confined to a relatively narrow range. We point out that, within the context of the FIMP DM model, one may thus determine $H_*$ reliably even in the absence of observable tensor perturbations.

We show how standard Newtonian N-body simulations can be interpreted in terms of the weak-field limit of general relativity by employing the recently developed Newtonian motion gauge. Our framework allows the inclusion of radiation perturbations and the non-linear evolution of matter. We show how to construct the weak-field metric by combining Newtonian simulations with results from Einstein-Boltzmann codes. We discuss observational effects on weak lensing and ray tracing, identifying important relativistic corrections.

We calculate the full probability density function (PDF) of inflationary curvature perturbations, even in the presence of large quantum backreaction. Making use of the stochastic-$\delta N$ formalism, two complementary methods are developed, one based on solving an ordinary differential equation for the characteristic function of the PDF, and the other based on solving a heat equation for the PDF directly. In the classical limit where quantum diffusion is small, we develop an expansion scheme that not only recovers the standard Gaussian PDF at leading order, but also allows us to calculate the first non-Gaussian corrections to the usual result. In the opposite limit where quantum diffusion is large, we find that the PDF is given by an elliptic theta function, which is fully characterised by the ratio between the squared width and height (in Planck mass units) of the region where stochastic effects dominate. We then apply these results to the calculation of the mass fraction of primordial black holes from inflation, and show that no more than $\sim 1$ $e$-fold can be spent in regions of the potential dominated by quantum diffusion. We explain how this requirement constrains inflationary potentials with two examples.

We study the stochastic distribution of spectator fields predicted in different slow-roll inflation backgrounds. Spectator fields have a negligible energy density during inflation but may play an important dynamical role later, even giving rise to primordial density perturbations within our observational horizon today. During de-Sitter expansion there is an equilibrium solution for the spectator field which is often used to estimate the stochastic distribution during slow-roll inflation. However slow roll only requires that the Hubble rate varies slowly compared to the Hubble time, while the time taken for the stochastic distribution to evolve to the de-Sitter equilibrium solution can be much longer than a Hubble time. We study both chaotic (monomial) and plateau inflaton potentials, with quadratic, quartic and axionic spectator fields. We give an adiabaticity condition for the spectator field distribution to relax to the de-Sitter equilibrium, and find that the de-Sitter approximation is never a reliable estimate for the typical distribution at the end of inflation for a quadratic spectator during monomial inflation. The existence of an adiabatic regime at early times can erase the dependence on initial conditions of the final distribution of field values. In these cases, spectator fields acquire sub-Planckian expectation values. Otherwise spectator fields may acquire much larger field displacements than suggested by the de-Sitter equilibrium solution. We quantify the information about initial conditions that can be obtained from the final field distribution. Our results may have important consequences for the viability of spectator models for the origin of structure, such as the simplest curvaton models.

Inflation in the early Universe is one of the most promising probes of gravity in the high-energy regime. However, observable scales give access to a limited window in the inflationary dynamics. In this essay, we argue that quantum corrections to the classical dynamics of cosmological fields allow us to probe much earlier epochs of the inflationary phase and extend this window by many orders of magnitude. We point out that both the statistics of cosmological fluctuations at observable scales, and the field displacements acquired by spectator fields that play an important role in many post-inflationary processes, are sensitive to a much longer phase of the inflationary epoch.

Bayesian inference techniques are used to investigate situations where an additional light scalar field is present during inflation and reheating. This includes (but is not limited to) curvaton-type models. We design a numerical pipeline where $\simeq 200$ inflaton setups $\times\, 10$ reheating scenarios $= 2000$ models are implemented and we present the results for a few prototypical potentials. We find that single-field models are remarkably robust under the introduction of light scalar degrees of freedom. Models that are ruled out at the single-field level are not improved in general, because good values of the spectral index and the tensor-to-scalar ratio can only be obtained for very fine-tuned values of the extra field parameters and/or when large non-Gaussianities are produced. The only exception is quartic large-field inflation, so that the best models after Planck are of two kinds: plateau potentials, regardless of whether an extra field is added or not, and quartic large-field inflation with an extra light scalar field, in some specific reheating scenarios. Using Bayesian complexity, we also find that more parameters are constrained for the models we study than for their single-field versions. This is because the added parameters not only contribute to the reheating kinematics but also to the cosmological perturbations themselves, to which the added field contributes. The interplay between these two effects lead to a suppression of degeneracies that is responsible for having more constrained parameters.

Initial conditions for (Newtonian) cosmological N-body simulations are usually set by re-scaling the present-day power spectrum obtained from linear (relativistic) Boltzmann codes to the desired initial redshift of the simulation. This back-scaling method can account for the effect of inhomogeneous residual thermal radiation at early times, which is absent in the Newtonian simulations. We analyse this procedure from a fully relativistic perspective, employing the recently-proposed Newtonian motion gauge framework. We find that N-body simulations for LambdaCDM cosmology starting from back-scaled initial conditions can be self-consistently embedded in a relativistic space-time with first-order metric potentials calculated using a linear Boltzmann code. This space-time coincides with a simple "N-body gauge" for z<50 for all observable modes. Care must be taken, however, when simulating non-standard cosmologies. As an example, we analyse the back-scaling method in a cosmology with decaying dark matter, and show that metric perturbations become large at early times in the back-scaling approach, indicating a breakdown of the perturbative description. We suggest a suitable "forwards approach" for such cases.

We derive the first systematic observational constraints on reheating in models of inflation where an additional light scalar field contributes to primordial density perturbations and affects the expansion history during reheating. This encompasses the original curvaton model but also covers a larger class of scenarios. We find that, compared to the single-field case, lower values of the energy density at the end of inflation and of the reheating temperature are preferred when an additional scalar field is introduced. For instance, if inflation is driven by a quartic potential, which is one of the most favoured models when a light scalar field is added, the upper bound $T_{\mathrm{reh}}<5\times 10^{4}\,\mathrm{GeV}$ on the reheating temperature $T_{\mathrm{reh}}$ is derived, and the implications of this value on post-inflationary physics are discussed. The information gained about reheating is also quantified and it is found that it remains modest in plateau inflation (though still larger than in the single-field version of the model) but can become substantial in quartic inflation. The role played by the vev of the additional scalar field at the end of inflation is highlighted, and opens interesting possibilities for exploring stochastic inflation effects that could determine its distribution.

Stochastic effects in generic scenarios of inflation with multiple fields are investigated. First passage time techniques are employed to calculate the statistical moments of the number of inflationary $e$-folds, which give rise to all correlation functions of primordial curvature perturbations through the stochastic $\delta N$ formalism. The number of fields is a critical parameter. The probability of exploring arbitrarily large-field regions of the potential becomes non-vanishing when more than two fields are driving inflation. The mean number of $e$-folds can be infinite, depending on the number of fields; for plateau potentials, this occurs even with one field. In such cases, correlation functions of curvature perturbations are infinite. They can, however, be regularised if a reflecting (or absorbing) wall is added at large energy or field value. The results are found to be independent of the exact location of the wall and this procedure is, therefore, well-defined for a wide range of cutoffs, above or below the Planck scale. Finally, we show that, contrary to single-field setups, multi-field models can yield large stochastic corrections even at sub-Planckian energy, opening interesting prospects for probing quantum effects on cosmological fluctuations.

The standard numerical tools for studying non-linear collapse of matter are Newtonian $N$-body simulations. Previous work has shown that these simulations are in accordance with General Relativity (GR) up to first order in perturbation theory, provided that the effects from radiation can be neglected. In this paper we show that the present day matter density receives more than 1$\%$ corrections from radiation on large scales if Newtonian simulations are initialised before $z=50$. We provide a relativistic framework in which \emph{unmodified} Newtonian simulations are compatible with linear GR even in the presence of radiation. Our idea is to use GR perturbation theory to keep track of the evolution of relativistic species and the relativistic space-time consistent with the Newtonian trajectories computed in $N$-body simulations. If metric potentials are sufficiently small, they can be computed using a first-order Einstein--Boltzmann code such as CLASS. We make this idea rigorous by defining a class of GR gauges, the \emph{Newtonian motion} gauges, which are defined such that matter particles follow Newtonian trajectories. We construct a simple example of a relativistic space-time within which unmodified Newtonian simulations can be interpreted.

Newtonian simulations are routinely used to examine the matter dynamics on non-linear scales. However, even on these scales, Newtonian gravity is not a complete description of gravitational effects. A post-Friedmann approach shows that the leading order correction to Newtonian theory is a vector potential in the metric. This vector potential can be calculated from N-body simulations, requiring a method for extracting the velocity field. Here, we present the full details of our calculation of the post-Friedmann vector potential, using the Delauney Tesselation Field Estimator (DTFE) code. We include a detailed examination of the robustness of our numerical result, including the effects of box size and mass resolution on the extracted fields. We present the power spectrum of the vector potential and find that the power spectrum of the vector potential is $\sim 10^5$ times smaller than the power spectrum of the fully non-linear scalar gravitational potential at redshift zero. Comparing our numerical results to perturbative estimates, we find that the fully non-linear result can be more than an order of magnitude larger than the perturbative estimate on small scales. We extend the analysis of the vector potential to multiple redshifts, showing that this ratio persists over a range of scales and redshifts. We also comment on the implications of our results for the validity and interpretation of Newtonian simulations.

We present a fully relativistic calculation of the matter bispectrum at second order in cosmological perturbation theory assuming a Gaussian primordial curvature perturbation. For the first time we perform a full numerical integration of the bispectrum for both baryons and cold dark matter using the second-order Einstein-Boltzmann code, SONG. We review previous analytical results and provide an improved analytic approximation for the second-order kernel in Poisson gauge which incorporates Newtonian nonlinear evolution, relativistic initial conditions, the effect of radiation at early times and the cosmological constant at late times. Our improved kernel provides a percent level fit to the full numerical result at late times for most configurations, including both equilateral shapes and the squeezed limit. We show that baryon acoustic oscillations leave an imprint in the matter bispectrum, making a significant impact on squeezed shapes.

Stochastic effects in multi-field inflationary scenarios are investigated. A hierarchy of diffusion equations is derived, the solutions of which yield moments of the numbers of inflationary $e$-folds. Solving the resulting partial differential equations in multi-dimensional field space is more challenging than the single-field case. A few tractable examples are discussed, which show that the number of fields is, in general, a critical parameter. When more than two fields are present for instance, the probability to explore arbitrarily large-field regions of the potential, otherwise inaccessible to single-field dynamics, becomes non-zero. In some configurations, this gives rise to an infinite mean number of $e$-folds, regardless of the initial conditions. Another difference with respect to single-field scenarios is that multi-field stochastic effects can be large even at sub-Planckian energy. This opens interesting new possibilities for probing quantum effects in inflationary dynamics, since the moments of the numbers of $e$-folds can be used to calculate the distribution of primordial density perturbations in the stochastic-$\delta N$ formalism.

Measurements of the non-Gaussianity of the primordial density field have the power to considerably improve our understanding of the physics of inflation. Indeed, if we can increase the precision of current measurements by an order of magnitude, a null-detection would rule out many classes of scenarios for generating primordial fluctuations. Large-scale galaxy redshift surveys represent experiments that hold the promise to realise this goal. Thus, we model the galaxy bispectrum and forecast the accuracy with which it will probe the parameter $f_{\rm NL}$, which represents the degree of primordial local-type non Gaussianity. Specifically, we address the problem of modelling redshift space distortions (RSD) in the tree-level galaxy bispectrum including $f_{\rm NL}$. We find novel contributions associated with RSD, with the characteristic large scale amplification induced by local-type non-Gaussianity. These RSD effects must be properly accounted for in order to obtain un-biased measurements of $f_{\rm NL}$ from the galaxy bispectrum. We propose an analytic template for the monopole which can be used to fit against data on large scales, extending models used in the recent measurements. Finally, we perform idealised forecasts on $\sigma_{f_{\rm NL}}$ -- the accuracy of the determination of local non-linear parameter $f_{\rm NL}$ -- from measurements of the galaxy bispectrum. Our findings suggest that current surveys can in principle provide $f_{\rm NL}$ constraints competitive with Planck, and future surveys could improve them further.

In General Relativity, the constraint equation relating metric and density perturbations is inherently nonlinear, leading to an effective non-Gaussianity in the dark matter density field on large scales - even if the primordial metric perturbation is Gaussian. Intrinsic non-Gaussianity in the large-scale dark matter overdensity in GR is real and physical. However, the variance smoothed on a local physical scale is not correlated with the large-scale curvature perturbation, so that there is no relativistic signature in the galaxy bias when using the simplest model of bias. It is an open question whether the observable mass proxies such as luminosity or weak lensing correspond directly to the physical mass in the simple halo bias model. If not, there may be observables that encode this relativistic signature.

An interaction between the vacuum energy and dark matter is an intriguing possibility which may offer a way of solving the cosmological constant problem. Adopting a general prescription for momentum exchange between the two dark components, we reconstruct $\alpha(a)$, the temporal evolution of the coupling strength between dark matter and vacuum energy, in a nonparametric Bayesian approach using combined observational data sets from the cosmic microwave background, supernovae and large scale structure. An evolving interaction between the vacuum energy and dark matter removes some of the tensions between different data sets. However, it is not preferred over $\Lambda$CDM in the Bayesian sense, as improvement in the fit is not sufficient to compensate for the increase in the volume of the parameter space.

The initial conditions for Newtonian $N$-body simulations are usually generated by applying the Zel'dovich approximation to the initial displacements of the particles using an initial power spectrum of density fluctuations generated by an Einstein-Boltzmann solver. We show that in most gauges the initial displacements generated in this way receive a first-order relativistic correction. We define a new gauge, the $N$-body gauge, in which this relativistic correction vanishes and show that a conventional Newtonian $N$-body simulation includes all first-order relativistic contributions (in the absence of radiation) if we identify the coordinates in Newtonian simulations with those in the relativistic $N$-body gauge.

We investigate whether the predictions of single-field models of inflation are robust under the introduction of additional scalar degrees of freedom, and whether these extra fields change the potentials for which the data show the strongest preference. We study the situation where an extra light scalar field contributes both to the total curvature perturbations and to the reheating kinematic properties. Ten reheating scenarios are identified, and all necessary formulas allowing a systematic computation of the predictions for this class of models are derived. They are implemented in the public library ASPIC, which contains more than 75 single-field potentials. This paves the way for a forthcoming full Bayesian analysis of the problem. A few representative examples are displayed and discussed.

In this paper we examine cosmological weak lensing on non-linear scales and show that there are Newtonian and relativistic contributions and that the latter can also be extracted from standard Newtonian simulations. We use the post-Friedmann formalism, a post-Newtonian type framework for cosmology, to derive the full weak-lensing deflection angle valid on non-linear scales for any metric theory of gravity. We show that the only contributing term that is quadratic in the first order deflection is the expected Born correction and lens-lens coupling term. We use this deflection angle to analyse the vector and tensor contributions to the E- and B- mode cosmic shear power spectra. In our approach, once the gravitational theory has been specified, the metric components are related to the matter content in a well-defined manner. Specifying General Relativity, we write down a complete set of equations for a GR$+\Lambda$CDM universe for computing all of the possible lensing terms from Newtonian N-body simulations. We illustrate this with the vector potential and show that, in a GR$+\Lambda$CDM universe, its contribution to the E-mode is negligible with respect to that of the conventional Newtonian scalar potential, even on non-linear scales. Thus, under the standard assumption that Newtonian N-body simulations give a good approximation of the matter dynamics, we show that the standard ray tracing approach gives a good description for a $\Lambda$CDM cosmology.

We discuss the question of gauge choice when analysing relativistic density perturbations at second order. We compare Newtonian and General Relativistic approaches. Some misconceptions in the recent literature are addressed. We show that the comoving-synchronous gauge is the unique gauge in General Relativity that corresponds to the Lagrangian frame and is entirely appropriate to describe the matter overdensity at second order. The comoving-synchronous gauge is the simplest gauge in which to describe Lagrangian bias at second order.

Primordial non-Gaussianity can lead to a scale-dependent bias in the density of collapsed halos relative to the underlying matter density. The galaxy power spectrum already provides constraints on local-type primordial non-Gaussianity complementary those from the cosmic microwave background (CMB), while the bispectrum contains additional shape information and has the potential to outperform CMB constraints in future. We develop the bias model for the halo density contrast in the presence of local-type primordial non-Gaussianity, deriving a bivariate expansion up to second order in terms of the local linear matter density contrast and the local gravitational potential in Lagrangian coordinates. Non-linear evolution of the matter density introduces a non-local tidal term in the halo model. Furthermore, the presence of local-type non-Gaussianity in the Lagrangian frame leads to a novel non-local convective term in the Eulerian frame, that is proportional to the displacement field when going beyond the spherical collapse approximation. We use an extended Press-Schechter approach to evaluate the halo mass function and thus the halo bispectrum. We show that including these non-local terms in the halo bispectra can lead to corrections of up to $25\%$ for some configurations, on large scales or at high redshift.

The cosmic microwave background (CMB) radiation provides a remarkable window onto the early universe, revealing its composition and structure. In these lectures we review and discuss the physics underlying the main features of the CMB.

We explore the hemispherical asymmetry predicted in cosmic microwave background polarization when there is an asymmetry in temperature anisotropies due to primordial perturbations. We consider the cases of asymmetries due to adiabatic and isocurvature modes, and tensor perturbations. We show that the asymmetry in the TE, EE and/or BB correlations can be substantially larger than those in the TT power spectrum in certain cases. The relative asymmetry in the different cross-correlations, as well as the angular scale dependence, can in principle distinguish between different origins for the asymmetry.

Using an effective field theory approach to inflation, we examine novel properties of the spectrum of inflationary tensor fluctuations, that arise when breaking some of the symmetries or requirements usually imposed on the dynamics of perturbations. During single-clock inflation, time-reparameterization invariance is broken by a time-dependent cosmological background. In order to explore more general scenarios, we consider the possibility that spatial diffeomorphism invariance is also broken by effective mass terms or by derivative operators for the metric fluctuations in the Lagrangian. We investigate the cosmological consequences of the breaking of spatial diffeomorphisms, focussing on operators that affect the power spectrum of fluctuations. We identify the operators for tensor fluctuations that can provide a blue spectrum without violating the null energy condition, and operators for scalar fluctuations that lead to non-conservation of the comoving curvature perturbation on superhorizon scales even in single-clock inflation. In the last part of our work, we also examine the consequences of operators containing more than two spatial derivatives, discussing how they affect the sound speed of tensor fluctuations, and showing that they can mimic some of the interesting effects of symmetry breaking operators, even in scenarios that preserve spatial diffeomorphism invariance.

We investigate whether non-adiabatic perturbations from inflation could produce an asymmetric distribution of temperature anisotropies on large angular scales in the cosmic microwave background (CMB). We use a generalised non-linear $\delta N$ formalism to calculate the non-Gaussianity of the primordial density and isocurvature perturbations due to the presence of non-adiabatic, but approximately scale-invariant field fluctuations during multi-field inflation. This local-type non-Gaussianity leads to a correlation between very long wavelength inhomogeneities, larger than our observable horizon, and smaller scale fluctuations in the radiation and matter density. Matter isocurvature perturbations contribute primarily to low CMB multipoles and hence can lead to a hemispherical asymmetry on large angular scales, with negligible asymmetry on smaller scales. In curvaton models, where the matter isocurvature perturbation is partly correlated with the primordial density perturbation, we are unable to obtain a significant asymmetry on large angular scales while respecting current observational constraints on the observed quadrupole. However in the axion model, where the matter isocurvature and primordial density perturbations are uncorrelated, we find it may be possible to obtain a significant asymmetry due to isocurvature modes on large angular scales. Such an isocurvature origin for the hemispherical asymmetry would naturally give rise to a distinctive asymmetry in the CMB polarisation on large scales.

We show that a general late-time interaction between cold dark matter and vacuum energy is favoured by current cosmological datasets. We characterize the strength of the coupling by a dimensionless parameter $q_V$ that is free to take different values in four redshift bins from the primordial epoch up to today. This interacting scenario is in agreement with measurements of cosmic microwave background temperature anisotropies from the Planck satellite, supernovae Ia from Union 2.1 and redshift space distortions from a number of surveys, as well as with combinations of these different datasets. We show that a non-zero interaction is very likely at late times. We then focus on the case $q_V\not=0$ in a single low-redshift bin, obtaining a nested one parameter extension of the standard $\Lambda$CDM model. We study the Bayesian evidence, with respect to $\Lambda$CDM, of this late-time interaction model, finding moderate evidence for an interaction starting at $z=0.9$, dependent upon the prior range chosen for the interaction strength parameter $q_V$. For this case the null interaction ($q_V=0$, i.e.$\Lambda$CDM) is excluded at 99% c.l..

We show how the non-linearity of general relativity generates a characteristic non-Gaussian signal in cosmological large-scale structure that we calculate at all perturbative orders in a large scale limit. Newtonian gravity and general relativity provide complementary theoretical frameworks for modelling large-scale structure in $\Lambda$CDM cosmology; a relativistic approach is essential to determine initial conditions which can then be used in Newtonian simulations studying the non-linear evolution of the matter density. Most inflationary models in the very early universe predict an almost Gaussian distribution for the primordial metric perturbation, $\zeta$. However, we argue that it is the Ricci curvature of comoving-orthogonal spatial hypersurfaces, $R$, that drives structure formation at large scales. We show how the non-linear relation between the spatial curvature, $R$, and the metric perturbation, $\zeta$, translates into a specific non-Gaussian contribution to the initial comoving matter density that we calculate for the simple case of an initially Gaussian $\zeta$. Our analysis shows the non-linear signature of Einstein's gravity in large-scale structure.

We compute the bispectrum induced in the cosmic microwave background (CMB) temperature and polarisation by the evolution of the primordial density perturbations using the second-order Boltzmann code SONG. We show that adding polarisation increases the signal-to-noise ratio by a factor four with respect to temperature alone and we estimate the observability of this intrinsic bispectrum and the bias it induces on measurements of primordial non-Gaussianity. When including all physical effects except the late-time non-linear evolution, we find for the intrinsic bispectrum a signal-to-noise of $S/N=3.8,\,2.9,\,1.6$ and $0.5$ for, respectively, an ideal experiment with an angular resolution of $\ell_\text{max}=3000$, the proposed CMB surveys PRISM and COrE, and Planck's polarised data; the bulk of this signal comes from the $E$-polarisation and from squeezed configurations. We discuss how CMB lensing is expected to reduce these estimates as it suppresses the bispectrum for squeezed configurations and contributes to the noise in the estimator. We find that the presence of the intrinsic bispectrum will bias a measurement of primordial non-Gaussianity of local type by $f_{NL}^\text{bias}=0.66$ for an ideal experiment with $\ell_\text{max}=3000$. Finally, we verify the robustness of our results by reproducing the analytical approximation for the squeezed-limit bispectrum in the general polarised case.

We present improved constraints on an interacting vacuum model using updated astronomical observations including the first data release from Planck. We consider a model with one dimensionless parameter, $\alpha$, describing the interaction between dark matter and vacuum energy (with fixed equation of state $w=-1$). The background dynamics correspond to a generalised Chaplygin gas cosmology, but the perturbations have a zero sound speed. The tension between the value of the Hubble constant, $H_0$, determined by Planck data plus WMAP polarisation (Planck+WP) and that determined by the Hubble Space Telescope (HST) can be alleviated by energy transfer from dark matter to vacuum ($\alpha>0$). A positive $\alpha$ increases the allowed values of $H_0$ due to parameter degeneracy within the model using only CMB data. Combining with additional datasets of including supernova type Ia (SN Ia) and baryon acoustic oscillation (BAO), we can significantly tighten the bounds on $\alpha$. Redshift-space distortions (RSD), which constrain the linear growth of structure, provide the tightest constraints on vacuum interaction when combined with Planck+WP, and prefer energy transfer from vacuum to dark matter ($\alpha<0$) which suppresses the growth of structure. Using the combined datasets of Planck+WP+Union2.1+BAO+RSD, we obtain the constraint on $\alpha$ to be $-0.083<\alpha<-0.006$ (95% C.L.), allowing low $H_0$ consistent with the measurement from 6dF Galaxy survey. This interacting vacuum model can alleviate the tension between RSD and Planck+WP in the $\Lambda$CDM model for $\alpha<0$, or between HST measurements of $H_0$ and Planck+WP for $\alpha>0$, but not both at the same time.

We estimate the B-polarisation induced in the Cosmic Microwave Background by the non-linear evolution of density perturbations. Using the second-order Boltzmann code SONG, our analysis incorporates, for the first time, all physical effects at recombination. We also include novel contributions from the redshift part of the Boltzmann equation and from the bolometric definition of the temperature in the presence of polarisation. The remaining line-of-sight terms (lensing and time-delay) have previously been studied and must be calculated non-perturbatively. The intrinsic B-mode polarisation is present independent of the initial conditions and might contaminate the signal from primordial gravitational waves. We find this contamination to be comparable to a primordial tensor-to-scalar ratio of $r\simeq10^{-7}$ at the angular scale $\ell\simeq100\,$, where the primordial signal peaks, and $r\simeq 5 \cdot 10^{-5}$ at $\ell\simeq700\,$, where the intrinsic signal peaks. Therefore, we conclude that the intrinsic B-polarisation from second-order effects is not likely to contaminate future searches of primordial gravitational waves.

The statistics of large-scale structure in the Universe can be used to probe non-Gaussianity of the primordial density field, complementary to existing constraints from the cosmic microwave background. In particular, the scale dependence of halo bias, which affects the halo distribution at large scales, represents a promising tool for analyzing primordial non-Gaussianity of local form. Future observations, for example, may be able to constrain the trispectrum parameter gNL that is difficult to study and constrain using the CMB alone. We investigate how galaxy and matter bispectra can distinguish between the two non-Gaussian parameters fNL and gNL, whose effects give nearly degenerate contributions to the power spectra. We use a generalization of the univariate bias approach, making the hypothesis that the number density of halos forming at a given position is a function of the local matter density contrast and of its local higher-order statistics. Using this approach, we calculate the halo-matter bispectra and analyze their properties. We determine a connection between the sign of the halo bispectrum on large scales and the parameter gNL. We also construct a combination of halo and matter bispectra that is sensitive to fNL, with little contamination from gNL. We study both the case of single and multiple sources to the primordial gravitational potential, discussing how to extend the concept of stochastic halo bias to the case of bispectra. We use a specific halo mass-function to calculate numerically the bispectra in appropriate squeezed limits, confirming our theoretical findings.

The goal of the present paper is to set initial conditions for structure formation at non-linear order, consistent with general relativity, while also allowing for primordial non-Gaussianity. We use the non-linear continuity and Raychauduri equations, which together with the non-linear energy constraint determine the evolution of the matter density fluctuation in general relativity. We solve this equations at first and second order in a perturbative expansion, recovering and extending previous results derived in the matter-dominated limit and in the Newtonian regime. We present a second-order solution for the comoving density contrast in a $\Lambda$CDM universe, identifying non-linear contributions coming from the Newtonian growing mode, primordial non-Gaussianity and intrinsic non-Gaussianity, due to the essential non-linearity of the relativistic constraint equations. We discuss the application of these results to initial conditions in N-body simulations, showing that relativistic corrections mimic a non-zero non-linear parameter $f_{\rm NL}$.

The new Planck constraints on the local bispectrum parameter fnl are about 10^5 times tighter than the current constraints on the trispectrum parameter gnl, which means that the allowed numerical values of the second and third order terms in the perturbative expansion of the curvature perturbation are comparable. We show that a consequence of this is that if gnl is large enough to be detectable, then it will induce a large variation between the observable value of fnl and its value in a larger inflated volume. Even if there were only a few extra efoldings between the beginning of inflation and horizon crossing of our Hubble horizon, an observably large gnl means that fnl is unlikely to be as small as its current constraint, regardless of its true background value. This result is very general, it holds regardless of how many fields contributed to the curvature perturbation. We also generalise this result to other shapes of non-Gaussianity, beyond the local model. We show that the variance of the 3-point function in the squeezed limit is bounded from below by the square of the squeezed limit of the 4-point function.

We present the first calculation of an intrinsically relativistic quantity in fully non-linear cosmolog- ical large-scale structure studies. Traditionally, non-linear structure formation in standard {\Lambda}CDM cosmology is studied using N-body simulations, based on Newtonian gravitational dynamics on an expanding background. When one derives the Newtonian regime in a way that is a consistent ap- proximation to the Einstein equations, a gravito-magnetic vector potential - giving rise to frame dragging - is present in the metric in addition to the usual Newtonian scalar potential. At leading order, this vector potential does not affect the matter dynamics, thus it can be computed from Newtonian N-body simulations. We explain how we compute the vector potential from simulations in {\Lambda}CDM and examine its magnitude relative to the scalar potential. We also discuss some possible observable effects.

We study a model of multiple-field DBI inflation leading to mixed form of primordial non-Gaussianity, including equilateral and local bispectrum shapes. We present a general formalism based on the Hamilton-Jacobi approach, allowing us to go beyond slow-roll, combining the three-point function for the fields at Hubble-exit with the non-linear evolution of super-Hubble scales. We are able to obtain analytic results by taking a separable Ansatz for the Hubble rate. We find general expressions for both the equilateral and local type non-Gaussianity parameter fNL. The equilateral non-Gaussianity includes the usual enhancement for small sound speeds, but multiplied by an analytic factor which can lead to a suppression. We illustrate our results with two scenarios. In the first model, previously found to have detectable local non-Gaussianity, we find that the equilateral signal is not sufficiently suppressed to evade current observational bounds. In our second scenario we construct a model which exhibits both a detectable equilateral fNL and a negative local fNL.

Any unified dark matter cosmology can be decomposed into dark matter interacting with vacuum energy, without introducing any additional degrees of freedom. We present observational constraints on an interacting vacuum plus dark energy corresponding to a generalised Chaplygin gas cosmology. We consider two distinct models for the interaction leading to either a barotropic equation of state or dark matter that follows geodesics, corresponding to a rest-frame sound speed equal to the adiabatic sound speed or zero sound speed, respectively. For the barotropic model, the most stringent constraint on $\alpha$ comes from the combination of CMB+SNIa+LSS(m) gives $\alpha<5.66\times10^{-6}$ at the 95% confidence level, which indicates that the barotropic model must be extremely close to the $\Lambda$CDM cosmology. For the case where the dark matter follows geodesics, perturbations have zero sound speed, and CMB+SNIa+gISW then gives the much weaker constraint $-0.15<\alpha<0.26$ at the 95% confidence level.

We develop a new, efficient code for solving the second-order Einstein-Boltzmann equations, and use it to estimate the intrinsic CMB non-Gaussianity arising from the non-linear evolution of density perturbations. The full calculation involves contributions from recombination and less tractable contributions from terms integrated along the line of sight. We investigate the bias that this intrinsic bispectrum implies for searches of primordial non-Gaussianity. We find that the inclusion or omission of certain line of sight terms can make a large impact. When including all physical effects but lensing and time-delay, we find that the local-type f_nl would be biased by f_nl ~ 0.5, below the expected sensitivity of the Planck satellite. The speed of our code allows us to confirm the robustness of our results with respect to a number of numerical parameters.

We analyze the dynamical system defined by a universe filled with a barotropic fluid plus a scalar field with modified kinetic term of the form L = F (X) - V (phi). After a suitable choice of variables that allows us to study the phase space of the system we obtain the critical points and their stability. We find that they reduce to the ones defined for the canonical case when F (X) = X. We also study the field energy conditions to have a nonsingular bounce.

Inflation is the leading paradigm for explaining the origin of primordial density perturbations and the observed temperature fluctuations of the cosmic microwave background. However many open questions remain, in particular whether one or more scalar fields were present during inflation and how they contributed to the primordial density perturbation. We propose a new observational test of whether multiple fields, or only one (not necessarily the inflaton) generated the perturbations. We show that our test, relating the bispectrum and trispectrum, is protected against loop corrections at all orders, unlike previous relations.

We study primordial density perturbations generated by the late decay of a curvaton field whose decay rate may be modulated by the local value of another isocurvature field, analogous to models of modulated reheating at the end of inflation. We calculate the primordial density perturbation and its local-type non-Gaussianity using the sudden-decay approximation for the curvaton field, recovering standard curvaton and modulated reheating results as limiting cases. We verify the Suyama-Yamaguchi inequality between bispectrum and trispectrum parameters for the primordial density field generated by multiple field fluctuations, and find conditions for the bound to be saturated.

We propose a method to probe higher-order correlators of the primordial density field through the inhomogeneity of local non-Gaussian parameters, such as f_NL, measured within smaller patches of the sky. Correlators between n-point functions measured in one patch of the sky and k-point functions measured in another patch depend upon the (n+k)-point functions over the entire sky. The inhomogeneity of non-Gaussian parameters may be a feasible way to detect or constrain higher-order correlators in local models of non-Gaussianity, as well as to distinguish between single and multiple-source scenarios for generating the primordial density perturbation, and more generally to probe the details of inflationary physics.

Vacuum energy is a simple model for dark energy driving an accelerated expansion of the universe. If the vacuum energy is inhomogeneous in spacetime then it must be interacting. We present the general equations for a spacetime-dependent vacuum energy in cosmology, including inhomogeneous perturbations. We show how any dark energy cosmology can be described by an interacting vacuum+matter. Different models for the interaction can lead to different behaviour (e.g., sound speed for dark energy perturbations) and hence could be distinguished by cosmological observations. As an example we present the cosmic microwave microwave background anisotropies and the matter power spectrum for two different versions of a generalised Chaplygin gas cosmology.

We study the primordial density perturbations and non-Gaussianities generated from the combined effects of an inhomogeneous end of inflation and curvaton decay in hybrid inflation. This dual role is played by a single isocurvature field which is massless during inflation but acquire a mass at the end of inflation via the waterfall phase transition. We calculate the resulting primordial non-Gaussianity characterized by the non-linearity parameter, $f_{NL}$, recovering the usual end-of-inflation result when the field decays promptly and the usual curvaton result if the field decays sufficiently late.

We study the effect of non-trivial sound speeds on local-type non-Gaussianity during multiple-field inflation. To this end, we consider a model of multiple-field DBI and use the deltaN formalism to track the super-horizon evolution of perturbations. By adopting a sum separable Hubble parameter we derive analytic expressions for the relevant quantities in the two-field case, valid beyond slow variation. We find that non-trivial sound speeds can, in principle, curve the trajectory in such a way that significant local-type non-Gaussianity is produced. Deviations from slow variation, such as rapidly varying sound speeds, enhance this effect. To illustrate our results we consider two-field inflation in the tip regions of two warped throats and find large local-type non-Gaussianity produced towards the end of the inflationary process.

We characterise the primordial perturbations produced due to both inflaton and curvaton fluctuations in models where the curvaton has a quadratic, cosine or hyperbolic potential, and the inflaton potential is characterised by the usual slow-roll parameters. Isocurvature curvaton field perturbations can produce significant non-Gaussianity in the primordial density field, in contrast with adiabatic inflaton field perturbations which produce negligible non-Gaussianity for canonical scalar fields. A non-self-interacting curvaton with quadratic potential produces a local-type non-Gaussianity that is well described by the non-linearity parameter fNL, which may be scale-dependent when the inflaton perturbations dominate the power spectrum. We show how observational bounds on non-linearity parameters and the tensor-scalar ratio can be used to constrain curvaton and inflaton parameters. We find a consistency relation between the bispectrum and trispectrum parameters in a mixed inflaton-curvaton model for a quadratic curvaton potential. Self-interaction terms in the curvaton potential can lead to both a large trispectrum parameter, gNL, and scale-dependence of the non-linearity parameters.

Vacuum energy remains the simplest model of dark energy which could drive the accelerated expansion of the Universe without necessarily introducing any new degrees of freedom. Inhomogeneous vacuum energy is necessarily interacting in general relativity. Although the four-velocity of vacuum energy is undefined, an interacting vacuum has an energy transfer and the vacuum energy defines a particular foliation of spacetime with spatially homogeneous vacuum energy in cosmological solutions. It is possible to give a consistent description of vacuum dynamics and in particular the relativistic equations of motion for inhomogeneous perturbations given a covariant prescription for the vacuum energy, or equivalently the energy transfer four-vector, and we construct gauge-invariant vacuum perturbations. We show that any dark energy cosmology can be decomposed into an interacting vacuum+matter cosmology whose inhomogeneous perturbations obey simple first-order equations.

Local non-Gaussianity, parametrized by $f_{\rm NL}$, introduces a scale-dependent bias that is strongest at large scales, precisely where General Relativistic (GR) effects also become significant. With future data, it should be possible to constrain $f_{\rm NL} = {\cal O}(1)$ with high redshift surveys. GR corrections to the power spectrum and ambiguities in the gauge used to define bias introduce effects similar to $f_{\rm NL}= {\cal O}(1)$, so it is essential to disentangle these effects. For the first time in studies of primordial non-Gaussianity, we include the consistent GR calculation of galaxy power spectra, highlighting the importance of a proper definition of bias. We present observable power spectra with and without GR corrections, showing that an incorrect definition of bias can mimic non-Gaussianity. However, these effects can be distinguished by their different redshift and scale dependence, so as to extract the true primordial non-Gaussianity.

We consider a simple model of cosmological collapse driven by canonical fields with exponential potentials. We generalise the two-field ekpyrotic collapse to consider non-orthogonal or tilted potentials and give the general condition for isocurvature field fluctuations to have a scale-invariant spectrum in this model. In particular we show that tilted potentials allow for a slightly red spectrum of perturbations as required by current observations. However a red spectrum of fluctuations implies that the two-field ekpyrotic phase must have a finite duration and requires a preceding phase which sets the initial conditions for what otherwise appears to be a fine-tuned trajectory in the phase space.

In this paper we consider how non-Gaussianity of the primordial density perturbation and the amplitude of gravitational waves from inflation can be used to determine parameters of the curvaton scenario for the origin of structure. We show that in the simplest quadratic model, where the curvaton evolves as a free scalar field, measurement of the bispectrum relative to the power spectrum, fNL, and the tensor-to-scalar ratio can determine both the expectation value of the curvaton field during inflation and its dimensionless decay rate relative to the curvaton mass. We show how these predictions are altered by the introduction of self-interactions, in models where higher-order corrections are determined by a characteristic mass scale and discuss how additional information about primordial non-Gaussianity and scale dependence may constrain curvaton interactions.

We calculate the scale dependence of the bispectrum and trispectrum in (quasi) local models of non-Gaussian primordial density perturbations, and characterize this scale dependence in terms of new observable parameters. They can help to discriminate between models of inflation, since they are sensitive to properties of the inflationary physics that are not probed by the standard observables. We find consistency relations between these parameters in certain classes of models. We apply our results to a scenario of modulated reheating, showing that the scale dependence of non-Gaussianity can be significant. We also discuss the scale dependence of the bispectrum and trispectrum, in cases where one varies the shape as well as the overall scale of the figure under consideration. We conclude providing a formulation of the curvature perturbation in real space, which generalises the standard local form by dropping the assumption that f_NL and g_NL are constants.

We estimate large-scale curvature perturbations from isocurvature fluctuations in the waterfall field during hybrid inflation, in addition to the usual inflaton field perturbations. The tachyonic instability at the end of inflation leads to an explosive growth of super-Hubble scale perturbations, but they retain the steep blue spectrum characteristic of vacuum fluctuations in a massive field during inflation. The power spectrum thus peaks around the Hubble-horizon scale at the end of inflation. We extend the usual delta-N formalism to include the essential role of these small fluctuations when estimating the large-scale curvature perturbation. The resulting curvature perturbation due to fluctuations in the waterfall field is second-order and the spectrum is expected to be of order 10^{-54} on cosmological scales.

We consider possible scale-dependence of the non-linearity parameter f_NL in local and quasi-local models of non-Gaussian primordial density perturbations. In the simplest model where the primordial perturbations are a quadratic local function of a single Gaussian field then f_NL is scale-independent by construction. However scale-dependence can arise due to either a local function of more than one Gaussian field, or due to non-linear evolution of modes after horizon-exit during inflation. We show that the scale dependence of f_NL is typically first order in slow-roll. For some models this may be observable with experiments such as Planck provided that f_NL is close to the current observational bounds.

We study perturbations of a scalar field cosmology in Horava-Lifshitz gravity, adopting the most general setup without detailed balance but with the projectability condition. We derive the generalized Klein-Gordon equation, which is sixth-order in spatial derivatives. Then we investigate scalar field perturbations coupled to gravity in a flat Friedmann-Robertson-Walker background. In the sub-horizon regime, the metric and scalar field modes have independent oscillations with different frequencies and phases except in particular cases. On super-horizon scales, the perturbations become adiabatic during slow-roll inflation driven by a single field, and the comoving curvature perturbation is constant.

The non-Gaussian distribution of primordial perturbations has the potential to reveal the physical processes at work in the very early Universe. Local models provide a well-defined class of non-Gaussian distributions that arise naturally from the non-linear evolution of density perturbations on super-Hubble scales starting from Gaussian field fluctuations during inflation. I describe the delta-N formalism used to calculate the primordial density perturbation on large scales and then review several models for the origin of local primordial non-Gaussianity, including the cuvaton, modulated reheating and ekpyrotic scenarios. I include an appendix with a table of sign conventions used in specific papers.

We consider the stochastic background of gravitational waves produced during the radiation-dominated hot big bang as a constraint on the primordial density perturbation on comoving length scales much smaller than those directly probed by the cosmic microwave background or large-scale structure. We place weak upper bounds on the primordial density perturbation from current data. Future detectors such as BBO and DECIGO will place much stronger constraints on the primordial density perturbation on small scales.

We compute the second-order matching conditions for tensor metric perturbations at an abrupt change in the equation of state. For adiabatic perturbations on large scales the matching hypersurface coincides with a uniform-density hypersurface. We show that in the uniform-density gauge both the tensor perturbation and its time-derivative are continuous in this case. For non-adiabatic perturbations, the matching hypersurface need not coincide with a uniform-density hypersurface and the tensor perturbation in the uniform-density gauge may be discontinuous. However, we show that in the Poisson gauge both the tensor perturbation and its time-derivative are continuous for adiabatic or non-adiabatic perturbations. As an application we solve the evolution equation for second-order tensor perturbations on large scales for a constant equation of state and we use the matching conditions to evolve the solutions through the transition from an inflationary era to a radiation era. We show that in the radiation era the resulting free part of the large-scale tensor perturbation (constant mode) is slow-roll suppressed in both the uniform-density and Poisson gauges. Thus, we conclude that second-order gravitational waves from slow-roll inflation are suppressed.

We review the study of inhomogeneous perturbations about a homogeneous and isotropic background cosmology. We adopt a coordinate based approach, but give geometrical interpretations of metric perturbations in terms of the expansion, shear and curvature of constant-time hypersurfaces and the orthogonal timelike vector field. We give the gauge transformation rules for metric and matter variables at first and second order. We show how gauge invariant variables are constructed by identifying geometric or matter variables in physically-defined coordinate systems, and give the relations between many commonly used gauge-invariant variables. In particular we show how the Einstein equations or energy-momentum conservation can be used to obtain simple evolution equations at linear order, and discuss extensions to non-linear order. We present evolution equations for systems with multiple interacting fluids and scalar fields, identifying adiabatic and entropy perturbations. As an application we consider the origin of primordial curvature and isocurvature perturbations from field perturbations during inflation in the very early universe.

In this note we examine the derivation of scale-dependent bias due to primordial non-Gaussianity of the local type in the context of general relativity. We justify the use of the Poisson equation in general relativistic perturbation theory and thus the derivation of scale-dependent bias as a test of primordial non-Gaussianity, using the spherical collapse model. The corollary is that the form of scale-dependent bias does not receive general relativistic corrections on scales larger than the Hubble radius. This leads to a formally divergent correlation function for biased tracers of the mass distribution which we discuss.

We compute the non-Gaussianity of the curvature perturbation generated by ekpyrotic collapse with multiple fields. The transition from the multi-field scaling solution to a single-field dominated regime converts initial isocurvature field perturbations to an almost scale-invariant comoving curvature perturbation. In the specific model of two fields, $\phi_1$ and $\phi_2$, with exponential potentials, $-V_i \exp (-c_i \phi_i)$, we calculate the bispectrum of the resulting curvature perturbation. We find that the non-Gaussianity is dominated by non-linear evolution on super-Hubble scales and hence is of the local form. The non-linear parameter of the curvature perturbation is given by $f_{NL} = 5 c_j^2 /12$, where $c_j$ is the exponent of the potential for the field which becomes sub-dominant at late times. Since $c_j^2$ must be large, in order to generate an almost scale invariant spectrum, the non-Gaussianity is inevitably large. By combining the present observational constraints on $f_{\rm NL}$ and the scalar spectral index, the specific model studied in this paper is thus ruled out by current observational data.

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 summarize the utility of precise cosmic microwave background (CMB) polarization measurements as probes of the physics of inflation. We focus on the prospects for using CMB measurements to differentiate various inflationary mechanisms. In particular, a detection of primordial B-mode polarization would demonstrate that inflation occurred at a very high energy scale, and that the inflaton traversed a super-Planckian distance in field space. We explain how such a detection or constraint would illuminate aspects of physics at the Planck scale. Moreover, CMB measurements can constrain the scale-dependence and non-Gaussianity of the primordial fluctuations and limit the possibility of a significant isocurvature contribution. Each such limit provides crucial information on the underlying inflationary dynamics. Finally, we quantify these considerations by presenting forecasts for the sensitivities of a future satellite experiment to the inflationary parameters.

A scale-invariant spectrum of isocurvature perturbations is generated during collapse in the ekpyrotic scaling solution in models where multiple fields have steep negative exponential potentials. The scale invariance of the spectrum is realized by a tachyonic instability in the isocurvature field. This instability drives the scaling solution to the late time attractor that is the old ekpyrotic collapse dominated by a single field. We show that the transition from the scaling solution to the single field dominated ekpyrotic collapse automatically converts the initial isocurvature perturbations about the scaling solution to comoving curvature perturbations about the late-time attractor. The final amplitude of the comoving curvature perturbation is determined by the Hubble scale at the transition.

We investigate the generation of gravitational waves due to the gravitational instability of primordial density perturbations in an early matter-dominated era which could be detectable by experiments such as LIGO and LISA. We use relativistic perturbation theory to give analytic estimates of the tensor perturbations generated at second order by linear density perturbations. We find that large enhancement factors with respect to the naive second-order estimate are possible due to the growth of density perturbations on sub-Hubble scales. However very large enhancement factors coincide with a breakdown of linear theory for density perturbations on small scales. To produce a primordial gravitational wave background that would be detectable with LIGO or LISA from density perturbations in the linear regime requires primordial comoving curvature perturbations on small scales of order 0.02 for Advanced LIGO or 0.005 for LISA, otherwise numerical calculations of the non-linear evolution on sub-Hubble scales are required.

We show that simple scalar field models can give rise to curvature singularities in the effective Friedmann dynamics of Loop Quantum Cosmology (LQC). We find singular solutions for spatially flat Friedmann-Robertson-Walker cosmologies with a canonical scalar field and a negative exponential potential, or with a phantom scalar field and a positive potential. While LQC avoids big bang or big rip type singularities, we find sudden singularities where the Hubble rate is bounded, but the Ricci curvature scalar diverges. We conclude that the effective equations of LQC are not in themselves sufficient to avoid the occurrence of singularities.

We study non-linear primordial adiabatic and isocurvature perturbations and their non-Gaussianity. After giving a general formulation in the context of an extended delta N-formalism, we analyse in detail two illustrative examples. The first is a mixed curvaton-inflaton scenario in which fluctuations of both the inflaton and a curvaton (a light isocurvature field during inflation) contribute to the primordial density perturbation. The second example is that of double inflation involving two decoupled massive scalar fields during inflation. In the mixed curvaton-inflaton scenario we find that the bispectrum of primordial isocurvature perturbations may be large and comparable to the bispectrum of adiabatic curvature perturbations.

Several scenarios have been proposed in which primordial perturbations could originate from quantum vacuum fluctuations in a phase corresponding to a collapse phase (in an Einstein frame) preceding the Big Bang. I briefly review three models which could produce scale-invariant spectra during collapse: (1) curvature perturbations during pressureless collapse, (2) axion field perturbations in a pre big bang scenario, and (3) tachyonic fields during multiple-field ekpyrotic collapse. In the separate universes picture one can derive generalised perturbation equations to describe the evolution of large scale perturbations through a semi-classical bounce, assuming a large-scale limit in which inhomogeneous perturbations can be described by locally homogeneous patches. For adiabatic perturbations there exists a conserved curvature perturbation on large scales, but isocurvature perturbations can change the curvature perturbation through the non-adiabatic pressure perturbation on large scales. Different models for the origin of large scale structure lead to different observational predictions, including gravitational waves and non-Gaussianity.

We consider linear perturbation equations for long-wavelength scalar metric perturbations in generalised gravity, applicable to non-singular cosmological models including a bounce from collapse to expansion in the very early universe. We present the general form for the perturbation equations which follows from requiring that the inhomogeneous universe on large scales obeys the same local equations as the homogeneous Friedmann-Robertson-Walker background cosmology (the separate universes approach). In a pseudo-longitudinal gauge this becomes a homogeneous second-order differential equation for adiabatic perturbations, which reduces to the usual equation for the longitudinal gauge metric perturbation in general relativity with vanishing anisotropic stress. As an application we show that the scale-invariant spectrum of perturbations in the longitudinal gauge generated during an ekpyrotic collapse are not transferred to the growing mode of adiabatic density perturbations in the expanding phase in a simple bounce model.

We study a model where two scalar fields, that are subdominant during inflation, decay into radiation some time after inflation has ended but before primordial nucleosynthesis. Perturbations of these two curvaton fields can be responsible for the primordial curvature perturbation. We write down the full non-linear equations that relate the primordial perturbation to the curvaton perturbations on large scales, and solve them in a sudden-decay approximation. We calculate the power spectrum of the primordial perturbation, and finally go to second order to find the non-linearity parameter, fNL. Not surprisingly, we find large positive values of fNL if the energy densities of the curvatons are sub-dominant when they decay, as in the single curvaton case. But we also find a novel effect, which can be present only in multi-curvaton models: fNL becomes large even if the curvatons dominate the total energy density in the case when the inhomogeneous radiation produced by the first curvaton decay is diluted by the decay of a second nearly homogeneous curvaton. The minimum value min(fNL)=-5/4 which we find is the same as in the single-curvaton case. Using (non-)Gaussianity observations, Planck can be able to distinguish between single-field inflation and curvaton model. Hence it is important to derive theoretical predictions for curvaton model. From particle physics point of view it is more natural to assume multiple scalar fields (rather than just one ``curvaton'' in addition to inflaton). Our work updates the theoretical predictions of curvaton model to this case.

We present Feynman type diagrams for calculating the n-point function of the primordial curvature perturbation in terms of scalar field perturbations during inflation. The diagrams can be used to evaluate the corresponding terms in the n-point function at tree level or any required loop level. Rules are presented for drawing the diagrams and writing down the corresponding terms in real space and Fourier space. We show that vertices can be renormalised to automatically account for diagrams with dressed vertices. We apply these rules to calculate the primordial power spectrum up to two loops, the bispectrum including loop corrections, and the trispectrum.

We quantify the slow-roll corrections to primordial density perturbations arising from inflation driven by a four-dimensional scalar field with a monomial potential in a five-dimensional non-compact bulk spacetime. Although the difference between the classical brane-world solutions and standard four-dimensional solutions is large at early times, the change to the amplitude at late times of perturbations generated from quantum fluctuations is first-order in slow-roll parameters, leading to second-order slow-roll corrections to the spectral index. This confirms that the leading-order effects are correctly given by previous work in the literature.

We study a model where two scalar fields, that are subdominant during inflation, decay into radiation some time after inflation has ended but before primordial nucleosynthesis. Perturbations of these two curvaton fields can be responsible for the primordial curvature perturbation. We write down the full non-linear equations that relate the primordial perturbation to the curvaton perturbations on large scales, calculate the power spectrum of the primordial perturbation, and finally go to second order to find the non-linearity parameter, fNL. We find large positive values of fNL if the energy densities of the curvatons are sub-dominant when they decay, as in the single curvaton case. But we also find a large fNL even if the curvatons dominate the total energy density in the case when the inhomogeneous radiation produced by the first curvaton decay is diluted by the decay of a second nearly homogeneous curvaton. The minimum value min(fNL)=-5/4 which we find is the same as in the single-curvaton case.

A scale invariant spectrum of isocurvature perturbations is generated during collapse in the scaling solution in models where two or more fields have steep negative exponential potentials. The scale invariance of the spectrum is realised by a tachyonic instability in the isocurvature field. We show that this instability is due to the fact that the scaling solution is a saddle point in the phase space. The late time attractor is identified with a single field dominated ekpyrotic collapse in which a steep blue spectrum for isocurvature perturbations is found. Although quantum fluctuations do not necessarily to disrupt the classical solution, an additional preceding stage is required to establish classical homogeneity.

We calculate the power spectra of primordial curvature and isocurvature perturbations from a general two field inflation model at next-to-leading order correction in a slow-roll expansion. In particular we calculate the spectral indices to second order in slow-roll parameters. We show that the cross-correlation of the curvature and isocurvature perturbations at the time of Hubble-exit during inflation is non-zero at first-order in slow-roll parameters. We apply our results to different classes of inflation, including inflaton and curvaton scenarios. The spectrum of primordial gravitational waves, curvature and isocurvature perturbations obey generalised consistency relations in two-field inflation models. We give the first two consistency relations in an infinite hierarchy.

Inflation offers a simple model for very early evolution of our Universe and the origin of primordial perturbations on large scales. Over the last 25 years we have become familiar with the predictions of single-field models, but inflation with more than one light scalar field can alter preconceptions about the inflationary dynamics and our predictions for the primordial perturbations. I will discuss how future observational data could distinguish between inflation driven by one field, or many fields. As an example, I briefly review the curvaton as an alternative to the inflaton scenario for the origin of structure.

In this paper we quantise scalar perturbations in a Randall-Sundrum-type model of inflation where the inflaton field is confined to a single brane embedded in five-dimensional anti-de Sitter space-time. In the high energy regime, small-scale inflaton fluctuations are strongly coupled to metric perturbations in the bulk and gravitational back-reaction has a dramatic effect on the behaviour of inflaton perturbations on sub-horizon scales. This is in contrast to the standard four-dimensional result where gravitational back-reaction can be neglected on small scales. Nevertheless, this does not give rise to significant particle production, and the correction to the power spectrum of the curvature perturbations on super-horizon scales is shown to be suppressed by a slow-roll parameter. We calculate the complete first order slow-roll corrections to the spectrum of primordial curvature perturbations.

We use the delta N-formalism to describe the leading order contributions to the primordial power spectrum, bispectrum and trispectrum in multiple-field models of inflation at leading order in a perturbative expansion. In slow-roll models where the initial field fluctuations at Hubble-exit are nearly Gaussian, any detectable non-Gaussianity is expected to come from super-Hubble evolution. We show that the contribution to the primordial trispectrum can be described by two non-linearity parameters, tau_{NL} and g_{NL}, which are dependent upon the second and third derivatives of the local expansion with respect to the field values during inflation.

We discuss the gravitational wave background generated by primordial density perturbations evolving during the radiation era. At second-order in a perturbative expansion, density fluctuations produce gravitational waves. We calculate the power spectra of gravitational waves from this mechanism, and show that, in principle, future gravitational wave detectors could be used to constrain the primordial power spectrum on scales vastly different from those currently being probed by large-scale structure. As examples we compute the gravitational wave background generated by both a power-law spectrum on all scales, and a delta-function power spectrum on a single scale.

We use the delta N -formalism to investigate the non-Gaussianity of the primordial curvature perturbation in the curvaton scenario for the origin of structure. We numerically calculate the full probability distribution function allowing for the non-instantaneous decay of the curvaton and compare this with analytic results derived in the sudden-decay approximation. We also present results for the leading-order contribution to the primordial bispectrum and trispectrum. In the sudden-decay approximation we derive a fully non-linear expression relating the primordial perturbation to the initial curvaton perturbation. As an example of how non-Gaussianity provides additional constraints on model parameters, we show how the primordial bispectrum on CMB scales can be used to constrain variance on much smaller scales in the curvaton field. Our analytical and numerical results allow for multiple tests of primordial non-Gaussianity, and thus they can offer consistency tests of the curvaton scenario.

Recently, the primordial non-Gaussianity in the curvaton model has been predicted assuming sudden decay of the curvaton. We extend the calculation to non-instantaneous decay by employing delta N -formalism. The difference between the sudden-decay approximation and our numerical result is larger than 1% only if the non-linearity parameter is small, -1.16 < f_NL < 60. Thus it is safe to use the sudden-decay approximation when deriving constraints for the curvaton model from WMAP3 (f_NL < 114), but with the Planck forecast |f_NL| <5 one should employ the fully numerical result. Often, the curvaton perturbations $\delta\sigma$ have been assumed to be small compared to the background value of the curvaton field $\sigma_0$. Consequently, the variance $\Delta^2 = <\delta\sigma^2> / \sigma_0^2$ has been assumed to be negligible. However, the measurements of CMB or large-scale structure perturbation amplitude do not constrain the variance if the main contribution to it comes from the ultraviolet (UV) scales, i.e., from smaller than observable scales. We discuss how, even in this case, observational constraints on non-Gaussianity set an upper bound to the small scale variance, Delta^2_UV < 90.

We investigate the evolution of scalar metric perturbations across a sudden cosmological transition, allowing for an inhomogeneous surface stress at the transition leading to a discontinuity in the local expansion rate, such as might be expected in a big crunch/big bang event. We assume that the transition occurs when some function of local matter variables reaches a critical value, and that the surface stress is also a function of local matter variables. In particular we consider the case of a single scalar field and show that a necessary condition for the surface stress tensor to be perturbed at the transition is the presence of a non-zero intrinsic entropy perturbation of the scalar field. We present the matching conditions in terms of gauge-invariant variables assuming a sudden transition to a fluid-dominated universe with barotropic equation of state. For adiabatic perturbations the comoving curvature perturbation is continuous at the transition, while the Newtonian potential may be discontinuous if there is a discontinuity in the background Hubble expansion.

We study the bispectrum of the curvature perturbation on uniform energy density hypersurfaces in models of inflation with two scalar fields evolving simultaneously. In the case of a separable potential, it is possible to compute the curvature perturbation up to second order in the perturbations, generated on large scales due to the presence of non-adiabatic perturbations, by employing the $\delta N$-formalism, in the slow-roll approximation. In this case, we provide an analytic formula for the nonlinear parameter $f_{NL}$. We apply this formula to double inflation with two massive fields, showing that it does not generate significant non-Gaussianity; the nonlinear parameter at the end of inflation is slow-roll suppressed. Finally, we develop a numerical method for generic two-field models of inflation, which allows us to go beyond the slow-roll approximation and confirms our analytic results for double inflation.

We investigate whether enhanced gravitational scattering on small scales (< 0.1mm), which becomes possible in models with large extra dimensions, can establish statistical equilibrium between different particle species in the early Universe. We calculate the classical relativistic energy transfer rate for two species with a large ratio between their masses for a general elastic scattering cross section. Although the classical calculation suggests that ultra-light WIMPs (e.g., axions) can be thermalized by gravitational scattering, such interactions are considerably less efficient once quantum effects are taken into account on scales below the Compton wavelength. However the energy transfer rate in models with several extra dimensions may still be sensitive to trans-Planckian physics.

We review the theory of inflation with single and multiple fields paying particular attention to the dynamics of adiabatic and entropy/isocurvature perturbations which provide the primary means of testing inflationary models. We review the theory and phenomenology of reheating and preheating after inflation providing a unified discussion of both the gravitational and nongravitational features of multi-field inflation. In addition we cover inflation in theories with extra dimensions and models such as the curvaton scenario and modulated reheating which provide alternative ways of generating large-scale density perturbations. Finally we discuss the interesting observational implications that can result from adiabatic-isocurvature correlations and non-Gaussianity.

Brane-world models, where observers are restricted to a brane in a higher-dimensional spacetime, offer a novel perspective on cosmology. I discuss some approaches to cosmology in extra dimensions and some interesting aspects of gravity and cosmology in brane-world models.

We show how the primordial bispectrum of density perturbations from inflation may be characterised in terms of manifestly gauge-invariant cosmological perturbations at second order. The primordial metric perturbation, zeta, describing the perturbed expansion of uniform-density hypersurfaces on large scales is related to scalar field perturbations on unperturbed (spatially-flat) hypersurfaces at first- and second-order. The bispectrum of the metric perturbation is thus composed of (i) a local contribution due to the second-order gauge-transformation, and (ii) the instrinsic bispectrum of the field perturbations on spatially flat hypersurfaces. We generalise previous results to allow for scale-dependence of the scalar field power spectra and correlations that can develop between fields on super-Hubble scales.

Vacuum fluctuations in the inflaton field driving chaotic inflation with a quadratic potential give a red spectrum of primordial density perturbations, n=0.97. However angular fluctuations in an O(N)-symmetric quadratic potential have a very nearly scale-invariant spectrum, n=0.9998. We investigate the possibility that these isocurvature field perturbations could give the dominant contribution to the primordial density perturbation after inflation.

We propose a geometrical model of brane inflation where inflation is driven by the flux generated by opposing brane charges and terminated by the collision of the branes, with charge annihilation. We assume the collision process is completely inelastic and the kinetic energy is transformed into the thermal energy after collision. Thereafter the two branes coalesce together and behave as a single brane universe with zero effective cosmological constant. In the Einstein frame, the 4-dimensional effective theory changes abruptly at the collision point. Therefore, our inflationary model is necessarily 5-dimensional in nature. As the collision process has no singularity in 5-dimensional gravity, we can follow the evolution of fluctuations during the whole history of the universe. It turns out that the radion field fluctuations have a steeply tilted, red spectrum, while the primordial gravitational waves have a flat spectrum. Instead, primordial density perturbations could be generated by a curvaton mechanism.

We propose a model which produces dark energy perturbations large enough to explain the lack of power seen at the quadrupole scale in the cosmic microwave background. If the dark energy is frozen from horizon exit during inflation until dark energy domination, then it is not possible to have perturbations in the dark energy which are large enough. We propose using a tachyonic amplification mechanism to overcome this. The dark energy is taken to be a complex scalar field, where the radial field has a Mexican hat potential. During inflation, the radial component is trapped near the maximum of its potential. At the end of inflation, it rolls down to the minimum. The dark energy today is taken to be a pseudo-Nambu-Goldstone boson. The perturbations generated during inflation are amplified by the rolling of the radial field. We also examine the use of the variable decay mechanism in order to generate an anti-correlation between the dark energy perturbations and the curvature perturbation. We show that using this mechanism then constrains the properties of the dark energy and its evolution from redshift one until today.

We develop a gauge-invariant formalism for the study of density perturbations in a Friedmann-Robertson-Walker universe with multiple interacting fluids and/or scalar fields. We show how N scalar fields may be described by N kinetic fluids (with maximally stiff equation of state) interacting with a non-dynamical potential (with vacuum equation of state). We split generic perturbations into adiabatic and entropic parts, and give the coupled first-order evolution equations on all scales, including energy and momentum exchange. We identify the non-adiabatic effects on large scales, and define adiabatic initial conditions in the presence of multiple fluids and fields.

We investigate the scalar metric perturbations about a de Sitter brane universe in a 5-dimensional anti de Sitter bulk. We compare the master-variable formalism, describing metric perturbations in a 5-dimensional longitudinal gauge, with results in a Gaussian normal gauge. For a vacuum brane (with constant brane tension) there is a continuum of normalizable Kaluza-Klein modes, with m>3H/2, which remain in the vacuum state. A light radion mode, with m=\sqrt{2}H, satisfies the boundary conditions for two branes but is not normalizable in the single-brane case. When matter is introduced (as a test field) on the brane, this mode, together with the zero-mode and an infinite ladder of discrete tachyonic modes, become normalizable. However, the boundary condition requires the self-consistent 4-dimensional evolution of scalar field perturbations on the brane and the dangerous growing modes are not excited. These normalizable discrete modes introduce corrections at first-order to the scalar field perturbations computed in a slow-roll expansion. On super-Hubble scales, the correction is smaller than slow-roll corrections to the de Sitter background. However on small scales the corrections can become significant.

We present a detailed study of a simple scalar field model that yields non-singular cosmological solutions. We study both the qualitative dynamics of the homogeneous and isotropic background and the evolution of inhomogeneous linear perturbations. We calculate the spectrum of perturbations generated on super-Hubble scales during the collapse phase from initial vacuum fluctuations on small scales and then evolve these numerically through the bounce. We show there is a gauge in which perturbations remain well-defined and small throughout the bounce, even though perturbations in other commonly used gauges become large or ill-defined. We show that the comoving curvature perturbation calculated during the collapse phase provides a good estimate of the resulting large scale adiabatic perturbation in the expanding phase while the Bardeen metric potential is dominated by what becomes a decaying mode after the bounce. We show that a power-law collapse phase with scale factor proportional to $(-t)^{2/3}$ can yield a scale-invariant spectrum of adiabatic scalar perturbations in the expanding phase, but the amplitude of tensor perturbations places important constraints on the allowed initial conditions.

We calculate the amplitude of gravitational waves produced by inflation on a de Sitter brane embedded in five-dimensional anti-de Sitter bulk spacetime, extending previous calculations in Randall-Sundrum type cosmology to include the effect of induced gravity corrections on the brane. These corrections arise via a term in the brane action that is proportional to the brane Ricci scalar. We find that, as in the Randall-Sundrum case, there is a mass gap between the discrete zero-mode and a continuum of massive bulk modes, which are too heavy to be excited during inflation. We give the normalization of the zero-mode as a function of the Hubble rate on the brane and are thus able to calculate the high energy correction to the spectrum of gravitational wave (tensor) modes excited on large scales during inflation from initial vacuum fluctuations on small scales. We also calculate the amplitude of density (scalar) perturbations expected due to inflaton fluctuations on the brane, and show that the usual four-dimensional consistency relation for the tensor/scalar ratio remains valid for brane inflation with induced gravity corrections.

We present a method for constructing gauge-invariant cosmological perturbations which are gauge-invariant up to second order. As an example we give the gauge-invariant definition of the second-order curvature perturbation on uniform density hypersurfaces. Using only the energy conservation equation we show that this curvature perturbation is conserved at second order on large scales for adiabatic perturbations.

We study the evolution of the cosmological perturbations after inflation in curvaton models where the non-relativistic curvaton decays into both radiation and a cold dark matter component. We calculate the primordial curvature and correlated isocurvature perturbations inherited by the radiation and cold dark matter after the curvaton has decayed. We give the transfer coefficient in terms of the initial curvaton density relative to the curvaton decay rate.

We investigate the evolution of gravitational wave perturbations about a brane cosmology embedded in a five-dimensional anti-de Sitter bulk. During slow-roll inflation in a Randall-Sundrum brane-world, the zero mode of the 5-dimensional graviton is generated, while the massive modes remain in their vacuum state. When the zero mode re-enters the Hubble radius during radiation domination, massive modes are generated. We show that modes decouple in the low-energy/near-brane limit and develop perturbative techniques to calculate the mode-mixing at finite energy.

I discuss how parameters describing inflation in the very early universe may be related to primordial perturbation spectra. Precision observations of anisotropies in the cosmic microwave background (CMB) such as those provided by the WMAP satellite offer an unprecedented window onto the physics of the very early universe. To theorists exploring speculative models of physics at high energies, the CMB offers us the chance to put our ideas to the test.

A conserved cosmological perturbation is associated with each quantity whose local evolution is determined entirely by the local expansion of the Universe. It may be defined as the appropriately normalised perturbation of the quantity, defined using a slicing of spacetime such that the expansion between slices is spatially homogeneous. To first order, on super-horizon scales, the slicing with unperturbed intrinsic curvature has this property. A general construction is given for conserved quantities, yielding the curvature perturbation $\zeta$ as well as more recently-considered conserved perturbations. The construction may be extended to higher orders in perturbation theory and even into the non-perturbative regime.

We discuss the residual isocurvature perturbations, fully-correlated with the curvature perturbation, that are automatic in the curvaton scenario if curvaton decay is sufficiently late. We contrast these residual isocurvature perturbations with the generally un-correlated `intrinsic' isocurvature perturbation generated by an additional field such as the axion. We present a general formula for the residual isocurvature perturbations, referring only to the generation of the relevant quantity (Cold Dark Matter, baryon number or lepton number) in an unperturbed universe. Specific formulas for the residual isocurvature CDM perturbation are given, for most of the commonly-considered CDM candidates.

We present a gauge-invariant formalism to study the evolution of curvature perturbations in a Friedmann-Robertson-Walker universe filled by multiple interacting fluids. We resolve arbitrary perturbations into adiabatic and entropy components and derive their coupled evolution equations. We demonstrate that perturbations obeying a generalised adiabatic condition remain adiabatic in the large-scale limit, even when one includes energy transfer between fluids. As a specific application we study the recently proposed curvaton model, in which the curvaton decays into radiation. We use the coupled evolution equations to show how an initial isocurvature perturbation in the curvaton gives rise to an adiabatic curvature perturbation after the curvaton decays.

We analyse the primordial density perturbation when it is generated by a `curvaton' field different from the inflaton. In some cases this perturbation may have large isocurvature components, fully correlated or anti-correlated with the adiabatic component. It may also have a significant non-Gaussian component. All of these effects are calculated in a form which will enable direct comparison with current and forthcoming observational data.

We study adiabatic and isocurvature perturbation spectra produced by a period of cosmological inflation driven by two scalar fields. We show that there exists a model-independent consistency condition for all two-field models of slow-roll inflation, despite allowing for model-dependent linear processing of curvature and isocurvature perturbations during and after inflation on super-horizon scales. The scale-dependence of all spectra are determined solely in terms of slow-roll parameters during inflation and the dimensionless cross-correlation between curvature and isocurvature perturbations. We present additional model-dependent consistency relations that may be derived in specific two-field models, such as the curvaton scenario.

We compare the latest cosmic microwave background data with theoretical predictions including correlated adiabatic and CDM isocurvature perturbations with a simple power-law dependence. We find that there is a degeneracy between the amplitude of correlated isocurvature perturbations and the spectral tilt. A negative (red) tilt is found to be compatible with a larger isocurvature contribution. Estimates of the baryon and CDM densities are found to be almost independent of the isocurvature amplitude. The main result is that current microwave background data do not exclude a dominant contribution from CDM isocurvature fluctuations on large scales.

We consider the collision of self-gravitating n-branes in a (n+2)-dimensional spacetime. We show that there is a geometrical constraint which can be expressed as a simple sum rule for angles characterizing Lorentz boosts between branes and the intervening spacetime regions. This constraint can then be re-interpreted as either energy or momentum conservation at the collision.

I discuss cosmological models either derived from, or inspired by, string theory or M-theory. In particular I discuss solutions in the low-energy effective theory and the role of the dilaton, moduli and antisymmetric form fields in the dimensionally reduced effective action. The pre big bang model is an attempt to use cosmological solutions to make observational predictions. I then discuss the effective theory of gravity found in recent brane-world models where we live on a 3-brane embedded in a five-dimensional spacetime and how the study of cosmological perturbations may enable us to test these ideas.

We present a mechanism for the origin of the large-scale curvature perturbation in our Universe by the late decay of a massive scalar field, the curvaton. The curvaton is light during a period of cosmological inflation, when it acquires a perturbation with an almost scale-invariant spectrum. This corresponds initially to an isocurvature density perturbation, which generates the curvature perturbation after inflation when the curvaton density becomes a significant fraction of the total. The isocurvature density perturbation disappears if the curvaton completely decays into thermalised radiation. Any residual isocurvature perturbation is 100% correlated with the curvature. The same mechanism can also generate the curvature perturbation in pre big bang/ekpyrotic models, provided that the curvaton has a suitable non-canonical kinetic term so as to generate a flat spectrum.

I review the standard analysis of adiabatic scalar and tensor perturbations produced by slow-roll inflation driven by a single scalar field, before going on to discuss recent work on the role of non-adiabatic modes during and after inflation. Isocurvature perturbations correlated with adiabatic modes may be produced during multi-field inflation and would give valuable information about the physics of the early universe. A significant contribution from correlated isocurvature perturbations is not ruled out by present data, but should be either detected or ruled out by future observations.

We study cosmological perturbations in a brane-world scenario where the matter fields live on a four-dimensional brane and gravity propagates in the five-dimensional bulk. We present the equations of motion in an arbitrary gauge for metric perturbations in the bulk and matter perturbations on the brane. Gauge-invariant perturbations are then constructed corresponding to perturbations in longitudinal and Gaussian normal gauges. Longitudinal gauge metric perturbations may be directly derived from three master variables (separately describing scalar, vector and tensor metric perturbations) which obey five-dimensional wave-equations. Gaussian normal gauge perturbations are directly related to the induced metric perturbations on the brane with the additional bulk degrees of freedom interpreted as an effective Weyl energy-momentum tensor on the brane. We construct gauge-invariant perturbations describing the effective density, momentum and pressures of this Weyl fluid at the brane and throughout the bulk. We show that there exist gauge-invariant curvature perturbations on the brane and in the bulk that are conserved on large-scales when three-dimensional spatial gradients are negligible.

We show that there exists a simple mechanism which can enhance the amplitude of curvature perturbations on superhorizon scales during inflation, relative to their amplitude at horizon crossing. The enhancement may occur even in a single-field inflaton model, and occurs if the quantity $a\dot\phi/H$ becomes sufficiently small, as compared to its value at horizon crossing, for some time interval during inflation. We give a criterion for this enhancement in general single-field inflation models.

In brane-world cosmologies of Randall-Sundrum type, we show that evolution of large-scale curvature perturbations may be determined on the brane, without solving the bulk perturbation equations. The influence of the bulk gravitational field on the brane is felt through a projected Weyl tensor which behaves effectively like an imperfect radiation fluid with anisotropic stress. We define curvature perturbations on uniform density surfaces for both the matter and Weyl fluids, and show that their evolution on large scales follows directly from the energy conservation equations for each fluid. The total curvature perturbation is not necessarily constant for adiabatic matter perturbations, but can change due to the Weyl entropy perturbation. To relate this curvature perturbation to the longitudinal gauge metric potentials requires knowledge of the Weyl anisotropic stress which is not determined by the equations on the brane. We discuss the implications for large-angle anisotropies on the cosmic microwave background sky.

We study adiabatic (curvature) and entropy (isocurvature) perturbations produced during a period of cosmological inflation that is driven by multiple scalar fields with an arbitrary interaction potential. A local rotation in field space is performed to separate out the adiabatic and entropy modes. The resulting field equations show explicitly how on large scales entropy perturbations can source adiabatic perturbations if the background solution follows a curved trajectory in field space, and how adiabatic perturbations cannot source entropy perturbations in the long-wavelength limit. It is the effective mass of the entropy field that determines the amplitude of entropy perturbations during inflation. We present two applications of the equations. First, we show why one in general expects the adiabatic and entropy perturbations to be correlated at the end of inflation, and calculate the cross-correlation in the context of a double inflation model with two non-interacting fields. Second, we consider two-field preheating after inflation, examining conditions under which entropy perturbations can alter the large-scale curvature perturbation and showing how our new formalism has advantages in numerical stability when the background solution follows a non-trivial trajectory in field space.

We discuss the evolution of linear perturbations about a Friedmann-Robertson-Walker background metric, using only the local conservation of energy-momentum. We show that on sufficiently large scales the curvature perturbation on spatial hypersurfaces of uniform-density is constant when the non-adiabatic pressure perturbation is negligible. We clarify the conditions under which super-horizon curvature perturbations may vary, using preheating as an example.

We discuss the evolution of gravitational waves in a brane-world cosmology embedded in five-dimensional anti-de Sitter spacetime. We show that during slow-roll inflation, modelled as a period of quasi-de Sitter expansion on the brane, there is a discrete normalizable massless graviton mode. There is a mass gap due to the expansion, above which there is a continuum of massive modes. Only the massless mode is `light' compared with the Hubble scale during inflation, leading to the production of classical perturbations on large scales from vacuum fluctuations on small scales. We calculate the amplitude of these fluctuations at horizon-crossing and show that the standard four-dimensional result is recovered at low energies, but the amplitude of the perturbations is enhanced at high energies.

We discuss the evolution of linear perturbations about a Friedmann-Robertson-Walker background metric, using only the local conservation of energy-momentum. We show that on sufficiently large scales the curvature perturbation on spatial hypersurfaces of uniform-density is conserved when the non-adiabatic pressure perturbation is negligible. This is the first time that this result has been demonstrated independently of the gravitational field equations. A physical picture of long-wavelength perturbations as being composed of separate Robertson-Walker universes gives a simple understanding of the possible evolution of the curvature perturbation, in particular clarifying the conditions under which super-horizon curvature perturbations may vary.

We consider slow-roll inflation in the context of recently proposed four-dimensional effective gravity induced on the world-volume of a three-brane in five-dimensional Einstein gravity. We find significant modifications of the simplest chaotic inflationary scenario when the five-dimensional Planck scale is below about 10^{17} GeV. We use the comoving curvature perturbation, which remains constant on super-Hubble scales, in order to calculate the spectrum of adiabatic density perturbations generated. Modifications to the Friedmann constraint equation lead to a faster Hubble expansion at high energies and a more strongly damped evolution of the scalar field. This assists slow-roll, enhances the amount of inflation obtained in any given model, and drives the perturbations towards an exactly scale-invariant Harrison-Zel'dovich spectrum. In chaotic inflation driven by a massive scalar field we show that inflation can occur at field values far below the four-dimensional Planck scale, though above the five-dimensional fundamental scale.

Aspects of superstring cosmology are reviewed with an emphasis on the cosmological implications of duality symmetries in the theory. The string effective actions are summarized and toroidal compactification to four dimensions reviewed. Global symmetries that arise in the compactification are discussed and the duality relationships between the string effective actions are then highlighted. Higher-dimensional Kasner cosmologies are presented and interpreted in both string and Einstein frames, and then given in dimensionally reduced forms. String cosmologies containing both non-trivial Neveu-Schwarz/Neveu-Schwarz and Ramond-Ramond fields are derived by employing the global symmetries of the effective actions. Anisotropic and inhomogeneous cosmologies in four-dimensions are also developed. The review concludes with a detailed analysis of the pre-big bang inflationary scenario. The generation of primordial spectra of cosmological perturbations in such a scenario is discussed. Possible future directions offered in the Horava-Witten theory are outlined.

It has recently been claimed by Bassett et al that preheating after inflation may affect the amplitude of curvature perturbations on large scales, undermining the usual inflationary prediction. We analyze the simplest model, and confirm the results of Jedamzik and Sigl and of Ivanov that in linear perturbation theory the effect is negligible. However the dominant effect is second-order in the field perturbation and we show that this too is negligible, and hence conclude that preheating has no significant influence on large-scale perturbations in this model. We briefly discuss the likelihood of an effect in other models.

We present self-similar cosmological solutions for a barotropic fluid plus scalar field with Brans-Dicke-type coupling to the spacetime curvature and an arbitrary power-law potential energy. We identify all the fixed points in the autonomous phase-plane, including a scaling solution where the fluid density scales with the scalar field's kinetic and potential energy. This is related by a conformal transformation to a scaling solution for a scalar field with exponential potential minimally coupled to the spacetime curvature, but non-minimally coupled to the barotropic fluid. Radiation is automatically decoupled from the scalar field, but energy transfer between the field and non-relativistic dark matter can lead to a change to an accelerated expansion at late times in the Einstein frame. The scalar field density can mimic a cosmological constant even for steep potentials in the strong coupling limit.

I show that cosmological perturbation spectra produced from quantum fluctuations in massless or self-interacting scalar fields during an inflationary era remain invariant under a two parameter family of transformations of the homogeneous background fields. This relates slow-roll inflation models to solutions which may be far from the usual slow-roll limit. For example, a scale-invariant spectrum of perturbations in a minimally coupled, massless field can be produced by an exponential expansion with $a\propto e^{Ht}$, or by a collapsing universe with $a\propto (-t)^{2/3}$.

We investigate the dynamics of the recently proposed model of assisted inflation. In this model an arbitrary number of scalar fields with exponential potentials evolve towards an inflationary scaling solution, even if each of the individual potentials is too steep to support inflation on its own. By choosing an appropriate rotation in field space we can write down explicitly the potential for the weighted mean field along the scaling solution and for fields orthogonal to it. This demonstrates that the potential has a global minimum along the scaling solution. We show that the potential close to this attractor in the rotated field space is analogous to a hybrid inflation model, but with the vacuum energy having an exponential dependence upon a dilaton field. We present analytic solutions describing homogeneous and inhomogeneous perturbations about the attractor solution without resorting to slow-roll approximations. We discuss the curvature and isocurvature perturbation spectra produced from vacuum fluctuations during assisted inflation.

We present a phase-plane analysis of cosmologies containing a barotropic fluid with equation of state $p_\gamma = (\gamma-1) \rho_\gamma$, plus a scalar field $\phi$ with an exponential potential $V \propto \exp(-\lambda \kappa \phi)$ where $\kappa^2 = 8\pi G$. In addition to the well-known inflationary solutions for $\lambda^2 < 2$, there exist scaling solutions when $\lambda^2 > 3\gamma$ in which the scalar field energy density tracks that of the barotropic fluid (which for example might be radiation or dust). We show that the scaling solutions are the unique late-time attractors whenever they exist. The fluid-dominated solutions, where $V(\phi)/\rho_\gamma \to 0$ at late times, are always unstable (except for the cosmological constant case $\gamma = 0$). The relative energy density of the fluid and scalar field depends on the steepness of the exponential potential, which is constrained by nucleosynthesis to $\lambda^2 > 20$. We show that standard inflation models are unable to solve this `relic density' problem.

We show how gauge-invariant cosmological perturbations may be constructed by an unambiguous choice of hypersurface-orthogonal time-like vector field (i.e., time-slicing). This may be defined either in terms of the metric quantities such as curvature or shear, or using some matter field. As an example, we show how linear perturbations in the covariant fluid-flow approach can then be presented in an explicitly gauge-invariant form in the coordinate based formalism.

Pre--big bang models of inflation based on string cosmology produce a stochastic gravitational wave background whose spectrum grows with decreasing wavelength, and which may be detectable using interferometers such as LIGO. We point out that the gravitational wave spectrum is closely tied to the density perturbation spectrum, and that the condition for producing observable gravitational waves is very similar to that for producing an observable density of primordial black holes. Detection of both would provide strong support to the string cosmology scenario.

We discuss the appropriate normalization of modes required to generate a homogeneous random field in an open Friedmann-Robertson-Walker universe. We consider scalar random fields and certain tensor random fields that can be obtained by covariantly differentiating a scalar. Modes of interest fall into three categories: the familiar sub-curvature modes, the more recently discussed super-curvature modes, and a set of discrete modes with positive eigenvalues which can be used to generate homogeneous tensor random fields even though the underlying scalar field is not homogeneous. A particular example of the last case which has been discussed in the literature is the bubble wall fluctuation in open inflationary universes.

We discuss metric perturbations produced during a period of inflation in the early universe where two scalar fields evolve. The final scalar perturbation spectrum can be calculated in terms of the perturbed expansion along neighbouring trajectories in field-space. In the usual single field case this is fixed by the values of the fields at horizon-crossing, but in the presence of more than one field there is no longer a unique slow-roll trajectory. The presence of entropy as well as adiabatic fluctuations means that the super-horizon-sized metric perturbation $\zeta$ may no longer be conserved and the evolution must be integrated along the whole of the subsequent trajectory. In general there is an inequality between the ratio of tensor to scalar perturbations and the tilt of the gravitational wave spectrum, which becomes an equality when only adiabatic perturbations are possible and $\zeta$ is conserved.

We study the amplitude and spectral tilt of density perturbations in the simplest hybrid inflation models. We give an exact expression for the amplitude of quantum fluctuations on all scales in the limit where we can neglect the backreaction on the metric. This is a very good approximation for values of the inflaton field well below the Planck scale and our results remain valid far from the usual massless limit. We confirm that the primordial density spectrum in this model has a constant spectral index n>1 over all observable scales. For the small values of the tilt (n<1.4) required by observations, the results remain close to those obtained using the quasi-massless approximation.

We investigate the recently proposed hybrid inflation models with two stages of inflation. We show that quantum fluctuations at the time corresponding to the phase transition between the two inflationary stages can trigger the formation of a large number of inflating topological defects. In order to study density perturbations in these models we develop a new method to calculate density perturbations in a system of two scalar fields. We show that density perturbations in hybrid inflation models of the new type can be very large on the scale corresponding to the phase transition. The resulting density inhomogeneities lead to a copious production of black holes. This could be an argument against hybrid inflation models with two stages of inflation. However, we find a class of models where this problem can be easily avoided. The number of black holes produced in these models can be made extremely small, but in general it could be sufficiently large to have important cosmological and astrophysical implications. In particular, for certain values of parameters these black holes may constitute the dark matter in the universe. It is also possible to have hybrid models with two stages of inflation where the black hole production is not suppressed, but where the typical masses of the black holes are very small. Such models lead to a completely different thermal history of the universe, where post-inflationary reheating occurs via black hole evaporation.

We study the metric perturbations produced during inflation in models with two scalar fields evolving simultaneously. In particular, we emphasize how the large-scale curvature perturbation $\zeta$ on fixed energy density hypersurfaces may not be conserved in general for multiple field inflation due to the presence of entropy as well as adiabatic fluctuations. We show that the usual method of solving the linearized perturbation equations is equivalent to the recently proposed analysis of Sasaki and Stewart in terms of the perturbed expansion along neighboring trajectories in field-space. In the case of a separable potential it is possible to compute in the slow-roll approximation the spectrum of density perturbations and gravitational waves at the end of inflation. In general there is an inequality between the ratio of tensor to scalar perturbations and the tilt of the gravitational wave spectrum, which becomes an equality when only adiabatic perturbations are possible and $\zeta$ is conserved.

The Grishchuk--Zel'dovich effect is the contribution to the microwave background anisotropy from an extremely large scale adiabatic density perturbation, on the standard hypothesis that this perturbation is a typical realization of a homogeneous Gaussian random field. We analyze this effect in open universes, corresponding to density parameter $\Omega_0<1$ with no cosmological constant, and concentrate on the recently discussed super-curvature modes. The effect is present in all of the low multipoles of the anisotropy, in contrast with the $\Omega_0=1$ case where only the quadrupole receives a contribution. However, for no value of $\Omega_0$ can a very large scale perturbation generate a spectrum capable of matching observations across a wide range of multipoles. We evaluate the magnitude of the effect coming from a given wavenumber as a function of the magnitude of the density perturbation, conveniently specified by the mean-square curvature perturbation. From the absence of the effect at the observed level, we find that for $0.25\leq\Omega_0\leq 0.8$, a curvature perturbation of order unity is permitted only for inverse wavenumbers more than one thousand times the size of the observable universe. As $\Omega_0$ tends to one, the constraint weakens to the flat space result that the inverse wavenumber be more than a hundred times the size of the observable universe, whereas for $\Omega_0 < 0.25$ it becomes stronger. We explain the physical meaning of these results, by relating them to the correlation length of the perturbation. Finally, in an Appendix we consider the dipole anisotropy and show that it always leads to weaker constraints.

We show how observations of the perturbation spectra produced during inflation may be used to constrain the parameters of general scalar-tensor theories of gravity, which include both an inflaton and dilaton field. An interesting feature of these models is the possibility that the curvature perturbations on super-horizon scales may not be constant due to non-adiabatic perturbations of the two fields. Within a given model, the tilt and relative amplitude of the scalar and tensor perturbation spectra gives constraints on the parameters of the gravity theory, which may be comparable with those from primordial nucleosynthesis and post-Newtonian experiments.

Quantum fluctuations of scalar fields during inflation could determine the very large-scale structure of the universe. In the case of general scalar-tensor gravity theories these fluctuations lead to the diffusion of fundamental constants like the Planck mass and the effective Brans--Dicke parameter, $\omega$. In the particular case of Brans--Dicke gravity, where $\omega$ is constant, this leads to runaway solutions with infinitely large values of the Planck mass. However, in a theory with variable $\omega$ we find stationary probability distributions with a finite value of the Planck mass peaked at exponentially large values of $\omega$ after inflation. We conclude that general relativity is an attractor during the quantum diffusion of the fields.

I give a general formulation of the constraints on models of inflation ended by a first order phase transition arising from the requirement that they do not produce too many large (observable) true vacuum voids -- the `big bubble problem'. It is shown that this constraint can be satisfied by a simple model in Einstein gravity -- a variant of `hybrid' or `false vacuum' inflation. (Talk presented at `Birth of the Universe' workshop Rome, May 1994)

We investigate chaotic inflation models with two scalar fields, such that one field (the inflaton) rolls while the other is trapped in a false vacuum state. The false vacuum becomes unstable when the inflaton field falls below some critical value, and a first or second order transition to the true vacuum ensues. Particular attention is paid to Linde's second-order `Hybrid Inflation'; with the false vacuum dominating, inflation differs from the usual true vacuum case both in its cosmology and in its relation to particle physics. The spectral index of the adiabatic density perturbation can be very close to 1, or it can be around ten percent higher. The energy scale at the end of inflation can be anywhere between $10^{16}$\,GeV and $10^{11}$\,GeV, though reheating is prompt so the reheat temperature can't be far below $10^{11}\,$GeV. Topological defects are almost inevitably produced at the end of inflation, and if the inflationary energy scale is near its upper limit they can have significant effects. Because false vacuum inflation occurs with the inflaton field far below the Planck scale, it is easier to implement in the context of supergravity than standard chaotic inflation. That the inflaton mass is small compared with the inflationary Hubble parameter is still a problem for generic supergravity theories, but remarkably this can be avoided in a natural way for a class of supergravity models which follow from orbifold compactification of superstrings. This opens up the prospect of a truly realistic, superstring