Primordial black holes, formed from rare peaks in the primordial fluctuations $\zeta$, are non-rotating at zeroth order in $\zeta_{\rm rms}$. We show that the spin also vanishes at first order in $\zeta_{\rm rms}$, suggesting the dimensionless spin parameter $a_{\rm rms} \sim \zeta_{\rm rms}^2$. We identify one quadratic contribution to the spin by calculating (and extrapolating to the formation time) the torque on a black hole due to ambient acoustic waves. For a reasonable density of primordial black holes this implies a percent level spin parameter.

We write down the Lagrangian bias expansion in general relativity up to 4th order in terms of operators describing the curvature of an early-time hypersurface for comoving observers. They can be easily expanded in synchronous or comoving gauges. This is necessary for the computation of the one-loop halo bispectrum, where relativistic effects can be degenerate with a primordial non-Gaussian signal. Since the bispectrum couples scales, an accurate prediction of the squeezed limit behavior needs to be both non-linear and relativistic. We then evolve the Lagrangian bias operators in time in comoving gauge, obtaining non-local operators analogous to what is known in the Newtonian limit. Finally, we show how to renormalize the bias expansion at an arbitrary time and find that this is crucial in order to cancel unphysical $1/k^2$ divergences in the large-scale power spectrum and bispectrum that could be mistaken for a contamination to the non-Gaussian signal.

The large scale structure bispectrum in the squeezed limit couples large with small scales. Since relativity is important at large scales and non-linear loop corrections are important at small scales, the proper calculation of the observed bispectrum in this limit requires a non-linear relativistic calculation. We compute the matter bispectrum in general relativity in the weak field approximation. The calculation is as involved as existing second-order results. We find several differences with the Newtonian calculation such as the non-cancellation of IR divergences, the need to renormalize the background, and the fact that initial conditions must be set at second order in perturbation theory. For the bispectrum, we find relativistic corrections to be as large as the newtonian result in the squeezed limit. In that limit relativistic one-loop contributions, which we compute for the first time, can be as large as tree level results and have the same $1/k^2$ dependence as a primordial local non-Gaussianity signal where $k$ is the momentum approaching zero. Moreover, we find the time dependence of the relavistic corrections to the bispectrum to be the same as that of a primordial non-Gaussianity signal.

Modes with wavelengths larger than the survey window can have significant impact on the covariance within the survey window. The supersample covariance has been recognized as an important source of covariance for the power spectrum on small scales, and it can potentially be important for the bispectrum covariance as well. In this paper, using the response function formalism, we model the supersample covariance contributions to the bispectrum covariance and the cross covariance between the power spectrum and the bispectrum. The supersample covariances due to the long wavelength density and tidal perturbations are investigated, and the tidal contribution is a few orders of magnitude smaller than the density one because in configuration space the bispectrum estimator involves angular averaging and the tidal response function is anisotropic. The impact of the super-survey modes is quantified using numerical measurements with periodic box and subbox setups. For the matter bispectrum, the ratio between the supersample covariance correction and the small scale covariance, which can be computed using a periodic box, is roughly an order of magnitude smaller than that for the matter power spectrum. This is because for the bispectrum, the small scale non-Gaussian covariance is significantly larger than that for the power spectrum. For the cross covariance, the supersample covariance is as important as for the power spectrum covariance. The supersample covariance prediction with the halo model response function is in good agreement with numerical results.

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.

High-precision constraints on primordial non-Gaussianity (PNG) will significantly improve our understanding of the physics of the early universe. Among all the subtleties in using large scale structure observables to constrain PNG, accounting for relativistic corrections to the clustering statistics is particularly important for the upcoming galaxy surveys covering progressively larger fraction of the sky. We focus on relativistic projection effects due to the fact that we observe the galaxies through the light that reaches the telescope on perturbed geodesics. These projection effects can give rise to an effective $f_{\rm NL}$ that can be misinterpreted as the primordial non-Gaussianity signal and hence is a systematic to be carefully computed and accounted for in modelling of the bispectrum. We develop the technique to properly account for relativistic effects in terms of purely observable quantities, namely angles and redshifts. We give some examples by applying this approach to a subset of the contributions to the tree-level bispectrum of the observed galaxy number counts calculated within perturbation theory and estimate the corresponding non-Gaussianity parameter, $f_{\rm NL}$, for the local, equilateral and orthogonal shapes. For the local shape, we also compute the local non-Gaussianity resulting from terms obtained using the consistency relation for observed number counts. Our goal here is not to give a precise estimate of $f_{\rm NL}$ for each shape but rather we aim to provide a scheme to compute the non-Gaussian contamination due to relativistic projection effects. For the terms considered in this work, we obtain contamination of $f_{\rm NL}^{\rm loc} \sim {\mathcal O}(1)$.

Axions play a central role in many realizations of large field models of inflation and in recent alternative mechanisms for generating primordial tensor modes in small field models. If these axions couple to gauge fields, the coupling produces a tachyonic instability that leads to an exponential enhancement of the gauge fields, which in turn can decay into observable scalar or tensor curvature perturbations. Thus, a fully self-consistent treatment of axions during inflation is important, and in this work we discuss the perturbative constraints on axions coupled to gauge fields. We show how the recent proposal of generating tensor modes through these alternative mechanisms is in tension with perturbation theory in the in-in formalism. Interestingly, we point out that the constraints are parametrically weaker than one would estimate based on naive power counting of propagators of the gauge field. In the case of non-Abelian gauge fields, we derive new constraints on the size of the gauge coupling, which apply also in certain models of natural large field inflation, such as alignment mechanisms.

We study the halo-matter cross bispectrum in the presence of primordial non-Gaussianity of the local type. We restrict ourselves to the squeezed limit, for which the calculation are straightforward, and perform the measurements in the initial conditions of N-body simulations, to mitigate the contamination induced by nonlinear gravitational evolution. Interestingly, the halo-matter cross bispectrum is not trivial even in this simple limit as it is strongly sensitive to the scale-dependence of the quadratic and third-order halo bias. Therefore, it can be used to test biasing prescriptions. We consider three different prescription for halo clustering: excursion set peaks (ESP), local bias and a model in which the halo bias parameters are explicitly derived from a peak-background split. In all cases, the model parameters are fully constrained with statistics other than the cross bispectrum. We measure the cross bispectrum involving one halo fluctuation field and two mass overdensity fields for various halo masses and collapse redshifts. We find that the ESP is in reasonably good agreement with the numerical data, while the other alternatives we consider fail in various cases. This suggests that the scale-dependence of halo bias also is a crucial ingredient to the squeezed limit of the halo bispectrum.

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 obtain a consistency relation for the observed three-point correlator of galaxies. It includes relativistic effects and it is valid in the squeezed limit. Furthermore, the consistency relation is non-perturbative and can be used at arbitrarily small scales for the short modes. Our results are also useful to compute the non-linear relativistic corrections which induce a signal in the observations that might be misinterpreted as primordial non-Gaussianity with a local shape. We estimate the effective local non-Gaussianity parameter from the relativistic corrections. The exact value depends on the redshift and the magnification bias. At redshift of $ z = 1$, in the absence of magnification bias, we get $\,\, f^{\rm GR}_{\rm NL} \simeq - 3.7 $.

We examine the cosmological implications of the measurements of the linear growth rate of cosmological structure obtained in a companion paper from the power spectrum and bispectrum monopoles of the Sloan Digital Sky Survey III Baryon Oscillation Spectroscopic Survey Data, Release 11, CMASS galaxies. This measurement was of $f^{0.43}\sigma_8$, where $\sigma_8$ is the amplitude of dark matter density fluctuations, and $f$ is the linear growth rate, at the effective redshift of the survey, $z_{\rm eff}=0.57$. In conjunction with Cosmic Microwave Background (CMB) data, interesting constraints can be placed on models with non-standard neutrino properties and models where gravity deviates from general relativity on cosmological scales. In particular, the sum of the masses of the three species of the neutrinos is constrained to $m_\nu<0.49\,{\rm eV}$ (at 95\% confidence level) when the $f^{0.43}\sigma_8$ measurement is combined with state-of-the-art CMB measurements. Allowing the effective number of neutrinos to vary as a free parameter does not significantly change these results. When we combine the measurement of $f^{0.43}\sigma_8$ with the complementary measurement of $f\sigma_8$ from the monopole and quadrupole of the two-point correlation function we are able to obtain an independent measurements of $f$ and $\sigma_8$. We obtain $f=0.63\pm0.16$ and $\sigma_8=0.710\pm0.086$ (68\% confidence level). This is the first time when these parameters have been able to be measured independently using the redshift-space power spectrum and bispectrum measurements from galaxy clustering data only.

We analyse the anisotropic clustering of the Baryon Oscillation Spectroscopic Survey (BOSS) CMASS Data Release 11 sample, which consists of $690 827$ galaxies in the redshift range $0.43 < z < 0.70$ and has a sky coverage of $8 498$ deg$^2$ corresponding to an effective volume of $\sim6\,\rm{Gpc}^3$. We fit the Fourier space statistics, the power spectrum and bispectrum monopoles to measure the linear and quadratic bias parameters, $b_1$ and $b_2$, for a non-linear non-local bias model, the growth of structure parameter $f$ and the amplitude of dark matter density fluctuations parametrised by $\sigma_8$. We obtain $b_1(z_{\rm eff})^{1.40}\sigma_8(z_{\rm eff})=1.672\pm 0.060$ and $b_2^{0.30}(z_{\rm eff})\sigma_8(z_{\rm eff})=0.579\pm0.082$ at the effective redshift of the survey, $z_{\rm eff}=0.57$. The main cosmological result is the constraint on the combination $f^{0.43}(z_{\rm eff})\sigma_8(z_{\rm eff})=0.582\pm0.084$, which is complementary to $f\sigma_8$ constraints obtained from 2-point redshift space distortion analyses. A less conservative analysis yields $f^{0.43}(z_{\rm eff})\sigma_8(z_{\rm eff})=0.584\pm0.051$. We ensure that our result is robust by performing detailed systematic tests using a large suite of survey galaxy mock catalogs and N-body simulations. The constraints on $f^{0.43}\sigma_8$ are useful for setting additional constrains on neutrino mass, gravity, curvature as well as the number of neutrino species from galaxy surveys analyses (as presented in a companion paper).

We present a phenomenological modification of the standard perturbation theory prediction for the bispectrum in redshift space that allows us to extend the model to mildly non-linear scales over a wide range of redshifts, $z\leq1.5$. We find that we can describe the bispectrum of dark matter particles with $\sim5%$ accuracy for $k_i\lesssim0.10\,h/{\rm Mpc}$ at $z=0$, for $k_i\lesssim0.15\,h/{\rm Mpc}$ at $z=0.5$, for $k_i\lesssim0.17\,h/{\rm Mpc}$ at $z=1.0$ and for $k_i\lesssim0.20\,h/{\rm Mpc}$ at $z=1.5$. We also test that the fitting formula is able to describe with similar accuracy the bispectrum of cosmologies with different $\Omega_m$, in the range $0.2\lesssim \Omega_m \lesssim 0.4$, and consequently with different values of the logarithmic grow rate $f$ at $z=0$, $0.4\lesssim f(z=0) \lesssim 0.6$. We apply this new formula to recover the bias parameters, $f$ and $\sigma_8$, by combining the redshift space power spectrum monopole and quadrupole with the bispectrum monopole for both dark matter particles and haloes. We find that the combination of these three statistics can break the degeneracy between $b_1$, $f$ and $\sigma_8$. For dark matter particles the new model can be used to recover $f$ and $\sigma_8$ with $\sim1%$ accuracy. For dark matter haloes we find that $f$ and $\sigma_8$ present larger systematic shifts, $\sim10%$. The systematic offsets arise because of limitations in the modelling of the interplay between bias and redshift space distortions, and represent a limitation as the statistical errors of forthcoming surveys reach this level. Conveniently, we find that these residual systematics are mitigated for combinations of parameters. The improvement on the modelling of the bispectrum presented in this paper will be useful for extracting information from current and future galaxy surveys. [abridged]

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.

A large-scale temperature perturbation has a non-zero correlation with the power spectrum of B-modes of cosmological origin on short scales while the corresponding correlation is expected to be zero for B-modes sourced by galactic foregrounds. We thus compute the squeezed limit of a three-point function in which one correlates the temperature fluctuations at large scales with two polarization modes at small scales. In the particular case of the B-mode polarization we obtain a relation that connects the squeezed limit of the $TBB$ three-point function with the cosmological B-mode power spectrum, which can be used as a consistency relation. This could in principle help to distinguish a primordial signal from that induced by inter-stellar dust.

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.

Consistency relations involving the soft limit of the (n + 1)-correlator functions of dark matter and galaxy overdensities can be obtained, both in real and redshift space, thanks to the symmetries enjoyed by the Newtonian equations of motion describing the dark matter and galaxy fluids coupled through gravity. We study the implications of such symmetries for the theory of galaxy bias and for the theories of modified gravity. We find that the invariance of the fluid equations under a coordinate transformation that induces a long-wavelength velocity constrain the bias to depend only on a set of invariants, while the symmetry of such equations under Lifshitz scalings in the case of matter domination allows one to compute the time-dependence of the coefficients in the bias expansion. We also find that in theories of modified gravity which violate the equivalence principle induce a violation of the consistency relation which may be a signature for their observation. Thus, given adiabatic Gaussian initial conditions, the observation of a deviation from the consistency relation for galaxies would signal a break-down of the so-called non-local Eulerian bias model or the violation of the equivalence principle in the underlying theory of gravity.

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 consider whether the non-Gaussian scale-dependent halo bias can be used not only to constrain the local form of non-Gaussianity but also to distinguish among different shapes. In particular, we ask whether it can constrain the behavior of the primordial three-point function in the squeezed limit where one of the momenta is much smaller than the other two. This is potentially interesting since the observation of a three-point function with a squeezed limit that does not go like the local nor equilateral templates would be a signal of non-trivial dynamics during inflation. To this end we use the quasi-single field inflation model of Chen and Wang as a representative two-parameter model, where one parameter governs the amplitude of non-Gaussianity and the other the shape. We also perform a model-independent analysis by parametrizing the scale-dependent bias as a power-law on large scales, where the power is to be constrained from observations. We find that proposed large-scale structure surveys (with characteristics similar to the dark energy task force stage IV surveys) have the potential to distinguish among the squeezed limit behavior of different bispectrum shapes for a wide range of fiducial model parameters. Thus the halo bias can help discriminate between different models of inflation.

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.

We implement Slow Roll Reconstruction -- an optimal solution to the inverse problem for inflationary cosmology -- within ModeCode, a publicly available solver for the inflationary dynamics. We obtain up-to-date constraints on the reconstructed inflationary potential, derived from the WMAP 7-year dataset and South Pole Telescope observations, combined with large scale structure data derived from SDSS Data Release 7. Using ModeCode in conjunction with the MultiNest sampler, we compute Bayesian evidence for the reconstructed potential at each order in the truncated slow roll hierarchy. We find that the data are well-described by the first two slow roll parameters, \epsilon and \eta, and that there is no need to include a nontrivial \xi parameter.

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.

We present a method to minimize, or even cancel out, the nuisance parameters affecting a measurement. Our approach is general and can be applied to any experiment or observation. We compare it with the bayesian technique used to deal with nuisance parameters: marginalization, and show how the method compares and improves by avoiding biases. We illustrate the method with several examples taken from the astrophysics and cosmology world: baryonic acoustic oscillations, cosmic clocks, Supernova Type Ia luminosity distance, neutrino oscillations and dark matter detection. By applying the method we recover some known results but also find some interesting new ones. For baryonic acoustic oscillation (BAO) experiments we show how to combine radial and angular BAO measurements in order to completely eliminate the dependence on the sound horizon at radiation drag. In the case of exploiting SN1a as standard candles we show how the uncertainty in the luminosity distance by a second parameter modeled as a metallicity dependence can be eliminated or greatly reduced. When using cosmic clocks to measure the expansion rate of the universe, we demonstrate how a particular combination of observables nearly removes the metallicity dependence of the galaxy on determining differential ages, thus removing the age-metallicity degeneracy in stellar populations. We hope that these findings will be useful in future surveys to obtain robust constraints on the dark energy equation of state.

The abundance of collapsed objects in the universe, or halo mass function, is an important theoretical tool in studying the effects of primordially generated non-Gaussianities on the large scale structure. The non-Gaussian mass function has been calculated by several authors in different ways, typically by exploiting the smallness of certain parameters which naturally appear in the calculation, to set up a perturbative expansion. We improve upon the existing results for the mass function by combining path integral methods and saddle point techniques (which have been separately applied in previous approaches). Additionally, we carefully account for the various scale dependent combinations of small parameters which appear. Some of these combinations in fact become of order unity for large mass scales and at high redshifts, and must therefore be treated non-perturbatively. Our approach allows us to do this, and to also account for multi-scale density correlations which appear in the calculation. We thus derive an accurate expression for the mass function which is based on approximations that are valid over a larger range of mass scales and redshifts than those of other authors. By tracking the terms ignored in the analysis, we estimate theoretical errors for our result and also for the results of others. We also discuss the complications introduced by the choice of smoothing filter function, which we take to be a top-hat in real space, and which leads to the dominant errors in our expression. Finally, we present a detailed comparison between the various expressions for the mass functions, exploring the accuracy and range of validity of each.

Primordial non-Gaussianity (NG) affects the large scale structure (LSS) of the universe by leaving an imprint on the distribution of matter at late times. Much attention has been focused on using the distribution of collapsed objects (i.e. dark matter halos and the galaxies and galaxy clusters that reside in them) to probe primordial NG. An equally interesting and complementary probe however is the abundance of extended underdense regions or voids in the LSS. The calculation of the abundance of voids using the excursion set formalism in the presence of primordial NG is subject to the same technical issues as the one for halos, which were discussed e.g. in arXiv:1005.1203. However, unlike the excursion set problem for halos which involved random walks in the presence of one barrier $\delta_c$, the void excursion set problem involves two barriers $\delta_v$ and $\delta_c$. This leads to a new complication introduced by what is called the "void-in-cloud" effect discussed in the literature, which is unique to the case of voids. We explore a path integral approach which allows us to carefully account for all these issues, leading to a rigorous derivation of the effects of primordial NG on void abundances. The void-in-cloud issue in particular makes the calculation conceptually rather different from the one for halos. However, we show that its final effect can be described by a simple yet accurate approximation. Our final void abundance function is valid on larger scales than the expressions of other authors, while being broadly in agreement with those expressions on smaller scales.

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

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.