Source author record

Rien van de Weygaert

Rien van de Weygaert appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

41works
8topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

41 published item(s)

preprint2025arXiv

Chemodynamic evidence of pristine gas accretion in the void galaxy VGS 12

Accretion of metal-poor gas is expected to be an important channel of gas replenishment in galaxy evolution studies. However, observational evidence of this process is still relatively scarce. The unusual polar disk galaxy VGS 12 was found in the Void Galaxy Survey. It appears to be isolated and resides in the cosmological wall between two large voids. The suggested formation scenario for this peculiar system is accretion of metal-poor gas from the void interior. We present the data on the gas-phase chemical abundance of VGS 12 obtained with the Russian 6m telescope BTA. We complement our analysis with HI data obtained with VLA and the data on the kinematics of the ionized gas. VGS 12 appears to be a strong outlier from the "metallicity - luminosity" relation, with gas oxygen abundance ~0.7dex lower than expected for its luminosity. The nitrogen abundance, on the other hand, is higher than what is typically observed in galaxies with similar metallicity, but is consistent with the metallicity expected given its luminosity. Such behavior is what is expected in the case of metal-poor gas accretion. The HI reveals clear morphological and kinematical asymmetry between the northern and southern parts of the disk, which are likely related to its unsettled state due to the recent accretion event. The kinematics of the ionized gas seen in Halpha reveal prolate rotation and follow closely the rotation of the HI disk, so we suggest this is accreted HI gas ionized by the stars in the central region of the galaxy. Together, our findings provide strong, multiwavelength evidence of ongoing cold gas accretion in a galaxy caught in the act of growing from the cosmic web. This is one of the very few individual galaxies where a convincing case can be made for such a process, and demonstrates the potential for cold accretion to contribute to galaxy growth even in the low-redshift universe.

preprint2023arXiv

Statistics of tidal and deformation eigenvalue fields in the primordial Gaussian matter distribution: the two-dimensional case

We study the statistical properties of the eigenvalues of the primordial tidal and deformation tensor for random Gaussian cosmic density fields. With the tidal and deformation tensors, Hessians of the gravitational and velocity potential, being Gaussian, the corresponding eigenvalue fields are distinctly non-Gaussian. Following the extension of the Doroshkevich formula for the joined distribution of eigenvalues to two-dimensional fields, we evaluate the two- and three-point correlation functions of the eigenvalue fields. In addition, we assess the number densities of singular points of the eigenvalue fields and find their corresponding two- and three-point correlation functions. The role of tidal forces and the resulting mass element deformation in shaping the prominent anisotropic wall-like and filamentary components of the cosmic web has since long been recognized based on the Zel'dovich approximation. Less well-known is that the weblike spatial pattern is already recognizable in the primordial tidal and deformation eigenvalue field, even while the corresponding Gaussian density and the potential field appear merely as a spatially incoherent and unstructured random field. Furthermore, against the background of a full phase-space assessment of structure formation in the Universe, the caustic skeleton theory entails a fully analytical framework for the nonlinear evolution of the cosmic web. It describes the folding of the dark matter sheet and the emerging caustic singularities, fully specified by the deformation eigenvalues and eigenvectors. Finally, tidal tensor eigenvalues are of central importance, and understanding their distribution is critical in predicting the resulting rotation and orientation. The current study applies to two-dimensional Gaussian random fields and will be generalized to a three-dimensional analysis in an upcoming study.

preprint2021arXiv

Cosmic Ballet III: halo spin evolution in the cosmic web

We explore the evolution of halo spins in the cosmic web using a very large sample of dark matter haloes in the $Λ$CDM Planck-Millennium N-body simulation. We use the NEXUS+ multiscale formalism to identify the hierarchy of filaments and sheets of the cosmic web at several redshifts. We find that at all times the magnitude of halo spins correlates with the web environment, being largest in filaments, and, for the first time, we show that it also correlates with filament thickness as well as the angle between spin-orientation and the spine of the host filament. For example, massive haloes in thick filaments spin faster than their counterparts in thin filaments, while for low-mass haloes the reverse is true. We also have studied the evolution of alignment between halo spin orientations and the preferential axes of filaments and sheets. The alignment varies with halo mass, with the spins of low-mass haloes being predominantly along the filament spine, while those of high-mass haloes being predominantly perpendicular to the filament spine. On average, for all halo masses, halo spins become more perpendicular to the filament spine at later times. At all redshifts, the spin alignment shows a considerable variation with filament thickness, with the halo mass corresponding to the transition from parallel to perpendicular alignment varying by more than one order of magnitude. The environmental dependence of halo spin magnitude shows little evolution for $z\leq2$ and is likely a consequence of the correlations in the initial conditions or high redshift effects

preprint2021arXiv

Deviations from tidal torque theory: evolution of the halo spin-filament alignment

The alignment between halo spins and the cosmic web is still poorly understood despite being a widely studied topic. Here, we study this alignment within the context of tidal torque theory (TTT) and deviations from it. To this end, we analyze the evolution of the shape and spin direction of proto-haloes, i.e. of all the volume elements associated to a $z=0$ halo, with respect to the present-day filaments. We find that the major axis of proto-haloes undergoes a major change, from being strongly perpendicular to the filament spine in the initial conditions, to being preferentially aligned at the present time. In comparison, the spin orientation shows only a mild evolution: it starts slightly parallel to the filament spine, but the subsequent evolution, up to $z{\sim}1$, gradually changes its orientation to preferentially perpendicular. In order to analyze these signals in the TTT framework, we split the haloes according to their net spin growth with respect to the median TTT expectation, finding a clear correlation with the spin--filament alignment. At the present time, haloes whose spin grew the most are the ones most perpendicular to the filament spine, while haloes whose spin grew below the median TTT expectation are typically more aligned. The dependence of spin directions on net spin growth is already present in the initial conditions, and gets further modified by late-time, $z<2$, evolution. Also, spin directions mildly deviate from the TTT predictions even at high redshift, indicating the need for extensions to the model.

preprint2016arXiv

The structural elements of the cosmic web

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

preprint2016arXiv

The Zeldovich & Adhesion approximations, and applications to the local universe

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

preprint2016arXiv

Voids and the Cosmic Web: cosmic depressions & spatial complexity

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

preprint2015arXiv

Understanding the cosmic web

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

preprint2014arXiv

Evolution of the cosmic web

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

preprint2014arXiv

Milky Way mass constraints from the Galactic satellite gap

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

preprint2014arXiv

Statistics of Caustics in Large-Scale Structure Formation

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

preprint2014arXiv

Sub-Megaparsec Individual Photometric Redshift Estimation from Cosmic Web Constraints

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.

preprint2014arXiv

Subhalo statistics of galactic halos: beyond the resolution limit

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

preprint2014arXiv

The cosmic web in CosmoGrid void regions

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

preprint2014arXiv

The Void Galaxy Survey: Galaxy Evolution and Gas Accretion in Voids

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.

preprint2013arXiv

Alpha, Betti and the Megaparsec Universe: on the Topology of the Cosmic Web

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

preprint2013arXiv

Assembly of filamentary void galaxy configurations

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

preprint2013arXiv

Betti numbers of Gaussian fields

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

preprint2013arXiv

Hierarchical clustering and the BAO signature

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

preprint2013arXiv

The Zeldovich approximation: key to understanding Cosmic Web complexity

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

preprint2012arXiv

Hot and cold spots counts as probes of non-Gaussianity in the CMB

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

preprint2012arXiv

Nexus of the Cosmic Web

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

preprint2012arXiv

NEXUS: Tracing the Cosmic Web Connection

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

preprint2012arXiv

Skewness as a probe of Baryon Acoustic Oscillations

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

preprint2012arXiv

The darkness that shaped the void: dark energy and cosmic voids

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

preprint2012arXiv

The Sticky Geometry of the Cosmic Web

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

preprint2011arXiv

Probing Dark Energy with Alpha Shapes and Betti Numbers

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

preprint2011arXiv

Structural Analysis of the SDSS Cosmic Web I.Nonlinear Density Field Reconstructions

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.

preprint2010arXiv

Alpha Shape Topology of the Cosmic Web

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

preprint2010arXiv

Fossil evidence for spin alignment of SDSS galaxies in filaments

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

preprint2010arXiv

Hierarchy of N-point functions in the LCDM and ReBEL cosmologies

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

preprint2010arXiv

The Spine of the Cosmic Web

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.

preprint2009arXiv

Cosmic Voids: structure, dynamics and galaxies

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

preprint2009arXiv

Cosmology and Cluster Halo Scaling Relations

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

preprint2009arXiv

Future Evolution of Bound Superclusters in an Accelerating Universe

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

preprint2008arXiv

Cosmic Order out of Primordial Chaos: a tribute to Nikos Voglis

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

preprint2007arXiv

A Cosmic Watershed: the WVF Void Detection Technique

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

preprint1995arXiv

Constraining Peaks in Gaussian Primordial Density Fields: an Application of the Hoffman-Ribak Method

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

preprint1995arXiv

The Two--Point Correlation Function and the Size of Voids

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