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 present the first quantitative detection of large-scale filamentary structure at $z \simeq 0.7$ in the large cosmological volume probed by the VIMOS Public Extragalactic Redshift Survey (VIPERS). We use simulations to show the capability of VIPERS to recover robust topological features in the galaxy distribution, in particular the filamentary network. We then investigate how galaxies with different stellar masses and stellar activities are distributed around the filaments and find a significant segregation, with the most massive or quiescent galaxies being closer to the filament axis than less massive or active galaxies. The signal persists even after down-weighting the contribution of peak regions. Our results suggest that massive and quiescent galaxies assemble their stellar mass through successive mergers during their migration along filaments towards the nodes of the cosmic web. On the other hand, low-mass star-forming galaxies prefer the outer edge of filaments, a vorticity rich region dominated by smooth accretion, as predicted by the recent spin alignment theory. This emphasizes the role of large scale cosmic flows in shaping galaxy properties.

Resolving numerically Vlasov-Poisson equations for initially cold systems can be reduced to following the evolution of a three-dimensional sheet evolving in six-dimensional phase-space. We describe a public parallel numerical algorithm consisting in representing the phase-space sheet with a conforming, self-adaptive simplicial tessellation of which the vertices follow the Lagrangian equations of motion. The algorithm is implemented both in six- and four-dimensional phase-space. Refinement of the tessellation mesh is performed using the bisection method and a local representation of the phase-space sheet at second order relying on additional tracers created when needed at runtime. In order to preserve in the best way the Hamiltonian nature of the system, refinement is anisotropic and constrained by measurements of local Poincar\'e invariants. Resolution of Poisson equation is performed using the fast Fourier method on a regular rectangular grid, similarly to particle in cells codes. To compute the density projected onto this grid, the intersection of the tessellation and the grid is calculated using the method of Franklin and Kankanhalli (1993) generalised to linear order. As preliminary tests of the code, we study in four dimensional phase-space the evolution of an initially small patch in a chaotic potential and the cosmological collapse of a fluctuation composed of two sinusoidal waves. We also perform a "warm" dark matter simulation in six-dimensional phase-space that we use to check the parallel scaling of the code.

We perform a detailed comparison of the phase-space density traced by the particle distribution in Gadget simulations to the result obtained with a spherical Vlasov solver using the splitting algorithm. The systems considered are apodized H\'enon spheres with two values of the virial ratio, R ~ 0.1 and 0.5. After checking that spherical symmetry is well preserved by the N-body simulations, visual and quantitative comparisons are performed. In particular we introduce new statistics, correlators and entropic estimators, based on the likelihood of whether N-body simulations actually trace randomly the Vlasov phase-space density. When taking into account the limits of both the N-body and the Vlasov codes, namely collective effects due to the particle shot noise in the first case and diffusion and possible nonlinear instabilities due to finite resolution of the phase-space grid in the second case, we find a spectacular agreement between both methods, even in regions of phase-space where nontrivial physical instabilities develop. However, in the colder case, R=0.1, it was not possible to prove actual numerical convergence of the N-body results after a number of dynamical times, even with N=10$^8$ particles.

The kinematic analysis of dark matter and hydrodynamical simulations suggests that the vorticity in large-scale structure is mostly confined to, and predominantly aligned with their filaments, with an excess of probability of 20 per cent to have the angle between vorticity and filaments direction lower than 60 degrees relative to random orientations. The cross sections of these filaments are typically partitioned into four quadrants with opposite vorticity sign, arising from multiple flows, originating from neighbouring walls. The spins of halos embedded within these filaments are consistently aligned with this vorticity for any halo mass, with a stronger alignment for the most massive structures up to an excess of probability of 165 per cent. On large scales, adiabatic/cooling hydrodynamical simulations display the same vorticity in the gas as in the dark matter. The global geometry of the flow within the cosmic web is therefore qualitatively consistent with a spin acquisition for smaller halos induced by this large-scale coherence, as argued in Codis et al. (2012). In effect, secondary anisotropic infall (originating from the vortex-rich filament within which these lower-mass halos form) dominates the angular momentum budget of these halos. The transition mass from alignment to orthogonality is related to the size of a given multi-flow region with a given polarity. This transition may be reconciled with the standard tidal torque theory if the latter is augmented so as to account for the larger scale anisotropic environment of walls and filaments.

A large-scale hydrodynamical cosmological simulation, Horizon-AGN, is used to investigate the alignment between the spin of galaxies and the cosmic filaments above redshift 1.2. The analysis of more than 150 000 galaxies per time step in the redshift range 1.2<z<1.8 with morphological diversity shows that the spin of low-mass blue galaxies is preferentially aligned with their neighbouring filaments, while high-mass red galaxies tend to have a perpendicular spin. The reorientation of the spin of massive galaxies is provided by galaxy mergers, which are significant in their mass build-up. We find that the stellar mass transition from alignment to misalignment happens around 3.10^10 M_sun. Galaxies form in the vorticity-rich neighbourhood of filaments, and migrate towards the nodes of the cosmic web as they convert their orbital angular momentum into spin. The signature of this process can be traced to the properties of galaxies, as measured relative to the cosmic web. We argue that a strong source of feedback such as active galactic nuclei is mandatory to quench in situ star formation in massive galaxies and promote various morphologies. It allows mergers to play their key role by reducing post-merger gas inflows and, therefore, keeping spins misaligned with cosmic filaments.

We performed five pointing observations with Suzaku to search for hot gases associated with the junctions of galaxy filaments where no significant diffuse X-ray sources were detected so far. We discovered X-ray sources successfully in all five regions and analyzed two bright sources in each field. Spectral analysis indicates that three sources originate from X-ray diffuse halos associated with optically bright galaxies or groups of galaxies with kT~0.6-0.8 keV. Other three sources are possibly group- and cluster-scale X-ray halos with temperatures of ~1 keV and ~4 keV, respectively while the others are compact object origins such as AGNs. All the three observed intracluster media within the junctions of the galaxy filaments previously found are involved in ongoing mergers. Thus, we demonstrate that deep X-ray observations at the filament junctions identified by galaxy surveys are a powerful mean to explore growing halos in a hierarchical structure undetected so far.

Pillars and globules are present in many high-mass star-forming regions, such as the Eagle nebula (M16) and the Rosette molecular cloud, and understanding their origin will help characterize triggered star formation. The formation mechanisms of these structures are still being debated. Recent numerical simulations have shown how pillars can arise from the collapse of the shell in on itself and how globules can be formed from the interplay of the turbulent molecular cloud and the ionization from massive stars. The goal here is to test this scenario through recent observations of two massive star-forming regions, M16 and Rosette. The column density structure of the interface between molecular clouds and H ii regions was characterized using column density maps obtained from far-infrared imaging of the Herschel HOBYS key programme. Then, the DisPerSe algorithm was used on these maps to detect the compressed layers around the ionized gas and pillars in different evolutionary states. Finally, their velocity structure was investigated using CO data, and all observational signatures were tested against some distinct diagnostics established from simulations. The column density profiles have revealed the importance of compression at the edge of the ionized gas. The velocity properties of the structures, i.e. pillars and globules, are very close to what we predict from the numerical simulations. We have identified a good candidate of a nascent pillar in the Rosette molecular cloud that presents the velocity pattern of the shell collapsing on itself, induced by a high local curvature. Globules have a bulk velocity dispersion that indicates the importance of the initial turbulence in their formation, as proposed from numerical simulations. Altogether, this study re-enforces the picture of pillar formation by shell collapse and globule formation by the ionization of highly turbulent clouds.

We present the DIScrete PERsistent Structures Extractor (DisPerSE), an open source software for the automatic and robust identification of structures in 2D and 3D noisy data sets. The software is designed to identify all sorts of topological structures, such as voids, peaks, sources, walls and filaments through segmentation, with a special emphasis put on the later ones. Based on discrete Morse theory, DisPerSE is able to deal directly with noisy datasets using the concept of persistence (a measure of the robustness of topological features) and can be applied indifferently to various sorts of data-sets defined over a possibly bounded manifold : 2D and 3D images, structured and unstructured grids, discrete point samples via the delaunay tesselation, Healpix tesselations of the sphere, ... Although it was initially developed with cosmology in mind, various I/O formats have been implemented and the current version is quite versatile. It should therefore be useful for any application where a robust structure identification is required as well as for studying the topology of sampled functions (e.g. computing persistent Betti numbers). DisPerSE can be downloaded directly from the website http://www2.iap.fr/users/sousbie/ and a thorough online documentation is also available at the same address.

We investigate the alignment of the spin of dark matter halos relative (i) to the surrounding large-scale filamentary structure, and (ii) to the tidal tensor eigenvectors using the Horizon 4pi dark matter simulation which resolves over 43 million dark matter halos at redshift zero. We detect a clear mass transition: the spin of dark matter halos above a critical mass tends to be perpendicular to the closest filament, and aligned with the intermediate axis of the tidal tensor, whereas the spin of low-mass halos is more likely to be aligned with the closest filament. Furthermore, this critical mass of 5 10^12 is redshift-dependent and scales as (1+z)^-2.5. We propose an interpretation of this signal in terms of large-scale cosmic flows. In this picture, most low-mass halos are formed through the winding of flows embedded in misaligned walls; hence they acquire a spin parallel to the axis of the resulting filaments forming at the intersection of these walls. On the other hand, more massive halos are typically the products of later mergers along such filaments, and thus they acquire a spin perpendicular to this direction when their orbital angular momentum is converted into spin. We show that this scenario is consistent with both the measured excess probabilities of alignment w.r.t. the eigen-directions of the tidal tensor, and halo merger histories. On a more qualitative level, it also seems compatible with 3D visualization of the structure of the cosmic web as traced by "smoothed" dark matter simulations or gas tracer particles. Finally, it provides extra support to the disc forming paradigm presented by Pichon et al (2011) as it extends it by characterizing the geometry of secondary infall at high redshift.

In order to characterise the cloud structures responsible for the formation of high-mass stars, we present Herschel observations of the DR21 environment. Maps of the column density and dust temperature unveil the structure of the DR21 ridge and several connected filaments. The ridge has column densities larger than 1e23/cm^2 over a region of 2.3 pc^2. It shows substructured column density profiles and branching into two major filaments in the north. The masses in the studied filaments range between 130 and 1400 Msun whereas the mass in the ridge is 15000 Msun. The accretion of these filaments onto the DR21 ridge, suggested by a previous molecular line study, could provide a continuous mass inflow to the ridge. In contrast to the striations seen in e.g., the Taurus region, these filaments are gravitationally unstable and form cores and protostars. These cores formed in the filaments potentially fall into the ridge. Both inflow and collisions of cores could be important to drive the observed high-mass star formation. The evolutionary gradient of star formation running from DR21 in the south to the northern branching is traced by decreasing dust temperature. This evolution and the ridge structure can be explained by two main filamentary components of the ridge that merged first in the south.

We present Herschel images from the HOBYS key program of the Eagle Nebula (M16) in the far-infrared and sub-millimetre, using the PACS and SPIRE cameras at 70{\mu}m, 160{\mu}m, 250{\mu}m, 350{\mu}m, 500{\mu}m. M16, home to the Pillars of Creation, is largely under the influence of the nearby NGC6611 high-mass star cluster. The Herschel images reveal a clear dust temperature gradient running away from the centre of the cavity carved by the OB cluster. We investigate the heating effect of NGC6611 on the entire M16 star-forming complex seen by Herschel including the diffuse cloud environment and the dense filamentary structures identified in this region. In addition, we interpret the three-dimensional geometry of M16 with respect to the nebula, its surrounding environment, and the NGC6611 cavity. The dust temperature and column density maps reveal a prominent eastern filament running north-south and away from the high-mass star-forming central region and the NGC6611 cluster, as well as a northern filament which extends around and away from the cluster. The dust temperature in each of these filaments decreases with increasing distance from the NGC6611 cluster, indicating a heating penetration depth of \sim 10 pc in each direction in 3 - 6 \times 10^{22} cm-2 column density filaments. We show that in high-mass star-forming regions OB clusters impact the temperature of future star-forming sites, modifying the initial conditions for collapse and effecting the evolutionary criteria of protostars developed from spectral energy distributions. Possible scenarios for the origin of the morphology seen in this region are discussed, including a western equivalent to the eastern filament, which was destroyed by the creation of the OB cluster and its subsequent winds and radiation.

For many years feedback processes generated by OB-stars in molecular clouds, including expanding ionization fronts, stellar winds, or UV-radiation, have been proposed to trigger subsequent star formation. However, hydrodynamic models including radiation and gravity show that UV-illumination has little or no impact on the global dynamical evolution of the cloud. The Rosette molecular cloud, irradiated by the NGC2244 cluster, is a template region for triggered star-formation, and we investigated its spatial and density structure by applying a curvelet analysis, a filament-tracing algorithm (DisPerSE), and probability density functions (PDFs) on Herschel column density maps, obtained within the HOBYS key program. The analysis reveals not only the filamentary structure of the cloud but also that all known infrared clusters except one lie at junctions of filaments, as predicted by turbulence simulations. The PDFs of sub-regions in the cloud show systematic differences. The two UV-exposed regions have a double-peaked PDF we interprete as caused by shock compression. The deviations of the PDF from the log-normal shape typically associated with low- and high-mass star-forming regions at Av~3-4m and 8-10m, respectively, are found here within the very same cloud. This shows that there is no fundamental difference in the density structure of low- and high-mass star-forming regions. We conclude that star-formation in Rosette - and probably in high-mass star-forming clouds in general - is not globally triggered by the impact of UV-radiation. Moreover, star formation takes place in filaments that arose from the primordial turbulent structure built up during the formation of the cloud. Clusters form at filament mergers, but star formation can be locally induced in the direct interaction zone between an expanding HII--region and the molecular cloud.

A growing body of evidence indicates that the formation of filaments in interstellar clouds is a key component of the star formation process. In this paper, we present new Herschel PACS and SPIRE observations of the B59 and Stem regions in the Pipe Nebula complex, revealing a rich, organized network of filaments. The asymmetric column density profiles observed for several filaments, along with the bow-like edge of B59, indicates that the Pipe Nebula is being compressed from its western side, most likely by the winds from the nearby Sco OB2 association. We suggest that this compressive flow has contributed to the formation of some of the observed filamentary structures. In B59, the only region of the entire Pipe complex showing star formation activity, the same compressive flow has likely enhanced the initial column density of the clump, allowing it to become globally gravitationally unstable. Although more speculative, we propose that gravity has also been responsible for shaping the converging filamentary pattern observed in B59. While the question of the relative impact of large-scale compression and gravity remains open in B59, large-scale compression appears to be a plausible mechanism for the initial formation of filamentary structures in the rest of the complex

We present the first Herschel PACS and SPIRE results of the Vela C molecular complex in the far-infrared and submillimetre regimes at 70, 160, 250, 350, and 500 um, spanning the peak of emission of cold prestellar or protostellar cores. Column density and multi-resolution analysis (MRA) differentiates the Vela C complex into five distinct sub-regions. Each sub-region displays differences in their column density and temperature probability distribution functions (PDFs), in particular, the PDFs of the `Centre-Ridge' and `South-Nest' sub-regions appear in stark contrast to each other. The Centre-Ridge displays a bimodal temperature PDF representative of hot gas surrounding the HII region RCW 36 and the cold neighbouring filaments, whilst the South-Nest is dominated by cold filamentary structure. The column density PDF of the Centre-Ridge is flatter than the South-Nest, with a high column density tail, consistent with formation through large-scale flows, and regulation by self-gravity. At small to intermediate scales MRA indicates the Centre-Ridge to be twice as concentrated as the South-Nest, whilst on larger scales, a greater portion of the gas in the South-Nest is dominated by turbulence than in the Centre-Ridge. In Vela C, high-mass stars appear to be preferentially forming in ridges, i.e., dominant high column density filaments.

We provide a first look at the results of the Herschel Gould Belt survey toward the IC5146 molecular cloud and present a preliminary analysis of the filamentary structure in this region. The column density map, derived from our 70-500 micron Herschel data, reveals a complex network of filaments, and confirms that these filaments are the main birth sites of prestellar cores. We analyze the column density profiles of 27 filaments and show that the underlying radial density profiles fall off as r^{-1.5} to r^{-2.5} at large radii. Our main result is that the filaments seem to be characterized by a narrow distribution of widths having a median value of 0.10 +- 0.03 pc, which is in stark contrast to a much broader distribution of central Jeans lengths. This characteristic width of ~0.1 pc corresponds to within a factor of ~2 to the sonic scale below which interstellar turbulence becomes subsonic in diffuse gas, supporting the argument that the filaments may form as a result of the dissipation of large-scale turbulence.

We report on a new merging group of galaxies, Suzaku J1552+2739 at z ~ 0.08 revealed by a SUZAKU observation. The group was found by observing a junction of galaxy filaments optically identified in the Sloan Digital Sky Survey spectroscopic data. Suzaku J1552+2739 exhibits an irregular morphology and presents several peaks in its X-ray image. A bright elliptical galaxy, observable in the central peak, allows the localization of the group at z = 0.083. We found a significant hot spot visible in the X-ray hardness map, close to the second peak. The spectroscopic temperature is T = 1.6+0.4-0.1 keV within R500 = 0.6 Mpc and T = 3 - 5 keV in the hot spot. We interpret those results as Suzaku J1552+2739 being located in the center of a major merging process. The observation of a galaxy group showing multiple X-ray peaks and a hot spot at the same time is rare and we believe in particular that the study of Suzaku J1552+2739 potentially presents a significant interest to better understand the dynamical and thermal evolution of the intragroup and intracluster medium, as well as its relation with surrounding environment.

The recently introduced discrete persistent structure extractor (DisPerSE, Soubie 2010, paper I) is implemented on realistic 3D cosmological simulations and observed redshift catalogues (SDSS); it is found that DisPerSE traces equally well the observed filaments, walls, and voids in both cases. In either setting, filaments are shown to connect onto halos, outskirt walls, which circumvent voids. Indeed this algorithm operates directly on the particles without assuming anything about the distribution, and yields a natural (topologically motivated) self-consistent criterion for selecting the significance level of the identified structures. It is shown that this extraction is possible even for very sparsely sampled point processes, as a function of the persistence ratio. Hence astrophysicists should be in a position to trace and measure precisely the filaments, walls and voids from such samples and assess the confidence of the post-processed sets as a function of this threshold, which can be expressed relative to the expected amplitude of shot noise. In a cosmic framework, this criterion is comparable to friend of friend for the identifications of peaks, while it also identifies the connected filaments and walls, and quantitatively recovers the full set of topological invariants (Betti numbers) {\sl directly from the particles} as a function of the persistence threshold. This criterion is found to be sufficient even if one particle out of two is noise, when the persistence ratio is set to 3-sigma or more. The algorithm is also implemented on the SDSS catalogue and used to locat interesting configurations of the filamentary structure. In this context we carried the identification of an ``optically faint'' cluster at the intersection of filaments through the recent observation of its X-ray counterpart by SUZAKU. The corresponding filament catalogue will be made available online.

We present DisPerSE, a novel approach to the coherent multi-scale identification of all types of astrophysical structures, and in particular the filaments, in the large scale distribution of matter in the Universe. This method and corresponding piece of software allows a genuinely scale free and parameter free identification of the voids, walls, filaments, clusters and their configuration within the cosmic web, directly from the discrete distribution of particles in N-body simulations or galaxies in sparse observational catalogues. To achieve that goal, the method works directly over the Delaunay tessellation of the discrete sample and uses the DTFE density computed at each tracer particle; no further sampling, smoothing or processing of the density field is required. The idea is based on recent advances in distinct sub-domains of computational topology, which allows a rigorous application of topological principles to astrophysical data sets, taking into account uncertainties and Poisson noise. Practically, the user can define a given persistence level in terms of robustness with respect to noise (defined as a "number of sigmas") and the algorithm returns the structures with the corresponding significance as sets of critical points, lines, surfaces and volumes corresponding to the clusters, filaments, walls and voids; filaments, connected at cluster nodes, crawling along the edges of walls bounding the voids. The method is also interesting as it allows for a robust quantification of the topological properties of a discrete distribution in terms of Betti numbers or Euler characteristics, without having to resort to smoothing or having to define a particular scale. In this paper, we introduce the necessary mathematical background and describe the method and implementation, while we address the application to 3D simulated and observed data sets to the companion paper.

The correlation between the large-scale distribution of galaxies and their spectroscopic properties at z=1.5 is investigated using the Horizon MareNostrum cosmological run. We have extracted a large sample of 10^5 galaxies from this large hydrodynamical simulation featuring standard galaxy formation physics. Spectral synthesis is applied to these single stellar populations to generate spectra and colours for all galaxies. We use the skeleton as a tracer of the cosmic web and study how our galaxy catalogue depends on the distance to the skeleton. We show that galaxies closer to the skeleton tend to be redder, but that the effect is mostly due to the proximity of large haloes at the nodes of the skeleton, rather than the filaments themselves. This effects translate into a bimodality in the colour distribution of our sample. The origin of this bimodality is investigated and seems to follow from the ram pressure stripping of satellite galaxies within the more massive clusters of the simulation. The virtual catalogues (spectroscopical properties of the MareNostrum galaxies at various redshifts) are available online at http://www.iap.fr/users/pichon/MareNostrum/catalogues

We report on two quantitative, morphological estimators of the filamentary structure of the Cosmic Web, the so-called global and local skeletons. The first, based on a global study of the matter density gradient flow, allows us to study the connectivity between a density peak and its surroundings, with direct relevance to the anisotropic accretion via cold flows on galactic halos. From the second, based on a local constraint equation involving the derivatives of the field, we can derive predictions for powerful statistics, such as the differential length and the relative saddle to extrema counts of the Cosmic web as a function of density threshold (with application to percolation of structures and connectivity), as well as a theoretical framework to study their cosmic evolution through the onset of gravity-induced non-linearities.

Within the context of upcoming full-sky lensing surveys, the edge-preserving non- linear algorithm Aski is presented. Using the framework of Maximum A Posteriori inversion, it aims at recovering the full-sky convergence map from surveys with masks. It proceeds in two steps: CCD images of crowded galactic fields are deblurred using automated edge-preserving deconvolution; once the reduced shear is estimated, the convergence map is also inverted via an edge- preserving method. For the deblurring, it is found that when the observed field is crowded, this gain can be quite significant for realistic ground-based surveys when both positivity and edge-preserving penalties are imposed during the iterative deconvolution. For the convergence inversion, the quality of the reconstruction is investigated on noisy maps derived from the horizon N-body simulation, with and without Galactic cuts, and quantified using one-point statistics, power spectra, cluster counts, peak patches and the skeleton. It is found that the reconstruction is able to interpolate and extrapolate within the Galactic cuts/non-uniform noise; its sharpness-preserving penalization avoids strong biasing near the clusters of the map; it reconstructs well the shape of the PDF as traced by its skewness and kurtosis; the geometry and topology of the reconstructed map is close to the initial map as traced by the peak patch distribution and the skeleton's differential length; the two-points statistics of the recovered map is consistent with the corresponding smoothed version of the initial map; the distribution of point sources is also consistent with the corresponding smoothing, with a significant improvement when edge preserving prior is applied. The contamination of B-modes when realistic Galactic cuts are present is also investigated. Leakage mainly occurs on large scales.

A method to compute the full hierarchy of the critical subsets of a density field is presented. It is based on a watershed technique and uses a probability propagation scheme to improve the quality of the segmentation by circumventing the discreteness of the sampling. It can be applied within spaces of arbitrary dimensions and geometry. This recursive segmentation of space yields, for a $d$-dimensional space, a $d-1$ succession of $n$-dimensional subspaces that fully characterize the topology of the density field. The final 1D manifold of the hierarchy is the fully connected network of the primary critical lines of the field : the skeleton. It corresponds to the subset of lines linking maxima to saddle points, and provides a definition of the filaments that compose the cosmic web as a precise physical object, which makes it possible to compute any of its properties such as its length, curvature, connectivity etc... When the skeleton extraction is applied to initial conditions of cosmological N-body simulations and their present day non linear counterparts, it is shown that the time evolution of the cosmic web, as traced by the skeleton, is well accounted for by the Zel'dovich approximation. Comparing this skeleton to the initial skeleton undergoing the Zel'dovich mapping shows that two effects are competing during the formation of the cosmic web: a general dilation of the larger filaments that is captured by a simple deformation of the skeleton of the initial conditions on the one hand, and the shrinking, fusion and disappearance of the more numerous smaller filaments on the other hand. Other applications of the N dimensional skeleton and its peak patch hierarchy are discussed.

The local theory of the critical lines of 2D and 3D Gaussian fields that underline the cosmic structures is presented. In the context of cosmological matter distribution the subset of critical lines of the 3D density field serves to delineate the skeleton of the observed filamentary structure at large scales. A stiff approximation used to quantitatively describe the filamentary skeleton shows that the flux of the skeleton lines is related to the average Gaussian curvature of the 1D (2D) sections of the field, much in the same way as the density of the peaks. The distribution of the length of the critical lines with threshold is analyzed in detail, while the extended descriptors of the skeleton - its curvature and its singular points, are introduced and briefly described. Theoretical predictions are compared to measurements of the skeleton in realizations of Gaussian random fields in 2D and 3D. It is found that the stiff approximation predicts accurately the shape of the differential length, allows for analytical insight, and explicit closed form solutions. Finally, it provides a simple classification of the singular points of the critical lines: i) critical points; ii) bifurcation points; iii) slopping plateaux.

The onset of stochasticity is measured in $\Lambda$CDM cosmological simulations using a set of classical observables. It is quantified as the local derivative of the logarithm of the dispersion of a given observable (within a set of different simulations differing weakly through their initial realization), with respect to the cosmic growth factor. In an Eulerian framework, it is shown here that chaos appears at small scales, where dynamic is non-linear, while it vanishes at larger scales, allowing the computation of a critical transition scale corresponding to ~ 3.5 Mpc/h. This picture is confirmed by Lagrangian measurements which show that the distribution of substructures within clusters is partially sensitive to initial conditions, with a critical mass upper bound scaling roughly like the perturbation's amplitude to the power 0.15. The corresponding characteristic mass, $M_{\rm crit}=2 10^{13} M_{\odot}$, is roughly of the order of the critical mass of non linearities at z=1 and accounts for the decoupling induced by the dark energy triggered acceleration. The sensitivity to detailed initial conditions spills to some of the overall physical properties of the host halo (spin and velocity dispersion tensor orientation) while other "global" properties are quite robust and show no chaos (mass, spin parameter, connexity and center of mass position). This apparent discrepancy may reflect the fact that quantities which are integrals over particles rapidly average out details of difference in orbits, while the other observables are more sensitive to the detailed environment of forming halos and reflect the non-linear scale coupling characterizing the environments of halos.

We investigate how well the 3D density field of neutral hydrogen in the Intergalactic Medium (IGM) can be reconstructed using the Lyman-alpha absorptions observed along lines of sight to quasars separated by arcmin distances in projection on the sky. We use cosmological hydrodynamical simulations to compare the topologies of different fields: dark matter, gas and neutral hydrogen optical depth and to investigate how well the topology of the IGM can be recovered from the Wiener interpolation method implemented by Pichon et al. (2001). The global statistical and topological properties of the recovered field are analyzed quantitatively through the power-spectrum, the probability distribution function (PDF), the Euler characteristics, its associated critical point counts and the filling factor of underdense regions. The local geometrical properties of the field are analysed using the local skeleton by defining the concept of inter-skeleton distance. At scales larger than ~1.4 <d_LOS>, where <d_LOS> is the mean separation between lines of sight, the reconstruction accurately recovers the topological features of the large scale density distribution of the gas, in particular the filamentary structures. At scales larger than the intrinsic smoothing length of the inversion procedure, the power spectrum of the recovered HI density field matches well that of the original one and the low order moments of the PDF are well recovered as well as the shape of the Euler characteristic. The integral errors on the PDF and the critical point counts are indeed small, less than 20% for <d_LOS>~2.5 arcmin. The small deviations between the reconstruction and the exact solution mainly reflect departures from the log-normal behaviour that are ascribed to highly non-linear objects in overdense regions.

The skeleton formalism aims at extracting and quantifying the filamentary structure of the universe is generalized to 3D density fields; a numerical method for computating a local approximation of the skeleton is presented and validated here on Gaussian random fields. This method manages to trace well the filamentary structure in 3D fields such as given by numerical simulations of the dark matter distribution on large scales and is insensitive to monotonic biasing. Two of its characteristics, namely its length and differential length, are analyzed for Gaussian random fields. Its differential length per unit normalized density contrast scales like the PDF of the underlying density contrast times the total length times a quadratic Edgeworth correction involving the square of the spectral parameter. The total length scales like the inverse square smoothing length, with a scaling factor given by 0.21 (5.28+ n) where n is the power index of the underlying field. This dependency implies that the total length can be used to constrain the shape of the underlying power spectrum, hence the cosmology. Possible applications of the skeleton to galaxy formation and cosmology are discussed. As an illustration, the orientation of the spin of dark halos and the orientation of the flow near the skeleton is computed for dark matter simulations. The flow is laminar along the filaments, while spins of dark halos within 500 kpc of the skeleton are preferentially orthogonal to the direction of the flow at a level of 25%.

This paper presents MoLUSC, a new method for generating mock galaxy catalogs from a large scale ($\approx 1000^3$ Mpc$^3$) dark matter simulation, that requires only modest CPU time and memory allocation. The method uses a small-scale ($\approx 256^3$ Mpc$^3$) dark matter simulation on which the \galics semi-analytic code has been run in order to define the transformation from dark matter density to galaxy density transformation using a probabilistic treatment. MoLUSC is then applied to a large-scale dark matter simulation in order to produce a realistic distribution of galaxies and their associated spectra. This permits the fast generation of large-scale mock surveys using relatively low-resolution simulations. We describe various tests which have been conducted to validate the method, and demonstrate a first application to generate a mock Sloan Digital Sky Survey redshift survey.

The length of the three-dimensional filaments observed in the fourth public data-release of the SDSS is measured using the {\em local skeleton} method. It consists in defining the set of points where the gradient of the smoothed density field is extremal along its isocontours, with some additional constraints on local curvature to probe actual ridges in the galaxy distribution. A good fit to the mean filament length per unit volume, $\cal{L}$, in the SDSS survey is found to be ${\cal{L}}=(52500\pm6500) (L/{\rm Mpc})^{-1.75\pm0.06}\rm{Mpc}/(100 \rm{Mpc})^{3}$ for $8.2 \leq L \leq 16.4$ Mpc, where $L$ is the smoothing length in Mpc. This result, which deviates only slightly, as expected, from the trivial behavior ${\cal{L}} \propto L^{-2}$, is in excellent agreement with a $\Lambda$CDM cosmology, as long as the matter density parameter remains in the range $0.25 < \Omega_{\rm matter} < 0.4$ at one sigma confidence level, considering the universe is flat. These measurements, which are in fact dominated by linear dynamics, are not significantly sensitive to observational biases such as redshift distortion, edge effects, incompleteness, and biasing between the galaxy distribution and the dark matter distribution. Hence it is argued that the local skeleton is a rather promising and discriminating tool for the analysis of filamentary structures in three-dimensional galaxy surveys.

We present predictions for the abundance of Ly-$\alpha$ emitters in hierarchical structure formation models. We use the {\tt GALFORM} semi-analytical model to explore the impact on the predicted counts of varying assumptions about the escape fraction of Ly-$\alpha$ photons, the redshift at which the universe reionised and the cosmological density parameter. A model with a fixed escape fraction gives a remarkably good match to the observed counts over a wide redshift interval. We present predictions for the expected counts in a typical observation with the Multi Unit Spectroscopic Explorer instrument proposed for the Very Large Telescope.

In order to investigate the properties of large-scale structures of galaxies in the universe, we present an analysis of their spatial distribution at z$<$0.033. We used the LEDA extragalactic database containing over 1 million of galaxies covering the all-sky and the SDSS data included in the public release DR1, yielding to a sample of around 134,000 galaxies having a measured redshift in two survey areas representing 690 sq. degrees. The results of the study are 2D, 3D maps and magnitude number counts of galaxies, drawn from B-band samples. Movies and high resolution figures are available on http://www-obs.univ-lyon1.fr/cosmologie/helene.html