We introduce a method for generating a continuous, mass-conserving and high-order differentiable density field from a discrete point distribution such as particles or halos from an N-body simulation or galaxies from a spectroscopic survey. The method consists on generating an ensemble of point realizations by perturbing the original point set following the geometric constraints imposed by the Delaunay tessellation in the vicinity of each point in the set. By computing the mean field of the ensemble we are able to significantly reduce artifacts arising from the Delaunay tessellation in poorly sampled regions while conserving the features in the point distribution. Our implementation is based on the Delaunay Tessellation Field Estimation (DTFE) method, however other tessellation techniques are possible. The method presented here shares the same advantages of the DTFE method such as self-adaptive scale, mass conservation and continuity, while being able to reconstruct even the faintest structures of the point distribution usually dominated by artifacts in Delaunay-based methods. Additionally, we also present preliminary results of an application of this method to image denoising and artifact removal, highlighting the broad applicability of the technique introduced here.

The standard explanation for galaxy spin starts with the tidal-torque theory (TTT), in which an ellipsoidal dark-matter protohalo, which comes to host the galaxy, is torqued up by the tidal gravitational field around it. We discuss a complementary picture, using the relatively familiar velocity field, instead of the tidal field, whose intuitive connection to the surrounding, possibly faraway matter arrangement is more obscure. In this 'spin from primordial inner motions' (SPIM) concept, implicit in TTT derivations but not previously emphasized, the angular momentum from the gravity-sourced velocity field inside a protohalo largely cancels out, but has some excess from the aspherical outskirts. At first, the net spin scales according to linear theory, a sort of comoving conservation of familiar angular momentum. Then, at collapse, it is conserved in physical coordinates. Small haloes are then typically subject to secondary exchanges of angular momentum. The TTT is useful for analytic estimates. But a literal interpretation of the TTT is inaccurate in detail, without some implicit concepts about smoothing of the velocity and tidal fields. This could lead to misconceptions, for those first learning about how galaxies come to spin. Protohaloes are not perfectly ellipsoidal and do not uniformly torque up, as in a naive interpretation of the TTT; their inner velocity fields retain substantial dispersion. Furthermore, quantitatively, given initial conditions and protohalo boundaries, SPIM is more direct and accurate than the TTT to predict halo spins. We also discuss how SPIM applies to rotating filaments, and the relation between halo mass and spin, in which the total spin of a halo can be thought of as a sum of random contributions.

We study the effects of the local environmental density and the cosmic web environment (filaments, walls, and voids) on key properties of dark matter halos using the Bolshoi-Planck LCDM cosmological simulation. The z = 0 simulation is analysed into filaments, walls, and voids using the SpineWeb method and also the VIDE package of tools, both of which use the watershed transform. The key halo properties that we study are the specific mass accretion rate, spin parameter, concentration, prolateness, scale factor of the last major merger, and scale factor when the halo had half of its z = 0 mass. For all these properties, we find that there is no discernible difference between the halo properties in filaments, walls, or voids when compared at the same environmental density. As a result, we conclude that environmental density is the core attribute that affects these properties. This conclusion is in line with recent findings that properties of galaxies in redshift surveys are independent of their cosmic web environment at the same environmental density at z ~ 0. We also find that the local web environment of the Milky Way and the Andromeda galaxies near the centre of a cosmic wall does not appear to have any effect on the properties of these galaxies' dark matter halos except for their orientation, although we find that it is rather rare to have such massive halos near the centre of a relatively small cosmic wall.

We propose the Cosmic Web Detachment (CWD) model, a framework to interpret the star-formation history of galaxies in a cosmological context. The CWD model unifies several starvation mechanisms known to disrupt or stop star formation into one single physical framework. Galaxies begin accreting star-forming gas at early times via a network of primordial filaments, simply related to the pattern of density fluctuations in the initial conditions. But when shell-crossing occurs on intergalactic scales, this pattern is disrupted, and the galaxy detaches from its primordial filaments, ending the accretion of cold gas. We argue that CWD encompasses known external processes halting star formation, such as harassment, strangulation and starvation. On top of these external processes, internal feedback processes such as AGN contribute to stop in star formation as well. By explicitly pointing out the non-linear nature of CWD events we introduce a simple formalism to identify CWD events in N-body simulations. With it we reproduce and explain, in the context of CWD, several observations including downsizing, the cosmic star formation rate history, the galaxy mass-color diagram and the dependence of the fraction of red galaxies with mass and local density.

Estimating the parameters of a model describing a set of observations using a neural network is in general solved in a supervised way. In cases when we do not have access to the model's true parameters this approach can not be applied. Standard unsupervised learning techniques on the other hand, do not produce meaningful or semantic representations that can be associated to the model's parameters. Here we introduce a self-supervised hybrid network that combines traditional neural network elements with analytic or numerical models which represent a physical process to be learned by the system. Self-supervised learning is achieved by generating an internal representation equivalent to the parameters of the physical model. This semantic representation is used to evaluate the model and compare it to the input data during training. The Semantic Autoencoder architecture described here shares the robustness of neural networks while including an explicit model of the data, learns in an unsupervised way and estimates, by construction, parameters with direct physical interpretation. As an illustrative application we perform unsupervised learning for 2D model fitting of exponential light profiles.

We present a new technique that allows us to compute ensemble statistics on a local basis, directly relating halo properties to their local environment. This is achieved by the use of a correlated ensemble in which the Large Scale Structure is common to all realizations while having each an independent halo population. The correlated ensemble can be stacked, effectively increasing the halo number density by an arbitrary factor, thus breaking the fundamental limit in the halo number density given by the haloe mass function. This technique allows us to compute \textit{local ensemble statistics} of the matter/haloe distribution at any position in the simulation box, while removing the intrinsic stochasticity in the halo formation process and directly relating halo properties to their environment. We introduce the \textit{Multum In Parvo} (MIP) correlated ensemble simulation consisting of 220 realizations on a 32 h$^{-1}$ Mpc box with $256^3$ particles each. This is equivalent in terms of effective volume and number of particles to a box of $\sim 193$ h$^{-1}$ Mpc of side with $\sim 1540^3$ particles containing $\sim 5\times 10^6$ haloes with a minimum mass of $3.25 \times 10^9$ h$^{-1}$ M$_{\odot}$. The potential of the technique presented here is illustrated by computing the local ensemble statistics of the halo ellipticity and halo shape-LSS alignment. We show that, while there are general trends in the ellipticity and alignment of haloes with their LSS, there are also significant spatial variations which has important implications for observational studies of galaxy shape and alignment.

We present the first application of deep neural networks to the semantic segmentation of cosmological filaments and walls in the Large Scale Structure of the Universe. Our results are based on a deep Convolutional Neural Network (CNN) with a U-Net architecture trained using an existing state-of-the-art manually-guided segmentation method. We successfully trained an tested an U-Net with a Voronoi model and an N-body simulation. The predicted segmentation masks from the Voronoi model have a Dice coefficient of 0.95 and 0.97 for filaments and mask respectively. The predicted segmentation masks from the N-body simulation have a Dice coefficient of 0.78 and 0.72 for walls and filaments respectively. The relatively lower Dice coefficient in the filament mask is the result of filaments that were predicted by the U-Net model but were not present in the original segmentation mask. Our results show that for a well-defined dataset such as the Voronoi model the U-Net has excellent performance. In the case of the N-body dataset the U-Net produced a filament mask of higher quality than the segmentation mask obtained from a state-of-the art method. The U-Net performs better than the method used to train it, being able to find even the tenuous filaments that the manually-guided segmentation failed to identify. The U-Net presented here can process a $512^3$ volume in a few minutes and without the need of complex pre-processing. Deep CNN have great potential as an efficient and accurate analysis tool for the next generation large-volume computer N-body simulations and galaxy surveys.

The cosmic web plays a major role in the formation and evolution of galaxies and defines, to a large extent, their properties. However, the relation between galaxies and environment is still not well understood. Here we present a machine learning approach to study imprints of environmental effects on the mass assembly of haloes. We present a galaxy-LSS machine learning classifier based on galaxy properties sensitive to the environment. We then use the classifier to assess the relevance of each property. Correlations between galaxy properties and their cosmic environment can be used to predict galaxy membership to void/wall or filament/cluster with an accuracy of $93\%$. Our study unveils environmental information encoded in properties of haloes not normally considered directly dependent on the cosmic environment such as merger history and complexity. Understanding the physical mechanism by which the cosmic web is imprinted in a halo can lead to significant improvements in galaxy formation models. This is accomplished by extracting features from galaxy properties and merger trees, computing feature scores for each feature and then applying support vector machine to different feature sets. To this end, we have discovered that the shape and depth of the merger tree, formation time and density of the galaxy are strongly associated with the cosmic environment. We describe a significant improvement in the original classification algorithm by performing LU decomposition of the distance matrix computed by the feature vectors and then using the output of the decomposition as input vectors for support vector machine.

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 introduce NebulOS, a Big Data platform that allows a cluster of Linux machines to be treated as a single computer. With NebulOS, the process of writing a massively parallel program for a datacenter is no more complicated than writing a Python script for a desktop computer. The platform enables most pre-existing data analysis software to be used, as scale, in a datacenter without modification. The shallow learning curve and compatibility with existing software greatly reduces the time required to develop distributed data analysis pipelines. The platform is built upon industry-standard, open-source Big Data technologies, from which it inherits several fault tolerance features. NebulOS enhances these technologies by adding an intuitive user interface, automated task monitoring, and other usability features. We present a summary of the architecture, provide usage examples, and discuss the system's performance scaling.

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.

Warm dark matter (WDM) has been proposed as an alternative to cold dark matter (CDM), to resolve issues such as the apparent lack of satellites around the Milky Way. Even if WDM is not the answer to observational issues, it is essential to constrain the nature of the dark matter. The effect of WDM on haloes has been extensively studied, but the small-scale initial smoothing in WDM also affects the present-day cosmic web and voids. It suppresses the cosmic "sub-web" inside voids, and the formation of both void haloes and subvoids. In $N$-body simulations run with different assumed WDM masses, we identify voids with the ZOBOV algorithm, and cosmic-web components with the ORIGAMI algorithm. As dark-matter warmth increases (i.e., particle mass decreases), void density minima grow shallower, while void edges change little. Also, the number of subvoids decreases. The density field in voids is particularly insensitive to baryonic physics, so if void density profiles and minima could be measured observationally, they would offer a valuable probe of the nature of dark matter. Furthermore, filaments and walls become cleaner, as the substructures in between have been smoothed out; this leads to a clear, mid-range peak in the density PDF.

It is well-known that galaxy environment has a fundamental effect in shaping its properties. We study the environmental effects on galaxy evolution, with an emphasis on the environment defined as the local number density of galaxies. The density field is estimated with different estimators (weighted adaptive kernel smoothing, 10$^{th}$ and 5$^{th}$ nearest neighbors, Voronoi and Delaunay tessellation) for a K$_{s}<$24 sample of $\sim$190,000 galaxies in the COSMOS field at 0.1$<$z$<$3.1. The performance of each estimator is evaluated with extensive simulations. We show that overall, there is a good agreement between the estimated density fields using different methods over $\sim$2 dex in overdensity values. However, our simulations show that adaptive kernel and Voronoi tessellation outperform other methods. Using the Voronoi tessellation method, we assign surface densities to a mass complete sample of quiescent and star-forming galaxies out to z$\sim$3. We show that at a fixed stellar mass, the median color of quiescent galaxies does not depend on their host environment out to z$\sim$3. We find that the number and stellar mass density of massive ($>$10$^{11}$M$_{\odot}$) star-forming galaxies have not significantly changed since z$\sim$3, regardless of their environment. However, for massive quiescent systems at lower redshifts (z$\lesssim$1.3), we find a significant evolution in the number and stellar mass densities in denser environments compared to lower density regions. Our results suggest that the relation between stellar mass and local density is more fundamental than the color-density relation and that environment plays a significant role in quenching star formation activity in galaxies at z$\lesssim$1.

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 star formation history of galaxies is a complex process usually considered to be stochastic in nature, for which we can only give average descriptions such as the color-density relation. In this work we follow star-forming gas particles in a hydrodynamical N-body simulation back in time in order to study their initial spatial configuration. By keeping record of the time when a gas particle started forming stars we can produce gas-star isochrone surfaces delineating the surfaces of accreting gas that begin producing stars at different times. These accretion surfaces are closely packed inside dense regions, intersecting each other, and as a result galaxies inside proto-clusters stop accreting gas early, naturally explaining the color dependence on density. The process described here has a purely gravitational / geometrical origin, arguably operating at a more fundamental level than complex processes such as AGN and supernovae, and providing a conceptual origin for the color-density relation.

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.

As an alternative way of describing the cosmological velocity field, we discuss the evolution of rotational invariants constructed from the velocity gradient tensor. Compared with the traditional divergence-vorticity decomposition, these invariants, defined as coefficients of characteristic equation of the velocity gradient tensor, enable a complete classification of all possible flow patterns in the dark-matter comoving frame, including both potential and vortical flows. We show that this tool, first introduced in turbulence two decades ago, proves to be very useful in understanding the evolution of the cosmic web structure, and in classifying its morphology. Before shell-crossing, different categories of potential flow are highly associated with cosmic web structure, because of the coherent evolution of density and velocity. This correspondence is even preserved at some level when vorticity is generated after shell-crossing. The evolution from the potential to vortical flow can be traced continuously by these invariants. With the help of this tool, we show that the vorticity is generated in a particular way that is highly correlated with the large-scale structure. This includes a distinct spatial distribution and different types of alignment between cosmic web and vorticity direction for various vortical flows. Incorporating shell-crossing into closed dynamical systems is highly non-trivial, but we propose a possible statistical explanation for some of these phenomena relating to the internal structure of the three-dimensional invariants space.

We study the relationship between dark-matter haloes and matter in the MIP $N$-body simulation ensemble, which allows precision measurements of this relationship, even deeply into voids. What enables this is a lack of discreteness, stochasticity, and exclusion, achieved by averaging over hundreds of possible sets of initial small-scale modes, while holding fixed large-scale modes that give the cosmic web. We find (i) that dark-matter-halo formation is greatly suppressed in voids; there is an exponential downturn at low densities in the otherwise power-law matter-to-halo density bias function. Thus, the rarity of haloes in voids is akin to the rarity of the largest clusters, and their abundance is quite sensitive to cosmological parameters. The exponential downturn appears both in an excursion-set model, and in a model in which fluctuations evolve in voids as in an open universe with an effective $\Omega_m$ proportional to a large-scale density. We also find that (ii) haloes typically populate the average halo-density field in a super-Poisson way, i.e. with a variance exceeding the mean; and (iii) the rank-order-Gaussianized halo and dark-matter fields are impressively similar in Fourier space. We compare both their power spectra and cross-correlation, supporting the conclusion that one is roughly a strictly-increasing mapping of the other. The MIP ensemble especially reveals how halo abundance varies with `environmental' quantities beyond the local matter density; (iv) we find a visual suggestion that at fixed matter density, filaments are more populated by haloes than clusters.

Dark matter haloes in cosmological filaments and walls have their spin vector aligned (in average) with their host structure. While haloes in walls are aligned with the plane of the wall independently of their mass, haloes in filaments present a mass dependent two-regime orientation. Here we show that the transition mass determining the change in the alignment regime (from parallel to perpendicular) depends on the hierarchical level in which the halo is located, reflecting the hierarchical nature of the Cosmic Web. By explicitly exposing this hierarchy we are able to identify the contributions of different components of the filament network to the spin alignment signal. We discuss a unifying picture to describe the alignment of haloes in filaments and walls consistent with previous results and our findings based on a two-phase angular momentum acquisition, first via tidal torquening and later via anisotropic mass accretion. The hierarchical identification and characterization of cosmic structures was done with a new implementation of the Multiscale Morphology Filter. The MMF-2 (introduced here) represents a significant improvement over its predecessor in terms of robustness and by explicitly describing the hierarchical relations between the elements of the Cosmic Web.

We explore the possibility of detecting radio emission in the \emph{cosmic web} by analyzing shock waves in the MareNostrum cosmological simulation. This requires a careful calibration of shock finding algorithms in Smoothed-Particle Hydrodynamics simulations, which we present here. Moreover, we identify the elements of the cosmic web, namely voids, walls, filaments and clusters with the use of the SpineWeb technique, a procedure that classifies the structure in terms of its topology. Thus, we are able to study the Mach number distribution as a function of its environment. We find that the median Mach number, for clusters is $\mathcal{M}_{\mathrm{clusters}}\approx1.8$, for filaments is $\mathcal{M}_{\mathrm{filaments}}\approx 6.2$, for walls is $\mathcal{M}_{\mathrm{walls}}\approx 7.5$, and for voids is $\mathcal{M}_{\mathrm{voids}}\approx 18$. We then estimate the radio emission in the cosmic web using the formalism derived in Hoeft & Br\"{u}ggen (2007). We also find that in order to match our simulations with observational data (e.g., NVSS radio relic luminosity function), a fraction of energy dissipated at the shock of $\xi_{\mathrm{e}}=0.0005$ is needed, in contrast with the $\xi_{\mathrm{e}}=0.005$ proposed by Hoeft et al. (2008). We find that 41% of clusters with $M \ge 10^{14} M_{\odot}$ host diffuse radio emission in the form of radio relics. Moreover, we predict that the radio flux from filaments should be $S_{150 MHz}\sim 0.12$ $\mu$Jy at a frequency of 150 MHz.

Contrary to the common view voids have very complex internal structure and dynamics. Here we show how the hierarchy of structures in the density field inside voids is reflected by a similar hierarchy of structures in the velocity field. Voids defined by dense filaments and clusters can de described as simple expanding domains with coherent flows everywhere except at their boundaries. At scales smaller that the void radius the velocity field breaks into expanding sub-domains corresponding to sub- voids. These sub-domains break into even smaller sub-sub domains at smaller scales resulting in a nesting hierarchy of locally expanding domains. The ratio between the magnitude of the velocity field responsible for the expansion of the void and the velocity field defining the sub voids is approximately one order of magnitude. The small-scale components of the velocity field play a minor role in the shaping of the voids but they define the local dynamics directly affecting the processes of galaxy formation and evolution. The super-Hubble expansion inside voids makes them cosmic magnifiers by stretching their internal primordial density fluctuations allowing us to probe the small scales in the primordial density field. Voids also act like time machines by "freezing" the development of the medium-scale density fluctuations responsible for the formation of the tenuous web of structures seen connecting proto galaxies in computer simulations. As a result of this freezing haloes in voids can remain "connected" to this tenuous web until the present time. This may have an important effect in the formation and evolution of galaxies in voids by providing an efficient gas accretion mechanism via coherent low-velocity streams that can keep a steady inflow of matter for extended periods of time.

We investigate the use of a logarithmic density variable in estimating the Lagrangian displacement field, motivated by the success of a logarithmic transformation in restoring information to the matter power spectrum. The logarithmic relation is an extension of the linear relation, motivated by the continuity equation, in which the density field is assumed to be proportional to the divergence of the displacement field; we compare the linear and logarithmic relations by measuring both of these fields directly in a cosmological N-body simulation. The relative success of the logarithmic and linear relations depends on the scale at which the density field is smoothed. Thus we explore several ways of measuring the density field, including Cloud-In-Cell smoothing, adaptive smoothing, and the (scale-independent) Delaunay tessellation, and we use both a Fourier space and a geometrical tessellation approach to measuring the divergence. We find that the relation between the divergence of the displacement field and the density is significantly tighter with a logarithmic density variable, especially at low redshifts and for very small (~2 Mpc/h) smoothing scales. We find that the grid-based methods are more reliable than the tessellation-based method of calculating both the density and the divergence fields, though in both cases the logarithmic relation works better in the appropriate regime, which corresponds to nonlinear scales for the grid-based methods and low densities for the tessellation-based method.

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.

[abridged] We present a detailed comparison of fundamental dark matter halo properties retrieved by a substantial number of different halo finders. These codes span a wide range of techniques including friends-of-friends (FOF), spherical-overdensity (SO) and phase-space based algorithms. We further introduce a robust (and publicly available) suite of test scenarios that allows halo finder developers to compare the performance of their codes against those presented here. This set includes mock haloes containing various levels and distributions of substructure at a range of resolutions as well as a cosmological simulation of the large-scale structure of the universe. All the halo finding codes tested could successfully recover the spatial location of our mock haloes. They further returned lists of particles (potentially) belonging to the object that led to coinciding values for the maximum of the circular velocity profile and the radius where it is reached. All the finders based in configuration space struggled to recover substructure that was located close to the centre of the host halo and the radial dependence of the mass recovered varies from finder to finder. Those finders based in phase space could resolve central substructure although they found difficulties in accurately recovering its properties. Via a resolution study we found that most of the finders could not reliably recover substructure containing fewer than 30-40 particles. However, also here the phase space finders excelled by resolving substructure down to 10-20 particles. By comparing the halo finders using a high resolution cosmological volume we found that they agree remarkably well on fundamental properties of astrophysical significance (e.g. mass, position, velocity, and peak of the rotation curve).

We show that the "cold" Hubble flow observed for galaxies around the Milky Way does not represent a problem in cosmology but is due to the particular geometry and dynamics of our local wall. The behavior of the perturbed Hubble flow around the Milky Way is the result of two main factors: at small scales (R < 1 Mpc) the inflow is dominated by the gravitational influence of the Milky Way. At large scales (R > 1 Mpc) the out flow reflects the expansion of our local wall which "cools down" the peculiar velocities. This is an intrinsic property of walls and is independent of cosmology. We find the dispersion of the local Hubble flow (1 < R < 3 Mpc) around simulated "Milky Way" haloes located at the centre of low-density cosmological walls to be {\sigma}_H ~ 30 km/s, in excellent agreement with observations. The expansion of our local wall is also reflected in the value of the measured local Hubble constant. For "Milky Way" haloes inside walls, we find super-Hubble flows with h_local \simeq 0.77 - 1.13. The radius of equilibrium (R_0) depends not only on the mass of the central halo and the Hubble expansion but also on the dynamics given by the local LSS geometry. The super-Hubble flow inside our local wall has the effect of reducing the radius at which the local expansion balances the gravitational influence of the Milky Way. By ignoring the dynamical effect of the local wall, the mass of the Milky Way estimated from R_0 can be underestimated by as much as ~ 30%.

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

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

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 present a novel method for computing the Minkowski Functionals from isodensity surfaces extracted directly from the Delaunay tessellation of a point distribution. This is an important step forward compared to the previous cosmological studies when the isodensity surface was built in the field on a uniform cubic grid and therefore having a uniform spatial resolution. The density field representing a particular interest in cosmology is the density of galaxies which is obtained from the highly nonuniform distribution of the galaxy positions. Therefore, the constraints caused by the spatially uniform grid put severe limitations on the studies of the geometry and shapes of the large-scale objects: superclusters and voids of galaxies. Our technique potentially is able to eliminate most of these limitations. The method is tested with some simple geometric models and an application to the density field from an N-body simulation is shown.

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.

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

We measure the projected spatial correlation function w_p(r_p) from a large sample combining GALEX ultraviolet imaging with the SDSS spectroscopic sample. We study the dependence of the clustering strength for samples selected on (NUV - r)_abs color, specific star formation rate (SSFR), and stellar mass. We find that there is a smooth transition in the clustering of galaxies as a function of this color from weak clustering among blue galaxies to stronger clustering for red galaxies. The clustering of galaxies within the "green valley" has an intermediate strength, and is consistent with that expected from galaxy groups. The results are robust to the correction for dust extinction. The comparison with simple analytical modeling suggests that the halo occupation number increases with older star formation epochs. When splitting according to SSFR, we find that the SSFR is a more sensitive tracer of environment than stellar mass.

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.}