Cosmic voids are biased tracers of the large-scale structure of the universe. Separate universe simulations (SUS) enable accurate measurements of this biasing relation by implementing the peak-background split (PBS). In this work, we apply the SUS technique to measure the void bias parameters. We confirm that the PBS argument works well for underdense tracers. The response of the void size distribution depends on the void radius. For voids larger (smaller) than the size at the peak of the distribution, the void abundance responds negatively (positively) to a long wavelength mode. The linear bias from the SUS is in good agreement with the cross power spectrum measurement on large scales. Using the SUS, we have detected the quadratic void bias for the first time in simulations. We find that $ b_2 $ is negative when the magnitude of $ b_1 $ is small, and that it becomes positive and increases rapidly when $ |b_1| $ increases. We compare the results from voids identified in the halo density field with those from the dark matter distribution, and find that the results are qualitatively similar, but the biases generally shift to the larger voids sizes.

Cosmic voids offer an extraordinary opportunity to study the effects of massive neutrinos on cosmological scales. Because they are freely streaming, neutrinos can penetrate the interior of voids more easily than cold dark matter or baryons, which makes their relative contribution to the mass budget in voids much higher than elsewhere in the Universe. In simulations it has recently been shown how various characteristics of voids in the matter distribution are affected by neutrinos, such as their abundance, density profiles, dynamics, and clustering properties. However, the tracers used to identify voids in observations (e.g., galaxies or halos) are affected by neutrinos as well, and isolating the unique neutrino signatures inherent to voids becomes more difficult. In this paper we make use of the DEMNUni suite of simulations to investigate the clustering bias of voids in Fourier space as a function of their core density and compensation. We find a clear dependence on the sum of neutrino masses that remains significant even for void statistics extracted from halos. In particular, we observe that the amplitude of the linear void bias increases with neutrino mass for voids defined in dark matter, whereas this trend gets reversed and slightly attenuated when measuring the relative void-halo bias using voids identified in the halo distribution. Finally, we argue how the original behaviour can be restored when considering observations of the total matter distribution (e.g. via weak lensing), and comment on scale-dependent effects in the void bias that may provide additional information on neutrinos in the future.

We test a theoretical description of the void size distribution function against direct estimates from halo catalogues of the DEMNUni suite of large cosmological simulations. Besides standard $\Lambda$CDM, we consider deviations of the dark energy equation of state from $w=-1$, corresponding to four combinations in the popular Chevallier-Polarski-Linder parametrisation: $w_0=-0.9;-1.1$, $w_a=-0.3;0.3$. The theoretical void size function model, relying on the Sheth & van de Weygaert double barrier excursion set formalism, provides an accurate description of the simulation measurements for the different dark energy models considered, within the statistical errors. The model remains accurate for any value of the threshold for void formation $\delta_\mathrm{v}$. Its robust consistency with simulations demonstrates that the theoretical void size function can be applied to real data as a sensitive tool to constrain dark energy.

What are the mass and galaxy profiles of cosmic voids? In this paper we use two methods to extract voids in the Dark Energy Survey (DES) Year 1 redMaGiC galaxy sample to address this question. We use either 2D slices in projection, or the 3D distribution of galaxies based on photometric redshifts to identify voids. For the mass profile, we measure the tangential shear profiles of background galaxies to infer the excess surface mass density. The signal-to-noise ratio for our lensing measurement ranges between 10.7 and 14.0 for the two void samples. We infer their 3D density profiles by fitting models based on N-body simulations and find good agreement for void radii in the range 15-85 Mpc. Comparison with their galaxy profiles then allows us to test the relation between mass and light at the 10%-level, the most stringent test to date. We find very similar shapes for the two profiles, consistent with a linear relationship between mass and light both within and outside the void radius. We validate our analysis with the help of simulated mock catalogues and estimate the impact of photometric redshift uncertainties on the measurement. Our methodology can be used for cosmological applications, including tests of gravity with voids. This is especially promising when the lensing profiles are combined with spectroscopic measurements of void dynamics via redshift-space distortions.

Cosmic voids gravitationally lens the cosmic microwave background (CMB) radiation, resulting in a distinct imprint on degree scales. We use the simulated CMB lensing convergence map from the MICE N-body simulation to calibrate our detection strategy for a given void definition and galaxy tracer density. We then identify cosmic voids in DES Year 1 data and stack the Planck 2015 lensing convergence map on their locations, probing the consistency of simulated and observed void lensing signals. When fixing the shape of the stacked convergence profile to that calibrated from simulations, we find imprints at the $3{\sigma}$ significance level for various analysis choices. The best measurement strategies based on the MICE calibration process yield $S/N \sim 4$ for DES Y1, and the best-fit amplitude recovered from the data is consistent with expectations from MICE ($A \sim 1$). Given these results as well as the agreement between them and N-body simulations, we conclude that the previously reported excess integrated Sachs-Wolfe (ISW) signal associated with cosmic voids in DES Y1 has no counterpart in the Planck CMB lensing map.

We study the large-scale bias parameter of cosmic voids with primordial non-Gaussian (PNG) initial conditions of the local type. In this scenario, the dark matter halo bias exhibits a characteristic scale dependence on large scales, which has been recognized as one of the most promising probes of the local PNG. Using a suite of $N$-body simulations with Gaussian and non-Gaussian initial conditions, we find that the void bias features scale-dependent corrections on large scales, similar to its halo counterpart. We find excellent agreement between the numerical measurement of the PNG void bias and the general peak-background split prediction. Contrary to halos, large voids anti-correlate with the dark matter density field, and the large-scale Gaussian void bias ranges from positive to negative values depending on void size and redshift. Thus, the information in the clustering of voids can be complementary to that of the halos. Using the Fisher matrix formalism for multiple tracers, we demonstrate that including the scale-dependent bias information from voids, constraints on the PNG parameter $f_{\rm NL}$ can be tightened by a factor of two compared to the accessible information from halos alone, when the sampling density of tracers reaches $4 \times 10^{-3} \, h^3 \mathrm{Mpc}^{-3} $.

The aim of this study is to distinguish genuine cosmic voids, found in a galaxy catalog by the void finder ZOBOV-VIDE, from under-dense regions in a Poisson distribution of objects. For this purpose, we perform two multivariate analyses using the following physical void characteristics: volume, redshift, density contrast, minimum density, contrast significance and number of member galaxies of the void. The multivariate analyses are trained on a catalog of voids obtained from a random Poisson distribution of points, used as background, and a catalog of voids identified in a mock galaxy catalog, used as signal. The classifications are then applied to voids extracted from the Data Release 12 sample of Luminous Red Galaxies in the redshift range 0.45 $\leq$ z $\leq$ 0.7 from the Sloan Digital Sky Survey Baryon Oscillation Spectroscopic Survey (SDSS BOSS DR12 CMASS). Our results show that the resulting void catalog is nearly free of contamination by Poisson noise. We also study the effect of tracer sparsity and bias on the classification efficiencies.

Cosmic voids, the less dense patches of the Universe, are promising laboratories to extract cosmological information. Thanks to their unique low density character, voids are extremely sensitive to diffuse components such as neutrinos and dark energy, and represent ideal environments to study modifications of gravity, where the effects of such modifications are expected to be more prominent. Robust void-related observables, including for example redshift-space distortions (RSD) and weak lensing around voids, are a promising way to chase and test new physics. Cosmological analysis of the large-scale structure of the Universe predominantly relies on the high density regions. Current and upcoming surveys are designed to optimize the extraction of cosmological information from these zones, but leave voids under-exploited. A dense, large area spectroscopic survey with imaging capabilities is ideal to exploit the power of voids fully. Besides helping illuminate the nature of dark energy, modified gravity, and neutrinos, this survey will give access to a detailed map of under-dense regions, providing an unprecedented opportunity to observe and study a so far under-explored galaxy population.

The largest structures in the cosmic web probe the dynamical nature of dark energy through their integrated Sachs-Wolfe imprints. In the strength of the signal, typical cosmic voids have shown good consistency with expectation $A_{\rm ISW}=\Delta T^{\rm data} / \Delta T^{\rm theory}=1$, given the substantial cosmic variance. Discordantly, large-scale hills in the gravitational potential, or supervoids, have shown excess signals. In this study, we mapped out 87 new supervoids in the total 5000 deg$^2$ footprint of the Dark Energy Survey at $0.2<z<0.9$ to probe these anomalous claims. We found an excess imprinted profile with $ A_{\rm ISW}\approx4.1\pm2.0$ amplitude. The combination with independent BOSS data reveals an ISW imprint of supervoids at the $3.3\sigma$ significance level with an enhanced $A_{\rm ISW}\approx5.2\pm1.6$ amplitude. The tension with $\Lambda$CDM predictions is equivalent to $2.6\sigma$ and remains unexplained.

We constrain the mass--richness scaling relation of redMaPPer galaxy clusters identified in the Dark Energy Survey Year 1 data using weak gravitational lensing. We split clusters into $4\times3$ bins of richness $\lambda$ and redshift $z$ for $\lambda\geq20$ and $0.2 \leq z \leq 0.65$ and measure the mean masses of these bins using their stacked weak lensing signal. By modeling the scaling relation as $\langle M_{\rm 200m}|\lambda,z\rangle = M_0 (\lambda/40)^F ((1+z)/1.35)^G$, we constrain the normalization of the scaling relation at the 5.0 per cent level as $M_0 = [3.081 \pm 0.075 ({\rm stat}) \pm 0.133 ({\rm sys})] \cdot 10^{14}\ {\rm M}_\odot$ at $\lambda=40$ and $z=0.35$. The richness scaling index is constrained to be $F=1.356 \pm 0.051\ ({\rm stat})\pm 0.008\ ({\rm sys})$ and the redshift scaling index $G=-0.30\pm 0.30\ ({\rm stat})\pm 0.06\ ({\rm sys})$. These are the tightest measurements of the normalization and richness scaling index made to date. We use a semi-analytic covariance matrix to characterize the statistical errors in the recovered weak lensing profiles. Our analysis accounts for the following sources of systematic error: shear and photometric redshift errors, cluster miscentering, cluster member dilution of the source sample, systematic uncertainties in the modeling of the halo--mass correlation function, halo triaxiality, and projection effects. We discuss prospects for reducing this systematic error budget, which dominates the uncertainty on $M_0$. Our result is in excellent agreement with, but has significantly smaller uncertainties than, previous measurements in the literature, and augurs well for the power of the DES cluster survey as a tool for precision cosmology and upcoming galaxy surveys such as LSST, Euclid and WFIRST.

Luminous tracers of large-scale structure are not entirely representative of the distribution of mass in our Universe. As they arise from the highest peaks in the matter density field, the spatial distribution of luminous objects is biased towards those peaks. On large scales, where density fluctuations are mild, this bias simply amounts to a constant offset in the clustering amplitude of the tracer, known as linear bias. In this work we focus on the relative bias between galaxies and galaxy clusters that are located inside and in the vicinity of cosmic voids, extended regions of relatively low density in the large-scale structure of the Universe. With the help of hydro-dynamical simulations we verify that the relation between galaxy and cluster overdensity around voids remains linear. Hence, the void-centric density profiles of different tracers can be linked by a single multiplicative constant. This amounts to the same value as the relative linear bias between tracers for the largest voids in the sample. For voids of small sizes, which typically arise in higher density regions, this constant has a higher value, possibly showing an environmental dependence similar to that observed for the linear bias itself. We confirm our findings by analysing mocks and data obtained during the first year of observations by the Dark Energy Survey. As a side product, we present the first catalogue of three-dimensional voids extracted from a photometric survey with a controlled photo-z uncertainty. Our results will be relevant in forthcoming analyses that attempt to use voids as cosmological probes.

The large-scale structure of the universe can only be observed via luminous tracers of the dark matter. However, the clustering statistics of tracers are biased and depend on various properties, such as their host-halo mass and assembly history. On very large scales this tracer bias results in a constant offset in the clustering amplitude, known as linear bias. Towards smaller nonlinear scales, this is no longer the case and tracer bias becomes a complicated function of scale and time. We focus on tracer bias centred on cosmic voids, depressions of the density field that spatially dominate the universe. We consider three types of tracers: galaxies, galaxy clusters and AGN, extracted from the hydrodynamical simulation Magneticum Pathfinder. In contrast to common clustering statistics that focus on auto-correlations of tracers, we find that void-tracer cross-correlations are successfully described by a linear-bias relation. The tracer-density profile of voids can thus be related to their matter-density profile by a single number. We show that it coincides with the linear tracer bias extracted from the large-scale auto-correlation function and expectations from theory, if sufficiently large voids are considered. For smaller voids we observe a shift towards higher values. This has important consequences on cosmological parameter inference, as the problem of unknown tracer bias is alleviated up to a constant number. The smallest scales in existing datasets become accessible to simpler models, providing numerous modes of the density field that have been disregarded so far, but may help to further reduce statistical errors in constraining cosmology.

We perform a comprehensive redshift-space distortion analysis based on cosmic voids in the large-scale distribution of galaxies observed with the Sloan Digital Sky Survey. To this end, we measure multipoles of the void-galaxy cross-correlation function and compare them with standard model predictions in cosmology. Merely considering linear-order theory allows us to accurately describe the data on the entire available range of scales and to probe void-centric distances down to about $2h^{-1}{\rm Mpc}$. Common systematics, such as the Fingers-of-God effect, scale-dependent galaxy bias, and nonlinear clustering do not seem to play a significant role in our analysis. We constrain the growth rate of structure via the redshift-space distortion parameter $\beta$ at two median redshifts, $\beta(\bar{z}=0.32)=0.599^{+0.134}_{-0.124}$ and $\beta(\bar{z}=0.54)=0.457^{+0.056}_{-0.054}$, with a precision that is competitive with state-of-the-art galaxy-clustering results. While the high-redshift constraint perfectly agrees with model expectations, we observe a mild $2\sigma$ deviation at $\bar{z}=0.32$, which increases to $3\sigma$ when the data is restricted to the lowest available redshift range of $0.15<z<0.33$.

We present measurements of angular cross power spectra between galaxies and optically-selected galaxy clusters in the final photometric sample of the Sloan Digital Sky Survey (SDSS). We measure the auto- and cross-correlations between galaxy and cluster samples, from which we extract the effective biases and study the shot noise properties. We model the non-Poissonian shot noise by introducing an effective number density of tracers and fit for this quantity. We find that we can only describe the cross-correlation of galaxies and galaxy clusters, as well as the auto-correlation of galaxy clusters, on the relevant scales using a non-Poissonian shot noise contribution. The values of effective bias we finally measure for a volume-limited sample are $b_{cc}=4.09 \pm 0.47$ for the cluster auto-correlation and $b_{gc}=2.15 \pm 0.09$ for the galaxy-cluster cross-correlation. We find that these results are consistent with expectations from the auto-correlations of galaxies and clusters and are in good agreement with previous studies. The main result is two-fold: firstly we provide a measurement of the cross-correlation of galaxies and clusters, which can be used for further cosmological analysis, and secondly we describe an effective treatment of the shot noise.

Galaxies and their dark matter halos populate a complicated filamentary network around large, nearly empty regions known as cosmic voids. Cosmic voids are usually identified in spectroscopic galaxy surveys, where 3D information about the large-scale structure of the Universe is available. Although an increasing amount of photometric data is being produced, its potential for void studies is limited since photometric redshifts induce line-of-sight position errors of $\sim50$ Mpc/$h$ or more that can render many voids undetectable. In this paper we present a new void finder designed for photometric surveys, validate it using simulations, and apply it to the high-quality photo-$z$ redMaGiC galaxy sample of the Dark Energy Survey Science Verification (DES-SV) data. The algorithm works by projecting galaxies into 2D slices and finding voids in the smoothed 2D galaxy density field of the slice. Fixing the line-of-sight size of the slices to be at least twice the photo-$z$ scatter, the number of voids found in these projected slices of simulated spectroscopic and photometric galaxy catalogs is within 20% for all transverse void sizes, and indistinguishable for the largest voids of radius $\sim 70$ Mpc/$h$ and larger. The positions, radii, and projected galaxy profiles of photometric voids also accurately match the spectroscopic void sample. Applying the algorithm to the DES-SV data in the redshift range $0.2<z<0.8$, we identify 87 voids with comoving radii spanning the range 18-120 Mpc/$h$, and carry out a stacked weak lensing measurement. With a significance of $4.4\sigma$, the lensing measurement confirms the voids are truly underdense in the matter field and hence not a product of Poisson noise, tracer density effects or systematics in the data. It also demonstrates, for the first time in real data, the viability of void lensing studies in photometric surveys.

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 $\Omega_\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.

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 $\Lambda$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 $\sigma_{f/b}\left/(f/b)\right.\sim2\%$ ($20\%$) and $\sigma_{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.

We show that the number of observed voids in galaxy redshift surveys is a sensitive function of the equation of state of dark energy. Using the Fisher matrix formalism we find the error ellipses in the $w_0-w_a$ plane when the equation of state of dark energy is assumed to be of the form $w_{CPL}(z)=w_0 +w_a z/(1+z)$. We forecast the number of voids to be observed with the ESA Euclid satellite and the NASA WFIRST mission, taking into account updated details of the surveys to reach accurate estimates of their power. The theoretical model for the forecast of the number of voids is based on matches between abundances in simulations and the analytical prediction. To take into account the uncertainties within the model, we marginalize over its free parameters when calculating the Fisher matrices. The addition of the void abundance constraints to the data from Planck, HST and supernova survey data noticeably tighten the $w_0-w_a$ parameter space. We thus quantify the improvement in the constraints due to the use of voids and demonstrate that the void abundance is a sensitive new probe for the dark energy equation of state.

What do we know about voids in the dark matter distribution given the Sloan Digital Sky Survey (SDSS) and assuming the $\Lambda\mathrm{CDM}$ model? Recent application of the Bayesian inference algorithm BORG to the SDSS Data Release 7 main galaxy sample has generated detailed Eulerian and Lagrangian representations of the large-scale structure as well as the possibility to accurately quantify corresponding uncertainties. Building upon these results, we present constrained catalogs of voids in the Sloan volume, aiming at a physical representation of dark matter underdensities and at the alleviation of the problems due to sparsity and biasing on galaxy void catalogs. To do so, we generate data-constrained reconstructions of the presently observed large-scale structure using a fully non-linear gravitational model. We then find and analyze void candidates using the VIDE toolkit. Our methodology therefore predicts the properties of voids based on fusing prior information from simulations and data constraints. For usual void statistics (number function, ellipticity distribution and radial density profile), all the results obtained are in agreement with dark matter simulations. Our dark matter void candidates probe a deeper void hierarchy than voids directly based on the observed galaxies alone. The use of our catalogs therefore opens the way to high-precision void cosmology at the level of the dark matter field. We will make the void catalogs used in this work available at http://www.cosmicvoids.net.

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.

We study the clustering of voids using $N$-body simulations and simple theoretical models. The excursion-set formalism describes fairly well the abundance of voids identified with the watershed algorithm, although the void formation threshold required is quite different from the spherical collapse value. The void cross bias $b_{\rm c} $ is measured and its large-scale value is found to be consistent with the peak background split results. A simple fitting formula for $b_{\rm c} $ is found. We model the void auto-power spectrum taking into account the void biasing and exclusion effect. A good fit to the simulation data is obtained for voids with radii $\gtrsim$ 30 Mpc/$h$, especially when the void biasing model is extended to 1-loop order. However, the best-fit bias parameters do not agree well with the peak-background split results. Being able to fit the void auto-power spectrum is particularly important not only because it is the direct observable in galaxy surveys, but also our method enables us to treat the bias parameters as nuisance parameters, which are sensitive to the techniques used to identify voids.

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.

Understanding the internal structure and spatial distribution of cosmic voids is crucial when considering them as probes of cosmology. We present recent advances in modeling void density- and velocity-profiles in real space, as well as void two-point statistics in redshift space, by examining voids identified via the watershed transform in state-of-the-art $\Lambda$CDM n-body simulations and mock galaxy catalogs. The simple and universal characteristics that emerge from these statistics indicate the self-similarity of large-scale structure and suggest cosmic voids to be among the most pristine objects to consider for future studies on the nature of dark energy, dark matter and modified gravity.

We discuss various applications of 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. Based on a substantially enhanced version of ZOBOV, VIDE not only finds voids, but also summarizes their properties, extracts statistical information, and provides a Python-based platform for more detailed analysis, such as manipulating void catalogs and particle members, filtering, plotting, computing clustering statistics, stacking, comparing catalogs, and fitting density profiles. VIDE also provides significant additional functionality for pre-processing inputs: 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 has been used for a wide variety of applications, from discovering a universal density profile to estimating primordial magnetic fields, and is publicly available at http://bitbucket.org/cosmicvoids/vide_public and http://www.cosmicvoids.net.

We investigate the formation, growth, merger history, movement, and destruction of cosmic voids detected via the watershed transform code VIDE in a cosmological N-body dark matter {\Lambda}CDM simulation. By adapting a method used to construct halo merger trees, we are able to trace individual voids back to their initial appearance and record the merging and evolution of their progenitors at high redshift. For the scales of void sizes captured in our simulation, we find that the void formation rate peaks at scale factor 0.3, which coincides with a growth in the void hierarchy and the emergence of dark energy. Voids of all sizes appear at all scale factors, though the median initial void size decreases with time. When voids become detectable they have nearly their present-day volumes. Almost all voids have relatively stable growth rates and suffer only infrequent minor mergers. Dissolution of a void via merging is very rare. Instead, most voids maintain their distinct identity as annexed subvoids of a larger parent. The smallest voids are collapsing at the present epoch, but void destruction ceases after scale factor 0.3. In addition, voids centers tend to move very little, less than 0.01 of their effective radii per ln a, over their lifetimes. Overall, most voids exhibit little radical dynamical evolution; their quiet lives make them pristine probes of cosmological initial conditions and the imprint of dark energy.

We present a simple empirical function for the average density profile of cosmic voids, identified via the watershed technique in $\Lambda$CDM N-body simulations. This function is universal across void size and redshift, accurately describing a large radial range of scales around void centers with only two free parameters. In analogy to halo density profiles, these parameters describe the scale radius and the central density of voids. While we initially start with a more general four-parameter model, we find two of its parameters to be redundant, as they follow linear trends with the scale radius in two distinct regimes of the void sample, separated by its compensation scale. Assuming linear theory, we derive an analytic formula for the velocity profile of voids and find an excellent agreement with the numerical data as well. In our companion paper [Sutter et al., Mon. Not. R. Astron. Soc. 442, 462 (2014)] the presented density profile is shown to be universal even across tracer type, properly describing voids defined in halo and galaxy distributions of varying sparsity, allowing us to relate various void populations by simple rescalings. This provides a powerful framework to match theory and simulations with observational data, opening up promising perspectives to constrain competing models of cosmology and gravity.

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.

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.

The clustering of galaxies in ongoing and upcoming galaxy surveys contains a wealth of cosmological information, but extracting this information is a non-trivial task since galaxies and their host haloes are stochastic tracers of the matter density field. This stochasticity is usually modeled as Poisson shot noise, which is constant as a function of wavenumber with amplitude given by 1/n, where n is the number density of galaxies. Here we use dark matter haloes in N-body simulations to show evidence for deviations from this simple behaviour and develop models that explain the behaviour of the large scale stochasticity. First, haloes are extended, non-overlapping objects, i.e., their correlation function needs to go to -1 on small scales. This leads to a negative correction to the stochasticity relative to the Poisson value at low wavenumber k, decreasing to zero for wavenumbers large compared to the inverse exclusion scale. Second, haloes show a non-linear enhancement of clustering outside the exclusion scale, leading to a positive stochasticity correction. Both of these effects go to zero for high-k, making the stochasticity scale dependent even for k<0.1 h/Mpc. We show that the corrections in the low-k regime are the same in Eulerian and Lagrangian space, but that the transition scale is pushed to smaller scales for haloes observed at present time, relative to the initial conditions. These corrections vary with halo mass and redshift. We also discuss simple applications of these effects to the galaxy samples with non-vanishing satellite fraction, where the stochasticity can again deviate strongly from the fiducial Poisson expectation. Overall these effects affect the clustering of galaxies at a level of a few percent even on very large scales and need to be modelled properly if we want to extract high precision cosmological information from the upcoming galaxy surveys.

Using multiple tracers of large-scale structure allows to evade the limitations imposed by sampling variance for some parameters of interest in cosmology. We demonstrate the optimal way of carrying out a multitracer analysis in a galaxy redshift survey by considering the principal components of the shot noise matrix from two-point clustering statistics. We show how to construct two tracers that maximize the benefits of sampling variance and shot noise cancellation using optimal weights. On the basis of high-resolution N-body simulations of dark matter halos we apply this technique to the analysis of redshift-space distortions and demonstrate how constraints on the growth rate of structure formation can be substantially improved. The primary limitation are nonlinear effects, which cause significant biases in the method already at scales of k<0.1h/Mpc, suggesting the need to develop nonlinear models of redshift-space distortions in order to extract the maximum information from future redshift surveys. Nonetheless we find gains of a factor of a few in constraints on the growth rate achievable when merely the linear regime of a galaxy survey like EUCLID is considered.

The proper general relativistic description of the observed galaxy power spectrum is substantially different from the standard Newtonian description on large scales, providing a unique opportunity to test general relativity on horizon scales. Using the Einstein equations, the general relativistic effects can be classified as two new terms that represent the velocity and the gravitational potential, coupling to the time evolution of galaxy number density and Hubble parameter. Compared to the dominant density and velocity redshift-space distortion terms, the former scales as H/k and correlates the real and imaginary parts of the Fourier modes, while the latter scales as (H/k)^2, where k is the comoving wave number and H is the conformal Hubble parameter. We use the recently developed methods to reduce the sampling variance and shot noise to show that in an all sky galaxy redshift survey at low redshift the velocity term can be measured at 10-sigma confidence level, if one can utilize halos of mass M>10^{10} Msun, while the gravitational potential term itself can only be marginally detected. We also demonstrate that the general relativistic effect is not degenerate with the primordial non-Gaussian signature in galaxy bias, and the ability to detect the primordial non-Gaussianity is little compromised.

Kaiser redshift-space distortion formula describes well the clustering of galaxies in redshift surveys on small scales, but there are numerous additional terms that arise on large scales. Some of these terms can be described using Newtonian dynamics and have been discussed in the literature, while the others require proper general relativistic description that was only recently developed. Accounting for these terms in galaxy clustering is the first step toward tests of general relativity on horizon scales. The effects can be classified as two terms that represent the velocity and the gravitational potential contributions. Their amplitude is determined by effects such as the volume and luminosity distance fluctuation effects and the time evolution of galaxy number density and Hubble parameter. We compare the Newtonian approximation often used in the redshift-space distortion literature to the fully general relativistic equation, and show that Newtonian approximation accounts for most of the terms contributing to velocity effect. We perform a Fisher matrix analysis of detectability of these terms and show that in a single tracer survey they are completely undetectable. To detect these terms one must resort to the recently developed methods to reduce sampling variance and shot noise. We show that in an all-sky galaxy redshift survey at low redshift the velocity term can be measured at a few sigma if one can utilize halos of mass M>10^12 Msun (this can increase to 10-sigma or more in some more optimistic scenarios), while the gravitational potential term itself can only be marginally detected. We also demonstrate that the general relativistic effect is not degenerate with the primordial non-Gaussian signature in galaxy bias, and the ability to detect primordial non-Gaussianity is little compromised.

One of the main signatures of primordial non-Gaussianity of the local type is a scale-dependent correction to the bias of large-scale structure tracers such as galaxies or clusters, whose amplitude depends on the bias of the tracers itself. The dominant source of noise in the power spectrum of the tracers is caused by sampling variance on large scales (where the non-Gaussian signal is strongest) and shot noise arising from their discrete nature. Recent work has argued that one can avoid sampling variance by comparing multiple tracers of different bias, and suppress shot noise by optimally weighting halos of different mass. Here we combine these ideas and investigate how well the signatures of non-Gaussian fluctuations in the primordial potential can be extracted from the two-point correlations of halos and dark matter. On the basis of large $N$-body simulations with local non-Gaussian initial conditions and their halo catalogs we perform a Fisher matrix analysis of the two-point statistics. Compared to the standard analysis, optimal weighting- and multiple-tracer techniques applied to halos can yield up to one order of magnitude improvements in $\fnl$-constraints, even if the underlying dark matter density field is not known. We compare our numerical results to the halo model and find satisfactory agreement. Forecasting the optimal $\fnl$-constraints that can be achieved with our methods when applied to existing and future survey data, we find that a survey of $50h^{-1}\mathrm{Gpc}^3$ volume resolving all halos down to $10^{11}\hMsun$ at $z=1$ will be able to obtain $\sigma_{\fnl}\sim1$ (68% cl), a factor of $\sim20$ improvement over the current limits. Decreasing the minimum mass of resolved halos, increasing the survey volume or obtaining the dark matter maps can further improve these limits, potentially reaching the level of $\sigma_{\fnl}\sim0.1$. (abridged)

In recent work (Seljak, Hamaus and Desjacques 2009) it was found that weighting central halo galaxies by halo mass can significantly suppress their stochasticity relative to the dark matter, well below the Poisson model expectation. In this paper we extend this study with the goal of finding the optimal mass-dependent halo weighting and use $N$-body simulations to perform a general analysis of halo stochasticity and its dependence on halo mass. We investigate the stochasticity matrix, defined as $C_{ij}\equiv<(\delta_i -b_i\delta_m)(\delta_j-b_j\delta_m)>$, where $\delta_m$ is the dark matter overdensity in Fourier space, $\delta_i$ the halo overdensity of the $i$-th halo mass bin and $b_i$ the halo bias. In contrast to the Poisson model predictions we detect nonvanishing correlations between different mass bins. We also find the diagonal terms to be sub-Poissonian for the highest-mass halos. The diagonalization of this matrix results in one large and one low eigenvalue, with the remaining eigenvalues close to the Poisson prediction $1/\bar{n}$, where $\bar{n}$ is the mean halo number density. The eigenmode with the lowest eigenvalue contains most of the information and the corresponding eigenvector provides an optimal weighting function to minimize the stochasticity between halos and dark matter. We find this optimal weighting function to match linear mass weighting at high masses, while at the low-mass end the weights approach a constant whose value depends on the low-mass cut in the halo mass function. Finally, we employ the halo model to derive the stochasticity matrix and the scale-dependent bias from an analytical perspective. It is remarkably successful in reproducing our numerical results and predicts that the stochasticity between halos and the dark matter can be reduced further when going to halo masses lower than we can resolve in current simulations.

Galaxy surveys are one of the most powerful means to extract the cosmological information and for a given volume the attainable precision is determined by the galaxy shot noise sigma_n^2 relative to the power spectrum P. It is generally assumed that shot noise is white and given by the inverse of the number density n. In this paper we argue one may be able to considerably improve upon this: in the halo picture of cosmological structure all of the dark matter is in halos of varying mass and galaxies are formed inside these halos, but for the dark matter mass and momentum conservation guarantee that nonlinear effects cannot develop a white noise in the dark matter power spectrum on large scales. This suggests that with a suitable weighting a similar effect may be achieved for galaxies, suppressing their shot noise. We explore this idea with N-body simulations by weighting central halo galaxies by halo mass and find that the resulting shot noise can be reduced dramatically relative to expectations, with a 10-30 suppression at the highest number density of n=4*10^-3 (Mpc/h)^3 resolved in our simulations. For specific applications other weighting schemes may achieve even better results and for n=3*10^-4(Mpc/h)^3 we can reduce sigma_n^2/P by up to a factor of 10 relative to uniform weighting. These results open up new opportunities to extract cosmological information in galaxy surveys, such as the recently proposed multi-tracer approach to cancel sampling variance, and may have important consequences for the planning of future redshift surveys. Taking full advantage of these findings may require better understanding of galaxy formation process to develop accurate tracers of the halo mass.

Sagittarius A*, the supermassive compact object at the center of the Galaxy, exhibits outbursts in the near infrared and X-ray domains. These flares are likely due to energetic events very close to the central object, on a scale of a few Schwarzschild radii. Optical interferometry will soon be able to provide astrometry with an accuracy of this order (~10 muas). In this article we use recent photometric near infrared data observed with the adaptive optics system NACO at the Very Large Telescope combined with simulations in order to deploy a method to test the nature of the flares and to predict the possible outcome of observations with the Very Large Telescope Interferometer. To accomplish this we implement a hot spot model and investigate its appearance for a remote observer in terms of light curves and centroid tracks, based on general relativistic ray tracing methods. First, we use a simplified model of a small steady source in order to investigate the relativistic effects qualitatively. A more realistic scenario is then being developed by fitting our model to existing flare data. While indications for the spin of the black hole and multiple images due to lensing effects are marginal in the light curves, astrometric measurements offer the possibility to reveal these high-order general relativistic effects. This study makes predictions on these astrometric measurements and leads us to the conclusion that future infrared interferometers will be able to detect proper motion of hot spots in the vicinity of Sagittarius A*.

We present the second-generation VLTI instrument GRAVITY, which currently is in the preliminary design phase. GRAVITY is specifically designed to observe highly relativistic motions of matter close to the event horizon of Sgr A*, the massive black hole at center of the Milky Way. We have identified the key design features needed to achieve this goal and present the resulting instrument concept. It includes an integrated optics, 4-telescope, dual feed beam combiner operated in a cryogenic vessel; near infrared wavefront sensing adaptive optics; fringe tracking on secondary sources within the field of view of the VLTI and a novel metrology concept. Simulations show that the planned design matches the scientific needs; in particular that 10 microarcsecond astrometry is feasible for a source with a magnitude of K=15 like Sgr A*, given the availability of suitable phase reference sources.