The cosmic web is one of the most striking features of the distribution of galaxies and dark matter on the largest scales in the Universe. It is composed of dense regions packed full of galaxies, long filamentary bridges, flattened sheets and vast low density voids. The study of the cosmic web has focused primarily on the identification of such features, and on understanding the environmental effects on galaxy formation and halo assembly. As such, a variety of different methods have been devised to classify the cosmic web -- depending on the data at hand, be it numerical simulations, large sky surveys or other. In this paper we bring twelve of these methods together and apply them to the same data set in order to understand how they compare. In general these cosmic web classifiers have been designed with different cosmological goals in mind, and to study different questions. Therefore one would not {\it a priori} expect agreement between different techniques however, many of these methods do converge on the identification of specific features. In this paper we study the agreements and disparities of the different methods. For example, each method finds that knots inhabit higher density regions than filaments, etc. and that voids have the lowest densities. For a given web environment, we find substantial overlap in the density range assigned by each web classification scheme. We also compare classifications on a halo-by-halo basis; for example, we find that 9 of 12 methods classify around a third of group-mass haloes (i.e. $M_{\rm halo}\sim10^{13.5}h^{-1}M_{\odot}$) as being in filaments. Lastly, so that any future cosmic web classification scheme can be compared to the 12 methods used here, we have made all the data used in this paper public.

We have carefully selected a sample of 60 galaxies that reside in the deepest underdensities of geometrically identified voids within the SDSS. HI imaging of 55 galaxies with the WSRT reveals morphological and kinematic signatures of ongoing interactions and gas accretion. We probe a total volume of 485 Mpc^3 within the voids, with an angular resolution of 8 kpc at an average distance of 85 Mpc. We reach column density sensitivities of 5 x 10^19 cm^-2, corresponding to an HI mass limit of 3 x 10^8 M_sun. We detect HI in 41 galaxies, with total masses ranging from 1.7 x 10^8 to 5.5 x 10^9 M_sun. The upper limits on the 14 non-detections are not inconsistent with their luminosities, given their expected HI mass to light ratios. We find that the void galaxies are generally gas rich, low luminosity, blue disk galaxies, with optical and HI properties that are not unusual for their luminosity and morphology. The sample spans a range of absolute magnitudes (-16.1 > M_r > -20.4) and colors (0.06 < g-r < 0.87), and includes disk and irregular galaxies. We also identify three as early type galaxies, all of which are not detected in HI. All galaxies have stellar masses less than 3 x 10^10 M_sun, and many have kinematic and morphological signs of ongoing gas accretion, suggesting that the void galaxy population is still in the process of assembling. The small scale clustering in the void, within 600 kpc and 200 km/s, is similar to that in higher density regions, and we identify 18 HI rich neighboring galaxies in the voids. Most are within 100 kpc and 100 km/s of the targeted galaxy, and we find no significant population of HI rich low luminosity galaxies filling the voids, contrary to what is predicted by simulations.

We investigate the ability of three reconstruction techniques to analyze and investigate weblike features and geometries in a discrete distribution of objects. The three methods are the linear Delaunay Tessellation Field Estimator (DTFE), its higher order equivalent Natural Neighbour Field Estimator (NNFE) and a version of Kriging interpolation adapted to the specific circumstances encountered in galaxy redshift surveys, the Natural Lognormal Kriging technique. DTFE and NNFE are based on the local geometry defined by the Voronoi and Delaunay tessellations of the galaxy distribution. The three reconstruction methods are analysed and compared using mock magnitude-limited and volume-limited SDSS redshift surveys, obtained on the basis of the Millennium simulation. We investigate error trends, biases and the topological structure of the resulting fields, concentrating on the void population identified by the Watershed Void Finder. Environmental effects are addressed by evaluating the density fields on a range of Gaussian filter scales. Comparison with the void population in the original simulation yields the fraction of false void mergers and false void splits. In most tests DTFE, NNFE and Kriging have largely similar density and topology error behaviour. Cosmetically, higher order NNFE and Kriging methods produce more visually appealing reconstructions. Quantitatively, however, DTFE performs better, even while computationally far less demanding. A successful recovery of the void population on small scales appears to be difficult, while the void recovery rate improves significantly on scales > 3 h-1Mpc. A study of small scale voids and the void galaxy population should therefore be restricted to the local Universe, out to at most 100 h-1Mpc.

The Void Galaxy Survey (VGS) is a multi-wavelength program to study $\sim$60 void galaxies. Each has been selected from the deepest interior regions of identified voids in the SDSS redshift survey on the basis of a unique geometric technique, with no a prior selection of intrinsic properties of the void galaxies. The project intends to study in detail the gas content, star formation history and stellar content, as well as kinematics and dynamics of void galaxies and their companions in a broad sample of void environments. It involves the HI imaging of the gas distribution in each of the VGS galaxies. Amongst its most tantalizing findings is the possible evidence for cold gas accretion in some of the most interesting objects, amongst which are a polar ring galaxy and a filamentary configuration of void galaxies. Here we shortly describe the scope of the VGS and the results of the full analysis of the pilot sample of 15 void galaxies.

We present the SpineWeb framework for the topological analysis of the Cosmic Web and the identification of its walls, filaments and cluster nodes. Based on the watershed segmentation of the cosmic density field, the SpineWeb method invokes the local adjacency properties of the boundaries between the watershed basins to trace the critical points in the density field and the separatrices defined by them. The separatrices are classified into walls and the spine, the network of filaments and nodes in the matter distribution. Testing the method with a heuristic Voronoi model yields outstanding results. Following the discussion of the test results, we apply the SpineWeb method to a set of cosmological N-body simulations. The latter illustrates the potential for studying the structure and dynamics of the Cosmic Web.

Void galaxies, residing within the deepest underdensities of the Cosmic Web, present an ideal population for the study of galaxy formation and evolution in an environment undisturbed by the complex processes modifying galaxies in clusters and groups, as well as provide an observational test for theories of cosmological structure formation. We have completed a pilot survey for the HI imaging aspects of a new Void Galaxy Survey (VGS), imaging 15 void galaxies in HI in local (d < 100 Mpc) voids. HI masses range from 3.5 x 10^8 to 3.8 x 10^9 M_sun, with one nondetection with an upper limit of 2.1 x 10^8 M_sun. Our galaxies were selected using a structural and geometric technique to produce a sample that is purely environmentally selected and uniformly represents the void galaxy population. In addition, we use a powerful new backend of the Westerbork Synthesis Radio Telescope that allows us to probe a large volume around each targeted galaxy, simultaneously providing an environmentally constrained sample of fore- and background control sample of galaxies while still resolving individual galaxy kinematics and detecting faint companions in HI. This small sample makes up a surprisingly interesting collection of perturbed and interacting galaxies, all with small stellar disks. Four galaxies have significantly perturbed HI disks, five have previously unidentified companions at distances ranging from 50 to 200 kpc, two are in interacting systems, and one was found to have a polar HI disk. Our initial findings suggest void galaxies are a gas-rich, dynamic population which present evidence of ongoing gas accretion, major and minor interactions, and filamentary alignment despite the surrounding underdense environment.

We study the topology of the Megaparsec Cosmic Web on the basis of the Alpha Shapes of the galaxy distribution. The simplicial complexes of the alpha shapes are used to determine the set of Betti numbers ($\beta_{\rm k},k=1,...,D$), which represent a complete characterization of the topology of a manifold. This forms a useful extension of the geometry and topology of the galaxy distribution by Minkowski functionals, of which three specify the geometrical structure of surfaces and one, the Euler characteristic, represents a key aspect of its topology. In order to develop an intuitive understanding for the relation between Betti numbers and the running $\alpha$ parameter of the alpha shapes, and thus in how far they may discriminate between different topologies, we study them within the context of simple heuristic Voronoi clustering models. These may be tuned to consist of a few or even only one specific morphological element of the Cosmic Web, ie. clusters, filaments or sheets.

We present a framework for the hierarchical identification and characterization of voids based on the Watershed Void Finder. The Hierarchical Void Finder is based on a generalization of the scale space of a density field invoked in order to trace the hierarchical nature and structure of cosmological voids. At each level of the hierarchy, the watershed transform is used to identify the voids at that particular scale. By identifying the overlapping regions between watershed basins in adjacent levels, the hierarchical void tree is constructed. Applications on a hierarchical Voronoi model and on a set of cosmological simulations illustrate its potential.

We review and discuss aspects of Cosmic Voids that form the background for our Void Galaxy Survey (see accompanying paper by Stanonik et al.). Following a sketch of the general characteristics of void formation and evolution, we describe the influence of the environment on their development and structure and the characteristic hierarchical buildup of the cosmic void population. In order to be able to study the resulting tenuous void substructure and the galaxies populating the interior of voids, we subsequently set out to describe our parameter free tessellation-based watershed void finding technique. It allows us to trace the outline, shape and size of voids in galaxy redshift surveys. The application of this technique enables us to find galaxies in the deepest troughs of the cosmic galaxy distribution, and has formed the basis of our void galaxy program.

We review the analysis of the Cosmic Web by means of an extensive toolset based on the use of Delaunay and Voronoi tessellations. The Cosmic Web is the salient and pervasive foamlike pattern in which matter has organized itself on scales of a few up to more than a hundred Megaparsec. First, we describe the Delaunay Tessellation Field Estimator (DTFE). The DTFE formalism is shown to recover the hierarchical nature and the anisotropic morphology of the cosmic matter distribution. The Multiscale Morphology Filter (MMF) uses the DTFE density field to extract the diverse morphological elements - filaments, sheets and clusters - on the basis of a ScaleSpace analysis which searches for these morphologies over a range of scales. Subsequently, we discuss the Watershed Voidfinder (WVF), which invokes the discrete watershed transform to identify voids in the cosmic matter distribution. The WVF is able to determine the location, size and shape of the voids. The watershed transform is also a key element in the SpineWeb analysis of the cosmic matter distribution. It allows the determination of the filamentary spine and connected walls in the cosmic matter density field through the identification of the singularities and corresponding separatrices. Finally, we describe the concept of Alphashapes for assessing the topology of the cosmic matter distribution.

In this review we discuss several aspects of Cosmic Voids. Voids are a major component of the large scale distribution of matter and galaxies in the Universe. They are of instrumental importance for understanding the emergence of the Cosmic Web. Their relatively simple shape and structure makes them into useful tools for extracting the value of a variety cosmic parameters, possibly including even that of the influence of dark energy. Perhaps most promising and challenging is the issue of the galaxies found within their realm. Not only does the pristine environment of voids provide a promising testing ground for assessing the role of environment on the formation and evolution of galaxies, the dearth of dwarf galaxies may even represent a serious challenge to the standard view of cosmic structure formation.

We have completed a pilot survey imaging 15 SDSS selected void galaxies in HI in local (d=50 to 100 Mpc) voids. This small sample makes up a surprisingly interesting collection of galaxies, consisting of galaxies with asymmetric and perturbed HI disks, previously unidentified companions, and ongoing interactions. One was found to have a polar HI disk with no stellar counterpart. While our small number statistics so far are limiting, results support past findings that most void galaxies are typically late type galaxies with gas rich disks and small scale clustering similar to field galaxies despite their large scale underdense environment.

We have found an isolated polar disk galaxy in what appears to be a cosmological wall situated between two voids. This void galaxy is unique as its polar disk was discovered serendipitously in an HI survey of SDSS void galaxies, with no optical counterpart to the HI polar disk. Yet the HI mass in the disk is comparable to the stellar mass in the galaxy. This suggests slow accretion of the HI material at a relatively recent time. There is also a hint of a warp in the outer parts of the HI disk. The central, stellar disk appears relatively blue, with faint near UV emission, and is oriented (roughly) parallel to the surrounding wall, implying gas accretion from out of the voids. The considerable gas mass and apparent lack of stars in the polar disk, coupled with the general underdensity of the environment, supports recent theories of cold flow accretion as an alternate formation mechanism for polar disk galaxies.

Despite a history that dates back at least a quarter of a century studies of voids in the large--scale structure of the Universe are bedevilled by a major problem: there exist a large number of quite different void--finding algorithms, a fact that has so far got in the way of groups comparing their results without worrying about whether such a comparison in fact makes sense. Because of the recent increased interest in voids, both in very large galaxy surveys and in detailed simulations of cosmic structure formation, this situation is very unfortunate. We here present the first systematic comparison study of thirteen different void finders constructed using particles, haloes, and semi--analytical model galaxies extracted from a subvolume of the Millennium simulation. The study includes many groups that have studied voids over the past decade. We show their results and discuss their differences and agreements. As it turns out, the basic results of the various methods agree very well with each other in that they all locate a major void near the centre of our volume. Voids have very underdense centres, reaching below 10 percent of the mean cosmic density. In addition, those void finders that allow for void galaxies show that those galaxies follow similar trends. For example, the overdensity of void galaxies brighter than $m_B = -20 $ is found to be smaller than about -0.8 by all our void finding algorithms.

We investigate the shapes and mutual alignment of voids in the large scale matter distribution of a LCDM cosmology simulation. The voids are identified using the novel WVF void finder technique. The identified voids are quite nonspherical and slightly prolate, with axis ratios in the order of c:b:a approx. 0.5:0.7:1. Their orientations are strongly correlated with significant alignments spanning scales >30 Mpc/h. We also find an intimate link between the cosmic tidal field and the void orientations. Over a very wide range of scales we find a coherent and strong alignment of the voids with the tidal field computed from the smoothed density distribution. This orientation-tide alignment remains significant on scales exceeding twice the typical void size, which shows that the long range external field is responsible for the alignment of the voids. This confirms the view that the large scale tidal force field is the main agent for the large scale spatial organization of the Cosmic Web.

On megaparsec scales the Universe is permeated by an intricate filigree of clusters, filaments, sheets and voids, the Cosmic Web. For the understanding of its dynamical and hierarchical history it is crucial to identify objectively its complex morphological components. One of the most characteristic aspects is that of the dominant underdense Voids, the product of a hierarchical process driven by the collapse of minor voids in addition to the merging of large ones. In this study we present an objective void finder technique which involves a minimum of assumptions about the scale, structure and shape of voids. Our void finding method, the Watershed Void Finder (WVF), is based upon the Watershed Transform, a well-known technique for the segmentation of images. Importantly, the technique has the potential to trace the existing manifestations of a void hierarchy. The basic watershed transform is augmented by a variety of correction procedures to remove spurious structure resulting from sampling noise. This study contains a detailed description of the WVF. We demonstrate how it is able to trace and identify, relatively parameter free, voids and their surrounding (filamentary and planar) boundaries. We test the technique on a set of Kinematic Voronoi models, heuristic spatial models for a cellular distribution of matter. Comparison of the WVF segmentations of low noise and high noise Voronoi models with the quantitatively known spatial characteristics of the intrinsic Voronoi tessellation shows that the size and shape of the voids are succesfully retrieved. WVF manages to even reproduce the full void size distribution function.

We present a model for the distribution of void sizes and its evolution within the context of hierarchical scenarios of gravitational structure formation. For a proper description of the hierarchical buildup of the system of voids in the matter distribution, not only the "void-in-void" problem should be taken into account, but also that of the "void-in-cloud" issue. Within the context of the excursion set formulation of an evolving void hierarchy is one involving a "two-barrier" excursion problem, unlike the "one-barrier" problem for the dark halo evolution. This leads to voids having a peaked size distribution at any cosmic epoch, centered on a characteristic void size that evolves self-similarly in time, this in distinct contrast to the distribution of virialized halo masses which do not have a small-scale cut-off.