BZPEERPapers, people, evidence and research paths
Menu

Discover

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

Research tools

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

Network

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

Opportunities

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

Account

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

Source author record

Guilhem Lavaux

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

ResearcherUnclaimed source record
astro-ph.COastro-ph.GAastro-ph.IMgr-qcApplicationsastro-ph.HEphysics.data-an

Catalog footprint

What is connected

39works
7topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

39 published item(s)

preprint2024arXiv

Towards Accurate Field-Level Inference of Massive Cosmic Structures

We investigate the accuracy requirements for field-level inference of cluster and void masses using data from galaxy surveys. We introduce a two-step framework that takes advantage of the fact that cluster masses are determined by flows on larger scales than the clusters themselves. First, we determine the integration accuracy required to perform field-level inference of cosmic initial conditions on these large scales, by fitting to late-time galaxy counts using the Bayesian Origin Reconstruction from Galaxies (BORG) algorithm. A 20-step COLA integrator is able to accurately describe the density field surrounding the most massive clusters in the Local Super-Volume ($<135\,h^{-1}\,\mathrm{Mpc}$), but does not by itself lead to converged virial mass estimates. Therefore we carry out `posterior resimulations', using full $N$-body dynamics while sampling from the inferred initial conditions, and thereby obtain estimates of masses for nearby massive clusters. We show that these are in broad agreement with existing estimates, and find that mass functions in the Local Super-Volume are compatible with $Λ$CDM.

0 citations0 reviewsNo code linkedOpen access
preprint2023arXiv

Higher-order statistics of the large-scale structure from photometric redshifts

The large-scale structure is a major source of cosmological information. However, next-generation photometric galaxy surveys will only provide a distorted view of cosmic structures due to large redshift uncertainties. To address the need for accurate reconstructions of the large-scale structure in presence of photometric uncertainties, we present a framework that constrains the three-dimensional dark matter density jointly with galaxy photometric redshift probability density functions (PDFs), exploiting information from galaxy clustering. Our forward model provides Markov Chain Monte Carlo realizations of the primordial and present-day dark matter density, inferred jointly from data. Our method goes beyond 2-point statistics via field-level inference. It accounts for all observational uncertainties and the survey geometry. We showcase our method using mock catalogs that emulate next-generation surveys with a worst-case redshift uncertainty, equivalent to ${\sim}300$ Mpc. On scales $150$ Mpc, we improve the cross-correlation of the photometric galaxy positions with the ground truth from $28\%$ to $86\%$. The improvement is significant down to $13$ Mpc. On scales $150$ Mpc, we achieve a cross-correlation of $80-90\%$ with the ground truth for the dark matter density, radial peculiar velocities, tidal shear and gravitational potential.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Field-level inference of galaxy intrinsic alignment from the SDSS-III BOSS survey

As a large-scale overdensity collapses, it affects the orientation and shape of galaxies that form, by exerting tidal shear along their axes. Therefore, the shapes of elliptical galaxies align with the tidal field of cosmic structures. This intrinsic alignment provides insights into galaxy formation and the primordial universe, complements late-time cosmological probes and constitutes a significant systematic effect for weak gravitational lensing observations. In the present study, we provide constraints on the linear alignment model using a fully Bayesian field-level approach, using galaxy shape measurements from the SDSS-III BOSS LOWZ sample and three-dimensional tidal fields constrained with the LOWZ and CMASS galaxy samples of the SDSS-III BOSS survey. We find 4$σ$ evidence of intrinsic alignment, with an amplitude of $A_I=2.9 \pm 0.7$ at 20$h^{-1}\;\mathrm{Mpc}$.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

SIBELIUS-DARK: a galaxy catalogue of the Local Volume from a constrained realisation simulation

We present SIBELIUS-DARK, a constrained realisation simulation of the local volume to a distance of 200~Mpc from the Milky Way. SIBELIUS-DARK is the first study of the \textit{Simulations Beyond The Local Universe} (SIBELIUS) project, which has the goal of embedding a model Local Group-like system within the correct cosmic environment. The simulation is dark-matter-only, with the galaxy population calculated using the semi-analytic model of galaxy formation, GALFORM. We demonstrate that the large-scale structure that emerges from the SIBELIUS constrained initial conditions matches well the observational data. The inferred galaxy population of SIBELIUS-DARK also match well the observational data, both statistically for the whole volume and on an object-by-object basis for the most massive clusters. For example, the $K$-band number counts across the whole sky, and when divided between the northern and southern Galactic hemispheres, are well reproduced by SIBELIUS-DARK. We find that the local volume is somewhat unusual in the wider context of $Λ$CDM: it contains an abnormally high number of supermassive clusters, as well as an overall large-scale underdensity at the level of $\approx 5$\% relative to the cosmic mean. However, whilst rare, the extent of these peculiarities does not significantly challenge the $Λ$CDM model. SIBELIUS-DARK is the most comprehensive constrained realisation simulation of the local volume to date, and with this paper we publicly release the halo and galaxy catalogues at $z=0$, which we hope will be useful to the wider astronomy community.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The Milky Way's plane of satellites: consistent with $Λ$CDM

The "plane of satellites problem" describes the arrangement of the Milky Way's 11 brightest satellite galaxies in a remarkably thin plane, possibly supported by rotation. This is in apparent contradiction to the standard cosmological model, wherein the Galaxy is surrounded by a dispersion-supported dark matter halo. Here, we show that the reported exceptional anisotropy of the satellite system is strongly contingent on a lopsided radial distribution, which earlier simulations have failed to reproduce, combined with the close but fleeting conjunction of the two most distant satellites, Leo I and Leo II. Using Gaia proper motions, we show that the orbital pole alignment is much more common than previously reported, and reveal the plane of satellites to be transient rather than rotationally supported. Comparing to new simulations, where such short-lived planes are common, we find the Milky Way satellites to be compatible with standard model expectations.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Bayesian forward modelling of cosmic shear data

We present a Bayesian hierarchical modelling approach to infer the cosmic matter density field, and the lensing and the matter power spectra, from cosmic shear data. This method uses a physical model of cosmic structure formation to infer physically plausible cosmic structures, which accounts for the non-Gaussian features of the gravitationally evolved matter distribution and light-cone effects. We test and validate our framework with realistic simulated shear data, demonstrating that the method recovers the unbiased matter distribution and the correct lensing and matter power spectrum. While the cosmology is fixed in this test, and the method employs a prior power spectrum, we demonstrate that the lensing results are sensitive to the true power spectrum when this differs from the prior. In this case, the density field samples are generated with a power spectrum that deviates from the prior, and the method recovers the true lensing power spectrum. The method also recovers the matter power spectrum across the sky, but as currently implemented, it cannot determine the radial power since isotropy is not imposed. In summary, our method provides physically plausible inference of the dark matter distribution from cosmic shear data, allowing us to extract information beyond the two-point statistics and exploiting the full information content of the cosmological fields.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Detectability of large correlation length inflationary magnetic field with Cherenkov telescopes

Magnetic fields occupying the voids of the large scale structure may be a relic from the Early Universe originating from either Inflation or from cosmological phase transitions. We explore the possibility of identifying the inflationary origin of the void magnetic fields and measuring its parameters with gamma-ray astronomy methods. The large correlation length inflationary field is expected to impose a characteristic asymmetry of extended gamma-ray emission that is correlated between different sources on the sky. We show that a set of nearby blazars for which the extended emission is observable in the 0.1-1 TeV band with CTA can be used for the test of inflationary origin of the void magnetic fields.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Optimal machine-driven acquisition of future cosmological data

We present maps classifying regions of the sky according to their information gain potential as quantified by the Fisher information. These maps can guide the optimal retrieval of relevant physical information with targeted cosmological searches. Specifically, we calculate the response of observed cosmic structures to perturbative changes in the cosmological model and chart their respective contributions to the Fisher information. Our physical forward modeling machinery transcends the limitations of contemporary analyses based on statistical summaries to yield detailed characterizations of individual 3D structures. We demonstrate this using galaxy counts data and showcase the potential of our approach by studying the information gain of the Coma cluster. We find that regions in the vicinity of the filaments and cluster core, where mass accretion ensues from gravitational infall, are the most informative about our physical model of structure formation in the Universe. Hence, collecting data in those regions would be most optimal for testing our model predictions. The results presented in this work are the first of their kind and elucidate the inhomogeneous distribution of cosmological information in the Universe. This study paves a new way forward to perform efficient targeted searches for the fundamental physics of the Universe, where search strategies are progressively refined with new cosmological data sets within an active learning framework.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Cosmology Inference from Biased Tracers using the EFT-based Likelihood

The effective-field-theory (EFT) approach to the clustering of galaxies and other biased tracers allows for an isolation of the cosmological information that is protected by symmetries, in particular the equivalence principle, and thus is robust to the complicated dynamics of dark matter, gas, and stars on small scales. All existing implementations proceed by making predictions for the lowest-order $n$-point functions of biased tracers, as well as their covariance, and comparing with measurements. Recently, we presented an EFT-based expression for the conditional probability of the density field of a biased tracer given the matter density field, which in principle combines information from arbitrarily high order $n$-point functions. Here, we report results based on this likelihood by applying it to halo catalogs in real space, specifically on the inference of the power spectrum normalization $σ_8$. We include bias terms up to second order as well as the leading higher-derivative term. For a cutoff value of $Λ= 0.1 h\,{\rm Mpc}^{-1}$, we recover the ground-truth value of $σ_8$ to within 95% CL for different halo samples and redshifts. We discuss possible sources for the remaining systematic bias in $σ_8$ as well as future developments.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Setting the Stage: Structures from Gaussian Random Fields

We study structure formation in a set of cosmological simulations to uncover the scales in the initial density field that gave rise to the formation of present-day structures. Our simulations share a common primordial power spectrum (here Lambda-CDM), but the introduction of hierarchical variations of the phase information allows us to systematically study the scales that determine the formation of structure at later times. We consider the variance in z=0 statistics such as the matter power spectrum and halo mass function. We also define a criterion for the existence of individual haloes across simulations, and determine what scales in the initial density field contain sufficient information for the non-linear formation of unique haloes. We study how the characteristics of individual haloes such as the mass and concentration, as well as the position and velocity, are affected by variations on different scales, and give scaling relations for haloes of different mass. Finally, we use the example of a cluster-mass halo to show how our hierarchical parametrisation of the initial density field can be used to create variants of particular objects. With properties such as mass, concentration, kinematics and substructure of haloes set on distinct and well-determined scales, and its unique ability to introduce variations localised in real space, our method is a powerful tool to study structure formation in cosmological simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Velocity correction for Hubble constant measurements from standard sirens

Gravitational wave (GW) sources are an excellent probe of the luminosity distance and offer a novel measure of the Hubble constant, $H_0$. This estimation of $H_0$ from standard sirens requires an accurate estimation of the cosmological redshift of the host galaxy of the GW source, after correcting for its peculiar velocity. Absence of an accurate peculiar velocity correction affects both the precision and accuracy of the measurement of $H_0$, particularly for nearby sources. We propose a framework to incorporate such a peculiar velocity correction for GW sources. A first implementation of our method to the event GW170817 combined with the Very Large Baseline Interferometry (VLBI) observation leads to a revised value of $H_0= 68.3^{+ 4.6}_{-4.5}$ km/s/Mpc. While this revision is minor, it demonstrates that our method makes it possible for obtaining an unbiased and accurate measurements of $H_0$ at the precision required for the standard siren cosmology.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Neural physical engines for inferring the halo mass distribution function

An ambitious goal in cosmology is to forward-model the observed distribution of galaxies in the nearby Universe today from the initial conditions of large-scale structures. For practical reasons, the spatial resolution at which this can be done is necessarily limited. Consequently, one needs a mapping between the density of dark matter averaged over ~Mpc scales, and the distribution of dark matter halos (used as a proxy for galaxies) in the same region. Here we demonstrate a method for determining the halo mass distribution function by learning the tracer bias between density fields and halo catalogues using a neural bias model. The method is based on the Bayesian analysis of simple, physically motivated, neural network-like architectures, which we denote as neural physical engines, and neural density estimation. As a result, we are able to sample the initial phases of the dark matter density field whilst inferring the parameters describing the halo mass distribution function, providing a fully Bayesian interpretation of both the initial dark matter density distribution and the neural bias model. We successfully run an upgraded BORG inference using our new likelihood and neural bias model with halo catalogues derived from full N-body simulations. We notice orders of magnitude improvement in modelling compared to classical biasing techniques.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Comparing cosmic web classifiers using information theory

We introduce a decision scheme for optimally choosing a classifier, which segments the cosmic web into different structure types (voids, sheets, filaments, and clusters). Our framework, based on information theory, accounts for the design aims of different classes of possible applications: (i) parameter inference, (ii) model selection, and (iii) prediction of new observations. As an illustration, we use cosmographic maps of web-types in the Sloan Digital Sky Survey to assess the relative performance of the classifiers T-web, DIVA and ORIGAMI for: (i) analyzing the morphology of the cosmic web, (ii) discriminating dark energy models, and (iii) predicting galaxy colors. Our study substantiates a data-supported connection between cosmic web analysis and information theory, and paves the path towards principled design of analysis procedures for the next generation of galaxy surveys. We have made the cosmic web maps, galaxy catalog, and analysis scripts used in this work publicly available.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Constraints on Cosmology and Gravity from the Dynamics of Voids

The Universe is mostly composed of large and relatively empty domains known as cosmic voids, whereas its matter content is predominantly distributed along their boundaries. The remaining material inside them, either dark or luminous matter, is attracted to these boundaries and causes voids to expand faster and to grow emptier over time. Using the distribution of galaxies centered on voids identified in the Sloan Digital Sky Survey and adopting minimal assumptions on the statistical motion of these galaxies, we constrain the average matter content $Ω_\mathrm{m}=0.281\pm0.031$ in the Universe today, as well as the linear growth rate of structure $f/b=0.417\pm0.089$ at median redshift $\bar{z}=0.57$, where $b$ is the galaxy bias ($68\%$ C.L.). These values originate from a percent-level measurement of the anisotropic distortion in the void-galaxy cross-correlation function, $\varepsilon = 1.003\pm0.012$, and are robust to consistency tests with bootstraps of the data and simulated mock catalogs within an additional systematic uncertainty of half that size. They surpass (and are complementary to) existing constraints by unlocking cosmological information on smaller scales through an accurate model of nonlinear clustering and dynamics in void environments. As such, our analysis furnishes a powerful probe of deviations from Einstein's general relativity in the low-density regime which has largely remained untested so far. We find no evidence for such deviations in the data at hand.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Bayesian inference of CMB gravitational lensing

The Planck satellite, along with several ground based telescopes, have mapped the cosmic microwave background (CMB) at sufficient resolution and signal-to-noise so as to allow a detection of the subtle distortions due to the gravitational influence of the intervening matter distribution. A natural modeling approach is to write a Bayesian hierarchical model for the lensed CMB in terms of the unlensed CMB and the lensing potential. So far there has been no feasible algorithm for inferring the posterior distribution of the lensing potential from the lensed CMB map. We propose a solution that allows efficient Markov Chain Monte Carlo sampling from the joint posterior of the lensing potential and the unlensed CMB map using the Hamiltonian Monte Carlo technique. The main conceptual step in the solution is a re-parameterization of CMB lensing in terms of the lensed CMB and the "inverse lensing" potential. We demonstrate a fast implementation on simulated data including noise and a sky cut, that uses a further acceleration based on a very mild approximation of the inverse lensing potential. We find that the resulting Markov Chain has short correlation lengths and excellent convergence properties, making it promising for application to high resolution CMB data sets of the future.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Cosmological parameters from the comparison of peculiar velocities with predictions from the 2M++ density field

Peculiar velocity measurements are the only tool available in the low-redshift Universe for mapping the large-scale distribution of matter and can thus be used to constrain cosmology. Using redshifts from the 2M++ redshift compilation, we reconstruct the density of galaxies within 200 Mpc/h, allowing for the first time good sampling of important superclusters such as the Shapley Concentration. We compare the predicted peculiar velocities from 2M++ to Tully-Fisher and SNe peculiar velocities. We find a value of $β^* \equiv Ω_{\rm{m}}^{0.55}/b^* = 0.431 \pm 0.021$, suggesting $Ω_{\rm{m}}^{0.55}σ_{\rm{8,lin}} = 0.401 \pm 0.024$, in good agreement with other probes. The predicted peculiar velocity of the Local Group arising from the 2M++ volume alone is $540 \pm 40$ km/s, towards $l = 268 \pm 4, b= 38 \pm 6$, only $10^\circ$ out of alignment with the Cosmic Microwave Background dipole. To account for velocity contributions arising from sources outside the 2M++ volume, we fit simultaneously for $β^*$ and an external bulk flow in our analysis. We find that an external bulk flow is preferred at the 5.1$σ$ level, and the best fit has a velocity of $159\pm23$ km/s towards $l=304 \pm 11, b = 6 \pm 13$. Finally, the predicted bulk flow of a 50 Mpc/h Gaussian-weighted volume centred on the Local Group is $230 \pm 30$ km/s, in the direction $l=293\pm 8, b = 14 \pm 10$, in agreement with predictions from $Λ$CDM.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Probabilistic cartography of the large-scale structure

The BORG algorithm is an inference engine that derives the initial conditions given a cosmological model and galaxy survey data, and produces physical reconstructions of the underlying large-scale structure by assimilating the data into the model. We present the application of BORG to real galaxy catalogs and describe the primordial and late-time large-scale structure in the considered volumes. We then show how these results can be used for building various probabilistic maps of the large-scale structure, with rigorous propagation of uncertainties. In particular, we study dynamic cosmic web elements and secondary effects in the cosmic microwave background.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Probing cosmology and gravity with redshift-space distortions around voids

Cosmic voids in the large-scale structure of the Universe affect the peculiar motions of objects in their vicinity. Although these motions are difficult to observe directly, the clustering pattern of their surrounding tracers in redshift space is influenced in a unique way. This allows to investigate the interplay between densities and velocities around voids, which is solely dictated by the laws of gravity. With the help of $N$-body simulations and derived mock-galaxy catalogs we calculate the average density fluctuations around voids identified with a watershed algorithm in redshift space and compare the results with the expectation from general relativity and the $Λ$CDM model. We find linear theory to work remarkably well in describing the dynamics of voids. Adopting a Bayesian inference framework, we explore the full posterior of our model parameters and forecast the achievable accuracy on measurements of the growth rate of structure and the geometric distortion through the Alcock-Paczynski effect. Systematic errors in the latter are reduced from $\sim15\%$ to $\sim5\%$ when peculiar velocities are taken into account. The relative parameter uncertainties in galaxy surveys with number densities comparable to the SDSS MAIN (CMASS) sample probing a volume of $1h^{-3}{\rm Gpc}^3$ yield $σ_{f/b}\left/(f/b)\right.\sim2\%$ ($20\%$) and $σ_{D_AH}/D_AH\sim0.2\%$ ($2\%$), respectively. At this level of precision the linear-theory model becomes systematics dominated, with parameter biases that fall beyond these values. Nevertheless, the presented method is highly model independent; its viability lies in the underlying assumption of statistical isotropy of the Universe.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Real-space density profile reconstruction of stacked voids

We present a non-parametric, model-independent method to reconstruct the spherical density profiles of void stacks in real space, without redshift-space distortions. Our method uses the expected spherical symmetry of stacked voids to build the shape of the spherical density profile of cosmic voids in real space without any assumption about the cosmological model. We test the reconstruction algorithm with both a toy model and a full dark matter simulation. We present the result for the simulation: the reconstruction of the spherical density profile for a simulated stacked void in real space. We also present a first application of the algorithm to reconstruct real cosmic void stacks density profiles in real space from the Sloan Digital Sky Survey (Sutter et al. 2012b). We discuss capabilities of the algorithm and possible future improvements. Reconstructed density profiles from real voids open the way to the study of the spherically averaged dynamical structure of voids.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Unmasking the Masked Universe: the 2M++ catalogue through Bayesian eyes

This work describes a full Bayesian analysis of the Nearby Universe as traced by galaxies of the 2M++ survey. The analysis is run in two sequential steps. The first step self-consistently derives the luminosity dependent galaxy biases, the power-spectrum of matter fluctuations and matter density fields within a Gaussian statistic approximation. The second step makes a detailed analysis of the three dimensional Large Scale Structures, assuming a fixed bias model and a fixed cosmology. This second step allows for the reconstruction of both the final density field and the initial conditions at z=1000 assuming a fixed bias model. From these, we derive fields that self-consistently extrapolate the observed large scale structures. We give two examples of these extrapolation and their utility for the detection of structures: the visibility of the Sloan Great Wall, and the detection and characterization of the Local Void using DIVA, a Lagrangian based technique to classify structures.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Cosmology with Void-Galaxy Correlations

Galaxy bias, the unknown relationship between the clustering of galaxies and the underlying dark matter density field is a major hurdle for cosmological inference from large-scale structure. While traditional analyses focus on the absolute clustering amplitude of high-density regions mapped out by galaxy surveys, we propose a relative measurement that compares those to the underdense regions, cosmic voids. On the basis of realistic mock catalogs we demonstrate that cross correlating galaxies and voids opens up the possibility to calibrate galaxy bias and to define a static ruler thanks to the observable geometric nature of voids. We illustrate how the clustering of voids is related to mass compensation and show that volume-exclusion significantly reduces the degree of stochasticity in their spatial distribution. Extracting the spherically averaged distribution of galaxies inside voids from their cross correlations reveals a remarkable concordance with the mass-density profile of voids.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Matrix-free Large Scale Bayesian inference in cosmology

In this work we propose a new matrix-free implementation of the Wiener sampler which is traditionally applied to high dimensional analysis when signal covariances are unknown. Specifically, the proposed method addresses the problem of jointly inferring a high dimensional signal and its corresponding covariance matrix from a set of observations. Our method implements a Gibbs sampling adaptation of the previously presented messenger approach, permitting to cast the complex multivariate inference problem into a sequence of uni-variate random processes. In this fashion, the traditional requirement of inverting high dimensional matrices is completely eliminated from the inference process, resulting in an efficient algorithm that is trivial to implement. Using cosmic large scale structure data as a showcase, we demonstrate the capabilities of our Gibbs sampling approach by performing a joint analysis of three dimensional density fields and corresponding power-spectra from Gaussian mock catalogues. These tests clearly demonstrate the ability of the algorithm to accurately provide measurements of the three dimensional density field and its power-spectrum and corresponding uncertainty quantification. Moreover, these tests reveal excellent numerical and statistical efficiency which will generally render the proposed algorithm a valuable addition to the toolbox of large scale Bayesian inference in cosmology and astrophysics.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Sparse sampling, galaxy bias, and voids

To study the impact of sparsity and galaxy bias on void statistics, we use a single large-volume, high-resolution N-body simulation to compare voids in multiple levels of subsampled dark matter, halo populations, and mock galaxies from a Halo Occupation Distribution model tuned to different galaxy survey densities. We focus our comparison on three key observational statistics: number functions, ellipticity distributions, and radial density profiles. We use the hierarchical tree structure of voids to interpret the impacts of sampling density and galaxy bias, and theoretical and empirical functions to describe the statistics in all our sample populations. We are able to make simple adjustments to theoretical expectations to offer prescriptions for translating from analytics to the void properties measured in realistic observations. We find that sampling density has a much larger effect on void sizes than galaxy bias. At lower tracer density, small voids disappear and the remaining voids are larger, more spherical, and have slightly steeper profiles. When a proper lower mass threshold is chosen, voids in halo distributions largely mimic those found in galaxy populations, except for ellipticities, where galaxy bias leads to higher values. We use the void density profile of Hamaus et al. (2014) to show that voids follow a self-similar and universal trend, allowing simple translations between voids studied in dark matter and voids identified in galaxy surveys. We have added the mock void catalogs used in this work to the Public Cosmic Void Catalog at http://www.cosmicvoids.net.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Testing cosmic geometry without dynamic distortions using voids

We propose a novel technique to probe the expansion history of the Universe based on the clustering statistics of cosmic voids. In particular, we compute their two-point statistics in redshift space on the basis of realistic mock galaxy catalogs and apply the Alcock-Paczynski test. In contrast to galaxies, we find void auto-correlations to be marginally affected by peculiar motions, providing a model-independent measure of cosmological parameters without systematics from redshift-space distortions. Because only galaxy-galaxy and void-galaxy correlations have been considered in these types of studies before, the presented method improves both statistical and systematic uncertainties on the product of angular diameter distance and Hubble rate, furnishing the potentially cleanest probe of cosmic geometry available to date.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The dark matter of galaxy voids

How do observed voids relate to the underlying dark matter distribution? To examine the spatial distribution of dark matter contained within voids identified in galaxy surveys, we apply Halo Occupation Distribution models representing sparsely and densely sampled galaxy surveys to a high-resolution N-body simulation. We compare these galaxy voids to voids found in the halo distribution, low-resolution dark matter, and high-resolution dark matter. We find that voids at all scales in densely sampled surveys - and medium- to large-scale voids in sparse surveys - trace the same underdensities as dark matter, but they are larger in radius by ~20%, they have somewhat shallower density profiles, and they have centers offset by ~0.4Rv rms. However, in void-to-void comparison we find that shape estimators are less robust to sampling, and the largest voids in sparsely sampled surveys suffer fragmentation at their edges. We find that voids in galaxy surveys always correspond to underdensities in the dark matter, though the centers may be offset. When this offset is taken into account, we recover almost identical radial density profiles between galaxies and dark matter. All mock catalogs used in this work are available at http://www.cosmicvoids.net.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

VIDE: The Void IDentification and Examination toolkit

We present VIDE, the Void IDentification and Examination toolkit, an open-source Python/C++ code for finding cosmic voids in galaxy redshift surveys and N-body simulations, characterizing their properties, and providing a platform for more detailed analysis. At its core, VIDE uses a substantially enhanced version of ZOBOV (Neyinck 2008) to calculate a Voronoi tessellation for estimating the density field and a performing a watershed transform to construct voids. Additionally, VIDE provides significant functionality for both pre- and post-processing: for example, vide can work with volume- or magnitude-limited galaxy samples with arbitrary survey geometries, or dark matter particles or halo catalogs in a variety of common formats. It can also randomly subsample inputs and includes a Halo Occupation Distribution model for constructing mock galaxy populations. VIDE uses the watershed levels to place voids in a hierarchical tree, outputs a summary of void properties in plain ASCII, and provides a Python API to perform many analysis tasks, such as loading and manipulating void catalogs and particle members, filtering, plotting, computing clustering statistics, stacking, comparing catalogs, and fitting density profiles. While centered around ZOBOV, the toolkit is designed to be as modular as possible and accommodate other void finders. VIDE has been in development for several years and has already been used to produce a wealth of results, which we summarize in this work to highlight the capabilities of the toolkit. VIDE is publicly available at http://bitbucket.org/cosmicvoids/vide public and http://www.cosmicvoids.net.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Voids in the SDSS DR9: observations, simulations, and the impact of the survey mask

We present and study cosmic voids identified using the watershed void finder VIDE in the Sloan Digital Sky Survey Data Release 9, compare these voids to ones identified in mock catalogs, and assess the impact of the survey mask on void statistics such as number functions, ellipticity distributions, and radial density profiles. The nearly 1,000 identified voids span three nearly volume-limited samples from redshift z = 0.43 to 0.7. For comparison we use 98 of the publicly available 2LPT-based mock galaxy catalogs of Manera et al., and also generate our own mock catalogs by applying a Halo Occupation Distribution model to an N-body simulation. We find that the mask reduces the number density of voids at all scales by a factor of three and slightly skews the relative size distributions. This engenders an increase in the mean ellipticity by roughly 30%. However, we find that radial density profiles are largely robust to the effects of the mask. We see excellent agreement between the data and both mock catalogs, and find no tension between the observed void properties and the properties derived from ΛCDM simulations. We have added the void catalogs from both data and mock galaxy populations discussed in this work to the Public Cosmic Void Catalog at http://www.cosmicvoids.net.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

A response to arXiv:1310.2791: A self-consistent public catalogue of voids and superclusters in the SDSS Data Release 7 galaxy surveys

Recently, Nadathur & Hotchkiss (2013) submitted a paper discussing a new cosmic void catalog. This paper includes claims about the void catalog described in Sutter et al. (2012). In this note, we respond to those claims, clarify some discrepancies between the text of Sutter et al. (2012) and the most recent version of the catalog, and provide some comments on the differences between our catalog and that of Nadathur & Hotchkiss (2013). All updates and documentation for our catalog are available at http://www.cosmicvoids.net.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

An Optimal and Model-Independent Measurement of the Intracluster Pressure Profile I: Methodology and First Applications

We present a statistically-optimal and model-independent method to extract the pressure profile of hot gas in the intracluster medium (ICM). Using the thermal Sunyaev-Zeldovich effect, we constrain the mean pressure profile of the ICM by appropriately considering all primary cosmic microwave background (CMB) and instrumental noise correlations, while using the maximum resolution and sensitivity of all frequency channels. As a first application, we analyze CMB maps of WMAP 9-year data through a study of the Meta-Catalogue of X-ray detected Clusters of galaxies (MCXC). We constrain the universal pressure profile out to 4R_500 with 15-sigma confidence, though our measurements are only significant out to R_200. Using a temperature profile constrained from X-ray observations, we measure the mean gas mass fraction out to R_200. Within statistical and systematic uncertainties, our constraints are compatible with the cosmic baryon fraction and the expected gas fraction in halos. While Planck multi-frequency CMB data are expected to reduce statistical uncertainties by a factor of 20, we argue that systematic errors in determining mass of clusters dominate the uncertainty in gas mass fraction measurements at the level of 20 percent.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

First measurement of the bulk flow of nearby galaxies using the cosmic microwave background

Peculiar velocities in the nearby Universe can be measured via the kinetic Sunyaev-Zel'dovich (kSZ) effect. Using a statistical method based on an optimised cross-correlation with nearby galaxies, we extract the kSZ signal generated by plasma halo of galaxies from the Cosmic Microwave Background (CMB) temperature anisotropies observed by the Wilkinson Microwave Anisotropy Probe (WMAP). Marginalising over the thermal Sunyaev-Zel'dovich contribution from clusters of galaxies, possible unresolved point source contamination, and Galactic foregrounds (dust, synchrotron and free-free emission), we report a kSZ bulk flow signal present at the 90% confidence level in the seven-year WMAP data. When only galaxies within 50 Mpc/h are included in the kSZ template we find a bulk flow in the CMB frame of |V|=533 +/- 263 km/s, in the direction l=324 +/- 27, b=-7 +/- 17, consistent with bulk flow measurements on a similar scale using classical distance indicators. We show how this comparison constrains the (ionised) baryonic budget in the local universe. On very large (~ 500 Mpc/h) scales, we find a 95% upper limit of 470 km/s, inconsistent with some analyses of bulk flow of clusters from the kSZ. We estimate that the significance of the bulk flow signal may increase to 3-5 sigma using data from the PLANCK probe.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A first application of the Alcock-Paczynski test to stacked cosmic voids

We report on the first application of the Alcock-Paczynski test to stacked voids in spectroscopic galaxy redshift surveys.We use voids from the Sutter et al. (2012) void catalog, which was derived from the Sloan Digital Sky Survey Data Release 7 main sample and luminous red galaxy catalogs. The construction of that void catalog removes potential shape measurement bias by using a modified version of the ZOBOV algorithm and by removing voids near survey boundaries and masks. We apply the shape-fitting procedure presented in Lavaux & Wandelt (2012) to ten void stacks out to redshift z=0.36. Combining these measurements, we determine the mean cosmologically induced "stretch" of voids in three redshift bins, with 1-sigma errors of 5-15%. The mean stretch is consistent with unity, providing no indication of a distortion induced by peculiar velocities. While the statistical errors are too large to detect the Alcock-Paczynski effect over our limited redshift range, this proof-of-concept analysis defines procedures that can be applied to larger spectroscopic galaxy surveys at higher redshifts to constrain dark energy using the expected statistical isotropy of structures that are minimally affected by uncertainties in galaxy velocity bias.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A public void catalog from the SDSS DR7 Galaxy Redshift Surveys based on the watershed transform

We produce the most comprehensive public void catalog to date using the Sloan Digital Sky Survey Data Release 7 main sample out to redshift z=0.2 and the luminous red galaxy sample out to z=0.44. Using a modified version of the parameter-free void finder ZOBOV, we fully take into account the presence of the survey boundary and masks. Our strategy for finding voids is thus appropriate for any survey configuration. We produce two distinct catalogs: a complete catalog including voids near any masks, which would be appropriate for void galaxy surveys, and a bias-free catalog of voids away from any masks, which is necessary for analyses that require a fair sampling of void shapes and alignments. Our discovered voids have effective radii from 5 to 135 h^-1 Mpc. We discuss basic catalog statistics such as number counts and redshift distributions and describe some additional data products derived from our catalog, such as radial density profiles and projected density maps. We find that radial profiles of stacked voids show a qualitatively similar behavior across nearly two decades of void radii and throughout the full redshift range.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Precision cosmography with stacked voids

We present a purely geometrical method for probing the expansion history of the Universe from the observation of the shape of stacked voids in spectroscopic redshift surveys. Our method is an Alcock-Paczynski (AP) test based on the average sphericity of voids posited on the local isotropy of the Universe. It works by comparing the temporal extent of cosmic voids along the line of sight with their angular, spatial extent. We describe the algorithm that we use to detect and stack voids in redshift shells on the light cone and test it on mock light cones produced from N-body simulations. We establish a robust statistical model for estimating the average stretching of voids in redshift space and quantify the contamination by peculiar velocities. Finally, assuming that the void statistics that we derive from N-body simulations is preserved when considering galaxy surveys, we assess the capability of this approach to constrain dark energy parameters. We report this assessment in terms of the figure of merit (FoM) of the dark energy task force and in particular of the proposed EUCLID mission which is particularly suited for this technique since it is a spectroscopic survey. The FoM due to stacked voids from the EUCLID wide survey may double that of all other dark energy probes derived from EUCLID data alone (combined with Planck priors). In particular, voids seem to outperform Baryon Acoustic Oscillations by an order of magnitude. This result is consistent with simple estimates based on mode-counting. The AP test based on stacked voids may be a significant addition to the portfolio of major dark energy probes and its potentialities must be studied in detail.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Robust, data-driven inference in non-linear cosmostatistics

We discuss two projects in non-linear cosmostatistics applicable to very large surveys of galaxies. The first is a Bayesian reconstruction of galaxy redshifts and their number density distribution from approximate, photometric redshift data. The second focuses on cosmic voids and uses them to construct cosmic spheres that allow reconstructing the expansion history of the Universe using the Alcock-Paczynski test. In both cases we find that non-linearities enable the methods or enhance the results: non-linear gravitational evolution creates voids and our photo-z reconstruction works best in the highest density (and hence most non-linear) portions of our simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Straightening the Density-Displacement Relation with a Logarithmic Transform

We investigate the use of a logarithmic density variable in estimating the Lagrangian displacement field, motivated by the success of a logarithmic transformation in restoring information to the matter power spectrum. The logarithmic relation is an extension of the linear relation, motivated by the continuity equation, in which the density field is assumed to be proportional to the divergence of the displacement field; we compare the linear and logarithmic relations by measuring both of these fields directly in a cosmological N-body simulation. The relative success of the logarithmic and linear relations depends on the scale at which the density field is smoothed. Thus we explore several ways of measuring the density field, including Cloud-In-Cell smoothing, adaptive smoothing, and the (scale-independent) Delaunay tessellation, and we use both a Fourier space and a geometrical tessellation approach to measuring the divergence. We find that the relation between the divergence of the displacement field and the density is significantly tighter with a logarithmic density variable, especially at low redshifts and for very small (~2 Mpc/h) smoothing scales. We find that the grid-based methods are more reliable than the tessellation-based method of calculating both the density and the divergence fields, though in both cases the logarithmic relation works better in the appropriate regime, which corresponds to nonlinear scales for the grid-based methods and low densities for the tessellation-based method.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The 2M++ galaxy redshift catalogue

Peculiar velocities arise from gravitational instability, and thus are linked to the surrounding distribution of matter. In order to understand the motion of the Local Group with respect to the Cosmic Microwave Background, a deep all-sky map of the galaxy distribution is required. Here we present a new redshift compilation of 69~160 galaxies, dubbed 2M++, to map large-scale structures of the Local Universe over nearly the whole sky, and reaching depths of K <= 12.5, or 200 Mpc/h. The target catalogue is based on the Two-Micron-All-Sky Extended Source Catalog (2MASS-XSC). The primary sources of redshifts are the 2MASS Redshift Survey, the 6dF galaxy redshift survey and the Sloan Digital Sky Survey (DR7). We assess redshift completeness in each region and compute the weights required to correct for redshift incompleteness and apparent magnitude limits, and discuss corrections for incompleteness in the Zone of Avoidance. We present the density field for this survey, and discuss the importance of large-scale structures such as the Shapley Concentration.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Fast CMB lensing using statistical interpolation on the sphere

We describe a accurate and fast pixel-based statistical method to interpolate fields of arbitrary spin on the sphere. We call this method Fast and Lean Interpolation on the Sphere (FLINTS). The method predicts the optimal interpolated values based on the theory of isotropic Gaussian random fields and provides an accurate error estimate at no additional cost. We use this method to compute lensed Cosmic Microwave Background (CMB) maps precisely and quickly, achieving a relative precision of 0.02% at a HEALPix resolution of Nside=4096, for a bandlimit of l_max=4096 in the same time it takes to simulate the original, unlensed CMB map. The method is suitable for efficient, distributed memory parallelization. The power spectra of our lensed maps are accurate to better than 0.5% at l=3000 for the temperature, the E and B mode of the polarization. As expected theoretically, we demonstrate that, on realistic cases, this method is between two to three orders of magnitude more precise than other known interpolation methods for the same computational cost.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Precision constrained simulation of the Local Universe

We use the formalism of constrained Gaussian random field to compute a precise large scale simulation of the 60 Mpc/h volume of our Local Universe. We derive the constraints from the reconstructed peculiar velocities of the 2MASS Redshift Survey. We obtain a correlation of 0.97 between the log-density field of the dark matter distribution of the simulation and the log-density of observed galaxies of the Local Universe. We achieve a good comparison of the simulated velocity field to the observed velocity field obtained from the galaxy distances of the NBG-3k. At the end, we compare the two-point correlation function of both the 2MRS galaxies and of the dark matter particles of the simulation. We conclude that this method is a very promising technique of exploring the dynamics and the particularities the Universe in our neighbourhood.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Precision cosmology with voids: definition, methods, dynamics

We propose a new definition of cosmic voids based on methods of Lagrangian orbit reconstruction as well as an algorithm to find them in actual data called DIVA. Our technique is intended to yield results which can be modeled sufficiently accurately to create a new probe of precision cosmology. We then develop an analytical model of the ellipticity of voids found by our method based on Zel'dovich approximation. We measure in N-body simulation that this model is precise at the 0.1% level for the mean ellipticity of voids of size greater than ~4 Mpc/h. We estimate that at this scale, we are able to predict the ellipticity with an accuracy of 0.02. Finally, we compare the distribution of void shapes in N-body simulation for two different equations of state w of the dark energy. We conclude that our method is far more accurate than Eulerian methods and is therefore promising as a precision probe of dark energy phenomenology.

0 citations0 reviewsNo code linkedOpen access