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

Roman Scoccimarro

Roman Scoccimarro 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-thastro-phgr-qcastro-ph.GA

Catalog footprint

What is connected

27works
5topics
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

27 published item(s)

preprint2023arXiv

Squeezing $f_{\rm NL}$ out of the matter bispectrum with consistency relations

We show how consistency relations can be used to robustly extract the amplitude of local primordial non-Gaussianity ($f_{\rm NL}$) from the squeezed limit of the matter bispectrum, well into the non-linear regime. First, we derive a non-perturbative relation between primordial non-Gaussianity and the leading term in the squeezed bispectrum, revising some results present in the literature. This relation is then used to successfully measure $f_{\rm NL}$ from $N$-body simulations. We discuss the dependence of our results on different scale cuts and redshifts. Specifically, the analysis is strongly dependent on the choice of the smallest soft momentum, $q_{\rm min}$, which is the most sensitive to primordial bispectrum contributions, but is largely independent of the choice of the largest hard momentum, $k_{\rm max}$, due to the non-Gaussian nature of the covariance. We also show how the constraints on $f_{\rm NL}$ improve at higher redshift, due to a reduced off-diagonal covariance. In particular, for a simulation with $f_{\rm NL} = 100$ and a volume of $(2.4 \text{ Gpc}/h)^3$, we measure $f_{\rm NL} = 98 \pm 12$ at redshift $z=0$ and $f_{\rm NL} = 97 \pm 8$ at $z=0.97$. Finally, we compare our results with a Fisher forecast, showing that the current version of the analysis is satisfactorily close to the Fisher error. We regard this as a first step towards the realistic application of consistency relations to constrain primordial non-Gaussianity using observations.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Cosmological implications of the full shape of anisotropic clustering measurements in BOSS and eBOSS

We present the analysis of the full shape of anisotropic clustering measurement from the extended Baryon Oscillation Spectroscopic Survey (eBOSS) quasar sample together with the combined galaxy sample from the Baryon Oscillation Spectroscopic Survey (BOSS), re-analysed using an updated recipe for the non-linear matter power spectrum and the non-local bias parameters. We obtain constraints for flat $Λ$CDM cosmologies, focusing on the cosmological parameters that are independent of the Hubble parameter $h$. Our recovered value for the RMS linear perturbation theory variance as measured on the scale of $12\,{\rm Mpc}$ is $σ_{12}=0.805\pm 0.049$, while using the traditional reference scale of $8\,h^{-1}{\rm Mpc}$ gives $σ_{8}=0.815\pm 0.044$. We quantify the agreement between our measurements and the latest CMB data from Planck using the suspiciousness metric, and find them to be consistent within $0.64 \pm 0.03σ$. Combining our clustering constraints with the $3\times2$pt data sample from the Dark Energy Survey (DES) Year 1 release slightly degrades this agreement to the level of $1.54 \pm 0.08σ$, while still showing an overall consistency with Planck. We furthermore study the effect of imposing a Planck - like prior on the parameters that define the shape of the linear matter power spectrum, and find significantly tighter constraints on the parameters that control the evolution of density fluctuations. In particular, the combination of low-redshift data sets prefers a value of the physical dark energy density $ω_{\rm DE}=0.335 \pm 0.011$, which is 1.7$σ$ higher than the one preferred by Planck.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Constraining $M_ν$ with the Bispectrum I: Breaking Parameter Degeneracies

Massive neutrinos suppress the growth of structure below their free-streaming scale and leave an imprint on large-scale structure. Measuring this imprint allows us to constrain the sum of neutrino masses, $M_ν$, a key parameter in particle physics beyond the Standard Model. However, degeneracies among cosmological parameters, especially between $M_ν$ and $σ_8$, limit the constraining power of standard two-point clustering statistics. In this work, we investigate whether we can break these degeneracies and constrain $\smnu$ with the next higher-order correlation function --- the bispectrum. We first examine the redshift-space halo bispectrum of $800$ $N$-body simulations from the HADES suite and demonstrate that the bispectrum helps break the $M_ν$--$σ_8$ degeneracy. Then using 22,000 $N$-body simulations of the Quijote suite, we quantify for the first time the full information content of the redshift-space halo bispectrum down to nonlinear scales using a Fisher matrix forecast of $\{Ω_m$, $Ω_b$, $h$, $n_s$, $σ_8$, $M_ν\}$. For $k_{\rm max}{=}0.5~h/{\rm Mpc}$, the bispectrum provides $Ω_m$, $Ω_b$, $h$, $n_s$, and $σ_8$ constraints 1.9, 2.6, 3.1, 3.6, and 2.6 times tighter than the power spectrum. For $M_ν$, the bispectrum improves the 1$σ$ constraint from 0.2968 to 0.0572 eV --- over 5 times tighter than the power spectrum. Even with priors from {\em Planck}, the bispectrum improves $M_ν$ constraints by a factor of 1.8. Although we reserve marginalizing over a more complete set of bias parameters to the next paper of the series, these constraints are derived for a $(1~h^{-1}{\rm Gpc})^3$ box, a substantially smaller volume than upcoming surveys. Thus, our results demonstrate that the bispectrum offers significant improvements over the power spectrum, especially for constraining $M_ν$.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Matter power spectrum and the challenge of percent accuracy

Future galaxy surveys require one percent precision in the theoretical knowledge of the power spectrum over a large range including very nonlinear scales. While this level of accuracy is easily obtained in the linear regime with perturbation theory, it represents a serious challenge for small scales where numerical simulations are required. In this paper we quantify the precision of present-day $N$-body methods, identifying main potential error sources from the set-up of initial conditions to the measurement of the final power spectrum. We directly compare three widely used $N$-body codes, Ramses, Pkdgrav3, and Gadget3 which represent three main discretisation techniques: the particle-mesh method, the tree method, and a hybrid combination of the two. For standard run parameters, the codes agree to within one percent at $k\leq1$ $h\,\rm Mpc^{-1}$ and to within three percent at $k\leq10$ $h\,\rm Mpc^{-1}$. We also consider the bispectrum and show that the reduced bispectra agree at the sub-percent level for $k\leq 2$ $h\,\rm Mpc^{-1}$. In a second step, we quantify potential errors due to initial conditions, box size, and resolution using an extended suite of simulations performed with our fastest code Pkdgrav3. We demonstrate that the simulation box size should not be smaller than $L=0.5$ $h^{-1}\rm Gpc$ to avoid systematic finite-volume effects (while much larger boxes are required to beat down the statistical sample variance). Furthermore, a maximum particle mass of $M_{\rm p}=10^{9}$ $h^{-1}\rm M_{\odot}$ is required to conservatively obtain one percent precision of the matter power spectrum. As a consequence, numerical simulations covering large survey volumes of upcoming missions such as DES, LSST, and Euclid will need more than a trillion particles to reproduce clustering properties at the targeted accuracy.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The clustering of galaxies in the completed SDSS-III Baryon Oscillation Spectroscopic Survey: cosmological implications of the configuration-space clustering wedges

We explore the cosmological implications of anisotropic clustering measurements in configuration space of the final galaxy samples from Data Release 12 of the SDSS-III Baryon Oscillation Spectroscopic Survey. We implement a new detailed modelling of the effects of non-linearities, galaxy bias and redshift-space distortions that can be used to extract unbiased cosmological information from our measurements for scales $s \gtrsim 20\,h^{-1}{\rm Mpc}$. We combined the galaxy clustering information from BOSS with the latest cosmic microwave background (CMB) observations and Type Ia supernovae samples and found no significant evidence for a deviation from the $Λ$CDM cosmological model. In particular, these data sets can constrain the dark energy equation of state parameter to $w_{\rm DE}=-0.996\pm0.042$ when assumed time-independent, the curvature of the Universe to $Ω_{k}=-0.0007\pm 0.0030$ and the sum of the neutrino masses to $\sum m_ν < 0.25\,{\rm eV}$ at 95 per cent CL. We explore the constraints on the growth rate of cosmic structures assuming $f(z)=Ω_{\rm m}(z)^γ$ and obtain $γ= 0.609\pm 0.079$, in good agreement with the predictions of general relativity of $γ=0.55$. We compress the information of our clustering measurements into constraints on the parameter combinations $D_{\rm V}(z)/r_{\rm d}$, $F_{\rm AP}(z)$ and $fσ_8(z)$ at the effective redshifts of $z=0.38$, $0.51$ and $0.61$ with their respective covariance matrices and find good agreement with the predictions for these parameters obtained from the best-fitting $Λ$CDM model to the CMB data from the Planck satellite. This paper is part of a set that analyses the final galaxy clustering dataset from BOSS. The measurements and likelihoods presented here are combined with others in Alam et al. (2016) to produce the final cosmological constraints from BOSS.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Fast Estimators for Redshift-Space Clustering

Redshift-space distortions in galaxy surveys happen along the radial direction, breaking statistical translation invariance. We construct estimators for radial distortions that, using only Fast Fourier Transforms (FFTs) of the overdensity field multipoles for a given survey geometry, compute the power spectrum monopole, quadrupole and hexadecapole, and generalize such estimators to the bispectrum. Using realistic mock catalogs we compare the signal to noise of two estimators for the power spectrum hexadecapole that require different number of FFTs and measure the bispectrum monopole, quadrupole and hexadecapole. The resulting algorithm is very efficient, e.g. for the BOSS survey requires about three minutes for $\ell=0,2,4$ power spectra for scales up to $k=0.3~h/$Mpc and about fifteen additional minutes for $\ell=0,2,4$ bispectra for all scales and triangle shapes up to $k=0.2~h/$Mpc on a single core. The speed of these estimators is essential as it makes possible to compute covariance matrices from large number of realizations of mock catalogs with realistic survey characteristics, and paves the way for improved constrains of gravity on cosmological scales, inflation and galaxy bias.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The Spatial Distribution of Satellite Galaxies Within Halos: Measuring the Very Small Scale Angular Clustering of SDSS Galaxies

We measure the angular clustering of galaxies from the Sloan Digital Sky Survey Data Release 7 in order to probe the spatial distribution of satellite galaxies within their dark matter halos. Specifically, we measure the angular correlation function on very small scales (7-320") in a range of luminosity threshold samples (absolute r-band magnitudes of -18 up to -21) that are constructed from the subset of SDSS that has been spectroscopically observed more than once (the so-called plate overlap region). We choose to measure angular clustering in this reduced survey footprint in order to minimize the effects of fiber collision incompleteness, which are otherwise substantial on these small scales. We model our clustering measurements using a fully numerical halo model that populates dark matter halos in N-body simulations to create realistic mock galaxy catalogs. The model has free parameters that specify both the number and spatial distribution of galaxies within their host halos. We adopt a flexible density profile for the spatial distribution of satellite galaxies that is similar to the dark matter Navarro-Frenk-White (NFW) profile, except that the inner slope is allowed to vary. We find that the angular clustering of our most luminous samples (Mr< -20 and -21) suggests that luminous satellite galaxies have substantially steeper inner density profiles than NFW. Lower luminosity samples are less constraining, however, and are consistent with satellite galaxies having shallow density profiles. Our results confirm the findings of Watson et al. 2012 while using different clustering measurements and modeling methodology.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Constraining Primordial Non-Gaussianity with Moments of the Large Scale Density Field

We use cosmological N-body simulations to investigate whether measurements of the moments of large-scale structure can yield constraints on primordial non-Gaussianity. We measure the variance, skewness, and kurtosis of the evolved density field from simulations with Gaussian and three different non-Gaussian initial conditions: a local model with f_NL=100, an equilateral model with f_NL=-400, and an orthogonal model with f_NL=-400. We show that the moments of the dark matter density field differ significantly between Gaussian and non-Gaussian models. We also make the measurements on mock galaxy catalogs that contain galaxies with clustering properties similar to those of luminous red galaxies (LRGs). We find that, in the case of skewness and kurtosis, galaxy bias reduces the detectability of non-Gaussianity, though we can still clearly discriminate between different models in our simulation volume. However, in the case of the variance, galaxy bias greatly amplifies the detectability of non-Gaussianity. In all cases we find that redshift distortions do not significantly affect the detectability. When we restrict our measurements to volumes equivalent to the Sloan Digital Sky Survey II (SDSS-II) or Baryon Oscillation Spectroscopic Survey (BOSS) samples, the probability of detecting a departure from the Gaussian model is high by using measurements of the variance, but very low by using only skewness and kurtosis measurements. We find that skewness and kurtosis measurements are never likely to yield useful constraints on primordial non-Gaussianity, but future surveys should be large enough to place meaningful constraints using measurements of the galaxy variance.(Abridged)

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Non-local Lagrangian bias

Halos are biased tracers of the dark matter distribution. It is often assumed that the patches from which halos formed are locally biased with respect to the initial fluctuation field, meaning that the halo-patch fluctuation field can be written as a Taylor series in that of the dark matter. If quantities other than the local density influence halo formation, then this Lagrangian bias will generically be nonlocal; the Taylor series must be performed with respect to these other variables as well. We illustrate the effect with Monte-Carlo simulations of a model in which halo formation depends on the local shear (the quadrupole of perturbation theory), and provide an analytic model which provides a good description of our results. Our model, which extends the excursion set approach to walks in more than one dimension, works both when steps in the walk are uncorrelated, as well as when there are correlations between steps. For walks with correlated steps, our model includes two distinct types of nonlocality: one is due to the fact that the initial density profile around a patch which is destined to form a halo must fall sufficiently steeply around it -- this introduces k-dependence to even the linear bias factor, but otherwise only affects the monopole of the clustering signal. The other is due to the surrounding shear field; this affects the quadratic and higher order bias factors, and introduces an angular dependence to the clustering signal. In both cases, our analysis shows that these nonlocal Lagrangian bias terms can be significant, particularly for massive halos; they must be accounted for in analyses of higher order clustering such as the halo bispectrum in Lagrangian or Eulerian space. Although we illustrate these effects using halos, our analysis and conclusions also apply to the other constituents of the cosmic web -- filaments, sheets and voids.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Novel Probes of Gravity and Dark Energy

The discovery of cosmic acceleration has stimulated theorists to consider dark energy or modifications to Einstein's General Relativity as possible explanations. The last decade has seen advances in theories that go beyond smooth dark energy -- modified gravity and interactions of dark energy. While the theoretical terrain is being actively explored, the generic presence of fifth forces and dark sector couplings suggests a set of distinct observational signatures. This report focuses on observations that differ from the conventional probes that map the expansion history or large-scale structure. Examples of such novel probes are: detection of scalar fields via lab experiments, tests of modified gravity using stars and galaxies in the nearby universe, comparison of lensing and dynamical masses of galaxies and clusters, and the measurements of fundamental constants at high redshift. The observational expertise involved is very broad as it spans laboratory experiments, high resolution astronomical imaging and spectroscopy and radio observations. In the coming decade, searches for these effects have the potential for discovering fundamental new physics. We discuss how the searches can be carried out using experiments that are already under way or with modest adaptations of existing telescopes or planned experiments. The accompanying paper on the Growth of Cosmic Structure describes complementary tests of gravity with observations of large-scale structure.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Constructing Regularized Cosmic Propagators

We present a new scheme for the general computation of cosmic propagators that allow to interpolate between standard perturbative results at low-k and their expected large-k resummed behavior. This scheme is applicable to any multi-point propagator and allows the matching of perturbative low-k calculations to any number of loops to their large-k behavior, and can potentially be applied in case of non-standard cosmological scenarios such as those with non-Gaussian initial conditions. The validity of our proposal is checked against previous prescriptions and measurements in numerical simulations showing a remarkably good agreement. Such a generic prescription for multi-point propagators provides the necessary building blocks for the computation of polyspectra in the context of the so-called Gamma-expansion introduced by Bernardeau et al. (2008). As a concrete application we present a consistent calculation of the matter bispectrum at one-loop order.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Gravity and Large-Scale Non-local Bias

The relationship between galaxy and matter overdensities, bias, is most often assumed to be local. This is however unstable under time evolution, we provide proofs under several sets of assumptions. In the simplest model galaxies are created locally and linearly biased at a single time, and subsequently move with the matter (no velocity bias) conserving their comoving number density (no merging). We show that, after this formation time, the bias becomes unavoidably non-local and non-linear at large scales. We identify the non-local gravitationally induced fields in which the galaxy overdensity can be expanded, showing that they can be constructed out of the invariants of the deformation tensor (Galileons). In addition, we show that this result persists if we include an arbitrary evolution of the comoving number density of tracers. We then include velocity bias, and show that new contributions appear, a dipole field being the signature at second order. We test these predictions by studying the dependence of halo overdensities in cells of fixed matter density: measurements in simulations show that departures from the mean bias relation are strongly correlated with the non-local gravitationally induced fields identified by our formalism. The effects on non-local bias seen in the simulations are most important for the most biased halos, as expected from our predictions. The non-locality seen in the simulations is not fully captured by assuming local bias in Lagrangian space. Accounting for these effects when modeling galaxy bias is essential for correctly describing the dependence on triangle shape of the galaxy bispectrum, and hence constraining cosmological parameters and primordial non-Gaussianity. We show that using our formalism we remove an important systematic in the determination of bias parameters from the galaxy bispectrum, particularly for luminous galaxies. (abridged)

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Halo Sampling, Local Bias and Loop Corrections

We develop a new test of local bias, by constructing a locally biased halo density field from sampling the dark matter-halo distribution. Our test differs from conventional tests in that it preserves the full scatter in the bias relation and it does not rely on perturbation theory. We put forward that bias parameters obtained using a smoothing scale R can only be applied to computing the halo power spectrum at scales k ~ 1/R. Our calculations can automatically include the running of bias parameters and give vanishingly small loop corrections at low-k. Our proposal results in much better agreement of the sampling and perturbation theory results with simulations. In particular, unlike the standard interpretation of local bias in the literature, our treatment of local bias does not generate a constant power in the low-k limit. We search for extra noise in the Poisson corrected halo power spectrum at wavenumbers below its turn-over and find no evidence of significant positive noise (as predicted by the standard interpretation) while we find evidence of negative noise coming from halo exclusion for very massive halos. Using perturbation theory and our non-perturbative sampling technique we also demonstrate that nonlocal bias effects discovered recently in simulations impact the power spectrum only at the few percent level in the weakly nonlinear regime.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Large-scale Bias and Efficient Generation of Initial Conditions for Non-Local Primordial Non-Gaussianity

We study the scale-dependence of halo bias in generic (non-local) primordial non-Gaussian (PNG) initial conditions of the type motivated by inflation, parametrized by an arbitrary quadratic kernel. We first show how to generate non-local PNG initial conditions with minimal overhead compared to local PNG models for a general class of primordial bispectra that can be written as linear combinations of separable templates. We run cosmological simulations for the local, and non-local equilateral and orthogonal models and present results on the scale-dependence of halo bias. We also derive a general formula for the Fourier-space bias using the peak-background split (PBS) in the context of the excursion set approach to halos and discuss the difference and similarities with the known corresponding result from local bias models. Our PBS bias formula generalizes previous results in the literature to include non-Markovian effects and non-universality of the mass function and are in better agreement with measurements in numerical simulations than previous results for a variety of halo masses, redshifts and halo definitions. We also derive for the first time quadratic bias results for arbitrary non-local PNG, and show that non-linear bias loops give small corrections at large-scales. The resulting well-behaved perturbation theory paves the way to constrain non-local PNG from measurements of the power spectrum and bispectrum in galaxy redshift surveys.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

MPTbreeze: A fast renormalized perturbative scheme

We put forward and test a simple description of multi-point propagators (MP), which serve as building-blocks to calculate the nonlinear matter power spectrum. On large scales these propagators reduce to the well-known kernels in standard perturbation theory, while at smaller scales they are suppresed due to nonlinear couplings. Through extensive testing with numerical simulations we find that this decay is characterized by the same damping scale for both two and three-point propagators. In turn this transition can be well modeled with resummation results that exponentiate one-loop computations. For the first time, we measure the four components of the non-linear (two-point) propagator using dedicated simulations started from two independent random Gaussian fields for positions and velocities, verifying in detail the fundamentals of propagator resummation. We use these results to develop an implementation of the MP-expansion for the nonlinear power spectrum that only requires seconds to evaluate at BAO scales. To test it we construct six suites of large numerical simulations with different cosmologies. From these and LasDamas runs we show that the nonlinear power spectrum can be described at the ~ 2% level at BAO scales for redshifts in the range [0-2.5]. We make a public release of the MPTbreeze code with the hope that it can be useful to the community.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The clustering of galaxies in the SDSS-III Baryon Oscillation Spectroscopic Survey: a large sample of mock galaxy catalogues

We present a fast method of producing mock galaxy catalogues that can be used to compute covariance matrices of large-scale clustering measurements and test the methods of analysis. Our method populates a 2nd-order Lagrangian Perturbation Theory (2LPT) matter field, where we calibrate masses of dark matter halos by detailed comparisons with N-body simulations. We demonstrate the clustering of halos is recovered at ~10 per cent accuracy. We populate halos with mock galaxies using a Halo Occupation Distribution (HOD) prescription, which has been calibrated to reproduce the clustering measurements on scales between 30 and 80 Mpc/h. We compare the sample covariance matrix from our mocks with analytic estimates, and discuss differences. We have used this method to make catalogues corresponding to Data Release 9 of the Baryon Oscillation Spectroscopic Survey (BOSS),producing 600 mock catalogues of the "CMASS" galaxy sample. These mocks enabled detailed tests of methods and errors that formed an integral part of companion analyses of these galaxy data.

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

Regarding the Line-of-Sight Baryonic Acoustic Feature in the Sloan Digital Sky Survey and Baryon Oscillation Spectroscopic Survey Luminous Red Galaxy Samples

We analyze the line-of-sight baryonic acoustic feature in the two-point correlation function ξ of the Sloan Digital Sky Survey (SDSS) luminous red galaxy (LRG) sample (0.16 < z < 0.47). By defining a narrow line-of-sight region, rp < 5.5 Mpc/h, where rp is the transverse separation component, we measure a strong excess of clustering at ~ 110 Mpc/h, as previously reported in the literature. We also test these results in an alternative coordinate system, by defining the line-of-sight as θ < 3°, where θ is the opening angle. This clustering excess appears much stronger than the feature in the better-measured monopole. A fiducial ΛCDM non-linear model in redshift-space predicts a much weaker signature. We use realistic mock catalogs to model the expected signal and noise. We find that the line-of-sight measurements can be explained well by our mocks as well as by a featureless ξ = 0. We conclude that there is no convincing evidence that the strong clustering measurement is the line-of-sight baryonic acoustic feature. We also evaluate how detectable such a signal would be in the upcoming Baryon Oscillation Spectroscopic Survey LRG volume (BOSS). Mock LRG catalogs (z < 0.6) suggest that: (i) the narrow line- of-sight cylinder and cone defined above probably will not reveal a detectable acoustic feature in BOSS; (ii) a clustering measurement as high as that in the current sample can be ruled out (or confirmed) at a high confidence level using a BOSS-sized data set; and (iii) an analysis with wider angular cuts, which provide better signal-to-noise ratios, can nevertheless be used to compare line-of-sight and transverse distances, and thereby constrain the expansion rate H(z) and diameter distance DA(z).

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Three-Point Correlation Functions of SDSS Galaxies: Constraining Galaxy-Mass Bias

We constrain the linear and quadratic bias parameters from the configuration dependence of the three-point correlation function (3PCF) in both redshift and projected space, utilizing measurements of spectroscopic galaxies in the Sloan Digital Sky Survey (SDSS) Main Galaxy Sample. We show that bright galaxies (M_r < -21.5) are biased tracers of mass, measured at a significance of 4.5 sigma in redshift space and 2.5 sigma in projected space by using a thorough error analysis in the quasi-linear regime (9-27 Mpc/h). Measurements on a fainter galaxy sample are consistent with an unbiased model. We demonstrate that a linear bias model appears sufficient to explain the galaxy-mass bias of our samples, although a model using both linear and quadratic terms results in a better fit. In contrast, the bias values obtained from the linear model appear in better agreement with the data by inspection of the relative bias, and yield implied values of sigma_8 that are more consistent with current constraints. We investigate the covariance of the 3PCF, which itself is a measurement of galaxy clustering. We assess the accuracy of our error estimates by comparing results from mock galaxy catalogs to jackknife re-sampling methods. We identify significant differences in the structure of the covariance. However, the impact of these discrepancies appears to be mitigated by an eigenmode analysis that can account for the noisy, unresolved modes. Our results demonstrate that using this technique is sufficient to remove potential systematics even when using less-than-ideal methods to estimate errors.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Three-Point Correlation Functions of SDSS Galaxies: Luminosity and Color Dependence in Redshift and Projected Space

The three-point correlation function (3PCF) provides an important view into the clustering of galaxies that is not available to its lower order cousin, the two-point correlation function (2PCF). Higher order statistics, such as the 3PCF, are necessary to probe the non-Gaussian structure and shape information expected in these distributions. We measure the clustering of spectroscopic galaxies in the Main Galaxy Sample of the Sloan Digital Sky Survey (SDSS), focusing on the shape or configuration dependence of the reduced 3PCF in both redshift and projected space. This work constitutes the largest number of galaxies ever used to investigate the reduced 3PCF, using over 220,000 galaxies in three volume-limited samples. We find significant configuration dependence of the reduced 3PCF at 3-27 Mpc/h, in agreement with LCDM predictions and in disagreement with the hierarchical ansatz. Below 6 Mpc/h, the redshift space reduced 3PCF shows a smaller amplitude and weak configuration dependence in comparison with projected measurements suggesting that redshift distortions, and not galaxy bias, can make the reduced 3PCF appear consistent with the hierarchical ansatz. The reduced 3PCF shows a weaker dependence on luminosity than the 2PCF, with no significant dependence on scales above 9 Mpc/h. On scales less than 9 Mpc/h, the reduced 3PCF appears more affected by galaxy color than luminosty. We demonstrate the extreme sensitivity of the 3PCF to systematic effects such as sky completeness and binning scheme, along with the difficulty of resolving the errors. Some comparable analyses make assumptions that do not consistently account for these effects.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Large scale bias and the inaccuracy of the peak-background split

The peak-background split argument is commonly used to relate the abundance of dark matter halos to their spatial clustering. Testing this argument requires an accurate determination of the halo mass function. We present a Maximum Likelihood method for fitting parametric functional forms to halo abundances which differs from previous work because it does not require binned counts. Our conclusions do not depend on whether we use our method or more conventional ones. In addition, halo abundances depend on how halos are defined. Our conclusions do not depend on the choice of link length associated with the friends-of-friends halo-finder, nor do they change if we identify halos using a spherical overdensity algorithm instead. The large scale halo bias measured from the matter-halo cross spectrum b_x and the halo autocorrelation function b_xi (on scales k~0.03h/Mpc and r ~50 Mpc/h) can differ by as much as 5% for halos that are significantly more massive than the characteristic mass M*. At these large masses, the peak background split estimate of the linear bias factor b1 is 3-5% smaller than b_xi, which is 5% smaller than b_x. We discuss the origin of these discrepancies: deterministic nonlinear local bias, with parameters determined by the peak-background split argument, is unable to account for the discrepancies we see. A simple linear but nonlocal bias model, motivated by peaks theory, may also be difficult to reconcile with our measurements. More work on such nonlocal bias models may be needed to understand the nature of halo bias at this level of precision.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Large-Scale Structure in Brane-Induced Gravity II. Numerical Simulations

We use N-body simulations to study the nonlinear structure formation in brane-induced gravity, developing a new method that requires alternate use of Fast Fourier Transforms and relaxation. This enables us to compute the nonlinear matter power spectrum and bispectrum, the halo mass function, and the halo bias. From the simulation results, we confirm the expectations based on analytic arguments that the Vainshtein mechanism does operate as anticipated, with the density power spectrum approaching that of standard gravity within a modified background evolution in the nonlinear regime. The transition is very broad and there is no well defined Vainshtein scale, but roughly this corresponds to k_*~ 2 at redshift z=1 and k_*~ 1 at z=0. We checked that while extrinsic curvature fluctuations go nonlinear, and the dynamics of the brane-bending mode C receives important nonlinear corrections, this mode does get suppressed compared to density perturbations, effectively decoupling from the standard gravity sector. At the same time, there is no violation of the weak field limit for metric perturbations associated with C. We find good agreement between our measurements and the predictions for the nonlinear power spectrum presented in paper I, that rely on a renormalization of the linear spectrum due to nonlinearities in the modified gravity sector. A similar prediction for the mass function shows the right trends. Our simulations also confirm the induced change in the bispectrum configuration dependence predicted in paper I.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

The Baryonic Acoustic Feature and Large-Scale Clustering in the SDSS LRG Sample

We examine the correlation function ξof the Sloan Digital Sky Survey (SDSS) Luminous Red Galaxy sample (LRG) at large scales (60<s<400 Mpc/h) using the final data release (DR7; 105,831 LRGs between 0.16<z<0.47). Using mock catalogs, we demonstrate that the observed baryonic acoustic peak and larger scale signal are consistent with LCDM at the 1.5σlevel. The signal at 155<s<200 Mpc/h tends to be high relative to theoretical expectations; this slight deviation can be attributed to a bright subsample of the LRGs. Fitting data to a non-linear, redshift-space, template based-model, we constrain the peak position at s_p=103.6+3.6-2.4 Mpc/h when fitting the range 60<s<150 Mpc/h (1σuncertainties measured from the mocks. This redshift-space distance s_p is related to the comoving sound horizon scale r_s after taking into account matter clustering non-linearities, redshift distortions and galaxy clustering bias. Mock catalogs show that the probability that a DR7-sized sample would not have an identifiable peak is at least 10%. As a consistency check of a fiducial cosmology, we use the observed s_p to obtain the distance D_V=[(1+z)^2D_A^2cz/H(z)]^(1/3) relative to the acoustic scale. We find r_s/D_V(z=0.278)=0.1394+-0.0049. This result is in excellent agreement with Percival et. al (2009), who examine roughly the same data set, but using the power spectrum. Comparison with other determinations in the literature are also in very good agreement. We have tested our results against a battery of possible systematic effects, finding all effects are smaller than our estimated sample variance.

0 citations0 reviewsNo code linkedOpen access
preprint2004arXiv

Galaxy Bias and Halo-Occupation Numbers from Large-Scale Clustering

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

0 citations0 reviewsNo code linkedOpen access
preprint2000arXiv

A New Angle on Gravitational Clustering

We describe a new approach to gravitational instability in large-scale structure, where the equations of motion are written and solved as in field theory in terms of Feynman diagrams. The basic objects of interest are the propagator (which propagates solutions forward in time), the vertex (which describes non-linear interactions between waves) and a source with prescribed statistics which describes the effect of initial conditions. We show that loop corrections renormalize these quantities, and discuss applications of this formalism to a better understanding of gravitational instability and to improving non-linear perturbation theory in the transition to the non-linear regime. We also consider the role of vorticity creation due to shell-crossing and show using N-body simulations that at small (virialized) scales the velocity field reaches equipartition, i.e. the vorticity power spectrum is about twice the divergence power spectrum.

0 citations0 reviewsNo code linkedOpen access
preprint1997arXiv

Non-Linear Evolution of the Bispectrum of Cosmological Perturbations

The bispectrum, the three-point function of density fluctuations in Fourier space, is the lowest order statistic that carries information about the spatial coherence of large-scale structures. For Gaussian initial conditions, when the density fluctuation amplitude is small, tree-level (leading order) perturbation theory predicts a characteristic dependence of the bispectrum on the shape of the triangle formed by the three wave vectors. This configuration dependence provides a signature of gravitational instability, and departures from it in galaxy catalogs can be interpreted as due to bias, that is, non-gravitational effects. On the other hand, N-body simulations indicate that the reduced three-point function becomes relatively shape-independent in the strongly non-linear regime. In order to understand this non-linear transition and assess the domain of reliability of shape-dependence as a probe of bias, we calculate the one-loop (next-to-leading order) corrections to the bispectrum in perturbation theory. We compare these results with measurements in numerical simulations with scale-free and CDM initial power spectra. We find that the one-loop corrections account very well for the departures from the tree-level results measured in numerical simulations on weakly non-linear scales. In this regime, the reduced bispectrum qualitatively retains its tree-level shape, but the amplitude can change significantly. At smaller scales, the reduced bispectrum in the simulations starts to flatten, an effect which can be partially understood from the one-loop results. In the strong clustering regime, where perturbation theory breaks down entirely, the simulation results confirm that the reduced bispectrum has almost no dependence on triangle shape, in rough agreement with the hierarchical ansatz.

0 citations0 reviewsNo code linkedOpen access
preprint1995arXiv

Loop Corrections in Non-Linear Cosmological Perturbation Theory

Using a diagrammatic approach to Eulerian perturbation theory, we analytically calculate the variance and skewness of the density and velocity divergence induced by gravitational evolution from Gaussian initial conditions, including corrections *beyond* leading order. Except for the power spectrum, previous calculations in cosmological perturbation theory have been confined to leading order (tree level)-we extend these to include loop corrections. For scale-free initial power spectra, the one-loop variance σ^2 = σ^2_l + 1.82 σ^4_l and the skewness S_3 = 34/7 + 9.8 σ^2_l, where σ_l is the rms fluctuation of the linear density field. We also compute loop corrections to the variance, skewness, and kurtosis for several non-linear approximation schemes, where the calculation can be easily generalized to 1-point cumulants of higher order and arbitrary number of loops. We find that the Zel'dovich approximation gives the best approximation to the loop corrections of exact perturbation theory, followed by the Linear Potential approximation (LPA) and the Frozen Flow approximation (FFA), in qualitative agreement with the relative behavior of tree-level results. In LPA and FFA, loop corrections are infrared divergent for spectral indices n < 0; this is related to the breaking of Galilean invariance in these schemes.

0 citations0 reviewsNo code linkedOpen access