Fluctuations with wavelengths larger than the volume of a galaxy survey affect the measurement of the galaxy power spectrum within the survey itself. In the presence of local Primordial Non- Gaussianities (PNG), in addition to the super-sample matter density and tidal fluctuations, the large-scale gravitational potential also induces a modulation of the observed power spectrum. In this work we investigate this modulation by computing for the first time the response of the redshift-space galaxy power spectrum to the presence of a long wavelength gravitational potential, fully accounting for the stochastic contributions. For biased tracers new response functions arise due to couplings between the small-scale fluctuations in the density, velocity and gravitational fields, the latter through scale dependent bias operators, and the large-scale gravitational potential. We study the impact of the super-sample modes on the measurement of the amplitude of the primordial bispectrum of the local-shape, $f_{\rm NL}^{\rm loc}$, accounting for modulations of both the signal and the covariance of the galaxy power spectrum by the long modes. Considering DESI-like survey specifications, we show that in most cases super-sample modes cause little or no degradation of the constraints, and could actually reduce the errorbars on $f_{\rm NL}^{\rm loc}$ by (10 - 30)\%, if external information on the bias parameters is available.

Future large-scale structure surveys will measure three-point statistics with high statistical significance. This will offer significant improvements on our understanding of gravity, provided we can model these statistics accurately. We assess the performance of several schemes for theoretical modelling of the matter bispectrum, including halo-model based approaches and fitting formulae. We compare the model predictions against N-body simulations, considering scales up to $k_{\rm max} = 4 h/{\rm Mpc}$, well into non-linear regime of structure formation. Focusing on the equilateral configuration, we conduct this analysis for three theories of gravity: general relativity, $f(R)$ gravity, and the DGP braneworld model. Additionally, we compute the lensing convergence bispectrum for these models. We find that all current modelling prescriptions in modified gravity, in particular for theories with scale-dependent linear growth, fail to attain the accuracy required by the precision of the Stage IV surveys such as \emph{Euclid}. Among these models, we find that a halo-model corrected fitting formula achieves the best overall performance.

The forthcoming generation of wide-field galaxy surveys will probe larger volumes and galaxy densities, thus allowing for a much larger signal-to-noise ratio for higher-order clustering statistics, in particular the galaxy bispectrum. Extracting this information, however, is more challenging than using the power spectrum due to more complex theoretical modeling, as well as significant computational cost of evaluating the bispectrum signal and the error budget. To overcome these challenges, several proxy statistics have been proposed in the literature, which partially or fully capture the information in the bispectrum, while being computationally less expensive than the bispectrum. One such statistics are {\it weighted skew-spectra}, which are cross-spectra of the density field and appropriately weighted quadratic fields. Using Fisher forecasts, we show that the information in these skew-spectra is equivalent to that in the bispectrum for parameters that appear as amplitudes in the bispectrum model, such as galaxy bias parameters or the amplitude of primordial non-Gaussianity. We consider three shapes of the primordial bispectrum: local, equilateral and that due to massive particles with spin two during inflation. To obtain constraints that match those from a measurement of the full bispectrum, we find that it is crucial to account for the full covariance matrix of the skew-spectra.

This paper discusses the science case for a sensitive spectro-polarimetric survey of the microwave sky. Such a survey would provide a tomographic and dynamic census of the three-dimensional distribution of hot gas, velocity flows, early metals, dust, and mass distribution in the entire Hubble volume, exploit CMB temperature and polarisation anisotropies down to fundamental limits, and track energy injection and absorption into the radiation background across cosmic times by measuring spectral distortions of the CMB blackbody emission. In addition to its exceptional capability for cosmology and fundamental physics, such a survey would provide an unprecedented view of microwave emissions at sub-arcminute to few-arcminute angular resolution in hundreds of frequency channels, a data set that would be of immense legacy value for many branches of astrophysics. We propose that this survey be carried-out with a large space mission featuring a broad-band polarised imager and a moderate resolution spectro-imager at the focus of a 3.5m aperture telescope actively cooled to about 8K, complemented with absolutely-calibrated Fourier Transform Spectrometer modules observing at degree-scale angular resolution in the 10-2000 GHz frequency range. We propose two observing modes: a survey mode to map the entire sky as well as a few selected wide fields, and an observatory mode for deeper observations of regions of specific interest.

This paper outlines the science case for line-intensity mapping with a space-borne instrument targeting the sub-millimeter (microwaves) to the far-infrared (FIR) wavelength range. Our goal is to observe and characterize the large-scale structure in the Universe from present times to the high redshift Epoch of Reionization. This is essential to constrain the cosmology of our Universe and form a better understanding of various mechanisms that drive galaxy formation and evolution. We argue that the proposed frequency range would make it possible to probe important metal cooling lines such as [CII] up to very high redshift as well as a large number of rotational lines of the CO molecule. These can be used to trace molecular gas and dust evolution and constrain the buildup in both the cosmic star formation rate density and the cosmic infrared background (CIB). Moreover, surveys at the highest frequencies will detect FIR lines which are used as diagnostics of galaxies and AGN. Tomography of these lines over a wide redshift range will enable invaluable measurements of the cosmic expansion history at epochs inaccessible to other methods, competitive constraints on the parameters of the standard model of cosmology, and numerous tests of dark matter, dark energy, modified gravity and inflation. To reach these goals, large-scale structure must be mapped over a wide range in frequency to trace its time evolution over a reasonable fraction of the volume of the observable Universe. In addition, the surveyed area needs to be very large to beat cosmic variance and to probe the largest scales where its easier to separate the astrophysical and cosmological contributions to the observed signal. Only, a space-borne mission can properly meet these requirements.

A common mechanism to form primordial black holes in the early universe is by enhancing at small-scales the scalar perturbations generated during inflation. If these fluctuations have a large enough amplitude, they may collapse into primordial black holes upon horizon re-entry. Such primordial black holes may comprise the totality of the dark matter. We offer some considerations about the formation and clustering of primordial black holes when the scalar perturbations are characterised by a broad spectrum. Using the excursion set method, as well as the supreme statistics, we show that the cloud-in-cloud phenomenon, for which small mass primordial black holes may be absorbed by bigger mass ones, is basically absent. This is due to the fact that the formation of a primordial black hole is an extremely rare event. We also show that, from the point of view of mass distribution, broad and narrow spectra give similar results in the sense that the mass distribution is tilted towards a single mass. Furthermore, we argue that primordial black holes from broad spectra are not clustered at formation, their distribution is dominantly Poissonian.

Our current understanding of the Universe is established through the pristine measurements of structure in the cosmic microwave background (CMB) and the distribution and shapes of galaxies tracing the large scale structure (LSS) of the Universe. One key ingredient that underlies cosmological observables is that the field that sources the observed structure is assumed to be initially Gaussian with high precision. Nevertheless, a minimal deviation from Gaussianityis perhaps the most robust theoretical prediction of models that explain the observed Universe; itis necessarily present even in the simplest scenarios. In addition, most inflationary models produce far higher levels of non-Gaussianity. Since non-Gaussianity directly probes the dynamics in the early Universe, a detection would present a monumental discovery in cosmology, providing clues about physics at energy scales as high as the GUT scale.

Constraints on primordial non-Gaussianity (PNG) will shed light on the origin of primordial fluctuations and the physics of the early universe. The intensity mapping technique is a promising probe of structure formation on large scales; at high redshifts, it can provide complementary information to other cosmological probes. We explore the potential of future wide-field [CII] and CO intensity mapping surveys in constraining PNG of the local shape, which induces a distinct, scale-dependent correction to the line bias. We explore the parameter space of CO and [CII] survey designs that can achieve the nominal target of $\sigma(f_{\rm NL}^{\rm loc}) =1$ and further calculate what constraints such surveys can place on PNG of the equilateral, orthogonal, and quasi-single field shapes. We further test the dependence of these constraints on various model assumptions; namely the halo mass function, modeling of line bias, and correlation between mass and [CII]/CO luminosity. We find that the ability of CO and [CII] intensity mapping surveys to constrain $f_{\rm NL}^{\rm loc}$ relies heavily on the spectral (redshift) coverage, requiring at least a coverage of an octave in frequency, to produce significant results, with the optimized surveys generally covering between $z\sim[2-8]$. As this redshift window partially overlaps with 21-cm EoR experiments like the Hydrogen Epoch of Reionization Array (HERA), we additionally explore the prospects for a hypothetical set of surveys to place constraints on $\sigma(f_{\rm NL}^{\rm loc})$ utilizing multi-tracer analysis between the lines of HI, CO, and [CII].

Primordial non-Gaussianity is an invaluable window into the physical processes that gave rise to the cosmological structure. The presence of local shape PNG imprints a distinct scale-dependent correction to the bias of dark matter tracers on large scales, which can be effectively probed via the technique of intensity mapping. Considering an upcoming generation of experiments, we demonstrate that intensity mapping of CO and [CII] emission can improve upon the current best constraints from the {\it Planck} satellite. We show that measurement of the CO intensity power spectrum by a hypothetical next stage of the ground-based COMAP experiment can achieve $\sigma(f_{\rm NL}^{\rm loc}) = 3.4$, and that the proposed CMB satellite mission PIXIE can achieve $\sigma(f_{\rm NL}^{\rm loc}) = 3.9$ via measurement of [CII] intensity power spectrum.

The presence of massive particles with spin during inflation induces distinct signatures on correlation functions of primordial curvature fluctuations. In particular, the bispectrum of primordial perturbations obtains an angular dependence determined by the spin of the particle, which can be used to set constraints on the presence of such particles. If these particles are long-lived on super-Hubble scales, as is the case for example for partially massless particles, their imprint on correlation functions of curvature perturbations would be unsuppressed. In this paper, we make a forecast for how well such angular dependence can be constrained by the upcoming EUCLID spectroscopic survey via the measurement of galaxy bispectrum.

Massive spinning particles, if present during inflation, lead to a distinctive bispectrum of primordial perturbations, the shape and amplitude of which depend on the masses and spins of the extra particles. This signal, in turn, leaves an imprint in the statistical distribution of galaxies; in particular, as a non-vanishing galaxy bispectrum, which can be used to probe the masses and spins of these particles. In this paper, we present for the first time a new theoretical template for the bispectrum generated by massive spinning particles, valid for a general triangle configuration. We then proceed to perform a Fisher-matrix forecast to assess the potential of two next-generation spectroscopic galaxy surveys, EUCLID and DESI, to constrain the primordial non-Gaussianity sourced by these extra particles. We model the galaxy bispectrum using tree-level perturbation theory, accounting for redshift-space distortions and the Alcock-Paczynski effect, and forecast constraints on the primordial non-Gaussianity parameters marginalizing over all relevant biases and cosmological parameters. Our results suggest that these surveys would potentially be sensitive to any primordial non-Gaussianity with an amplitude larger than $f_{\rm NL}\approx 1$, for massive particles with spins 2, 3, and 4. Interestingly, if non-Gaussianities are present at that level, these surveys will be able to infer the masses of these spinning particles to within tens of percent. If detected, this would provide a very clear window into the particle content of our Universe during inflation.

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.

The presence of additional particles during inflation leads to non-Gaussianity in late-time correlators of primordial curvature perturbations. The shape and amplitude of this signal depend on the mass and spin of the extra particles. Constraints on this distinct form of primordial non-Gaussianity, therefore, provide a wealth of information on the particle content during inflation. We investigate the potential of upcoming galaxy surveys in constraining such a signature through its impact on the observed galaxy power spectrum. Primordial non-Gaussianity of various shapes induces a scale-dependent bias on tracers of large-scale structure, such as galaxies. Using this signature we obtain constraints on massive particles during inflation, which can have non-zero spins. In particular, we show that the prospects for constraining particles with spins 0 and 1 are promising, while constraining particles with spin 2 from power spectrum alone seems challenging. We show that the multi-tracer technique can significantly improve the constraints from the power spectrum by at least an order of magnitude. Furthermore, we analyze the effect of non-linearities due to gravitational evolution on the forecasted constraints on the masses of the extra particles and the amplitudes of the imprinted non-Gaussian signal. We find that gravitational evolution affects the constraints by less than a factor of 2.

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 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 study the impact of the often neglected lensing contribution to galaxy number counts on the $E_g$ statistics which is used to constrain deviations from GR. This contribution affects both the galaxy-galaxy and the convergence-galaxy spectra, while it is larger for the latter. At higher redshifts probed by upcoming surveys, neglecting this term induces a substantial error in the spectra and therefore on the $E_g$ statistics which is constructed from the combination of the two. Moreover, including it, renders the $E_g$ statistics scale and bias-dependent and hence puts into question its very objective.

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 develop a count-in-cells approach to the distribution of ultraviolet background fluctuations that includes source clustering. We demonstrate that an exact expression can be obtained if the clustering of ionizing sources follows the hierarchical ansatz. In this case, the intensity distribution depends solely on their 2-point correlation function. We show that the void scaling function of high redshift mock quasars is consistent with the Negative Binomial form, before applying our formalism to the description of HeII-ionizing fluctuations at the end of helium reionization. The model inputs are the observed quasar luminosity function and 2-point correlation at redshift $z\sim 3$. We find that, for an (comoving) attenuation length $\lesssim $ 55 Mpc, quasar clustering contributes less than 30% of the variance of intensity fluctuations so long as the quasar correlation length does not exceed 15 Mpc. We investigate also the dependence of the intensity distribution on the large-scale environment. Differences in the mean HeII-ionizing intensity between low- and high-density regions could be a factor of few if the sources are highly clustered. An accurate description of quasar demographics and their correlation with strong absorption systems is required to make more precise predictions.

We point out that the ability of some models of inflation, such as Higgs inflation and the universal attractor models, in reproducing the available data is due to their relation to the Starobinsky model of inflation. For large field values, where the inflationary phase takes place, all these classes of models are indeed identical to the Starobinsky model. Nevertheless, the inflaton is just an auxiliary field in the Jordan frame of the Starobinsky model and this leads to two important consequences: first, the inflationary predictions of the Starobinsky model and its descendants are slightly different (albeit not measurably); secondly the theories have different small-field behaviour, leading to different ultra-violet cut-off scales. In particular, one interesting descendant of the Starobinsky model is the non-minimally- coupled quadratic chaotic inflation. Although the standard quadratic chaotic inflation is ruled out by the recent Planck data, its non-minimally coupled version is in agreement with observational data and valid up to Planckian scales.

We study the impact of primordial non-Gaussianity on the density profile of dark matter halos by using the semi-analytical model introduced recently by Dalal {\it et al.} which relates the peaks of the initial linear density field to the final density profile of dark matter halos. Models with primordial non-Gaussianity typically produce an initial density field that differs from that produced in Gaussian models. We use the path-integral formulation of excursion set theory to calculate the non-Gaussian corrections to the peak profile and derive the statistics of the peaks of non-Gaussian density field. In the context of the semi-analytic model for halo profiles, currently allowed values for primordial non-Gaussianity would increase the shapes of the inner dark matter profiles, but only at the sub-percent level except in the very innermost regions.

We study how uncertainty in the reionization history of the universe affects estimates of other cosmological parameters from the Cosmic Microwave Background. We analyze WMAP7 data and synthetic Planck-quality data generated using a realistic scenario for the reionization history of the universe obtained from high-resolution numerical simulation. We perform parameter estimation using a simple sudden reionization approximation, and using the Principal Component Analysis (PCA) technique proposed by Mortonson and Hu. We reach two main conclusions: (1) Adopting a simple sudden reionization model does not introduce measurable bias into values for other parameters, indicating that detailed modeling of reionization is not necessary for the purpose of parameter estimation from future CMB data sets such as Planck. (2) PCA analysis does not allow accurate reconstruction of the actual reionization history of the universe in a realistic case.

We consider Cosmic Microwave Background constraints on inflation models for which the primordial power spectrum is a mixture of perturbations generated by inflaton fluctuations and fluctuations in a curvaton field. If future experiments do not detect isocurvature modes or large non-Gaussianity, it will not be possible to directly distinguish inflaton and curvaton contributions. We investigate whether current and future data can instead constrain the relative contributions of the two sources. We model the spectrum with a bimodal form consisting of a sum of two independent power laws, with different spectral indices. We quantify the ability of current and upcoming data sets to constrain the difference $\Delta n$ in spectral indices, and relative fraction $f$ of the subdominant power spectrum at a pivot scale of $k = 0.017\ {{\rm Mpc^{-1}}\ h}$. Data sets selected are the WMAP 7-year data, alone and in conjunction with South Pole Telescope data, and a synthetic data set comparable to the upcoming Planck data set. We find that current data show no increase in quality of fit for a mixed inflaton/curvaton power spectrum, and a pure power-law spectrum is favored. The ability to constrain independent parameters such as the tensor/scalar ratio is not substantially affected by the additional parameters in the fit. Planck will be capable of placing significant constraints on the parameter space for a bimodal spectrum.

We generalize the cosmological duality between inflation and cyclic contraction under the interchange $a \leftrightarrow H$ to the case of non-canonical scalar field theories with varying speed of sound. The single duality in the canonical case generalizes to a family of three dualities constructed to leave the cosmological acoustic horizon invariant. We find three classes of models: (I) DBI inflation, (II) the non-canonical generalization of cyclic contraction, and (III) a new cosmological solution with rapidly decreasing speed of sound and relatively slowly growing scale factor, which we dub {\it stalled} cosmology. We construct dual analogs to the inflationary slow roll approximation, and solve for the curvature perturbation in all three cases. Both cyclic contraction and stalled cosmology predict a strongly blue spectrum for the curvature perturbations inconsistent with observations.

We investigate the general properties of expanding cosmological models which generate scale-invariant curvature perturbations in the presence of a variable speed of sound. We show that in an expanding universe, generation of a super-Hubble, nearly scale-invariant spectrum of perturbations over a range of wavelengths consistent with observation requires at least one of three conditions: (1) accelerating expansion, (2) a speed of sound faster than the speed of light, or (3) super-Planckian energy density.

In this paper, we use a known duality between expanding and contracting cosmologies to construct a dual of the inflationary flow hierarchy applicable to contracting cosmologies such as Ekpyrotic and Cyclic models. We show that the inflationary flow equations are invariant under the duality and therefore apply equally well to inflation or to cyclic cosmology. We construct a self-consistent small-parameter approximation dual to the slow-roll approximation in inflation, and calculate the power spectrum of perturbations in this limit. We also recover the matter-dominated contracting solution of Wands, and the recently proposed Adiabatic Ekpyrosis solution.