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

We study the prediction of the dark matter power spectrum at two-loop order in the Effective Field Theory of Large Scale Structures (EFTofLSS) using high precision numerical simulations. In our universe, short distance non-linear perturbations, not under perturbative control, affect long distance fluctuations through an effective stress tensor that needs to be parametrized in terms of counterterms that are functions of the long distance fluctuating fields. We find that at two-loop order it is necessary to include three counterterms: a linear term in the over density, $\delta$, a quadratic term, $\delta^2$, and a higher derivative term, $\partial^2\delta$. After the inclusion of these three terms, the EFTofLSS at two-loop order matches simulation data up to $k\simeq 0.34 \,h\, {\rm Mpc}^{-1}$ at redshift $z=0$, up to $k\simeq 0.55\,h\, {\rm Mpc}^{-1}$ at $z=1$, and up to $k\simeq 1.1\,h\, {\rm Mpc}^{-1}$ at $z=2$. At these wavenumbers, the cosmic variance of the simulation is at least as small as $10^{-3}$, providing a high precision comparison between theory and data. The actual reach of the theory is affected by theoretical uncertainties associated to not having included higher order terms in perturbation theory, for which we provide an estimate, and by potentially overfitting the data, which we also try to address. Since in the EFTofLSS the coupling constants associated with the counterterms are unknown functions of time, we show how a simple parametrization gives a sensible description of their time-dependence. Overall, the $k$-reach of the EFTofLSS is much larger than previous analytical techniques, showing that the amount of cosmological information amenable to high-precision analytical control might be much larger than previously believed.

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

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 discuss the consistency relations involving the soft limit of the (n + 1)-correlator functions of dark matter at equal times and their consequences for the halo model.

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 analyze the class of models where a suitable coupling between the inflaton field and the vector field gives rise to scale-invariant vector perturbations. We exploit the fact that the de Sitter isometry group acts as conformal group on the three-dimensional Euclidean space for the super-Hubble fluctuations in order to characterize the correlators involving the inflaton and the vector fields.

The primordial non-Gaussianity parameters fNL and tauNL may be scale-dependent. We investigate the capability of future measurements of the CMB mu-distortion, which is very sensitive to small scales, and of the large-scale halo bias to test the running of local non-Gaussianity. We show that, for an experiment such as PIXIE, a measurement of the mu-temperature correlation can pin down the spectral indices n_fNL and n_tauNL to values of the order of 0.3 if fNL = 20 and tauNL = 5000. A similar value can be achieved with an all-sky survey extending to redshift z ~ 1. In the particular case in which the two spectral indices are equal, as predicted in models where the cosmological perturbations are generated by a single-field other than the inflaton, then the 1-sigma error on the scale-dependence of the non-linearity parameters goes down to 0.2.

We have recently proposed that the Standard Model Higgs might be responsible for generating the cosmological perturbations of the universe by acting as an isocurvature mode during a de Sitter inflationary stage. In this paper we study the level of non-Gaussianity in the cosmological perturbations which are inevitably generated due to the non-linearities of the Standard Model Higgs potential. In particular, for the current central value of the top mass, we find that a future detection of non-Gaussianity would exclude the detection of tensor modes by the PLANCK satellite.

We consider a model of inflation which has recently been proposed in the literature and where inflation is induced by corrections to the energy density coming from the non-commutativity of spacetime. We show that the very rapid inflationary expansion typical of this model is responsible for a burst of particle production which ends inflation and leads to a radiation-dominated phase. We analytically estimate the energy density of these particles and we confront the results with more precise numerical calculations. We estimate the number of inflationary e-folds before the back-reaction of the radiation energy density overcomes the non-commutative effects and terminate inflation naturally.