The forthcoming generation of galaxy redshift surveys will sample the large-scale structure of the Universe over unprecedented volumes with high-density tracers. This advancement will make robust measurements of three-point clustering statistics possible. In preparation for this improvement, we investigate how several methodological choices can influence inferences based on the bispectrum about galaxy bias and shot noise. We first measure the real-space bispectrum of dark-matter haloes extracted from 298 N-body simulations covering a volume of approximately $1000 h^{-3} \mathrm{Gpc}^3$. We then fit a series of theoretical models based on tree-level perturbation theory to the numerical data. To achieve this, we estimate the covariance matrix of the measurement errors by using 10,000 mock catalogues generated with the Pinocchio code. We study how the model constraints are influenced by the binning strategy for the bispectrum configurations and by the form of the likelihood function. We also use Bayesian model-selection techniques to single out the optimal theoretical description of our data. We find that a three-parameter bias model combined with Poissonian shot noise is necessary to model the halo bispectrum up to scales of $k_\mathrm{max}\lesssim 0.08 h \mathrm{Mpc}^{-1}$, although fitting formulae that relate the bias parameters can be helpful to reduce the freedom of the model without compromising accuracy. Our data clearly disfavour local Eulerian and local Lagrangian bias models and do not require corrections to Poissonian shot noise. We anticipate that model-selection diagnostics will be particularly useful to extend the analysis to smaller scales as, in this case, the number of model parameters will grow significantly large.

This paper is the first in a set that analyses the covariance matrices of clustering statistics obtained from several approximate methods for gravitational structure formation. We focus here on the covariance matrices of anisotropic two-point correlation function measurements. Our comparison includes seven approximate methods, which can be divided into three categories: predictive methods that follow the evolution of the linear density field deterministically (ICE-COLA, Peak Patch, and Pinocchio), methods that require a calibration with N-body simulations (Patchy and Halogen), and simpler recipes based on assumptions regarding the shape of the probability distribution function (PDF) of density fluctuations (log-normal and Gaussian density fields). We analyse the impact of using covariance estimates obtained from these approximate methods on cosmological analyses of galaxy clustering measurements, using as a reference the covariances inferred from a set of full N-body simulations. We find that all approximate methods can accurately recover the mean parameter values inferred using the N-body covariances. The obtained parameter uncertainties typically agree with the corresponding N-body results within 5% for our lower mass threshold, and 10% for our higher mass threshold. Furthermore, we find that the constraints for some methods can differ by up to 20% depending on whether the halo samples used to define the covariance matrices are defined by matching the mass, number density, or clustering amplitude of the parent N-body samples. The results of our configuration-space analysis indicate that most approximate methods provide similar results, with no single method clearly outperforming the others.

Massive neutrinos, due to their free streaming, produce a suppression in the matter power spectrum at intermediate and small scales which could be probed by galaxy clustering and/or weak lensing observables. This effect happens at scales that are also influenced by baryon feedback, i.e. galactic winds or Active Galactic Nuclei (AGN) feedback, which in realistic hydrodynamic simulations has also been shown to produce a suppression of power. Leaving aside, for the moment, the complex issue of galaxy bias, we focus here on matter clustering and tomographic weak lensing, we investigate the possible degeneracy between baryon feedback and neutrinos showing that it is not likely to degrade significantly the measurement of neutrino mass in future surveys. To do so, we generate mock data sets and fit them using the Markov Chain Monte Carlo (MCMC) technique and explore degeneracies between feedback parameters and neutrino mass. We model baryon feedback through fitting functions, while massive neutrinos are accounted for, also in the non-linear regime, using Halofit calibrated against accurate N-body neutrino simulations. In the error budget, we include the uncertainty in the modelling of non-linearities. For both matter clustering and weak lensing, we always recover the input neutrino mass within $\sim 0.25\sigma$ confidence level. Finally, we also take into account the intrinsic alignment effect in the weak lensing mock data. Even in this case, we are able to recover the right parameters: in particular, we find a significant degeneracy pattern between $M_\nu$ and the intrinsic alignment parameter $A_\mathrm{IA}$ .

We present a blind method to determine the properties of a foreground contamination, given by a visibility mask, that affects a deep galaxy survey. Angular cross correlations of density fields in different redshift bins are expected to vanish (apart from a contribution due to lensing), but are sensitive to the presence of a foreground that modulates the flux limit across the sky. After formalizing the expected effect of a foreground mask on the measured galaxy density, under a linear, luminosity-dependent bias model for galaxies, we construct two estimators that single out the mask contribution if a sufficient number of independent redshift bins is available. These estimators are combined to give a reconstruction of the mask. We use Milky-Way reddening as a prototype for the mask. Using a set of 20 large mock catalogs covering $1/4$-th of the sky and number-matched to $H\alpha$ emitters to mimic an Euclid-like sample, we demonstrate that our method can reconstruct the mask and its angular clustering at scales $\ell<100$, beyond which the cosmological signal becomes dominant. The uncertainty of this reconstruction is quantified to be $1/3$-rd of the sample variance of the signal. Such a reconstruction requires knowledge of the average and square average of the mask, but we show that it is possible to recover this information either from external models or internally from the data. It also relies on knowledge of how the impact of the foreground changes with redshift (due to the extinction curve in our case), but this can be tightly constrained by cross correlations of different redshift bins. The strong points of this blind reconstruction technique lies in the ability to find "unknown unknowns" that affect a survey, and in the facility to quantify, using sets of mock catalogs, how its uncertainty propagates to clustering measurements. [Abridged]

We study the accuracy of several approximate methods for gravitational dynamics in terms of halo power spectrum multipoles and their estimated covariance matrix. We propagate the differences in covariances into parameter constrains related to growth rate of structure, Alcock-Paczynski distortions and biasing. We consider seven methods in three broad categories: algorithms that solve for halo density evolution deterministically using Lagrangian trajectories (ICE-COLA, Pinocchio and PeakPatch), methods that rely on halo assignment schemes onto dark-matter overdensities calibrated with a target N-body run (Halogen, Patchy) and two standard assumptions about the full density PDF (Gaussian and Lognormal). We benchmark their performance against a set of three hundred N-body simulations, running similar sets of approximate simulations with matched initial conditions, for each method. We find that most methods reproduce the monopole to within $5\%$, while residuals for the quadrupole are sometimes larger and scale dependent. The variance of the multipoles is typically reproduced within $10\%$. Overall, we find that covariances built from approximate simulations yield errors on model parameters within $10\%$ of those from the N-body based covariance.

We estimate the velocity field in a large set of $N$-body simulations including massive neutrino particles, and measure the auto-power spectrum of the velocity divergence field as well as the cross-power spectrum between the cold dark matter density and the velocity divergence. We perform these measurements at four different redshifts and within four different cosmological scenarios, covering a wide range in neutrino masses. We find that the nonlinear correction to the velocity power spectra largely depend on the degree of nonlinear evolution with no specific dependence on the value of neutrino mass. We provide a fitting formula, based on the value of the r.m.s. of the matter fluctuations in spheres of $8h^{-1}$Mpc, describing the nonlinear corrections with 3\% accuracy on scales below $k=0.7\; h$ Mpc$^{-1}$.

We compare the measurements of the bispectrum and the estimate of its covariance obtained from a set of different methods for the efficient generation of approximate dark matter halo catalogs to the same quantities obtained from full N-body simulations. To this purpose we employ a large set of three-hundred realisations of the same cosmology for each method, run with matching initial conditions in order to reduce the contribution of cosmic variance to the comparison. In addition, we compare how the error on cosmological parameters such as linear and nonlinear bias parameters depends on the approximate method used for the determination of the bispectrum variance. As general result, most methods provide errors within 10% of the errors estimated from N-body simulations. Exceptions are those methods requiring calibration of the clustering amplitude but restrict this to two-point statistics. Finally we test how our results are affected by being limited to a few hundreds measurements from N-body simulation, and therefore to the bispectrum variance, by comparing with a larger set of several thousands realisations performed with one approximate method.

The main effect of massive neutrinos on the large-scale structure consists in a few percent suppression of matter perturbations on all scales below their free-streaming scale. Such effect is of particular importance as it allows to constraint the value of the sum of neutrino masses from measurements of the galaxy power spectrum. In this work, we present the first measurements of the next higher-order correlation function, the bispectrum, from N-body simulations that include massive neutrinos as particles. This is the simplest statistics characterising the non-Gaussian properties of the matter and dark matter halos distributions. We investigate, in the first place, the suppression due to massive neutrinos on the matter bispectrum, comparing our measurements with the simplest perturbation theory predictions, finding the approximation of neutrinos contributing at quadratic order in perturbation theory to provide a good fit to the measurements in the simulations. On the other hand, as expected, a linear approximation for neutrino perturbations would lead to O($f_{\nu}$) errors on the total matter bispectrum at large scales. We then attempt an extension of previous results on the universality of linear halo bias in neutrino cosmologies, to non-linear and non-local corrections finding consistent results with the power spectrum analysis.

Efficient estimators of Fourier-space statistics for large number of objects rely on Fast Fourier Transforms (FFTs), which are affected by aliasing from unresolved small scale modes due to the finite FFT grid. Aliasing takes the form of a sum over images, each of them corresponding to the Fourier content displaced by increasing multiples of the sampling frequency of the grid. These spurious contributions limit the accuracy in the estimation of Fourier-space statistics, and are typically ameliorated by simultaneously increasing grid size and discarding high-frequency modes. This results in inefficient estimates for e.g. the power spectrum when desired systematic biases are well under per-cent level. We show that using interlaced grids removes odd images, which include the dominant contribution to aliasing. In addition, we discuss the choice of interpolation kernel used to define density perturbations on the FFT grid and demonstrate that using higher-order interpolation kernels than the standard Cloud in Cell algorithm results in significant reduction of the remaining images. We show that combining fourth-order interpolation with interlacing gives very accurate Fourier amplitudes and phases of density perturbations. This results in power spectrum and bispectrum estimates that have systematic biases below 0.01% all the way to the Nyquist frequency of the grid, thus maximizing the use of unbiased Fourier coefficients for a given grid size and greatly reducing systematics for applications to large cosmological datasets.

We present a test to quantify how well some approximate methods, designed to reproduce the mildly non-linear evolution of perturbations, are able to reproduce the clustering of DM halos once the grouping of particles into halos is defined and kept fixed. The following methods have been considered: Lagrangian Perturbation Theory (LPT) up to third order, Truncated LPT, Augmented LPT, MUSCLE and COLA. The test runs as follows: halos are defined by applying a friends-of-friends (FoF) halo finder to the output of an N-body simulation. The approximate methods are then applied to the same initial conditions of the simulation, producing for all particles displacements from their starting position and velocities. The position and velocity of each halo are computed by averaging over the particles that belong to that halo, according to the FoF halo finder. This procedure allows us to perform a well-posed test of how clustering of the matter density and halo density fields are recovered, without asking to the approximate method an accurate reconstruction of halos. We have considered the results at $z=0,0.5,1$, and we have analysed power spectrum in real and redshift space, object-by-object difference in position and velocity, density Probability Distribution Function (PDF) and its moments, phase difference of Fourier modes. We find that higher LPT orders are generally able to better reproduce the clustering of halos, while little or no improvement is found for the matter density field when going to 2LPT and 3LPT. Augmentation provides some improvement when coupled with 2LPT, while its effect is limited when coupled with 3LPT. Little improvement is brought by MUSCLE with respect to Augmentation. The more expensive particle-mesh code COLA outperforms all LPT methods [abridged]

The forecasted accuracy of upcoming surveys of large-scale structure cannot be achieved without a proper quantification of the error induced by foreground removal (or other systematics like 0-point photometry offset). Because these errors are highly correlated on the sky, their influence is expected to be especially important at very large scales. In this work we quantify how the uncertainty in the visibility mask of a survey influences the measured power spectrum of a sample of tracers of the density field and its covariance matrix. We start from a very large set of 10,000 catalogs of dark matter halos in periodic cosmological boxes, produced with the PINOCCHIO approximate method. To make an analytic approach feasible, we assume luminosity-independent halo bias and an idealized geometry for the visibility mask. We find that the power spectrum of these biased tracers can be expressed as the sum of a cosmological term, a mask term and a term involving their convolution. The mask and convolution terms scale like $P\propto l^2\sigma_A^2$, where $\sigma_A^2$ is the variance of the uncertainty on the visibility mask. With $l=30-100$ Mpc$/h$ and $\sigma_A=5-20$\%, the mask term can be significant at $k\sim0.01-0.1\ h/$Mpc, and the convolution term can amount to $\sim 1-10$\% of the total. For the power spectrum covariance, the coupling of the convolution term with the other two gives rise to several mixed terms, that we quantify by difference using the mock catalogs. These are found to be of the same order of the mask covariance, and to introduce non-diagonal terms at large scales. Then, the power spectrum covariance matrix cannot be expressed as the sum of a cosmological and of a mask term. Our results lie down the theoretical bases to quantify the impact that uncertainties in the mask calibration have on the derivation of cosmological constraints from large spectroscopic surveys. [Abridged]

The set-up of the initial conditions in cosmological N-body simulations is usually implemented by rescaling the desired low-redshift linear power spectrum to the required starting redshift consistently with the Newtonian evolution of the simulation. The implementation of this practical solution requires more care in the context of massive neutrino cosmologies, mainly because of the non-trivial scale-dependence of the linear growth that characterises these models. In this work we consider a simple two-fluid, Newtonian approximation for cold dark matter and massive neutrinos perturbations that can reproduce the cold matter linear evolution predicted by Boltzmann codes such as CAMB or CLASS with a 0.1% accuracy or below for all redshift relevant to nonlinear structure formation. We use this description, in the first place, to quantify the systematic errors induced by several approximations often assumed in numerical simulations, including the typical set-up of the initial conditions for massive neutrino cosmologies adopted in previous works. We then take advantage of the flexibility of this approach to rescale the late-time linear power spectra to the simulation initial redshift, in order to be as consistent as possible with the dynamics of the N-body code and the approximations it assumes. We implement our method in a public code providing the initial displacements and velocities for cold dark matter and neutrino particles that will allow accurate, i.e. one-percent level, numerical simulations for this cosmological scenario.

We present a method that extends the capabilities of the PINpointing Orbit-Crossing Collapsed HIerarchical Objects (PINOCCHIO) code, allowing it to generate accurate dark matter halo mock catalogues in cosmological models where the linear growth factor and the growth rate depend on scale. Such cosmologies comprise, among others, models with massive neutrinos and some classes of modified gravity theories. We validate the code by comparing the halo properties from PINOCCHIO against N-body simulations, focusing on cosmologies with massive neutrinos: $\nu\Lambda$CDM. We analyse the halo mass function, halo two-point correlation function, halo power spectrum and the moments of the halo density field, showing that PINOCCHIO reproduces the results from simulations with the same level of precision as the original code ($\sim5-10\%$). We demonstrate that the abundance of halos in cosmologies with massless and massive neutrinos from PINOCCHIO matches very well the outcome of simulations, and point out that PINOCCHIO can reproduce the $\Omega_\nu-\sigma_8$ degeneracy that affects the halo mass function. We show that the clustering properties of the halos from PINOCCHIO matches accurately those from simulations both in real and redshift-space, in the latter case up to $k=0.3~h~{\rm Mpc}^{-1}$. We finally point out that the first moments of the halo density field from simulations are precisely reproduced by PINOCCHIO. We emphasize that the computational time required by PINOCCHIO to generate mock halo catalogues is orders of magnitude lower than the one needed for N-body simulations. This makes this tool ideal for applications like covariance matrix studies within the standard $\Lambda$CDM model but also in cosmologies with massive neutrinos or some modified gravity theories.

We present the latest version of Pinocchio, a code that generates catalogues of DM haloes in an approximate but fast way with respect to an N-body simulation. This code version extends the computation of particle and halo displacements up to 3rd-order Lagrangian Perturbation Theory (LPT), in contrast with previous versions that used Zeldovich approximation (ZA). We run Pinocchio on the same initial configuration of a reference N-body simulation, so that the comparison extends to the object-by-object level. We consider haloes at redshifts 0 and 1, using different LPT orders either for halo construction - where displacements are needed to decide particle accretion onto a halo or halo merging - or to compute halo final positions. We compare the clustering properties of Pinocchio haloes with those from the simulation by computing the power spectrum and 2-point correlation function (2PCF) in real and redshift space (monopole and quadrupole), the bispectrum and the phase difference of halo distributions. We find that 2LPT and 3LPT give noticeable improvement. 3LPT provides the best agreement with N-body when it is used to displace haloes, while 2LPT gives better results for constructing haloes. At the highest orders, linear bias is typically recovered at a few per cent level. In Fourier space and using 3LPT for halo displacements, the halo power spectrum is recovered to within 10 per cent up to $k_{max}\sim0.5\ h/$Mpc. The results presented in this paper have interesting implications for the generation of large ensemble of mock surveys aimed at accurately compute covariance matrices for clustering statistics.

(abridged) We analyse the clustering features of Large Scale Structures (LSS) in the presence of massive neutrinos, employing a set of large-volume, high-resolution cosmological N-body simulations, where neutrinos are treated as a separate collisionless fluid. The volume of 8$\cGpc$, combined with a resolution of about $8\times 10^{10}\Ms$ for the cold dark matter (CDM) component, represents a significant improvement over previous N-body simulations in massive neutrino cosmologies. We show that most of the nonlinear evolution is generated exclusively by the CDM component. We find that accounting only for the nonlinear evolution of the CDM power spectrum allows to recover the total matter power spectrum with the same accuracy as the massless case. Indeed, we show that, the most recent version of the \halofit\ formula calibrated on $\Lambda$CDM simulations can be applied directly to the linear CDM power spectrum without requiring additional fitting parameters in the massive case. As a second step, we study the abundance and clustering properties of CDM halos, confirming that, in massive neutrino cosmologies, the proper definition of the halo bias should be made with respect to the {\em cold} rather than the {\em total} matter distribution, as recently shown in the literature. Here we extend these results to the redshift space, finding that, when accounting for massive neutrinos, an improper definition of the linear bias can lead to a systematic error of about 1-$2 \%$ in the determination of the linear growth rate from anisotropic clustering. This result is quite important if we consider that future spectroscopic galaxy surveys, as \eg\ Euclid, are expected to measure the linear growth-rate with statistical errors less than about $3 \%$ at $z\lesssim1$.

We re-evaluate the extragalactic gamma-ray flux prediction from dark matter annihilation in the approach of integrating over the nonlinear matter power spectrum, extrapolated to the free-streaming scale. We provide an estimate of the uncertainty based entirely on available N-body simulation results and minimal theoretical assumptions. We illustrate how an improvement in the simulation resolution, exemplified by the comparison between the Millennium and Millennium II simulations, affects our estimate of the flux uncertainty and we provide a "best guess" value for the flux multiplier, based on the assumption of stable clustering for the dark matter perturbations described as a collision-less fluid. We achieve results comparable to traditional Halo Model calculations, but with a much simpler procedure and a more general approach, as it relies only on one, directly measurable quantity. In addition we discuss the extension of our calculation to include baryonic effects as modeled in hydrodynamical cosmological simulations and other possible sources of uncertainty that would in turn affect indirect dark matter signals. Upper limit on the integrated power spectrum from supernovae lensing magnification are also derived and compared with theoretical expectations.

We extend the resummation method of Anselmi & Pietroni (2012) to compute the total density power spectrum in models of quintessence characterized by a vanishing speed of sound. For standard $\Lambda$CDM cosmologies, this resummation scheme allows predictions with an accuracy at the few percent level beyond the range of scales where acoustic oscillations are present, therefore comparable to other, common numerical tools. In addition, our theoretical approach indicates an approximate but valuable and simple relation between the power spectra for standard quintessence models and models where scalar field perturbations appear at all scales. This, in turn, provides an educated guess for the prediction of nonlinear growth in models with generic speed of sound, particularly valuable since no numerical results are yet available.

We use a suite of N-body simulations that incorporate massive neutrinos as an extra-set of particles to investigate their effect on the halo mass function. We show that for cosmologies with massive neutrinos the mass function of dark matter haloes selected using the spherical overdensity (SO) criterion is well reproduced by the fitting formula of Tinker et al. (2008) once the cold dark matter power spectrum is considered instead of the total matter power, as it is usually done. The differences between the two implementations, i.e. using $P_{\rm cdm}(k)$ instead of $P_{\rm m}(k)$, are more pronounced for large values of the neutrino masses and in the high end of the halo mass function: in particular, the number of massive haloes is higher when $P_{\rm cdm}(k)$ is considered rather than $P_{\rm m}(k)$. As a quantitative application of our findings we consider a Planck-like SZ-clusters survey and show that the differences in predicted number counts can be as large as $30\%$ for $\sum m_\nu = 0.4$ eV. Finally, we use the Planck-SZ clusters sample, with an approximate likelihood calculation, to derive Planck-like constraints on cosmological parameters. We find that, in a massive neutrino cosmology, our correction to the halo mass function produces a shift in the $\sigma_8(\Omega_{\rm m}/0.27)^\gamma$ relation which can be quantified as $\Delta \gamma \sim 0.05$ and $\Delta \gamma \sim 0.14$ assuming one ($N_\nu=1$) or three ($N_\nu=3$) degenerate massive neutrino, respectively. The shift results in a lower mean value of $\sigma_8$ with $\Delta \sigma_8 = 0.01$ for $N_\nu=1$ and $\Delta \sigma_8 = 0.02$ for $N_\nu=3$, respectively. Such difference, in a cosmology with massive neutrinos, would increase the tension between cluster abundance and Planck CMB measurements.

We use a large suite of N-body simulations to study departures from universality in halo abundances and clustering in cosmologies with non-vanishing neutrino masses. To this end, we study how the halo mass function and halo bias factors depend on the scaling variable $\sigma^2(M,z)$, the variance of the initial matter fluctuation field, rather than on halo mass $M$ and redshift $z$ themselves. We show that using the variance of the cold dark matter rather than the total mass field, i.e., $\sigma^2_{cdm}(M,z)$ rather than $\sigma^2_{m}(M,z)$, yields more universal results. Analysis of halo bias yields similar conclusions: When large-scale halo bias is defined with respect to the cold dark matter power spectrum, the result is both more universal, and less scale- or $k$-dependent. These results are used extensively in Papers I and III of this series.

By using a suite of large box-size N-body simulations that incorporate massive neutrinos as an extra set of particles, we investigate the impact of neutrino masses on the spatial distribution of dark matter haloes and galaxies. We compute the bias between the spatial distribution of dark matter haloes and the overall matter and cold dark matter distributions using statistical tools such as the power spectrum and the two-point correlation function. Overall we find a scale-dependent bias on large scales for the cosmologies with massive neutrinos. However, our results indicate that the scale-dependence in the bias is reduced if the latter is computed with respect to the cold dark matter distribution only. We find that the value of the bias on large scales is reasonably well reproduced by the Tinker fitting formula once the linear cold dark matter power spectrum is used, instead of the total matter power spectrum. We investigate whether scale-dependent bias really comes from purely neutrino's effect or from nonlinear gravitational collapse of haloes. For this purpose, we address the $\Omega_\nu$-$\sigma_8$ degeneracy and find that such degeneracy is not perfect, implying that neutrinos imprint a slight scale dependence on the large-scale bias. Finally, by using a simple halo occupation distribution (HOD) model, we investigate the impact of massive neutrinos on the distribution of galaxies within dark matter haloes. We use the main galaxy sample in the Sloan Digital Sky Survey II Data Release 7 to investigate if the small-scale galaxy clustering alone can be used to discriminate among different cosmological models with different neutrino masses. Our results suggest that different choices of the HOD parameters can reproduce the observational measurements relatively well, and we quantify the difference between the values of the HOD parameters between massless and massive neutrino cosmologies.

We explore the scale dependence of halo bias using real space cross-correlation measurements in N-body simulations and in Pinocchio, an algorithm based on Lagrangian Perturbation Theory. Recent work has shown how to interpret such real space measurements in terms of k-dependent bias in Fourier space, and how to remove the k-dependence to reconstruct the k-independent peak-background split halo bias parameters. We compare our reconstruction of the linear bias, which requires no free parameters, with previous estimates from N-body simulations which were obtained directly in Fourier space at large scales, and find very good agreement. Our reconstruction of the quadratic bias is similarly parameter-free, although in this case there are no previous Fourier space measurements to compare with. Our analysis of N-body simulations explicitly tests the predictions of the excursion set peaks (ESP) formalism of Paranjape et al. (2013) for the scale dependence of bias; we find that the ESP predictions accurately describe our measurements. In addition, our measurements in Pinocchio serve as a useful, successful consistency check between Pinocchio and N-body simulations that is not accessible to traditional measurements.

We present a new parallel implementation of the PINpointing Orbit Crossing-Collapsed HIerarchical Objects (PINOCCHIO) algorithm, a quick tool, based on Lagrangian Perturbation Theory, for the hierarchical build-up of Dark Matter halos in cosmological volumes. To assess its ability to predict halo correlations on large scales, we compare its results with those of an N-body simulation of a 3 Gpc/h box sampled with 2048^3 particles taken from the MICE suite, matching the same seeds for the initial conditions. Thanks to the FFTW libraries and to the relatively simple design, the code shows very good scaling properties. The CPU time required by PINOCCHIO is a tiny fraction (~1/2000) of that required by the MICE simulation. Varying some of PINOCCHIO numerical parameters allows one to produce a universal mass function that lies in the range allowed by published fits, although it underestimates the MICE mass function of FoF halos in the high mass tail. We compare the matter-halo and the halo-halo power spectra with those of the MICE simulation and find that these 2-point statistics are well recovered on large scales. In particular, when catalogs are matched in number density, agreement within ten per cent is achieved for the halo power spectrum. At scales k>0.1 h/Mpc, the inaccuracy of the Zel'dovich approximation in locating halo positions causes an underestimate of the power spectrum that can be modeled as a Gaussian factor with a damping scale of d=3 Mpc/h at z=0, decreasing at higher redshift. Finally, a remarkable match is obtained for the reduced halo bispectrum, showing a good description of nonlinear halo bias. Our results demonstrate the potential of PINOCCHIO as an accurate and flexible tool for generating large ensembles of mock galaxy surveys, with interesting applications for the analysis of large galaxy redshift surveys.

We compute the matter bispectrum in the presence of primordial local non-Gaussianity over a wide range of scales, including the very small nonlinear ones. We use the Halo Model approach, considering non-Gaussian corrections to the halo profiles, the halo mass function and the bias functions. We compare our results in the linear and mildly nonlinear scales to a large ensemble of Gaussian and non-Gaussian numerical simulations. We consider both squeezed and equilateral configurations, at redshift z = 0 and z = 1. For z = 0, the deviations between the Halo Model and the simulations are smaller than 10% in the squeezed limit, both in the Gaussian and non-Gaussian cases. The Halo Model allows to make predictions on scales much smaller than those reached by numerical simulations. For local non-Gaussian initial conditions with a parameter fNL = 100, we find an enhancement of the bispectrum in the squeezed configuration k = k3 = k2 >> k1 \sim 0.01 h^{-1} Mpc, of \sim 15% and \sim 25% on scales k \sim 1 h^{-1} Mpc, at z = 0 and z = 1 respectively. This is mainly due to the non-Gaussian corrections in the linear bias. Finally we provide a very simple expression valid for any scenario, i.e. for any choice of the halo profile, mass and bias functions, which allow for a fast evaluation of the bispectrum on squeezed configurations.

Quasi-single field inflation predicts a peculiar momentum dependence in the squeezed limit of the primordial bispectrum which smoothly interpolates between the local and equilateral models. This dependence is directly related to the mass of the isocurvatons in the theory which is determined by the supersymmetry. Therefore, in the event of detection of a non-zero primordial bispectrum, additional constraints on the parameter controlling the momentum-dependence in the squeezed limit becomes an important question. We explore the effects of these non-Gaussian initial conditions on large-scale structure and the cosmic microwave background, with particular attention to the galaxy power spectrum at large scales and scale-dependence corrections to galaxy bias. We determine the simultaneous constraints on the two parameters describing the QSF bispectrum that we can expect from upcoming large-scale structure and cosmic microwave background observations. We find that for relatively large values of the non-Gaussian amplitude parameters, but still well within current uncertainties, galaxy power spectrum measurements will be able to distinguish the QSF scenario from the predictions of the local model. A CMB likelihood analysis, as well as Fisher matrix analysis, shows that there is also a range of parameter values for which Planck data may be able distinguish between QSF models and the related local and equilateral shapes. Given the different observational weightings of the CMB and LSS results, degeneracies can be significantly reduced in a joint analysis.

We revisit the computation of the extragalactic gamma-ray signal from cosmological dark matter annihilations. The prediction of this signal is notoriously model dependent, due to different descriptions of the clumpiness of the dark matter distribution at small scales, responsible for an enhancement with respect to the smoothly distributed case. We show how a direct computation of this "flux multiplier" in terms of the nonlinear power spectrum offers a conceptually simpler approach and may ease some problems, such as the extrapolation issue. In fact very simple analytical recipes to construct the power spectrum yield results similar to the popular Halo Model expectations, with a straightforward alternative estimate of errors. For this specific application, one also obviates to the need of identifying (often literature-dependent) concepts entering the Halo Model, to compare different simulations.

We present measurements of the bispectrum of dark matter halos in numerical simulations with non-Gaussian initial conditions of the local type. We show, in the first place, that the overall effect of primordial non-Gaussianity on the halo bispectrum is larger than on the halo power spectrum when all measurable configurations are taken into account. We then compare our measurements with a tree-level perturbative prediction finding good agreement at large scale when the constant Gaussian bias parameter, both linear and quadratic, and their constant non-Gaussian corrections are fitted for. The best-fit values of the Gaussian bias factors and their non-Gaussian, scale-independent corrections are in qualitative agreement with the peak-background split expectations. In particular, we show that the effect of non-Gaussian initial conditions on squeezed configurations is fairly large (up to 30% for f_NL=100 at redshift z=0.5) and results from contributions of similar amplitude induced by the initial matter bispectrum, scale-dependent bias corrections as well as from nonlinear matter bispectrum corrections. We show, in addition, that effects at second order in f_NL are irrelevant for the range of values allowed by CMB and galaxy power spectrum measurements, at least on the scales probed by our simulations. Finally, we present a Fisher matrix analysis to assess the possibility of constraining primordial non-Gaussianity with future measurements of the galaxy bispectrum. We find that a survey with a volume of about 10 cubic Gpc at mean redshift z ~ 1 could provide an error on f_NL of the order of a few. This shows the relevance of a joint analysis of galaxy power spectrum and bispectrum in future redshift surveys.

We study the nonlinear evolution of density perturbations in cosmologies where the late-time accelerated expansion is driven by a quintessence field with vanishing speed of sound. For these models matter and quintessence perturbations are comoving and it is possible to write a single continuity equation for the total density fluctuations given by a weighted sum of the two components. Including the Euler equation for the common velocity field we solve the evolution equations for the nonlinear, total density power spectrum in the Time-Renormalization Group approach. In fact any cosmological observable is directly related by gravity only to the total density perturbations, with the two components being individually unobservable. We show that the clustering of quintessence perturbations induces small corrections with respect to the nonlinear evolution of power spectrum in smooth quintessence models described by the same equation of state. Such small corrections, however, contrast with the large effect of a vanishing speed of sound on the linear growth function at low redshift. For this reason, models with the same normalization of the linear density power spectrum can present significantly different nonlinear corrections depending on the value of the sound speed. Although such differences vanish in the w \rightarrow -1, {\Lambda}CDM limit, we argue that the relation between linear and nonlinear growth of structures should be properly taken into account in constraining models with inhomogeneous dark energy.

Estimators for weak lensing observables such as shear and convergence generally have non-linear corrections, which, in principle, make weak lensing power spectra sensitive to primordial non-Gaussianity. In this paper, we quantitatively evaluate these contributions for weak lensing auto- and cross-correlation power spectra, and show that they are strongly suppressed by projection effects. This is a consequence of the central limit theorem, which suppresses departures from Gaussianity when the projection reaches over several correlation lengths of the density field, L_P~55 [Mpc/h]. Furthermore, the typical scales that contribute to projected bispectra are generally smaller than those that contribute to projected power spectra. Both of these effects are not specific to lensing, and thus affect the statistics of non-linear tracers (e.g., peaks) of any projected density field. Thus, the clustering of biased tracers of the three-dimensional density field is generically more sensitive to non-Gaussianity than observables constructed from projected density fields.

We study large-scale structure formation in the presence of a quintessence component with zero speed of sound in the framework of Eulerian Perturbation Theory. Due to the absence of pressure gradients, quintessence and dark matter are comoving and can be studied as a unique fluid in terms of the total energy density contrast and the common velocity. In this description the clustering of quintessence enhances the linear term proportional to the velocity divergence in the continuity equation by a factor (1+w) Omega_Q / Omega_m. This is responsible for a rapid evolution of the growth rate at low redshifts, and modifies the standard relation between the velocity divergence and the growth factor. For the total fluid, the solutions for the linear growth function and growth rate can be written in integral forms and admit simple fitting formulae, as in the LambdaCDM case. At second order in perturbation theory, we derive an explicit expression for the kernels F_2 and G_2. They receive modifications of the order of the ratio between quintessence and total energy density perturbations, which affect the corresponding tree-level bispectra. We finally compute the cumulative signal-to-noise in the power spectrum, bispectrum and reduced bispectrum, expected for departures from a LambdaCDM cosmology both in the clustering and smooth quintessence scenarios. The reduced bispectrum, in particular, receives sensible modifications only in the clustering case and can potentially be used to detect or rule out the model.

We show here how Renormalized Perturbation Theory (RPT) calculations applied to the quasi-linear growth of the large-scale structure can be carried on in presence of primordial non-Gaussian (PNG) initial conditions. It is explicitly demonstrated that the series reordering scheme proposed in Bernardeau, Crocce and Scoccimarro (2008) is preserved for non-Gaussian initial conditions. This scheme applies to the power spectrum and higher order spectra and is based on a reorganization of the contributing terms into sum of products of multi-point propagators. In case of PNG new contributing terms appear, the importance of which is discussed in the context of current PNG models. The properties of the building blocks of such resummation schemes, the multi-point propagators, are then investigated. It is first remarked that their expressions are left unchanged at one-loop order irrespectively of statistical properties of the initial field. We furthermore show that the high-momemtum limit of each of these propagators can be explicitly computed even for arbitrary initial conditions. They are found to be damped by an exponential cutoff whose expression is directly related to the moment generating function of the one-dimensional displacement field. This extends what had been established for multi-point propagators for Gaussian initial conditions. Numerical forms of the cut-off are shown for the so-called local model of PNG.

We present measurements of the dark matter bispectrum in N-body simulations with non-Gaussian initial conditions of the local kind for a large variety of triangular configurations and compare them with predictions from Eulerian Perturbation Theory up to one-loop corrections. We find that the effects of primordial non-Gaussianity at large scales, when compared to Perturbation Theory, are well described by the initial component of the matter bispectrum, linearly extrapolated at the redshift of interest. In addition, we find that, for f_NL=100, the nonlinear corrections due to non-Gaussian initial conditions are of the order of ~3, 4% for generic triangles up to ~20% for squeezed configurations, at any redshift. We show that the predictions of Perturbation Theory at tree-level fail to describe the simulation results at redshift z=0 already at scales corresponding to k ~ 0.02 - 0.08 h/Mpc, depending on the triangle, while one-loop corrections can significantly extend their validity to smaller scales. At higher redshift, one-loop Perturbation Theory provides indeed quite accurate predictions, particularly with respect to the relative correction due to primordial non-Gaussianity.

The most direct probe of non-Gaussian initial conditions has come from bispectrum measurements of temperature fluctuations in the Cosmic Microwave Background and of the matter and galaxy distribution at large scales. Such bispectrum estimators are expected to continue to provide the best constraints on the non-Gaussian parameters in future observations. We review and compare the theoretical and observational problems, current results and future prospects for the detection of a non-vanishing primordial component in the bispectrum of the Cosmic Microwave Background and large-scale structure, and the relation to specific predictions from different inflationary models.

We derive the expressions for the 1-loop corrections in cosmological, Eulerian, perturbation theory to the matter bispectrum and to the galaxy bispectrum, assuming local galaxy bias, in presence of non-Gaussian initial conditions. We compute them explicitly for the particular case of non-vanishing initial bispectrum and trispectrum in the local model and for a non-vanishing initial bispectrum alone for the equilateral model of primordial non-Gaussianity. While the primordial contribution to the matter bispectrum for values compatible with CMB observations is dominant over the component due to gravitational instability at large scales, 1-loop perturbative corrections due to non-Gaussian initial conditions correspond to just a few percent of the gravity-induced bispectrum at mildly non-linear scales, similarly to what happens for the matter power spectrum. However, in the perturbative expansion for the galaxy bispectrum, 1-loop diagrams arising from non-linear bias are responsible for significant large-scale contributions, indeed exceeding the primordial component, both for the local and equilateral model. We study the peculiar dependence on scale and on the shape of the triangular configurations of such additional terms, similar in their origin to the large-scale corrections to the halo and galaxy power spectra that raised significant interest in the recent literature.

The primordial non-Gaussian parameter fNL has been shown to be scale-dependent in several models of inflation with a variable speed of sound. Starting from a simple ansatz for a scale-dependent amplitude of the primordial curvature bispectrum for two common phenomenological models of primordial non-Gaussianity, we perform a Fisher matrix analysis of the bispectra of the temperature and polarization of the Cosmic Microwave Background (CMB) radiation and derive the expected constraints on the parameter nNG that quantifies the running of fNL(k) for current and future CMB missions such as WMAP, Planck and CMBPol. We find that CMB information alone, in the event of a significant detection of the non-Gaussian component, corresponding to fNL = 50 for the local model and fNL = 100 for the equilateral model of non-Gaussianity, is able to determine nNG with a 1-sigma uncertainty of Delta nNG = 0.1 and Delta nNG = 0.3, respectively, for the Planck mission. In addition, we consider a Fisher matrix analysis of the galaxy power spectrum to determine the expected constraints on the running parameter nNG for the local model and of the galaxy bispectrum for the equilateral model from future photometric and spectroscopic surveys. We find that, in both cases, large-scale structure observations should achieve results comparable to or even better than those from the CMB, while showing some complementarity due to the different distribution of the non-Gaussian signal over the relevant range of scales. Finally, we compare our findings to the predictions on the amplitude and running of non-Gaussianity of DBI inflation, showing how the constraints on a scale-dependent fNL(k) translate into constraints on the parameter space of the theory.

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

We measure the two-point spatial correlation function for clusters selected from the photometric MaxBCG galaxy cluster catalog for the Sloan Digital Sky Survey (SDSS). We evaluate the correlation function for several cluster samples using different cuts in cluster richness. Fitting the results to power-laws, $\xi_{cc}(r) = (r/R_0)^{-\gamma}$, the estimated correlation length $R_0$ as a function of richness is broadly consistent with previous cluster observations and with expectations from N-body simulations. We study how the linear bias parameter scales with richness and compare our results to theoretical predictions. Since these measurements extend to very large scales, we also compare them to models that include the baryon acoustic oscillation feature and that account for the smoothing effects induced by errors in the cluster photometric redshift estimates. For the largest cluster sample, corresponding to a richness threshold of $\Ng\ge 10$, we find only weak evidence, of about $1.4-1.7\sigma$ significance, for the baryonic acoustic oscillation signature in the cluster correlation function.

We study the corrections to the galaxy three-point correlation function (3PCF) induced by weak lensing magnification due to the matter distribution along the line of sight. We consistently derive all the correction terms arising up to second order in perturbation theory and provide analytic expressions as well as order of magnitude estimates for their relative importance. The magnification contributions depend on the geometry of the projected triangle on the sky plane, and scale with different powers of the number count slope and redshift of the galaxy sample considered. We evaluate all terms numerically and show that, depending on the triangle configuration as well as the galaxy sample considered, weak lensing can in general significantly contribute to and alter the three-point correlation function observed through galaxy and quasar catalogs.

The greatest challenge in the interpretation of galaxy clustering data from any surveys is galaxy bias. Using a simple Fisher matrix analysis, we show that the bispectrum provides an excellent determination of linear and non-linear bias parameters of intermediate and high-z galaxies, when all measurable triangle configurations down to mildly non-linear scales, where perturbation theory is still valid, are included. The bispectrum is also a powerful probe of primordial non-Gaussianity. The planned galaxy surveys at z>2 should yield constraints on non-Gaussian parameters, f_{NL}^{loc.} and f_{NL}^{eq.}, that are comparable to, or even better than, those from CMB experiments. We study how these constraints improve with volume, redshift range, as well as the number density of galaxies. Finally, we show that a halo occupation distribution may be used to improve these constraints further by lifting degeneracies between gravity, bias, and primordial non-Gaussianity.

Galaxy cluster surveys will be a powerful probe of dark energy. At the same time, cluster abundance is sensitive to any non-Gaussianity of the primordial density field. It is therefore possible that non-Gaussian initial conditions might be misinterpreted as a sign of dark energy or at least degrade the expected constraints on dark energy parameters. To address this issue, we perform a likelihood analysis of an ideal cluster survey similar in size and depth to the upcoming South Pole Telescope/Dark Energy Survey (SPT-DES). We analyze a model in which the strength of the non-Gaussianity is parameterized by the constant fNL; this model has been used extensively to derive Cosmic Microwave Background (CMB) anisotropy constraints on non-Gaussianity, allowing us to make contact with those works. We find that the constraining power of the cluster survey on dark energy observables is not significantly diminished by non-Gaussianity provided that cluster redshift information is included in the analysis. We also find that even an ideal cluster survey is unlikely to improve significantly current and future CMB constraints on non-Gaussianity. However, when all systematics are under control, it could constitute a valuable cross check to CMB observations.

The present spatial distribution of galaxies in the Universe is non-Gaussian, with 40% skewness in 50 Mpc/h spheres, and remarkably little is known about the information encoded in it about cosmological parameters beyond the power spectrum. In this work we present an attempt to bridge this gap by studying the bispectrum, paying particular attention to a joint analysis with the power spectrum and their combination with CMB data. We address the covariance properties of the power spectrum and bispectrum including the effects of beat coupling that lead to interesting cross-correlations, and discuss how baryon acoustic oscillations break degeneracies. We show that the bispectrum has significant information on cosmological parameters well beyond its power in constraining galaxy bias, and when combined with the power spectrum is more complementary than combining power spectra of different samples of galaxies, since non-Gaussianity provides a somewhat different direction in parameter space. In the framework of flat cosmological models we show that most of the improvement of adding bispectrum information corresponds to parameters related to the amplitude and effective spectral index of perturbations, which can be improved by almost a factor of two. Moreover, we demonstrate that the expected statistical uncertainties in sigma8 of a few percent are robust to relaxing the dark energy beyond a cosmological constant.

We show that current surveys have at least as much signal to noise in higher-order statistics as in the power spectrum at weakly nonlinear scales. We discuss how one can use this information to determine the mean of the galaxy halo occupation distribution (HOD) using only large-scale information, through galaxy bias parameters determined from the galaxy bispectrum and trispectrum. After introducing an averaged, reasonably fast to evaluate, trispectrum estimator, we show that the expected errors on linear and quadratic bias parameters can be reduced by at least 20-40%. Also, the inclusion of the trispectrum information, which is sensitive to "three-dimensionality" of structures, helps significantly in constraining the mass dependence of the HOD mean. Our approach depends only on adequate modeling of the abundance and large-scale clustering of halos and thus is independent of details of how galaxies are distributed within halos. This provides a consistency check on the traditional approach of using two-point statistics down to small scales, which necessarily makes more assumptions. We present a detailed forecast of how well our approach can be carried out in the case of the SDSS.

We consider primordial non-Gaussianity due to quadratic corrections in the gravitational potential parametrized by a non-linear coupling parameter fnl. We study constraints on fnl from measurements of the galaxy bispectrum in redshift surveys. Using estimates for idealized survey geometries of the 2dF and SDSS surveys and realistic ones from SDSS mock catalogs, we show that it is possible to probe |fnl|~100, after marginalization over bias parameters. We apply our methods to the galaxy bispectrum measured from the PSCz survey, and obtain a 2sigma-constraint |fnl|< 1800. We estimate that an all sky redshift survey up to z~1 can probe |fnl|~1. We also consider the use of cluster abundance to constrain fnl and find that in order to be sensitive to |fnl|~100, cluster masses need to be determined with an accuracy of a few percent, assuming perfect knowledge of the mass function and cosmological parameters.