Poincar\'e invariance is a well-tested symmetry of nature and sits at the core of our description of relativistic particles and gravity. At the same time, in most systems Poincar\'e invariance is not a symmetry of the ground state and is hence broken spontaneously. This phenomenon is ubiquitous in cosmology where Lorentz boosts are spontaneously broken by the existence of a preferred reference frame in which the universe is homogeneous and isotropic. This motivates us to study scattering amplitudes without requiring invariance of the interactions under Lorentz boosts. In particular, using on-shell methods and assuming massless, relativistic and luminal particles of any spin, we show that the allowed interactions around Minkowski spacetime are severely constrained by unitarity and locality in the form of consistent factorization. The existence of an interacting massless spin-2 particle enforces (analytically continued) three-particle amplitudes to be Lorentz invariant, even those that do not involve a graviton, such as cubic scalar couplings. We conjecture this to be true for all n-particle amplitudes. Also, particles of spin S > 2 cannot self-interact nor can be minimally coupled to gravity, while particles of spin S > 1 cannot have electric charge. Given the growing evidence that free gravitons are well described by massless, luminal relativistic particles, our results imply that cubic graviton interactions in Minkowski must be those of general relativity up to a unique Lorentz-invariant higher derivative correction of mass dimension 9. Finally, we discuss a surprising and general lesson revealed by applying our analysis to the sub-Hubble limit of inflation and Dark Energy models.

The space of inflationary models is vast, containing wide varieties of mechanisms, symmetries, and spectra of particles. Consequently, the space of observational signatures is similarly complex. Hence, it is natural to look for boundaries of the space of models and their signatures. In this paper, we explore the possible symmetries associated with the primordial cosmological perturbations and their correlators in the asymptotic future. Assuming the observed homogeneity, isotropy and (approximate) scale invariance, we prove three main results. First, correlation functions of scalar metric fluctuations are uniquely characterized by soft theorems and are free from ambiguity under field redefinitions. Second, whatever the particle content and interactions, when the standard soft theorems apply, invariance under de Sitter boosts (linearly realized conformal invariance) is only possible if all connected correlators vanish identically, i.e. if the theory is free. Third, conformal invariance is the largest set of linearly realized (bosonic) symmetries of the correlators of any single scalar, irrespectively of any soft theorems or particle content.

The effect of spatial curvature on primordial perturbations is controlled by $ \Omega_{K,0}/c_{s}^{2} $, where $ \Omega_{K,0} $ is today's fractional density of spatial curvature and $ c_{s} $ is the speed of sound during inflation. Here we study these effects in the limit $ c_{s}\ll 1 $. First, we show that the standard cosmological soft theorems in flat universes are violated in curved universes and the soft limits of correlators can have non-universal contributions even in single-clock inflation. This is a consequence of the fact that, in the presence of spatial curvature, there is a gap between the spectrum of residual diffeomorphisms and that of physical modes. Second, there are curvature corrections to primordial correlators, which are not scale invariant. We provide explicit formulae for these corrections to the power spectrum and the bispectrum to linear order in curvature in single-clock inflation. We show that the large-scale CMB anisotropies could provide interesting new constraints on these curvature effects, and therefore on $ \Omega_{K,0}/c_{s}^{2} $, but it is necessary to go beyond our linear-order treatment.

We provide a complete classification of Poincar\'e-invariant scalar field theories with an enlarged set of classical symmetries to leading order in derivatives, namely for the so-called $P(X,\phi)$ theories, in two or more spacetime dimensions. We find only three possibilities: Dirac-Born-Infeld, Cuscuton and Scaling theories. The latter two classes of actions involve an arbitrary function of the scalar field. As an application, we use the scaling symmetry to derive an infinite set of constraints on the Wilsonian coefficients of the low-energy Effective Field Theory. Furthermore, we study the extension of these results to cosmological (FLRW) and (Anti-)de Sitter spacetimes. We find in particular that the Cuscuton action has a generic set of symmetries around any background spacetime that possesses Killing vector fields, while the DBI actions have well-known analogues that we summarize explicitly.

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.

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.

We derive cosmological soft theorems for solids coupled to gravity. To this end, we first derive all cosmological adiabatic modes for solids, which display the interesting novelty of non-vanishing anisotropic stresses on large scales. Then, from the corresponding symmetries of the action of perturbations we compute the leading order related soft theorems using the operator product expansion. For the scalar bispectrum, we re-derive the result that Maldacena's consistency relation is recovered only upon angular averaging over the long mode direction. In addition, we find theorems for soft tensor and vector perturbations. In passing, we also clarify the derivation of these soft theorems in gauges where no residual diffeomorphisms exist.

Adiabatic modes are cosmological perturbations that are locally indistinguishable from a (large) change of coordinates. At the classical level, they provide model independent solutions. At the quantum level, they lead to soft theorems for cosmological correlators. We present a systematic derivation of adiabatic modes in spatially-flat cosmological backgrounds with asymptotically-perfect fluids. We find several new adiabatic modes including vector, time-dependent tensor and time-dependent scalar modes. The new vector and tensor modes decay with time in standard cosmologies but are the leading modes in contracting universes. We present a preliminary derivation of the related soft theorems. In passing, we discuss a distinction between classical and quantum adiabatic modes, we clarify the subtle nature of Weinberg second adiabatic mode and point out that the adiabatic nature of a perturbation is a gauge dependent statement.

Distortions of the Cosmic Microwave Background energy spectrum of the $\mu$ type are sensitive to the primordial power spectrum through the dissipation of curvature perturbations on scales $k\simeq 50$ - $10^4\,\mathrm{Mpc}^{-1}$. Their angular correlation with large-scale temperature anisotropies is then sensitive to the squeezed limit of the primordial bispectrum. For inflationary models obeying the single-field consistency relation, we show that the observed $\mu T$ angular correlation that would correspond to the local shape vanishes exactly. All leading non-primordial contributions, including all non-linear production and projection effects, are of the "equilateral shape", namely suppressed by $k^2/{\cal H}_f^2$, where ${\cal H}_f\simeq 10^{-1}\,{\rm Mpc}^{-1}$ is the Hubble radius at the end of the $\mu$-era. Therefore, these non-primordial contributions are orthogonal to a potential local primordial signal (e.g. from multi-field inflation). Moreover, they are very small in amplitude. Our results strengthen the position of $\mu$ distortions as the ultimate probe of local primordial non-Gaussianity.

Long-short mode coupling during inflation, encoded in the squeezed bispectrum of curvature perturbations, induces a dependence of the local, small-scale power spectrum on long-wavelength perturbations, leading to a scale-dependent halo bias. While this scale dependence is absent in the large-scale limit for single-field inflation models that satisfy the consistency relation, certain models such as resonant non-Gaussianity show a peculiar behavior on intermediate scales. We reconsider the predictions for the halo bias in this model by working in Conformal Fermi Coordinates, which isolate the physical effects of long-wavelength perturbations on short-scale physics. We find that the bias oscillates with scale with an envelope similar to that of equilateral non-Gaussianity. Moreover, the bias shows a peculiar modulation with the halo mass. Unfortunately, we find that upcoming surveys will be unable to detect the signal because of its very small amplitude. We also discuss non-Gaussianity due to interactions between the inflaton and massive fields: our results for the bias agree with those in the literature.

We make progress towards an analytical understanding of the regime of validity of perturbation theory for large scale structures and the nature of some non-perturbative corrections. We restrict ourselves to 1D gravitational collapse, for which exact solutions before shell crossing are known. We review the convergence of perturbation theory for the power spectrum, recently proven by McQuinn and White, and extend it to non-Gaussian initial conditions and the bispectrum. In contrast, we prove that perturbation theory diverges for the real space two-point correlation function and for the probability density function (PDF) of the density averaged in cells and all the cumulants derived from it. We attribute these divergences to the statistical averaging intrinsic to cosmological observables, which, even on very large and "perturbative" scales, gives non-vanishing weight to all extreme fluctuations. Finally, we discuss some general properties of non-perturbative effects in real space and Fourier space.

We derive soft theorems for single-clock cosmologies that enjoy a shift symmetry. These so-called consistency conditions arise from a combination of a large diffeomorphism and the internal shift-symmetry and fix the squeezed limit of all correlators with a soft scalar mode. As an application, we show that our results reproduce the squeezed bispectrum for Ultra-slow-roll inflation, a particular shift-symmetric, non-attractor model which is known to violate Maldacena's consistency relation. Similar results have been previously obtained by Mooij and Palma using background-wave methods. Our results shed new light on the infrared structure of single-clock cosmological spacetimes.

A shift symmetry is a ubiquitous ingredient in inflationary models, both in effective constructions and in UV-finite embeddings such as string theory. It has also been proposed to play a key role in certain Dark Energy and Dark Matter models. Despite the crucial role it plays in cosmology, the observable, model independent consequences of a shift symmetry are yet unknown. Here, assuming an exact shift symmetry, we derive these consequences for single-clock cosmologies within the framework of the Effective Field Theory of Inflation. We find an infinite set of relations among the otherwise arbitrary effective coefficients, which relate non-Gaussianity to their time dependence. For example, to leading order in derivatives, these relations reduce the infinitely many free functions in the theory to just a single one. Our Effective Theory of shift-symmetric cosmologies describes, among other systems, perfect and imperfect superfluids coupled to gravity and driven superfluids in the decoupling limit. Our results are the first step to determine observationally whether a shift symmetry is at play in the laws of nature and whether it is broken by quantum gravity effects.

Non-Gaussianities of dynamical origin are disentangled from primordial ones using the formalism of large deviation statistics with spherical collapse dynamics. This is achieved by relying on accurate analytical predictions for the one-point probability distribution function (PDF) and the two-point clustering of spherically-averaged cosmic densities (sphere bias). Sphere bias extends the idea of halo bias to intermediate density environments and voids as underdense regions. In the presence of primordial non-Gaussianity, sphere bias displays a strong scale dependence relevant for both high and low density regions, which is predicted analytically. The statistics of densities in spheres are built to model primordial non-Gaussianity via an initial skewness with a scale-dependence that depends on the bispectrum of the underlying model. The analytical formulas with the measured nonlinear dark matter variance as input are successfully tested against numerical simulations. For local non-Gaussianity with a range from $f_{\rm NL}=-100$ to $+100$ they are found to agree within 2\% or better for densities $\rho\in[0.5,3]$ in spheres of radius 15 Mpc$/h$ down to $z=0.35$. The validity of the large deviation statistics formalism is thereby established for all observationally relevant local-type departures from perfectly Gaussian initial conditions. The corresponding estimators for the amplitude of the nonlinear variance $\sigma_8$ and primordial skewness $f_{\rm NL}$ are validated using a fiducial joint maximum likelihood experiment. The influence of observational effects and the prospects for a future detection of primordial non-Gaussianity from joint one- and two-point densities-in-spheres statistics are discussed.

Future observations of cosmic microwave background (CMB) polarisation have the potential to answer some of the most fundamental questions of modern physics and cosmology. In this paper, we list the requirements for a future CMB polarisation survey addressing these scientific objectives, and discuss the design drivers of the CORE space mission proposed to ESA in answer to the "M5" call for a medium-sized mission. The rationale and options, and the methodologies used to assess the mission's performance, are of interest to other future CMB mission design studies. CORE is designed as a near-ultimate CMB polarisation mission which, for optimal complementarity with ground-based observations, will perform the observations that are known to be essential to CMB polarisation scienceand cannot be obtained by any other means than a dedicated space mission.

We forecast the scientific capabilities to improve our understanding of cosmic inflation of CORE, a proposed CMB space satellite submitted in response to the ESA fifth call for a medium-size mission opportunity. The CORE satellite will map the CMB anisotropies in temperature and polarization in 19 frequency channels spanning the range 60-600 GHz. CORE will have an aggregate noise sensitivity of $1.7 \mu$K$\cdot \,$arcmin and an angular resolution of 5' at 200 GHz. We explore the impact of telescope size and noise sensitivity on the inflation science return by making forecasts for several instrumental configurations. This study assumes that the lower and higher frequency channels suffice to remove foreground contaminations and complements other related studies of component separation and systematic effects, which will be reported in other papers of the series "Exploring Cosmic Origins with CORE." We forecast the capability to determine key inflationary parameters, to lower the detection limit for the tensor-to-scalar ratio down to the $10^{-3}$ level, to chart the landscape of single field slow-roll inflationary models, to constrain the epoch of reheating, thus connecting inflation to the standard radiation-matter dominated Big Bang era, to reconstruct the primordial power spectrum, to constrain the contribution from isocurvature perturbations to the $10^{-3}$ level, to improve constraints on the cosmic string tension to a level below the presumptive GUT scale, and to improve the current measurements of primordial non-Gaussianities down to the $f_{NL}^{\rm local} < 1$ level. For all the models explored, CORE alone will improve significantly on the present constraints on the physics of inflation. Its capabilities will be further enhanced by combining with complementary future cosmological observations.

Gravity is a non-linear theory, and hence, barring cancellations, the initial super-horizon perturbations produced by inflation must contain some minimum amount of mode coupling, or primordial non-Gaussianity. In single-field slow-roll models, where this lower bound is saturated, non-Gaussianity is controlled by two observables: the tensor-to-scalar ratio, which is uncertain by more than fifty orders of magnitude; and the scalar spectral index, or tilt, which is relatively well measured. It is well known that to leading and next-to-leading order in derivatives, the contributions proportional to the tilt disappear from any local observable, and suspicion has been raised that this might happen to all orders, allowing for an arbitrarily low amount of primordial non-Gaussianity. Employing Conformal Fermi Coordinates, we show explicitly that this is not the case. Instead, a contribution of order the tilt appears in local observables. In summary, the floor of physical primordial non-Gaussianity in our Universe has a squeezed-limit scaling of $k_\ell^2/k_s^2$, similar to equilateral and orthogonal shapes, and a dimensionless amplitude of order $0.1\times(n_\mathrm{s}-1)$.

We argue that the non-detection of primordial tensor modes has taught us a great deal about the primordial universe. In single-field slow-roll inflation, the current upper bound on the tensor-to-scalar ratio, $r < 0.07$ $(95 \% ~CL)$, implies that the Hubble slow-roll parameters obey $\varepsilon \ll \eta$, and therefore establishes the existence of a new hierarchy. We dub this regime the conformal limit of (slow-roll) inflation, and show that it includes Starobinsky-like inflation as well as all viable single-field models with a sub-Planckian field excursion. In this limit, all primordial correlators are constrained by the full conformal group to leading non-trivial order in slow-roll. This fixes the power spectrum and the full bispectrum, and leads to the "conformal" shape of non-Gaussianity. The size of non-Gaussianity is related to the running of the spectral index by a consistency condition, and therefore it is expected to be small. In passing, we clarify the role of boundary terms in the $\zeta$ action, the order to which constraint equations need to be solved, and re-derive our results using the Wheeler-deWitt formalism.

We use the recent observations of Cosmic Microwave Background temperature and polarization anisotropies provided by the Planck satellite experiment to place constraints on the running $\alpha_\mathrm{s} = \mathrm{d}n_{\mathrm{s}} / \mathrm{d}\log k$ and the running of the running $\beta_{\mathrm{s}} = \mathrm{d}\alpha_{\mathrm{s}} / \mathrm{d}\log k$ of the spectral index $n_{\mathrm{s}}$ of primordial scalar fluctuations. We find $\alpha_\mathrm{s}=0.011\pm0.010$ and $\beta_\mathrm{s}=0.027\pm0.013$ at $68\%\,\mathrm{CL}$, suggesting the presence of a running of the running at the level of two standard deviations. We find no significant correlation between $\beta_{\mathrm{s}}$ and foregrounds parameters, with the exception of the point sources amplitude at $143\,\mathrm{GHz}$, $A^{PS}_{143}$, which shifts by half sigma when the running of the running is considered. We further study the cosmological implications of such preference for $\alpha_\mathrm{s},\beta_\mathrm{s}\sim0.01$ by including in the analysis the lensing amplitude $A_L$, the curvature parameter $\Omega_k$, and the sum of neutrino masses $\sum m_{\nu}$. We find that when the running of the running is considered, Planck data are more compatible with the standard expectations of $A_L = 1$ and $\Omega_k = 0$ but still hint at possible deviations. The indication for $\beta_\mathrm{s} > 0$ survives at two standard deviations when external datasets such as BAO and CFHTLenS are included in the analysis, and persists at $\sim 1.7$ standard deviations when CMB lensing is considered. We discuss the possibility of constraining $\beta_\mathrm{s}$ with current and future measurements of CMB spectral distortions, showing that an experiment like PIXIE could provide strong constraints on $\alpha_\mathrm{s}$ and $\beta_\mathrm{s}$.

Primordial non-Gaussianity (PNG) in Large Scale Structures is obfuscated by the many additional sources of non-linearity. Within the Effective Field Theory approach to Standard Perturbation Theory, we show that matter non-linearities in the bispectrum can be modeled sufficiently well to strengthen current bounds with near future surveys, such as Euclid. We find that the EFT corrections are crucial to this improvement in sensitivity. Yet, our understanding of non-linearities is still insufficient to reach important theoretical benchmarks for equilateral PNG, while, for local PNG, our forecast is more optimistic. We consistently account for the theoretical error intrinsic to the perturbative approach and discuss the details of its implementation in Fisher forecasts.

We discuss the implications of a PIXIE-like experiment, which would measure $\mu$-type spectral distortions of the CMB at a level of $\sigma_{\mu}=(1/n)\times 10^{-8}$, with $n\geq1$ representing an improved sensitivity (e.g. $n=10$ corresponds to PRISM). Using Planck data and considering the six-parameter $\Lambda$CDM model, we compute the posterior for $\mu_8\equiv\mu\times 10^{8}$ and find $\mu_8=1.57^{+0.11}_{-0.13}$ ($68\%\,\mathrm{CL}$). This becomes $\mu_{8} = 1.28^{+0.30}_{-0.52}$ ($68\%\,\mathrm{CL}$) when the running $\alpha_\mathrm{s}$ of the spectral index is included. We point out that a sensitivity of about $3\times$ PIXIE implies a guaranteed discovery: $\mu$-distortion is detected or $\alpha_\mathrm{s}\geq 0$ is excluded (both at $95\%\,\mathrm{CL}$ or higher). This threshold sensitivity sets a clear benchmark for CMB spectrometry. For a combined analysis of PIXIE and current Planck data, we discuss the improvement on measurements of the tilt $n_\mathrm{s}$ and the running $\alpha_\mathrm{s}$ and the dependence on the choice of the pivot. A fiducial running of $\alpha_\mathrm{s}=-0.01$ (close to the Planck best-fit) leads to a detection of negative running at $2\sigma$ for $5\times$ PIXIE. A fiducial running of $\alpha_\mathrm{s}=-0.02$, still compatible with Planck, requires $3\times$ PIXIE to rule out $\alpha_\mathrm{s} = 0$ (at $95\%\,\mathrm{CL}$). We propose a convenient and compact visualization of the improving constraints on the tilt, running and tensor-to-scalar ratio.

A perturbative description of Large Scale Structure is a cornerstone of our understanding of the observed distribution of matter in the universe. Renormalization is an essential and defining step to make this description physical and predictive. Here we introduce a systematic renormalization procedure, which neatly associates counterterms to the UV-sensitive diagrams order by order, as it is commonly done in quantum field theory. As a concrete example, we renormalize the one-loop power spectrum and bispectrum of both density and velocity. In addition, we present a series of results that are valid to all orders in perturbation theory. First, we show that while systematic renormalization requires temporally non-local counterterms, in practice one can use an equivalent basis made of local operators. We give an explicit prescription to generate all counterterms allowed by the symmetries. Second, we present a formal proof of the well-known general argument that the contribution of short distance perturbations to large scale density contrast $\delta$ and momentum density $\mathbf\pi(\mathbf k)$ scale as $k^2$ and $k$, respectively. Third, we demonstrate that the common practice of introducing counterterms only in the Euler equation when one is interested in correlators of $ \delta$ is indeed valid to all orders.

We develop the effective theory of large-scale structure for non-Gaussian initial conditions. The effective stress tensor in the dark matter equations of motion contains new operators, which originate from the squeezed limit of the primordial bispectrum. Parameterizing the squeezed limit by a scaling and an angular dependence, captures large classes of primordial non-Gaussianity. Within this parameterization, we classify the possible contributions to the effective theory. We show explicitly how all terms consistent with the symmetries arise from coarse graining the dark matter equations of motion and its initial conditions. We also demonstrate that the system is closed under renormalization and that the basis of correction terms is therefore complete. The relevant corrections to the matter power spectrum and bispectrum are computed numerically and their relative importance is discussed.

Cosmological observations have provided us with the measurement of just three numbers that characterize the very early universe: $ 1-n_{s} $, $ N $ and $\ln\Delta_R^2$. Although each of the three numbers individually carries limited information about the physics of inflation, one may hope to extract non-trivial information from relations among them. Invoking minimality, namely the absence of ad hoc large numbers, we find two viable and mutually exclusive inflationary scenarios. The first is the well-known inverse relation between $1- n_{s} $ and $ N $. The second implies a new relation between $ 1-n_{s} $ and $\ln\Delta_R^2$, which might provide us with a handle on the beginning of inflation and predicts the intriguing $\textit{lower}$ bound on the tensor-to-scalar ratio $ r> 0.006 $ ($ 95\% $ CL).

(abridged version) The separate universe conjecture states that in General Relativity a density perturbation behaves locally (i.e. on scales much smaller than the wavelength of the mode) as a separate universe with different background density and curvature. We prove this conjecture for a spherical compensated tophat density perturbation of arbitrary amplitude and radius in $\Lambda$CDM. We then use Conformal Fermi Coordinates to generalize this result to scalar perturbations of arbitrary configuration and scale. In this case, the separate universe conjecture holds for the isotropic part of the perturbations. The anisotropic part on the other hand is exactly captured by a tidal field in the Newtonian form. We show that the separate universe picture is restricted to scales larger than the sound horizons of all fluid components. We then derive an expression for the locally measured matter bispectrum induced by a long-wavelength mode of arbitrary wavelength. We show that nonlinear gravitational dynamics does not generate observable contributions that scale like local-type non-Gaussianity $f_{\rm NL}^{\rm loc}$, and hence does not contribute to a scale-dependent galaxy bias $\Delta b \propto k^{-2}$ on large scales; rather, the locally measurable long-short mode coupling assumes a form essentially identical to subhorizon perturbation theory results, once the long-mode density perturbation is replaced by the synchronous-comoving gauge density perturbation. Apparent $f_{\rm NL}^{\rm loc}$-type contributions arise through projection effects on photon propagation, which depend on the specific large-scale structure tracer and observable considered, and are in principle distinguishable from the local mode coupling induced by gravity. We conclude that any observation of $f_{\rm NL}^{\rm loc}$ beyond these projection effects signals a departure from standard single-clock inflation.

Fermi Normal Coordinates (FNC) are a useful frame for isolating the locally observable, physical effects of a long-wavelength spacetime perturbation. Their cosmological application, however, is hampered by the fact that they are only valid on scales much smaller than the horizon. We introduce a generalization that we call Conformal Fermi Coordinates (CFC). CFC preserve all the advantages of FNC, but in addition are valid outside the horizon. They allow us to calculate the coupling of long- and short-wavelength modes on all scales larger than the sound horizon of the cosmological fluid, starting from the epoch of inflation until today, by removing the complications of the second order Einstein equations to a large extent, and eliminating all gauge ambiguities. As an application, we present a calculation of the effect of long-wavelength tensor modes on small scale density fluctuations. We recover previous results, but clarify the physical content of the individual contributions in terms of locally measurable effects and "projection" terms.

We study the bispectrum in the Effective Field Theory of Large Scale Structure, consistently accounting for the effects of short-scale dynamics. We begin by proving that, as long as the theory is perturbative, it can be formulated to arbitrary order using only operators that are local in time. We then derive all the new operators required to cancel the UV-divergences and obtain a physically meaningful prediction for the one-loop bispectrum. In addition to new, subleading stochastic noises and the viscosity term needed for the one-loop power spectrum, we find three new effective operators. The three new parameters can be constrained by comparing with N-body simulations. The best fit is precisely what is suggested by the structure of UV-divergences, hence justifying a formula for the EFTofLSS bispectrum whose only fitting parameter is already fixed by the power spectrum. This result predicts the bispectrum of N-body simulations up to $k \approx 0.22\, h\, \text{Mpc}^{-1}$ at $z=0$, an improvement by nearly a factor of two as compared to one-loop standard perturbation theory.

We construct a class of random potentials for N >> 1 scalar fields using non-equilibrium random matrix theory, and then characterize multifield inflation in this setting. By stipulating that the Hessian matrices in adjacent coordinate patches are related by Dyson Brownian motion, we define the potential in the vicinity of a trajectory. This method remains computationally efficient at large N, permitting us to study much larger systems than has been possible with other constructions. We illustrate the utility of our approach with a numerical study of inflation in systems with up to 100 coupled scalar fields. A significant finding is that eigenvalue repulsion sharply reduces the duration of inflation near a critical point of the potential: even if the curvature of the potential is fine-tuned to be small at the critical point, small cross-couplings in the Hessian cause the curvature to grow in the neighborhood of the critical point.

We compute the renormalized two-point functions of density, divergence and vorticity of the velocity in the Effective Field Theory of Large Scale Structures. Because of momentum and mass conservation, the corrections from short scales to the large-scale power spectra of density, divergence and vorticity must start at order $k^{4}$. For the vorticity this constitutes one of the two leading terms. Exact (approximated) self-similarity of an Einstein-de Sitter ($\Lambda$CDM) background fixes the time dependence so that the vorticity power spectrum at leading order is determined by the symmetries of the problem and the power spectrum around the non-linear scale. We show that to cancel all divergences in the velocity correlators one needs new counterterms. These fix the definition of velocity and do not represent new properties of the system. For an Einstein-de Sitter universe, we show that all three renormalized cross- and auto-correlation functions have the same structure but different numerical coefficients, which we compute. We elucidate the differences between using momentum and velocity.

The leading locally observable effect of a long-wavelength metric perturbation corresponds to a tidal field. We derive the tidal field induced by scalar, vector, and tensor perturbations, and use second order perturbation theory to calculate the effect on the locally measured small-scale density fluctuations. For sub-horizon scalar perturbations, we recover the standard perturbation theory result ($F_2$ kernel). For tensor modes of wavenumber $k_L$, we find that effects persist for $k_L\tau \gg 1$, i.e. even long after the gravitational wave has entered the horizon and redshifted away, i.e. it is a "fossil" effect. We then use these results, combined with the "ruler perturbations" of arXiv:1204.3625, to predict the observed distortion of the small-scale matter correlation function induced by a long-wavelength tensor mode. We also estimate the observed signal in the B mode of the cosmic shear from a gravitational wave background, including both tidal (intrinsic alignment) and projection (lensing) effects. The non-vanishing tidal effect in the $k_L\tau \gg 1$ limit significantly increases the intrinsic alignment contribution to shear B modes, especially at low redshifts $z \lesssim 2$.

The squeezed limit of the three-point function of cosmological perturbations is a powerful discriminant of different models of the early Universe. We present a conceptually simple and complete framework to relate any primordial bispectrum in this limit to late time observables, such as the CMB temperature bispectrum and the scale-dependent halo bias. We employ a series of convenient coordinate transformations to capture the leading non-linear effects of cosmological perturbation theory on these observables. This makes crucial use of Fermi Normal Coordinates and their conformal generalization, which we introduce here and discuss in detail. As an example, we apply our formalism to standard slow-roll single-field inflation. We show explicitly that Maldacena's results for the squeezed limits of the scalar bispectrum [proportional to (ns-1) in comoving gauge] and the tensor-scalar-scalar bispectrum lead to no deviations from a Gaussian universe, except for projection effects. In particular, the primordial contributions to the squeezed CMB bispectrum and scale dependent halo bias vanish, and there are no primordial "fossil" correlations between long-wavelength tensor perturbations and small-scale perturbations. The contributions to observed correlations are then only due to projection effects such as gravitational lensing and redshift perturbations.

We carry out a comprehensive analysis of the possible constraints on cosmological and astrophysical parameters achievable with measurements of the thermal Sunyaev-Zel'dovich (tSZ) power spectrum from upcoming full-sky CMB observations, with a particular focus on one-parameter extensions to the LCDM standard model involving local primordial non-Gaussianity (described by fNL) and massive neutrinos (described by Mnu). We include all of the relevant physical effects due to these additional parameters, including the change to the halo mass function and the scale-dependent halo bias induced by local primordial non-Gaussianity. [...] We compute forecasts for Planck, PIXIE, and a cosmic variance (CV)-limited experiment, using multifrequency subtraction to remove foregrounds and implementing two masking criteria based on the ROSAT and eROSITA cluster catalogs to reduce the significant CV errors at low multipoles. We find that Planck can detect the tSZ power spectrum with >30-sigma significance, regardless of the masking scenario. However, neither Planck or PIXIE is likely to provide competitive constraints on fNL from the tSZ power spectrum due to CV noise at low-ell overwhelming the unique signature of the scale-dependent bias. A future CV-limited experiment could provide a 3-sigma detection of fNL~37, which is the WMAP9 maximum-likelihood result. The outlook for neutrino masses is more optimistic: Planck can reach levels comparable to the current upper bounds ~<0.3 eV with conservative assumptions about the intracluster medium (ICM); stronger ICM priors could allow Planck to provide 1-2-sigma evidence for massive neutrinos from the tSZ power spectrum, depending on the true value of the sum of the neutrino masses. We also forecast a ~<10% constraint on the outer slope of the ICM pressure profile using the unmasked Planck tSZ power spectrum.

Because the inflationary mechanism is extremely sensitive to UV-physics, the construction of theoretically robust models of inflation provides a unique window on Planck-scale physics. We review efforts to use an axion with a shift symmetry to ensure a prolonged slow-roll background evolution. The symmetry dictates which operators are allowed, and these in turn determine the observational predictions of this class of models, which include observable gravitational waves (potentially chiral), oscillations in all primordial correlators, specific deviations from scale invariance and Gaussianity and primordial black holes. We discuss the constraints on this class of models in light of the recent Planck results and comment on future perspectives. The shift symmetry is very useful in models of large-field inflation, which typically have monomial potentials, but it cannot explain why two or more terms in the potential are fine-tuned against each other, as needed for typical models of small-field inflation. Therefore some additional symmetries or fine-tuning will be needed if forthcoming experiments will constrain the tensor-to-scalar ratio to be r < 0.01.

When inflation is driven by a pseudo-scalar field \chi coupled to vectors as \alpha/4 \chi F \tilde F, this coupling may lead to a copious production of gauge quanta, which in turns induces non-Gaussian and non-scale invariant corrections to curvature perturbations. We point out that this mechanism is generically at work in a broad class of inflationary models in supergravity hence providing them with a rich set of observational predictions. When the gauge fields are massless, significant effects on CMB scales emerge only for relatively large \alpha. We show that in this regime, the curvature perturbations produced at the last stages of inflation have a relatively large amplitude that is of the order of the upper bound set by the possible production of primordial black holes by non-Gaussian perturbations. On the other hand, within the supergravity framework described in our paper, the gauge fields can often acquire a mass through a coupling to additional light scalar fields. Perturbations of these fields modulate the duration of inflation, which serves as a source for non-Gaussian perturbations of the metric. In this regime, the bounds from primordial black holes are parametrically satisfied and non-Gaussianity of the local type can be generated at the observationally interesting level f_NL =O(10).

Standard perturbation theory (SPT) for large-scale matter inhomogeneities is unsatisfactory for at least three reasons: there is no clear expansion parameter since the density contrast is not small on all scales; it does not fully account for deviations at large scales from a perfect pressureless fluid induced by short-scale non-linearities; for generic initial conditions, loop corrections are UV-divergent, making predictions cutoff dependent and hence unphysical. The Effective Field Theory of Large Scale Structures successfully addresses all three issues. Here we focus on the third one and show explicitly that the terms induced by integrating out short scales, neglected in SPT, have exactly the right scale dependence to cancel all UV-divergences at one loop, and this should hold at all loops. A particularly clear example is an Einstein deSitter universe with no-scale initial conditions P_in=A k^n. After renormalizing the theory, we use self-similarity to derive a very simple result for the final power spectrum for any n, excluding two-loop corrections and higher. We show how the relative importance of different corrections depend on n. For n=-1.5, relevant for our universe, pressure and dissipative corrections are more important than the two-loop corrections.

We know very little about primordial curvature perturbations on scales smaller than about a Mpc. Measurements of the mu-type distortion of the CMB spectrum provide the unique opportunity to probe these scales over the unexplored range from 50 to 10^4 Mpc^-1. This is a very clean probe, in that it relies only on well-understood linear evolution. We point out that correlations between mu-distortion and temperature anisotropies can be used to test Gaussianity at these very small scales. In particular the mu-T cross correlation is proportional to the very squeezed limit of the primordial bispectrum and hence measures fNL^loc, while mu-mu is proportional to the primordial trispectrum and measures tauNL. We present a Fisher matrix forecast of the observational constraints.

Spectral distortion of the cosmic microwave background provides a unique opportunity to probe primordial perturbations on very small scales by performing large-scale measurements. We discuss in a systematic and pedagogic way all the relevant physical phenomena involved in the production and evolution of the mu-type spectral distortion. Our main results agree with previous estimates (in particular we show that a recently found factor of 3/4 arises from relativistic corrections to the wave energy). We also discuss several subleading corrections such as adiabatic cooling and the effects of bulk viscosity, baryon loading and photon heat conduction. Finally we calculate the transfer function for mu-distortions between the end of the mu-era and now.

Models of axion inflation are particularly interesting since they provide a natural justification for the flatness of the potential over a super-Planckian distance, namely the approximate shift-symmetry of the inflaton. In addition, most of the observational consequences are directly related to this symmetry and hence are correlated. Large tensor modes can be accompanied by the observable effects of a the shift-symmetric coupling $\phi F\tilde F$ to a gauge field. During inflation this coupling leads to a copious production of gauge quanta and consequently a very distinct modification of the primordial curvature perturbations. In this work we compare these predictions with observations. We find that the leading constraint on the model comes from the CMB power spectrum when considering both WMAP 7-year and ACT data. The bispectrum generated by the non-Gaussian inverse-decay of the gauge field leads to a comparable but slightly weaker constraint. There is also a constraint from mu-distortion using TRIS plus COBE/FIRAS data, but it is much weaker. Finally we comment on a generalization of the model to massive gauge fields. When the mass is generated by some light Higgs field, observably large local non-Gaussianity can be produced.

Models of inflation based on axions, which owe their popularity to the robustness against UV corrections, have also a very distinct class of signatures. The relevant interactions of the axion are a non-perturbative oscillating contribution to the potential and a shift-symmetric coupling to gauge fields. We review how these couplings affect the cosmological perturbations via a unified study based on the in-in formalism. We then note that, when the inflaton coupling to gauge fields is high enough to lead to interesting observational results, the backreaction of the produced gauge quanta on the inflaton dynamics becomes relevant during the final stage of inflation, and prolongs its duration by about 10 e-foldings. We extend existing results on gravity wave production in these models to account for this late inflationary phase. The strong backreaction phase results in an enhancement of the gravity wave signal at the interferometer scales. As a consequence, the signal is potentially observable at Advanced LIGO/VIRGO for the most natural duration of inflation in such models. Finally, we explicitly compute the axion couplings to gauge fields in string theory construction of axion monodromy inflation and identify cases where they can trigger interesting phenomenological effects.

We compute all N-point primordial curvature correlation functions from inflation at tree-level up to N of order ten or more depending on the choice of parameters. This is achieved for resonant inflationary models in which the inflaton potential has a periodic modulation on top of a slow-roll flat term. These models find a natural UV completion in string theory implementation of axion monodromy. Key to the success of our computation is the observation that gravitational interactions among the perturbations can be neglected, which we argue is justified for any model of inflation with parametrically large non-Gaussianity. We provide a comprehensive review and detailed derivations of known consistency relations for squeezed and collinear limits, and generalize them to any N-point function.

We provide a derivation from first principles of the primordial bispectrum of scalar perturbations produced during inflation driven by a canonically normalized scalar field whose potential exhibits small sinusoidal modulations. A potential of this type has been derived in a class of string theory models of inflation based on axion monodromy. We use this model as a concrete example, but we present our derivations and results for a general slow-roll potential with superimposed modulations. We show analytically that a resonance between the oscillations of the background and the oscillations of the fluctuations is responsible for the production of an observably large non-Gaussian signal. We provide an explicit expression for the shape of this resonant non-Gaussianity. We show that there is essentially no overlap between this shape and the local, equilateral, and orthogonal shapes, and we stress that resonant non-Gaussianity is not captured by the simplest version of the effective field theory of inflation. We hope our analytic expression will be useful to further observationally constrain this class of models.

The primordial non-Gaussian parameter fNL has been shown to be scale-dependent in several models of inflation with a variable speed of sound. Starting from a simple ansatz for a scale-dependent amplitude of the primordial curvature bispectrum for two common phenomenological models of primordial non-Gaussianity, we perform a Fisher matrix analysis of the bispectra of the temperature and polarization of the Cosmic Microwave Background (CMB) radiation and derive the expected constraints on the parameter nNG that quantifies the running of fNL(k) for current and future CMB missions such as WMAP, Planck and CMBPol. We find that CMB information alone, in the event of a significant detection of the non-Gaussian component, corresponding to fNL = 50 for the local model and fNL = 100 for the equilateral model of non-Gaussianity, is able to determine nNG with a 1-sigma uncertainty of Delta nNG = 0.1 and Delta nNG = 0.3, respectively, for the Planck mission. In addition, we consider a Fisher matrix analysis of the galaxy power spectrum to determine the expected constraints on the running parameter nNG for the local model and of the galaxy bispectrum for the equilateral model from future photometric and spectroscopic surveys. We find that, in both cases, large-scale structure observations should achieve results comparable to or even better than those from the CMB, while showing some complementarity due to the different distribution of the non-Gaussian signal over the relevant range of scales. Finally, we compare our findings to the predictions on the amplitude and running of non-Gaussianity of DBI inflation, showing how the constraints on a scale-dependent fNL(k) translate into constraints on the parameter space of the theory.