Emiliano Sefusatti

The bispectrum of galaxies from high-redshift galaxy surveys: Primordial non-Gaussianity and non-linear galaxy bias

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 {approx}> 2 should yield constraints on non-Gaussian parameters, f{sub NL}{sup loc.} and f{sub NL}{sup 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.

Primordial Non-Gaussianity and Bispectrum Measurements in the Cosmic Microwave Background and Large-Scale Structure

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 nonvanishing primordial component in the bispectrum of the Cosmic Microwave Background and large-scale structure, and the relation to specific predictions from different inflationary models.

Constraining running non-gaussianity

The primordial non-Gaussian parameter fNL has been shown to be scale-dependent in several models of inflation with a variable speed of sound, such as Dirac-Born-Infeld (DBI) models. 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. 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-σ uncertainty of nNG 0.1 and ΔnNG 0.3, respectively, for the Planck mission and a factor of two better for CMBPol. In addition, we show how future 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.

The matter bispectrum in N-body simulations with non-Gaussian initial conditions

We present measurements of the bispectrum of dark matter haloes in numerical simulations with non-Gaussian initial conditions of 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 scales when the constant Gaussian bias parameter, both linear and quadratic, and their constant non-Gaussian corrections are fitted for. The best-fitting 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 per cent for fNL = 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 non-linear matter bispectrum corrections. We show, in addition, that effects at second order in fNL are irrelevant for the range of values allowed by cosmic microwave background and galaxy power spectrum measurements, at least on the scales probed by our simulations (k > 0.01 h Mpc−1). 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 h−3 Gpc3 at mean redshift z ≃ 1 could provide an error on fNL of the order of a few. This shows the relevance of a joint analysis of galaxy power spectrum and bispectrum in future redshift surveys.

The Correlation Function of Optically Selected Galaxy Clusters in the Sloan Digital Sky Survey

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, ?cc(r) = (r/R 0)??, 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 N 200 ? 10, we find only weak evidence, of about 1.4?-1.7? significance, for the baryonic acoustic oscillation signature in the cluster correlation function.

Primordial non-gaussianity and dark energy constraints from cluster surveys

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 survey and Dark Energy Survey. 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 significantly improve current and future CMB constraints on non-Gaussianity. However, when all systematics are under control, such surveys could constitute a valuable cross-check on CMB observations.

Initial Conditions for Accurate N-Body Simulations of Massive Neutrino Cosmologies

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 characterizes 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 per cent accuracy or below for all redshift relevant to non-linear 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 (REPS rescaled power spectra for initial conditions with massive neutrinos https://github.com/matteozennaro/reps) providing the initial displacements and velocities for cold dark matter and neutrino particles that will allow accurate, i.e. 1 per cent level, numerical simulations for this cosmological scenario.

Improving fast generation of halo catalogues with higher order Lagrangian perturbation theory

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

Weak Lensing Effects on the Galaxy Three-Point Correlation Function

k_{max}\sim0.5\ h/

DEMNUni: Massive neutrinos and the bispectrum of large scale structures

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.