We present a general method for accelerating by more than an order of magnitude the convolution of pixelated functions on the sphere with a radially-symmetric kernel. Our method splits the kernel into a compact real-space component and a compact sphe
rical harmonic space component. These components can then be convolved in parallel using an inexpensive commodity GPU and a CPU. We provide models for the computational cost of both real-space and Fourier space convolutions and an estimate for the approximation error. Using these models we can determine the optimum split that minimizes the wall clock time for the convolution while satisfying the desired error bounds. We apply this technique to the problem of simulating a cosmic microwave background (CMB) anisotropy sky map at the resolution typical of the high resolution maps produced by the Planck mission. For the main Planck CMB science channels we achieve a speedup of over a factor of ten, assuming an acceptable fractional rms error of order 1.e-5 in the power spectrum of the output map.
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 a
nalysis. 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.
Taking N-body simulations with volumes and particle densities tuned to match the SDSS DR7 spectroscopic main sample, we assess the ability of current void catalogs (e.g., Sutter et al. 2012b) to distinguish a model of coupled dark matter-dark energy
from {Lambda}CDM cosmology using properties of cosmic voids. Identifying voids with the VIDE toolkit, we find no statistically significant differences in the ellipticities, but find that coupling produces a population of significantly larger voids, possibly explaining the recent result of Tavasoli et al. (2013). In addition, we use the universal density profile of Hamaus et al. (2014) to quantify the relationship between coupling and density profile shape, finding that the coupling produces broader, shallower, undercompensated profiles for large voids by thinning the walls between adjacent medium-scale voids. We find that these differences are potentially measurable with existing void catalogs once effects from survey geometries and peculiar velocities are taken into account.
We perform an Alcock-Paczynski test using stacked cosmic voids identified in the SDSS Data Release 7 main sample and Data Release 10 LOWZ and CMASS samples. We find ~1,500 voids out to redshift $0.6$ using a heavily modified and extended version of t
he watershed algorithm ZOBOV, which we call VIDE (Void IDentification and Examination). To assess the impact of peculiar velocities we use the mock void catalogs presented in Sutter et al. (2013). We find a constant uniform flattening of 14% along the line of sight when peculiar velocities are included. This flattening appears universal for all void sizes at all redshifts and for all tracer densities. We also use these mocks to identify an optimal stacking strategy. After correcting for systematic effects we find that our Alcock-Paczynski measurement leads to a preference of our best-fit value of $Omega_{rm M}sim 0.15$ over $Omega_{rm M} = 1.0$ by a likelihood ratio of 10. Likewise, we find a factor of $4.5$ preference of the likelihood ratio for a $Lambda$CDM $Omega_{rm M} = 0.3$ model and a null measurement. Taken together, we find substantial evidence for the Alcock-Paczynski signal in our sample of cosmic voids. Our assessment using realistic mocks suggests that measurements with future SDSS releases and other surveys will provide tighter cosmological parameter constraints. The void-finding algorithm and catalogs used in this work will be made publicly available at http://www.cosmicvoids.net.
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 suc
h 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 {Lambda}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.
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 Occupat
ion 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.
