Cosmic voids, the under-dense regions of the cosmic web, are widely used to constrain cosmology. Voids contain few, isolated galaxies, presumably expected to be less evolved and preserving memory of the pristine Universe. We use the cosmological hydrodynamical simulation Horizon-AGN coupled to the void finder {\sc \texttt{VIDE}} to investigate properties of galaxies in voids at z=0. We find that, closer to void centers, low-mass galaxies are more common than their massive counterparts. At fixed dark matter halo mass, they have smaller stellar masses than in denser regions. The star formation rate of void galaxies diminishes when approaching void centers, but their sSFR slightly increases, suggesting that void galaxies form stars more efficiently with respect to their stellar mass. We find that this can not only be attributed to the prevalence of low-mass galaxies. The inner region of voids also predominantly host low-mass BHs. However, the BH mass to galaxy mass ratios resemble those of the whole simulation at z=0. Our results suggest that even if the growth channels in cosmic voids are different than in denser environments, voids grow their galaxies and BHs in a similar way. While a large fraction of the BHs have low Eddington ratios, we find that 20\% could be observed as AGN with log10 L=41.5-42.5 erg/s in hard X-ray (2-10 keV). These results pave the way to future work with larger next-generation hydro simulations, aiming to confirm our findings and prepare the application on data from upcoming large surveys such as PFS, Euclid and WFIRST.

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.

Cosmic voids, the underdense regions in the universe, are particularly sensitive to diffuse density components such as cosmic neutrinos. This sensitivity is enhanced by the match between void sizes and the free-streaming scale of massive neutrinos. Using the massive neutrino simulations \texttt{MassiveNuS}, we investigate the effect of neutrino mass on dark matter halos as a function of environment. We find that the halo mass function depends strongly on neutrino mass and that this dependence is more pronounced in voids than in high-density environments. An observational program that measured the characteristic mass of the most massive halos in voids should be able to place novel constraints on the sum of the masses of neutrinos $\sum m_\nu$. The neutrino mass effect in the simulations is quite strong: In a 512$^3$ $h^{-3}$ Mpc$^3$ survey, the mean mass of the 1000 most massive halos in the void interiors is $(4.82 \pm 0.11) \times 10^{12} h^{-1}M_{\odot}$ for $\sum m_\nu = 0.6$ eV and $(8.21 \pm 0.13) \times 10^{12} h^{-1}M_{\odot}$ for $\sum m_\nu = 0.1$ eV. Subaru (SuMIRe), Euclid and WFIRST will have both spectroscopic and weak lensing surveys. Covering volumes at least 50 times larger than our simulations, they should be sensitive probes of neutrino mass through void substructure.

The Quijote simulations are a set of 43100 full N-body simulations spanning more than 7000 cosmological models in the $\{\Omega_{\rm m}, \Omega_{\rm b}, h, n_s, \sigma_8, M_\nu, w \}$ hyperplane. At a single redshift the simulations contain more than 8.5 trillions of particles over a combined volume of 43100 $(h^{-1}{\rm Gpc})^3$. Billions of dark matter halos and cosmic voids have been identified in the simulations, whose runs required more than 35 million core hours. The Quijote simulations have been designed for two main purposes: 1) to quantify the information content on cosmological observables, and 2) to provide enough data to train machine learning algorithms. In this paper we describe the simulations and show a few of their applications. We also release the Petabyte of data generated, comprising hundreds of thousands of simulation snapshots at multiple redshifts, halo and void catalogs, together with millions of summary statistics such as power spectra, bispectra, correlation functions, marked power spectra, and estimated probability density functions.

We present catalogues of cosmic voids identified in the distribution of Luminous Red Galaxies (LRGs) and Quasi Stellar Objects (QSOs) in the fourteenth data release (DR14) of the extended Baryon Oscillation Spectroscopic Survey (eBOSS). We perform a multivariate analysis to assess the level of contamination in these catalogues by spurious Poisson underdensities. We find that the LRG void catalogue is largely free from contamination but that the QSO catalogue may be heavily contaminated. We analyse the multipoles of the void-galaxy cross-correlation function in these catalogues to obtain constraints on the growth rate of structure around voids. We find a value of $\beta(z=0.703)=0.58^{+0.33}_{-0.28}$ for the LRG voids and $\beta(z=1.53)=0.15^{+0.13}_{-0.12}$ for the QSO voids.

Do void statistics contain information beyond the tracer 2-point correlation function? Yes! As we vary the sum of the neutrino masses, we find void statistics contain information absent when using just tracer 2-point statistics. Massive neutrinos uniquely affect cosmic voids. We explore their impact on void clustering using both the DEMNUni and MassiveNuS simulations. For voids, neutrino effects depend on the observed void tracers. As the neutrino mass increases, the number of small voids traced by cold dark matter particles increases and the number of large voids decreases. Surprisingly, when massive, highly biased, halos are used as tracers, we find the opposite effect. The scale at which voids cluster, as well as the void correlation, is similarly sensitive to the sum of neutrino masses and the tracers. This scale dependent trend is not due to simulation volume or halo density. The interplay of these signatures in the void abundance and clustering leaves a distinct fingerprint that could be detected with observations and potentially help break degeneracies between different cosmological parameters. This paper paves the way to exploit cosmic voids in future surveys to constrain the mass of neutrinos.

Two decades after its discovery, cosmic acceleration remains the most profound mystery in cosmology and arguably in all of physics. Either the Universe is dominated by a form of dark energy with exotic physical properties not predicted by standard model physics, or General Relativity is not an adequate description of gravity over cosmic distances. WFIRST emerged as a top priority of Astro2010 in part because of its ability to address the mystery of cosmic acceleration through both high precision measurements of the cosmic expansion history and the growth of cosmic structures with multiple and redundant probes. We illustrate in this white paper how mission design changes since Astro2010 have made WFIRST an even more powerful dark energy facility and have improved the ability of WFIRST to respond to changes in the experimental landscape. WFIRST is the space-based probe of DE the community needs in the mid-2020s.

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 Wide Field Infrared Survey Telescope (WFIRST) is a 2.4m space telescope with a 0.281 deg^2 field of view for near-IR imaging and slitless spectroscopy and a coronagraph designed for > 10^8 starlight suppresion. As background information for Astro2020 white papers, this article summarizes the current design and anticipated performance of WFIRST. While WFIRST does not have the UV imaging/spectroscopic capabilities of the Hubble Space Telescope, for wide field near-IR surveys WFIRST is hundreds of times more efficient. Some of the most ambitious multi-cycle HST Treasury programs could be executed as routine General Observer (GO) programs on WFIRST. The large area and time-domain surveys planned for the cosmology and exoplanet microlensing programs will produce extraordinarily rich data sets that enable an enormous range of Archival Research (AR) investigations. Requirements for the coronagraph are defined based on its status as a technology demonstration, but its expected performance will enable unprecedented observations of nearby giant exoplanets and circumstellar disks. WFIRST is currently in the Preliminary Design and Technology Completion phase (Phase B), on schedule for launch in 2025, with several of its critical components already in production.

SPHEREx is a proposed NASA MIDEX mission selected for Phase A study. SPHEREx would carry out the first all-sky spectral survey in the near infrared. At the end of its two-year mission, SPHEREx would obtain 0.75-to-5$\mu$m spectra of every 6.2 arcsec pixel on the sky, with spectral resolution R>35 and a 5-$\sigma$ sensitivity AB$>$19 per spectral/spatial resolution element. More details concerning SPHEREx are available at http://spherex.caltech.edu. The SPHEREx team has proposed three specific science investigations to be carried out with this unique data set: cosmic inflation, interstellar and circumstellar ices, and the extra-galactic background light. Though these three themes are undoubtedly compelling, they are far from exhausting the scientific output of SPHEREx. Indeed, SPHEREx would create a unique all-sky spectral database including spectra of very large numbers of astronomical and solar system targets, including both extended and diffuse sources. These spectra would enable a wide variety of investigations, and the SPHEREx team is dedicated to making the data available to the community to enable these investigations, which we refer to as Legacy Science. To that end, we have sponsored two workshops for the general scientific community to identify the most interesting Legacy Science themes and to ensure that the SPHEREx data products are responsive to their needs. In February of 2016, some 50 scientists from all fields met in Pasadena to develop these themes and to understand their implications for the SPHEREx mission. The 2016 workshop highlighted many synergies between SPHEREx and other contemporaneous astronomical missions, facilities, and databases. Consequently, in January 2018 we convened a second workshop at the Center for Astrophysics in Cambridge to focus specifically on these synergies. This white paper reports on the results of the 2018 SPHEREx workshop.

The fourth generation of the Sloan Digital Sky Survey (SDSS-IV) has been in operation since July 2014. This paper describes the second data release from this phase, and the fourteenth from SDSS overall (making this, Data Release Fourteen or DR14). This release makes public data taken by SDSS-IV in its first two years of operation (July 2014-2016). Like all previous SDSS releases, DR14 is cumulative, including the most recent reductions and calibrations of all data taken by SDSS since the first phase began operations in 2000. New in DR14 is the first public release of data from the extended Baryon Oscillation Spectroscopic Survey (eBOSS); the first data from the second phase of the Apache Point Observatory (APO) Galactic Evolution Experiment (APOGEE-2), including stellar parameter estimates from an innovative data driven machine learning algorithm known as "The Cannon"; and almost twice as many data cubes from the Mapping Nearby Galaxies at APO (MaNGA) survey as were in the previous release (N = 2812 in total). This paper describes the location and format of the publicly available data from SDSS-IV surveys. We provide references to the important technical papers describing how these data have been taken (both targeting and observation details) and processed for scientific use. The SDSS website (www.sdss.org) has been updated for this release, and provides links to data downloads, as well as tutorials and examples of data use. SDSS-IV is planning to continue to collect astronomical data until 2020, and will be followed by SDSS-V.

Cosmic acceleration is the most surprising cosmological discovery in many decades. Testing and distinguishing among possible explanations requires cosmological measurements of extremely high precision probing the full history of cosmic expansion and structure growth and, ideally, compare and contrast matter and relativistic tracers of the gravity potential. This program is one of the defining objectives of the Wide-Field Infrared Survey Telescope (WFIRST), as set forth in the New Worlds, New Horizons report (NWNH) in 2010. The WFIRST mission has the ability to improve these measurements by 1-2 orders of magnitude compared to the current state of the art, while simultaneously extending their redshift grasp, greatly improving control of systematic effects, and taking a unified approach to multiple probes that provide complementary physical information and cross-checks of cosmological results. We describe in this annual report the activities of the Science Investigation Team (SIT) "Cosmology with the High Latitude Survey (HLS)" during the year 2017. This team was selected by NASA in December 2015 in order to address the stringent challenges of the WFIRST dark energy (DE) program through the Project's formulation phase. This SIT has elected to jointly address Galaxy Redshift Survey, Weak Lensing and Cluster Growth and thus fully embrace the fact that the imaging and spectroscopic elements of the HLS will be realized as an integrated observing program, and they jointly impose requirements on performance and operations. WFIRST is designed to be able to deliver a definitive result on the origin of cosmic acceleration. It is not optimized for Figure of Merit sensitivity but for control of systematic uncertainties and for having multiple techniques each with multiple cross-checks. Our SIT work focuses on understanding the potential systematics in the WFIRST DE measurements.

The fourth generation of the Sloan Digital Sky Survey (SDSS-IV) began observations in July 2014. It pursues three core programs: APOGEE-2, MaNGA, and eBOSS. In addition, eBOSS contains two major subprograms: TDSS and SPIDERS. This paper describes the first data release from SDSS-IV, Data Release 13 (DR13), which contains new data, reanalysis of existing data sets and, like all SDSS data releases, is inclusive of previously released data. DR13 makes publicly available 1390 spatially resolved integral field unit observations of nearby galaxies from MaNGA, the first data released from this survey. It includes new observations from eBOSS, completing SEQUELS. In addition to targeting galaxies and quasars, SEQUELS also targeted variability-selected objects from TDSS and X-ray selected objects from SPIDERS. DR13 includes new reductions of the SDSS-III BOSS data, improving the spectrophotometric calibration and redshift classification. DR13 releases new reductions of the APOGEE-1 data from SDSS-III, with abundances of elements not previously included and improved stellar parameters for dwarf stars and cooler stars. For the SDSS imaging data, DR13 provides new, more robust and precise photometric calibrations. Several value-added catalogs are being released in tandem with DR13, in particular target catalogs relevant for eBOSS, TDSS, and SPIDERS, and an updated red-clump catalog for APOGEE. This paper describes the location and format of the data now publicly available, as well as providing references to the important technical papers that describe the targeting, observing, and data reduction. The SDSS website, http://www.sdss.org, provides links to the data, tutorials and examples of data access, and extensive documentation of the reduction and analysis procedures. DR13 is the first of a scheduled set that will contain new data and analyses from the planned ~6-year operations of SDSS-IV.

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 describe the Sloan Digital Sky Survey IV (SDSS-IV), a project encompassing three major spectroscopic programs. The Apache Point Observatory Galactic Evolution Experiment 2 (APOGEE-2) is observing hundreds of thousands of Milky Way stars at high resolution and high signal-to-noise ratio in the near-infrared. The Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey is obtaining spatially-resolved spectroscopy for thousands of nearby galaxies (median redshift of z = 0.03). The extended Baryon Oscillation Spectroscopic Survey (eBOSS) is mapping the galaxy, quasar, and neutral gas distributions between redshifts z = 0.6 and 3.5 to constrain cosmology using baryon acoustic oscillations, redshift space distortions, and the shape of the power spectrum. Within eBOSS, we are conducting two major subprograms: the SPectroscopic IDentification of eROSITA Sources (SPIDERS), investigating X-ray AGN and galaxies in X-ray clusters, and the Time Domain Spectroscopic Survey (TDSS), obtaining spectra of variable sources. All programs use the 2.5-meter Sloan Foundation Telescope at Apache Point Observatory; observations there began in Summer 2014. APOGEE-2 also operates a second near-infrared spectrograph at the 2.5-meter du Pont Telescope at Las Campanas Observatory, with observations beginning in early 2017. Observations at both facilities are scheduled to continue through 2020. In keeping with previous SDSS policy, SDSS-IV provides regularly scheduled public data releases; the first one, Data Release 13, was made available in July 2016.

Despite the ability of the cosmological concordance model ($\Lambda$CDM) to describe the cosmological observations exceedingly well, power law expansion of the Universe scale radius, $R(t)\propto t^n$, has been proposed as an alternative framework. We examine here these models, analyzing their ability to fit cosmological data using robust model comparison criteria. Type Ia supernovae (SNIa), baryonic acoustic oscillations (BAO) and acoustic scale information from the cosmic microwave background (CMB) have been used. We find that SNIa data either alone or combined with BAO can be well reproduced by both $\Lambda$CDM and power law expansion models with $n\sim 1.5$, while the constant expansion rate model $(n=1)$ is clearly disfavored. Allowing for some redshift evolution in the SNIa luminosity essentially removes any clear preference for a specific model. The CMB data are well known to provide the most stringent constraints on standard cosmological models, in particular, through the position of the first peak of the temperature angular power spectrum, corresponding to the sound horizon at recombination, a scale physically related to the BAO scale. Models with $n\geq 1$ lead to a divergence of the sound horizon and do not naturally provide the relevant scales for the BAO and the CMB. We retain an empirical footing to overcome this issue: we let the data choose the preferred values for these scales, while we recompute the ionization history in power law models, to obtain the distance to the CMB. In doing so, we find that the scale coming from the BAO data is not consistent with the observed position of the first peak of the CMB temperature angular power spectrum for any power law cosmology. Therefore, we conclude that when the three standard probes are combined, the $\Lambda$CDM model is very strongly favored over any of these alternative models, which are then essentially ruled out.

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.

We present a fast and robust alternative method to compute covariance matrix in case of cosmology studies. Our method is based on the jackknife resampling applied on simulation mock catalogues. Using a set of 600 BOSS DR11 mock catalogues as a reference, we find that the jackknife technique gives a similar galaxy clustering covariance matrix estimate by requiring a smaller number of mocks. A comparison of convergence rates show that $\sim$7 times fewer simulations are needed to get a similar accuracy on variance. We expect this technique to be applied in any analysis where the number of available N-body simulations is low.

The Extended Baryon Oscillation Spectroscopic Survey (eBOSS) will conduct novel cosmological observations using the BOSS spectrograph at Apache Point Observatory. Observations will be simultaneous with the Time Domain Spectroscopic Survey (TDSS) designed for variability studies and the Spectroscopic Identification of eROSITA Sources (SPIDERS) program designed for studies of X-ray sources. eBOSS will use four different tracers to measure the distance-redshift relation with baryon acoustic oscillations (BAO). Using more than 250,000 new, spectroscopically confirmed luminous red galaxies at a median redshift z=0.72, we project that eBOSS will yield measurements of $d_A(z)$ to an accuracy of 1.2% and measurements of H(z) to 2.1% when combined with the z>0.6 sample of BOSS galaxies. With ~195,000 new emission line galaxy redshifts, we expect BAO measurements of $d_A(z)$ to an accuracy of 3.1% and H(z) to 4.7% at an effective redshift of z= 0.87. A sample of more than 500,000 spectroscopically-confirmed quasars will provide the first BAO distance measurements over the redshift range 0.9<z<2.2, with expected precision of 2.8% and 4.2% on $d_A(z)$ and H(z), respectively. Finally, with 60,000 new quasars and re-observation of 60,000 quasars known from BOSS, we will obtain new Lyman-alpha forest measurements at redshifts z>2.1; these new data will enhance the precision of $d_A(z)$ and H(z) by a factor of 1.44 relative to BOSS. Furthermore, eBOSS will provide improved tests of General Relativity on cosmological scales through redshift-space distortion measurements, improved tests for non-Gaussianity in the primordial density field, and new constraints on the summed mass of all neutrino species. Here, we provide an overview of the cosmological goals, spectroscopic target sample, demonstration of spectral quality from early data, and projected cosmological constraints from eBOSS.

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.

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.

How do peculiar velocities affect observed voids? To answer this question we use the VIDE toolkit to identify voids in mock galaxy populations embedded within an N-body simulation both with and without peculiar velocities included. We compare the resulting void populations to assess the impact on void properties. We find that void abundances and spherically-averaged radial density profiles are mildly affected by peculiar velocities. However, peculiar velocities can distort by up to 10% the shapes for a particular subset of voids depending on the void size and density contrast, which can lead to increased variance in Alcock-Paczy\'nski test. We offer guidelines for performing optimal cuts on the void catalogue to reduce this variance by removing the most severely affected voids while preserving the unaffected ones. In addition, since this shape distortion is largely limited to the line of sight, we show that the void radii are only affected at the $\sim$ 10% level and the macrocenter positions at the $\sim$ 20% (even before performing cuts), meaning that cosmological probes based on the Integrated Sachs-Wolfe and gravitational lensing are not severely impacted by peculiar velocities.

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.

Cosmology has come a long way from being based on a small number of observations to being a data-driven precision science. We discuss the questions "What is observable?", "What in the Universe is knowable?" and "What are the fundamental limits to cosmological knowledge?". We then describe the methodology for investigation: theoretical hypotheses are used to model, predict and anticipate results; data is used to infer theory. We illustrate with concrete examples of principled analysis approaches from the study of cosmic microwave background anisotropies and surveys of large-scale structure, culminating in a summary of the highest precision probe to date of the physical origin of cosmic structures: the Planck 2013 constraints on primordial non-Gaussianity.

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 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 the 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 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 {\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.

We search for and characterize substructures in the projected distribution of galaxies observed in the wide field CCD images of the 77 nearby clusters of the WIde-field Nearby Galaxy-cluster Survey (WINGS). This sample is complete in X-ray flux in the redshift range 0.04<z<0.07. We search for substructures in WINGS clusters with DEDICA, an adaptive-kernel procedure. We test the procedure on Monte-Carlo simulations of the observed frames and determine the reliability for the detected structures. DEDICA identifies at least one reliable structure in the field of 55 clusters. 40 of these clusters have a total of 69 substructures at the same redshift of the cluster (redshift estimates of substructures are from color-magnitude diagrams). The fraction of clusters with subclusters (73%) is higher than in most studies. The presence of subclusters affects the relative luminosities of the brightest cluster galaxies (BCGs). Down to L ~ 10^11.2 L_Sun, our observed differential distribution of subcluster luminosities is consistent with the theoretical prediction of the differential mass function of substructures in cosmological simulations.

We measure the distribution of velocity dispersions of groups of galaxies identified in the UZC catalog; we use the distribution to derive the group mass function. We introduce a new method which makes efficient use of the entire magnitude limited catalog. Our determination of \nst includes a significant contribution from low luminosity systems that would be missing in a volume limited sample. We start from a model for the probability density function of the total number of group members and reproduce the observed distribution. We take several effects like the local fluctuations in volume density, limited sampling and group selection into account. We estimate the relation between total number of members, total luminosity and true velocity dispersion. We can then reproduce not only the observed distribution of \sv but also the distributions of the number of group members, the total groups luminosity, and virial mass. The best fit to the data in the true velocity dispersion range 100 km/s < \sigma < 750 km/s is a power law model with a slope -3.4 and a normalization (1.27 +/- 0.21) x 10^{-5} h^3 Mpc^{-3}. The predictions of CDM models are consistent with the mass function we derive from this distribution of velocity dispersions.

We apply a friends-of-friends algorithm to the combined UZC and SSRS2 redshift surveys to construct a catalog of 1168 groups of galaxies; 411 of these groups have 5 or more members within the redshift survey. The group catalog covers 4.69 steradians and all groups exceed a number density contrast threshold of 80. We demonstrate that the groups catalog is homogeneous across the two underlying redshift surveys; the catalog of groups and their members thus provides a basis for other statistical studies of the large-scale distribution of groups and their physical properties. The median physical properties of the groups are similar to those for groups derived from independent surveys including the ESO Key Programme and the Las Campanas Redshift Survey. We include tables of groups and their members.

In this paper we describe the Nearby Optical Galaxy (NOG) sample, which is a complete, distance-limited ($cz\leq$6000 km/s) and magnitude-limited (B$\leq$14) sample of $\sim$7000 optical galaxies. The sample covers 2/3 (8.27 sr) of the sky ($|b|>20^{\circ}$) and appears to have a good completeness in redshift (98%). We select the sample on the basis of homogenized corrected total blue magnitudes in order to minimize systematic effects in galaxy sampling. We identify the groups in this sample by means of both the hierarchical and the percolation {\it friends of friends} methods. The resulting catalogs of loose groups appear to be similar and are among the largest catalogs of groups presently available. Most of the NOG galaxies ($\sim$60%) are found to be members of galaxy pairs ($\sim$580 pairs for a total of $\sim$15% of objects) or groups with at least three members ($\sim$500 groups for a total of $\sim$45% of objects). About 40% of galaxies are left ungrouped (field galaxies). We illustrate the main features of the NOG galaxy distribution. Compared to previous optical and IRAS galaxy samples, the NOG provides a denser sampling of the galaxy distribution in the nearby universe. Given its large sky coverage, the identification of groups, and its high-density sampling, the NOG is suited for the analysis of the galaxy density field of the nearby universe, especially on small scales.

We discuss the 3D real-space reconstruction of the optical galaxy density field in the local Universe as derived from the galaxies of the Nearby Optical Galaxy (NOG) sample. NOG is a distance-limited (cz_{LG} < 6000 km/s) and magnitude--limited (B<14 mag) sample of 7076 optical galaxies which covers 2/3 (8.29 sr) of the sky (|b|>20). We have replaced ``true distances'' measurements for all the objects in order to correct for redshif distortions. Using homogenized photometric information for the whole sample, NOG is meant to be an approximation to a homogeneous all-sky 3D optically selected and statistically well-controlled galaxy sample that probes in great detail volumes of cosmological interest. Our goal is to construct a reliable, robust and unbiased field of density contrasts covering interesting regions of galaxy and mass overdensities of the local universe.

In order to map the galaxy density field in the local universe, we select the Nearby Optical Galaxy (NOG) sample, which is a distance-limited (cz < 6000 km/s) and magnitude--limited (B < 14 mag) sample of 7076 optical galaxies which covers 2/3 (8.29 sr) of the sky (|b|>20^{\circ}) and has a good completeness in redshift (98%). In order to trace the galaxy density field on small scales, we identify the NOG galaxy systems by means of both the hierarchical and the percolation (friends of friends) methods. The NOG provides high resolution in both spatial sampling of the nearby universe and morphological galaxy classification. The NOG is meant to be the first step towards the construction of a statistically well-controlled galaxy sample with homogenized photometric data covering most of the celestial sphere.

In this paper we identify objectively and analyze groups of galaxies in the recently completed ESP survey. We find 231 groups above the number overdensity threshold delta(rho)/rho=80 in the redshift range 5000<cz <60000 km/s. These groups contain 1250 members, 40.5% of the 3085 ESP galaxies within the same redshift range. The median velocity dispersion (corrected for measurement errors and computed at the redshift of the group) is sigma_{ESP,median} = 194 km/s. We show that our result is reliable in spite of the particular geometry of the ESP survey (two rows of tangent circular fields of radius 15 arcmin), which causes most systems to be only partially surveyed. In general, we find that the properties of ESP groups are consistent with those of groups in shallower (and wider) catalogs (e.g. CfA2N and SSRS2). As in shallower catalogs, ESP groups trace very well the geometry of the large scale structure. Our results are of particular interest because the depth of the ESP survey allows us to sample group properties over a large number of structures. We also compare luminosity function and spectral properties of galaxies that are members of groups with those of isolated galaxies. We find that galaxies in groups have a brighter M* with respect to non--member galaxies; the slope \alpha is the same, within the errors, in the two cases. We find that 34% of ESP galaxies with detectable emission lines are members of groups. The fraction of galaxies without detectable emission lines in groups is significantly higher: 45%. More generally, we find a gradual decrease of the fraction of emission line galaxies among members of systems of increasing richness. This result confirms that the morphology-density relation found for clusters also extends toward systems of lower density.

The "algorithm driven by the density estimate for the identification of clusters" (DEDICA, Pisani 1993, 1996) is applied to the A3558 cluster complex in order to find substructures. This complex, located at the center of the Shapley Concentration supercluster, is a chain formed by the ACO clusters A3556, A3558 and A3562 and the two poor clusters SC 1327-312 and SC 1329-313. We find a large number of clumps, indicating that strong dynamical processes are active. In particular, it is necessary to use a fully three-dimensional sample(i.e. using the galaxy velocity as third coordinate) in order to recover also the clumps superimposed along the line of sight. Even if a great number of detected substructures were already found in a previous analysis (Bardelli et al. 1998), this method is more efficient and faster when compared with the use of a wide battery of tests and permits the direct estimate of the detection significance. Almost all subclusters previously detected by the wavelet analyses found in the literature are recognized by DEDICA. On the basis of the substructure analysis, we also briefly discuss the origin of the A3558 complex by comparing two hypotheses: 1) the structure is a cluster-cluster collision seen just after the first core-core encounter; 2) this complex is the result of a series of incoherent group-group and cluster-group mergings, focused in that region by the presence of the surrounding supercluster. We studied the fraction of blue galaxies in the detected substructures and found that the bluest groups reside between A3562 and A3558, i.e. in the expected position in the scenario of the cluster-cluster collision.

A general method is described for detecting and analysing galaxy systems. The multivariate geometrical structure of the sample is studied by using an extension of the method which we introduced in a previous paper. The method is based on an estimate of the probability density underlying a data sample. The density is estimated by using an iterative and adaptive kernel estimator. The used kernels have spherical symmetry, however we describe a method in order to estimate the locally optimal shape of the kernels. We use the results of the geometrical structure analysis in order to study the effects that is has on the cluster parameter estimate. This suggests a possible way to distinguish between structure and substructure within a sample. The method is tested by using simulated numerical models and applied to two galaxy samples taken from the literature. The results obtained for the Coma cluster suggest a core-halo structure formed by a large number of geometrically independent systems. A different conclusion is suggested by the results for the Cancer cluster indicating the presence of at least two independent structures both containing substructure. The dynamical consequences of the results obtained from the geometrical analysis will be described in a later paper. Further applications of the method are suggested and are currently in progress.