Source author record

Martin Crocce

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

30works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

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

30 published item(s)

preprint2022arXiv

Clustering with general photo-$z$ uncertainties: Application to Baryon Acoustic Oscillations

Photometric data can be analyzed using the three-dimensional correlation function $ξ_{\rm p}$ to extract cosmological information via e.g., measurement of the Baryon Acoustic Oscillations (BAO). Previous studies modeled $ξ_{\rm p} $ assuming a Gaussian photo-$z$ approximation. In this work we improve the modeling by incorporating realistic photo-$z$ distribution. We show that the position of the BAO scale in $ξ_{\rm p}$ is determined by the photo-$z$ distribution and the Jacobian of the transformation. The latter diverges at the transverse scale of the separation $s_\perp $, and it explains why $ξ_{\rm p } $ traces the underlying correlation function at $s_\perp $, rather than $s$, when the photo-$z$ uncertainty $ σ_z / (1+ z) \gtrsim 0.02$. We also obtain the Gaussian covariance for $ξ_{\rm p}$. Due to photo-$z$ mixing, the covariance of $ξ_{\mathrm{p}}$ shows strong off-diagonal elements. The high correlation of the data causes some issues to the data fitting. Nonetheless, we find that either it can be solved by suppressing the largest eigenvalues of the covariance or it is not directly related to the BAO. We test our BAO fitting pipeline using a set of mock catalogs. The data set is dedicated for Dark Energy Survey Year 3 (DES Y3) BAO analyses and includes realistic photo-$z$ distributions. The theory template is in good agreement with mock measurement. Based on the DES Y3 mocks, $ξ_{\rm p}$ statistic is forecast to constrain the BAO shift parameter $α$ to be $1.001 \pm 0.023$, which is well consistent with the corresponding constraint derived from the angular correlation function measurements. Thus $ξ_{\rm p}$ offers a competitive alternative for the photometric data analyses.

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.

preprint2022arXiv

Modeling Intrinsic Galaxy Alignment in the MICE Simulation

The intrinsic alignment (IA) of galaxies is potentially a major limitation in deriving cosmological constraints from weak lensing surveys. In order to investigate this effect we assign intrinsic shapes and orientations to galaxies in the light-cone output of the MICE simulation, spanning $\sim5000\,{\rm deg}^2$ and reaching redshift $z=1.4$. This assignment is based on a 'semi-analytic' IA model that uses photometric properties of galaxies as well as the spin and shape of their host halos. Advancing on previous work, we include more realistic distributions of galaxy shapes and a luminosity dependent galaxy-halo alignment. The IA model parameters are calibrated against COSMOS and BOSS LOWZ observations. The null detection of IA in observations of blue galaxies is accounted for by setting random orientations for these objects. We compare the two-point alignment statistics measured in the simulation against predictions from the analytical IA models NLA and TATT over a wide range of scales, redshifts and luminosities for red and blue galaxies separately. We find that both models fit the measurements well at scales above $8\,h^{-1}{\rm Mpc}$, while TATT outperforms NLA at smaller scales. The IA parameters derived from our fits are in broad agreement with various observational constraints from red galaxies. Lastly, we build a realistic source sample, mimicking DES Year 3 observations and use it to predict the IA contamination to the observed shear statistics. We find this prediction to be within the measurement uncertainty, which might be a consequence of the random alignment of blue galaxies in the simulation.

preprint2022arXiv

New catalogue of dark-matter halo properties identified in MICE-GC -- I. Analysis of density profile distributions

Constraints on dark matter halo masses from weak gravitational lensing can be improved significantly by using additional information about the morphology of their density distribution, leading to tighter cosmological constraints derived from the halo mass function. This work is the first of two in which we investigate the accuracy of halo morphology and mass measurements in 2D and 3D. To this end, we determine several halo physical properties in the MICE-Grand Challenge dark matter only simulation. We present a public catalogue of these properties that includes density profiles and shape parameters measured in 2D and 3D, the halo centre at the peak of the 3D density distribution as well as the gravitational and kinetic energies and angular momentum vectors. The density profiles are computed using spherical and ellipsoidal radial bins, taking into account the halo shapes. We also provide halo concentrations and masses derived from fits to 2D and 3D density profiles using NFW and Einasto models for halos with more than $1000$ particles ($\gtrsim 3 \times 10^{13} h^{-1} M_{\odot}$). We find that the Einasto model provides better fits compared to NFW, regardless of the halo relaxation state and shape. The mass and concentration parameters of the 3D density profiles derived from fits to the 2D profiles are in general biased. Similar biases are obtained when constraining mass and concentrations using a weak-lensing stacking analysis. We show that these biases depend on the radial range and density profile model adopted in the fitting procedure, but not on the halo shape.

preprint2021arXiv

Testing one-loop galaxy bias: Power spectrum

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.

preprint2020arXiv

Dark Energy Survey Year 1 Results: Cosmological Constraints from Cluster Abundances and Weak Lensing

We perform a joint analysis of the counts and weak lensing signal of redMaPPer clusters selected from the Dark Energy Survey (DES) Year 1 dataset. Our analysis uses the same shear and source photometric redshifts estimates as were used in the DES combined probes analysis. Our analysis results in surprisingly low values for $S_8 =σ_8(Ω_{\rm m}/0.3)^{0.5}= 0.65\pm 0.04$, driven by a low matter density parameter, $Ω_{\rm m}=0.179^{+0.031}_{-0.038}$, with $σ_8-Ω_{\rm m}$ posteriors in $2.4σ$ tension with the DES Y1 3x2pt results, and in $5.6σ$ with the Planck CMB analysis. These results include the impact of post-unblinding changes to the analysis, which did not improve the level of consistency with other data sets compared to the results obtained at the unblinding. The fact that multiple cosmological probes (supernovae, baryon acoustic oscillations, cosmic shear, galaxy clustering and CMB anisotropies), and other galaxy cluster analyses all favor significantly higher matter densities suggests the presence of systematic errors in the data or an incomplete modeling of the relevant physics. Cross checks with X-ray and microwave data, as well as independent constraints on the observable--mass relation from SZ selected clusters, suggest that the discrepancy resides in our modeling of the weak lensing signal rather than the cluster abundance. Repeating our analysis using a higher richness threshold ($λ\ge 30$) significantly reduces the tension with other probes, and points to one or more richness-dependent effects not captured by our model.

preprint2020arXiv

KiDS-1000 Cosmology: Multi-probe weak gravitational lensing and spectroscopic galaxy clustering constraints

We present a joint cosmological analysis of weak gravitational lensing observations from the Kilo-Degree Survey (KiDS-1000), with redshift-space galaxy clustering observations from the Baryon Oscillation Spectroscopic Survey (BOSS), and galaxy-galaxy lensing observations from the overlap between KiDS-1000, BOSS and the spectroscopic 2-degree Field Lensing Survey (2dFLenS). This combination of large-scale structure probes breaks the degeneracies between cosmological parameters for individual observables, resulting in a constraint on the structure growth parameter $S_8=σ_8 \sqrt{Ω_{\rm m}/0.3} = 0.766^{+0.020}_{-0.014}$, that has the same overall precision as that reported by the full-sky cosmic microwave background observations from Planck. The recovered $S_8$ amplitude is low, however, by $8.3 \pm 2.6$ % relative to Planck. This result builds from a series of KiDS-1000 analyses where we validate our methodology with variable depth mock galaxy surveys, our lensing calibration with image simulations and null-tests, and our optical-to-near-infrared redshift calibration with multi-band mock catalogues and a spectroscopic-photometric clustering analysis. The systematic uncertainties identified by these analyses are folded through as nuisance parameters in our cosmological analysis. Inspecting the offset between the marginalised posterior distributions, we find that the $S_8$-difference with Planck is driven by a tension in the matter fluctuation amplitude parameter, $σ_8$. We quantify the level of agreement between the CMB and our large-scale structure constraints using a series of different metrics, finding differences with a significance ranging between $\sim\! 3\,σ$, when considering the offset in $S_{8}$, and $\sim\! 2\,σ$, when considering the full multi-dimensional parameter space.

preprint2016arXiv

A new method to measure galaxy bias by combining the density and weak lensing fields

We present a new method to measure the redshift-dependent galaxy bias by combining information from the galaxy density field and the weak lensing field. This method is based on Amara et al. (2012), where they use the galaxy density field to construct a bias-weighted convergence field kg. The main difference between Amara et al. (2012) and our new implementation is that here we present another way to measure galaxy bias using tomography instead of bias parameterizations. The correlation between kg and the true lensing field k allows us to measure galaxy bias using different zero-lag correlations, such as <kgk>/<kk> or <kgkg>/<kgk>. Our method measures the linear bias factor on linear scales under the assumption of no stochasticity between galaxies and matter. We use the MICE simulation to measure the linear galaxy bias for a flux-limited sample (i < 22.5) in tomographic redshift bins using this method. This paper is the first that studies the accuracy and systematic uncertainties associated with the implementation of the method, and the regime where it is consistent with the linear galaxy bias defined by projected 2-point correlation functions (2PCF). We find that our method is consistent with linear bias at the percent level for scales larger than 30 arcmin, while nonlinearities appear at smaller scales. This measurement is a good complement to other measurements of bias, since it does not depend strongly on sigma8 as the 2PCF measurements. We apply this method to the Dark Energy Survey Science Verification data in a follow-up paper.

preprint2016arXiv

Cosmology from large scale galaxy clustering and galaxy-galaxy lensing with Dark Energy Survey Science Verification data

We present cosmological constraints from the Dark Energy Survey (DES) using a combined analysis of angular clustering of red galaxies and their cross-correlation with weak gravitational lensing of background galaxies. We use a 139 square degree contiguous patch of DES data from the Science Verification (SV) period of observations. Using large scale measurements, we constrain the matter density of the Universe as Omega_m = 0.31 +/- 0.09 and the clustering amplitude of the matter power spectrum as sigma_8 = 0.74 +/- 0.13 after marginalizing over seven nuisance parameters and three additional cosmological parameters. This translates into S_8 = sigma_8(Omega_m/0.3)^{0.16} = 0.74 +/- 0.12 for our fiducial lens redshift bin at 0.35 <z< 0.5, while S_8 = 0.78 +/- 0.09 using two bins over the range 0.2 <z< 0.5. We study the robustness of the results under changes in the data vectors, modelling and systematics treatment, including photometric redshift and shear calibration uncertainties, and find consistency in the derived cosmological parameters. We show that our results are consistent with previous cosmological analyses from DES and other data sets and conclude with a joint analysis of DES angular clustering and galaxy-galaxy lensing with Planck CMB data, Baryon Accoustic Oscillations and Supernova type Ia measurements.

preprint2016arXiv

ICE-COLA: Towards fast and accurate synthetic galaxy catalogues optimizing a quasi $N$-body method

Next generation galaxy surveys demand the development of massive ensembles of galaxy mocks to model the observables and their covariances, what is computationally prohibitive using $N$-body simulations. COLA is a novel method designed to make this feasible by following an approximate dynamics but with up to 3 orders of magnitude speed-ups when compared to an exact $N$-body. In this paper we investigate the optimization of the code parameters in the compromise between computational cost and recovered accuracy in observables such as two-point clustering and halo abundance. We benchmark those observables with a state-of-the-art $N$-body run, the MICE Grand Challenge simulation (MICE-GC). We find that using 40 time steps linearly spaced since $z_i \sim 20$, and a force mesh resolution three times finer than that of the number of particles, yields a matter power spectrum within $1\%$ for $k \lesssim 1\,h {\rm Mpc}^{-1}$ and a halo mass function within $5\%$ of those in the $N$-body. In turn the halo bias is accurate within $2\%$ for $k \lesssim 0.7\,h {\rm Mpc}^{-1}$ whereas, in redshift space, the halo monopole and quadrupole are within $4\%$ for $k \lesssim 0.4\,h {\rm Mpc}^{-1}$. These results hold for a broad range in redshift ($0 < z < 1$) and for all halo mass bins investigated ($M > 10^{12.5} \, h^{-1} \, {\rm M_{\odot}}$). To bring accuracy in clustering to one percent level we study various methods that re-calibrate halo masses and/or velocities. We thus propose an optimized choice of COLA code parameters as a powerful tool to optimally exploit future galaxy surveys.

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.

preprint2016arXiv

The DESI Experiment Part I: Science,Targeting, and Survey Design

DESI (Dark Energy Spectroscopic Instrument) is a Stage IV ground-based dark energy experiment that will study baryon acoustic oscillations (BAO) and the growth of structure through redshift-space distortions with a wide-area galaxy and quasar redshift survey. To trace the underlying dark matter distribution, spectroscopic targets will be selected in four classes from imaging data. We will measure luminous red galaxies up to $z=1.0$. To probe the Universe out to even higher redshift, DESI will target bright [O II] emission line galaxies up to $z=1.7$. Quasars will be targeted both as direct tracers of the underlying dark matter distribution and, at higher redshifts ($ 2.1 < z < 3.5$), for the Ly-$α$ forest absorption features in their spectra, which will be used to trace the distribution of neutral hydrogen. When moonlight prevents efficient observations of the faint targets of the baseline survey, DESI will conduct a magnitude-limited Bright Galaxy Survey comprising approximately 10 million galaxies with a median $z\approx 0.2$. In total, more than 30 million galaxy and quasar redshifts will be obtained to measure the BAO feature and determine the matter power spectrum, including redshift space distortions.

preprint2016arXiv

The DESI Experiment Part II: Instrument Design

DESI (Dark Energy Spectropic Instrument) is a Stage IV ground-based dark energy experiment that will study baryon acoustic oscillations and the growth of structure through redshift-space distortions with a wide-area galaxy and quasar redshift survey. The DESI instrument is a robotically-actuated, fiber-fed spectrograph capable of taking up to 5,000 simultaneous spectra over a wavelength range from 360 nm to 980 nm. The fibers feed ten three-arm spectrographs with resolution $R= λ/Δλ$ between 2000 and 5500, depending on wavelength. The DESI instrument will be used to conduct a five-year survey designed to cover 14,000 deg$^2$. This powerful instrument will be installed at prime focus on the 4-m Mayall telescope in Kitt Peak, Arizona, along with a new optical corrector, which will provide a three-degree diameter field of view. The DESI collaboration will also deliver a spectroscopic pipeline and data management system to reduce and archive all data for eventual public use.

preprint2015arXiv

nIFTy Cosmology: Galaxy/halo mock catalogue comparison project on clustering statistics

We present a comparison of major methodologies of fast generating mock halo or galaxy catalogues. The comparison is done for two-point and the three-point clustering statistics. The reference catalogues are drawn from the BigMultiDark N-body simulation. Both friend-of-friends (including distinct halos only) and spherical overdensity (including distinct halos and subhalos) catalogs have been used with the typical number density of a large-volume galaxy surveys. We demonstrate that a proper biasing model is essential for reproducing the power spectrum at quasilinear and even smaller scales. With respect to various clustering statistics a methodology based on perturbation theory and a realistic biasing model leads to very good agreement with N-body simulations. However, for the quadrupole of the correlation function or the power spectrum, only the method based on semi-N-body simulation could reach high accuracy (1% level) at small scales, i.e., r<25 Mpc/h or k>0.15 h/Mpc. Full N-body solutions will remain indispensable to produce reference catalogues. Nevertheless, we have demonstrated that the far more efficient approximate solvers can reach a few percent accuracy in terms of clustering statistics at the scales interesting for the large-scale structure analysis after calibration with a few reference N-body calculations. This makes them useful for massive production aimed at covariance studies, to scan large parameter spaces, and to estimate uncertainties in data analysis techniques, such as baryon acoustic oscillation reconstruction, redshift distortion measurements, etc.

preprint2014arXiv

Measuring the growth of matter fluctuations with third-order galaxy correlations

Measurements of the linear growth factor $D$ at different redshifts $z$ are key to distinguish among cosmological models. One can estimate the derivative $dD(z)/d\ln(1+z)$ from redshift space measurements of the 3D anisotropic galaxy two-point correlation $ξ(z)$, but the degeneracy of its transverse (or projected) component with galaxy bias $b$, i.e. $ξ_{\perp}(z) \propto\ D^2(z) b^2(z)$, introduces large errors in the growth measurement. Here we present a comparison between two methods which break this degeneracy by combining second- and third-order statistics. One uses the shape of the reduced three-point correlation and the other a combination of third-order one- and two-point cumulants. These methods use the fact that, for Gaussian initial conditions and scales larger than $20$ $h^{-1}$Mpc, the reduced third-order matter correlations are independent of redshift (and therefore of the growth factor) while the third-order galaxy correlations depend on $b$. We use matter and halo catalogs from the MICE-GC simulation to test how well we can recover $b(z)$ and therefore $D(z)$ with these methods in 3D real space. We also present a new approach, which enables us to measure $D$ directly from the redshift evolution of second- and third-order galaxy correlations without the need of modelling matter correlations. For haloes with masses lower than $10^{14}$ $h^{-1}$M$_\odot$, we find $10%$ deviations between the different estimates of $D$, which are comparable to current observational errors. At higher masses we find larger differences that can probably be attributed to the breakdown of the bias model and non-Poissonian shot noise.

preprint2014arXiv

Precision cosmology in muddy waters: Cosmological constraints and N-body codes

Future large-scale structure surveys of the Universe will aim to constrain the cosmological model and the true nature of dark energy with unprecedented accuracy. In order for these surveys to achieve their designed goals, they will require predictions for the nonlinear matter power spectrum to sub-percent accuracy. Through the use of a large ensemble of cosmological N-body simulations, we demonstrate that if we do not understand the uncertainties associated with simulating structure formation, i.e. knowledge of the `true' simulation parameters, and simply seek to marginalize over them, then the constraining power of such future surveys can be significantly reduced. However, for the parameters {n_s, h, Om_b, Om_m}, this effect can be largely mitigated by adding the information from a CMB experiment, like Planck. In contrast, for the amplitude of fluctuations sigma8 and the time-evolving equation of state of dark energy {w_0, w_a}, the mitigation is mild. On marginalizing over the simulation parameters, we find that the dark-energy figure of merit can be degraded by ~2. This is likely an optimistic assessment, since we do not take into account other important simulation parameters. A caveat is our assumption that the Hessian of the likelihood function does not vary significantly when moving from our adopted to the 'true' simulation parameter set. This paper therefore provides strong motivation for rigorous convergence testing of N-body codes to meet the future challenges of precision cosmology.

preprint2014arXiv

Redshift-space distortions from the cross-correlation of photometric populations

Several papers have recently highlighted the possibility of measuring redshift space distortions from angular auto-correlations of galaxies in photometric redshift bins. In this work we extend this idea to include as observables the cross-correlations between redshift bins, as an additional way of measuring radial information. We show that this extra information allows to reduce the recovered error in the growth rate index γby a factor of ~2. Although the final error in γdepends on the bias and the mean photometric accuracy of the galaxy sample, the improvement from adding cross-correlations is robust in different settings. Another factor of 2-3 improvement in the determination of γcan be achieved by considering two galaxy populations over the same photometric sky area but with different biases. This additional gain is shown to be much larger than the one from the same populations when observed over different areas of the sky (with twice the combined area). The total improvement of ~5 implies that a photometric survey such as the Dark Energy Survey should be able to recover γat the 5-10% from the angular clustering in linear scales of two different tracers. It can also constrain the evolution of f(z)xσ_8(z) in few bins beyond z~0.8-0.9 at the 10-15% level per-bin, compatible with recent constrains from lower-z spectroscopic surveys. We also show how further improvement can be achieved by reducing the photometric redshift error.

preprint2014arXiv

The 6dF Galaxy Velocity Survey: Cosmological constraints from the velocity power spectrum

We present scale-dependent measurements of the normalised growth rate of structure $fσ_{8}(k, z=0)$ using only the peculiar motions of galaxies. We use data from the 6-degree Field Galaxy Survey velocity sample (6dFGSv) together with a newly-compiled sample of low-redshift $(z < 0.07)$ type Ia supernovae. We constrain the growth rate in a series of $Δk \sim 0.03 h{\rm Mpc^{-1}}$ bins to $\sim35\%$ precision, including a measurement on scales $>300 h^{-1}{\rm Mpc}$, which represents one of the largest-scale growth rate measurement to date. We find no evidence for a scale dependence in the growth rate, or any statistically significant variation from the growth rate as predicted by the {\it Planck} cosmology. Bringing all the scales together, we determine the normalised growth rate at $z=0$ to $\sim15\%$ in a manner {\it independent} of galaxy bias and in excellent agreement with the constraint from the measurements of redshift-space distortions from 6dFGS. We pay particular attention to systematic errors. We point out that the intrinsic scatter present in Fundamental-Plane and Tully-Fisher relations is only Gaussian in logarithmic distance units; wrongly assuming it is Gaussian in linear (velocity) units can bias cosmological constraints. We also analytically marginalise over zero-point errors in distance indicators, validate the accuracy of all our constraints using numerical simulations, and demonstrate how to combine different (correlated) velocity surveys using a matrix `hyper-parameter' analysis. Current and forthcoming peculiar velocity surveys will allow us to understand in detail the growth of structure in the low-redshift universe, providing strong constraints on the nature of dark energy.

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.

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.

preprint2012arXiv

Recovering 3D clustering information with angular correlations

We study how to recover the full 3D clustering information of P(\vec{k},z), including redshift space distortions (RSD), from 2D tomography using the angular auto and cross spectra of different redshift bins C_\ell(z,z'). We focus on quasilinear scales where the minimum scale λ_{min} or corresponding maximum wavenumber k_{max}= 2π/λ_{min} is targeted to be between k_{max}={0.05-0.2} h/Mpc. For spectroscopic surveys, we find that we can recover the full 3D clustering information when the redshift bin width Δz used in the 2D tomography is similar to the targeted minimum scale, i.e. Δz ~ {0.6-0.8} λ_{min} H(z)/c which corresponds to Δz ~ 0.01-0.05 for z<1. This value of Δz is optimal in the sense that larger values of Δz lose information, while smaller values violate our minimum scale requirement. For a narrow-band photometric survey, with photo-z error σ_z=0.004, we find almost identical results to the spectroscopic survey because the photo-z error is smaller than the optimal bin width σ_z<Δz. For a typical broad-band photometric survey with σ_z=0.1, we have that σ_z>Δz and most radial information is intrinsically lost. The remaining information can be recovered from the 2D tomography if we use Δz ~ 2σ_z. While 3D and 2D analysis are shown here to be equivalent, the advantage of using angular positions and redshifts is that we do not need a fiducial cosmology to convert to 3D coordinates. This avoids assumptions and marginalization over the fiducial model. In addition, it becomes straight forward to combine RSD, clustering and weak lensing in 2D space.

preprint2011arXiv

Clustering of photometric luminous red galaxies I : Growth of Structure and Baryon Acoustic Feature

The possibility of measuring redshift space (RSD) distortions using photometric data have been recently highlighted. This effect complements and significantly alters the detectability of baryon acoustic oscillations (BAO) in photometric surveys. In this paper we present measurements of the angular correlation function of luminous red galaxies (LRGs) in the photometric catalog of the final data release (DR7) of the Sloan Digital Sky Survey II (SDSS). The sample compromise ~ 1.5 x 10^6 LRGs distributed in 0.45 < z < 0.65, with a characteristic photometric error of ~ 0.05. Our measured correlation centered at z=0.55 is in very good agreement with predictions from standard LCDM in a broad range of angular scales, $0.5^\circ < θ< 6^\circ$. We find that the growth of structure can indeed be robustly measured, with errors matching expectations. The velocity growth rate is recovered as $f σ_8 = 0.53 \pm 0.42$ when no prior is imposed on the growth factor and the background geometry follows a LCDM model with WMAP7+SNIa priors. This is compatible with the corresponding General Relativity (GR) prediction $f σ_8 = 0.45$ for our fiducial cosmology. If we adopt a parametrization such that $f=Ω^γ_m(z)$, with $γ\approx 0.55$ in GR, and combine our $fσ_8$ measurement with the corresponding ones from spectroscopic LRGs at lower redshifts we obtain $γ=0.54 \pm 0.17$. In addition we find evidence for the presence of the baryon acoustic feature matching the amplitude, location and shape of LCDM predictions. The photometric BAO feature is detected with 98 % confidence level at z=0.55.

preprint2011arXiv

Measuring Redshift-Space Distortions using Photometric Surveys

We outline how redshift-space distortions (RSD) can be measured from the angular correlation function w(θ), of galaxies selected from photometric surveys. The natural degeneracy between RSD and galaxy bias can be minimized by comparing results from bins with top-hat galaxy selection in redshift, and bins based on the radial position of galaxy pair centres. This comparison can also be used to test the accuracy of the photometric redshifts. The presence of RSD will be clearly detectable with the next generation of photometric redshift surveys. We show that the Dark Energy Survey (DES) will be able to measure f(z)σ_8(z) to a 1σ accuracy of (17 {\times} b)%, using galaxies drawn from a single narrow redshift slice centered at z = 1. Here b is the linear bias, and f is the logarithmic rate of change of the linear growth rate with respect to the scale factor. Extending to measurements of w(θ) for a series of bins of width 0.02(1 + z) over 0.5 < z < 1.4 will measure γ to a 1σ accuracy of 25%, given the model f = Ω_m(z)^γ, and assuming a linear bias model that evolves such that b = 0.5 + z (and fixing other cosmological parameters). The accuracy of our analytic predictions is confirmed using mock catalogs drawn from simulations conducted by the MICE collaboration.

preprint2011arXiv

Modeling the angular correlation function and its full covariance in Photometric Galaxy Surveys

Near future cosmology will see the advent of wide area photometric galaxy surveys, like the Dark Energy Survey (DES), that extent to high redshifts (z ~ 1 - 2) but with poor radial distance resolution. In such cases splitting the data into redshift bins and using the angular correlation function $w(θ)$, or the $C_{\ell}$ power spectrum, will become the standard approach to extract cosmological information or to study the nature of dark energy through the Baryon Acoustic Oscillations (BAO) probe. In this work we present a detailed model for $w(θ)$ at large scales as a function of redshift and bin width, including all relevant effects, namely nonlinear gravitational clustering, bias, redshift space distortions and photo-z uncertainties. We also present a model for the full covariance matrix characterizing the angular correlation measurements, that takes into account the same effects as for $w(θ)$ and also the possibility of a shot-noise component and partial sky coverage. Provided with a large volume N-body simulation from the MICE collaboration we built several ensembles of mock redshift bins with a sky coverage and depth typical of forthcoming photometric surveys. The model for the angular correlation and the one for the covariance matrix agree remarkably well with the mock measurements in all configurations. The prospects for a full shape analysis of $w(θ)$ at BAO scales in forthcoming photometric surveys such as DES are thus very encouraging.

preprint2011arXiv

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

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

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.

preprint2010arXiv

Multi-Point Propagators for Non-Gaussian Initial Conditions

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

preprint2010arXiv

On the radiative efficiencies, Eddington ratios, and duty cycles of luminous high-redshift quasars

We investigate the characteristic radiative efficiency ε, Eddington ratio λ, and duty cycle P_0 of high-redshift active galactic nuclei (AGN), drawing on measurements of the AGN luminosity function at z=3-6 and, especially, on recent measurements of quasar clustering at z=3-4.5 from the Sloan Digital Sky Survey. The free parameters of our models are ε, λ, and the normalization, scatter, and redshift evolution of the relation between black hole mass \mbh and halo virial velocity V_vir. We compute the luminosity function from the implied growth of the black hole mass function and the quasar correlation length from the bias of the host halos. We test our adopted formulae for the halo mass function and halo bias against measurements from the large N-body simulation developed by the MICE collaboration. The strong clustering of AGNs observed at z=3 and, especially, at z=4 implies that massive black holes reside in rare, massive dark matter halos. Reproducing the observed luminosity function then requires high efficiency εand/or low Eddington ratio λ, with a lower limit (based on 2σagreement with the measured z=4 correlation length) ε> 0.7λ/(1+0.7λ), implying ε> 0.17 for λ> 0.25. Successful models predict high duty cycles, P_0~0.2, 0.5, and 0.9 at z=3.1, 4.5 and 6, respectively, and they require that the fraction of halo baryons locked in the central black hole is much larger than the locally observed value. The rapid drop in the abundance of the massive and rare host halos at z>7 implies a proportionally rapid decline in the number density of luminous quasars, much stronger than simple extrapolations of the z=3-6 luminosity function would predict. (abridged)

preprint2010arXiv

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

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

preprint2009arXiv

Simulating the Universe with MICE: The abundance of massive clusters

We introduce a new set of large N-body runs, the MICE simulations, that provide a unique combination of very large cosmological volumes with good mass resolution. They follow the gravitational evolution of ~ 8.5 billion particles (2048^3) in volumes covering up to 450 (Gpc/h)^3. Our main goal is to accurately model and calibrate basic cosmological probes that will be used by upcoming astronomical surveys. Here we take advantage of the very large volumes of MICE to make a robust sampling of the high-mass tail of the halo mass function (MF). We discuss and avoid possible systematic effects in our study, and do a detailed analysis of different error estimators. We find that available fits to the local abundance of halos (Warren et al. (2006)) match well the abundance in MICE up to M ~ 10^{14}\Msun, but significantly deviate for larger masses, underestimating the mass function by 10% (30%) at M = 3.16 x 10^{14}\Msun (10^{15}\Msun). Similarly, the widely used Sheth & Tormen (1999) fit, if extrapolated to high redshift assuming universality, leads to an underestimation of the cluster abundance by 30%, 20% and 15% at z=0, 0.5, 1 for M ~ [7 - 2.5 - 0.8] x 10^{14}\Msun respectively ($ν= δ_c/σ~ 3$). We provide a re-calibration of the halo MF valid over 5 orders of magnitude in mass, 10^{10} < M/(\Msun) < 10^{15}, that accurately describes its redshift evolution up to z=1. We explore the impact of this re-calibration on the determination of dark-energy, and conclude that using available fits may systematically bias the estimate of w by as much as 50% for medium-depth (z <= 1) surveys. MICE halo catalogues are publicly available at http://www.ice.cat/mice