It is commonly assumed that the energy density of the Universe was dominated by radiation between reheating after inflation and the onset of matter domination 54,000 years later. While the abundance of light elements indicates that the Universe was radiation dominated during Big Bang Nucleosynthesis (BBN), there is scant evidence that the Universe was radiation dominated prior to BBN. It is therefore possible that the cosmological history was more complicated, with deviations from the standard radiation domination during the earliest epochs. Indeed, several interesting proposals regarding various topics such as the generation of dark matter, matter-antimatter asymmetry, gravitational waves, primordial black holes, or microhalos during a nonstandard expansion phase have been recently made. In this paper, we review various possible causes and consequences of deviations from radiation domination in the early Universe - taking place either before or after BBN - and the constraints on them, as they have been discussed in the literature during the recent years.

We have calculated the detectable merger rate of primordial black holes, as a function of the redshift, as well as the binary's mass ratio, total mass and chirp mass (observables that have not previously been explored in great detail for PBHs). We consider both the current and design sensitivity of LIGO and five different primordial black hole mass functions, as well as showing a comparison to a predicted astrophysical black hole merger rate. We show that the empirical preference for nearly equal-mass binaries in current LIGO/Virgo data can be consistent with a PBH hypothesis once observational selection effects are taken into account. However, current data do exclude some PBH mass distributions, and future data may be able to rule out the possibility that all observed BH mergers had a primordial origin.

The formation and abundance of primordial black holes (PBHs) arising from the curvature perturbation $\zeta$ is studied. The non-linear relation between $\zeta$ and the density contrast $\delta$ means that, even when $\zeta$ has an exactly Gaussian distribution, significant non-Gaussianities affecting PBH formation must be considered. Numerical simulations are used to investigate the critical value and the mass of PBHs which form, and peaks theory is used to calculate the mass fraction of the universe collapsing to form PBHs at the time of formation. A formalism to calculate the total present day PBH abundance and mass function is also derived. It is found that the abundance of PBHs is very sensitive to the non-linear effects, and that the power spectrum $\mathcal{P}_\zeta$ must be a factor of $\sim2$ larger to produce the same number of PBHs as the linear model (where the exact value depends on the critical value for a region to collapse and form a PBH). This also means that the derived constraints on the small-scale power spectrum from constraints on the abundance of PBHs are weaker by the same factor.

If the primordial curvature perturbation followed a Gaussian distribution, primordial black holes (PBHs) will be Poisson distributed with no additional clustering. We consider local non-Gaussianity and its impact on the initial PBH clustering and mass function due to mode coupling between long and short wavelength modes. We show that even a small amount of non-Gaussianity results in a significant enhancement on the PBH initial clustering and subsequent merger rate and that the PBH mass function shifts to higher mass PBHs. However, as the clustering becomes strong, the local number density of PBHs becomes large, leading to a large theoretical uncertainty in the merger rate.

We recently showed that postulated ultracompact minihalos with a steep density profile do not form in realistic simulations with enhanced initial perturbations. In this paper we assume that a small fraction of the dark matter consists of primordial black holes (PBHs) and simulate the formation of structures around them. We find that in this scenario halos with steep density profiles do form, consistent with theoretical predictions. If the rest of the dark matter consists of weakly interacting massive particles (WIMPs), we also show that WIMPs in the dense innermost part of these halos would annihilate and produce a detectable gamma-ray signal. The non-detection of this signal implies that PBHs make up at most one billionth of the dark matter, provided that their mass is greater than one millionth of the mass of the Sun. Similarly, a detection of PBHs would imply that the remaining dark matter could not be WIMPs.

We derive analytic bounds on the shape of the primordial power spectrum in the context of single-field inflation. In particular, the steepest possible growth has a spectral index of $n_s - 1 = 4$ once transients have died down. Its primary implication is that any constraint on the power spectrum at a particular scale can be extrapolated to an upper bound over an extended range of scales. This is important for models which generate relics due to an enhanced amplitude of the primordial scalar perturbations, such as primordial black holes. In order to generate them, the power spectrum needs to grow many orders of magnitude larger than its observed value on CMB scales - typically achieved through a phase of ultra slow-roll inflation - and is thus subject to additional constraints at small scales. We plot all relevant constraints including CMB spectral distortions and gravitational waves sourced by scalar perturbations at second order. We show how this limits the allowed mass of PBHs, especially for the large masses of interest following recent detections by LIGO and prospects for constraining them further with future observations. We show that any transition from approximately constant $\epsilon$ slow-roll inflation to a phase where the power spectrum rapidly rises necessarily implies an intervening dip in power. We also show how to reconstruct a potential that can reproduce an arbitrary time-varying $\epsilon$, offering a complementary perspective on how ultra slow-roll can be achieved.

Following the pioneering observations with COBE in the early 1990s, studies of the cosmic microwave background (CMB) have focused on temperature and polarization anisotropies. CMB spectral distortions - tiny departures of the CMB energy spectrum from that of a perfect blackbody - provide a second, independent probe of fundamental physics, with a reach deep into the primordial Universe. The theoretical foundation of spectral distortions has seen major advances in recent years, which highlight the immense potential of this emerging field. Spectral distortions probe a fundamental property of the Universe - its thermal history - thereby providing additional insight into processes within the cosmological standard model (CSM) as well as new physics beyond. Spectral distortions are an important tool for understanding inflation and the nature of dark matter. They shed new light on the physics of recombination and reionization, both prominent stages in the evolution of our Universe, and furnish critical information on baryonic feedback processes, in addition to probing primordial correlation functions at scales inaccessible to other tracers. In principle the range of signals is vast: many orders of magnitude of discovery space could be explored by detailed observations of the CMB energy spectrum. Several CSM signals are predicted and provide clear experimental targets, some of which are already observable with present-day technology. Confirmation of these signals would extend the reach of the CSM by orders of magnitude in physical scale as the Universe evolves from the initial stages to its present form. The absence of these signals would pose a huge theoretical challenge, immediately pointing to new physics.

The expansion of the Universe is understood to have accelerated during two epochs: in its very first moments during a period of Inflation and much more recently, at $z < 1$, when Dark Energy is hypothesized to drive cosmic acceleration. The undiscovered mechanisms behind these two epochs represent some of the most important open problems in fundamental physics. The large cosmological volume at $2 < z < 5$, together with the ability to efficiently target high-$z$ galaxies with known techniques, enables large gains in the study of Inflation and Dark Energy. A future spectroscopic survey can test the Gaussianity of the initial conditions up to a factor of ~50 better than our current bounds, crossing the crucial theoretical threshold of $\sigma(f_{NL}^{\rm local})$ of order unity that separates single field and multi-field models. Simultaneously, it can measure the fraction of Dark Energy at the percent level up to $z = 5$, thus serving as an unprecedented test of the standard model and opening up a tremendous discovery space.

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.

The LIGO discoveries have rekindled suggestions that primordial black holes (BHs) may constitute part to all of the dark matter (DM) in the Universe. Such suggestions came from 1) the observed merger rate of the BHs, 2) their unusual masses, 3) their low/zero spins, and 4) also from the independently uncovered cosmic infrared background (CIB) fluctuations signal of high amplitude and coherence with unresolved cosmic X-ray background (CXB). Here we summarize the prospects to resolve this important issue with electromagnetic observations using the instruments and tools expected in the 2020's. These prospects appear promising to make significant, and potentially critical, advances. We demonstrate that in the next decade, new space- and ground-borne electromagnetic instruments, combined with concurrent theoretical efforts, should shed critical light on the long-considered link between primordial BHs and DM. Specifically the new data and methodologies under this program will involve: I) Probing with high precision the spatial spectrum of source-subtracted CIB with Euclid and WFIRST, and its coherence with unresolved cosmic X-ray background using eROSITA and Athena, II) Advanced searches for microlensing of Galactic stars by the intervening Galactic Halo BHs with OGLE, Gaia, LSST and WFIRST, III) Supernovae (SNe) lensing in the upcoming surveys with WFIRST, LSST and also potentially with Euclid and JWST, IV) Advanced theoretical work to understand the details of PBH accretion and evolution and their influence on cosmic microwave background (CMB) anisotropies in light of the next generation CMB experiments, V) Better new samples and theoretical understanding involving stability and properties of ultra faint dwarf galaxies, pulsar timing, and cosmological quasar lensing.

The grand challenges of contemporary fundamental physics---dark matter, dark energy, vacuum energy, inflation and early universe cosmology, singularities and the hierarchy problem---all involve gravity as a key component. And of all gravitational phenomena, black holes stand out in their elegant simplicity, while harbouring some of the most remarkable predictions of General Relativity: event horizons, singularities and ergoregions. The hitherto invisible landscape of the gravitational Universe is being unveiled before our eyes: the historical direct detection of gravitational waves by the LIGO-Virgo collaboration marks the dawn of a new era of scientific exploration. Gravitational-wave astronomy will allow us to test models of black hole formation, growth and evolution, as well as models of gravitational-wave generation and propagation. It will provide evidence for event horizons and ergoregions, test the theory of General Relativity itself, and may reveal the existence of new fundamental fields. The synthesis of these results has the potential to radically reshape our understanding of the cosmos and of the laws of Nature. The purpose of this work is to present a concise, yet comprehensive overview of the state of the art in the relevant fields of research, summarize important open problems, and lay out a roadmap for future progress.

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.

Making use of definitive new lattice computations of the Standard Model thermodynamics during the quantum chromodynamic (QCD) phase transition, we calculate the enhancement in the mass distribution of primordial black holes (PBHs) due to the softening of the equation of state. We find that the enhancement peaks at approximately $0.7M_\odot$, with the formation rate increasing by at least two orders of magnitude due to the softening of the equation of state at this time, with a range of approximately $0.3M_\odot<M<1.4M_\odot$ at full width half-maximum. PBH formation is increased by a smaller amount for PBHs with masses spanning a large range, $10^{-3}M_\odot<M_{\rm PBH}<10^{3}M_\odot$, which includes the masses of the BHs that LIGO detected. The most significant source of uncertainty in the number of PBHs formed is now due to unknowns in the formation process, rather than from the phase transition. A near scale-invariant density power spectrum tuned to generate a population with mass and merger rate consistent with that detected by LIGO should also produce a much larger energy density of PBHs with solar mass. The existence of BHs below the Chandresekhar mass limit would be a smoking gun for a primordial origin and they could arguably constitute a significant fraction of the cold dark matter density. They also pose a challenge to inflationary model building which seek to produce the LIGO BHs without overproducing lighter PBHs.

Observational constraints on the abundance of primordial black holes (PBHs) constrain the allowed amplitude of the primordial power spectrum on both the smallest and the largest ranges of scales, covering over 20 decades from $1-10^{20}/ \rm{Mpc}$. Despite tight constraints on the allowed fraction of PBHs at their time of formation near horizon entry in the early universe, the corresponding constraints on the primordial power spectrum are quite weak, typically ${\cal P}_\mathcal{R}\lesssim 10^{-2}$ assuming Gaussian perturbations. Motivated by recent claims that the evaporation of just one PBH would destabilise the Higgs vacuum and collapse the universe, we calculate the constraints which follow from assuming there are zero PBHs within the observable universe. This extends the constraints right down to the horizon scale at the end of inflation, but does not significantly tighten the existing power spectrum constraints, even though the constraint on PBH abundance can decrease by up to 46 orders of magnitude. This shows that no future improvement in observational constraints can ever lead to a significant tightening in constraints on inflation (via the power spectrum amplitude). The power spectrum constraints are weak because an order unity perturbation is required in order to overcome pressure forces. We therefore consider an early matter dominated era, during which exponentially more PBHs form for the same initial conditions. We show this leads to far tighter constraints, which approach ${\cal P}_\mathcal{R}\lesssim10^{-9}$, albeit over a smaller range of scales and are very sensitive to when the early matter dominated era ends. Finally, we show that an extended early matter era is incompatible with the argument that an evaporating PBH would destroy the universe, unless the power spectrum amplitude decreases by up to ten orders of magnitude.

We perform three-dimensional simulations of structure formation in the early Universe, when boosting the primordial power spectrum on approximately kpc scales. We demonstrate that our simulations are capable of producing power-law profiles close to the steep $\rho\propto r^{-9/4}$ halo profiles that are commonly assumed to be a good approximation to ultracompact minihalos (UCMHs). However, we show that for more realistic initial conditions in which halos are neither perfectly symmetric nor isolated, the steep power-law profile is disrupted and we find that the Navarro-Frenk-White profile is a better fit to most halos. In the presence of background fluctuations even extreme, nearly spherical initial conditions do not remain exceptional. Nonetheless, boosting the amplitude of initial fluctuations causes all structures to form earlier and thus at larger densities. With sufficiently large amplitude of fluctuations we find that values for the concentration of typical halos in our simulations can become very large. However, despite the signal coming from dark matter annihilation inside the cores of these halos being enhanced, it is still orders-of-magnitude smaller compared to the usually assumed UCMH profile. The upper bound on the primordial power spectrum from the non-observation of UCMHs should therefore be re-evaluated.

Negative absolute temperatures (NAT) are an exotic thermodynamical consequence of quantum physics which has been known since the 1950's (having been achieved in the lab on a number of occasions). Recently, the work of Braun et al (2013) has rekindled interest in negative temperatures and hinted at a possibility of using NAT systems in the lab as dark energy analogues. This paper goes one step further, looking into the cosmological consequences of the existence of a NAT component in the Universe. NAT-dominated expanding Universes experience a borderline phantom expansion ($w<-1$) with no Big Rip, and their contracting counterparts are forced to bounce after the energy density becomes sufficiently large. Both scenarios might be used to solve horizon and flatness problems analogously to standard inflation and bouncing cosmologies. We discuss the difficulties in obtaining and ending a NAT-dominated epoch, and possible ways of obtaining density perturbations with an acceptable spectrum.

The spectral index of scalar perturbations is an important observable that allows us to learn about inflationary physics. In particular, a detection of a significant deviation from a constant spectral index could enable us to rule out the simplest class of inflation models. We investigate whether future observations could rule out canonical single-field slow-roll inflation given the parameters allowed by current observational constraints. We find that future measurements of a constant running (or running of the running) of the spectral index over currently available scales are unlikely to achieve this. However, there remains a large region of parameter space (especially when considering the running of the running) for falsifying the assumed class of slow-roll models if future observations accurately constrain a much wider range of scales.

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.

If the Standard Model (SM) Higgs is weakly coupled to the inflationary sector, the Higgs is expected to be universally in the form of a condensate towards the end of inflation. The Higgs decays rapidly after inflation - via non-perturbative effects - into an out-of-equilibrium distribution of SM species, which thermalize soon afterwards. If the post-inflationary equation of state of the universe is stiff, $w \simeq +1$, the SM species eventually dominate the total energy budget. This provides a natural origin for the relativistic thermal plasma of SM species, required for the onset of the `hot Big Bang' era. The viability of this scenario requires the inflationary Hubble scale $H_*$ to be lower than the instability scale for Higgs vacuum decay, the Higgs not to generate too large curvature perturbations at cosmological scales, and the SM dominance to occur before Big Bang Nucleosynthesis. We show that successful reheating into the SM can only be obtained in the presence of a non-minimal coupling to gravity $\xi \gtrsim 1$, with a reheating temperature of $T_{\rm RH} \gtrsim \mathcal{O}(10^{10})\xi^{3/2}(H_*/10^{14}{\rm GeV})^2~{\rm GeV}$.

We consider the possible decay of the inflaton into curvaton particles during reheating and analyse its effect on curvaton scenarios. Typical decay curvatons are initially relativistic then become non-relativistic, changing the background history of the Universe. We show that this change to the background is the only way in which observational predictions of the scenario are modified. Moreover, once the required amplitude of perturbations is fixed by observation there are no signatures of such decays in other cosmological observables. The decay curvatons can prevent the Universe from becoming dominated by the curvaton condensate, making it impossible to match observations in parts of parameter space. This constrains the branching ratio of the inflaton to curvaton to be less than of order $0.1$ typically. If the branching ratio is below about $10^{-4}$ it has negligible impact on the model parameter space and can be ignored.

We reconsider the observed CMB dipolar asymmetry in the context of open inflation, where a supercurvature mode might survive the bubble nucleation. If such a supercurvature mode modulates the amplitude of the curvature power spectrum, it would easily produce an asymmetry in the power spectrum. We show that current observational data can be accommodated in a three-field model, with simple quadratic potentials and a non-trivial field-space metric. Despite the presence of three fields, we believe this model is so far the simplest that can match current observations. We are able to match the observed strong scale dependence of the dipolar asymmetry, without a fine tuning of initial conditions, breaking slow roll or adding a feature to the evolution of any field.

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

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

Primordial black holes represent a unique probe to constrain the early universe on small scales - providing the only constraints on the primordial power spectrum on the majority of scales. However, these constraints are strongly dependent on even small amounts of non-Gaussianity, which is unconstrained on scales significantly smaller than those visible in the CMB. This paper goes beyond previous considerations to consider the effects of a bispectrum of the equilateral, orthogonal and local shapes with arbitrary magnitude upon the abundance of primordial black holes. Non-Gaussian density maps of the early universe are generated from a given bispectrum and used to place constraints on the small scale power spectrum. When small, we show that the skewness provides an accurate estimate for how the constraint depends on non-Gaussianity, independently of the shape of the bispectrum. We show that the orthogonal template of non-Gaussianity has an order of magnitude weaker effect on the constraints than the local and equilateral templates.

We introduce an alternative parametrisation for the scale dependence of the non-linearity parameter $f_{\rm NL}$ in quasi-local models of non-Gaussianity. Our parametrisation remains valid when $f_{\rm NL}$ changes sign, unlike the commonly adopted power law ansatz $f_{\rm NL}(k) \propto k^{ n_{f_{\rm NL}} }$. We motivate our alternative parametrisation by appealing to the self-interacting curvaton scenario, and as an application, we apply it to the CMB power asymmetry. Explaining the power asymmetry requires a strongly scale dependent non-Gaussianity. We show that regimes of model parameter space where $f_{\rm NL}$ is strongly scale dependent are typically associated with a large $g_{\rm NL}$ and quadrupolar power asymmetry, which can be ruled out by existing observational constraints.

We perform a Bayesian model comparison for scenarios within the quadratic curvaton model, determining the degree to which both are disfavoured with respect to the $\Lambda$CDM concordance model and single-field quadratic inflation, using the recent \emph{Planck} data release. Despite having three additional model parameters, the simplest curvaton scenario is not disfavoured relative to single-field quadratic inflation, and it becomes favoured against this single-field model when we include the joint BICEP/Keck/\emph{Planck} analysis. In all cases we assume an instantaneous inflaton decay and no surviving isocurvature perturbations. Despite the success of \emph{Planck} reaching its forecast measurement accuracy, we show that the current constraints on local non-Gaussianity are insufficiently precise to have any significant impact on the evidence ratios so far. We also determine the precision $\sigma(f_{\mathrm{NL}})$ required by future measurements assuming a fiducial value of $f_{\mathrm{NL}}=-5/4$ or $10.8$ to no longer disfavour the curvaton against the $\Lambda$CDM parametrisation, and we discuss the effect that the predicted increase in precision from future measurements on $f_{\mathrm{NL}}$ may have. We show that our results are not very sensitive to our choice of priors.

Primordial black holes (PBHs) are black holes which may have formed very early on during the radiation dominated era in the early universe. We present here a method by which the large scale perturbations in the density of primordial black holes may be used to place tight constraints on non-gaussianity if PBHs account for dark matter (DM). The presence of local-type non-gaussianity is known to have a significant effect on the abundance of primordial black holes, and modal coupling from the observed CMB scale modes can significantly alter the number density of PBHs that form within different regions of the universe, which appear as DM isocurvature modes. Using the recent \emph{Planck} constraints on isocurvature perturbations, we show that PBHs are excluded as DM candidates for even very small local-type non-gaussianity, $|f_{NL}|\approx0.001$ and remarkably the constraint on $g_{NL}$ is almost as strong. Even small non-gaussianity is excluded if DM is composed of PBHs. If local non-Gaussianity is ever detected on CMB scales, the constraints on the fraction of the universe collapsing into PBHs (which are massive enough to have not yet evaporated) will become much tighter.

The effects of non-gaussianity on the constraints on the primordial curvature perturbation power spectrum from primordial black holes (PBHs) are considered. We extend previous analyses to include the effects of coupling between the modes of the horizon scale at the time the PBH forms and super-horizon modes. We consider terms of up to third order in the Gaussian perturbation. For the weakest constraints on the abundance of PBHs in the early universe (corresponding to a fractional energy density of PBHs of $10^{-5}$ at the time of formation), in the case of gaussian perturbations, constraints on the power spectrum are $\mathcal{P}_{\zeta}<0.05$ but can significantly tighter when even a small amount of non-gaussianity is considered, to $\mathcal{P}_{\zeta}<0.01$, and become approximately $\mathcal{P}_{\zeta}<0.003$ in more special cases. Surprisingly, even when there is negative skew (which naively would suggest fewer areas of high density, leading to weaker constraints), we find that the constraints on the power spectrum become tighter than the purely gaussian case - in strong contrast with previous results. We find that the constraints are highly sensitive to both the non-gaussianity parameters as well as the amplitude of super-horizon perturbations.

We reinspect the calculation for the mass fraction of primordial black holes (PBHs) which are formed from primordial perturbations, finding that performing the calculation using the comoving curvature perturbation $\mathcal{R}_{c}$ in the standard way vastly overestimates the number of PBHs, by many orders of magnitude. This is because PBHs form shortly after horizon entry, meaning modes significantly larger than the PBH are unobservable and should not affect whether a PBH forms or not - this important effect is not taken into account by smoothing the distribution in the standard fashion. We discuss alternative methods and argue that the density contrast, $\Delta$, should be used instead as super-horizon modes are damped by a factor $k^{2}$. We make a comparison between using a Press-Schechter approach and peaks theory, finding that the two are in close agreement in the region of interest. We also investigate the effect of varying the spectral index, and the running of the spectral index, on the abundance of primordial black holes.

We discuss how primordial non-Gaussianity of the curvature perturbation helps to constrain models of the early universe. Observations are consistent with Gaussian initial conditions, compatible with the predictions of the simplest models of inflation. Deviations are constrained to be at the sub percent level, constraining alternative models such as those with multiple fields, non-canonical kinetic terms or breaking the slow-roll conditions. We introduce some of the most important models of inflation which generate non-Gaussian perturbations and provide practical tools on how to calculate the three-point correlation function for a popular class of non-Gaussian models. The current state of the field is summarised and an outlook is given.

We investigate the weak lensing signature of primordial non-Gaussianities of the local type by constraining the magnitude of the weak convergence bi- and trispectra expected for the EUCLID weak lensing survey. Starting from expressions for the weak convergence spectra, bispectra and trispectra, whose relative magnitudes we investigate as a function of scale, we compute their respective signal to noise ratios by relating the polyspectra's amplitude to their Gaussian covariance using a Monte-Carlo technique for carrying out the configuration space integrations. In computing the Fisher-matrix on the non-Gaussianity parameters f_nl, g_nl and tau_nl with a very similar technique, we can derive Bayesian evidences for a violation of the Suyama-Yamaguchi relation tau_nl>=(6 f_nl/5)^2 as a function of the true f_nl and tau_nl-values and show that the relation can be probed down to levels of f_nl~10^2 and tau_nl~10^5. In a related study, we derive analytical expressions for the probability density that the SY-relation is exactly fulfilled, as required by models in which any one field generates the perturbations. We conclude with an outlook on the levels of non-Gaussianity that can be probed with tomographic lensing surveys.

We carry out a comprehensive analysis of the simplest curvaton model, which is based on two non-interacting massive fields. Our analysis encompasses cases where the inflaton and curvaton both contribute to observable perturbations, and where the curvaton itself drives a second period of inflation. We consider both power spectrum and non-Gaussianity observables, and focus on presenting constraints in model parameter space. The fully curvaton-dominated regime is in some tension with observational data, while an admixture of inflaton-generated perturbations improves the fit. The inflating curvaton regime mimics the predictions of Nflation. Some parts of parameter space permitted by power spectrum data are excluded by non-Gaussianity constraints. The recent BICEP2 results [1] require that the inflaton perturbations provide a significant fraction of the total perturbation, ruling out the usual curvaton scenario in which the inflaton perturbations are negligible, though not the admixture regime where both inflaton and curvaton contribute to the spectrum.

Primordial black holes (PBHs) can form in the early Universe from the collapse of rare, large density fluctuations. They have never been observed, but this fact is enough to constrain the amplitude of fluctuations on very small scales which cannot be otherwise probed. Because PBHs form only in very rare large fluctuations, the number of PBHs formed is extremely sensitive to changes in the shape of the tail of the fluctuation distribution - which depends on the amount of non-Gaussianity present. We first study how local non-Gaussianity of arbitrary size up to fifth order affects the abundance and constraints from PBHs, finding that they depend strongly on even small amounts of non-Gaussianity and the upper bound on the allowed amplitude of the power spectrum can vary by several orders of magnitude. The sign of the non-linearity parameters (f_{NL}, g_{NL}, etc) are particularly important. We also study the abundance and constraints from PBHs in the curvaton scenario, in which case the complete non-linear probability distribution is known, and find that truncating to any given order (i.e. to order f_{NL} or g_{NL}, etc) does not give accurate results.

We study the evolution of the non-linear curvature perturbation during perturbative reheating, and hence how observables evolve to their final values which we may compare against observations. Our study includes the evolution of the two trispectrum parameters, $\gnl$ and $\taunl$, as well as the scale dependence of both $\fnl$ and $\taunl$. In general the evolution is significant and must be taken into account, which means that models of multifield inflation cannot be compared to observations without specifying how the subsequent reheating takes place. If the trispectrum is large at the end of inflation, it normally remains large at the end of reheating. In the classes of models we study, it is very hard to generate $\taunl\gg\fnl^2$, regardless of the decay rates of the fields. Similarly, for the classes of models in which $\gnl\simeq\taunl$ during slow--roll inflation, we find the relation typically remains valid during reheating. Therefore it is possible to observationally test such classes of models without specifying the parameters of reheating, even though the individual observables are sensitive to the details of reheating. It is hard to generate an observably large $\gnl$ however. The runnings, $\nfnl$ and $\ntaunl$, tend to satisfy a consistency relation $\ntaunl=(3/2)\nfnl$, but are in general too small to be observed for the class of models considered regardless of reheating timescale.

Primordial Black Holes (PBH's) can form in the early Universe from the collapse of large density fluctuations. Tight observational limits on their abundance constrain the amplitude of the primordial fluctuations on very small scales which can not otherwise be constrained, with PBH's only forming from the extremely rare large fluctuations. The number of PBH's formed is therefore sensitive to small changes in the shape of the tail of the fluctuation distribution, which itself depends on the amount of non-Gaussianity present. We study, for the first time, how quadratic and cubic local non-Gaussianity of arbitrary size (parameterised by f_nl and g_nl respectively) affects the PBH abundance and the resulting constraints on the amplitude of the fluctuations on very small scales. Intriguingly we find that even non-linearity parameters of order unity have a significant impact on the PBH abundance. The sign of the non-Gaussianity is particularly important, with the constraint on the allowed fluctuation amplitude tightening by an order of magnitude as f_nl changes from just -0.5 to 0.5. We find that if PBH's are observed in the future, then regardless of the amplitude of the fluctuations, non-negligible negative f_nl would be ruled out. Finally we show that g_nl can have an even larger effect on the number of PBH's formed than f_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.

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.

Primordial perturbations with wavelengths greater than the observable universe shift the effective background fields in our observable patch from their global averages over the inflating space. This leads to a landscape picture where the properties of our observable patch depend on its location and may significantly differ from the expectation values predicted by the underlying fundamental inflationary model. We show that if multiple fields are present during inflation, this may happen even if our horizon exit would be preceded by only a few e-foldings of inflation. Non-Gaussian statistics are especially affected: for example models of local non-Gaussianity predicting |f_NL|>> 10 over the entire inflating volume can have a probability up to a few tens of percent to generate a non-detectable bispectrum in our observable patch |fNL^{obs.}|<10. In this work we establish systematic connections between the observable local properties of primordial perturbations and the global properties of the inflating space which reflect the underlying high energy physics. We study in detail the implications of both a detection and non-detection of primordial non-Gaussianity by Planck, and discover novel ways of characterising the naturalness of different observational configurations.

We compute the scale dependence of fNL for models of multi-field inflation, allowing for an arbitrary field space metric. We show that, in addition to multi-field effects and self interactions, the curved field space metric provides another source of scale dependence, which arises from the field-space Riemann curvature tensor and its derivatives. The scale dependence may be detectable within the near future if the amplitude of fNL is not too far from the current observational bounds.

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.

We study the effects of perturbative reheating on the evolution of the curvature perturbation \zeta, in two-field inflation models. We use numerical methods to explore the sensitivity of f_NL, n_s and r to the reheating process, and present simple qualitative arguments to explain our results. In general, if a large non-Gaussian signal exists at the start of reheating, it will remain non zero at the end of reheating. Unless all isocurvature modes have completely decayed before the start of reheating, we find that the non-linearity parameter, f_NL, can be sensitive to the reheating timescale, and that this dependence is most appreciable for `runaway' inflationary potentials that only have a minimum in one direction. For potentials with a minimum in both directions, f_NL can also be sensitive to reheating if a mild hierarchy exists between the decay rates of each field. Within the class of models studied, we find that the spectral index n_s, is fairly insensitive to large changes in the field decay rates, indicating that n_s is a more robust inflationary observable, unlike the non-linearity parameter f_NL. Our results imply that the statistics of \zeta, especially f_NL, can only be reliably used to discriminate between models of two-field inflation if the physics of reheating are properly accounted for.

We propose a novel scenario to generate primordial magnetic fields during inflation induced by an oscillating coupling of the electromagnetic field to the inflaton. This resonant mechanism has two key advantages over previous proposals. First of all, it generates a narrow band of magnetic fields at any required wavelength, thereby allaying the usual problem of a strongly blue spectrum and its associated backreaction. Secondly, it avoids the need for a strong coupling as the coupling is oscillating rather than growing or decaying exponentially. Despite these major advantages, we find that the backreaction is still far too large during inflation if the generated magnetic fields are required to have a strength of order 10^{-15} Gauss today on observationally interesting scales. We provide a more general no-go argument, proving that this problem will apply to any model in which the magnetic fields are generated on subhorizon scales and freeze after horizon crossing.

We study the scale dependence of the non-linearity parameters f_NL and g_NL in curvaton models with self-interactions. We show that the spectral indices n_fNL=d ln|f_NL|/(d ln k) and n_gNL=d ln |g_NL|/(d ln k) can take values much greater than the slow--roll parameters and the spectral index of the power spectrum. This means that the scale--dependence of the bi and trispectrum could be easily observable in this scenario with Planck, which would lead to tight additional constraints on the model. Inspite of the highly non-trivial behaviour of f_NL and g_NL in the curvaton models with self-interactions, we find that the model can be falsified if g_NL(k) is also observed.

We study the non-Gaussianity generated during multiple-field inflation. We provide an exact expression for the bispectrum parameter f_NL which is valid beyond the slow-roll regime, valid for certain classes of inflationary models. We then study a new, exact multi-field inflationary model considering a case where the bispectrum grows to observable values at the end of inflation. We show that in this case the trispectrum is also large and may even provide the dominant signal of non-Gaussianity.

We investigate the signature of primordial non-Gaussianities in the weak lensing bispectrum, in particular the signals generated by local, orthogonal and equilateral non-Gaussianities. The questions we address include the signal-to-noise ratio generated in the Euclid weak lensing survey (we find the 1sigma-errors for fNL are 200, 575 and 1628 for local, orthogonal and equilateral non-Gaussianities, respectively), misestimations of fNL if one chooses the wrong non-Gaussianity model (misestimations by up to a factor of +/-3 in fNL are possible, depending on the choice of the model), the probability of noticing such a mistake (improbably large values for the chi^2-functional occur from fNL 200 on), degeneracies of the primordial bispectrum with other cosmological parameters (only the matter density Omega_m plays a significant role), and the subtraction of the much larger, structure-formation generated bispectrum. If a prior on a standard wCDM-parameter set is available from Euclid and Planck, the structure formation bispectrum can be predicted accurately enough for subtraction, and any residual structure formation bispectrum would influence the estimation of fNL to a minor degree. Configuration-space integrations which appear in the evaluation of chi^2-functionals and related quantities can be carried out very efficiently with Monte-Carlo techniques, which reduce the complexity by a factor of O(10^4) while delivering sub-percent accuracies. Weak lensing probes smaller scales than the CMB and hence provide an additional constraint on non-Gaussianities, even though they are not as sensitive to primordial non-Gaussianities as the CMB.

We consider a two-dimensional model of inflation, where the inflationary trajectory is "deformed" by a grazing encounter with an Extra Species/Symmetry Point (ESP) after the observable cosmological scales have left the Hubble radius. The encounter entails a sudden production of particles, whose backreaction causes a bending of the trajectory and a temporary decrease in speed, both of which are sensitive to initial conditions. This "modulated" effect leads to an additional contribution to the curvature perturbation, which can be dominant if the encounter is close. We compute associated non-Gaussianities, the bispectrum and its scale dependence as well as the trispectrum, which are potentially detectable in many cases. In addition, we consider a direct modulation of the coupling to the light field at the ESP via a modulaton field, a mixed scenario whereby the modulaton is identified with a second inflaton, and an extended Extra Species Locus (ESL); all of these scenarios lead to similar additional contributions to observables. We conclude that inflaton interactions throughout inflation are strongly constrained if primordial non-Gaussianities remain unobserved in current experiments such as PLANCK. If they are observed, an ESP encounter leaves additional signatures on smaller scales which may be used to identify the model.

We investigate the scale-dependence of f_NL in the self-interacting curvaton model. We show that the scale-dependence, encoded in the spectral index n_{f_NL}, can be observable by future cosmic microwave background observations, such as CMBpol, in a significant part of the parameter space of the model. We point out that together with information about the trispectrum g_NL, the self-interacting curvaton model parameters could be completely fixed by observations. We also discuss the scale-dependence of g_NL and its implications for the curvaton model, arguing that it could provide a complementary probe in cases where the theoretical value of n_{f_NL} is below observational sensitivity.

In this paper we analyse three models of the early universe, for which the respective mechanisms for generating the curvature perturbation are considered disparate. We find that in fact the mechanisms are very similar, and hence explain why they give rise to a large non-gaussianity. We show that the mechanism for generating the primordial curvature perturbation, and hence the observable non-gaussianity, is similar in both the Curvaton and Modulated Reheating models. In both cases the model can be written in terms of an energy transfer between the constituting fluids. We then show that this is also true for the mechanism of generating the curvature perturbation by symmetry breaking the end of inflation. We then relate this to the non-gaussian contribution to the curvature perturbation and find that it is inversely proportional to the efficiency with which the curvature perturbation is transferred between the fluids. For the first time, we generalise models of modulated reheating to allow for a non-linear energy transfer rate.

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 describe a simple way of incorporating fluctuations of the Hubble scale during the horizon exit of scalar perturbations into the delta N formalism. The dominant effect comes from the dependence of the Hubble scale on low-frequency modes of the inflaton. This modifies the coefficient of the log-enhanced term appearing in the curvature spectrum at second order in field fluctuations. With this modification, the relevant coefficient turns out to be proportional to the second derivative of the tree-level spectrum with respect to the inflaton phi at horizon exit. A logarithm with precisely the same coefficient appears in a calculation of the log-enhancement of the curvature spectrum based purely on the geometry of the reheating surface. We take this agreement as strong support for the proposed implementation of the delta N formalism. Moreover, our analysis makes it apparent that the log-enhancement of the inflationary power-spectrum is indeed physical if this quantity is defined using a global coordinate system on the reheating surface (or any other post-inflationary surface of constant energy density). However, it can be avoided by defining the spectrum using invariant distances on this surface.

We outline the expected constraints on non-Gaussianity from the cosmic microwave background (CMB) with current and future experiments, focusing on both the third (f_{NL}) and fourth-order (g_{NL} and \tau_{NL}) amplitudes of the local configuration or non-Gaussianity. The experimental focus is the skewness (two-to-one) and kurtosis (two-to-two and three-to-one) power spectra from weighted maps. In adition to a measurement of \tau_{NL} and g_{NL} with WMAP 5-year data, our study provides the first forecasts for future constraints on g_{NL}. We describe how these statistics can be corrected for the mask and cut-sky through a window function, bypassing the need to compute linear terms that were introduced for the previous-generation non-Gaussianity statistics, such as the skewness estimator. We discus the ratio A_{NL} = \tau_{NL}/(6f_{NL}/5)^2 as an additional test of single-field inflationary models and discuss the physical significance of each statistic. Using these estimators with WMAP 5-Year V+W-band data out to l_{max}=600 we constrain the cubic order non-Gaussianity parameters \tau_{NL}, and g_{NL} and find -7.4 < g_{NL}/10^5 < 8.2 and -0.6 < \tau_{NL}/10^4 < 3.3 improving the previous COBE-based limit on \tau_{NL} < 10^8 nearly four orders of magnitude with WMAP.

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

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 study the generation of non-Gaussianity in models of hybrid inflation with two inflaton fields, (2-brid inflation). We analyse the region in the parameter and the initial condition space where a large non-Gaussianity may be generated during slow-roll inflation which is generally characterised by a large f_NL, tau_NL and a small g_NL. For certain parameter values we can satisfy tau_NL>>f_NL^2. The bispectrum is of the local type but may have a significant scale dependence. We show that the loop corrections to the power spectrum and bispectrum are suppressed during inflation, if one assume that the fields follow a classical background trajectory. We also include the effect of the waterfall field, which can lead to a significant change in the observables after the waterfall field is destabilised, depending on the couplings between the waterfall and inflaton fields.

We analyse models of inflation in which isocurvature perturbations present during inflation are converted into the primordial curvature perturbation during instant preheating. This can be due to an asymmetry between the fields present either during inflation or during preheating. We consider all the constraints that the model must satisfy in order to be theoretically valid and to satisfy observations. We show that the constraints are very tight in all of the models proposed and special initial conditions are required for the models to work. In the case where the symmetry is strongly broken during inflation the non-Gaussianity parameter f_NL is generally large and negative.

We study the level of primordial non-Gaussianity in slow-roll two-field inflation. Using an analytic formula for the nonlinear parameter f_nl in the case of a sum or product separable potential, we find that it is possible to generate significant non-Gaussianity even during slow-roll inflation with Gaussian perturbations at Hubble exit. In this paper we give the general conditions to obtain large non-Gaussianity and calculate the level of fine-tuning required to obtain this. We present explicit models in which the non-Gaussianity at the end of inflation can exceed the current observational bound of |f_nl|<100.

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 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.

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.

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.