Sharp features in the primordial power spectrum are a powerful window into the inflationary epoch. To date, the cosmic microwave background (CMB) has offered the most sensitive avenue to search for these signatures. In this paper, we demonstrate the power of large-scale structure observations to surpass the CMB as a probe of primordial features. We show that the signatures in galaxy surveys can be separated from the broadband power spectrum and are as robust to the nonlinear evolution of matter as the standard baryon acoustic oscillations. As a result, analyses can exploit a significant range of scales beyond the linear regime available in the datasets. We develop a feature search for large-scale structure, apply it to the final data release of the Baryon Oscillation Spectroscopic Survey and find new bounds on oscillatory features that exceed the sensitivity of Planck for a significant range of frequencies. Moreover, we forecast that the next generation of galaxy surveys, such as DESI and Euclid, will be able to improve current constraints by up to an order of magnitude over an expanded frequency range.

Cosmic voids are biased tracers of the large-scale structure of the universe. Separate universe simulations (SUS) enable accurate measurements of this biasing relation by implementing the peak-background split (PBS). In this work, we apply the SUS technique to measure the void bias parameters. We confirm that the PBS argument works well for underdense tracers. The response of the void size distribution depends on the void radius. For voids larger (smaller) than the size at the peak of the distribution, the void abundance responds negatively (positively) to a long wavelength mode. The linear bias from the SUS is in good agreement with the cross power spectrum measurement on large scales. Using the SUS, we have detected the quadratic void bias for the first time in simulations. We find that $ b_2 $ is negative when the magnitude of $ b_1 $ is small, and that it becomes positive and increases rapidly when $ |b_1| $ increases. We compare the results from voids identified in the halo density field with those from the dark matter distribution, and find that the results are qualitatively similar, but the biases generally shift to the larger voids sizes.

Galaxy shapes have been observed to align with external tidal fields generated by the large-scale structures of the Universe. While the main source for these tidal fields is provided by long-wavelength density perturbations, tensor perturbations also contribute with a non-vanishing amplitude at linear order. We show that parity-breaking gravitational waves produced during inflation leave a distinctive imprint in the galaxy shape power spectrum which is not hampered by any scalar-induced tidal field. We also show that a certain class of tensor non-Gaussianities produced during inflation can leave a signature in the density-weighted galaxy shape power spectrum. We estimate the possibility of observing such imprints in future galaxy surveys.

The understanding of the primordial mechanism that seeded the cosmic structures we observe today in the sky is one of the major goals in cosmology. The leading paradigm for such a mechanism is provided by the inflationary scenario, a period of violent accelerated expansion in the very early stages of evolution of the Universe. While our current knowledge of the physics of inflation is limited to phenomenological models which fit observations, an exquisite understanding of the particle content and interactions taking place during inflation would provide breakthroughs in our understanding of fundamental physics at high energies. In this review, we summarize recent theoretical progress in the modelling of the imprint of primordial interactions in the large scale structures of the Universe. We focus specifically on the effects of such interactions on the statistical distribution of dark matter halos, providing a consistent treatment of the steps required to connect the correlations generated among fields during inflation all the way to the late-time correlations of halos.

We study the large-scale bias parameter of cosmic voids with primordial non-Gaussian (PNG) initial conditions of the local type. In this scenario, the dark matter halo bias exhibits a characteristic scale dependence on large scales, which has been recognized as one of the most promising probes of the local PNG. Using a suite of $N$-body simulations with Gaussian and non-Gaussian initial conditions, we find that the void bias features scale-dependent corrections on large scales, similar to its halo counterpart. We find excellent agreement between the numerical measurement of the PNG void bias and the general peak-background split prediction. Contrary to halos, large voids anti-correlate with the dark matter density field, and the large-scale Gaussian void bias ranges from positive to negative values depending on void size and redshift. Thus, the information in the clustering of voids can be complementary to that of the halos. Using the Fisher matrix formalism for multiple tracers, we demonstrate that including the scale-dependent bias information from voids, constraints on the PNG parameter $f_{\rm NL}$ can be tightened by a factor of two compared to the accessible information from halos alone, when the sampling density of tracers reaches $4 \times 10^{-3} \, h^3 \mathrm{Mpc}^{-3} $.

The expansion of the Universe is understood to have accelerated during two epochs: in its very first moments during a period of Inflation and much more recently, at $z < 1$, when Dark Energy is hypothesized to drive cosmic acceleration. The undiscovered mechanisms behind these two epochs represent some of the most important open problems in fundamental physics. The large cosmological volume at $2 < z < 5$, together with the ability to efficiently target high-$z$ galaxies with known techniques, enables large gains in the study of Inflation and Dark Energy. A future spectroscopic survey can test the Gaussianity of the initial conditions up to a factor of ~50 better than our current bounds, crossing the crucial theoretical threshold of $\sigma(f_{NL}^{\rm local})$ of order unity that separates single field and multi-field models. Simultaneously, it can measure the fraction of Dark Energy at the percent level up to $z = 5$, thus serving as an unprecedented test of the standard model and opening up a tremendous discovery space.

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.

Primordial black holes as dark matter may be generated in single-field models of inflation thanks to the enhancement at small scales of the comoving curvature perturbation. This mechanism requires leaving the slow-roll phase to enter a non-attractor phase during which the inflaton travels across a plateau and its velocity drops down exponentially. We argue that quantum diffusion has a significant impact on the primordial black hole mass fraction making the classical standard prediction not trustable.

Using cosmological N-body simulations, we study the abundance of local maxima (peaks) and minima (dips) identified in the smoothed distribution of halos and dark matter (DM) on scales of $10-100$s Mpcs. The simulations include Gaussian and local-type $f_{\rm NL}$ non-Gaussian initial conditions. The expression derived in the literature for the abundance (irrespective of height) of peaks for Gaussian fields is surprisingly accurate for the evolved halo and DM density fields for all initial conditions considered. Furthermore, the height distribution is very well fitted by a log-normal on quasi-linear scales. The abundance as a function of scale depends on the cosmological parameters ($H_0$ and background matter densities) through the shape of the power spectrum, but it is insensitive to the clustering amplitude. Further, the abundance in the smoothed halo distribution is substantially different in the non-Gaussian from the Gaussian simulations. The interpretation of this effect is straightforward in terms of the scale dependence of halo bias in non-Gaussian models. The abundance of extrema extracted from three-dimensional large galaxy redshift surveys could be a competitive probe of the cosmological parameters and initial non-Gaussianity. It breaks the degeneracy between $f_{\rm NL}$ and the clustering amplitude, making it complementary to counts of galaxy clusters and peaks in weak-lensing maps.

We study the imprints of a massive spin-2 field on inflationary observables, and in particular on the breaking of consistency relations. In this setup, the minimal inflationary field content interacts with the massive spin-2 field through dRGT interactions, thus guaranteeing the absence of Boulware-Deser ghostly degrees of freedom. The unitarity requirement on spinning particles, known as Higuchi bound, plays a crucial role for the size of the observable signal.

We measure the large-scale bias of dark matter halos in simulations with non-Gaussian initial conditions of the local type, and compare this bias to the response of the mass function to a change in the primordial amplitude of fluctuations. The two are found to be consistent, as expected from physical arguments, for three halo-finder algorithms which use different Spherical Overdensity (SO) and Friends-of-Friends (FoF) methods. On the other hand, we find that the commonly used prediction for universal mass functions, that the scale-dependent bias is proportional to the first-order Gaussian Lagrangian bias, does not yield a good agreement with the measurements. For all halo finders, high-mass halos show a non-Gaussian bias suppressed by 10-15% relative to the universal mass function prediction. For SO halos, this deviation changes sign at low masses, where the non-Gaussian bias becomes larger than the universal prediction.

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.

Dark matter halos are the building blocks of the universe as they host galaxies and clusters. The knowledge of the clustering properties of halos is therefore essential for the understanding of the galaxy statistical properties. We derive an effective halo Boltzmann equation which can be used to describe the halo clustering statistics. In particular, we show how the halo Boltzmann equation encodes a statistically biased gravitational force which generates a bias in the peculiar velocities of virialized halos with respect to the underlying dark matter, as recently observed in N-body simulations.

In this paper we operate under the assumption that no tensors from inflation will be measured in the future by the dedicated experiments and argue that, while for single-field slow-roll models of inflation the running of the spectral index will be hard to be detected, in multi-field models the running can be large due to its strong correlation with non-Gaussianity. A detection of the running might therefore be related to the presence of more than one active scalar degree of freedom during inflation.

With the advent of large scale galaxy surveys, constraints on primordial non-Gaussianity (PNG) are expected to reach ${\cal O}(f_\text{NL}) \sim 1$. In order to fully exploit the potential of these future surveys, a deep theoretical understanding of the signatures imprinted by PNG on the large scale structure of the Universe is necessary. In this paper, we explore the effect of a stochastic moving barrier on the amplitude of the non-Gaussian bias induced by local quadratic PNG. We show that, in the peak approach to halo clustering, the amplitude of the non-Gaussian bias will generally differ from the peak-background split prediction unless the barrier is flat and deterministic. For excursion set peaks with a square-root barrier, which reproduce reasonably well the linear bias $b_1$ and mass function $\bar{n}_\text{h}$ of SO haloes, the non-Gaussian bias amplitude is $\sim 40$% larger than the peak-background split expectation $d\ln\bar{n}_\text{h}/d\ln\sigma_8$ for haloes of mass $\sim 10^{13} {\it h}^{-1}M_\odot$ at $z=0$. Furthermore, we argue that the effect of PNG on squeezed configurations of the halo bispectrum differs significantly from that predicted by standard local bias approaches. Our predictions can be easily confirmed, or invalidated, with N-body simulations.

It has been recently shown that any halo velocity bias present in the initial conditions does not decay to unity, in agreement with predictions from peak theory. However, this is at odds with the standard formalism based on the coupled fluids approximation for the coevolution of dark matter and halos. Starting from conservation laws in phase space, we discuss why the fluid momentum conservation equation for the biased tracers needs to be modified in accordance with the change advocated in Baldauf, Desjacques & Seljak (2014). Our findings indicate that a correct description of the halo properties should properly take into account peak constraints when starting from the Vlasov-Boltzmann equation.

The most conventional mechanism for gravitational waves (gw) production during inflation is the amplification of vacuum metric fluctuations. In this case the gw production can be uniquely related to the inflationary expansion rate $H$. For example, a gw detection close to the present experimental limit (tensor-to-scalar ratio $r \sim 0.1$) would indicate an inflationary expansion rate close to $10^{14} \, {\rm GeV}$. This conclusion, however, would be invalid if the observed gw originated from a different source. We construct and study one of the possible covariant formulations of the mechanism suggested in [43], where a spectator field $\sigma$ with a sound speed $c_{s} \ll 1$ acts as a source for gw during inflation. In our formulation $\sigma$ is described by a so-called $P(X)$ Lagrangian and a non-minimal coupling to gravity. This field interacts only gravitationally with the inflaton, which has a standard action. We compute the amount of scalar and tensor density fluctuations produced by $\sigma$ and find that, in our realization, $r$ is not enhanced with respect to the standard result but it is strongly sensitive to $c_s$, thus breaking the direct $r \leftrightarrow H$ connection.

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.

Understanding the biasing between the clustering properties of halos and the underlying dark matter distribution is important for extracting cosmological information from ongoing and upcoming galaxy surveys. While on sufficiently larges scales the halo overdensity is a local function of the mass density fluctuations, on smaller scales the gravitational evolution generates non-local terms in the halo density field. We characterize the magnitude of these contributions at third-order in perturbation theory by identifying the coefficients of the non-local invariant operators, and extend our calculation to include non-local (Lagrangian) terms induced by a peak constraint. We apply our results to describe the scale-dependence of halo bias in cosmologies with massive neutrinos. The inclusion of gravity-induced non-local terms and, especially, a Lagrangian $k^2$-contribution is essential to reproduce the numerical data accurately. We use the peak-background split to derive the numerical values of the various bias coefficients from the excursion set peak mass function. For neutrino masses in the range $0\leq \sum_i m_{\nu_i} \leq 0.6$ eV, we are able to fit the data with a precision of a few percents up to $k=0.3\, h {\rm \,Mpc^{-1}}$ without any free parameter.

We investigate nonlocal Lagrangian bias contributions involving gradients of the linear density field, for which we have predictions from the excursion set peak formalism. We begin by writing down a bias expansion which includes all the bias terms, including the nonlocal ones. Having checked that the model furnishes a reasonable fit to the halo mass function, we develop a 1-point cross-correlation technique to measure bias factors associated with 2-distributed quantities. We validate the method with numerical realizations of peaks of Gaussian random fields before we apply it to N-body simulations. We focus on the lowest (quadratic) order nonlocal contributions. We can reproduce our measurement of \chi_{10} if we allow for an offset between the Lagrangian halo center-of-mass and the peak position. The sign and magnitude of \chi_{10} is consistent with Lagrangian haloes sitting near linear density maxima. The resulting contribution to the halo bias can safely be ignored for M = 10^13 Msun/h, but could become relevant at larger halo masses. For the second nonlocal bias \chi_{01} however, we measure a much larger magnitude than predicted by our model. We speculate that some of this discrepancy might originate from nonlocal Lagrangian contributions induced by nonspherical collapse.

We show that a large contribution to tensor modes during inflation can be generated by a spectator scalar field with speed of sound lower than unity.

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.

Multi-field models of inflation predict an inequality between the amplitude tauNL of the collapsed limit of the four-point correlator of the primordial curvature perturbation and the amplitude fNL of the squeezed limit of its three-point correlator. While a convincing detection of non-Gaussianity through the squeezed limit of the three-point correlator would rule out all single-field models, a robust confirmation or disproval of the inequality between tauNL and fNL would provide crucial information about the validity of multi-field models of inflation. In this paper, we discuss to which extent future measurements of the scale-dependence of galaxy bias can test multi-field inflationary scenarios. The strong degeneracy between the effect of a non-vanishing fNL and tauNL on halo bias can be broken by considering multiple tracer populations of the same surveyed volume. If halos down to 1e13 Msun/h are resolved in a survey of volume 25(Gpc/h)^3, then testing multi-field models of inflation at the 3-\sigma level would require, for instance, a detection of tauNL at the level of tauNL~1e5 given a measurement of a local bispectrum with amplitude fNL~10. However, we find that disproving multi-field models of inflation with measurements of the non-Gaussian bias only will be very challenging, unless |fNL| > 80 and one can achieve a halo mass resolution of 1e10 Msun/h.