We present the DTFE public software, a code for reconstructing fields from a discrete set of samples/measurements using the maximum of information contained in the point distribution. The code is written in C++ using the CGAL library and is parallelized using OpenMP. The software was designed for the analysis of cosmological data but can be used in other fields where one must interpolate quantities given at a discrete point set. The software comes with a wide suite of options to facilitate the analysis of 2- and 3-dimensional data and of both numerical simulations and galaxy redshift surveys. For comparison purposes, the code also implements the TSC and SPH grid interpolation methods. The code comes with an extensive user guide detailing the program options, examples and the inner workings of the code. The DTFE public software and further information can be found at http://www.astro.rug.nl/~voronoi/DTFE/dtfe.html .

The topology and geometry of random fields - in terms of the Euler characteristic and the Minkowski functionals - has received a lot of attention in the context of the Cosmic Microwave Background (CMB), as the detection of primordial non-Gaussianities would form a valuable clue on the physics of the early Universe. The virtue of both the Euler characteristic and the Minkowski functionals in general, lies in the fact that there exist closed form expressions for their expectation values for Gaussian random fields. However, the Euler characteristic and Minkowski functionals are summarizing characteristics of topology and geometry. Considerably more topological information is contained in the homology of the random field, as it completely describes the creation, merging and disappearance of topological features in superlevel set filtrations. In the present study we extend the topological analysis of the superlevel set filtrations of two-dimensional Gaussian random fields by analysing the statistical properties of the Betti numbers - counting the number of connected components and loops - and the persistence diagrams - describing the creation and mergers of homological features. Using the link between homology and the critical points of a function - as illustrated by the Morse-Smale complex - we derive a one-parameter fitting formula for the expectation value of the Betti numbers and forward this formalism to the persistent diagrams. We, moreover, numerically demonstrate the sensitivity of the Betti numbers and persistence diagrams to the presence of non-Gaussianities.

We present a self-consistent Bayesian formalism to sample the primordial density fields compatible with a set of dark matter density tracers after cosmic evolution observed in redshift space. Previous works on density reconstruction did not self-consistently consider redshift space distortions or included an additional iterative distortion correction step. We present here the analytic solution of coherent flows within a Hamiltonian Monte Carlo posterior sampling of the primordial density field. We test our method within the Zel'dovich approximation, presenting also an analytic solution including tidal fields and spherical collapse on small scales using augmented Lagrangian perturbation theory. Our resulting reconstructed fields are isotropic and their power spectra are unbiased compared to the true one defined by our mock observations. Novel algorithmic implementations are introduced regarding the mass assignment kernels when defining the dark matter density field and optimization of the time step in the Hamiltonian equations of motions. Our algorithm, dubbed barcode, promises to be specially suited for analysis of the dark matter cosmic web down to scales of a few Megaparsecs. This large scale structure is implied by the observed spatial distribution of galaxy clusters --- such as obtained from X-ray, SZ or weak lensing surveys --- as well as that of the intergalactic medium sampled by the Lyman alpha forest or perhaps even by deep hydrogen intensity mapping. In these cases, virialized motions are negligible, and the tracers cannot be modeled as point-like objects. It could be used in all of these contexts as a baryon acoustic oscillation reconstruction algorithm.

We introduce the 4-metre Multi-Object Spectroscopic Telescope (4MOST), a new high-multiplex, wide-field spectroscopic survey facility under development for the four-metre-class Visible and Infrared Survey Telescope for Astronomy (VISTA) at Paranal. Its key specifications are: a large field of view (FoV) of 4.2 square degrees and a high multiplex capability, with 1624 fibres feeding two low-resolution spectrographs ($R = \lambda/\Delta\lambda \sim 6500$), and 812 fibres transferring light to the high-resolution spectrograph ($R \sim 20\,000$). After a description of the instrument and its expected performance, a short overview is given of its operational scheme and planned 4MOST Consortium science; these aspects are covered in more detail in other articles in this edition of The Messenger. Finally, the processes, schedules, and policies concerning the selection of ESO Community Surveys are presented, commencing with a singular opportunity to submit Letters of Intent for Public Surveys during the first five years of 4MOST operations.

We investigate the alignment of galaxies and haloes relative to cosmic web filaments using the EAGLE hydrodynamical simulation. We identify filaments by applying the NEXUS+ method to the mass distribution and the Bisous formalism to the galaxy distribution. Both web finders return similar filamentary structures that are well aligned and that contain comparable galaxy populations. EAGLE haloes have an identical spin alignment with filaments as their counterparts in dark matter only simulations: a complex mass dependent trend with low mass haloes spinning preferentially parallel to and high mass haloes spinning preferentially perpendicular to filaments. In contrast, galaxy spins do not show such a spin transition and have a propensity for perpendicular alignments at all masses, with the degree of alignment being largest for massive galaxies. This result is valid for both NEXUS+ and Bisous filaments. When splitting by morphology, we find that elliptical galaxies show a stronger orthogonal spin--filament alignment than spiral galaxies of similar mass. The same is true of their haloes, with the host haloes of elliptical galaxies having a larger degree of orthogonal alignment than the host haloes of spirals. Due to the misalignment between galaxy shape and spin, galaxy minor axes are oriented differently with filaments than galaxy spins. We find that the galaxies whose minor axis is perpendicular to a filament are much better aligned with their host haloes. This suggests that many of the same physical processes determine both the galaxy--filament and the galaxy--halo alignments.

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.

This study presents a numerical analysis of the topology of a set of cosmologically interesting three-dimensional Gaussian random fields in terms of their Betti numbers $\beta_0$, $\beta_1$ and $\beta_2$. We show that Betti numbers entail a considerably richer characterization of the topology of the primordial density field. Of particular interest is that Betti numbers specify which topological features - islands, cavities or tunnels - define its spatial structure. A principal characteristic of Gaussian fields is that the three Betti numbers dominate the topology at different density ranges. At extreme density levels, the topology is dominated by a single class of features. At low levels this is a \emph{Swiss-cheeselike} topology, dominated by isolated cavities, at high levels a predominantly \emph{Meatball-like} topology of isolated objects. At moderate density levels, two Betti number define a more \emph{Sponge-like} topology. At mean density, the topology even needs three Betti numbers, quantifying a field consisting of several disconnected complexes, not of one connected and percolating overdensity. A {\it second} important aspect of Betti number statistics is that they are sensitive to the power spectrum. It reveals a monotonic trend in which at a moderate density range a lower spectral index corresponds to a considerably higher (relative) population of cavities and islands. We also assess the level of complementary information that Betti numbers represent, in addition to conventional measures such as Minkowski functionals. To this end, we include an extensive description of the Gaussian Kinematic Formula (GKF), which represents a major theoretical underpinning for this discussion.

We study the topology generated by the temperature fluctuations of the Cosmic Microwave Background (CMB) radiation, as quantified by the number of components and holes, formally given by the Betti numbers, in the growing excursion sets. We compare CMB maps observed by the Planck satellite with a thousand simulated maps generated according to the LCDM paradigm with Gaussian distributed fluctuations. The survey of the CMB over $\mathbb{S}^2$ is incomplete due to obfuscation effects by bright point sources and other extended foreground objects like our own galaxy. To deal with such situations, where analysis in the presence of "masks" is of importance, we introduce the concept of relative homology. The parametric $\chi^2$-test shows differences between observations and simulations, yielding $p$-values at per-cent to less than per-mil levels roughly between 2 to 7 degrees. The highest observed deviation for $b_0$ and $b_1$ is approximately between $3\sigma$-4$\sigma$ at scales of 3 to 7 degrees. There are reports of mildly unusual behaviour of the Euler characteristic at 3.66 degrees in the literature, computed from independent measurements of the CMB temperature fluctuations by Planck's predecessor WMAP satellite. The mildly anomalous behaviour of Euler characteristic is related to the strongly anomalous behaviour of components and holes. These are also the scales at which the observed maps exhibit low variance compared to the simulations. Non-parametric tests show even stronger differences at almost all scales. Regardless, beyond the trivial possibility that this may still be a manifestation of an extreme Gaussian case, these observations, along with the super-horizon scales involved, may motivate to look at primordial non-Gaussianity. Alternative scenarios worth exploring may be models with non-trivial topology.

We present a new formalism for studying the topology of HII regions during the Epoch of Reionisation, based on persistent homology theory. With persistent homology, it is possible to follow the evolution of topological features over time. We introduce the notion of a persistence field as a statistical summary of persistence data and we show how these fields can be used to identify different stages of reionisation. We identify two new stages common to all bubble ionisation scenarios. Following an initial pre-overlap and subsequent overlap stage, the topology is first dominated by neutral filaments (filament stage) and then by enclosed patches of neutral hydrogen undergoing outside-in ionisation (patch stage). We study how these stages are affected by the degree of galaxy clustering. We also show how persistence fields can be used to study other properties of the ionisation topology, such as the bubble size distribution and the fractal-like topology of the largest ionised region.

We investigate the alignment of haloes with the filaments of the cosmic web using an unprecedently large sample of dark matter haloes taken from the P-Millennium $\Lambda$CDM cosmological N-body simulation. We use the state-of-the-art NEXUS morphological formalism which, due to its multiscale nature, simultaneously identifies structures at all scales. We find strong and highly significant alignments, with both the major axis of haloes and their peculiar velocity tending to orient along the filament. However, the spin - filament alignment displays a more complex trend changing from preferentially parallel at low masses to preferentially perpendicular at high masses. This "spin flip" occurs at an average mass of $5\times10^{11}~h^{-1}M_\odot$. This mass increases with increasing filament diameter, varying by more than an order of magnitude between the thinnest and thickest filament samples. We also find that the inner parts of haloes have a spin flip mass that is several times smaller than that of the halo as a whole. These results confirm that recent accretion is responsible for the complex behaviour of the halo spin - filament alignment. Low-mass haloes mainly accrete mass along directions perpendicular to their host filament and thus their spins tend to be oriented along the filaments. In contrast, high-mass haloes mainly accrete along their host filaments and have their spins preferentially perpendicular to them. Furthermore, haloes located in thinner filaments are more likely to accrete along their host filaments than haloes of the same mass located in thicker filaments.

For over twenty years, the term 'cosmic web' has guided our understanding of the large-scale arrangement of matter in the cosmos, accurately evoking the concept of a network of galaxies linked by filaments. But the physical correspondence between the cosmic web and structural-engineering or textile 'spiderwebs' is even deeper than previously known, and extends to origami tessellations as well. Here we explain that in a good structure-formation approximation known as the adhesion model, threads of the cosmic web form a spiderweb, i.e. can be strung up to be entirely in tension. The correspondence is exact if nodes sampling voids are included, and if structure is excluded within collapsed regions (walls, filaments and haloes), where dark-matter multistreaming and baryonic physics affect the structure. We also suggest how concepts arising from this link might be used to test cosmological models: for example, to test for large-scale anisotropy and rotational flows in the cosmos.

We present a general formalism for identifying the caustic structure of an evolving mass distribution in an arbitrary dimensional space. For the class of Hamiltonian fluids the identification corresponds to the classification of singularities in Lagrangian catastrophe theory. Based on this we develop a theoretical framework for the formation of the cosmic web, and specifically those aspects that characterize its unique nature: its complex topological connectivity and multiscale spinal structure of sheetlike membranes, elongated filaments and compact cluster nodes. The present work represents an extension of the work by Arnol'd et al., who classified the caustics for the 1- and 2-dimensional Zel'dovich approximation. His seminal work established the role of emerging singularities in the formation of nonlinear structures in the universe. At the transition from the linear to nonlinear structure evolution, the first complex features emerge at locations where different fluid elements cross to establish multistream regions. The classification and characterization of these mass element foldings can be encapsulated in caustic conditions on the eigenvalue and eigenvector fields of the deformation tensor field. We introduce an alternative and transparent proof for Lagrangian catastrophe theory, and derive the caustic conditions for general Lagrangian fluids, with arbitrary dynamics, including dissipative terms and vorticity. The new proof allows us to describe the full 3-dimensional complexity of the gravitationally evolving cosmic matter field. One of our key findings is the significance of the eigenvector field of the deformation field for outlining the spatial structure of the caustic skeleton. We consider the caustic conditions for the 3-dimensional Zel'dovich approximation, extending earlier work on those for 1- and 2-dimensional fluids towards the full spatial richness of the cosmic web.

Future observations of cosmic microwave background (CMB) polarisation have the potential to answer some of the most fundamental questions of modern physics and cosmology. In this paper, we list the requirements for a future CMB polarisation survey addressing these scientific objectives, and discuss the design drivers of the CORE space mission proposed to ESA in answer to the "M5" call for a medium-sized mission. The rationale and options, and the methodologies used to assess the mission's performance, are of interest to other future CMB mission design studies. CORE is designed as a near-ultimate CMB polarisation mission which, for optimal complementarity with ground-based observations, will perform the observations that are known to be essential to CMB polarisation scienceand cannot be obtained by any other means than a dedicated space mission.

We discuss the potential of a next generation space-borne Cosmic Microwave Background (CMB) experiment for studies of extragalactic sources. Our analysis has particular bearing on the definition of the future space project, CORE, that has been submitted in response to ESA's call for a Medium-size mission opportunity as the successor of the Planck satellite. Even though the effective telescope size will be somewhat smaller than that of Planck, CORE will have a considerably better angular resolution at its highest frequencies, since, in contrast with Planck, it will be diffraction limited at all frequencies. The improved resolution implies a considerable decrease of the source confusion, i.e. substantially fainter detection limits. In particular, CORE will detect thousands of strongly lensed high-z galaxies distributed over the full sky. The extreme brightness of these galaxies will make it possible to study them, via follow-up observations, in extraordinary detail. Also, the CORE resolution matches the typical sizes of high-z galaxy proto-clusters much better than the Planck resolution, resulting in a much higher detection efficiency; these objects will be caught in an evolutionary phase beyond the reach of surveys in other wavebands. Furthermore, CORE will provide unique information on the evolution of the star formation in virialized groups and clusters of galaxies up to the highest possible redshifts. Finally, thanks to its very high sensitivity, CORE will detect the polarized emission of thousands of radio sources and, for the first time, of dusty galaxies, at mm and sub-mm wavelengths, respectively.

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 forecast the scientific capabilities to improve our understanding of cosmic inflation of CORE, a proposed CMB space satellite submitted in response to the ESA fifth call for a medium-size mission opportunity. The CORE satellite will map the CMB anisotropies in temperature and polarization in 19 frequency channels spanning the range 60-600 GHz. CORE will have an aggregate noise sensitivity of $1.7 \mu$K$\cdot \,$arcmin and an angular resolution of 5' at 200 GHz. We explore the impact of telescope size and noise sensitivity on the inflation science return by making forecasts for several instrumental configurations. This study assumes that the lower and higher frequency channels suffice to remove foreground contaminations and complements other related studies of component separation and systematic effects, which will be reported in other papers of the series "Exploring Cosmic Origins with CORE." We forecast the capability to determine key inflationary parameters, to lower the detection limit for the tensor-to-scalar ratio down to the $10^{-3}$ level, to chart the landscape of single field slow-roll inflationary models, to constrain the epoch of reheating, thus connecting inflation to the standard radiation-matter dominated Big Bang era, to reconstruct the primordial power spectrum, to constrain the contribution from isocurvature perturbations to the $10^{-3}$ level, to improve constraints on the cosmic string tension to a level below the presumptive GUT scale, and to improve the current measurements of primordial non-Gaussianities down to the $f_{NL}^{\rm local} < 1$ level. For all the models explored, CORE alone will improve significantly on the present constraints on the physics of inflation. Its capabilities will be further enhanced by combining with complementary future cosmological observations.

We forecast the main cosmological parameter constraints achievable with the CORE space mission which is dedicated to mapping the polarisation of the Cosmic Microwave Background (CMB). CORE was recently submitted in response to ESA's fifth call for medium-sized mission proposals (M5). Here we report the results from our pre-submission study of the impact of various instrumental options, in particular the telescope size and sensitivity level, and review the great, transformative potential of the mission as proposed. Specifically, we assess the impact on a broad range of fundamental parameters of our Universe as a function of the expected CMB characteristics, with other papers in the series focusing on controlling astrophysical and instrumental residual systematics. In this paper, we assume that only a few central CORE frequency channels are usable for our purpose, all others being devoted to the cleaning of astrophysical contaminants. On the theoretical side, we assume LCDM as our general framework and quantify the improvement provided by CORE over the current constraints from the Planck 2015 release. We also study the joint sensitivity of CORE and of future Baryon Acoustic Oscillation and Large Scale Structure experiments like DESI and Euclid. Specific constraints on the physics of inflation are presented in another paper of the series. In addition to the six parameters of the base LCDM, which describe the matter content of a spatially flat universe with adiabatic and scalar primordial fluctuations from inflation, we derive the precision achievable on parameters like those describing curvature, neutrino physics, extra light relics, primordial helium abundance, dark matter annihilation, recombination physics, variation of fundamental constants, dark energy, modified gravity, reionization and cosmic birefringence. (ABRIDGED)

We examine the cosmological constraints that can be achieved with a galaxy cluster survey with the future CORE space mission. Using realistic simulations of the millimeter sky, produced with the latest version of the Planck Sky Model, we characterize the CORE cluster catalogues as a function of the main mission performance parameters. We pay particular attention to telescope size, key to improved angular resolution, and discuss the comparison and the complementarity of CORE with ambitious future ground-based CMB experiments that could be deployed in the next decade. A possible CORE mission concept with a 150 cm diameter primary mirror can detect of the order of 50,000 clusters through the thermal Sunyaev-Zeldovich effect (SZE). The total yield increases (decreases) by 25% when increasing (decreasing) the mirror diameter by 30 cm. The 150 cm telescope configuration will detect the most massive clusters ($>10^{14}\, M_\odot$) at redshift $z>1.5$ over the whole sky, although the exact number above this redshift is tied to the uncertain evolution of the cluster SZE flux-mass relation; assuming self-similar evolution, CORE will detect $\sim 500$ clusters at redshift $z>1.5$. This changes to 800 (200) when increasing (decreasing) the mirror size by 30 cm. CORE will be able to measure individual cluster halo masses through lensing of the cosmic microwave background anisotropies with a 1-$\sigma$ sensitivity of $4\times10^{14} M_\odot$, for a 120 cm aperture telescope, and $10^{14} M_\odot$ for a 180 cm one. [abridged]

We introduce a multiscale topological description of the Megaparsec weblike cosmic matter distribution. Betti numbers and topological persistence offer a powerful means of describing the rich connectivity structure of the cosmic web and of its multiscale arrangement of matter and galaxies. Emanating from algebraic topology and Morse theory, Betti numbers and persistence diagrams represent an extension and deepening of the cosmologically familiar topological genus measure, and the related geometric Minkowski functionals. In addition to a description of the mathematical background, this study presents the computational procedure for computing Betti numbers and persistence diagrams for density field filtrations. The field may be computed starting from a discrete spatial distribution of galaxies or simulation particles. The main emphasis of this study concerns an extensive and systematic exploration of the imprint of different weblike morphologies and different levels of multiscale clustering in the corresponding computed Betti numbers and persistence diagrams. To this end, we use Voronoi clustering models as templates for a rich variety of weblike configurations, and the fractal-like Soneira-Peebles models exemplify a range of multiscale configurations. We have identified the clear imprint of cluster nodes, filaments, walls, and voids in persistence diagrams, along with that of the nested hierarchy of structures in multiscale point distributions. We conclude by outlining the potential of persistent topology for understanding the connectivity structure of the cosmic web, in large simulations of cosmic structure formation and in the challenging context of the observed galaxy distribution in large galaxy surveys.

Voids form a prominent aspect of the Megaparsec distribution of galaxies and matter. Not only do they represent a key constituent of the Cosmic Web, they also are one of the cleanest probes and measures of global cosmological parameters. The shape and evolution of voids are highly sensitive to the nature of dark energy, while their substructure and galaxy population provides a direct key to the nature of dark matter. Also, the pristine environment of void interiors is an important testing ground for our understanding of environmental influences on galaxy formation and evolution. In this paper, we review the key aspects of the structure and dynamics of voids, with a particular focus on the hierarchical evolution of the void population. We demonstrate how the rich structural pattern of the Cosmic Web is related to the complex evolution and buildup of voids.

In 1970 Zel'dovich published a far-reaching paper presenting a simple equation describing the nonlinear growth of primordial density inhomogeneities. The equation was remarkably successful in explaining the large scale structure in the Universe that we observe: a Universe in which the structure appears to be delineated by filaments and clusters of galaxies surrounding huge void regions. In order to concretise this impression it is necessary to define these structural elements through formal techniques with which we can compare the Zel'dovich model and N-body simulations with the observational data. We present an overview of recent efforts to identify voids, filaments and clusters in both the observed galaxy distribution and in numerical simulations of structure formation. We focus, in particular, on methods that involve no fine-tuning of parameters and that handle scale dependence automatically. It is important that these techniques should result in finding structures that relate directly to the dynamical mechanism of structure formation.

We describe the Bayesian BARCODE formalism that has been designed towards the reconstruction of the Cosmic Web in a given volume on the basis of the sampled galaxy cluster distribution. Based on the realization that the massive compact clusters are responsible for the major share of the large scale tidal force field shaping the anisotropic and in particular filamentary features in the Cosmic Web. Given the nonlinearity of the constraints imposed by the cluster configurations, we resort to a state-of-the-art constrained reconstruction technique to find a proper statistically sampled realization of the original initial density and velocity field in the same cosmic region. Ultimately, the subsequent gravitational evolution of these initial conditions towards the implied Cosmic Web configuration can be followed on the basis of a proper analytical model or an N-body computer simulation. The BARCODE formalism includes an implicit treatment for redshift space distortions. This enables a direct reconstruction on the basis of observational data, without the need for a correction of redshift space artifacts. In this contribution we provide a general overview of the the Cosmic Web connection with clusters and a description of the Bayesian BARCODE formalism. We conclude with a presentation of its successful workings with respect to test runs based on a simulated large scale matter distribution, in physical space as well as in redshift space.

The Zeldovich approximation (ZA) predicts the formation of a web of singularities. While these singularities may only exist in the most formal interpretation of the ZA, they provide a powerful tool for the analysis of initial conditions. We present a novel method to find the skeleton of the resulting cosmic web based on singularities in the primordial deformation tensor and its higher order derivatives. We show that the A_3-lines predict the formation of filaments in a two-dimensional model. We continue with applications of the adhesion model to visualise structures in the local (z < 0.03) universe.

We analyze photometry from deep B-band images of 59 void galaxies in the Void Galaxy Survey (VGS), together with their near-infrared 3.6$\mu$m and 4.5$\mu$m Spitzer photometry. The VGS galaxies constitute a sample of void galaxies that were selected by a geometric-topological procedure from the SDSS DR7 data release, and which populate the deep interior of voids. Our void galaxies span a range of absolute B-magnitude from $\rm{M_B=-15.5}$ to $\rm{M_B=-20}$, while at the 3.6$\mu$m band their magnitudes range from $\rm{M_{3.6}=-18}$ to $\rm{M_{3.6}=-24}$. Their B-[3.6] colour and structural parameters indicate these are star forming galaxies. A good reflection of the old stellar population, the near-infrared band photometry also provide a robust estimate of the stellar mass, which for the VGS galaxies we confirm to be smaller than $3 \times 10^{10}$ M$_\odot$. In terms of the structural parameters and morphology, our findings align with other studies in that our VGS galaxy sample consists mostly of small late-type galaxies. Most of them are similar to Sd-Sm galaxies, although a few are irregularly shaped galaxies. The sample even includes two early-type galaxies, one of which is an AGN. Their S\'{e}rsic indices are nearly all smaller than $n=2$ in both bands and they also have small half-light radii. In all, we conclude that the principal impact of the void environment on the galaxies populating them mostly concerns their low stellar mass and small size.

We study the star formation properties of 59 void galaxies as part of the Void Galaxy Survey (VGS). Current star formation rates are derived from $\rm{H\alpha}$ and recent star formation rates from near-UV imaging. In addition, infrared 3.4 $\rm{\mu m}$, 4.6 $\rm{\mu m}$, 12 $\rm{\mu m}$ and 22 $\rm{\mu m}$ WISE emission is used as star formation and mass indicator. Infrared and optical colours show that the VGS sample displays a wide range of dust and metallicity properties. We combine these measurements with stellar and HI masses to measure the specific SFRs ($\rm{SFR/M_{*}}$) and star formation efficiencies ($\rm{SFR/M_{HI}}$). We compare the star formation properties of our sample with galaxies in the more moderate density regions of the cosmic web, 'the field'. We find that specific SFRs of the VGS galaxies as a function of stellar and HI mass are similar to those of the galaxies in these field regions. Their $\rm{SFR\alpha}$ is slightly elevated than the galaxies in the field for a given total HI mass. In the global star formation picture presented by Kennicutt-Schmidt, VGS galaxies fall into the regime of low average star formation and correspondingly low HI surface density. Their mean $\rm{SFR\alpha/M_{HI}}$ and $\rm{SFR\alpha/M_{*}}$ are of the order of $\rm{10^{-9.9}}$ $\rm{yr^{-1}}$. We conclude that while the large scale underdense environment must play some role in galaxy formation and growth through accretion, we find that even with respect to other galaxies in the more mildly underdense regions, the increase in star formation rate is only marginal.

The large-scale structure of the universe is comprised of virialized blob-like clusters, linear filaments, sheet-like walls and huge near empty three-dimensional voids. Characterizing the large scale universe is essential to our understanding of the formation and evolution of galaxies. The density range of clusters, walls and voids are relatively well separated, when compared to filaments, which span a relatively larger range. The large scale filamentary network thus forms an intricate part of the cosmic web. In this paper, we describe Felix, a topology based framework for visual exploration of filaments in the cosmic web. The filamentary structure is represented by the ascending manifold geometry of the 2-saddles in the Morse-Smale complex of the density field. We generate a hierarchy of Morse-Smale complexes and query for filaments based on the density ranges at the end points of the filaments. The query is processed efficiently over the entire hierarchical Morse-Smale complex, allowing for interactive visualization. We apply Felix to computer simulations based on the heuristic Voronoi kinematic model and the standard $\Lambda$CDM cosmology, and demonstrate its usefulness through two case studies. First, we extract cosmic filaments within and across cluster like regions in Voronoi kinematic simulation datasets. We demonstrate that we produce similar results to existing structure finders. Filaments that form the spine of the cosmic web, which exist in high density regions in the current epoch, are isolated using Felix. Also, filaments present in void-like regions are isolated and visualized. These filamentary structures are often over shadowed by higher density range filaments and are not easily characterizable and extractable using other filament extraction methodologies.

We present the structural and star formation properties of 59 void galaxies as part of the Void Galaxy Survey (VGS). Our aim is to study in detail the physical properties of these void galaxies and study the effect of the void environment on galaxy properties. We use Spitzer 3.6 $\rm{\mu m}$ and B-band imaging to study the morphology and color of the VGS galaxies. For their star formation properties, we use Halpha and GALEX near-UV imaging. We compare our results to a range of galaxies of different morphologies in higher density environments. We find that the VGS galaxies are in general disk dominated and star forming galaxies. Their star formation rates are, however, often less than 1 $\rm{M_{\odot}}$ $\rm{yr^{-1}}$. There are two early-type galaxies in our sample as well. In $\rm{r_{e}}$ versus $\rm{M_{B}}$ parameter space, VGS galaxies occupy the same space as dwarf irregulars and spirals.

We investigate the characteristics and the time evolution of the cosmic web from redshift, z=2, to present time, within the framework of the NEXUS+ algorithm. This necessitates the introduction of new analysis tools optimally suited to describe the very intricate and hierarchical pattern that is the cosmic web. In particular, we characterize filaments (walls) in terms of their linear (surface) mass density. This is very good in capturing the evolution of these structures. At early times the cosmos is dominated by tenuous filaments and sheets, which, during subsequent evolution, merge together, such that the present day web is dominated by fewer, but much more massive, structures. We also show that voids are more naturally described in terms of their boundaries and not their centres. We illustrate this for void density profiles, which, when expressed as a function of the distance from void boundary, show a universal profile in good qualitative agreement with the theoretical shell-crossing framework of expanding underdense regions.

We present a method, PhotoWeb, for estimating photometric redshifts of individual galaxies, and their equivalent distance, with megaparsec and even sub-megaparsec accuracy using the Cosmic Web as a constraint over photo-z estimates. PhotoWeb redshift errors for individual galaxies are of the order of delta_z = 0.0007, compared to errors of delta_z = 0.02 for current photo-z techniques. The mean redshift error is of the order of 0.00005-0.0004 compared to mean errors in the range delta_z =z 0.001-0.01 for the best available photo-z estimates in the literature. Current photo-z techniques based on the spectral energy distribution of galaxies and projected clustering produce redshift estimates with large errors due to the poor constraining power the galaxy's spectral energy distribution and projected clustering can provide. The Cosmic Web, on the other hand, provides the strongest constraints on the position of galaxies. The network of walls, filaments and voids occupy ~%10 of the volume of the Universe, yet they contain ~%95 of galaxies. The cosmic web, being a cellular system with well-defined boundaries, sets a restricted set of intermittent positions a galaxy can occupy along a given line-of-sight. Using the information in the density field computed from spectroscopic redshifts we can narrow the possible locations of a given galaxy along the line of sight from a single broad probability distribution (from photo-z) to one or a few narrow peaks. Our first results improve previous photo-z errors by more than one order of magnitude allowing sub-megaparsec errors in some cases. Such accurate estimates for tens of millions of galaxies will allow unprecedented galaxy-LSS studies. In this work we apply our technique to the SDSS photo-z galaxy sample and discuss its performance and future improvements.

The cosmic web is a complex spatial pattern of walls, filaments, cluster nodes and underdense void regions. It emerged through gravitational amplification from the Gaussian primordial density field. Here we infer analytical expressions for the spatial statistics of caustics in the evolving large-scale mass distribution. In our analysis, following the quasi-linear Zeldovich formalism and confined to the 1D and 2D situation, we compute number density and correlation properties of caustics in cosmic density fields that evolve from Gaussian primordial conditions. The analysis can be straightforwardly extended to the 3D situation. We moreover, are currently extending the approach to the non-linear regime of structure formation by including higher order Lagrangian approximations and Lagrangian effective field theory.

The current Lambda CDM cosmological model predicts that galaxy evolution proceeds more slowly in lower density environments, suggesting that voids are a prime location to search for relatively pristine galaxies that are representative of the building blocks of early massive galaxies. To test the assumption that void galaxies are more pristine, we compare the evolutionary properties of a sample of dwarf galaxies selected specifically to lie in voids with a sample of similar isolated dwarf galaxies in average density environments. We measure gas-phase oxygen abundances and gas fractions for eight dwarf galaxies (M_r > -16.2), carefully selected to reside within the lowest density environments of seven voids, and apply the same calibrations to existing samples of isolated dwarf galaxies. We find no significant difference between these void dwarf galaxies and the isolated dwarf galaxies, suggesting that dwarf galaxy chemical evolution proceeds independent of the large-scale environment. While this sample is too small to draw strong conclusions, it suggests that external gas accretion is playing a limited role in the chemical evolution of these systems, and that this evolution is instead dominated mainly by the internal secular processes that are linking the simultaneous growth and enrichment of these galaxies.

We study the formation and evolution of the cosmic web, using the high-resolution CosmoGrid \$\Lambda\$CDM simulation. In particular, we investigate the evolution of the large-scale structure around void halo groups, and compare this to observations of the VGS-31 galaxy group, which consists of three interacting galaxies inside a large void. The structure around such haloes shows a great deal of tenuous structure, with most of such systems being embedded in intra-void filaments and walls. We use the Nexus+ algorithm to detect walls and filaments in CosmoGrid, and find them to be present and detectable at every scale. The void regions embed tenuous walls, which in turn embed tenuous filaments. We hypothesize that the void galaxy group of VGS-31 formed in such an environment.

Voids represent a unique environment for the study of galaxy evolution, as the lower density environment is expected to result in shorter merger histories and slower evolution of galaxies. This provides an ideal opportunity to test theories of galaxy formation and evolution. Imaging of the neutral hydrogen, central in both driving and regulating star formation, directly traces the gas reservoir and can reveal interactions and signs of cold gas accretion. For a new Void Galaxy Survey (VGS), we have carefully selected a sample of 59 galaxies that reside in the deepest underdensities of geometrically identified voids within the SDSS at distances of ~100 Mpc, and pursued deep UV, optical, Halpha, IR, and HI imaging to study in detail the morphology and kinematics of both the stellar and gaseous components. This sample allows us to not only examine the global statistical properties of void galaxies, but also to explore the details of the dynamical properties. We present an overview of the VGS, and highlight key results on the HI content and individually interesting systems. In general, we find that the void galaxies are gas rich, low luminosity, blue disk galaxies, with optical and HI properties that are not unusual for their luminosity and morphology. We see evidence of both ongoing assembly, through the gas dynamics between interacting systems, and significant gas accretion, seen in extended gas disks and kinematic misalignments. The VGS establishes a local reference sample to be used in future HI surveys (CHILES, DINGO, LADUMA) that will directly observe the HI evolution of void galaxies over cosmic time.

We study the substructure population of Milky Way (MW)-mass halos in the $\Lambda$CDM cosmology using a novel procedure to extrapolate subhalo number statistics beyond the resolution limit of N-body simulations. The technique recovers the mean and the variance of the subhalo abundance, but not its spatial distribution. It extends the dynamic range over which precise statistical predictions can be made by the equivalent of performing a simulation with 50 times higher resolution, at no additional computational cost. We apply this technique to MW-mass halos, but it can easily be applied to halos of any mass. We find up to $20\%$ more substructures in MW-mass halos than found in previous studies. Our analysis lowers the mass of the MW halo required to accommodate the observation that the MW has only three satellites with a maximum circular velocity $V_{max}\ge30 km/s$ in the $\Lambda$CDM cosmology. The probability of having a subhalo population similar to that in the MW is $20\%$ for a virial mass, $M_{200}=1\times10^{12} M_\odot$ and practically zero for halos more massive than $M_{200}=2\times10^{12} M_\odot$.

We use the distribution of maximum circular velocities, $V_{max}$, of satellites in the Milky Way (MW) to constrain the virial mass, $M_{200}$, of the Galactic halo under an assumed prior of a $\Lambda$CDM universe. This is done by analysing the subhalo populations of a large sample of halos found in the Millennium II cosmological simulation. The observation that the MW has at most three subhalos with $V_{max}\ge30 km/s$ requires a halo mass $M_{200}\le1.4\times10^{12} M_\odot$, while the existence of the Magellanic Clouds (assumed to have $V_{max}\ge60 km/s$) requires $M_{200}\ge1.0\times10^{12} M_\odot$. The first of these conditions is necessary to avoid the "too-big-to-fail" problem highlighted by Boylan-Kolchin et al., while the second stems from the observation that massive satellites like the Magellanic Clouds are rare. When combining both requirements, we find that the MW halo mass must lie in the range $0.25 \le M_{200}/(10^{12} M_\odot) \le 1.4$ at $90\%$ confidence. The gap in the abundance of Galactic satellites between $30 km/s\le V_{max} \le 60 km/s$ places our galaxy in the tail of the expected satellite distribution.

The cosmic web is the largest scale manifestation of the anisotropic gravitational collapse of matter. It represents the transitional stage between linear and non-linear structures and contains easily accessible information about the early phases of structure formation processes. Here we investigate the characteristics and the time evolution of morphological components since. Our analysis involves the application of the NEXUS Multiscale Morphology Filter (MMF) technique, predominantly its NEXUS+ version, to high resolution and large volume cosmological simulations. We quantify the cosmic web components in terms of their mass and volume content, their density distribution and halo populations. We employ new analysis techniques to determine the spatial extent of filaments and sheets, like their total length and local width. This analysis identifies cluster and filaments as the most prominent components of the web. In contrast, while voids and sheets take most of the volume, they correspond to underdense environments and are devoid of group-sized and more massive haloes. At early times the cosmos is dominated by tenuous filaments and sheets, which, during subsequent evolution, merge together, such that the present day web is dominated by fewer, but much more massive, structures. The analysis of the mass transport between environments clearly shows how matter flows from voids into walls, and then via filaments into cluster regions, which form the nodes of the cosmic web. We also study the properties of individual filamentary branches, to find long, almost straight, filaments extending to distances larger than 100Mpc/h. These constitute the bridges between massive clusters, which seem to form along approximatively straight lines.

We describe how the dynamics of cosmic structure formation defines the intricate geometric structure of the spine of the cosmic web. The Zeldovich approximation is used to model the backbone of the cosmic web in terms of its singularity structure. The description by Arnold et al. (1982) in terms of catastrophe theory forms the basis of our analysis. This two-dimensional analysis involves a profound assessment of the Lagrangian and Eulerian projections of the gravitationally evolving four-dimensional phase-space manifold. It involves the identification of the complete family of singularity classes, and the corresponding caustics that we see emerging as structure in Eulerian space evolves. In particular, as it is instrumental in outlining the spatial network of the cosmic web, we investigate the nature of spatial connections between these singularities. The major finding of our study is that all singularities are located on a set of lines in Lagrangian space. All dynamical processes related to the caustics are concentrated near these lines. We demonstrate and discuss extensively how all 2D singularities are to be found on these lines. When mapping this spatial pattern of lines to Eulerian space, we find a growing connectedness between initially disjoint lines, resulting in a percolating network. In other words, the lines form the blueprint for the global geometric evolution of the cosmic web.

We study the formation and evolution of filamentary configurations of dark matter haloes in voids. Our investigation uses the high-resolution LambdaCDM simulation CosmoGrid to look for void systems resembling the VGS_31 elongated system of three interacting galaxies that was recently discovered by the Void Galaxy Survey (VGS) inside a large void in the SDSS galaxy redshift survey. HI data revealed these galaxies to be embedded in a common elongated envelope, possibly embedded in intravoid filament. In the CosmoGrid simulation we look for systems similar to VGS_31 in mass, size and environment. We find a total of eight such systems. For these systems, we study the distribution of neighbour haloes, the assembly and evolution of the main haloes and the dynamical evolution of the haloes, as well as the evolution of the large-scale structure in which the systems are embedded. The spatial distribution of the haloes follows that of the dark matter environment. We find that VGS_31-like systems have a large variation in formation time, having formed between 10 Gyr ago and the present epoch. However, the environments in which the systems are embedded evolved resemble each other substantially. Each of the VGS_31-like systems is embedded in an intra-void wall, that no later than z = 0.5 became the only prominent feature in its environment. While part of the void walls retain a rather featureless character, we find that around half of them are marked by a pronounced and rapidly evolving substructure. Five haloes find themselves in a tenuous filament of a few Mpc/h long inside the intra-void wall. Finally, we compare the results to observed data from VGS_31. Our study implies that the VGS_31 galaxies formed in the same (proto)filament, and did not meet just recently. The diversity amongst the simulated halo systems indicates that VGS_31 may not be typical for groups of galaxies in voids.

We present the relation between the genus in cosmology and the Betti numbers for excursion sets of three- and two-dimensional smooth Gaussian random fields, and numerically investigate the Betti numbers as a function of threshold level. Betti numbers are topological invariants of figures that can be used to distinguish topological spaces. In the case of the excursion sets of a three-dimensional field there are three possibly non-zero Betti numbers; $\beta_0$ is the number of connected regions, $\beta_1$ is the number of circular holes, and $\beta_2$ is the number of three-dimensional voids. Their sum with alternating signs is the genus of the surface of excursion regions. It is found that each Betti number has a dominant contribution to the genus in a specific threshold range. $\beta_0$ dominates the high-threshold part of the genus curve measuring the abundance of high density regions (clusters). $\beta_1$ dominates the genus near the median thresholds which measures the topology of negatively curved iso-density surfaces, and $\beta_2$ corresponds to the low-threshold part measuring the void abundance. We average the Betti number curves (the Betti numbers as a function of the threshold level) over many realizations of Gaussian fields and find that both the amplitude and shape of the Betti number curves depend on the slope of the power spectrum $n$ in such a way that their shape becomes broader and their amplitude drops less steeply than the genus as $n$ decreases. This behaviour contrasts with the fact that the shape of the genus curve is fixed for all Gaussian fields regardless of the power spectrum. Even though the Gaussian Betti number curves should be calculated for each given power spectrum, we propose to use the Betti numbers for better specification of the topology of large scale structures in the universe.

We study the topology of the Megaparsec Cosmic Web in terms of the scale-dependent Betti numbers, which formalize the topological information content of the cosmic mass distribution. While the Betti numbers do not fully quantify topology, they extend the information beyond conventional cosmological studies of topology in terms of genus and Euler characteristic. The richer information content of Betti numbers goes along the availability of fast algorithms to compute them. For continuous density fields, we determine the scale-dependence of Betti numbers by invoking the cosmologically familiar filtration of sublevel or superlevel sets defined by density thresholds. For the discrete galaxy distribution, however, the analysis is based on the alpha shapes of the particles. These simplicial complexes constitute an ordered sequence of nested subsets of the Delaunay tessellation, a filtration defined by the scale parameter, $\alpha$. As they are homotopy equivalent to the sublevel sets of the distance field, they are an excellent tool for assessing the topological structure of a discrete point distribution. In order to develop an intuitive understanding for the behavior of Betti numbers as a function of $\alpha$, and their relation to the morphological patterns in the Cosmic Web, we first study them within the context of simple heuristic Voronoi clustering models. Subsequently, we address the topology of structures emerging in the standard LCDM scenario and in cosmological scenarios with alternative dark energy content. The evolution and scale-dependence of the Betti numbers is shown to reflect the hierarchical evolution of the Cosmic Web and yields a promising measure of cosmological parameters. We also discuss the expected Betti numbers as a function of the density threshold for superlevel sets of a Gaussian random field.

Cosmological voids provide a unique environment for the study of galaxy formation and evolution. The galaxy population in their interior have significantly different properties than average field galaxies. As part of our Void Galaxy Survey (VGS), we have found a system of three interacting galaxies (VGS_31) inside a large void. VGS_31 is a small elongated group whose members are embedded in a common HI envelope. The HI picture suggests a filamentary structure with accretion of intergalactic cold gas from the filament onto the galaxies. We present deep optical and narrow band H_alpha data, optical spectroscopy, near-UV and far-UV GALEX and CO(1-0) data. We find that one of the galaxies, a Markarian object, has a ring-like structure and a tail evident both in optical and HI. While all three galaxies form stars in their central parts, the tail and the ring of the Markarian object are devoid of star formation. We discuss these findings in terms of a gravitational interaction and ongoing growth of galaxies out of a filament. VGS_31 is one of the first observed examples of a filamentary structure in a void. It is an important prototype for understanding the formation of substructure in a void. This system also shows that the galaxy evolution in voids can be as dynamic as in high density environments.

In this contribution we present the preliminary results regarding the non-linear BAO signal in higher-order statistics of the cosmic density field. We use ensembles of N-body simulations to show that the non-linear evolution changes the amplitudes of the BAO signal, but has a negligible effect on the scale of the BAO feature. The latter observation accompanied by the fact that the BAO feature amplitude roughly doubles as one moves to higher orders, suggests that the higher-order correlation amplitudes can be used as probe of the BAO signal.

In this study we show that the skewness S_3 of the cosmic density field contains a significant and potentially detectable and clean imprint of Baryonic Acoustic Oscillations. Although the BAO signal in the skewness has a lower amplitude than second order measures like the two-point correlation function and power spectrum, it has the advantage of a considerably lower sensitivity to systematic influences. Because it lacks a direct dependence on bias if this concerns simple linear bias, skewness will be considerably less beset by uncertainties due to galaxy bias. Also, it has a weaker sensitivity to redshift distortion effects. We use perturbation theory to evaluate the magnitude of the effect on the volume-average skewness, for various cosmological models. One important finding of our analysis is that the skewness BAO signal occurs at smaller scales than that in second order statistics. For an LCDM spectrum with WMAP7 normalization, the BAO feature has a maximum wiggle amplitude of ~3% and appears at a scale of ~82Mpc/h. We conclude that the detection of BAO wiggles in future extensive galaxy surveys via the skewness of the observed galaxy distribution may provide us with a useful, and potentially advantageous, measure of the nature of Dark Energy.

We present a new way to formulate the geometry of the Cosmic Web in terms of Lagrangian space. The Adhesion model has an ingenious geometric interpretation out of which the spine of the Cosmic Web emerges naturally. Within this context we demonstrate a deep connection of the relation between Eulerian and Lagrangian space with that between Voronoi and Delaunay tessellations.

One of the important unknowns of current cosmology concerns the effects of the large scale distribution of matter on the formation and evolution of dark matter haloes and galaxies. One main difficulty in answering this question lies in the absence of a robust and natural way of identifying the large scale environments and their characteristics. This work summarizes the NEXUS+ formalism which extends and improves our multiscale scale-space MMF method. The new algorithm is very successful in tracing the Cosmic Web components, mainly due to its novel filtering of the density in logarithmic space. The method, due to its multiscale and hierarchical character, has the advantage of detecting all the cosmic structures, either prominent or tenuous, without preference for a certain size or shape. The resulting filamentary and wall networks can easily be characterized by their direction, thickness, mass density and density profile. These additional environmental properties allows to us to investigate not only the effect of environment on haloes, but also how it correlates with the environment characteristics.

We introduce the NEXUS algorithm for the identification of Cosmic Web environments: clusters, filaments, walls and voids. This is a multiscale and automatic morphological analysis tool that identifies all the cosmic structures in a scale free way, without preference for a certain size or shape. We develop the NEXUS method to incorporate the density, tidal field, velocity divergence and velocity shear as tracers of the Cosmic Web. We also present the NEXUS+ procedure which, taking advantage of a novel filtering of the density in logarithmic space, is very successful at identifying the filament and wall environments in a robust and natural way. To asses the algorithms we apply them to an N-body simulation. We find that all methods correctly identify the most prominent filaments and walls, while there are differences in the detection of the more tenuous structures. In general, the structures traced by the density and tidal fields are clumpier and more rugged than those present in the velocity divergence and velocity shear fields. We find that the NEXUS+ method captures much better the filamentary and wall networks and is successful in detecting even the fainter structures. We also confirm the efficiency of our methods by examining the dark matter particle and halo distributions.

We showed how the shape of cosmic voids can be used to distinguish between different models of dark energy using galaxy positions.

Aims: We assess the sensitivity of void shapes to the nature of dark energy that was pointed out in recent studies. We investigate whether or not void shapes are useable as an observational probe in galaxy redshift surveys. We focus on the evolution of the mean void ellipticity and its underlying physical cause. Methods: We analyse the morphological properties of voids in five sets of cosmological N-body simulations, each with a different nature of dark energy. Comparing voids in the dark matter distribution to those in the halo population, we address the question of whether galaxy redshift surveys yield sufficiently accurate void morphologies. Voids are identified using the parameter free Watershed Void Finder. The effect of redshift distortions is investigated as well. Results: We confirm the statistically significant sensitivity of voids in the dark matter distribution. We identify the level of clustering as measured by \sigma_8(z) as the main cause of differences in mean void shape <\epsilon>. We find that in the halo and/or galaxy distribution it is practically unfeasible to distinguish at a statistically significant level between the various cosmologies due to the sparsity and spatial bias of the sample.

We introduce the numbers of hot and cold spots, $n_h$ and $n_c$, of excursion sets of the CMB temperature anisotropy maps as statistical observables that can discriminate different non-Gaussian models. We numerically compute them from simulations of non-Gaussian CMB temperature fluctuation maps. The first kind of non-Gaussian model we study is the local type primordial non-Gaussianity. The second kind of models have some specific form of the probability distribution function from which the temperature fluctuation value at each pixel is drawn, obtained using HEALPIX. We find the characteristic non-Gaussian deviation shapes of $n_h$ and $n_c$, which is distinct for each of the models under consideration. We further demonstrate that $n_h$ and $n_c$ carry additional information compared to the genus, which is just their linear combination, making them valuable additions to the Minkowski Functionals in constraining non-Gaussianity.

In this video we highlight the application of Computational Geometry to our understanding of the formation and dynamics of the Cosmic Web. The emergence of this intricate and pervasive weblike structure of the Universe on Megaparsec scales can be approximated by a well-known equation from fluid mechanics, the Burgers' equation. The solution to this equation can be obtained from a geometrical formalism. We have extended and improved this method by invoking weighted Delaunay and Voronoi tessellations. The duality between these tessellations finds a remarkable and profound reflection in the description of physical systems in Eulerian and Lagrangian terms. The resulting Adhesion formalism provides deep insight into the dynamics and topology of the Cosmic Web. It uncovers a direct connection between the conditions in the very early Universe and the complex spatial patterns that emerged out of these under the influence of gravity.

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 introduce a new descriptor of the weblike pattern in the distribution of galaxies and matter: the scale dependent Betti numbers which formalize the topological information content of the cosmic mass distribution. While the Betti numbers do not fully quantify topology, they extend the information beyond conventional cosmological studies of topology in terms of genus and Euler characteristic used in earlier analyses of cosmological models. The richer information content of Betti numbers goes along with the availability of fast algorithms to compute them. When measured as a function of scale they provide a "Betti signature" for a point distribution that is a sensitive yet robust discriminator of structure. The signature is highly effective in revealing differences in structure arising in different cosmological models, and is exploited towards distinguishing between different dark energy models and may likewise be used to trace primordial non-Gaussianities. In this study we demonstrate the potential of Betti numbers by studying their behaviour in simulations of cosmologies differing in the nature of their dark energy.

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.

We have found that KK 246, the only confirmed galaxy located within the nearby Tully Void, is a dwarf galaxy with an extremely extended H I disk and signs of an H I cloud with anomalous velocity. It also exhibits clear misalignment between the kinematical major and minor axes, indicative of an oval distortion, and a general misalignment between the H I and optical major axes. We measure a H I mass of 1.05 +- 0.08 x 10^8 M_sun, and a H I extent 5 times that of the stellar disk, one of the most extended H I disks known. We estimate a dynamical mass of 4.1 x 10^9 M_sun, making this also one of the darkest galaxies known, with a mass-to-light ratio of 89. The relative isolation and extreme underdense environment make this an interesting case for examining the role of gas accretion in galaxy evolution.

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.

In this work we investigate higher order statistics for the $\lcdm$ and ReBEL scalar-interacting dark matter models by analyzing $180\hmpc$ dark matter N-body simulation ensembles. The N-point correlation functions and the related hierarchical amplitudes, such as skewness and kurtosis, are computed using the Count-In-Cells method. Our studies demonstrate that the hierarchical amplitudes $S_n$ of the scalar-interacting dark matter model significantly deviate from the values in the $\lcdm$ cosmology on scales comparable and smaller then the screening length $r_s$ of a given scalar-interacting model. The corresponding additional forces that enhance the total attractive force exerted on dark matter particles at galaxy scales lowers the values of the hierarchical amplitudes $S_n$. We conclude that hypothetical additional exotic interactions in the dark matter sector should leave detectable markers in the higher-order correlation statistics of the density field. We focussed in detail on the redshift evolution of the dark matter field's skewness and kurtosis. From this investigation we find that the deviations from the canonical $\lcdm$ model introduced by the presence of the ``fifth'' force attain a maximum value at redshifts $0.5<z<2$. We therefore conclude that moderate redshift data are better suited for setting observational constraints on the investigated ReBEL models.

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 analyze the structure and connectivity of the distinct morphologies that define the Cosmic Web. With the help of our Multiscale Morphology Filter (MMF), we dissect the matter distribution of a cosmological $\Lambda$CDM N-body computer simulation into cluster, filaments and walls. The MMF is ideally suited to adress both the anisotropic morphological character of filaments and sheets, as well as the multiscale nature of the hierarchically evolved cosmic matter distribution. The results of our study may be summarized as follows: i).- While all morphologies occupy a roughly well defined range in density, this alone is not sufficient to differentiate between them given their overlap. Environment defined only in terms of density fails to incorporate the intrinsic dynamics of each morphology. This plays an important role in both linear and non linear interactions between haloes. ii).- Most of the mass in the Universe is concentrated in filaments, narrowly followed by clusters. In terms of volume, clusters only represent a minute fraction, and filaments not more than 9%. Walls are relatively inconspicous in terms of mass and volume. iii).- On average, massive clusters are connected to more filaments than low mass clusters. Clusters with $M \sim 10^{14}$ M$_{\odot}$ h$^{-1}$ have on average two connecting filaments, while clusters with $M \geq 10^{15}$ M$_{\odot}$ h$^{-1}$ have on average five connecting filaments. iv).- Density profiles indicate that the typical width of filaments is 2$\Mpch$. Walls have less well defined boundaries with widths between 5-8 Mpc h$^{-1}$. In their interior, filaments have a power-law density profile with slope ${\gamma}\approx -1$, corresponding to an isothermal density profile.

We search for and find fossil evidence that the distribution of the spin axes of galaxies in cosmic web filaments relative to their host filaments are not randomly distributed. This would indicate that the action of large scale tidal torques effected the alignments of galaxies located in cosmic filaments. To this end, we constructed a catalogue of clean filaments containing edge-on galaxies. We started by applying the Multiscale Morphology Filter (MMF) technique to the galaxies in a redshift-distortion corrected version of the Sloan Digital Sky Survey DR5. From that sample we extracted those 426 filaments that contained edge-on galaxies (b/a < 0.2). These filaments were then visually classified relative to a variety of quality criteria. Statistical analysis using "feature measures" indicates that the distribution of orientations of these edge-on galaxies relative to their parent filament deviate significantly from what would be expected on the basis of a random distribution of orientations. The interpretation of this result may not be immediately apparent, but it is easy to identify a population of 14 objects whose spin axes are aligned perpendicular to the spine of the parent filament (\cos \theta < 0.2). The candidate objects are found in relatively less dense filaments. This might be expected since galaxies in such locations suffer less interaction with surrounding galaxies, and consequently better preserve their tidally induced orientations relative to the parent filament. The technique of searching for fossil evidence of alignment yields relatively few candidate objects, but it does not suffer from the dilution effects inherent in correlation analysis of large samples.

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.

The evolution of marginally bound supercluster-like objects in an accelerating LambdaCDM Universe is followed, by means of cosmological simulations, from the present time to an expansion factor a = 100. The objects are identified on the basis of the binding density criterion introduced by Dunner et al. (2006). superclusters are identified with the ones whose mass M > 10^15 M_sun/h, the most massive one with M ~ 8x10^15 M_sun/h, comparable to the Shapley supercluster. The spatial distribution of the superclusters remains essentially the same after the present epoch, reflecting the halting growth of the Cosmic Web as Lambda gets to dominate the expansion of the Universe. The same trend can be seen in the stagnation of the development of the mass function of virialized haloes and bound objects. The situation is considerably different when looking at the internal evolution, quantified in terms of their shape, compactness and density profile, and substructure in terms of their multiplicity function. We find a continuing evolution from a wide range of triaxial shapes at a = 1 to almost perfect spherical shapes at a = 100. We also find a systematic trend towards a higher concentration. Meanwhile, we see their substructure gradually disappearing, as the surrounding subclumps fall in and merge to form one coherent, virialized system.

We explore the effects of dark matter and dark energy on the dynamical scaling properties of galaxy clusters. We investigate the cluster Faber-Jackson (FJ), Kormendy and Fundamental Plane (FP) relations between the mass, radius and velocity dispersion of cluster size halos in cosmological $N$-body simulations. The simulations span a wide range of cosmological parameters, representing open, flat and closed Universes. Independently of the cosmology, we find that the simulated clusters are close to a perfect virial state and do indeed define a Fundamental Plane. The fitted parameters of the FJ, Kormendy and FP relationships do not show any significant dependence on $\Omega_m$ and/or $\Omega_{\Lambda}$. The one outstanding effect is the influence of $\Omega_{m}$ on the thickness of the Fundamental Plane. Following the time evolution of our models, we find slight changes of FJ and Kormendy parameters in high $\Omega_m$ universe, along with a slight decrease of FP fitting parameters. We also see an initial increase of the FP thickness followed by a convergence to a nearly constant value. The epoch of convergence is later for higher values of $\Omega_m$ while the thickness remains constant in the low $\Omega_m$ $\Lambda$-models. We also find a continuous increase of the FP thickness in the Standard CDM (SCDM) cosmology. There is no evidence that these differences are due to the different power spectrum slope at cluster scales. From the point of view of the FP, there is little difference between clusters that quietly accreted their mass and those that underwent massive mergers. The principal effect of strong mergers is to change significantly the ratio of the half-mass radius $r_{half}$ to the harmonic mean radius $r_h$.

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.

Nikos Voglis had many astronomical interests, among them was the question of the origin of galactic angular momentum. In this short tribute we review how this subject has changed since the 1970's and how it has now become evident that gravitational tidal forces have not only caused galaxies to rotate, but have also acted to shape the very cosmic structure in which those galaxies are found. We present recent evidence for this based on data analysis techniques that provide objective catalogues of clusters, filaments and voids.

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.

The lecture notes describe the Delaunay Tessellation Field Estimator for Cosmic Web analysis. The high sensitivity of Voronoi/Delaunay tessellations to the local point distribution is used to obtain estimates of density and related quantities. The adaptive and minimum triangulation properties of Delaunay tessellations are utilized for interpolation to continuous volume-covering density and velocity fields. The adaptivity to the local density and geometry preserves three key characteristics of the cosmic matter distribution: 1) its hierarchical nature,2) its anisotropic weblike morphology, 3) the presence and structure of voids. The lecture notes discuss the mathematical context and background of the DTFE method, in essence a first order version of Natural Neighbour methods, smooth and higher order spatial interpolation techniques. This is followed by an outline of the DTFE procedure. A series of tests on spatially complex point distributions provides quantitative information on its performance and noise characteristics. We conclude with a few cosmological applications. The reconstruction of the spatial weblike patterns in N-body simulations and in the 2dFGRS provide a visual impression of DTFE's performance. The analysis of the velocity divergence pdf and maps obtained of the velocity flow in the Local Universe (implied by PSCz) illustrate its promise for probing the dynamics of cosmic structure. The Watershed Void Finder (Platen et al. 2007), for the study of voids, and Multiscale Morphology Filter (Aragon-Calvo et al. 2007), for identifying and analyzing weblike networks of filaments, sheets and clusters, are based upon DTFE density fields. We also mention SimpleX, a related random lattice radiative transfer formalism (Ritzerveld & Icke 2006).

The spatial cosmic matter distribution on scales of a few up to more than a hundred Megaparsec displays a salient and pervasive foamlike pattern. Voronoi tessellations are a versatile and flexible mathematical model for such weblike spatial patterns. They would be the natural asymptotic result of an evolution in which low-density expanding void regions dictate the spatial organization of the Megaparsec Universe, while matter assembles in high-density filamentary and wall-like interstices between the voids. We describe the results of ongoing investigations of a variety of aspects of cosmologically relevant spatial distributions and statistics within the framework of Voronoi tessellations. Particularly enticing is the finding of a profound scaling of both clustering strength and clustering extent for the distribution of tessellation nodes, suggestive for the clustering properties of galaxy clusters. Cellular patterns may be the source of an intrinsic ``geometrically biased'' clustering.

We apply the Delaunay Tessellation Field Estimator (DTFE) to reconstruct and analyze the matter distribution and cosmic velocity flows in the Local Universe on the basis of the PSCz galaxy survey. The prime objective of this study is the production of optimal resolution three-dimensional maps fully volume-covering of the volume-weighted velocity and density fields throughout the nearby Universe, out to a distance of 150 Mpc/h. Based on the Voronoi and Delaunay tessellation defined by the spatial galaxy sample, DTFE involves the estimate of density values on the basis of the volume of the related Delaunay tetrahedra and the subsequent use of theDelaunay tessellation as natural multidimensional (linear) interpolation grid for the corresponding density and velocity fields throughout the sample volume. The linearized model of the spatial galaxy distribution and the corresponding peculiar velocities of the PSCz galaxy sample, produced by Brachini et al. (1999), forms the input sample for the DTFE study. The DTFE maps reproduce the high-density supercluster regions in optimal detail, both their internal structure as well as their elongated or flattened shape. The corresponding velocity flows trace the bulk and shear flows marking the region extending from the Pisces-Perseus supercluster, via the Local superclusters, towards the Hydra-Centaurus and the Shapley concentration. The most outstanding and unique feature of the DTFE maps is the sharply defined radial outflow regions in and around underdense voids, marking the dynamical importance of voids in the Local Universe. The maximum expansion rate of voids defines a sharp cutoff in the DTFE velocity divergence pdf. We found that on the basis of this cutoff DTFE manages to consistently reproduce the value of O_m ~ 0.35 underlying the linearized velocity dataset.

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 here a new method, MMF, for automatically segmenting cosmic structure into its basic components: clusters, filaments, and walls. Importantly, the segmentation is scale independent, so all structures are identified without prejudice as to their size or shape. The method is ideally suited for extracting catalogues of clusters, walls, and filaments from samples of galaxies in redshift surveys or from particles in cosmological N-body simulations: it makes no prior assumptions about the scale or shape of the structures.}

The MMF technique is used to segment the cosmic web as seen in a cosmological N-body simulation into wall-like and filament-like structures. We find that the spins and shapes of dark matter haloes are significantly correlated with each other and with the orientation of their host structures. The shape orientation is such that the halo minor axes tend to lie perpendicular to the host structure, be it a wall or filament. The orientation of the halo spin vector is mass dependent. Low mass haloes in walls and filaments have a tendency to have their spins oriented within the parent structure, while higher mass haloes in filaments have spins that tend to lie perpendicular to the parent structure.

We discuss the intimate relationship between the filamentary features and the rare dense compact cluster nodes in this network, via the large scale tidal field going along with them, following the cosmic web theory developed Bond et al. The Megaparsec scale tidal shear pattern is responsible for the contraction of matter into filaments, and its link with the cluster locations can be understood through the implied quadrupolar mass distribution in which the clusters are to be found at the sites of the overdense patches. We present a new technique for tracing the cosmic web, identifying planar walls, elongated filaments and cluster nodes in the galaxy distribution. This will allow the practical exploitation of the concept of the cosmic web towards identifying and tracing the locations of the gaseous WHIM. These methods, the Delaunay Tessellation Field Estimator (DTFE) and the Morphology Multiscale Filter (MMF) find their basis in computational geometry and visualization.

We study the clustering properties of the recently compiled HIPASS neutral hydrogen (HI) sources catalogue using the two point correlation function in redshift space. If the two point correlation is modelled as a power law, $\xi(r)=(r_{0}/r)^{\gamma}$, then the best-fitting parameters for the HI selected galaxies are $r_{0}=4.3 \pm 0.3 h^{-1}$ Mpc with $\gamma=1.42 \pm 0.14$. Fixing the slope to its universal value $\gamma=1.8$, we obtain $r_{0}= 4.1^{+0.2}_{-0.3} h^{-1}$ Mpc. Comparing the measured two point correlation function with the predictions of the concordance cosmological model, we find that at the present epoch the HI selected galaxies are anti-biased with respect to the underline matter fluctuation field with bias value being $b_{0}\simeq 0.63$. Also, we investigate the evolution of the linear bias factor, $b(z)$, and we find to be a strongly dependent function of redshift. Finally, we estimate the HI selected galaxies redshift space distortion parameter $\beta \simeq \Omega_{\rm m}^{0.6}/b_{0}$ in order to correct the correlation function for the peculiar motions and we find $\beta\simeq 0.77$. Taking into account the corrections for the redshift space distortions the HI correlation length in real space is $r^{\rm re}_{0}=3^{+0.4}_{-0.3}h^{-1}$Mpc for $\gamma=1.42$ while using $\gamma=1.8$ we get $r^{\rm re}_{0}= 3.1 \pm 0.3 h^{-1}$Mpc.

We assess the extent to which the PSCz survey encapsulates the complete or major share of matter inhomogeneities responsible for the external tidal forces affecting the peculiar velocity flow within the Local Supercluster [LS]. We investigate this issue on the basis of mock galaxy catalogs drawn from two large Nbody simulations of cosmic structure formation in two different cosmological scenarios. From these a variety of datasets is selected imitating the observational conditions of either the volume-limited LS mimicking NBG catalog or the deeper magnitude-limited PSCz catalog. The mildly nonlinear dynamics in the mock NBG and PSCz velocities are analyzed by means of an optimized implementation of the Least Action Principle, the Fast Action Method [FAM]. By comparing the modeled and true velocities of the corresponding mock catalogs we assess the extent and nature of the external tidal influence on the LS volume. We find that the dynamics in the inner 30Mpc/h volume is strongly affected by the external forces. Most of the external forces can be traced back to a depth of no more than 100Mpc/h. This is concluded from the fact that the FAM reconstructions of the 100Mpc/h PSCz volume appear to have included most gravitational influences. In addition, we demonstrate that for all considered cosmological models the bulk flow and shear components of the tidal velocity field generated by the external distribution of PSCz galaxies provides sufficient information for representing the full external tidal force field.

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.

We present a model for the distribution of void sizes and its evolution in the context of hierarchical scenarios of gravitational structure formation. We find that at any cosmic epoch the voids have a size distribution which is well-peaked about a characteristic void size which evolves self-similarly in time. This is in distinct contrast to the distribution of virialized halo masses which does not have a small-scale cut-off. In our model, the fate of voids is ruled by two processes. The first process affects those voids which are embedded in larger underdense regions: the evolution is effectively one in which a larger void is made up by the mergers of smaller voids, and is analogous to how massive clusters form from the mergers of less massive progenitors. The second process is unique to voids, and occurs to voids which happen to be embedded within a larger scale overdensity: these voids get squeezed out of existence as the overdensity collapses around them. It is this second process which produces the cut-off at small scales. In the excursion set formulation of cluster abundance and evolution, solution of the cloud-in-cloud problem, i.e., counting as clusters only those objects which are not embedded in larger clusters, requires study of random walks crossing one barrier. We show that a similar formulation of void evolution requires study of a two-barrier problem: one barrier is required to account for voids-in-voids, and the other for voids-in-clouds. Thus, in our model, the void size distribution is a function of two parameters, one of which reflects the dynamics of void formation, and the other the formation of collapsed objects.

Within the context of hierarchical scenarios of gravitational structure formation we describe how an evolving hierarchy of voids evolves on the basis of {\em two} processes, the {\em void-in-void} process and the {\it void-in-cloud} process. The related analytical formulation in terms of a {\it two-barrier} excursion problem leads to a self-similarly evolving peaked void size distribution.

We present the results of a study confronting density maps reconstructed by the Delaunay Tessellation Field Estimator (DTFE) and by regular SPH kernel-based techniques. The comparison between the two methods clearly demonstrates the superior performance of the DTFE with respect to conventional SPH methods, in particular at locations where SPH appears to fail. The DTFE is a fully self-adaptive technique for reconstructing continuous density fields from discrete particle distributions, and is based upon the corresponding Delaunay tessellation. Its principal asset is its complete independence of arbitrary smoothing functions and parameters specifying the properties of these. As a result it manages to faithfully reproduce the anisotropies of the local particle distribution and through its adaptive and local nature proves to be optimally suited for uncovering the full structural richness in the density distribution. Through the improvement in local density estimates, calculations invoking the DTFE will yield a much better representation of physical processes which depend on density. The presented results form an encouraging step towards the application and insertion of the DTFE in astrophysical hydrocodes. We describe an outline for the construction of a particle hydrodynamics code in which the DTFE replaces kernel-based methods. Further discussion addresses the issue and possibilities for a moving grid based hydrocode invoking the DTFE, and Delaunay tessellations, in an attempt to combine the virtues of the Eulerian and Lagrangian approaches.

A review on the properties and dynamical origin and nature of the cosmic foam, the tenuous space-filling frothy network permeating the interior of the Universe. We discuss the properties of this striking and intriguing pattern, describing its observational appearance, and seeking to elucidate its dynamical origin and nature. An extensive discussion on the gravitational formation and dynamical evolution of weblike patterns puts particular emphasis on the formative role of the generic anisotropy of the cosmic gravitational force fields. These tidal fields play an essential role in shaping the pattern of the large scale cosmic matter distribution. Special attention is put on a geometrical and stochastic framework for a systematic evaluation of its fossil information content on the cosmic structure formation process. Its distinct geometric character and the stochastic nature provides the cosmic web with some unique and at first unexpected properties. The implications for galaxy clustering are discussed on the basis of its relevant branch of mathematics, stochastic geometry. Central within this context are Voronoi tessellations, which have been found to represent a surprisingly versatile model for spatial cellular distributions.

The presentations in the 2nd Hellenic Cosmology Workshop covered a broad range of cosmological subjects. Following the analogy with the architectural masters of Ancient Greece, the contributions are presented in ``architectural'' sessions. The session on the ``Architects'' involves ``fundamental cosmology'', focussing on the topology and geometry of our Universe and the fundamental physical laws and processes at work in the first decisive moments of the (Very) Early Universe. Following these cosmic blueprints, the ``Constructors'' set out to present the ultimate realization of the cosmic framework. This concerns the search for the best available estimates of the parameters characterizing the Friedmann-Robertson-Walker metric. The specification of these cosmic measures paves the way for the ``Interior Decorators'' to take care of the cosmic ``Infrastructure''. This involves the patterns and morphology of Megaparsec scale cosmic structure and its emergence through gravitational instability, as well as aspects concerning the rise of galaxies through the complex interplay of gravitational, hydrodynamic, radiative, stellar and other dissipative processes.

We stress the importance of stochastic geometry as a branch of mathematical statistics particularly suited to model and investigate nontrivial spatial patterns. One of its key concepts, Voronoi tessellations, represents a versatile and flexible mathematical model for foamlike patterns. Based on a seemingly simple definition, Voronoi tessellations define a wealthy stochastic network of interconnected anisotropic components, each of which can be identified with the various structural elements of the cosmic galaxy distribution. Here we describe results of an ongoing investigation of a variety of aspects of cosmologically relevant spatial distributions and statistics within the framework of Voronoi tessellations. Interesting is the recent finding of a profound scaling of both clustering strength and clustering extent for the distribution of tessellation nodes, representing galaxy clusters. It suggests a hitherto unexpected fundamental and profound property of foamlike geometries. In a sense, cellular networks may be the source of an intrinsic ``geometrically biased'' clustering.

We assess the extent to which the flux-limited PSCz redshift sample is capable of accounting for the major share of mass concentrations inducing the external tidal forces affecting the peculiar velocities within the Local Supercluster (LS). The investigation is based upon a comparison of the ``true'' velocities in 2 large N-body simulations and their reconstruction from ``observation-mimicking'' mock catalogues. The mildly nonlinear ``mock'' LS and PSCz velocities are analyzed by means of the Least Action Principle technique in its highly optimized implementation of Nusser & Branchini's Fast Action Method (FAM). For both model N-body Universes, we conclude that the dipolar and quadrupolar force field implied by the PSCz galaxy distribution would indeed be sufficiently representing the full external tidal force field.

Optimally resolved one-dimensional density and velocity profiles through cosmological N-body simulations are constructed by means of the Voronoi-Delaunay tessellation reconstruction technique. In a fully self-adaptive fashion a strikingly detailed view of the density features and the corresponding cosmic motions is recovered.

Astronomical observations, physical experiments as well as computer simulations often involve discrete data sets supposed to represent a fair sample of an underlying smooth and continuous field. Reconstructing the underlying fields from a set of irregularly sampled data is therefore a recurring key issue in operations on astronomical data sets. Conventional methods involve artificial filtering through a grid or a smoothing kernel and fail to achieve an optimal result. Here we describe a fully self-adaptive geometric method which does not make use of artificial filtering, and which makes optimal use of the available information.

The application of Voronoi and Delaunay tessellation based methods for reconstructing continuous fields from discretely sampled data sets is discussed. The succesfull operation as ``multidimensional interpolation'' method is corroborated through their ability to reproduce even intricate statistical aspects of the analytical predictions of perturbation theory of cosmic velocity field evolution. The newly developed and fully self-adaptive technique for density field estimation by means of Delaunay tessellations, basically exploiting their ``minimum triangulation'' quality, is shown to succesfully reproduce the morphology of the foamlike structure in N-body simulations of cosmic structure formation. The full hierarchy of structure is implicitly and directly reproduced at every spatial resolution scale present in the particle distribution, at the same time automatically and sharply rendering its characteristic anisotropic filamentary and wall-like features.

Here we introduce the Delaunay Density Estimator Method. Its purpose is rendering a fully volume-covering reconstruction of a density field from a set of discrete data points sampling this field. Reconstructing density or intensity fields from a set of irregularly sampled data is a recurring key issue in operations on astronomical data sets, both in an observational context as well as in the context of numerical simulations. Our technique is based upon the stochastic geometric concept of the Delaunay tessellation generated by the point set. We shortly describe the method, and illustrate its virtues by means of an application to an N-body simulation of cosmic structure formation. The presented technique is a fully adaptive method: automatically it probes high density regions at maximum possible resolution, while low density regions are recovered as moderately varying regions devoid of the often irritating shot-noise effects. Of equal importance is its capability to sharply and undilutedly recover anisotropic density features like filaments and walls. The prominence of such features at a range of resolution levels within a hierarchical clustering scenario as the example of the standard CDM scenario is shown to be impressively recovered by our scheme.

Voronoi Tessellations form an attractive and versatile geometrical asymptotic model for the foamlike cosmic distribution of matter and galaxies. In the Voronoi model the vertices are identified with clusters of galaxies. For a substantial range out to a scale in the order of the cellsize, their spatial two-point correlation function is a power-law with a slope $\gamma \approx 1.95$. This study presents recent results showing that subsets of vertices selected on the basis of their ``richness'', i.e. inflow volume, retain this power-law correlation behaviour. Interestingly, they do so with a ``clustering length'' $r_{\rm o}$ that is exactly linearly proportional to the average inter-vertex distance in the sample, thus forming a realization of the Szalay-Schramm prescription. For the geometry and structural patterns even more significant is the finding tessellation vertices display a similar linear increase for their correlation length $r_{\rm a}$, the coherence length at which $\xi(r_{\rm a})=0$. Such patterns therefore exhibit positive correlations out to distances considerably in excess of the cellsize. Most intriguing is the implication of self-similar scaling, while these results may be regarded as the presentation of a ``geometrical bias'' effect. A seemingly rigid structure appears to represent a flexible and useful geometric model for exploring the statistical and dynamical repercussions of the nontrivial cellular patterns in Megaparsec scale cosmic structure.

Unlike the substantial coherent bulk motion in which our local patch of the Cosmos is participating, the amplitude of the random motions around this large scale flow seems to be surprisingly low. Attempts to invoke global explanations to account for this coldness of the local cosmic velocity field have not yet been succesfull. Here we propose a different view on this cosmic dilemma, stressing the repercussions of our cosmic neighbourhood embodying a rather uncharacteristic region of the Cosmos. Suspended between two huge mass concentrations, the Great Attractor region and the Perseus-Pisces chain, we find ourselves in a region of relatively low density yet with a very strong tidal shear. By means of constrained realizations of our local Universe, based on Wiener-filtered reconstructions inferred from the Mark III catalogue of galaxy peculiar velocities, we show that indeed this configuration may induce locally cold regions. Hence, the coldness of the local flow may be a cosmic variance effect.

Several studies have indicated that the local cosmic velocity field is rather cold, in particular in the regions outside the massive, virialized clusters of galaxies. If our local cosmic environment is taken to be a representative volume of the Universe, the repercussion of this finding is that either we live in a low-$\Omega$ Universe and/or that the galaxy distribution is a biased reflection of the underlying mass distribution. Otherwise, the pronounced nature of the observed galaxy distribution would be irreconcilable with the relatively quiet flow of the galaxies. Here we propose a different view on this cosmic dilemma, stressing the fact that our cosmic neighbourhood embodies a region of rather particular dynamical properties, and henceforth we are apt to infer flawed conclusions with respect to the global Universe. Suspended between two huge mass concentrations, the Great Attractor region and the Perseus-Pisces chain, we find ourselves in a region of relatively low density yet with a very strong tidal shear. This tidal field induces a local velocity field with a significant large-scale bulk flow but a low small-scale velocity dispersion. By means of constrained realizations of our local Universe, consisting of Wiener-filtered reconstructions inferred from the Mark III catalogue of galaxy peculiar velocities in combination with appropriate spectrally determined fluctuations, we study the implications for our local velocity field. We find that we live near a local peak in the distribution of the cosmic Mach number, $|v_{bulk}|/\sigma_v$, and that our local cosmic niche is located in the tail of the Mach number distribution function.

Over the past three decades the established view of a nearly homogeneuous, featureless Universe on scales larger than a few Megaparsec has been completely overhauled. In particular through the advent of ever larger galaxy redshift surveys we were revealed a galaxy distribution displaying an intriguing cellular pattern in which filamentary and wall-like structures, as well as huge regions devoid of galaxies, are amongst the most conspicuous morphological elements. In this contribution we will provide an overview of the present observational state of affairs concerning the distribution of galaxies and the structure traced out by the matter distribution in our Universe. In conjunction with the insight on the dynamics of the structure formation process obtained through the mapping of the peculiar velocities of galaxies in our local Universe and the information on the embryonic circumstances that prevailed at the epoch of Recombination yielded by the various Cosmic Microwave Background experiments, we seek to arrive at a more or less compelling theoretical framework of structure formation.The main aspects of this framework of the rise of structure through gravitational instability can probably be most readily appreciated through illustrative examples of various scenarios, as for instance provided by some current state-of-the-art N-body simulations. We will subsequently wrap up the observational and theoretical evidence for the emergence and evolution of structure in the Universe by sketching the stage for the ultimate Holy Grail of late 20th century astrophysics, understanding the saga of the formation of what arguably are the most prominent and at the same time intoxicatingly beautiful and intriguing denizens of our Cosmos, the {\it galaxies}.

We describe two new -- stochastic-geometrical -- methods to obtain reliable velocity field statistics from N-body simulations and from any general density and velocity fluctuation field sampled at a discrete set of locations. These methods, the Voronoi tessellation method and Delaunay tessellation method, are based on the use of the Voronoi and Delaunay tessellations of the point distribution defined by the locations at which the velocity field is sampled. Adjusting themselves automatically to the density of sampling points, they represent the optimal estimator for volume-averaged quantities. They are therefore particularly suited for checking the validity of the predictions of quasi-linear analytical density and velocity field perturbation theory through the results of N-body simulations of structure formation. We illustrate the subsequent succesfull application of the two methods to estimate the bias-independent value of $\Omega$ in the N-body simulations on the basis of the predictions of perturbation theory for the $\Omega$-dependence of the moments and PDF of the velocity divergence in gravitational instability structure formation scenarios with Gaussian initial conditions. We will also shortly discuss practical and complicating issues involved in the obvious extension of the Voronoi and Delaunay method to the analysis of observational samples of galaxy peculiar velocities.

Analytical studies based on perturbative theory have shown that the moments of the Probability Distribution Function (PDF) of the local smoothed velocity divergence are expected to have a very specific dependence on the density parameter Omega in the quasi-linear regime. This dependence is particularly interesting as it does not involve the possible bias between the galaxy spatial distribution and the underlying mass distribution. This implies a new and promising method for determining a bias-independent value of Omega based on a reliable determination of the velocity divergence PDF. In this paper we study the Omega dependence of the velocity divergence PDF and its first moments in a set of N-body simulations, using the so-called Voronoi and Delaunay methods. We show that this dependence is in agreement with the theoretical prediction, even while the number density of velocity field tracers has been diluted to a value comparable to that available in current galaxy catalogues. In addition, we demonstrate that a sufficiently reliable determination of these statistical quantities is also possible when the measurement of the galaxy peculiar velocities is restricted to the one component along the line-of-sight. Under ideal, noise-free circumstances we can successfully discriminate between low and high Omega.

The cosmological case for a next generation radio observatory, the Square Kilometer Array, is discussed and reviewed. An instrument like the SKA would be able to measure galaxy redshifts of normal late-type galaxies, via the 21 cm line of HI, out to redshifts of $\sim 3$. Not only would such very deep redshift surveys enable us to map the large scale galaxy distribution and probe the large scale structure of the universe out to previously unexplored scales, it would also allow for the first time to obtain direct observational data on the evolution of this structure. Other promising applications concern the mapping of the local velocity field of the universe, study of the formation and evolution of galaxies, and determining the global cosmological parameters $H_0$, $q_0$ and $\Lambda$ through the application of classical cosmological tests like source counts. Particularly emphasized is the redshift survey capability of the SKA. A review is given of the current knowledge of the galaxy distribution, starting from an inventarisation of nearby cosmic structures, through a discussion of how it all fits together in a coherent ``foamlike pattern''. After providing a short overview of the basics of theories of structure formation, a description is provided of different observational strategies to probe the structure of the universe out to larger depths, ranging from pencil-beam surveys and cluster surveys out to the new and ambitious complete and deep galaxy redshift surveys like the 2dF and the Sloan survey. It is argued that a survey with the SKA would be a natural and complementary follow-up. We finally conclude with a specification of the technical requirements for the SKA to make it into an instrument ideally suited for these purposes.

The role of tidal shear in the formation of structure in the Universe is explored. To illustrate the possible and sometimes dramatic impact of tidal fields we focus on the evolution of voids. We firstly analyze the role of tidal fields both in the highly symmetric situation of an initially homogeneous ellipsoidal underdensity embedded in an artificially imposed tidal field. In addition, we present selfconsistent case studies consisting of N-body simulations that start from constrained Gaussian initial conditions in which the matter distribution is appropriately sculpted to yield an a priori specified tidal field. We conclude that voids may indeed be induced to collapse. Also, we present evidence for the strong relation between tidal fields and filaments in the mass distribution.

We introduce two new methods to obtain reliable velocity field statistics from N-body simulations, or indeed from any general density and velocity fluctuation field sampled by discrete points. These methods, the {\it Voronoi tessellation method} and {\it Delaunay tessellation method}, are based on the use of the Voronoi and Delaunay tessellations of the point distribution defined by the locations at which the velocity field is sampled. In the Voronoi method the velocity is supposed to be uniform within the Voronoi polyhedra, whereas the Delaunay method constructs a velocity field by linear interpolation between the four velocities at the locations defining each Delaunay tetrahedron. The most important advantage of these methods is that they provide an optimal estimator for determining the statistics of volume-averaged quantities, as opposed to the available numerical methods that mainly concern mass-averaged quantities. As the major share of the related analytical work on velocity field statistics has focussed on volume-averaged quantities, the availability of appropriate numerical estimators is of crucial importance for checking the validity of the analytical perturbation calculations. In addition, it allows us to study the statistics of the velocity field in the highly nonlinear clustering regime. Although in most cases both the Voronoi and the Delaunay method lead to equally good velocity field estimates, the Delaunay method may be slightly preferable. In particular it performs considerably better at small radii. Note that it is more CPU-time intensive while its requirement for memory space is almost a factor 8 lower than the Voronoi method.

We develop an algorithm for setting up initial Gaussian random density and velocity fields containing one or more peaks or dips, in an arbitrary cosmological scenario. The intention is to generate appropriate initial conditions for cosmological N-body simulations that focus on the evolution of the progenitors of the present-day galaxies and clusters. The procedure is an application of the direct and accurate prescription of Hoffman & Ribak (1991) for generating constrained random fields. For each peak a total of 21 physical characteristics can be specified, including its scale, position, density Hessian, velocity, and velocity gradient. The velocity (or, equivalently, gravity) field constrants are based on a generalization of the formalism developed by Bardeen et al. (1986). The resulting density field is sculpted such that it induces the desired amount of net gravitational and tidal forces. We provide a detailed mathematical presentation of the formalism. Afterwards we provide analytical estimates of the likelihood of the imposed constraints. Amongst others, it is shown that the tidal field has a strong tendency to align itself along the principal axes of the mass tensor. The method is illustrated by means of some concrete examples. In addition to the illustration of constraint-field correlation functions and how they add up to the mean fields, followed by illustrations of the variance characteristics of field realizations, we concentrate in particular on the consequences of imposing gravitational field constraints (or, equivalent in the linear regime for growing mode fluctuations, peculiar velocity field constraints).

Under the assumption of a void-filled Universe we investigate if the characteristic scale of voids can be determined from existing surveys. We use the Voronoi tessellation to create mock surveys and study the properties of the first zero-crossing of the two-point correlation function for various survey geometries. Our main conclusion is that the available data sets should be able to discriminate between 5000 \kms and 12000 \kms voids, if one of these scales actually characterizes the distribution of galaxies.

We present preliminary results of our investigation into the influence of shear fields on the evolution of galactic scale fluctuations in a primordial Gaussian random density field. Specifically, we study how the matter associated with a galaxy-scale peak evolves, to determine whether the shear can affect the peak's ability to form a virialized structure. We find that the evolution of the mass distribution in the immediate surroundings of $\geq 1\sigma$ initial density peaks is sensitive to the nature and magnitude of the shear. Its final fate and configuration spans a plethora of possibilities, ranging from accretion onto the peak to complete disruption. On the other hand, the mass defining the peak itself always tends to form a virialized object, though a given galaxy-scale peak need not necessarily form only one halo. Furthermore, galaxy-size halos need not necessarily be associated with initial galaxy-scale density peaks. Under certain conditions, the shear field is capable of breaking up a single primordial peak into two (and perhaps, more) distinct halos, or of promoting the growth of smaller-scale peaks into galaxy-size halos.