We point out that dark-energy perturbations may become unstable in the presence of a gravitational wave of sufficiently large amplitude. We study this effect for the cubic Horndeski operator (braiding), proportional to $\alpha_{\rm B}$. The scalar that describes dark-energy fluctuations features ghost and/or gradient instabilities for gravitational-wave amplitudes that are produced by typical binary systems. Taking into account the populations of binary systems, we conclude that the instability is triggered in the whole Universe for $|\alpha_{\rm B} |\gtrsim 10^{-2}$, i.e. when the modification of gravity is sizeable. The instability is triggered by massive black-hole binaries down to frequencies corresponding to $10^{10}$ km: the instability is thus robust, unless new physics enters on even longer wavelengths. The fate of the instability and the subsequent time-evolution of the system depend on the UV completion, so that the theory may end up in a state very different from the original one. The same kind of instability is present in beyond-Horndeski theories for $|\alpha_{\rm H}| \gtrsim 10^{-20}$. In conclusion, the only dark-energy theories with sizeable cosmological effects that avoid these problems are $k$-essence models, with a possible conformal coupling with matter.

We study, using Mean Curvature Flow methods, 3+1 dimensional cosmologies with a positive cosmological constant, matter satisfying the dominant and the strong energy conditions, and with spatial slices that can be foliated by 2-dimensional surfaces that are the closed orbits of a symmetry group. If these surfaces have non-positive Euler characteristic (or in the case of 2-spheres, if the initial 2-spheres are large enough) and also if the initial spatial slice is expanding everywhere, then we prove that asymptotically the spacetime becomes physically indistinguishable from de Sitter space on arbitrarily large regions of spacetime. This holds true notwithstanding the presence of initial arbitrarily-large density fluctuations.

Scalar hair of black holes in theories with a shift symmetry are constrained by the no-hair theorem of Hui and Nicolis, assuming spherical symmetry, time-independence of the scalar field and asymptotic flatness. The most studied counterexample is a linear coupling of the scalar with the Gauss-Bonnet invariant. However, in this case the norm of the shift-symmetry current $J^2$ diverges at the horizon casting doubts on whether the solution is physically sound. We show that this is not an issue since $J^2$ is not a scalar quantity, since $J^\mu$ is not a diff-invariant current in the presence of Gauss-Bonnet. The same theory can be written in Horndeski form with a non-analytic function $G_5 \sim \log X$. In this case the shift-symmetry current is diff-invariant, but contains powers of $X$ in the denominator, so that its divergence at the horizon is again immaterial. We confirm that other hairy solutions in the presence of non-analytic Horndeski functions are pathological, featuring divergences of physical quantities as soon as one departs from time-independence and spherical symmetry. We generalise the no-hair theorem to Beyond Horndeski and DHOST theories, showing that the coupling with Gauss-Bonnet is necessary to have hair.

We study, using Mean Curvature Flow methods, 2+1 dimensional cosmologies with a positive cosmological constant and matter satisfying the dominant and the strong energy conditions. If the spatial slices are compact with non-positive Euler characteristic and are initially expanding everywhere, then we prove that the spatial slices reach infinite volume, asymptotically converge on average to de Sitter and they become, almost everywhere, physically indistinguishable from de Sitter. This holds true notwithstanding the presence of initial arbitrarily-large density fluctuations and the formation of black holes.

We study the decay of gravitational waves into dark energy fluctuations $\pi$, taking into account the large occupation numbers. We describe dark energy using the effective field theory approach, in the context of generalized scalar-tensor theories. When the $m_3^3$ (cubic Horndeski) and $\tilde m_4^2$ (beyond Horndeski) operators are present, the gravitational wave acts as a classical background for $\pi$ and modifies its dynamics. In particular, $\pi$ fluctuations are described by a Mathieu equation and feature instability bands that grow exponentially. Focusing on the regime of small gravitational-wave amplitude, corresponding to narrow resonance, we calculate analytically the produced $\pi$, its energy and the change of the gravitational-wave signal. The resonance is affected by $\pi$ self-interactions in a way that we cannot describe analytically. This effect is very relevant for the operator $m_3^3$ and it limits the instability. In the case of the $\tilde m_4^2$ operator self-interactions can be neglected, at least in some regimes. The modification of the gravitational-wave signal is observable for $3 \times 10^{-20} \lesssim \alpha_{\rm H} \lesssim 10^{-17}$ with a LIGO/Virgo-like interferometer and for $10^{-16} \lesssim \alpha_{\rm H} \lesssim 10^{-10}$ with a LISA-like one.

We study the decay of gravitational waves into dark energy fluctuations $\pi$, through the processes $\gamma \to \pi\pi$ and $\gamma \to \gamma \pi$, made possible by the spontaneous breaking of Lorentz invariance. Within the EFT of Dark Energy (or Horndeski/beyond Horndeski theories) the first process is large for the operator $\frac12 \tilde m_4^2(t) \, \delta g^{00}\, \left( {}^{(3)}\! R + \delta K_\mu^\nu \delta K^\mu_\nu -\delta K^2 \right)$, so that the recent observations force $\tilde m_4 =0$ (or equivalently $\alpha_{\rm H}=0$). This constraint, together with the requirement that gravitational waves travel at the speed of light, rules out all quartic and quintic GLPV theories. Additionally, we study how the same couplings affect the propagation of gravitons at loop order. The operator proportional to $\tilde m_4^2$ generates a calculable, non-Lorentz invariant higher-derivative correction to the graviton propagation. The modification of the dispersion relation provides a bound on $\tilde m_4^2$ comparable to the one of the decay. Conversely, operators up to cubic Horndeski do not generate sizeable higher-derivative corrections.

The existence of light particles with spin during inflation is prohibited by the Higuchi bound. This conclusion can be evaded if one considers states with a sizeable coupling with the inflaton foliation, since this breaks the de Sitter isometries. The action for these states can be constructed within the Effective Field Theory of Inflation, or using a CCWZ procedure. Light particles with spin have prescribed couplings with soft inflaton perturbations, which are encoded in consistency relations. We study the phenomenology of light states with spin 2. These mix with the graviton changing the tensor power spectrum and can lead to sizeable tensor non-Gaussianities. They also give rise to a scalar bispectrum and trispectrum with a characteristic angle-dependent non-Gaussianity.

The observation of GW170817 and its electromagnetic counterpart implies that gravitational waves travel at the speed of light, with deviations smaller than a few $\times 10^{-15}$. We discuss the consequences of this experimental result for models of dark energy and modified gravity characterized by a single scalar degree of freedom. To avoid tuning, the speed of gravitational waves must be unaffected not only for our particular cosmological solution, but also for nearby solutions obtained by slightly changing the matter abundance. For this to happen the coefficients of various operators must satisfy precise relations that we discuss both in the language of the Effective Field Theory of Dark Energy and in the covariant one, for Horndeski, beyond Horndeski and degenerate higher-order theories. The simplification is dramatic: of the three functions describing quartic and quintic beyond Horndeski theories, only one remains and reduces to a standard conformal coupling to the Ricci scalar for Horndeski theories. We show that the deduced relations among operators do not introduce further tuning of the models, since they are stable under quantum corrections.

We study generalized disformal transformations, including derivatives of the metric, in the context of the Effective Field Theory of Inflation. All these transformations do not change the late-time cosmological observables but change the coefficients of the operators in the action: some couplings are effectively redundant. At leading order in derivatives and up to cubic order in perturbations, one has 6 free functions that can be used to set to zero 6 of the 17 operators at this order. This is used to show that the tensor three-point function cannot be modified at leading order in derivatives, while the scalar-tensor-tensor correlator can only be modified by changing the scalar dynamics. At higher order in derivatives there are transformations that do not affect the Einstein-Hilbert action: one can find 6 additional transformations that can be used to simplify the inflaton action, at least when the dynamics is dominated by the lowest derivative terms. We also identify the leading higher-derivative corrections to the tensor power spectrum and bispectrum.

We argue that $isotropic$ scalar fluctuations in solid inflation are adiabatic in the super-horizon limit. During the solid phase this adiabatic mode has peculiar features: constant energy-density slices and comoving slices do not coincide, and their curvatures, parameterized respectively by $\zeta$ and $\mathcal R$, both evolve in time. The existence of this adiabatic mode implies that Maldacena's squeezed limit consistency relation holds after angular average over the long mode. The correlation functions of a long-wavelength spherical scalar mode with several short scalar or tensor modes is fixed by the scaling behavior of the correlators of short modes, independently of the solid inflation action or dynamics of reheating.

We propose a technically natural scenario whereby an initially large cosmological constant (c.c.) is relaxed down to the observed value due to the dynamics of a scalar evolving on a very shallow potential. The model crucially relies on a sector that violates the null energy condition (NEC) and gets activated only when the Hubble rate becomes sufficiently small --- of the order of the present one. As a result of NEC violation, this low-energy universe evolves into inflation, followed by reheating and the standard Big Bang cosmology. The symmetries of the theory force the c.c. to be the same before and after the NEC-violating phase, so that a late-time observer sees an effective c.c. of the correct magnitude. Importantly, our model allows neither for eternal inflation nor for a set of possible values of dark energy, the latter fixed by the parameters of the theory.

We study the stability of spatially flat FRW solutions which are geodesically complete, i.e. for which one can follow null (graviton) geodesics both in the past and in the future without ever encountering singularities. This is the case of NEC-violating cosmologies such as smooth bounces or solutions which approach Minkowski in the past. We study the EFT of linear perturbations around a solution of this kind, including the possibility of multiple fields and fluids. One generally faces a gradient instability which can be avoided only if the operator $~^{(3)}{R} \delta N~$ is present and its coefficient changes sign along the evolution. This operator (typical of beyond-Horndeski theories) does not lead to extra degrees of freedom, but cannot arise starting from any theory with second-order equations of motion. The change of sign of this operator prevents to set it to zero with a generalised disformal transformation.

We point out that tensor consistency relations-i.e. the behavior of primordial correlation functions in the limit a tensor mode has a small momentum-are more universal than scalar consistency relations. They hold in the presence of multiple scalar fields and as long as anisotropies are diluted exponentially fast. When de Sitter isometries are approximately respected during inflation this is guaranteed by the Higuchi bound, which forbids the existence of light particles with spin: De Sitter space can support scalar hair but no curly hair. We discuss two indirect ways to look for the violation of tensor con- sistency relations in observations, as a signature of models in which inflation is not a strong isotropic attractor, such as solid inflation: (a) Graviton exchange contribution to the scalar four-point function; (b) Quadrupolar anisotropy of the scalar power spectrum due to super-horizon tensor modes. This anisotropy has a well-defined statistics which can be distinguished from cases in which the background has a privileged direction.

We study how the coupling with gravity of theories with non-linearly realized space-time symmetries is modified when one changes the parametrization of the coset. As an example, we focus on the so-called Galileon duality: a reparametrization which maps a Galilean invariant action into another one which enjoys the same symmetry. Starting with a standard coupling with gravity, with a parametric separation between the Planck scale and the typical scale of the coset, one ends up with a theory without such a separation. In particular an infinite set of higher-dimension operators are relevant when the superluminality of the Galileon is measurable in the effective theory. This addresses an apparent paradox since superluminality arises in the dual theory even when absent in the original one.

We investigate the possible implications of the measured value of the scalar tilt $n_s$ for the tensor-to-scalar ratio $r$ in slow-roll, single-field inflationary models. The measured value of the tilt satisfies $n_s -1\sim 1/N_*$, where $N_* \sim 60$ is the number of $e$-folds for observationally relevant scales. If this is not a coincidence and the scaling holds for different values of $N$, it strongly suggests that either $r$ is as big as $10^{-1}$ (a possibility in tension with the recent data), or smaller than $10^{-2}$ and exponentially dependent on $n_s$. A large region of the ($n_s$,$r$) plane is not compatible with this scaling.

We update the forecasts for the measurement of the tensor-to-scalar ratio $r$ for various ground-based experiments (AdvACT, CLASS, Keck/BICEP3, Simons Array, SPT-3G), balloons (EBEX 10k and Spider) and satellites (CMBPol, COrE and LiteBIRD), taking into account the recent Planck data on polarized dust and using a component separation method. The forecasts do not change significantly with respect to previous estimates when at least three frequencies are available, provided foregrounds can be accurately described by few parameters. We argue that a theoretically motivated goal for future experiments is $r\sim2\times10^{-3}$, and that this is achievable if the noise is reduced to $\sim1\,\mu$K-arcmin and lensing is reduced to $10\%$ in power. We study the constraints experiments will be able to put on the frequency and $\ell$-dependence of the tensor signal as a check of its primordial origin. Futuristic ground-based and balloon experiments can have good constraints on these parameters, even for $r\sim2\times10^{-3}$. For the same value of $r$, satellites will marginally be able to detect the presence of the recombination bump, the most distinctive feature of the primordial signal.

Non-linear realizations of spacetime symmetries can be obtained by a generalization of the coset construction valid for internal ones. The physical equivalence of different representations for spacetime symmetries is not obvious, since their relation involves not only a redefinition of the fields but also a field-dependent change of coordinates. A simple and relevant spacetime symmetry is obtained by the contraction of the 4D conformal group that leads to the Galileon group. We analyze two non-linear realizations of this group, focusing in particular on the propagation of signals around non-trivial backgrounds. The aperture of the lightcone is in general different in the two representations and in particular a free (luminal) massless scalar is mapped in a Galileon theory which admits superluminal propagation. We show that in this theory, if we consider backgrounds that vanish at infinity, there is no asymptotic effect: the displacement of the trajectory integrates to zero, as can be expected since the S-matrix is trivial. Regarding local measurements, we show that the puzzle is solved taking into account that a local coupling with fixed sources in one theory is mapped into a non-local coupling and we show that this effect compensates the different lightcone. Therefore the two theories have a different notion of locality. The same applies to the different non-linear realizations of the conformal group and we study the particular case of a cosmologically interesting background: the Galilean Genesis scenarios.

We show that the prediction for the primordial tensor power spectrum cannot be modified at leading order in derivatives. Indeed, one can always set to unity the speed of propagation of gravitational waves during inflation by a suitable disformal transformation of the metric, while a conformal one can make the Planck mass time-independent. Therefore, the tensor amplitude unambiguously fixes the energy scale of inflation. Using the Effective Field Theory of Inflation, we check that predictions are independent of the choice of frame, as expected. The first corrections to the standard prediction come from two parity violating operators with three derivatives. Also the correlator $<\gamma\gamma\gamma>$ is standard and only receives higher derivative corrections. These results hold also in multifield models of inflation and in alternatives to inflation and make the connection between a (quasi) scale-invariant tensor spectrum and inflation completely robust.

The recently derived consistency relations for Large Scale Structure do not hold if the Equivalence Principle (EP) is violated. We show it explicitly in a toy model with two fluids, one of which is coupled to a fifth force. We explore the constraints that galaxy surveys can set on EP violation looking at the squeezed limit of the 3-point function involving two populations of objects. We find that one can explore EP violations of order 10^{-3} - 10^{-4} on cosmological scales. Chameleon models are already very constrained by the requirement of screening within the Solar System and only a very tiny region of the parameter space can be explored with this method. We show that no violation of the consistency relations is expected in Galileon models.

The simplest inflationary model $V=\frac12 m^2\phi^2$ represents the benchmark for future constraints. For a quadratic potential, the quantity $(n_s-1)+r/4+11 (n_s-1)^2/24$ vanishes (up to corrections which are cubic in slow roll) and can be used to parametrize small deviations from the minimal scenario. Future constraints on this quantity will be able to distinguish a quadratic potential from a pseudo-Nambu-Goldstone boson with $f \lesssim 30 M_{\rm pl}$ and set limits on the deviation from unity of the speed of sound $| c_s-1| \lesssim 3\times 10^{-2}$ (corresponding to an energy scale $\Lambda\gtrsim 2\times 10^{16}\, \mathrm{GeV}$), and on the contribution of a second field to perturbations ($\lesssim 6 \times 10^{-2}$). The limiting factor for these bounds will be the uncertainty on the spectral index. The error on the number of e-folds will be $\Delta N \simeq 0.4$, corresponding to an error on the reheating temperature $\Delta T_\mathrm{rh}/T_\mathrm{rh}\simeq 1.2$. We comment on the relevance of non-Gaussianity after BICEP2 results.

In the simplest inflationary model $V=\frac12 m^2\phi^2$, we provide a prediction accurate up to $1\%$ for the spectral index $n_s$ and the tensor-to-scalar ratio $r$ assuming instantaneous reheating and a standard thermal history: $n_s = 0.9668\pm0.0003$ and $r=0.131\pm 0.001$. This represents the simplest and most informative point in the $(n_s,r)$ plane. The result is independent of the details of reheating (or preheating) provided the conversion to radiation is sufficiently fast. A slower reheating or a modified post-inflationary evolution push towards smaller $n_s$ (and larger $r$), so that our prediction corresponds to the maximum $n_s$ (and minimum $r$) for the quadratic potential. We also derive similar results for a general $V \propto \phi^p$ potential.

We generalize the recently derived single-field consistency relations of Large Scale Structure in two directions. First, we treat the effect of the long modes (with momentum q) on the short ones (with momentum k) non-perturbatively, by writing resummed consistency relations which do not require k/q \delta_q << 1. These relations do not make any assumptions on the short-scales physics and are extended to include (an arbitrary number of) multiple long modes, internal lines with soft momenta and soft loops. We do several checks of these relations in perturbation theory and we verify that the effect of soft modes always cancels out in equal-time correlators. Second, we write the relations directly in redshift space, without assuming the single-stream approximation: not only the long mode affects the short scales as a homogeneous gravitational field, but it also displaces them by its velocity along the line-of-sight. Redshift space consistency relations still vanish when short modes are taken at equal time: an observation of a signal in the squeezed limit would point towards multifield inflation or a violation of the equivalence principle.

In DBI inflation the cubic action is a particular linear combination of the two, otherwise independent, cubic operators \dot \pi^3 and \dot \pi (\partial_i \pi)^2. We show that in the Effective Field Theory (EFT) of inflation this is a consequence of an approximate 5D Poincar\'e symmetry, ISO(4,1), non-linearly realized by the Goldstone \pi. This symmetry uniquely fixes, at lowest order in derivatives, all correlation functions in terms of the speed of sound c_s. In the limit c_s \to 1, the ISO(4,1) symmetry reduces to the Galilean symmetry acting on \pi. On the other hand, we point out that the non-linear realization of SO(4,2), the isometry group of 5D AdS space, does not fix the cubic action in terms of c_s.

We derive consistency relations for the late universe (CDM and \Lambda CDM): relations between an n-point function of the density contrast \delta and an (n+1)-point function in the limit in which one of the (n+1) momenta becomes much smaller than the others. These are based on the observation that a long mode, in single-field models of inflation, reduces to a diffeomorphism since its freezing during inflation all the way until the late universe, even when the long mode is inside the horizon (but out of the sound horizon). These results are derived in Newtonian gauge, at first and second order in the small momentum q of the long mode and they are valid non-perturbatively in the short-scale \delta. In the non-relativistic limit our results match with (Kehagias and Riotto '13) and (Peloso and Pietroni '13). These relations are a consequence of diffeomorphism invariance; they are not satisfied in the presence of extra degrees of freedom during inflation or violation of the Equivalence Principle (extra forces) in the late universe.

In single-field models of inflation the effect of a long mode with momentum q reduces to a diffeomorphism at zeroth and first order in q. This gives the well-known consistency relations for the n-point functions. At order q^2 the long mode has a physical effect on the short ones, since it induces curvature, and we expect that this effect is the same as being in a curved FRW universe. In this paper we verify this intuition in various examples of the three-point function, whose behaviour at order q^2 can be written in terms of the power spectrum in a curved universe. This gives a simple alternative understanding of the level of non-Gaussianity in single-field models. Non-Gaussianity is always parametrically enhanced when modes freeze at a physical scale k_{ph, f} shorter than H: f_{NL} \sim (k_{ph, f}/H)^2.

There are two common non-linear realizations of the 4D conformal group: in the first, the dilaton is the conformal factor of the effective metric \eta_{\mu\nu} e^{-2 \pi}; in the second it describes the fluctuations of a brane in AdS_5. The two are related by a complicated field redefinition, found by Bellucci, Ivanov and Krivonos (2002) to all orders in derivatives. We show that this field redefinition can be understood geometrically as a change of coordinates in AdS_5. In one gauge the brane is rigid at a fixed radial coordinate with a conformal factor on the AdS_5 boundary, while in the other one the brane bends in an unperturbed AdS_5. This geometrical picture illuminates some aspects of the mapping between the two representations. We show that the conformal Galileons in the two representations are mapped into each other in a quite non-trivial way: the DBI action, for example, is mapped into a complete linear combination of all the five Galileons in the other representation. We also verify the equivalence of the dilaton S-matrix in the two representations and point out that the aperture of the dilaton light-cone around non-trivial backgrounds is not the same in the two representations.

We systematically derive the consistency relations associated to the non-linearly realized symmetries of theories with spontaneously broken conformal symmetry but with a linearly-realized de Sitter subalgebra. These identities relate (N+1)-point correlation functions with a soft external Goldstone to N-point functions. These relations have direct implications for the recently proposed conformal mechanism for generating density perturbations in the early universe. We study the observational consequences, in particular a novel one-loop contribution to the four-point function, relevant for the stochastic scale-dependent bias and CMB mu-distortion.

We study the possibility that the approximate time shift symmetry during inflation is promoted to the full invariance under time reparametrization t \to \tilde t(t), or equivalently under field redefinition of the inflaton \phi \to \tilde\phi(\phi). The symmetry allows only two operators at leading order in derivatives, so that all n-point functions of scalar perturbations are fixed in terms of the power spectrum normalization and the speed of sound. During inflation the decaying mode only decays as 1/a and this opens up the possibility to violate some of the consistency relations in the squeezed limit, although this violation is suppressed by the (small) breaking of the field reparametrization symmetry. In particular one can get terms in the 3-point function that are only suppressed by 1/k_L in the squeezed limit k_L \to 0 compared to the local shape.

We generalize the single-field consistency relations to capture not only the leading term in the squeezed limit---going as 1/q^3, where q is the small wavevector---but also the subleading one, going as 1/q^2. This term, for an (n+1)-point function, is fixed in terms of the variation of the n-point function under a special conformal transformation; this parallels the fact that the 1/q^3 term is related with the scale dependence of the n-point function. For the squeezed limit of the 3-point function, this conformal consistency relation implies that there are no terms going as 1/q^2. We verify that the squeezed limit of the 4-point function is related to the conformal variation of the 3-point function both in the case of canonical slow-roll inflation and in models with reduced speed of sound. In the second case the conformal consistency conditions capture, at the level of observables, the relation among operators induced by the non-linear realization of Lorentz invariance in the Lagrangian. These results mean that, in any single-field model, primordial correlation functions of \zeta are endowed with an SO(4,1) symmetry, with dilations and special conformal transformations non-linearly realized by \zeta. We also verify the conformal consistency relations for any n-point function in models with a modulation of the inflaton potential, where the scale dependence is not negligible. Finally, we generalize (some of) the consistency relations involving tensors and soft internal momenta.

In inflationary models where the source of scalar perturbations is not the inflaton, but one or more scalars with negligible coupling with the inflaton, the resulting perturbations are not only scale invariant, but fully conformally invariant with conformal dimension close to zero. This is closely related to the fact that correlation functions can only depend on the de Sitter invariant distances. These properties follow from the isometries of the inflationary de Sitter space and are thus completely independent of the dynamics. The 3-point function is fixed in terms of two constants, while the 4-point function is a function of two parameters (instead of five as in the absence of conformal invariance). The conformal invariance of correlators can be directly checked in Fourier space, as we show in an explicit example. A detection of a non-conformal correlation function, for example an equilateral 3-point function, would imply that the source of perturbations is not decoupled from the inflaton.

We prove that, in a generic single-field model, the consistency relation for the 3-point function in the squeezed limit receives corrections that vanish quadratically in the ratio of the momenta, i.e. as (k_L/k_S)^2. This implies that a detection of a bispectrum signal going as 1/k_L^2 in the squeezed limit, that is suppressed only by one power of k_L compared with the local shape, would rule out all single-field models. The absence of this kind of terms in the bispectrum holds also for multifield models, but only if all the fields have a mass much smaller than H. The detection of any scale dependence of the bias, for scales much larger than the size of the haloes, would disprove all single-field models. We comment on the regime of squeezing that can be probed by realistic surveys.

The CMB bispectrum generated by second-order effects at recombination can be calculated analytically when one of the three modes has a wavelength much longer than the other two and is outside the horizon at recombination. This was pointed out in \cite{Creminelli:2004pv} and here we correct their results. We derive a simple formula for the bispectrum, $f_{NL}^{loc} = - (1/6+ \cos 2 \theta) \cdot (1- 1/2 \cdot d \ln (l_S^2 C_{S})/d \ln l_S)$, where $C_S$ is the short scale spectrum and $\theta$ the relative orientation between the long and the short modes. This formula is exact and takes into account all effects at recombination, including recombination-lensing, but neglects all late-time effects such as ISW-lensing. The induced bispectrum in the squeezed limit is small and will negligibly contaminate the Planck search for a local primordial signal: this will be biased only by $f_{NL}^{loc}\approx-0.4$. The above analytic formula includes the primordial non-Gaussianity of any single-field model. It also represents a consistency check for second-order Boltzmann codes: we find substantial agreement with the CMBquick code.

We study the consequences of imposing an approximate Galilean symmetry on the Effective Theory of Inflation, the theory of small perturbations around the inflationary background. This approach allows us to study the effect of operators with two derivatives on each field, which can be the leading interactions due to non-renormalization properties of the Galilean Lagrangian. In this case cubic non-Gaussianities are given by three independent operators, containing up to six derivatives, two with a shape close to equilateral and one peaking on flattened isosceles triangles. The four-point function is larger than in models with small speed of sound and potentially observable with the Planck satellite.

We propose a novel cosmological scenario, in which standard inflation is replaced by an expanding phase with a drastic violation of the Null Energy Condition (NEC): \dot H >> H^2. The model is based on the recently introduced Galileon theories, that allow NEC violating solutions without instabilities. The unperturbed solution describes a Universe that is asymptotically Minkowski in the past, expands with increasing energy density until it exits the regime of validity of the effective field theory and reheats. This solution is a dynamical attractor and the Universe is driven to it, even if it is initially contracting. The study of perturbations of the Galileon field reveals some subtleties, related to the gross violation of the NEC and it shows that adiabatic perturbations are cosmologically irrelevant. The model, however, suggests a new way to produce a scale invariant spectrum of isocurvature perturbations, which can later be converted to adiabatic: the Galileon is forced by symmetry to couple to the other fields as a dilaton; the effective metric it yields on the NEC violating solution is that of de Sitter space, so that all light scalars will automatically acquire a nearly scale-invariant spectrum of perturbations.

We study the spherical collapse model in the presence of quintessence with negligible speed of sound. This case is particularly motivated for w<-1 as it is required by stability. As pressure gradients are negligible, quintessence follows dark matter during the collapse. The spherical overdensity behaves as a separate closed FLRW universe, so that its evolution can be studied exactly. We derive the critical overdensity for collapse and we use the extended Press-Schechter theory to study how the clustering of quintessence affects the dark matter mass function. The effect is dominated by the modification of the linear dark matter growth function. A larger effect occurs on the total mass function, which includes the quintessence overdensities. Indeed, here quintessence constitutes a third component of virialized objects, together with baryons and dark matter, and contributes to the total halo mass by a fraction ~ (1+w) Omega_Q / Omega_m. This gives a distinctive modification of the total mass function at low redshift.

We calculate the Cosmic Microwave Background anisotropy bispectrum on large angular scales in the absence of primordial non-Gaussianities, assuming exact matter dominance and extending at second order the classic Sachs-Wolfe result \delta T/T=\Phi/3. The calculation is done in Poisson gauge. Besides intrinsic contributions calculated at last scattering, one must consider integrated effects. These are associated to lensing, and to the time dependence of the potentials (Rees-Sciama) and of the vector and tensor components of the metric generated at second order. The bispectrum is explicitly computed in the flat-sky approximation. It scales as l^(-4) in the scale invariant limit and the shape dependence of its various contributions is represented in 3d plots. Although all the contributions to the bispectrum are parametrically of the same order, the full bispectrum is dominated by lensing. In the squeezed limit it corresponds to f_NL^local = -1/6 - cos(2 \theta), where \theta is the angle between the short and the long modes; the angle dependent contribution comes from lensing. In the equilateral limit it corresponds to f_NL^equil ~ 3.13.

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

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

Modern cosmology has sharpened questions posed for millennia about the origin of our cosmic habitat. The age-old questions have been transformed into two pressing issues primed for attack in the coming decade: How did the Universe begin? and What physical laws govern the Universe at the highest energies? The clearest window onto these questions is the pattern of polarization in the Cosmic Microwave Background (CMB), which is uniquely sensitive to primordial gravity waves. A detection of the special pattern produced by gravity waves would be not only an unprecedented discovery, but also a direct probe of physics at the earliest observable instants of our Universe. Experiments which map CMB polarization over the coming decade will lead us on our first steps towards answering these age-old questions.

We study generic single-field dark energy models, by a parametrization of the most general theory of their perturbations around a given background, including higher derivative terms. In appropriate limits this approach reproduces standard quintessence, k-essence and ghost condensation. We find no general pathology associated to an equation of state w_Q < -1 or in crossing the phantom divide w_Q = -1. Stability requires that the w_Q < -1 side of dark energy behaves, on cosmological scales, as a k-essence fluid with a virtually zero speed of sound. This implies that one should set the speed of sound to zero when comparing with data models with w_Q < -1 or crossing the phantom divide. We summarize the theoretical and stability constraints on the quintessential plane (1+w_Q) vs. speed of sound squared.

For slow-roll inflation we study the phase transition to the eternal regime. Starting from a finite inflationary volume, we consider the volume of the universe at reheating as order parameter. We show that there exists a critical value for the classical inflaton speed, \dot\phi^2/H^4 = 3/(2 \pi^2), where the probability distribution for the reheating volume undergoes a sharp transition. In particular, for sub-critical inflaton speeds all distribution moments become infinite. We show that at the same transition point the system develops a non-vanishing probability of having a strictly infinite reheating volume, while retaining a finite probability for finite values. Our analysis represents the exact quantum treatment of the system at lowest order in the slow-roll parameters and H^2/M_Pl^2.

We study nonlinear cosmological perturbations during the post-inflationary evolution, using the equivalence between a perfect barotropic fluid and a derivatively coupled scalar field with Lagrangian [-(\partial \phi)^2]^[(1+w)/2w]. Since this Lagrangian is just a special case of k-inflation, this approach is analogous to the one employed in the study of non-Gaussianities from inflation. We use this method to derive the second order metric during matter dominance in the comoving gauge directly as function of the primordial inflationary perturbation \zeta. Going to Poisson gauge, we recover the metric previously derived in the literature.

We study the effective field theory of inflation, i.e. the most general theory describing the fluctuations around a quasi de Sitter background, in the case of single field models. The scalar mode can be eaten by the metric by going to unitary gauge. In this gauge, the most general theory is built with the lowest dimension operators invariant under spatial diffeomorphisms, like g^{00} and K_{mu nu}, the extrinsic curvature of constant time surfaces. This approach allows us to characterize all the possible high energy corrections to simple slow-roll inflation, whose sizes are constrained by experiments. Also, it describes in a common language all single field models, including those with a small speed of sound and Ghost Inflation, and it makes explicit the implications of having a quasi de Sitter background. The non-linear realization of time diffeomorphisms forces correlation among different observables, like a reduced speed of sound and an enhanced level of non-Gaussianity.

We present a bouncing cosmology which evolves from the contracting to the expanding phase in a smooth way, without developing instabilities or pathologies and remaining in the regime of validity of 4d effective field theory. A nearly scale invariant spectrum of perturbations is generated during the contracting phase by an isocurvature scalar with a negative exponential potential and then converted to adiabatic. The model predicts a slightly blue spectrum, n_S >~ 1, no observable gravitational waves and a high (but model dependent) level of non-Gaussianities with local shape. The model represents an explicit and predictive alternative to inflation, although, at present, it is clearly less compelling.

We study the Likelihood function of data given f_NL for the so-called local type of non-Gaussianity. In this case the curvature perturbation is a non-linear function, local in real space, of a Gaussian random field. We compute the Cramer-Rao bound for f_NL and show that for small values of f_NL the 3-point function estimator saturates the bound and is equivalent to calculating the full Likelihood of the data. However, for sufficiently large f_NL, the naive 3-point function estimator has a much larger variance than previously thought. In the limit in which the departure from Gaussianity is detected with high confidence, error bars on f_NL only decrease as 1/ln Npix rather than Npix^-1/2 as the size of the data set increases. We identify the physical origin of this behavior and explain why it only affects the local type of non-Gaussianity, where the contribution of the first multipoles is always relevant. We find a simple improvement to the 3-point function estimator that makes the square root of its variance decrease as Npix^-1/2 even for large f_NL, asymptotically approaching the Cramer-Rao bound. We show that using the modified estimator is practically equivalent to computing the full Likelihood of f_NL given the data. Thus other statistics of the data, such as the 4-point function and Minkowski functionals, contain no additional information on f_NL. In particular, we explicitly show that the recent claims about the relevance of the 4-point function are not correct. By direct inspection of the Likelihood, we show that the data do not contain enough information for any statistic to be able to constrain higher order terms in the relation between the Gaussian field and the curvature perturbation, unless these are orders of magnitude larger than the size suggested by the current limits on f_NL.

We analyze the 3-year WMAP data and look for a deviation from Gaussianity in the form of a 3-point function that has either of the two theoretically motivated shapes: local and equilateral. There is no evidence of departure from Gaussianity and the analysis gives the presently tightest bounds on the parameters f_ NL^local and f_NL^equil., which define the amplitude of respectively the local and the equilateral non-Gaussianity: -36 < f_NL^local < 100, -256 < f_NL^equil. < 332 at 95% C.L.

We present a consistent effective theory that violates the null energy condition (NEC) without developing any instabilities or other pathological features. The model is the ghost condensate with the global shift symmetry softly broken by a potential. We show that this system can drive a cosmological expansion with dH/dt > 0. Demanding the absence of instabilities in this model requires dH/dt <~ H^2. We then construct a general low-energy effective theory that describes scalar fluctuations about an arbitrary FRW background, and argue that the qualitative features found in our model are very general for stable systems that violate the NEC. Violating the NEC allows dramatically non-standard cosmological histories. To illustrate this, we construct an explicit model in which the expansion of our universe originates from an asymptotically flat state in the past, smoothing out the big-bang singularity within control of a low-energy effective theory. This gives an interesting alternative to standard inflation for solving the horizon problem. We also construct models in which the present acceleration has w < -1; a periodic ever-expanding universe and a model with a smooth ``bounce'' connecting a contracting and expanding phase.

We study the possibility that cosmological density perturbations are generated by the inhomogeneous decay of right-handed neutrinos. This will occur if a scalar field whose fluctuations are created during inflation is coupled to the neutrino sector. Robust predictions of the model are a detectable level of non-Gaussianity and, if standard leptogenesis is the source of the baryon asymmetry, a baryon isocurvature perturbations at the level of the present experimental constraints.

We develop a method to constrain the level of non-Gaussianity of density perturbations when the 3-point function is of the "equilateral" type. Departures from Gaussianity of this form are produced by single field models such as ghost or DBI inflation and in general by the presence of higher order derivative operators in the effective Lagrangian of the inflaton. We show that the induced shape of the 3-point function can be very well approximated by a factorizable form, making the analysis practical. We also show that, unless one has a full sky map with uniform noise, in order to saturate the Cramer-Rao bound for the error on the amplitude of the 3-point function, the estimator must contain a piece that is linear in the data. We apply our technique to the WMAP data obtaining a constraint on the amplitude f_NL^equil of "equilateral" non-Gaussianity: -366 < f_NL^equil < 238 at 95% C.L. We also apply our technique to constrain the so-called "local" shape, which is predicted for example by the curvaton and variable decay width models. We show that the inclusion of the linear piece in the estimator improves the constraint over those obtained by the WMAP team, to -27 < f_NL^local < 121 at 95% C.L.

In the context of Lorentz-invariant massive gravity we show that classical solutions around heavy sources are plagued by ghost instabilities. The ghost shows up in the effective field theory at huge distances from the source, much bigger than the Vainshtein radius. Its presence is independent of the choice of the non-linear terms added to the Fierz-Pauli Lagrangian. At the Vainshtein radius the mass of the ghost is of order of the inverse radius, so that the theory cannot be trusted inside this region, not even at the classical level.

For bouncing cosmologies such as the ekpyrotic/cyclic scenarios we show that it is possible to make predictions for density perturbations which are independent of the details of the bouncing phase. This can be achieved, as in inflationary cosmology, thanks to the existence of a dynamical attractor, which makes local observables equal to the unperturbed solution up to exponentially small terms. Assuming that the physics of the bounce is not extremely sensitive to these corrections, perturbations can be evolved even at non-linear level. The resulting spectrum is not scale invariant and thus incompatible with experimental data. This can be explicitly shown in synchronous gauge where, contrary to what happens in the commonly used Newtonian gauge, all perturbations remain small going towards the bounce and the existence of the attractor is manifest.

We point out the existence of a consistency relation involving the 3-point function of scalar perturbations which is valid in any inflationary model, independently of the inflaton Lagrangian under the assumption that the inflaton is the only dynamical field. The 3-point function in the limit in which one of the momenta is much smaller than the other two is fixed in terms of the power spectrum and its tilt. This relation, although very hard to verify experimentally, could be easily proved wrong by forecoming data, thus ruling out any scenario with a single dynamical field in a model independent way.

We study the 3-point functions generated at recombination in the squeezed triangle limit, when one mode has a wavelength much larger than the other two and is outside the horizon. The presence of the long wavelength mode cannot change the physics inside the horizon but modifies how a late time observer sees the anisotropies. The effect of the long wavelength mode can be divided into a redefinition of time and spatial scales, a Shapiro time delay and gravitational lensing. The separation is gauge dependent but helps develop intuition. We show that the resulting 3-point function corresponds to an f_NL < 1 and that its shape is different from that created by the f_NL (or local) model.

We study the dependence on configuration in momentum space of the primordial 3-point function of density perturbations in several different scenarios: standard slow-roll inflation, curvaton and variable decay models, ghost inflation, models with higher derivative operators and the DBI model of inflation. We define a cosine between the distributions using a measure based on the ability of experiments to distinguish between them. We find that models fall into two broad categories with fairly orthogonal distributions. Models where non-Gaussianity is created at horizon-crossing during inflation and models in which the evolution outside the horizon dominates. In the first case the 3-point function is largest for equilateral triangles, while in the second the dominant contribution to the signal comes from the influence of long wavelength modes on small wavelength ones. We show that, because the distributions in these two cases are so different, translating constraints on parameters of one model to those of another based on the normalization of the 3-point function for equilateral triangles can be very misleading.

We propose a new scenario for early cosmology, where an inflationary de Sitter phase is obtained with a ghost condensate. The transition to radiation dominance is triggered by the ghost itself, without any slow-roll potential. Density perturbations are generated by fluctuations around the ghost condensate and can be reliably computed in the effective field theory. The fluctuations are scale invariant as a consequence of the de Sitter symmetries, however, the size of the perturbations are parametrically different from conventional slow-roll inflation, and the inflation happens at far lower energy scales. The model makes definite predictions that distinguish it from standard inflation, and can be sharply excluded or confirmed by experiments in the near future. The tilt in the scalar spectrum is predicted to vanish (n_s=1), and the gravity wave signal is negligible. The non-Gaussianities in the spectrum are predicted to be observable: the 3-point function is determined up to an overall order 1 constant, and its magnitude is much bigger than in conventional inflation, with an equivalent f_NL ~ 100, not far from the present WMAP bounds.

We study the impact of higher dimension operators in the inflaton Lagrangian on the non-gaussianity of the scalar spectrum. These terms can strongly enhance the effect without spoiling slow-roll, though it is difficult to exceed f_NL ~ 1, because the scale which suppresses the operators cannot be too low, if we want the effective field theory description to make sense. In particular we explicitly calculate the 3-point function given by an higher derivative interaction of the form (\nabla\phi)^4, which is expected to give the most important contribution. The angular dependence of the result turns out to be quite different from the minimal case without higher dimension operators.

We present a new model of inflation in which the inflaton is the extra component of a gauge field in a 5d theory compactified on a circle. The chief merit of this model is that the potential comes only from non-local effects so that its flatness is not spoiled by higher dimensional operators or quantum gravity corrections. The model predicts a red spectrum (n ~ 0.96) and a significant production of gravitational waves (r ~ 0.11). We also comment on the relevance of this idea to quintessence.

We study how to obtain a sufficiently flat inflaton potential that is natural from the particle physics point of view. Supersymmetry, which is broken during inflation, cannot protect the potential against non-renormalizable operators violating slow-roll. We are therefore led to consider models based on non-linearly realized symmetries. The basic scenario with a single four-dimensional pseudo Nambu Goldstone boson requires the spontaneous breaking scale to be above the Planck scale, which is beyond the range of validity of the field theory description, so that quantum gravity corrections are not under control. A nice way to obtain consistent models with large field values is to consider simple extensions in extra-dimensional setups. We also consider the minimal structures necessary to obtain purely four-dimensional models with spontaneous breaking scale below M_P; we show that they require an approximate symmetry that is supplemented by either the little-Higgs mechanism or supersymmetry to give trustworthy scenarios.