Search

We present non-linear solutions of Vlasov Perturbation Theory (VPT), describing gravitational clustering of collisionless dark matter with dispersion and higher cumulants induced by orbit crossing. We show that VPT can be cast into a form that is formally analogous to standard perturbation theory (SPT), but including additional perturbation variables, non-linear interactions, and a more complex propagation. VPT non-linear kernels have a crucial decoupling property: for fixed total momentum, the kernels becomes strongly suppressed when any of the individual momenta cross the dispersion scale into the non-linear regime. This screening of UV modes allows us to compute non-linear corrections to power spectra even for cosmologies with very blue power-law input spectra, for which SPT diverges. We compare predictions for the density and velocity divergence power spectra as well as the bispectrum at one-loop order to N-body results in a scaling universe with spectral indices $-1\leq n_s\leq +2$. We find a good agreement up to the non-linear scale for all cases, with a reach that increases with the spectral index $n_s$. We discuss the generation of vorticity as well as vector and tensor modes of the velocity dispersion, showing that neglecting vorticity when including dispersion would lead to a violation of momentum conservation. We verify momentum conservation when including vorticity, and compute the vorticity power spectrum at two-loop order, necessary to recover the correct large-scale limit with slope $n_w=2$. Comparing to our N-body measurements confirms the cross-over from $k^4$ to $k^2$ scaling on large scales. Our results provide a proof-of-principle that perturbative techniques for dark matter clustering can be systematically improved based on the known underlying collisionless dynamics., Comment: 37 pages + appendices, 23 figures

The standard perturbation theory (SPT) approach to gravitational clustering is based on a fluid approximation of the underlying Vlasov-Poisson dynamics, taking only the zeroth and first cumulant of the phase-space distribution function into account (density and velocity fields). This assumption breaks down when dark matter particle orbits cross and leads to well-known problems, e.g. an anomalously large backreaction of small-scale modes onto larger scales that compromises predictivity. We extend SPT by incorporating second and higher cumulants generated by orbit crossing. For collisionless matter, their equations of motion are completely fixed by the Vlasov-Poisson system, and thus we refer to this approach as Vlasov Perturbation Theory (VPT). Even cumulants develop a background value, and they enter the hierarchy of coupled equations for the fluctuations. The background values are in turn sourced by power spectra of the fluctuations. The latter can be brought into a form that is formally analogous to SPT, but with an extended set of variables and linear as well as non-linear terms, that we derive explicitly. In this paper, we focus on linear solutions, which are far richer than in SPT, showing that modes that cross the dispersion scale set by the second cumulant are highly suppressed. We derive stability conditions on the background values of even cumulants from the requirement that exponential instabilities be absent. We also compute the expected magnitude of averaged higher cumulants for various halo models and show that they satisfy the stability conditions. Finally, we derive self-consistent solutions of perturbations and background values for a scaling universe and study the convergence of the cumulant expansion. The VPT framework provides a conceptually straightforward and deterministic extension of SPT that accounts for the decoupling of small-scale modes., Comment: 31 pages + appendices, 10 figures

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., Comment: 17 pages, 8 figures. Minor changes. Matches version published in PRD

In this paper we present COMET, a Gaussian process emulator of the galaxy power spectrum multipoles in redshift-space. The model predictions are based on one-loop perturbation theory and we consider two alternative descriptions of redshift-space distortions: one that performs a full expansion of the real- to redshift-space mapping, as in recent effective field theory models, and another that preserves the non-perturbative impact of small-scale velocities by means of an effective damping function. The outputs of COMET can be obtained at arbitrary redshifts (up to $z \sim 3$), for arbitrary fiducial background cosmologies, and for a large parameter space that covers the shape parameters $\omega_c$, $\omega_b$, and $n_s$, as well as the evolution parameters $h$, $A_s$, $\Omega_K$, $w_0$, and $w_a$. This flexibility does not impair COMET's accuracy, since we exploit an exact degeneracy between the evolution parameters that allows us to train the emulator on a significantly reduced parameter space. While the predictions are sped up by at least two orders of magnitude, validation tests reveal an accuracy of $0.1\,\%$ for the monopole and quadrupole ($0.3\,\%$ for the hexadecapole), or alternatively, better than $0.25\,\sigma$ for all three multipoles in comparison to statistical uncertainties expected for the Euclid survey with a tenfold increase in volume. We show that these differences translate into shifts in mean posterior values that are at most of the same size, meaning that COMET can be used with the same confidence as the exact underlying models. COMET is a publicly available Python package that also provides the tree-level bispectrum multipoles in redshift-space and Gaussian covariance matrices., Comment: 19 pages, 11 figures; for the COMET Python package, see https://gitlab.com/aegge/comet-emu; v2: added discussion on binning and discreteness effects, minor improvements to match published version

Galaxy power spectrum and bispectrum signals are distorted by peculiar velocities and other relativistic effects arising from a perturbed spacetime background. In addition, study of correlation functions of tracers in Fourier space is often done in the plane-parallel approximation under which it is assumed that line-of-sight (LOS) vectors are parallel. In this work we show that a simple perturbative procedure can be employed for a fast evaluation of beyond plane-parallel (wide-angle) corrections to the power spectrum and bispectrum. We also show that evolution of linear matter density fluctuations in a relativistic context can be found from a simple method. For the power spectrum at linear level, we compare leading order wide-angle contributions to multipoles of the galaxy power spectrum with those from non-integrated and integrated relativistic corrections and estimate their possible contamination on local fNL measurements to be of order a few. We also compute wide-angle corrections in the presence of nonlinear terms at one-loop order. For the bispectrum, we show that wide-angle effects alone, even with fully symmetric choices of LOS, give rise to imaginary, odd-parity multipoles of the galaxy bispectrum (dipole, octupole, etc.) which are in many cases larger than previously known ones of relativistic origin. We calculate these contributions and provide an estimator for measuring the leading order bispectrum dipole from data, using a symmetric LOS definition. Finally, we calculate the leading order corrections to multipoles of real plane-parallel bispectrum multipoles and estimate the apparent local fNL induced to be of order unity., Comment: Minor improvements and corrections, discussion added about inclusion of window function, matches the published version

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 $\Lambda$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 $\sigma_{12}=0.805\pm 0.049$, while using the traditional reference scale of $8\,h^{-1}{\rm Mpc}$ gives $\sigma_{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\sigma$. 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\sigma$, 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 $\omega_{\rm DE}=0.335 \pm 0.011$, which is 1.7$\sigma$ higher than the one preferred by Planck., Comment: 14 pages, 10 figures, matches the version accepted by MNRAS

We investigate the impact of different assumptions in the modeling of one-loop galaxy bias on the recovery of cosmological parameters, as a follow up of the analysis done in the first paper of the series at fixed cosmology. We use three different synthetic galaxy samples whose clustering properties match the ones of the CMASS and LOWZ catalogues of BOSS and the SDSS Main Galaxy Sample. We investigate the relevance of allowing for either short range non-locality or scale-dependent stochasticity by fitting the real-space galaxy auto power spectrum or the combination of galaxy-galaxy and galaxy-matter power spectrum. From a comparison among the goodness-of-fit ($\chi^2$), unbiasedness of cosmological parameters (FoB), and figure-of-merit (FoM), we find that a four-parameter model (linear, quadratic, cubic non-local bias, and constant shot-noise) with fixed quadratic tidal bias provides a robust modelling choice for the auto power spectrum of the three samples, up to $k_{\rm max}=0.3\,h\,\mathrm{Mpc}^{-1}$ and for an effective volume of $6\,h^{-3}\,\mathrm{Gpc}^3$. Instead, a joint analysis of the two observables fails at larger scales, and a model extension with either higher derivatives or scale-dependent shot-noise is necessary to reach a similar $k_{\rm max}$, with the latter providing the most stable results. These findings are obtained with three, either hybrid or perturbative, prescriptions for the matter power spectrum, \texttt{RESPRESSO}, gRPT and EFT. In all cases, the inclusion of scale-dependent shot-noise increases the range of validity of the model in terms of FoB and $\chi^2$. Interestingly, these model extensions with additional free parameters do not necessarily lead to an increase in the maximally achievable FoM for the cosmological parameters $\left(h,\,\Omega_ch^2,\,A_s\right)$, which are generally consistent to those of the simpler model at smaller $k_{\rm max}$., Comment: 24 pages, 8 figures

We present a joint likelihood analysis of the real-space power spectrum and bispectrum measured from a variety of halo and galaxy mock catalogs. A novel aspect of this work is the inclusion of nonlinear triangle configurations for the bispectrum, made possible by a complete next-to-leading order ("one-loop") description of galaxy bias, as is already common practice for the power spectrum. Based on the goodness-of-fit and the unbiasedness of the parameter posteriors, we accomplish a stringent validation of this model compared to the leading order ("tree-level") bispectrum. Using measurement uncertainties that correspond to an effective survey volume of $6\,(\mathrm{Gpc}/h)^3$, we determine that the one-loop corrections roughly double the applicable range of scales, from $\sim 0.17\,h/\mathrm{Mpc}$ (tree-level) to $\sim 0.3\,h/\mathrm{Mpc}$. This converts into a $1.5 - 2$x improvement on constraints of the linear bias parameter at fixed cosmology, and a $1.5 - 2.4$x shrinkage of uncertainties on the amplitude of fluctuations $A_s$, which clearly demonstrates the benefit of extracting information from nonlinear scales despite having to marginalize over a larger number of bias parameters. Besides, our precise measurements of galaxy bias parameters up to fourth order allow for thorough comparisons to coevolution relations, showing excellent agreement for all contributions generated by the nonlocal action of gravity. Using these relations in the likelihood analysis does not compromise the model validity and is crucial for obtaining the quoted improvements on $A_s$. We also analyzed the impact of higher-derivative and scale-dependent stochastic terms, finding that for a subset of our tracers the former can boost the performance of the tree-level model with constraints on $A_s$ that are only slightly degraded compared to the one-loop model., Comment: 33 pages, 13 figures, supporting Mathematica notebook available at https://doi.org/10.5281/zenodo.4529886, Python package for Delaunay binning correction available at https://github.com/aegge/BispTools

We use analytic covariance matrices to carry out a full-shape analysis of the galaxy power spectrum multipoles from the Baryon Oscillation Spectroscopic Survey (BOSS). We obtain parameter estimates that agree well with those based on the sample covariance from two thousand galaxy mock catalogs, thus validating the analytic approach and providing substantial reduction in computational cost. We also highlight a number of additional advantages of analytic covariances. First, the analysis does not suffer from sampling noise, which biases the constraints and typically requires inflating parameter error bars. Second, it allows us to study convergence of the cosmological constraints when recomputing the analytic covariances to match the best-fit power spectrum, which can be done at a negligible computational cost, unlike when using mock catalogs. These effects reduce the systematic error budget of cosmological constraints, which suggests that the analytic approach may be an important tool for upcoming high-precision galaxy redshift surveys such as DESI and Euclid. Finally, we study the impact of various ingredients in the power spectrum covariance matrix and show that the non-Gaussian part, which includes the regular trispectrum and super-sample covariance, has a marginal effect ($\lesssim 10 \%$) on the cosmological parameter error bars. We also suggest improvements to analytic covariances that are commonly used in Fisher forecasts., Comment: 24 pages, 12 figures. Version appearing in PRD. Likelihoods available at https://github.com/Michalychforever/lss_montepython . CovaPT code for calculating analytic covariances is available at https://github.com/JayWadekar/CovaPT

We measure the anisotropic clustering of the quasar sample from Data Release 16 (DR16) of the Sloan Digital Sky Survey IV extended Baryon Oscillation Spectroscopic Survey (eBOSS). A sample of $343,708$ spectroscopically confirmed quasars between redshift $0.8

We test the regime of validity of one-loop galaxy bias for a wide variety of biased tracers. Our most stringent test asks the bias model to simultaneously match the galaxy-galaxy and galaxy-mass spectrum, using the measured nonlinear matter spectrum from the simulations to test one-loop effects from the bias expansion alone. In addition, we investigate the relevance of short-range nonlocality and halo exclusion through higher-derivative and scale-dependent noise terms, as well as the impact of using co-evolution relations to reduce the number of free fitting parameters. From comparing validity and merit of these assumptions we find that a four-parameter model (linear, quadratic, cubic nonlocal bias, and constant shot noise) with fixed quadratic tidal bias provides a robust modeling choice for the auto power spectrum of the less massive halos in our set of samples and their galaxy populations (up to $k_{\mathrm{max}} = 0.35\,h/\mathrm{Mpc}$ for a sample volume of $6\,(\mathrm{Gpc}/h)^3$). For the more biased tracers it is most beneficial to include scale-dependent noise. This is also the preferred option when considering combinations of the auto and cross power spectrum, which might be relevant in joint studies of galaxy clustering and weak lensing. We also test the use of perturbation theory to account for matter loops through gRPT, EFT and the hybrid approach RESPRESSO. While all these have similar performance, we find the latter to be the best in terms of validity and recovered mean posterior values, in accordance with it being based partially on simulations., Comment: 30 pages, 12 figures; v2: includes gRPT expressions and some minor edits, matches published version

We compute the covariance of the galaxy power spectrum multipoles in perturbation theory, including the effects of nonlinear evolution, nonlinear and nonlocal bias, radial redshift-space distortions, arbitrary survey window and shot noise. We rewrite the power spectrum FKP estimator in terms of the usual windowed galaxy fluctuations and the fluctuations in the number of galaxies inside the survey volume. We show that this leads to a stronger super-sample covariance than assumed in the literature and causes a substantial leakage of Gaussian information. We decompose the covariance matrix into several contributions that provide an insight into its behavior for different biased tracers. We show that for realistic surveys, the covariance of power spectrum multipoles is already dominated by shot noise and super survey mode coupling in the weakly non-linear regime. Both these effects can be accurately modeled analytically, making a perturbative treatment of the covariance very compelling. Our method allows for the covariance to be varied as a function of cosmology and bias parameters very efficiently, with survey geometry entering as fixed kernels that can be computed separately using fast fourier transforms (FFTs). We find excellent agreement between our analytic covariance and that estimated from BOSS DR12 Patchy mock catalogs in the whole range we tested, up to $k=0.6$ h/Mpc. This bodes well for application to future surveys such as DESI and Euclid. The CovaPT code that accompanies this paper is available at https://github.com/JayWadekar/CovaPT, Comment: 28 pages + appendices; 14 figures. Version appearing in PRD. The busy reader can focus on the sections: 5C, 6 and 7 which contain the main results. CovaPT code is available at https://github.com/JayWadekar/CovaPT

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_\nu$, a key parameter in particle physics beyond the Standard Model. However, degeneracies among cosmological parameters, especially between $M_\nu$ and $\sigma_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_\nu$--$\sigma_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 $\{\Omega_m$, $\Omega_b$, $h$, $n_s$, $\sigma_8$, $M_\nu\}$. For $k_{\rm max}{=}0.5~h/{\rm Mpc}$, the bispectrum provides $\Omega_m$, $\Omega_b$, $h$, $n_s$, and $\sigma_8$ constraints 1.9, 2.6, 3.1, 3.6, and 2.6 times tighter than the power spectrum. For $M_\nu$, the bispectrum improves the 1$\sigma$ 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_\nu$ 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_\nu$., Comment: 33 pages, 13 figures

The Quijote simulations are a set of 44,100 full N-body simulations spanning more than 7,000 cosmological models in the $\{\Omega_{\rm m}, \Omega_{\rm b}, h, n_s, \sigma_8, M_\nu, w \}$ hyperplane. At a single redshift the simulations contain more than 8.5 trillions of particles over a combined volume of 44,100 $(h^{-1}{\rm Gpc})^3$; each simulation follow the evolution of $256^3$, $512^3$ or $1024^3$ particles in a box of $1~h^{-1}{\rm Gpc}$ length. Billions of dark matter halos and cosmic voids have been identified in the simulations, whose runs required more than 35 million core hours. The Quijote simulations have been designed for two main purposes: 1) to quantify the information content on cosmological observables, and 2) to provide enough data to train machine learning algorithms. In this paper we describe the simulations and show a few of their applications. We also release the Petabyte of data generated, comprising hundreds of thousands of simulation snapshots at multiple redshifts, halo and void catalogs, together with millions of summary statistics such as power spectra, bispectra, correlation functions, marked power spectra, and estimated probability density functions., Comment: 20 pages, 15 figures. Matches published version. Simulations publicly available at https://github.com/franciscovillaescusa/Quijote-simulations

In this paper, we verify the large scale structure consistency relations using N-body simulations, including modes in the highly non-linear regime. These relations (pointed out by Kehagias & Riotto and Peloso & Pietroni) follow from the symmetry of the dynamics under a shift of the Newtonian potential by a constant and a linear gradient, and predict the absence of certain poles in the ratio between the (equal time) squeezed bispectrum and power spectrum. The consistency relations, as symmetry statements, are exact, but have not been previously checked beyond the perturbative regime. Our test using N-body simulations not only offers a non-perturbative check, but also serves as a warm-up exercise for applications to observational data. A number of subtleties arise when taking the squeezed limit of the bispectrum--we show how to circumvent or address them. An interesting by-product of our investigation is an explicit demonstration that the linear-gradient symmetry is unaffected by the periodic boundary condition of the simulations. Lastly, we verify using simulations that the consistency relations are violated when the initial conditions are non-gaussian (of the local fNL type). The methodology developed here paves the way for constraining primordial non-gaussianity using large scale structure data, including (numerous) highly non-linear modes that are otherwise hard to interpret and utilize., Comment: 10 pages, 5 figures, 1 table

Combination of the power spectrum and bispectrum is a powerful way of breaking degeneracies between galaxy bias and cosmological parameters, enabling us to maximize the constraining power from galaxy surveys. Recent cosmological constraints treat the power spectrum and bispectrum on an uneven footing: they include one-loop bias corrections for the power spectrum but not the bispectrum. To bridge this gap, we develop the galaxy bias description up to fourth order in perturbation theory, conveniently expressed through a basis of Galilean invariants that clearly split contributions that are local and nonlocal in the second derivatives of the linear gravitational potential. In addition, we consider relevant contributions from short-range nonlocality (higher-derivative terms), stress-tensor corrections and stochasticity. To sidestep the usual renormalization of bias parameters that complicates predictions beyond leading order, we recast the bias expansion in terms of multipoint propagators, which take a simple form in our split-basis with loop corrections depending only on bias parameters corresponding to nonlocal operators. We show how to take advantage of Galilean invariance to compute the time evolution of bias and present results for the fourth-order parameters for the first time. We also discuss the possibilities of reducing the bias parameter space by using the evolution of bias and exploiting degeneracies between various bias contributions in the large-scale limit. Our baseline model allows to verify these simplifications for any application to large-scale structure data sets., Comment: 36 pages, 6 figures; added two figures and various references, some minor changes throughout

We explore the cosmological implications of anisotropic clustering measurements of the quasar sample from Data Release 14 of the Sloan Digital Sky Survey IV Extended Baryon Oscillation Spectroscopic Survey (eBOSS) in configuration space. The $\sim 147,000$ quasar sample observed by eBOSS offers a direct tracer of the density field and bridges the gap of previous BAO measurements between redshift $0.8

We compare reduced three-point correlations $Q$ of matter, haloes (as proxies for galaxies) and their cross correlations, measured in a total simulated volume of $\sim100 \ (h^{-1} \text{Gpc})^{3}$, to predictions from leading order perturbation theory on a large range of scales in configuration space. Predictions for haloes are based on the non-local bias model, employing linear ($b_1$) and non-linear ($c_2$, $g_2$) bias parameters, which have been constrained previously from the bispectrum in Fourier space. We also study predictions from two other bias models, one local ($g_2=0$) and one in which $c_2$ and $g_2$ are determined by $b_1$ via an approximately universal relation. Overall, measurements and predictions agree when $Q$ is derived for triangles with $(r_1r_2r_3)^{1/3} \gtrsim 60 h^{-1}\text{Mpc}$, where $r_{1-3}$ are the sizes of the triangle legs. Predictions for $Q_{matter}$, based on the linear power spectrum, show significant deviations from the measurements at the BAO scale (given our small measurement errors), which strongly decrease when adding a damping term or using the non-linear power spectrum, as expected. Predictions for $Q_{halo}$ agree best with measurements at large scales when considering non-local contributions. The universal bias model works well for haloes and might therefore be also useful for tightening constraints on $b_1$ from $Q$ in galaxy surveys. Such constraints are independent of the amplitude of matter density fluctuation ($\sigma_8$) and hence break the degeneracy between $b_1$ and $\sigma_8$, present in galaxy two-point correlations., Comment: 18 pages, 20 figures, 2 tables, added references and systematic tests, accepted for publication in MNRAS

Interpreting the small-scale clustering of galaxies with halo models can elucidate the connection between galaxies and dark matter halos. Unfortunately, the modelling is typically not sufficiently accurate for ruling out models statistically. It is thus difficult to use the information encoded in small scales to test cosmological models or probe subtle features of the galaxy-halo connection. In this paper, we attempt to push halo modelling into the "accurate" regime with a fully numerical mock-based methodology and careful treatment of statistical and systematic errors. With our forward-modelling approach, we can incorporate clustering statistics beyond the traditional two-point statistics. We use this modelling methodology to test the standard $\Lambda\mathrm{CDM}$ + halo model against the clustering of SDSS DR7 galaxies. Specifically, we use the projected correlation function, group multiplicity function and galaxy number density as constraints. We find that while the model fits each statistic separately, it struggles to fit them simultaneously. Adding group statistics leads to a more stringent test of the model and significantly tighter constraints on model parameters. We explore the impact of varying the adopted halo definition and cosmological model and find that changing the cosmology makes a significant difference. The most successful model we tried (Planck cosmology with Mvir halos) matches the clustering of low luminosity galaxies, but exhibits a 2.3$\sigma$ tension with the clustering of luminous galaxies, thus providing evidence that the "standard" halo model needs to be extended. This work opens the door to adding interesting freedom to the halo model and including additional clustering statistics as constraints., Comment: Replaced to match the published version

(abridged) We investigate the signatures left by the cosmic neutrino background on the clustering of matter, CDM+baryons and halos in redshift-space using a set of more than 1000 N-body and hydrodynamical simulations with massless and massive neutrinos. We find that the effect neutrinos induce on the clustering of CDM+baryons in redshift-space on small scales is almost entirely due to the change in $\sigma_8$. Neutrinos imprint a characteristic signature in the quadrupole of the matter (CDM+baryons+neutrinos) field on small scales, that can be used to disentangle the effect of $\sigma_8$ and $M_\nu$. We show that the effect of neutrinos on the clustering of halos is very different, on all scales, to the one induced by $\sigma_8$. We find that the effects of neutrinos of the growth rate of CDM+baryons ranges from $\sim0.3\%$ to $2\%$ on scales $k\in[0.01, 0.5]~h{\rm Mpc}^{-1}$ for neutrinos with masses $M_\nu \leqslant 0.15$ eV. We compute the bias between the momentum of halos and the momentum of CDM+baryon and find it to be 1 on large scales for all models with massless and massive neutrinos considered. This point towards a velocity bias between halos and total matter on large scales that it is important to account for in order to extract unbiased neutrino information from velocity/momentum surveys such as kSZ observations. We show that baryonic effects can affect the clustering of matter and CDM+baryons in redshift-space by up to a few percent down to $k=0.5~h{\rm Mpc}^{-1}$. We find that hydrodynamics and astrophysical processes, as implemented in our simulations, only distort the relative effect that neutrinos induce on the anisotropic clustering of matter, CDM+baryons and halos in redshift-space by less than $1\%$. Thus, the effect of neutrinos in the fully non-linear regime can be written as a transfer function with very weak dependence on astrophysics., Comment: 19 pages, 12 figures

The protohalo patches from which halos form are defined by a number of constraints imposed on the Lagrangian dark matter density field. Each of these constraints contributes to biasing the spatial distribution of the protohalos relative to the matter. We show how measurements of this spatial distribution -- linear combinations of protohalo bias factors -- can be used to make inferences about the physics of halo formation. Our analysis exploits the fact that halo bias factors satisfy consistency relations which encode this physics, and that these relations are the same even for sub-populations in which assembly bias has played a role. We illustrate our methods using a model in which three parameters matter: a density threshold, the local slope and the curvature of the smoothed density field. The latter two are nearly degenerate; our approach naturally allows one to build an accurate effective two-parameter model for which the consistency relations still apply. This, with an accurate description of the smoothing window, allows one to describe the protohalo-matter cross-correlation very well, both in Fourier and configuration space. We then use our determination of the large scale bias parameters together with the consistency relations, to estimate the enclosed density and mean slope on the Lagrangian radius scale of the protohalos. Direct measurements of these quantities, made on smaller scales than those on which the bias parameters are typically measured, are in good agreement., Comment: 17 pages, 19 figures, match the published version

Fiber-fed multi-object spectroscopic surveys, with their ability to collect an unprecedented number of redshifts, currently dominate large-scale structure studies. However, physical constraints limit these surveys from successfully collecting redshifts from galaxies too close to each other on the focal plane. This ultimately leads to significant systematic effects on galaxy clustering measurements. Using simulated mock catalogs, we demonstrate that fiber collisions have a significant impact on the power spectrum, $P(k)$, monopole and quadrupole that exceeds sample variance at scales smaller than $k\sim0.1~h/Mpc$. We present two methods to account for fiber collisions in the power spectrum. The first, statistically reconstructs the clustering of fiber collided galaxy pairs by modeling the distribution of the line-of-sight displacements between them. It also properly accounts for fiber collisions in the shot-noise correction term of the $P(k)$ estimator. Using this method, we recover the true $P(k)$ monopole of the mock catalogs with residuals of $<0.5\%$ at $k=0.3~h/Mpc$ and $<4\%$ at $k=0.83~h/Mpc$ -- a significant improvement over existing correction methods. The quadrupole, however, does not improve significantly. The second method models the effect of fiber collisions on the power spectrum as a convolution with a configuration space top-hat function that depends on the physical scale of fiber collisions. It directly computes theoretical predictions of the fiber-collided $P(k)$ multipoles and reduces the influence of smaller scales to a set of nuisance parameters. Using this method, we reliably model the effect of fiber collisions on the monopole and quadrupole down to the scale limits of theoretical predictions. The methods we present in this paper will allow us to robustly analyze galaxy power spectrum multipole measurements to much smaller scales than previously possible., Comment: 18 pages, 10 figures

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 $\Lambda$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 $\Omega_{k}=-0.0007\pm 0.0030$ and the sum of the neutrino masses to $\sum m_{\nu} < 0.25\,{\rm eV}$ at 95 per cent CL. We explore the constraints on the growth rate of cosmic structures assuming $f(z)=\Omega_{\rm m}(z)^\gamma$ and obtain $\gamma = 0.609\pm 0.079$, in good agreement with the predictions of general relativity of $\gamma=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\sigma_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 $\Lambda$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., Comment: 20 pages, 16 figures

We extract cosmological information from the anisotropic power spectrum measurements from the recently completed Baryon Oscillation Spectroscopic Survey (BOSS), extending the concept of clustering wedges to Fourier space. Making use of new FFT-based estimators, we measure the power spectrum clustering wedges of the BOSS sample by filtering out the information of Legendre multipoles l > 4. Our modelling of these measurements is based on novel approaches to describe non-linear evolution, bias, and redshift-space distortions, which we test using synthetic catalogues based on large-volume N-body simulations. We are able to include smaller scales than in previous analyses, resulting in tighter cosmological constraints. Using three overlapping redshift bins, we measure the angular diameter distance, the Hubble parameter, and the cosmic growth rate, and explore the cosmological implications of our full shape clustering measurements in combination with CMB and SN Ia data. Assuming a {\Lambda}CDM cosmology, we constrain the matter density to {\Omega}_m = 0.311 -0.010 +0.009 and the Hubble parameter to H_0 = 67.6 -0.6 +0.7 km s^-1 Mpc^-1, at a confidence level (CL) of 68 per cent. We also allow for non-standard dark energy models and modifications of the growth rate, finding good agreement with the {\Lambda}CDM paradigm. For example, we constrain the equation-of-state parameter to w = -1.019 -0.039 +0.048. 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., Comment: 23 pages, 16 figures in the main text, appendix of 6 pages and 5 figures; revision submitted to MNRAS. The data used in this analysis is publicly available at https://www.sdss3.org/science/boss_publications.php

We investigate the cosmological implications of studying galaxy clustering using a tomographic approach applied to the final BOSS DR12 galaxy sample, including both auto- and cross-correlation functions between redshift shells. We model the signal of the full shape of the angular correlation function, $\omega(\theta)$, in redshift bins using state-of-the-art modelling of non-linearities, bias and redshift-space distortions. We present results on the redshift evolution of the linear bias of BOSS galaxies, which cannot be obtained with traditional methods for galaxy-clustering analysis. We also obtain constraints on cosmological parameters, combining this tomographic analysis with measurements of the cosmic microwave background (CMB) and type Ia supernova (SNIa). We explore a number of cosmological models, including the standard $\Lambda$CDM model and its most interesting extensions, such as deviations from $w_\rm{DE} = -1$, non-minimal neutrino masses, spatial curvature and deviations from general relativity using the growth-index $\gamma$ parametrisation. These results are, in general, comparable to the most precise present-day constraints on cosmological parameters, and show very good agreement with the standard model. In particular, combining CMB, $\omega(\theta)$ and SNIa, we find a value of $w_\rm{DE}$ consistent with $-1$ to a precision better than 5\% when it is assumed to be constant in time, and better than 6\% when we also allow for a spatially-curved Universe., Comment: 21 pages, 18 figures, accepted for publication MNRAS. The data used in this analysis is publicly available at https://sdss3.org/science/boss_publications.php

We apply the Alcock-Paczynski (AP) test to the stacked voids identified using the large-scale structure galaxy catalog from the Baryon Oscillation Spectroscopic Survey (BOSS). This galaxy catalog is part of the Sloan Digital Sky Survey (SDSS) Data Release 12 and is the final catalog of SDSS-III. We also use 1000 mock galaxy catalogs that match the geometry, density, and clustering properties of the BOSS sample in order to characterize the statistical uncertainties of our measurements and take into account systematic errors such as redshift space distortions. For both BOSS data and mock catalogs, we use the ZOBOV algorithm to identify voids, we stack together all voids with effective radii of 30-100Mpc/h in the redshift range 0.43-0.7, and we accurately measure the shape of the stacked voids. Our tests with the mock catalogs show that we measure the stacked void ellipticity with a statistical precision of 2.6%. We find that the stacked voids in redshift space are slightly squashed along the line of sight, which is consistent with previous studies. We repeat this measurement of stacked void shape in the BOSS data assuming several values of Omega_m within the flat LCDM model, and we compare to the mock catalogs in redshift space in order to perform the AP test. We obtain a constraint of $\Omega_m = 0.38^{+0.18}_{-0.15}$ at the 68% confidence level from the AP test. We discuss the various sources of statistical and systematic noise that affect the constraining power of this method. In particular, we find that the measured ellipticity of stacked voids scales more weakly with cosmology than the standard AP prediction, leading to significantly weaker constraints. We discuss how AP constraints will improve in future surveys with larger volumes and densities., Comment: 10 pages, 7 figures. Submitted to the ApJ. Replaced with updated version

We present a cosmic void catalog using the large-scale structure galaxy catalog from the Baryon Oscillation Spectroscopic Survey (BOSS). This galaxy catalog is part of the Sloan Digital Sky Survey (SDSS) Data Release 12 and is the final catalog of SDSS-III. We take into account the survey boundaries, masks, and angular and radial selection functions, and apply the ZOBOV void finding algorithm to the galaxy catalog. We identify a total of 10,643 voids. After making quality cuts to ensure that the voids represent real underdense regions, we obtain 1,228 voids with effective radii spanning the range 20-100Mpc/h and with central densities that are, on average, 30% of the mean sample density. We release versions of the catalogs both with and without quality cuts. We discuss the basic statistics of voids, such as their size and redshift distributions, and measure the radial density profile of the voids via a stacking technique. In addition, we construct mock void catalogs from 1000 mock galaxy catalogs, and find that the properties of BOSS voids are in good agreement with those in the mock catalogs. We compare the stellar mass distribution of galaxies living inside and outside of the voids, and find no significant difference. These BOSS and mock void catalogs are useful for a number of cosmological and galaxy environment studies., Comment: 7 pages, 7 figures, submitted to the ApJ

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

We report a measurement of the large-scale 3-point correlation function of galaxies using the largest dataset for this purpose to date, 777, 202 Luminous Red Galaxies in the Sloan Digital Sky Survey Baryon Acoustic Oscillation Spectroscopic Survey (SDSS BOSS) DR12 CMASS sample. This work exploits the novel algorithm of Slepian & Eisenstein (2015b) to compute the multipole moments of the 3PCF in $\mathcal{O}(N^2)$ time, with $N$ the number of galaxies. Leading-order perturbation theory models the data well in a compressed basis where one triangle side is integrated out. We also present an accurate and computationally efficient means of estimating the covariance matrix. With these techniques the redshift-space linear and non-linear bias are measured, with 2.6% precision on the former if $\sigma_8$ is fixed. The data also indicates a $2.8\sigma$ preference for the BAO, confirming the presence of BAO in the 3-point function., Comment: 16 pages, 10 figures, submitted MNRAS

The window function for protohalos in Lagrangian space is often assumed to be a tophat in real space. We measure this profile directly and find that it is more extended than a tophat but less extended than a Gaussian; its shape is well-described by rounding the edges of the tophat by convolution with a Gaussian that has a scale length about 5 times smaller. This effective window $W_{\rm eff}$ is particularly simple in Fourier space, and has an analytic form in real space. Together with the excursion set bias parameters, $W_{\rm eff}$ describes the scale-dependence of the Lagrangian halo-matter cross correlation up to $kR_{\rm Lag} \sim 10 $, where $R_{\rm Lag}$ is the Lagrangian size of the protohalo. Moreover, with this $W_{\rm eff}$, all the spectral moments of the power spectrum are finite, allowing a straightforward estimate of the excursion set peak mass function. This estimate requires a prescription of the critical overdensity enclosed within a protohalo if it is to collapse, which we calibrate from simulations. We find that the resulting estimate of halo abundances is only accurate to about 20%, and we discuss why: A tophat in `infall time' towards the protohalo center need not correspond to a tophat in the initial spatial distribution, so models in which infall rather than smoothed overdensity is the relevant variable may be more accurate., Comment: 15 pages, 14 figures, matched the published version, discussion and presentation significantly improved, conclusions unchanged

We present a study of the clustering and halo occupation distribution of BOSS CMASS galaxies in the redshift range 0.43 < z < 0.7 drawn from the Final SDSS-III Data Release. We compare the BOSS results with the predictions of a Halo Abundance Matching (HAM) clustering model that assigns galaxies to dark matter halos selected from the large BigMultiDark $N$-body simulation of a flat $\Lambda$CDM Planck cosmology. We compare the observational data with the simulated ones on a light-cone constructed from 20 subsequent outputs of the simulation. Observational effects such as incompleteness, geometry, veto masks and fiber collisions are included in the model, which reproduces within 1-$\sigma$ errors the observed monopole of the 2-point correlation function at all relevant scales: from the smallest scales, 0.5 $h^{-1}$ Mpc, up to scales beyond the Baryonic Acoustic Oscillation feature. This model also agrees remarkably well with the BOSS galaxy power spectrum (up to $k\sim1$ $h$ Mpc$^{-1}$), and the Three-point correlation function. The quadrupole of the correlation function presents some tensions with observations. We discuss possible causes that can explain this disagreement, including target selection effects. Overall, the standard HAM model describes remarkably well the clustering statistics of the CMASS sample. We compare the stellar to halo mass relation for the CMASS sample measured using weak lensing in the CFHT Stripe 82 Survey with the prediction of our clustering model, and find a good agreement within 1-$\sigma$. The BigMD-BOSS light-cone including properties of BOSS galaxies and halo properties is made publicly available., Comment: 16 pages, 16 figures. See companion papers that share the "The clustering of galaxies in the SDSS-III Baryon Oscillation Spectroscopic Survey:" title

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., Comment: 9 pages, 3 figures

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., Comment: matches published version

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., Comment: 13 pages, 8 figures. Published in the Astrophysical Journal

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)

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., Comment: Report from the "Dark Energy and CMB" working group for the American Physical Society's Division of Particles and Fields long-term planning exercise ("Snowmass")

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., Comment: 13 pages, 6 figures; new section shows agreement with bispectrum measurements in simulations

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., Comment: 16 pages, 10 figures. Accepted for publication in MNRAS (minor comments included to match accepted version). Public code available at http://maia.ice.cat/crocce/mptbreeze/

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., Comment: 25 pages, 14 figures; V2: significant revision including more details about halo exclusion and low-k noise. Conclusions unchanged

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., Comment: 21 pages, 18 figures, 2 tables. See companion papers that share the "The clustering of galaxies in the SDSS-III Baryon Oscillation Spectroscopic Survey:" title. Mocks and covariance matrices will be available at http://www.marcmanera.net/mocks when the paper is accepted and DR9 released. Minor edits. Amended author's names. Submitted to MNRAS

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), Comment: 26 pages, 9 figures. v2: improved appendix A

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., Comment: 21 pages, 14 figures. v2: corrections in response to referee report

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., Comment: 43 pages, 10 figures. v2: references added. 2LPT parallel code for generating non-local PNG initial conditions available at http://cosmo.nyu.edu/roman/2LPT

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., Comment: 18 pages, 17 figures. Submitted to ApJ

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., Comment: (v1): 23 pages text, 8 figures + appendix (v2): typos fixed, references added, accepted for publication in PRD

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., Comment: 27 pages, 21 figures. Updated to match accepted version. Published in ApJ

We analyze the line-of-sight baryonic acoustic feature in the two-point correlation function {\xi} 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 {\theta} < 3{\deg}, where {\theta} is the opening angle. This clustering excess appears much stronger than the feature in the better-measured monopole. A fiducial {\Lambda}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 {\xi} = 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)., Comment: 12 double-column pages of Text, 6 Figures, 2 Tables

We examine the correlation function \xi of the Sloan Digital Sky Survey (SDSS) Luminous Red Galaxy sample (LRG) at large scales (60

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., Comment: 19 pages, 13 figures. v2: corrected typos, added more simulations, better test of predictions in large mass regime. v3: minor changes, published version