BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Vincent Desjacques

Vincent Desjacques appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
astro-ph.COhep-phgr-qchep-thastro-ph.GAastro-ph.HEastro-phastro-ph.SR

Catalog footprint

What is connected

58works
8topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

58 published item(s)

preprint2023arXiv

Assembly bias in the local PNG halo bias and its implication for $f_{\rm NL}$ constraints

We use $N$-body simulations to study halo assembly bias (i.e., the dependence of halo clustering on properties beyond total mass) in the density and primordial non-Gaussianity (PNG) linear bias parameters $b_1$ and $b_ϕ$, respectively. We consider concentration, spin and sphericity as secondary halo properties, for which we find a clear detection of assembly bias for $b_1$ and $b_ϕ$. At fixed total mass, halo spin and sphericity impact $b_1$ and $b_ϕ$ in a similar manner, roughly preserving the shape of the linear $b_ϕ(b_1)$ relation satisfied by the global halo population. Halo concentration, however, drives $b_1$ and $b_ϕ$ in opposite directions. This induces significant changes to the $b_ϕ(b_1)$ relation, with higher concentration halos having higher amplitude of $b_ϕ(b_1)$. For $z=0.5$ and $b_1 \approx 2$ in particular, the population comprising either all halos, those with the $33\%$ lowest or those with the $33\%$ highest concentrations have a PNG bias of $b_ϕ\approx 3$, $b_ϕ\approx -1$ and $b_ϕ\approx 9$, respectively. Varying the halo concentration can make $b_ϕ$ very small and even change its sign. These results have important ramifications for galaxy clustering constraints of the local PNG parameter $f_{\rm NL}$ that assume fixed forms for the $b_ϕ(b_1)$ relation. We illustrate the significant impact of halo assembly bias in actual data using the BOSS DR12 galaxy power spectrum: assuming that BOSS galaxies are representative of all halos, the $33\%$ lowest or the $33\%$ highest concentration halos yields $σ_{f_{\rm NL}} = 44, 165, 19$, respectively. Our results suggest taking host halo concentration into account in galaxy selection strategies to maximize the signal-to-noise on $f_{\rm NL}$. They also motivate more simulation-based efforts to study the $b_ϕ(b_1)$ relation of halos and galaxies.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Analytic solution to the dynamical friction acting on circularly moving perturbers

We present an analytic approach to the dynamical friction (DF) acting on a circularly moving point mass perturber in a gaseous medium. We demonstrate that, when the perturber is turned on at $t=0$, steady-state (infinite time perturbation) is achieved after exactly one sound-crossing time. At low Mach number $\mathcal{M}~\ll~1$, the circular-motion steady-state DF converges to the linear-motion, finite time perturbation expression. The analytic results describe both the radial and tangential forces on the perturbers caused by the backreaction of the wake propagating in the medium. The radial force is directed inward, toward the motion centre, and is dominant at large Mach numbers. For subsonic motion, this component is negligible. For moderate and low Mach numbers, the tangential force is stronger and opposes the motion of the perturber. The analytic solution to the circular-orbit DF suffers from a logarithmic divergence in the supersonic regime. This divergence appears at short distances from the perturber solely (unlike the linear motion result which is also divergent at large distances) and can be encoded in a maximum multipole. This is helpful to assess the resolution dependence of numerical simulations implementing DF at the level of Liénard-Wiechert potentials. We also show how our approach can be generalised to calculate the DF acting on a compact circular binary.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Axion core - halo mass and the black hole - halo mass relation: constraints on a few parsec scales

If the dark matter is made of ultra-light axions, stable solitonic cores form at the centers of virialized halos. In some range for the mass $m$ of the axion particle, these cores are sufficiently compact and can mimic supermassive black holes (SMBH) residing at galactic nuclei. We use the solitonic core--halo mass relation, validated in numerical simulations, to constrain a new range of allowed axion mass from measurements of the SMBH mass in (pseudo)bulge and bulgeless galaxies. These limits are based on observations of galactic nuclei on scales smaller than 10 pc. Our analysis suggests that $m < 10^{-18}$ eV is ruled out by the data. We briefly discuss whether an attractive self-interaction among axions could alleviate this constraint.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Axion resonances in binary pulsar systems

We investigate the extent to which resonances between an oscillating background of ultra-light axion and a binary Keplerian system can affect the motion of the latter. These resonances lead to perturbations in the instantaneous time-of-arrivals, and to secular variations in the period of the binary. While the secular changes at exact resonance have recently been explored, the instantaneous effects have been overlooked. In this paper, we examine the latter using N-body simulations including the external oscillatory forcing induced by the axion background. While the secular effects are restricted to a narrow width near the resonance, the instantaneous changes, albeit strongest close to resonances, are apparent for wide range of configurations. We compute the signal-to-noise ratio (SNR) as a function of semi-major axis for a detection of axion oscillations through the R\{o} mer delay. The latter can be extracted from the time-of-arrivals of binary pulsars. The SNR broadly increases with increasing binary eccentricity in agreement with the secular expectation. However, we find that it differs significantly from the scaling a^{5/2} around the lowest orders of resonance. Future observations could probe these effects away from resonances and, therefore, constrain a much broader range of axion masses provided that binary pulsar systems are found near the central region of our Galaxy, and that the time-or-arrival measurement accuracy reaches < 10 ns

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Shot noise in multi-tracer constraints on $f_\text{NL}$ and relativistic projections: Power Spectrum

Multiple tracers of the same surveyed volume can enhance the signal-to-noise on a measurement of local primordial non-Gaussianity and the relativistic projections. Increasing the number of tracers comparably increases the number of shot noise terms required to describe the stochasticity of the data. Although the shot noise is white on large scales, it is desirable to investigate the extent to which it can degrade constraints on the parameters of interest. In a multi-tracer analysis of the power spectrum, a marginalization over shot noise does not degrade the constraints on $f_\text{NL}$ by more than $\sim 30$% so long as halos of mass $M\lesssim 10^{12}M_\odot$ are resolved. However, ignoring cross shot noise terms induces large systematics on a measurement of $f_\text{NL}$ at redshift $z<1$ when small mass halos are resolved. These effects are less severe for the relativistic projections, especially for the dipole term. In the case of a low and high mass tracer, the optimal sample division maximizes the signal-to-noise on $f_\text{NL}$ and the projection effects simultaneously, reducing the errors to the level of $\sim 10$ consecutive mass bins of equal number density. We also emphasize that the non-Poissonian noise corrections that arise from small-scale clustering effects cannot be measured with random dilutions of the data. Therefore, they must either be properly modeled or marginalized over.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Gravitational waves from in-spirals of compact objects in binary common-envelope evolution

Detection of gravitational-wave (GW) sources enables the characterisation of binary compact objects and of their in-spiral. However, other dissipative processes can affect the in-spiral. Here we show that the in-spiral of compact objects through a gaseous common-envelope (CE) arising from an evolved stellar companion produces a novel type of GW-sources, whose evolution is dominated by the dissipative gas dynamical friction effects from the CE, rather than the GW-emission itself. The evolution and properties of the GW-signals differ from those of isolated gas-poor mergers significantly. We find characteristic strains of $\sim10^{-23}$-$10^{-21}$ ($10{\rm kpc}/{D}$) for such sources -- observable by next-generation space-based GW-detectors. The evolution of the GW-signal can serve as a probe of the interior regions of the evolved star, and the final stages of CE-evolution, otherwise inaccessible through other observational means. Moreover, such CE-mergers are frequently followed by observable explosive electromagnetic counterparts and/or the formation of exotic stars.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

The information content of Cosmic Infrared Background anisotropies

We use analytic computations to predict the power spectrum as well as the bispectrum of Cosmic Infrared Background (CIB) anisotropies. Our approach is based on the halo model and takes into account the mean luminosity-mass relation. The model is used to forecast the possibility to simultaneously constrain cosmological, CIB and halo occupation distribution (HOD) parameters in the presence of foregrounds. For the analysis we use wavelengths in eight frequency channels between 200 and 900$\;\mathrm{GHz}$ with survey specifications given by Planck and LiteBird. We explore the sensitivity to the model parameters up to multipoles of $\ell =1000$ using auto- and cross-correlations between the different frequency bands. With this setting, cosmological, HOD and CIB parameters can be constrained to a few percent. Galactic dust is modeled by a power law and the shot noise contribution as a frequency dependent amplitude which are marginalized over. We find that dust residuals in the CIB maps only marginally influence constraints on standard cosmological parameters. Furthermore, the bispectrum yields tighter constraints (by a factor four in $1σ$ errors) on almost all model parameters while the degeneracy directions are very similar to the ones of the power spectrum. The increase in sensitivity is most pronounced for the sum of the neutrino masses. Due to the similarity of degeneracies a combination of both analysis is not needed for most parameters. This, however, might be due to the simplified bias description generally adopted in such halo model approaches.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Cosmic Infrared Background anisotropies as a window into primordial non-Gaussianity

The angular power spectrum of the cosmic infrared background (CIB) is a sensitive probe of the local primordial bispectrum. CIB measurements are integrated over a large volume so that the scale dependent bias from the primordial non-Gaussianity leaves a strong signal in the CIB power spectrum. Although galactic dust dominates over the non-Gaussian CIB signal, it is possible to mitigate the dust contamination with enough frequency channels, especially if high frequencies such as the Planck 857 GHz channel are available. We show that, in this case, measurements of the cosmic microwave background from future space missions should be able to probe the local bispectrum shape down to an amplitude |f_nl| < 1.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Impacts of biasing schemes in the one-loop integrated perturbation theory

The impact of biasing schemes on the clustering of tracers of the large-scale structure is analytically studied in the weakly nonlinear regime. For this purpose, we use the one-loop approximation of the integrated perturbation theory together with the renormalized bias functions of various, physically motivated Lagrangian bias schemes. These include the halo, peaks and excursion set peaks model, for which we derive useful formulae for the evaluation of their renormalized bias functions. The shapes of the power spectra and correlation functions are affected by the different bias models at the level of a few percent on weakly nonlinear scales. These effects are studied quantitatively both in real and redshift space. The amplitude of the scale-dependent bias in the presence of primordial non-Gaussianity also depends on the details of the bias models. If left unaccounted for, these theoretical uncertainties could affect the robustness of the cosmological constraints extracted from galaxy clustering data.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Lagrangian bias of generic large-scale structure tracers

The dark matter halos that host galaxies and clusters form out of initial high-density patches, providing a biased tracer of the linear matter density field. In the simplest local bias approximation, the halo field is treated as a perturbative series in the average overdensity of the Lagrangian patch. In more realistic models, however, additional quantities will affect the clustering of halo-patches, and this expansion becomes a function of several stochastic variables. In this paper, we present a general multivariate expansion scheme that can parametrize the clustering of any biased Lagrangian tracer, given only the variables involved and their symmetry (in our case rotational invariance). This approach is based on an expansion in the orthonormal polynomials associated with the relevant variables, so that no renormalization of the coefficients ever occurs. We provide explicit expression for the series coefficients, or Lagrangian bias parameters, in the case of peaks of the linear density field. As an application of our formalism, we present a simple derivation of the original BBKS formula, and compute the non-Gaussian bias in the presence of a primordial trispectrum of the local shape.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Squeezing the halo bispectrum: a test of bias models

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The Halo Boltzmann Equation

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Velocity bias in the distribution of dark matter halos

The standard formalism for the co-evolution of halos and dark matter predicts that any initial halo velocity bias rapidly decays to zero. We argue that, when the purpose is to compute statistics like power spectra etc., the coupling in the momentum conservation equation for the biased tracers must be modified. Our new formulation predicts the constancy in time of any statistical halo velocity bias present in the initial conditions, in agreement with peak theory. We test this prediction by studying the evolution of a conserved halo population in N-body simulations. We establish that the initial simulated halo density and velocity statistics show distinct features of the peak model and, thus, deviate from the simple local Lagrangian bias. We demonstrate, for the first time, that the time evolution of their velocity is in tension with the rapid decay expected in the standard approach.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Detecting the cosmological recombination signal from space

Spectral distortions of the CMB have recently experienced an increased interest. One of the inevitable distortion signals of our cosmological concordance model is created by the cosmological recombination process, just a little before photons last scatter at redshift $z\simeq 1100$. These cosmological recombination lines, emitted by the hydrogen and helium plasma, should still be observable as tiny deviation from the CMB blackbody spectrum in the cm--dm spectral bands. In this paper, we present a forecast for the detectability of the recombination signal with future satellite experiments. We argue that serious consideration for future CMB experiments in space should be given to probing spectral distortions and, in particular, the recombination line signals. The cosmological recombination radiation not only allows determination of standard cosmological parameters, but also provides a direct observational confirmation for one of the key ingredients of our cosmological model: the cosmological recombination history. We show that, with present technology, such experiments are futuristic but feasible. The potential rewards won by opening this new window to the very early universe could be considerable.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Galaxy power spectrum in redshift space: combining perturbation theory with the halo model

Theoretical modeling of the redshift-space power spectrum of galaxies is crucially important to correctly extract cosmological information from redshift surveys. The task is complicated by the nonlinear biasing and redshift space distortion effects, which change with halo mass, and by the wide distribution of halo masses and their occupations by galaxies. One of the main modeling challenges is the existence of satellite galaxies that have both radial distribution and large virial velocities inside halos, a phenomenon known as the Finger-of-God effect. We present a model for the galaxy power spectrum of in which we decompose a given galaxy sample into central and satellite galaxies and relate different contributions to 1- and 2-halo terms in a halo model. Our primary goal is to ensure that any parameters that we introduce have physically meaningful values, and are not just fitting parameters. For the 2-halo terms we use the previously developed RSD modeling of halos in the context of distribution function and perturbation theory approach. This term needs to be multiplied by the effect of radial distances and velocities of satellites inside the halo. To this one needs to add the 1-halo terms, which are non-perturbative. We show that the real space 1-halo terms can be modeled as almost constant, with the finite extent of the satellites inside the halo inducing a small k^2R^2 term, where R is related to the size of the halo. We adopt a similar model for FoG in redshift space, ensuring that FoG velocity dispersion is related to the halo mass. For FoG k^2 type expansions do not work and FoG resummation must be used instead. We test several damping functions to model the velocity dispersion FoG effect. Applying the formalism to mock galaxies modeled after the "CMASS" sample of the BOSS survey, we find that our predictions for the redshift-space power spectra are accurate up to k~0.4Mpc/h.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Intrinsic alignment of simulated galaxies in the cosmic web: implications for weak lensing surveys

The intrinsic alignment of galaxy shapes (by means of their angular momentum) and their cross-correlation with the surrounding dark matter tidal field are investigated using the 160 000, z=1.2 synthetic galaxies extracted from the high-resolution cosmological hydrodynamical simulation Horizon-AGN. One- and two-point statistics of the spin of the stellar component are measured as a function of mass and colour. For the low-mass galaxies, this spin is locally aligned with the tidal field `filamentary' direction while, for the high-mass galaxies, it is perpendicular to both filaments and walls. The bluest galaxies of our synthetic catalog are more strongly correlated with the surrounding tidal field than the reddest galaxies, and this correlation extends up to 10 Mpc/h comoving distance. We also report a correlation of the projected ellipticities of blue, intermediate mass galaxies on a similar scale at a level of 10^(-4) which could be a concern for cosmic shear measurements. We do not report any measurable intrinsic alignments of the reddest galaxies of our sample. This work is a first step toward the use of very realistic catalog of synthetic galaxies to evaluate the contamination of weak lensing measurement by the intrinsic galactic alignments.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Peak exclusion, stochasticity and convergence of perturbative bias expansions in 1+1 gravity

The Lagrangian peaks of a 1D cosmological random field representing dark matter are used as a proxy for a catalogue of biased tracers in order to investigate the small-scale exclusion in the two-halo term. The two-point correlation function of peaks of a given height is numerically estimated and analytical approximations that are valid inside the exclusion zone are derived. The resulting power spectrum of these tracers is investigated and shows clear deviations from Poisson noise at low frequencies. On large scales, the convergence of a perturbative bias expansion is discussed. Finally, we go beyond Gaussian statistics for the initial conditions and investigate the subsequent evolution of the two-point clustering of peaks through their Zel'dovich ballistic displacement, to clarify how exclusion effects mix up with scale-dependencies induced by nonlinear gravitational evolution. While the expected large-scale separation limit is recovered, significant deviations are found in the exclusion zone that tends in particular to be reduced at later times. Even though these findings apply to the clustering of one-dimensional tracers, they provide useful insights into halo exclusion and its impact on the two-halo term.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Dancing in the dark: galactic properties trace spin swings along the cosmic web

A large-scale hydrodynamical cosmological simulation, Horizon-AGN, is used to investigate the alignment between the spin of galaxies and the cosmic filaments above redshift 1.2. The analysis of more than 150 000 galaxies per time step in the redshift range 1.2<z<1.8 with morphological diversity shows that the spin of low-mass blue galaxies is preferentially aligned with their neighbouring filaments, while high-mass red galaxies tend to have a perpendicular spin. The reorientation of the spin of massive galaxies is provided by galaxy mergers, which are significant in their mass build-up. We find that the stellar mass transition from alignment to misalignment happens around 3.10^10 M_sun. Galaxies form in the vorticity-rich neighbourhood of filaments, and migrate towards the nodes of the cosmic web as they convert their orbital angular momentum into spin. The signature of this process can be traced to the properties of galaxies, as measured relative to the cosmic web. We argue that a strong source of feedback such as active galactic nuclei is mandatory to quench in situ star formation in massive galaxies and promote various morphologies. It allows mergers to play their key role by reducing post-merger gas inflows and, therefore, keeping spins misaligned with cosmic filaments.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Large-Scale Clustering of Cosmic Voids

We study the clustering of voids using $N$-body simulations and simple theoretical models. The excursion-set formalism describes fairly well the abundance of voids identified with the watershed algorithm, although the void formation threshold required is quite different from the spherical collapse value. The void cross bias $b_{\rm c} $ is measured and its large-scale value is found to be consistent with the peak background split results. A simple fitting formula for $b_{\rm c} $ is found. We model the void auto-power spectrum taking into account the void biasing and exclusion effect. A good fit to the simulation data is obtained for voids with radii $\gtrsim$ 30 Mpc/$h$, especially when the void biasing model is extended to 1-loop order. However, the best-fit bias parameters do not agree well with the peak-background split results. Being able to fit the void auto-power spectrum is particularly important not only because it is the direct observable in galaxy surveys, but also our method enables us to treat the bias parameters as nuisance parameters, which are sensitive to the techniques used to identify voids.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Measuring nonlocal Lagrangian peak bias

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 χ_{10} if we allow for an offset between the Lagrangian halo center-of-mass and the peak position. The sign and magnitude of χ_{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 χ_{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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Non-local halo bias with and without massive neutrinos

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_{ν_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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

On the halo velocity bias

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Peculiar velocities in redshift space: formalism, N-body simulations and perturbation theory

Direct measurements of peculiar velocities of galaxies and clusters of galaxies can in principle provide explicit information on the three dimensional mass distribution, but this information is modulated by the fact that velocity field is sampled at galaxy positions, and is thus probing galaxy momentum. We derive expressions for the cross power spectrum between the density and momentum field and the auto spectrum of the momentum field in redshift space, by extending the distribution function method to these statistics. The resulting momentum cross and auto power spectra in redshift space are expressed as infinite sums over velocity moment correlators in real space, as is the case for the density power spectrum in redshift space. We compute each correlator using Eulerian perturbation theory (PT) and halo biasing model and compare the resulting redshift-space velocity statistics to those measured from N-body simulations for both dark matter and halos. We find that in redshift space linear theory predictions for the density-momentum cross power spectrum as well as for the momentum auto spectrum fail to predict the N-body results at very large scales. On the other hand, our nonlinear PT prediction for these velocity statistics, together with real-space power spectrum for dark matter from simulations, improves the accuracy for both dark matter and halos. We also present the same analysis in configuration space, computing the redshift-space pairwise mean infall velocities and velocity correlation function and compare to nonlinear PT.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Ultraviolet background fluctuations with clustered sources

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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

A local bias approach to the clustering of discrete density peaks

Maxima of the linear density field form a point process that can be used to understand the spatial distribution of virialized halos that collapsed from initially overdense regions. However, owing to the peak constraint, clustering statistics of discrete density peaks are difficult to evaluate. For this reason, local bias schemes have received considerably more attention in the literature thus far. In this paper, we show that the 2-point correlation function of maxima of a homogeneous and isotropic Gaussian random field can be thought of, up to second order at least, as arising from a local bias expansion formulated in terms of rotationally invariant variables. This expansion relies on a unique smoothing scale, which is the Lagrangian radius of dark matter halos. The great advantage of this local bias approach is that it circumvents the difficult computation of joint probability distributions. We demonstrate that the bias factors associated with these rotational invariants can be computed using a peak-background split argument, in which the background perturbation shifts the corresponding probability distribution functions. Consequently, the bias factors are orthogonal polynomials averaged over those spatial locations that satisfy the peak constraint. In particular, asphericity in the peak profile contributes to the clustering at quadratic and higher order, with bias factors given by generalized Laguerre polynomials. We speculate that our approach remains valid at all orders, and that it can be extended to describe clustering statistics of any point process of a Gaussian random field. Our results will be very useful to model the clustering of discrete tracers with more realistic collapse prescriptions involving the tidal shear for instance.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

All-Sky Analysis of the General Relativistic Galaxy Power Spectrum

We perform an all-sky analysis of the general relativistic galaxy power spectrum using the well-developed spherical Fourier decomposition. Spherical Fourier analysis expresses the observed galaxy fluctuation in terms of the spherical harmonics and spherical Bessel functions that are angular and radial eigenfunctions of the Helmholtz equation, providing a natural orthogonal basis for all-sky analysis of the large-scale mode measurements. Accounting for all the relativistic effects in galaxy clustering, we compute the spherical power spectrum and its covariance matrix and compare it to the standard three-dimensional power spectrum to establish a connection. The spherical power spectrum recovers the three-dimensional power spectrum at each wavenumber k with its angular dependence mu_k encoded in angular multipole l, and the contributions of the line-of-sight projection to galaxy clustering such as the gravitational lensing effect can be readily accommodated in the spherical Fourier analysis. A complete list of formulas for computing the relativistic spherical galaxy power spectrum is also presented.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Bias deconstructed: Unravelling the scale dependence of halo bias using real space measurements

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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Distribution function approach to redshift space distortions. Part V: perturbation theory applied to dark matter halos

Numerical simulations show that redshift space distortions (RSD) introduce strong scale dependence in the power spectra of halos, with ten percent deviations relative to linear theory predictions even on relatively large scales (k<0.1h/Mpc) and even in the absence of satellites (which induce Fingers-of-God, FoG, effects). If unmodeled these effects prevent one from extracting cosmological information from RSD surveys. In this paper we use perturbation theory (PT) and halo biasing model and apply it to the distribution function approach to RSD, in which RSD is decomposed into several correlators of density weighted velocity moments. We model each of these correlators using PT and compare the results to simulations over a wide range of halo masses and redshifts. We find that with an introduction of a physically motivated halo biasing, and using dark matter power spectra from simulations, we can reproduce the simulation results at a percent level on scales up to k~0.15h/Mpc at z=0, without the need to have free FoG parameters in the model.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Excursion set peaks: a self-consistent model of dark halo abundances and clustering

We describe how to extend the excursion set peaks framework so that its predictions of dark halo abundances and clustering can be compared directly with simulations. These extensions include: a halo mass definition which uses the TopHat filter in real space; the mean dependence of the critical density for collapse delta_c on halo mass m; and the scatter around this mean value. All three of these are motivated by the physics of triaxial rather than spherical collapse. A comparison of the resulting mass function with N-body results shows that, if one uses delta_c(m) and its scatter as determined from simulations, then all three are necessary ingredients for obtaining ~10% accuracy. E.g., assuming a constant value of delta_c with no scatter, as motivated by the physics of spherical collapse, leads to many more massive halos than seen in simulations. The same model is also in excellent agreement with N-body results for the linear halo bias, especially at the high mass end where the traditional peak-background split argument applied to the mass function fit is known to underpredict the measured bias by ~10%. In the excursion set language, our model is about walks centered on special positions (peaks) in the initial conditions -- we discuss what it implies for the usual calculation in which all walks contribute to the statistics.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Halo Stochasticity from Exclusion and non-linear Clustering

The clustering of galaxies in ongoing and upcoming galaxy surveys contains a wealth of cosmological information, but extracting this information is a non-trivial task since galaxies and their host haloes are stochastic tracers of the matter density field. This stochasticity is usually modeled as Poisson shot noise, which is constant as a function of wavenumber with amplitude given by 1/n, where n is the number density of galaxies. Here we use dark matter haloes in N-body simulations to show evidence for deviations from this simple behaviour and develop models that explain the behaviour of the large scale stochasticity. First, haloes are extended, non-overlapping objects, i.e., their correlation function needs to go to -1 on small scales. This leads to a negative correction to the stochasticity relative to the Poisson value at low wavenumber k, decreasing to zero for wavenumbers large compared to the inverse exclusion scale. Second, haloes show a non-linear enhancement of clustering outside the exclusion scale, leading to a positive stochasticity correction. Both of these effects go to zero for high-k, making the stochasticity scale dependent even for k<0.1 h/Mpc. We show that the corrections in the low-k regime are the same in Eulerian and Lagrangian space, but that the transition scale is pushed to smaller scales for haloes observed at present time, relative to the initial conditions. These corrections vary with halo mass and redshift. We also discuss simple applications of these effects to the galaxy samples with non-vanishing satellite fraction, where the stochasticity can again deviate strongly from the fiducial Poisson expectation. Overall these effects affect the clustering of galaxies at a level of a few percent even on very large scales and need to be modelled properly if we want to extract high precision cosmological information from the upcoming galaxy surveys.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

High Performance P3M N-body code: CUBEP3M

This paper presents CUBEP3M, a publicly-available high performance cosmological N-body code and describes many utilities and extensions that have been added to the standard package. These include a memory-light runtime SO halo finder, a non-Gaussian initial conditions generator, and a system of unique particle identification. CUBEP3M is fast, its accuracy is tuneable to optimize speed or memory, and has been run on more than 27,000 cores, achieving within a factor of two of ideal weak scaling even at this problem size. The code can be run in an extra-lean mode where the peak memory imprint for large runs is as low as 37 bytes per particles, which is almost two times leaner than other widely used N-body codes. However, load imbalances can increase this requirement by a factor of two, such that fast configurations with all the utilities enabled and load imbalances factored in require between 70 and 120 bytes per particles. CUBEP3M is well designed to study large scales cosmological systems, where imbalances are not too large and adaptive time-stepping not essential. It has already been used for a broad number of science applications that require either large samples of non-linear realizations or very large dark matter N-body simulations, including cosmological reionization, halo formation, baryonic acoustic oscillations, weak lensing or non-Gaussian statistics. We discuss the structure, the accuracy, known systematic effects and the scaling performance of the code and its utilities, when applicable.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Non-Gaussian bias: insights from discrete density peaks

Corrections induced by primordial non-Gaussianity to the linear halo bias can be computed from a peak-background split or the widespread local bias model. However, numerical simulations clearly support the prediction of the former, in which the non-Gaussian amplitude is proportional to the linear halo bias. To understand better the reasons behind the failure of standard Lagrangian local bias, in which the halo overdensity is a function of the local mass overdensity only, we explore the effect of a primordial bispectrum on the 2-point correlation of discrete density peaks. We show that the effective local bias expansion to peak clustering vastly simplifies the calculation. We generalize this approach to excursion set peaks and demonstrate that the resulting non-Gaussian amplitude, which is a weighted sum of quadratic bias factors, precisely agrees with the peak-background split expectation, which is a logarithmic derivative of the halo mass function with respect to the normalisation amplitude. We point out that statistics of thresholded regions can be computed using the same formalism. Our results suggest that halo clustering statistics can be modelled consistently (in the sense that the Gaussian and non-Gaussian bias factors agree with peak-background split expectations) from a Lagrangian bias relation only if the latter is specified as a set of constraints imposed on the linear density field. This is clearly not the case of standard Lagrangian local bias. Therefore, one is led to consider additional variables beyond the local mass overdensity.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Peak-Background Split, Renormalization, and Galaxy Clustering

We present a derivation of two-point correlations of general tracers in the peak-background split (PBS) framework by way of a rigorous definition of the PBS argument. Our expressions only depend on connected matter correlators and "renormalized" bias parameters with clear physical interpretation, and are independent of any coarse-graining scale. This result should be contrasted with the naive expression derived from a local bias expansion of the tracer number density with respect to the matter density perturbation δ_L coarse-grained on a scale R_L. In the latter case, the predicted tracer correlation function receives contributions of order <δ_L^n> at each perturbative order n, whereas, in our formalism, these are absorbed in the PBS bias parameters at all orders. Further, this approach naturally predicts both a scale-dependent bias ~ k^2 such as found for peaks of the density field, and the scale-dependent bias induced by primordial non-Gaussianity in the initial conditions. The only assumption made about the tracers is that their abundance at a given position depends solely on the matter distribution within a finite region around that position.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Testing the Running of non-Gaussianity through the CMB mu-distortion and the Halo Bias

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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The mass of the Milky Way and M31 using the method of least action

We constrain the most likely range of masses for the Milky Way and M31 using an application of the Numerical Action Method (NAM) that optimizes the fit to observed parameters over a large ensemble of NAM-generated solutions. Our 95% confidence level mass ranges, 1.5-4.5 x 10^12 m_sun for MW and 1.5-5.5 x 10^12 m_sun for M31, are consistent with the upper range of estimates from other methods and suggests that a larger proportion of the total mass becomes detectable when the peculiar motions of many nearby satellites are taken into account in the dynamical analysis. We test the method against simulated Local Group catalogs extracted from the Millennium Run to confirm that mass predictions are consistent with actual galaxy halo masses.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Distribution function approach to redshift space distortions, Part III: halos and galaxies

It was recently shown that the power spectrum in redshift space can be written as a sum of cross-power spectra between number weighted velocity moments. We investigate the properties of these power spectra for simulated galaxies and dark matter halos and compare them to the dark matter power spectra, generalizing the concept of the bias. Because all of the quantities are number weighted this approach is well defined even for sparse systems such as massive halos, in contrasts to the previous approaches to RSD where velocity correlations have been explored. We find that the number density weighting leads to a strong scale dependence of the bias terms for momentum density auto-correlation and cross-correlation with density. This trend becomes more significant for the more biased halos and leads to an enhancement of RSD power relative to the linear theory. Fingers-of-god effects, which in this formalism come from the correlations of the higher order moments beyond the momentum density, lead to smoothing of the power spectrum and can reduce this enhancement of power, but are relatively small for halos with no small-scale velocity dispersion. In comparison, for a more realistic galaxy sample with satellites the velocity dispersion generated by satellite motions inside the halos leads to a larger power suppression on small scales, but this depends on the satellite fraction. We investigate several statistics such as the two-dimensional power spectrum, its multipole moments, its powers of mu^2, and configuration space statistics. Overall we find that the nonlinear effects in realistic galaxy samples such as luminous red galaxies affect the redshift space clustering on very large scales: for example, the quadrupole moment is affected by 10% for k<0.1h/mpc, which means that these effects need to be understood if we want to extract cosmological information from the redshift space distortions.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Distribution function approach to redshift space distortions. Part II: N-body simulations

Measurement of redshift-space distortions (RSD) offers an attractive method to directly probe the cosmic growth history of density perturbations. A distribution function approach where RSD can be written as a sum over density weighted velocity moment correlators has recently been developed. We use Nbody simulations to investigate the individual contributions and convergence of this expansion for dark matter. If the series is expanded as a function of powers of mu, cosine of the angle between the Fourier mode and line of sight, there are a finite number of terms contributing at each order. We present these terms and investigate their contribution to the total as a function of wavevector k. For mu^2 the correlation between density and momentum dominates on large scales. Higher order corrections, which act as a Finger-of-God (FoG) term, contribute 1% at k~0.015h/Mpc, 10% at k~0.05h/Mpc at z=0, while for k>0.15h/Mpc they dominate and make the total negative. These higher order terms are dominated by density-energy density correlations which contribute negatively to the power, while the contribution from vorticity part of momentum density auto-correlation is an order of magnitude lower. For mu^4 term the dominant term on large scales is the scalar part of momentum density auto-correlation, while higher order terms dominate for k>0.15h/Mpc. For mu^6 and mu^8 we find it has very little power for k<0.15h/Mpc. We also compare the expansion to the full 2D P^ss(k,mu) as well as to their multipoles. For these statistics an infinite number of terms contribute and we find that the expansion achieves percent level accuracy for kmu<0.15h/Mpc at 6th order, but breaks down on smaller scales because the series is no longer perturbative. We explore resummation of the terms into FoG kernels, which extend the convergence up to a factor of 2 in scale. We find that the FoG kernels are approximately Lorentzian.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Evidence for Quadratic Tidal Tensor Bias from the Halo Bispectrum

The relation between the clustering properties of luminous matter in the form of galaxies and the underlying dark matter distribution is of fundamental importance for the interpretation of ongoing and upcoming galaxy surveys. The so called local bias model, where galaxy density is a function of local matter density, is frequently discussed as a means to infer the matter power spectrum or correlation function from the measured galaxy correlation. However, gravitational evolution generates a term quadratic in the tidal tensor and thus non-local in the density field, even if this term is absent in the initial conditions (Lagrangian space). Because the term is quadratic, it contributes as a loop correction to the power spectrum, so the standard linear bias picture still applies on large scales, however, it contributes at leading order to the bispectrum for which it is significant on all scales. Such a term could also be present in Lagrangian space if halo formation were influenced by the tidal field. We measure the corresponding coupling strengths from the matter-matter-halo bispectrum in numerical simulations and find a non-vanishing coefficient for the tidal tensor term. We find no scale dependence of the bias parameters up to k=0.1 h/Mpc and that the tidal effect is increasing with halo mass. While the Lagrangian bias picture is a better description of our results than the Eulerian bias picture, our results suggest that there might be a tidal tensor bias already in the initial conditions. We also find that the coefficients of the quadratic density term deviate quite strongly from the theoretical predictions based on the spherical collapse model and a universal mass function. Both quadratic density and tidal tensor bias terms must be included in the modeling of galaxy clustering of current and future surveys if one wants to achieve the high precision cosmology promise of these datasets.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Measuring primordial non-Gaussianity through weak lensing peak counts

We explore the possibility of detecting primordial non-Gaussianity of the local type using weak lensing peak counts. We measure the peak abundance in sets of simulated weak lensing maps corresponding to three models f_NL={0, +100, -100}. Using survey specifications similar to those of Euclid and without assuming any knowledge of the lens and source redshifts, we find the peak functions of the non-Gaussian models with f_NL=+-100 to differ by up to 15% from the Gaussian peak function at the high-mass end. For the assumed survey parameters, the probability of fitting an f_NL=0 peak function to the f_NL=+-100 peak functions is less than 0.1%. Assuming the other cosmological parameters known, f_NL can be measured with an error Δf_NL ~ 13. It is therefore possible that future weak lensing surveys like Euclid and LSST may detect primordial non-Gaussianity from the abundance of peak counts, and provide complementary information to that obtained from the cosmic microwave background.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Measuring primordial non-Gaussianity with weak-lensing surveys

We study the ability of future weak lensing (WL) surveys to constrain primordial non-Gaussianity of the local type. We use a large ensemble of simulated WL maps with survey specifications relevant to Euclid and LSST. The simulations assume Cold Dark Matter cosmologies that vary certain parameters around fiducial values: the non-Gaussianity parameter f_NL, the matter density parameter Omega_m, the amplitude of the matter power spectrum sigma_8, the spectral index of the primordial power spectrum n_s, and the dark-energy equation-of-state parameter w_0. We assess the sensitivity of the cosmic shear correlation functions, the third-order aperture mass statistics, and the abundance of shear peaks to these parameters. We find that each of the considered probes provides unmarginalized constraints of Delta f_NL ~ 20 on f_NL. Marginalized constraints from any individual WL probe are much weaker due to strong correlations between parameters. However, the parameter errors can be substantially reduced by combining information from different WL probes. Combining all WL probes yields the following marginal (68% confidence level) uncertainties: Delta f_NL ~ 50, Delta Omega_m ~ 0.002, Delta sigma_8 ~ 0.004, Delta n_s ~ 0.007, and Delta w_0 ~ 0.03. We examine the bias induced by neglecting f_NL on the constraints on the other parameters. We find sigma_8 and w_0 to be the most affected. Moreover, neglecting non-Gaussianity leads to a severe underestimation of the uncertainties in the other cosmological parameters.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Optimal Constraints on Local Primordial Non-Gaussianity from the Two-Point Statistics of Large-Scale Structure

One of the main signatures of primordial non-Gaussianity of the local type is a scale-dependent correction to the bias of large-scale structure tracers such as galaxies or clusters, whose amplitude depends on the bias of the tracers itself. The dominant source of noise in the power spectrum of the tracers is caused by sampling variance on large scales (where the non-Gaussian signal is strongest) and shot noise arising from their discrete nature. Recent work has argued that one can avoid sampling variance by comparing multiple tracers of different bias, and suppress shot noise by optimally weighting halos of different mass. Here we combine these ideas and investigate how well the signatures of non-Gaussian fluctuations in the primordial potential can be extracted from the two-point correlations of halos and dark matter. On the basis of large $N$-body simulations with local non-Gaussian initial conditions and their halo catalogs we perform a Fisher matrix analysis of the two-point statistics. Compared to the standard analysis, optimal weighting- and multiple-tracer techniques applied to halos can yield up to one order of magnitude improvements in $\fnl$-constraints, even if the underlying dark matter density field is not known. We compare our numerical results to the halo model and find satisfactory agreement. Forecasting the optimal $\fnl$-constraints that can be achieved with our methods when applied to existing and future survey data, we find that a survey of $50h^{-1}\mathrm{Gpc}^3$ volume resolving all halos down to $10^{11}\hMsun$ at $z=1$ will be able to obtain $σ_{\fnl}\sim1$ (68% cl), a factor of $\sim20$ improvement over the current limits. Decreasing the minimum mass of resolved halos, increasing the survey volume or obtaining the dark matter maps can further improve these limits, potentially reaching the level of $σ_{\fnl}\sim0.1$. (abridged)

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Optimal Weighting in Galaxy Surveys: Application to Redshift-Space Distortions

Using multiple tracers of large-scale structure allows to evade the limitations imposed by sampling variance for some parameters of interest in cosmology. We demonstrate the optimal way of carrying out a multitracer analysis in a galaxy redshift survey by considering the principal components of the shot noise matrix from two-point clustering statistics. We show how to construct two tracers that maximize the benefits of sampling variance and shot noise cancellation using optimal weights. On the basis of high-resolution N-body simulations of dark matter halos we apply this technique to the analysis of redshift-space distortions and demonstrate how constraints on the growth rate of structure formation can be substantially improved. The primary limitation are nonlinear effects, which cause significant biases in the method already at scales of k<0.1h/Mpc, suggesting the need to develop nonlinear models of redshift-space distortions in order to extract the maximum information from future redshift surveys. Nonetheless we find gains of a factor of a few in constraints on the growth rate achievable when merely the linear regime of a galaxy survey like EUCLID is considered.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Testing Multi-Field Inflation with Galaxy Bias

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

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Accurate Predictions for the Scale-Dependent Galaxy Bias from Primordial Non-Gaussianity

The large-scale clustering of galaxies can serve as a probe of primordial non-Gaussianity in the Universe competitive with the anisotropies of the CMB. Here, we present improved theoretical predictions which include an important, previously overlooked correction to the bias. We demonstrate that the new predictions are able to reproduce the results of N-body simulations, explaining the significant departures seen from previous theoretical results. These refined predictions open the way to accurate constraints on primordial physics with large-scale structure surveys.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Non-Gaussian Halo Bias Re-examined: Mass-dependent Amplitude from the Peak-Background Split and Thresholding

Recent results of N-body simulations have shown that current theoretical models are not able to correctly predict the amplitude of the scale-dependent halo bias induced by primordial non-Gaussianity, for models going beyond the simplest, local quadratic case. Motivated by these discrepancies, we carefully examine three theoretical approaches based on (1) the statistics of thresholded regions, (2) a peak-background split method based on separation of scales, and (3) a peak-background split method using the conditional mass function. We first demonstrate that the statistics of thresholded regions, which is shown to be equivalent at leading order to a local bias expansion, cannot explain the mass-dependent deviation between theory and N-body simulations. In the two formulations of the peak-background split on the other hand, we identify an important, but previously overlooked, correction to the non-Gaussian bias that strongly depends on halo mass. This new term is in general significant for any primordial non-Gaussianity going beyond the simplest local fNL model. In a separate paper, we compare these new theoretical predictions with N-body simulations, showing good agreement for all simulated types of non-Gaussianity.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Nonlinear clustering in models with primordial non-Gaussianity: the halo model approach

We develop the halo model of large-scale structure as an accurate tool for probing primordial non-Gaussianity. In this study we focus on understanding the matter clustering at several redshifts. The primordial non-Gaussianity is modeled as a quadratic correction to the local Gaussian potential, and is characterized by the parameter f_NL. In our formulation of the halo model we pay special attention to the effect of halo exclusion, and show that this can potentially solve the long standing problem of excess power on large scales in this model. The model depends on the mass function, clustering and density profiles of halos. We test these ingredients using a large ensemble of high-resolution Gaussian and non-Gaussian numerical simulations. In particular, we provide a first exploration of how density profiles change in the presence of primordial non-Gaussianities. We find that for f_NL positive/negative high mass halos have an increased/decreased core density, so being more/less concentrated than in the Gaussian case. We also examine the halo bias and show that, if the halo model is correct, then there is a small asymmetry in the scale-dependence of the bias on very large scales, which arises because the Gaussian bias must be renormalized. We show that the matter power spectrum is modified by ~2.5% and ~3.5% on scales k~1.0 h/Mpc at z=0 and z=1, respectively. Our halo model calculation reproduces the absolute amplitude to within 10% and the ratio of non-Gaussian to Gaussian spectra to within 1%. We also measure the matter correlation functions and find similarly good agreement between the model and the data. We anticipate that this modeling will be useful for constraining f_NL from measurements of the shear correlation function in future weak lensing surveys such as Euclid.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The Halo Bispectrum in N-body Simulations with non-Gaussian Initial Conditions

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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Minimizing the stochasticity of halos in large-scale structure surveys

In recent work (Seljak, Hamaus and Desjacques 2009) it was found that weighting central halo galaxies by halo mass can significantly suppress their stochasticity relative to the dark matter, well below the Poisson model expectation. In this paper we extend this study with the goal of finding the optimal mass-dependent halo weighting and use $N$-body simulations to perform a general analysis of halo stochasticity and its dependence on halo mass. We investigate the stochasticity matrix, defined as $C_{ij}\equiv<(δ_i -b_iδ_m)(δ_j-b_jδ_m)>$, where $δ_m$ is the dark matter overdensity in Fourier space, $δ_i$ the halo overdensity of the $i$-th halo mass bin and $b_i$ the halo bias. In contrast to the Poisson model predictions we detect nonvanishing correlations between different mass bins. We also find the diagonal terms to be sub-Poissonian for the highest-mass halos. The diagonalization of this matrix results in one large and one low eigenvalue, with the remaining eigenvalues close to the Poisson prediction $1/\bar{n}$, where $\bar{n}$ is the mean halo number density. The eigenmode with the lowest eigenvalue contains most of the information and the corresponding eigenvector provides an optimal weighting function to minimize the stochasticity between halos and dark matter. We find this optimal weighting function to match linear mass weighting at high masses, while at the low-mass end the weights approach a constant whose value depends on the low-mass cut in the halo mass function. Finally, we employ the halo model to derive the stochasticity matrix and the scale-dependent bias from an analytical perspective. It is remarkably successful in reproducing our numerical results and predicts that the stochasticity between halos and the dark matter can be reduced further when going to halo masses lower than we can resolve in current simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Modeling scale-dependent bias on the baryonic acoustic scale with the statistics of peaks of Gaussian random fields

Models of galaxy and halo clustering commonly assume that the tracers can be treated as a continuous field locally biased with respect to the underlying mass distribution. In the peak model pioneered by BBKS, one considers instead density maxima of the initial, Gaussian mass density field as an approximation to the formation site of virialized objects. In this paper, the peak model is extended in two ways to improve its predictive accuracy. Firstly, we derive the two-point correlation function of initial density peaks up to second order and demonstrate that a peak-background split approach can be applied to obtain the k-independent and k-dependent peak bias factors at all orders. Secondly, we explore the gravitational evolution of the peak correlation function within the Zel'dovich approximation. We show that the local (Lagrangian) bias approach emerges as a special case of the peak model, in which all bias parameters are scale-independent and there is no statistical velocity bias. We apply our formulae to study how the Lagrangian peak biasing, the diffusion due to large scale flows and the mode-coupling due to nonlocal interactions affect the scale dependence of bias from small separations up to the baryon acoustic oscillation (BAO) scale. For 2-sigma density peaks collapsing at z=0.3, our model predicts a ~ 5% residual scale-dependent bias around the acoustic scale that arises mostly from first-order Lagrangian peak biasing (as opposed to second-order gravity mode-coupling). We also search for a scale dependence of bias in the large scale auto-correlation of massive halos extracted from a very large N-body simulation provided by the MICE collaboration. For halos with mass M>10^{14}Msun/h, our measurements demonstrate a scale-dependent bias across the BAO feature which is very well reproduced by a prediction based on the peak model.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Primordial non-Gaussianity from the large scale structure

Primordial non-Gaussianity is a potentially powerful discriminant of the physical mechanisms that generated the cosmological fluctuations observed today. Any detection of non-Gaussianity would have profound implications for our understanding of cosmic structure formation. In this paper, we review past and current efforts in the search for primordial non-Gaussianity in the large scale structure of the Universe.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Primordial non-Gaussianity in the large scale structure of the Universe

Primordial non-Gaussianity is a potentially powerful discriminant of the physical mechanisms that generated the cosmological fluctuations observed today. Any detection of significant non-Gaussianity would thus have profound implications for our understanding of cosmic structure formation. The large scale mass distribution in the Universe is a sensitive probe of the nature of initial conditions. Recent theoretical progress together with rapid developments in observational techniques will enable us to critically confront predictions of inflationary scenarios and set constraints as competitive as those from the Cosmic Microwave Background. In this paper, we review past and current efforts in the search for primordial non-Gaussianity in the large scale structure of the Universe.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Redshift space correlations and scale-dependent stochastic biasing of density peaks

We calculate the redshift space correlation function and the power spectrum of density peaks of a Gaussian random field. In the linear regime k < 0.1 h/Mpc, the redshift space power spectrum is P^s_{pk}(k,u) = exp(-f^2 s_{vel}^2 k^2 u^2) * [b_{pk}(k) + b_{vel}(k) f u^2]^2 * P_m(k), where u is the angle with respect to the line of sight, s_{vel} is the one-dimensional velocity dispersion, f is the growth rate, and b_{pk}(k) and b_{vel}(k) are k-dependent linear spatial and velocity bias factors. For peaks, the value of s_{vel} depends upon the functional form of b_{vel}. The peaks model is remarkable because it has unbiased velocities -- peak motions are driven by dark matter flows -- but, in order to achieve this, b_{vel} is k-dependent. We speculate that this is true in general: k-dependence of the spatial bias will lead to k-dependence of b_{vel} even if the biased tracers flow with the dark matter. Because of the k-dependence of the linear bias parameters, standard manipulations applied to the peak model will lead to k-dependent estimates of the growth factor that could erroneously be interpreted as a signature of modified dark energy or gravity. We use the Fisher formalism to show that the constraint on the growth rate f is degraded by a factor of two if one allows for a k-dependent velocity bias of the peak type. We discuss a simple estimate of nonlinear evolution and illustrate the effect of the peak bias on the redshift space multipoles. For k < 0.1 h/Mpc, the peak bias is deterministic but k-dependent, so the configuration space bias is stochastic and scale dependent, both in real and redshift space. We provide expressions for this stochasticity and its evolution (abridged).

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The large-scale clustering of massive dark matter haloes

The statistics of peaks of the initial, Gaussian density field can be used to interpret the abundance and clustering of massive dark matter haloes. I discuss some recent theoretical results related to their clustering and its redshift evolution. Predictions from the peak model are qualitatively consistent with measurements of the linear bias of high mass haloes, which also show some evidence for a dependence on the halo mass M at fixed peak height. The peak approach also predicts a distinctive scale-dependence in the bias of haloes across the baryon acoustic feature, a measurement of which would provide strong support for its validity. For 2-sigma density peaks collapsing at z=0.3, this residual scale-dependent bias is at the 5-10 percent level and should thus be within reach of very large simulations of structure formation.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The Matter Bispectrum in N-body Simulations with non-Gaussian Initial Conditions

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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The nonlinear redshift space probability distribution function in models with local primordial non-Gaussianity

We use the ellipsoidal collapse approximation to investigate the nonlinear redshift space evolution of the density field with primordial non-Gaussianity of the local f_{nl}-type. We utilize the joint distribution of eigenvalues of the initial non-Gaussian shear field and evaluate the evolved redshift space probability distribution function (PDF). It is shown that, similar to the real space analysis, the underdense tail of the nonlinear redshift space PDF differs significantly from that for Gaussian initial conditions. We also derive the lowest order correction of the Kaiser's formulain the presence of a non-zero f_{nl}.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Scale-dependent bias induced by local non-Gaussianity: A comparison to N-body simulations

We investigate the effect of primordial non-Gaussianity of the local f_NL type on the auto- and cross-power spectrum of dark matter haloes using simulations of the LCDM cosmology. We perform a series of large N-body simulations of both positive and negative f_NL, spanning the range between 10 and 100. Theoretical models predict a scale-dependent bias correction Δb(k,f_NL) that depends on the linear halo bias b(M). We measure the power spectra for a range of halo mass and redshifts covering the relevant range of existing galaxy and quasar populations. We show that auto and cross-correlation analyses of bias are consistent with each other. We find that for low wavenumbers with k<0.03 h/Mpc the theory and the simulations agree well with each other for biased haloes with b(M)>1.5. We show that a scale-independent bias correction improves the comparison between theory and simulations on smaller scales, where the scale-dependent effect rapidly becomes negligible. The current limits on f_NL from Slosar et al. (2008) come mostly from very large scales k<0.01 h/Mpc and, therefore, remain valid. For the halo samples with b(M)<1.5-2 we find that the scale- dependent bias from non-Gaussianity actually exceeds the theoretical predictions. Our results are consistent with the bias correction scaling linearly with f_NL.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Signature of primordial non-Gaussianity of phi^3-type in the mass function and bias of dark matter haloes

We explore the effect of a cubic correction gnl*phi^3 on the mass function and bias of dark matter haloes extracted from a series of large N-body simulations and compare it to theoretical predictions. Such cubic terms can be motivated in scenarios like the curvaton model, in which a large cubic correction can be produced while simultaneously keeping the quadratic fnl*phi^2 correction small. The deviation from the Gaussian halo mass function is in reasonable agreement with the theoretical predictions. The scale-dependent bias correction Delta b_kappa(k,gnl) measured from the auto- and cross-power spectrum of haloes, is similar to the correction in fnl models, but the amplitude is lower than theoretical expectations. Using the compilation of LSS data in Slosar et al. (2008), we obtain for the first time a limit on gnl of -3.5*10^5 < gnl < +8.2*10^5 (at 95% CL). This limit will improve with the future LSS data by 1-2 orders of magnitude, which should test many of the scenarios of this type.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

The initial shear field in models with primordial local non-Gaussianity and implications for halo and void abundances

We generalize Doroshkevich's celebrated formulae for the eigenvalues of the initial shear field associated with Gaussian statistics to the local non-Gaussian f_{nl} model. This is possible because, to at least second order in f_{nl}, distributions at fixed overdensity are unchanged from the case f_{nl}=0. We use this generalization to estimate the effect of f_{nl}\ne 0 on the abundance of virialized halos. Halo abundances are expected to be related to the probability that a certain quantity in the initial fluctuation field exceeds a threshold value, and we study two choices for this variable: it can either be the sum of the eigenvalues of the initial deformation tensor (the initial overdensity), or its smallest eigenvalue. The approach based on a critical overdensity yields results which are in excellent agreement with numerical measurements. We then use these same methods to develop approximations describing the sensitivity of void abundances on f_{nl}. While a positive f_{nl} produces more extremely massive halos, it makes fewer extremely large voids. Its effect thus is qualitatively different from a simple rescaling of the normalisation of the density fluctuation field σ_8. Therefore, void abundances furnish complementary information to cluster abundances, and a joint comparison of both might provide interesting constraints on primordial non-Gaussianity.

0 citations0 reviewsNo code linkedOpen access