-
Identifying Dark Matter Haloes by the Caustic Boundary
Authors:
Segei F. Shandarin
Abstract:
Dark matter density is formally infinite at the location of caustic surfaces, where dark matter sheet folds in phase space. The caustics separate multi-stream regions with different number of streams. Volume elements change the parity by turning inside out when passing through the caustic stage. Being measure-zero structures, identification of caustics via matter density fields is usually restrict…
▽ More
Dark matter density is formally infinite at the location of caustic surfaces, where dark matter sheet folds in phase space. The caustics separate multi-stream regions with different number of streams. Volume elements change the parity by turning inside out when passing through the caustic stage. Being measure-zero structures, identification of caustics via matter density fields is usually restricted to fine-grained simulations. Instead a generic purely geometric algorithm can be employed to identify caustics directly by using triangulation of Lagrangian sub-manifold x(q, t) where x and q are Eulerian and Lagrangian coordinates obtained in N-body simulations. The caustic surfaces are approximated by a set of triangles with vertices being the particles in the simulation. It is demonstrated that finding a dark matter halo is quite feasible by building its owtermost convex caustic. Neither more specific assumptions about the geometry of the boundary nor ad hoc parameters are needed. The halo boundary in our idealized but undoubtably generic simulation is neither spherical nor ellipsoidal but rather remarkably asymmetrical. The analysis of the kinetic and potential energies of individual particles and the halo as a whole along with an examination of the two-dimensional phase space has shown that the halo is gravitationally bound.
△ Less
Submitted 27 May, 2020;
originally announced May 2020.
-
The Caustic Design of the Dark Matter Web
Authors:
Sergei F. Shandarin,
Nesar S. Ramachandra
Abstract:
Matter density is formally infinite at the location of caustic surfaces, where dark matter sheet folds in phase-space. The caustics separate regions with different number of streams and the volume elements change the parity by turning inside out when passing through the caustic stage. Being measure-zero structures, identification of caustics via matter density fields is usually restricted to fine-…
▽ More
Matter density is formally infinite at the location of caustic surfaces, where dark matter sheet folds in phase-space. The caustics separate regions with different number of streams and the volume elements change the parity by turning inside out when passing through the caustic stage. Being measure-zero structures, identification of caustics via matter density fields is usually restricted to fine-grained simulations. Here we employ a generic algorithm to identify caustics directly using the triangulation of the Lagrangian sub-manifold x(q,t) obtained in N-body simulations. In our approach the caustic surfaces are approximated by a set of triangles whose vertices are particles of the simulation. The major obstacle we encountered was insufficient sampling of small scale perturbations. We overcame it by a brute force approach. We continued to raise the scale of the cutoff in the initial power spectrum until obtained the reliable resolution of the caustics shells up to seven layers. Although quite modest, our result is the first reliable direct construction of caustic surfaces in N-body simulation. It reveals a number of unexpected geometrical features. In particular shapes of some of them are contrastingly different from the known shapes of the caustics formed in the Zeldovich approximation.
△ Less
Submitted 13 June, 2019;
originally announced June 2019.
-
Dark matter haloes: a multistream view
Authors:
Nesar S. Ramachandra,
Sergei F. Shandarin
Abstract:
Mysterious dark matter constitutes about 85% of all mass in the Universe. Clustering of dark matter plays the dominant role in the formation of all observed structures on scales from a fraction to a few hundreds of Mega-parsecs. Galaxies play a role of lights illuminating these structures so they can be observed. The observations in the last several decades have unveiled opulent geometry of these…
▽ More
Mysterious dark matter constitutes about 85% of all mass in the Universe. Clustering of dark matter plays the dominant role in the formation of all observed structures on scales from a fraction to a few hundreds of Mega-parsecs. Galaxies play a role of lights illuminating these structures so they can be observed. The observations in the last several decades have unveiled opulent geometry of these structures currently known as the cosmic web. Haloes are the highest concentrations of dark matter and host luminous galaxies. Currently the most accurate modeling of dark matter haloes is achieved in cosmological N-body simulations. Identifying the haloes from the distribution of particles in N-body simulations is one of the problems attracting both considerable interest and efforts. We propose a novel framework for detecting potential dark matter haloes using the field unique for dark matter -- multistream field. The multistream field emerges at the nonlinear stage of the growth of perturbations because the dark matter is collisionless. Counting the number of velocity streams in gravitational collapses supplements our knowledge of spatial clustering. We assume that the virialized haloes have convex boundaries. Closed and convex regions of the multistream field are hence isolated by imposing a positivity condition on all three eigenvalues of the Hessian estimated on the smoothed multistream field. In a single-scale analysis of high multistream field resolution and low softening length, the halo substructures with local multistream maxima are isolated as individual halo sites.
△ Less
Submitted 9 June, 2017;
originally announced June 2017.
-
Tracing the cosmic web
Authors:
Noam I Libeskind,
Rien van de Weygaert,
Marius Cautun,
Bridget Falck,
Elmo Tempel,
Tom Abel,
Mehmet Alpaslan,
Miguel A. Aragoon-Calvo,
Jaime E. Forero-Romero,
Roberto Gonzalez,
Stefan Gottloober,
Oliver Hahn,
Wojciech A. Hellwing,
Yehuda Hoffman,
Bernard J. T. Jones,
Francisco Kitaura,
Alexander Knebe,
Serena Manti,
Mark Neyrinck,
Sebastiaan E. Nuza,
Nelson Padilla,
Erwin Platen,
Nesar Ramachandra,
Aaron Robotham,
Enn Saar
, et al. (5 additional authors not shown)
Abstract:
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 effe…
▽ More
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.
△ Less
Submitted 8 May, 2017;
originally announced May 2017.
-
The Zeldovich & Adhesion approximations, and applications to the local universe
Authors:
Johan Hidding,
Rien van de Weygaert,
Sergei Shandarin
Abstract:
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 deriva…
▽ More
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.
△ Less
Submitted 3 November, 2016;
originally announced November 2016.
-
The features of the Cosmic Web unveiled by the flip-flop field
Authors:
Sergei F. Shandarin,
Mikhail V. Medvedev
Abstract:
Understanding of the observed structure in the universe can be reached only in the theoretical framework of dark matter. N-body simulations are indispensable for the analysis of the formation and evolution of the dark matter web. Two primary fields - density and velocity fields - are used in most of studies. However dark matter provides two additional fields which are unique for collisionless medi…
▽ More
Understanding of the observed structure in the universe can be reached only in the theoretical framework of dark matter. N-body simulations are indispensable for the analysis of the formation and evolution of the dark matter web. Two primary fields - density and velocity fields - are used in most of studies. However dark matter provides two additional fields which are unique for collisionless media only. These are the multi- stream field in Eulerian space and flip-flop field in Lagrangian space. The flip-flop field represents the number of sign reversals of an elementary volume of each collisionless fluid element. This field can be estimated by counting the sign reversals of the Jacobian at each particle at every time step of the simulation. The Jacobian is evaluated by numerical differentiation of the Lagrangian submanifold, i.e., the three-dimensional dark matter sheet in the six-dimensional space formed by three Lagrangian and three Eulerian coordinates. We present the results of the statistical study of the evolution of the flip-flop field from z = 50 to the present time z = 0. A number of statistical characteristics show that the pattern of the flip-flop field remains remarkably stable from z = 30 to the present time. As a result the flip-flop field evaluated at z = 0 stores a wealth of information about the dynamical history of the dark matter web. In particular one of the most intriguing properties of the flip-flop is a unique capability to preserve the information about the merging history of dark matter haloes.
△ Less
Submitted 27 September, 2016;
originally announced September 2016.
-
Topology and geometry of the dark matter web: a multistream view
Authors:
Nesar S. Ramachandra,
Sergei F. Shandarin
Abstract:
Topological connections in the single-streaming voids and multistreaming filaments and walls reveal a cosmic web structure different from traditional mass density fields. A single void structure not only percolates the multistream field in all the directions, but also occupies over 99 per cent of all the single-streaming regions. Sub-grid analyses on scales smaller than simulation resolution revea…
▽ More
Topological connections in the single-streaming voids and multistreaming filaments and walls reveal a cosmic web structure different from traditional mass density fields. A single void structure not only percolates the multistream field in all the directions, but also occupies over 99 per cent of all the single-streaming regions. Sub-grid analyses on scales smaller than simulation resolution reveal tiny pockets of voids that are isolated by membranes of the structure. For the multistreaming excursion sets, the percolating structure is significantly thinner than the filaments in over-density excursion approach.
Hessian eigenvalues of the multistream field are used as local geometrical indicators of dark matter structures. Single-streaming regions have most of the zero eigenvalues. Parameter-free conditions on the eigenvalues in the multistream region may be used to delineate primitive geometries with concavities corresponding to filaments, walls and haloes.
△ Less
Submitted 4 May, 2017; v1 submitted 18 August, 2016;
originally announced August 2016.
-
Multi-stream portrait of the Cosmic web
Authors:
Nesar S. Ramachandra,
Sergei F. Shandarin
Abstract:
We report the results of the first study of the multi-stream environment of dark matter haloes in cosmological N-body simulations in the LCDM cosmology. The full dynamical state of dark matter can be described as a three-dimensional sub-manifold in six dimensional phase space - the dark matter sheet. In our study we use a Lagrangian sub-manifold x = x(q,t) (where x and q are co-moving Eulerian and…
▽ More
We report the results of the first study of the multi-stream environment of dark matter haloes in cosmological N-body simulations in the LCDM cosmology. The full dynamical state of dark matter can be described as a three-dimensional sub-manifold in six dimensional phase space - the dark matter sheet. In our study we use a Lagrangian sub-manifold x = x(q,t) (where x and q are co-moving Eulerian and Lagrangian coordinates respectively), which is dynamically equivalent to the dark matter sheet but is more convenient for numerical analysis. Our major results can be summarized as follows. At the resolution of the simulation i.e. without additional smoothing, the cosmic web represents a hierarchical structure: each halo is embedded in the filamentary framework of the web predominantly at the filament crossings, and each filament is embedded in the wall like fabric of the web at the wall crossings. Locally, each halo or sub-halo is a peak in the number of streams field. The number of streams in the neighbouring filaments is higher than in the neighbouring walls. The walls are regions where number of streams is equal to three or a few. Voids are uniquely defined by the local condition requiring to be a single-stream flow region. The shells of streams around haloes are quite thin and the closest void region is typically within one and a half FOF radius from the center of the halo.
△ Less
Submitted 19 August, 2016; v1 submitted 24 December, 2014;
originally announced December 2014.
-
Tracing the Cosmic Web substructure with Lagrangian submanifold
Authors:
Sergei F. Shandarin,
Mikhail V. Medvedev
Abstract:
A new computational paradigm for the analysis of substructure of the Cosmic Web in cosmological cold dark matter simulations is proposed. We introduce a new data-field --- the flip-flop field ---which carries wealth of information about the history and dynamics of the structure formation in the universe. The flip-flop field is an ordered data set in Lagrangian space representing the number of turn…
▽ More
A new computational paradigm for the analysis of substructure of the Cosmic Web in cosmological cold dark matter simulations is proposed. We introduce a new data-field --- the flip-flop field ---which carries wealth of information about the history and dynamics of the structure formation in the universe. The flip-flop field is an ordered data set in Lagrangian space representing the number of turns inside out sign reversals of an elementary volume of each collisionless fluid element represented by a computational particle in a N-body simulation. This field is computed using the Lagrangian submanifold, i.e. the three-dimensional dark matter sheet in the six-dimensional space formed by three Lagrangian and three Eulerian coordinates of the simulation particles. It is demonstrated that the very rich substructure of dark matter haloes and the void regions can be reliably and unambiguously recovered from the flip-flop field.
△ Less
Submitted 26 September, 2014;
originally announced September 2014.
-
The Zeldovich approximation: key to understanding Cosmic Web complexity
Authors:
Johan Hidding,
Sergei F. Shandarin,
Rien van de Weygaert
Abstract:
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 profo…
▽ More
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.
△ Less
Submitted 20 November, 2013;
originally announced November 2013.
-
The Phase-Space Structure of Cold Dark Matter in the Universe
Authors:
Sergei Shandarin
Abstract:
A novel method allowing to compute density, velocity and other fields in cosmological N--body simulations with unprecedentedly high spatial resolution is described. It is based on the tessellation of the three-dimensional manifold representing cold dark matter in six-dimensional phase space. The density, velocity and other fields are computed by projecting the tessellation on configuration space.…
▽ More
A novel method allowing to compute density, velocity and other fields in cosmological N--body simulations with unprecedentedly high spatial resolution is described. It is based on the tessellation of the three-dimensional manifold representing cold dark matter in six-dimensional phase space. The density, velocity and other fields are computed by projecting the tessellation on configuration space. The application of this technique to cosmological N--body simulations in LambdaCDM cosmology reveals a far more elaborate cosmic web then dot plots or self--adaptive SPH. In addition, this method allows to uniquely define physical voids and identify and study the caustic surfaces directly.
△ Less
Submitted 1 February, 2013;
originally announced February 2013.
-
Tessellating the cosmological dark-matter sheet: origami creases in the universe and ways to find them
Authors:
Mark C. Neyrinck,
Sergei F. Shandarin
Abstract:
Tessellations are valuable both conceptually and for analysis in the study of the large-scale structure of the universe. They provide a conceptual model for the 'cosmic web,' and are of great use to analyze cosmological data. Here we describe tessellations in another set of coordinates, of the initially flat sheet of dark matter that gravity folds up in rough analogy to origami. The folds that dev…
▽ More
Tessellations are valuable both conceptually and for analysis in the study of the large-scale structure of the universe. They provide a conceptual model for the 'cosmic web,' and are of great use to analyze cosmological data. Here we describe tessellations in another set of coordinates, of the initially flat sheet of dark matter that gravity folds up in rough analogy to origami. The folds that develop are called caustics, and they tessellate space into stream regions. Tessellations of the dark-matter sheet are also useful in simulation analysis, for instance for density measurement, and to identify structures where streams overlap.
△ Less
Submitted 18 July, 2012;
originally announced July 2012.
-
The Cosmic Web, Multi-Stream Flows, and Tessellations
Authors:
Sergei Shandarin,
Salman Habib,
Katrin Heitmann
Abstract:
Understanding the structure of the matter distribution in the Universe due to the action of the gravitational instability -- the cosmic web -- is complicated by lack of direct analytic access to the nonlinear domain of structure formation. Here, we suggest and apply a novel tessellation method designed for cold dark matter (CDM) N-body cosmological simulations. The method is based on the fact that…
▽ More
Understanding the structure of the matter distribution in the Universe due to the action of the gravitational instability -- the cosmic web -- is complicated by lack of direct analytic access to the nonlinear domain of structure formation. Here, we suggest and apply a novel tessellation method designed for cold dark matter (CDM) N-body cosmological simulations. The method is based on the fact that the initial CDM state can be described by a 3-D manifold (in a 6-D phase space) that remains continuous under evolution. Our technique uses the full phase space information and has no free parameters; it can be used to compute multi-stream and density fields, the main focus of this paper. Using a large-box LCDM simulation we carry out a variety of initial analyses with the technique. These include studying the correlation between multi-streaming and density, the identification of structures such as Zel'dovich pancakes and voids, and statistical measurements of quantities such as the volume fraction as a function of the number of streams -- where we find a remarkable scaling relation. Cosmological implications are briefly discussed.
△ Less
Submitted 9 November, 2011;
originally announced November 2011.
-
The multi-stream flows and the dynamics of the cosmic web
Authors:
Sergei F. Shandarin
Abstract:
A new numerical technique to identify the cosmic web is proposed. It is based on locating multi-stream flows, i.e. the places where the velocity field is multi-valued. The method is local in Eulerian space, simple and computaionally efficient. This technique uses the velocities of particles and thus takes into account the dynamical information. This is in contrast with the majority of standard met…
▽ More
A new numerical technique to identify the cosmic web is proposed. It is based on locating multi-stream flows, i.e. the places where the velocity field is multi-valued. The method is local in Eulerian space, simple and computaionally efficient. This technique uses the velocities of particles and thus takes into account the dynamical information. This is in contrast with the majority of standard methods that use the coordinates of particles only. Two quantities are computed in every mesh cell: the mean and variance of the velocity field. In the cells where the velocity is single-valued the variance must be equal to zero exactly, therefore the cells with non-zero variance are identified as multi-stream flows. The technique has been tested in a N-body simulation of the ŁCDM model. The preliminary analysis has shown that numerical noise does not pose a significant problem. The web identified by the new method has been compared with the web identified by the standard technique using only the particle coordinates. The comparison has shown overall similarity of two webs as expected, however they by no means are identical. For example, the isocontours of the corresponding fields have significantly different shapes and some density peaks of similar heights exhibit significant differences in the velocity variance and vice versa. This suggest that the density and velocity variance have a significant degree of independence. The shape of the two-dimensional pdf of density and velocity variance confirms this proposition. Thus, we conclude that the dynamical information probed by this technique introduces an additional dimension into analysis of the web.
△ Less
Submitted 10 November, 2010; v1 submitted 8 November, 2010;
originally announced November 2010.
-
Geometry of the Cosmic Web: Minkowski Functionals from the Delaunay Tessellation
Authors:
Miguel A. Aragon-Calvo,
Sergei F. Shandarin,
Alexander Szalay
Abstract:
We present a novel method for computing the Minkowski Functionals from isodensity surfaces extracted directly from the Delaunay tessellation of a point distribution. This is an important step forward compared to the previous cosmological studies when the isodensity surface was built in the field on a uniform cubic grid and therefore having a uniform spatial resolution. The density field representi…
▽ More
We present a novel method for computing the Minkowski Functionals from isodensity surfaces extracted directly from the Delaunay tessellation of a point distribution. This is an important step forward compared to the previous cosmological studies when the isodensity surface was built in the field on a uniform cubic grid and therefore having a uniform spatial resolution. The density field representing a particular interest in cosmology is the density of galaxies which is obtained from the highly nonuniform distribution of the galaxy positions. Therefore, the constraints caused by the spatially uniform grid put severe limitations on the studies of the geometry and shapes of the large-scale objects: superclusters and voids of galaxies. Our technique potentially is able to eliminate most of these limitations. The method is tested with some simple geometric models and an application to the density field from an N-body simulation is shown.
△ Less
Submitted 21 June, 2010;
originally announced June 2010.
-
Tessellating the Universe: the Zel'dovich and Adhesion tiling of space
Authors:
Sergei F. Shandarin
Abstract:
The adhesion approximation is a simple analytical model suggested for explanation of the major geometrical features of the observed structure in the galaxy distribution on scales from 1 to (a few)x100/h Mpc. It is based on Burgers' equation and therefore allows analysis in considerable detail. A particular version of the model that assumes the infinitesimal viscosity naturally results in irregul…
▽ More
The adhesion approximation is a simple analytical model suggested for explanation of the major geometrical features of the observed structure in the galaxy distribution on scales from 1 to (a few)x100/h Mpc. It is based on Burgers' equation and therefore allows analysis in considerable detail. A particular version of the model that assumes the infinitesimal viscosity naturally results in irregular tessellation of the universe. Generic elements of the tessellation: vertices, edges, faces and three-dimensional tiles can be associated with astronomical objects of different kinds: clusters, superclusters and voids of galaxies. Point-like vertices contain the most of the mass and one-dimensional edges (filaments) are the second massive elements. The least massive are the two-dimensional faces and tiles (voids). The evolution of the large-scale structure can be viewed as a continuous process that transports mass predominantly from the high- to low-dimensional elements of the tessellation. For instance, the mass from the cells flows into faces, edges and vertices, in turn the mass from faces flows into edges and vertices, etc. At the same time, the elements of the tessellation themselves are in continuous motion resulting in mergers of some vertices, growth of some tiles and shrinking and disappearance of the others as well as other metamorphoses.
△ Less
Submitted 22 December, 2009;
originally announced December 2009.
-
Origin of the Cosmic Network: Nature vs Nurture
Authors:
Sergei Shandarin,
Salman Habib,
Katrin Heitmann
Abstract:
The large-scale structure of the Universe, as traced by the distribution of galaxies, is now being revealed by large-volume cosmological surveys. The structure is characterized by galaxies distributed along filaments, the filaments connecting in turn to form a percolating network. Our objective here is to quantitatively specify the underlying mechanisms that drive the formation of the cosmic net…
▽ More
The large-scale structure of the Universe, as traced by the distribution of galaxies, is now being revealed by large-volume cosmological surveys. The structure is characterized by galaxies distributed along filaments, the filaments connecting in turn to form a percolating network. Our objective here is to quantitatively specify the underlying mechanisms that drive the formation of the cosmic network: By combining percolation-based analyses with N-body simulations of gravitational structure formation, we elucidate how the network has its origin in the properties of the initial density field (nature) and how its contrast is then amplified by the nonlinear mapping induced by the gravitational instability (nurture).
△ Less
Submitted 22 December, 2009;
originally announced December 2009.
-
The origin of 'Great Walls'
Authors:
Sergei Shandarin
Abstract:
A new semianalytical model that explains the formation and sizes of the 'great walls' - the largest structures observed in the universe is suggested. Although the basis of the model is the Zel'dovich approximation it has been used in a new way very different from the previous studies. Instead of traditional approach that evaluates the nonlinear density field it has been utilized for identificati…
▽ More
A new semianalytical model that explains the formation and sizes of the 'great walls' - the largest structures observed in the universe is suggested. Although the basis of the model is the Zel'dovich approximation it has been used in a new way very different from the previous studies. Instead of traditional approach that evaluates the nonlinear density field it has been utilized for identification of the regions in Lagrangian space that after the mapping to real or redshift space (depending on the kind of structure is studied) end up in the regions where shell-crossing occurs. The set of these regions in Lagrangian space form the progenitor of the structure and after the mapping it determines the pattern of the structure in real or redshift space. The particle trajectories have crossed in such regions and the mapping is no longer unique there. The progenitor after mapping makes only one stream in the multi-stream flow regions therefore it does not comprise all the mass. Nevertheless, it approximately retains the shape of the structure. The progenitor of the structure in redshift space depends on a few non-Gaussian fields and also it is strongly affected by two anisotropic fields that determine the pattern of great walls as well as their huge sizes. All the fields used in the mappings are derived from the linear potential smoothed at the current scale of nonlinearity which is $R_{nl} = 2.7$ {\hmpc} for the adopted parameters of the \lcdm universe normalized to $σ_8 = 0.8$. The model predicts the existence of walls with sizes significantly greater than 500 {\hmpc} that may be found in sufficiently large redshift surveys.
△ Less
Submitted 29 January, 2009; v1 submitted 27 December, 2008;
originally announced December 2008.
-
The Aspen--Amsterdam Void Finder Comparison Project
Authors:
Joerg M. Colberg,
Frazer Pearce,
Caroline Foster,
Erwin Platen,
Riccardo Brunino,
Mark Neyrinck,
Spyros Basilakos,
Anthony Fairall,
Hume Feldman,
Stefan Gottloeber,
Oliver Hahn,
Fiona Hoyle,
Volker Mueller,
Lorne Nelson,
Manolis Plionis,
Cristiano Porciani,
Sergei Shandarin,
Michael S. Vogeley,
Rien van de Weygaert
Abstract:
Despite a history that dates back at least a quarter of a century studies of voids in the large--scale structure of the Universe are bedevilled by a major problem: there exist a large number of quite different void--finding algorithms, a fact that has so far got in the way of groups comparing their results without worrying about whether such a comparison in fact makes sense. Because of the recen…
▽ More
Despite a history that dates back at least a quarter of a century studies of voids in the large--scale structure of the Universe are bedevilled by a major problem: there exist a large number of quite different void--finding algorithms, a fact that has so far got in the way of groups comparing their results without worrying about whether such a comparison in fact makes sense. Because of the recent increased interest in voids, both in very large galaxy surveys and in detailed simulations of cosmic structure formation, this situation is very unfortunate. We here present the first systematic comparison study of thirteen different void finders constructed using particles, haloes, and semi--analytical model galaxies extracted from a subvolume of the Millennium simulation. The study includes many groups that have studied voids over the past decade. We show their results and discuss their differences and agreements. As it turns out, the basic results of the various methods agree very well with each other in that they all locate a major void near the centre of our volume. Voids have very underdense centres, reaching below 10 percent of the mean cosmic density. In addition, those void finders that allow for void galaxies show that those galaxies follow similar trends. For example, the overdensity of void galaxies brighter than $m_B = -20 $ is found to be smaller than about -0.8 by all our void finding algorithms.
△ Less
Submitted 10 April, 2008; v1 submitted 6 March, 2008;
originally announced March 2008.
-
Dark matter caustics and the enhancement of self-annihilation flux
Authors:
Roya Mohayaee,
Sergei Shandarin,
Joseph Silk
Abstract:
Cold dark matter haloes are populated by caustics, which are yet to be resolved in N-body simulations or observed in the Universe. Secondary infall model provides a paradigm for the study of caustics in "typical" haloes assuming that they have had no major mergers and have grown only by smooth accretion. This is a particular characteristic of the smallest dark matter haloes of about 10^{-5} Mo,…
▽ More
Cold dark matter haloes are populated by caustics, which are yet to be resolved in N-body simulations or observed in the Universe. Secondary infall model provides a paradigm for the study of caustics in "typical" haloes assuming that they have had no major mergers and have grown only by smooth accretion. This is a particular characteristic of the smallest dark matter haloes of about 10^{-5} Mo, which although "atypical" contain no substructures and could have survived until now with no major mergers. Thus using this model as the first guidline, we evaluate the neutralino self-annihilation flux for these haloes. Our results show that caustics could leave a distinct sawteeth signature on the differential and cumulative fluxes coming from the outer regions of these haloes. The total annihilation signal from the regions away from the centre can be boosted by about forty percents.
△ Less
Submitted 16 April, 2007;
originally announced April 2007.
-
Morphology and Evolution of Simulated and Optical Clusters: A Comparative Analysis
Authors:
Nurur Rahman,
Janusz Krywult,
Patrick M. Motl,
Piotr Flin,
Sergei F. Shandarin
Abstract:
We have made a comparative study of morphological evolution in simulated DM halos and X-ray brightness distribution, and in optical clusters. Samples of simulated clusters include star formation with supernovae feedback, radiative cooling, and simulation in the adiabatic limit at three different redshifts, z = 0.0, 0.10, and 0.25. The optical sample contains 208 ACO clusters within redshift,…
▽ More
We have made a comparative study of morphological evolution in simulated DM halos and X-ray brightness distribution, and in optical clusters. Samples of simulated clusters include star formation with supernovae feedback, radiative cooling, and simulation in the adiabatic limit at three different redshifts, z = 0.0, 0.10, and 0.25. The optical sample contains 208 ACO clusters within redshift, $z \leq 0.25$. Cluster morphology, within 0.5 and 1.0 h$^{-1}$ Mpc from cluster center, is quantified by multiplicity and ellipticity.
We find that the distribution of the dark matter halos in the adiabatic simulation appear to be more elongated than the galaxy clusters. Radiative cooling brings halo shapes in excellent agreement with observed clusters, however, cooling along with feedback mechanism make the halos more flattened.
Our results indicate relatively stronger structural evolution and more clumpy distributions in observed clusters than in the structure of simulated clusters, and slower increase in simulated cluster shapes compared to those in the observed one.
Within $z \leq 0.1$, we notice an interesting agreement in the shapes of clusters obtained from the cooling simulations and observation. We also notice that the different samples of observed clusters differ significantly in morphological evolution with redshift. We highlight a few possibilities responsible for the discrepancy in morphological evolution of simulated and observed clusters.
△ Less
Submitted 29 December, 2005;
originally announced December 2005.
-
Caustics in Dark Matter Haloes
Authors:
Roya Mohayaee,
Stephane Colombi,
Bernard Fort,
Raphael Gavazzi,
Sergei Shandarin,
Jihad Touma
Abstract:
Caustics are formally singular structures, with infinite density, that form in collisionless media. The non-negligible velocity dispersion of dark matter particles renders their density finite. We evaluate the maximum density of the caustics within the framework of secondary infall model of formation of dark matter haloes. The result is then used to demonstrate that caustics can be probed by pro…
▽ More
Caustics are formally singular structures, with infinite density, that form in collisionless media. The non-negligible velocity dispersion of dark matter particles renders their density finite. We evaluate the maximum density of the caustics within the framework of secondary infall model of formation of dark matter haloes. The result is then used to demonstrate that caustics can be probed by properly stacking the weak-lensing signal of about 600 haloes. CFHTLS accompanied by X-ray observations and the space-based experiments like SNAP or DUNE can provide us with the required statistics and hence a way of distinguishing between the viable dark matter particle candidates. The extension of our results to more realistic models including the effects of mergers of haloes is briefly outlined.
△ Less
Submitted 24 October, 2005; v1 submitted 19 October, 2005;
originally announced October 2005.
-
Shapes and Sizes of Voids in the LCDM Universe: Excursion Set Approach
Authors:
Sergei Shandarin,
Hume A. Feldman,
Katrin Heitmann,
Salman Habib
Abstract:
We study the global distribution and morphology of dark matter voids in a LCDM universe using density fields generated by N-body simulations. Voids are defined as isolated regions of the low-density excursion set specified via density thresholds, the density thresholds being quantified by the corresponding filling factors, i.e., the fraction of the total volume in the excursion set. Our work enc…
▽ More
We study the global distribution and morphology of dark matter voids in a LCDM universe using density fields generated by N-body simulations. Voids are defined as isolated regions of the low-density excursion set specified via density thresholds, the density thresholds being quantified by the corresponding filling factors, i.e., the fraction of the total volume in the excursion set. Our work encompasses a systematic investigation of the void volume function, the volume fraction in voids, and the fitting of voids to corresponding ellipsoids and spheres. We emphasize the relevance of the percolation threshold to the void volume statistics of the density field both in the high redshift, Gaussian random field regime, as well as in the present epoch. By using measures such as the Inverse Porosity, we characterize the quality of ellipsoidal fits to voids, finding that such fits are a poor representation of the larger voids that dominate the volume of the void excursion set.
△ Less
Submitted 28 September, 2005;
originally announced September 2005.
-
Gravitational cooling and density profile near caustics in collisionless dark matter haloes
Authors:
Roya Mohayaee,
Sergei Shandarin
Abstract:
Cold dark matter haloes are populated by high-density structures with sharply-peaked profiles known as caustics which have not yet been resolved by 3-dimensional numerical simulations. Here, we derive semi-analytic expressions for the density profiles near caustics in haloes which form by self-similar accretions of dark matter with infinitesimal velocity dispersion. A simple rescaling shows that…
▽ More
Cold dark matter haloes are populated by high-density structures with sharply-peaked profiles known as caustics which have not yet been resolved by 3-dimensional numerical simulations. Here, we derive semi-analytic expressions for the density profiles near caustics in haloes which form by self-similar accretions of dark matter with infinitesimal velocity dispersion. A simple rescaling shows that similarly to the case of absolutely cold medium these profiles are universal: they are valid for all caustics and irrespective of physical parameters of the halo. We derive the maximum density of the caustics and show that it depends on the velocity dispersion and the caustic location. We show that both the absolute and relative thickness of the caustic monotonically decrease toward the center of the halo while the maximum density grows. This indicates that the radial component of the thermal velocities decreases in the inner streams, i.e. the collisionless medium cools down in the radial direction descending to the center of the halo. Finally, we demonstrate that there can be a significant contribution to the emission measure from dark matter particle annihilation in the caustics.
△ Less
Submitted 18 October, 2005; v1 submitted 8 March, 2005;
originally announced March 2005.
-
Morphological Statistics of the Cosmic Web
Authors:
Sergei F. Shandarin
Abstract:
We report the {\em first} systematic study of the supercluster-void network in the $Λ$CDM concordance cosmology treating voids and superclusters on an equal footing. We study the dark matter density field in real space smoothed with the $L_s = 5$ \hm1 Mpc Gaussian window. Superclusters and voids are defined as individual members of over-dense and under-dense excursion sets respectively. We deter…
▽ More
We report the {\em first} systematic study of the supercluster-void network in the $Λ$CDM concordance cosmology treating voids and superclusters on an equal footing. We study the dark matter density field in real space smoothed with the $L_s = 5$ \hm1 Mpc Gaussian window. Superclusters and voids are defined as individual members of over-dense and under-dense excursion sets respectively. We determine the morphological properties of the cosmic web at a large number of dark matter density levels by computing Minkowski functionals for every supercluster and void. At the adopted smoothing scale individual superclusters totally occupy no more than about 5% of the total volume and contain no more than 20% of mass if the largest supercluster is excluded. Likewise, individual voids totally occupy no more than 14% of volume and contain no more than 4% of mass if the largest void is excluded.
The genus of individual superclusters can be $\sim 5$ while the genus of individual voids reaches $\sim 55$, implying significant amount of substructure in superclusters and especially in voids. Large voids are typically distinctly non-spherical.
△ Less
Submitted 28 May, 2004; v1 submitted 16 May, 2004;
originally announced May 2004.
-
Morphology and Evolution in Galaxy Clusters I: Simulated Clusters in the Adiabatic limit and with Radiative Cooling
Authors:
Nurur Rahman,
Sergei F. Shandarin,
Patrick M. Motl,
Adrian L. Melott,
;
Abstract:
We have studied morphological evolution in clusters simulated in the adiabatic limit and with radiative cooling. Cluster morphology in the redshift range, $0 < z < 0.5$, is quantified by multiplicity and ellipticity. In terms of ellipticity, our result indicates slow evolution in cluster shapes compared to those observed in the X-ray and optical wavelengths. The result is consistent with Floor,…
▽ More
We have studied morphological evolution in clusters simulated in the adiabatic limit and with radiative cooling. Cluster morphology in the redshift range, $0 < z < 0.5$, is quantified by multiplicity and ellipticity. In terms of ellipticity, our result indicates slow evolution in cluster shapes compared to those observed in the X-ray and optical wavelengths. The result is consistent with Floor, Melott & Motl (2003). In terms of multiplicity, however, the result indicate relatively stronger evolution (compared to ellipticity but still weaker than observation) in the structure of simulated clusters suggesting that for comparative studies of simulation and observation, sub-structure measures are more sensitive than the shape measures. We highlight a few possibilities responsible for the discrepancy in the shape evolution of simulated and real clusters.
△ Less
Submitted 5 May, 2004;
originally announced May 2004.
-
Morphology of the supercluster-void network in LCDM cosmology
Authors:
S. F. Shandarin,
J. V. Sheth,
V. Sahni
Abstract:
We report here the first systematic study of the supercluster-void network in the LCDM concordance cosmology in which voids and superclusters are treated on an equal footing. Superclusters are defined as individual members of an over-dense excursion set and voids are defined as individual members of a complementary under-dense excursion set at the same density threshold. We determine the geometr…
▽ More
We report here the first systematic study of the supercluster-void network in the LCDM concordance cosmology in which voids and superclusters are treated on an equal footing. Superclusters are defined as individual members of an over-dense excursion set and voids are defined as individual members of a complementary under-dense excursion set at the same density threshold. We determine the geometric, topological and morphological properties of the cosmic web at a large set of density levels by computing Minkowski functionals for every supercluster and void using SURFGEN. The properties of the largest (percolating) supercluster and the complementary void are found to be very different from properties of individual superclusters and voids. Individual superclusters totally occupy no more than about 5% of the total volume and contain no more than 20% of mass if the largest supercluster is excluded. Likewise, individual voids totally occupy no more than 14% of volume and contain no more than 4% of mass if the largest void is excluded. Although superclusters are more massive and voids are more voluminous the difference in maximum volumes is not greater than by an order of magnitude. The genus value of individual superclusters can be 5 while the genus of individual voids can reach 40, implying significant amount of substructure in superclusters and especially in voids. One of our main results is that large voids, as defined through the density field (read dark matter distribution) can be distinctly non-spherical.
△ Less
Submitted 3 December, 2003;
originally announced December 2003.
-
Measuring Shapes of Galaxy Images II: Morphology of 2MASS Galaxies
Authors:
Nurur Rahman,
Sergei F. Shandarin
Abstract:
We study a sample of 112 galaxies of various Hubble types imaged in the Two Micron All Sky Survey (2MASS) in the Near-Infra Red (NIR; 1-2 $μ$m) $J$, $H$, and $K_s$ bands. The sample contains (optically classified) 32 elliptical, 16 lenticulars, and 64 spirals acquired from the 2MASS Extended Source Catalogue. We use a set of non-parametric shape measures constructed from the Minkowski Functional…
▽ More
We study a sample of 112 galaxies of various Hubble types imaged in the Two Micron All Sky Survey (2MASS) in the Near-Infra Red (NIR; 1-2 $μ$m) $J$, $H$, and $K_s$ bands. The sample contains (optically classified) 32 elliptical, 16 lenticulars, and 64 spirals acquired from the 2MASS Extended Source Catalogue. We use a set of non-parametric shape measures constructed from the Minkowski Functionals (MFs) for galaxy shape analysis.
We use ellipticity ($ε$) and orientation angle ($Φ$) as shape diagnostics. With these parameters as functions of area within the isophotal contour, we note that the NIR elliptical galaxies with $ε> 0.2$ show a trend of being centrally spherical and increasingly flattened towards the edge, a trend similar to images in optical wavelengths. The highly flattened elliptical galaxies show strong change in ellipticity between the center and the edge. The lenticular galaxies show morphological properties resembling either ellipticals or disk galaxies. Our analysis shows that almost half of the spiral galaxies appear to have bar like features while the rest are likely to be non-barred. Our results also indicate that almost one-third of spiral galaxies have optically hidden bars.
The isophotal twist noted in the orientations of elliptical galaxies decreases with the flattening of these galaxies indicating that twist and flattening are also anti-correlated in the NIR, as found in optical wavelengths. The orientations of NIR lenticular and spiral galaxies show a wide range of twists.
△ Less
Submitted 6 July, 2004; v1 submitted 9 October, 2003;
originally announced October 2003.
-
Quantifying the Bull's Eye Effect
Authors:
Brian C. Thomas,
Adrian L. Melott,
Hume A. Feldman,
Sergei F. Shandarin
Abstract:
We have used N-body simulations to develop two independent methods to quantify redshift distortions known as the Bull's Eye effect (large scale infall plus small scale virial motion). This effect depends upon the mass density, $Ω_0$, so measuring it can in principle give an estimate of this important cosmological parameter. We are able to measure the effect and distinguish between its strength f…
▽ More
We have used N-body simulations to develop two independent methods to quantify redshift distortions known as the Bull's Eye effect (large scale infall plus small scale virial motion). This effect depends upon the mass density, $Ω_0$, so measuring it can in principle give an estimate of this important cosmological parameter. We are able to measure the effect and distinguish between its strength for high and low values of $Ω_0$. Unlike other techniques which utilize redshift distortions, one of our methods is relatively insensitive to bias. In one approach, we use path lengths between contour crossings of the density field. The other is based upon percolation. We have found both methods to be successful in quantifying the effect and distinguishing between values of $Ω_0$. However, only the path lengths method exhibits low sensitivity to bias.
△ Less
Submitted 9 October, 2003; v1 submitted 23 May, 2003;
originally announced May 2003.
-
Measuring Shapes of Galaxy Images I: Ellipticity and Orientation
Authors:
Nurur Rahman,
Sergei F. Shandarin
Abstract:
We suggest a set of morphological measures that we believe can help in quantifying the shapes of two-dimensional cosmological images such as galaxies, clusters, and superclusters of galaxies. The method employs non-parametric morphological descriptors known as the Minkowski functionals in combination with geometric moments widely used in the image analysis. For the purpose of visualization of th…
▽ More
We suggest a set of morphological measures that we believe can help in quantifying the shapes of two-dimensional cosmological images such as galaxies, clusters, and superclusters of galaxies. The method employs non-parametric morphological descriptors known as the Minkowski functionals in combination with geometric moments widely used in the image analysis. For the purpose of visualization of the morphological properties of image contour lines we introduce three auxiliary ellipses representing the vector and tensor Minkowski functionals. We study the discreteness, seeing, and noise effects on elliptic contours as well as their morphological characteristics such as the ellipticity and orientation. In order to reduce the effect of noise we employ a technique of contour smoothing. We test the method by studying simulated elliptic profiles of toy spheroidal galaxies ranging in ellipticity from E0 to E7. We then apply the method to real galaxies, including eight spheroidals, three disk spirals and one peculiar galaxy, as imaged in the near-infrared $K_s$-band (2.2 microns) with the Two Micron All Sky Survey (2MASS). The method is numerically very efficient and can be used in the study of hundreds of thousands images obtained in modern surveys.
△ Less
Submitted 9 May, 2003; v1 submitted 17 March, 2003;
originally announced March 2003.
-
Measuring the Geometry and Topology of Large Scale Structure using SURFGEN: Methodology and Preliminary Results
Authors:
Jatush V. Sheth,
Varun Sahni,
Sergei F. Shandarin,
B. S. Sathyaprakash
Abstract:
We present a new ansatz which can successfully be used to determine the morphological properties of the supercluster-void network. The ansatz is based on a surface modelling scheme SURFGEN, which generates a triangulated surface from a discrete data set representing (say) the distribution of galaxies in real (or redshift) space. Four Minkowski functionals -- surface area, volume, extrinsic curva…
▽ More
We present a new ansatz which can successfully be used to determine the morphological properties of the supercluster-void network. The ansatz is based on a surface modelling scheme SURFGEN, which generates a triangulated surface from a discrete data set representing (say) the distribution of galaxies in real (or redshift) space. Four Minkowski functionals -- surface area, volume, extrinsic curvature and genus -- describe the geometry and topology of the supercluster-void network. Ratio's of Minkowski functionals -- Shapefinders -- provide us with an excellent diagnostic of three dimensional shapes of clusters, superclusters and voids. Minkowski functionals and Shapefinders are determined for a triangulated iso-density surface using SURFGEN. SURFGEN is tested against both simply and multiply connected eikonal surfaces such as triaxial ellipsoids and tori. Remarkably, the first three Minkowski functionals are computed to better than 1% accuracy while the fourth (genus) is known exactly. SURFGEN also gives excellent results when applied to Gaussian random fields. Our results indicate that the surface modelling scheme SURFGEN is accurate and robust and can successfully be used to quantify the topology and morphology of the supercluster-void network in the universe. We apply SURFGEN to three cosmological models, $Ł$CDM, $\T$CDM and SCDM and obtain interesting new results pertaining to the geometry, morphology and topology of large scale structure.
△ Less
Submitted 22 March, 2003; v1 submitted 6 October, 2002;
originally announced October 2002.
-
Effects of Noise on Galaxy Isophotes
Authors:
Nurur Rahman,
Sergei F. Shandarin
Abstract:
The study of shapes of the images of objects is an important issue not only because it reveals its dynamical state but also it helps to understand the object's evolutionary history. We discuss a new technique in cosmological image analysis which is based on a set of non-parametric shape descriptors known as the Minkowski Functionals (MFs). These functionals are extremely versatile and under some…
▽ More
The study of shapes of the images of objects is an important issue not only because it reveals its dynamical state but also it helps to understand the object's evolutionary history. We discuss a new technique in cosmological image analysis which is based on a set of non-parametric shape descriptors known as the Minkowski Functionals (MFs). These functionals are extremely versatile and under some conditions give a complete description of the geometrical properties of objects. We believe that MFs could be a useful tool to extract information about the shapes of galaxies, clusters of galaxies and superclusters. The information revealed by MFs can be utilized along with the knowledge obtained from currently popular methods and thus could improve our understanding of the true shapes of cosmological objects.
△ Less
Submitted 25 September, 2002;
originally announced September 2002.
-
Morphology of the Secondary CMB Anisotropies: the Key to "Smoldering" Reionization
Authors:
Nickolay Y. Gnedin,
Sergei F. Shandarin
Abstract:
We show how the morphological analysis of the maps of the secondary CMB anisotropies can detect an extended period of ``smoldering'' reionization, during which the universe remains partially ionized. Neither radio observations of the redshifted 21cm line nor IR observations of the redshifted Lyman-alpha forest will be able to detect such a period. The most sensitive to this kind of non-gaussiani…
▽ More
We show how the morphological analysis of the maps of the secondary CMB anisotropies can detect an extended period of ``smoldering'' reionization, during which the universe remains partially ionized. Neither radio observations of the redshifted 21cm line nor IR observations of the redshifted Lyman-alpha forest will be able to detect such a period. The most sensitive to this kind of non-gaussianity parameters are the number of regions in the excursion set, the perimeter of the excursion set, and the genus of the largest (by area) region. For example, if the universe reionized fully at z=6, but maintained about 1/3 mean ionized fraction since z=20, then a 2 arcmin map with 500x500 pixel resolution and a signal-to-noise ratio S/N=1/2 allows to detect the non-gaussianity due to reionization with better than 99% confidence level.
△ Less
Submitted 13 June, 2002;
originally announced June 2002.
-
Testing non-Gaussianity in CMB Maps by Morphological Statistic
Authors:
Sergei F. Shandarin
Abstract:
The assumption of Gaussianity of the primordial perturbations plays an important role in modern cosmology. The most direct test of this hypothesis consists in testing Gaussianity of the CMB maps. Counting the pixels with the temperatures in given ranges and thus estimating the one point probability function of the field is the simplest of all the tests. Other usually more complex tests of Gaussi…
▽ More
The assumption of Gaussianity of the primordial perturbations plays an important role in modern cosmology. The most direct test of this hypothesis consists in testing Gaussianity of the CMB maps. Counting the pixels with the temperatures in given ranges and thus estimating the one point probability function of the field is the simplest of all the tests. Other usually more complex tests of Gaussianity generally use a great deal of the information already contained in the probability function. However, the most interesting outcome of such a test would be the signal of non-Gaussianity independent of the probability function. It is shown that the independent information has purely morphological character i.e. it depends on the geometry and topology of the level contours only. As an example we discuss in detail the quadratic model $v=u+α(u^2-1)$ ($u$ is a Gaussian field with $\bar{u}=0$ and $<u^2>=1$, $α$ is a parameter) which may arise in slow-roll or two-field inflation models. We show that in the limit of small amplitude $α$ the full information about the non-Gaussianity is contained in the probability function. If other tests are performed on this model they simply recycle the same information. A simple procedure allowing to assess the sensitivity of any statistics to the morphological information is suggested. We provide an analytic estimate of the statistical limit for detecting the quadratic non-Gaussianity $\a_c$ as a function of the map size in the ideal situation when the scale of the field is resolved. This estimate is in a good agreement with the results of the Monte Carlo simulations of $256^2$ and $1024^2$ maps. The effect of resolution on the detection quadratic non-Gaussianity is also briefly discussed.
△ Less
Submitted 31 May, 2002; v1 submitted 17 July, 2001;
originally announced July 2001.
-
Morphological Measures of nonGaussianity in CMB Maps
Authors:
Sergei F. Shandarin,
Hume A. Feldman,
Yongzhong Xu,
Max Tegmark
Abstract:
We discuss the tests of nonGaussianity in CMB maps using morphological statistics known as Minkowski functionals. As an example we test degree-scale cosmic microwave background (CMB) anisotropy for Gaussianity by studying the \qmask map that was obtained from combining the QMAP and Saskatoon data. We compute seven morphological functions $M_i(\dt)$, $i=1,...,7$: six \mf and the number of regions…
▽ More
We discuss the tests of nonGaussianity in CMB maps using morphological statistics known as Minkowski functionals. As an example we test degree-scale cosmic microwave background (CMB) anisotropy for Gaussianity by studying the \qmask map that was obtained from combining the QMAP and Saskatoon data. We compute seven morphological functions $M_i(\dt)$, $i=1,...,7$: six \mf and the number of regions $N_c$ at a hundred $\dt$ levels. We also introduce a new parameterization of the morphological functions $M_i(A)$ in terms of the total area $A$ of the excursion set. We show that the latter considerably decorrelates the morphological statistics and makes them more robust because they are less sensitive to the measurements at extreme levels. We compare these results with those from 1000 Gaussian Monte Carlo maps with the same sky coverage, noise properties and power spectrum, and conclude that the \qmask map is neither a very typical nor a very exceptional realization of a Gaussian field. At least about 20% of the 1000 Gaussian Monte Carlo maps differ more than the \qmask map from the mean morphological parameters of the Gaussian fields.
△ Less
Submitted 1 June, 2002; v1 submitted 6 July, 2001;
originally announced July 2001.
-
Constraining Cosmological Models by the Cluster Mass Function
Authors:
Nurur Rahman,
Sergei F. Shandarin
Abstract:
We present a comparison between two observational and three theoretical mass functions for eight cosmological models suggested by the data from the recently completed BOOMERANG-98 and MAXIMA-1 cosmic microwave background (CMB) anisotropy experiments as well as peculiar velocities (PVs) and type Ia supernovae (SN) observations. The cosmological models have been proposed as the best fit models by…
▽ More
We present a comparison between two observational and three theoretical mass functions for eight cosmological models suggested by the data from the recently completed BOOMERANG-98 and MAXIMA-1 cosmic microwave background (CMB) anisotropy experiments as well as peculiar velocities (PVs) and type Ia supernovae (SN) observations. The cosmological models have been proposed as the best fit models by several groups. We show that no model is in agreement with the abundances of X-ray clusters at $\sim 10^{14.7} h^{-1}M_{\odot}$.On the other hand, we find that the BOOM+MAX+{\sl COBE}:I, Refined Concordance and $Λ$MDM are in a good agreement with the abundances of optical clusters. The P11 and especially Concordance models predict a slightly lower abundances than observed at $\sim 10^{14.6} h^{-1}M_{\odot}$. The BOOM+MAX+{\sl COBE}:II and PV+CMB+SN models predict a slightly higher abundances than observed at $\sim 10^{14.9} h^{-1}M_{\odot}$. The nonflat MAXIMA-1 is in a fatal conflict with the observational cluster abundances and can be safely ruled out.
△ Less
Submitted 2 April, 2001; v1 submitted 11 October, 2000;
originally announced October 2000.
-
The Cosmological Mass Function in the Zel'dovich Approximation
Authors:
Sergei F. Shandarin
Abstract:
The Press-Schechter theory of the cosmological mass function and its modifications allow to constraint cosmological scenarios of the structure formation. Recently a few new models have been suggested that explored the influence of anisotropic collapse on the shape of the mass function. I discuss in more detail a particular model that assumes that a fluid particle becomes a part of a gravitationa…
▽ More
The Press-Schechter theory of the cosmological mass function and its modifications allow to constraint cosmological scenarios of the structure formation. Recently a few new models have been suggested that explored the influence of anisotropic collapse on the shape of the mass function. I discuss in more detail a particular model that assumes that a fluid particle becomes a part of a gravitationally bound halo when the smallest eigenvalue of the deformation tensor of the filtered initial density field reaches a certain threshold (like the filtered density contrast reaches the threshold in the Press-Schechter formalism). Choosing the smallest eigenvalue guarantees that the fluid particle in question experiences collapse along all three axes. The model shows a better agreement with the N-body simulations than the standard Press-Schechter model.
△ Less
Submitted 31 July, 2000;
originally announced July 2000.
-
Evidence for Filamentarity in the Las Campanas Redshift Survey
Authors:
Somnath Bharadwaj,
Varun Sahni,
B. S. Sathyaprakash,
Sergei F. Shandarin,
Capp Yess
Abstract:
We apply Shapefinders, statistical measures of `shape' constructed from two dimensional partial Minkowski functionals, to study the degree of filamentarity in the Las Campanas Redshift Survey (LCRS). In two dimensions, three Minkowski functionals characterise the morphology of an object, they are: its perimeter (L), area (S), and genus. Out of L and S a single dimensionless Shapefinder Statistic…
▽ More
We apply Shapefinders, statistical measures of `shape' constructed from two dimensional partial Minkowski functionals, to study the degree of filamentarity in the Las Campanas Redshift Survey (LCRS). In two dimensions, three Minkowski functionals characterise the morphology of an object, they are: its perimeter (L), area (S), and genus. Out of L and S a single dimensionless Shapefinder Statistic, F can be constructed (0 <=F <=1). F acquires extreme values on a circle (F = 0) and a filament (F = 1). Using F, we quantify the extent of filamentarity in the LCRS by comparing our results with a Poisson distribution with similar geometrical properties and having the same selection function as the survey. Our results unambiguously demonstrate that the LCRS displays a high degree of filamentarity both in the Northern and Southern galactic sections a result that is in general agreement with the visual appearance of the catalogue. It is well known that gravitational clustering from Gaussian initial conditions gives rise to the development of non-Gaussianity reflected in the formation of a network-like filamentary structure on supercluster scales. Consequently the fact that the smoothed LCRS catalogue shows properties consistent with those of a Gaussian random field (Colley 1997) whereas the unsmoothed catalogue demonstrates the presence of filamentarity lends strong support to the conjecture that the large scale clustering of galaxies is driven by gravitational instability.
△ Less
Submitted 23 November, 1999; v1 submitted 29 April, 1999;
originally announced April 1999.
-
Disentangling the Cosmic Web I: Morphology of Isodensity Contours
Authors:
Jens Schmalzing,
Thomas Buchert,
Adrian L. Melott,
Varun Sahni,
B. S. Sathyaprakash,
Sergei F. Shandarin
Abstract:
We apply Minkowski functionals and various derived measures to decipher the morphological properties of large-scale structure seen in simulations of gravitational evolution. Minkowski functionals of isodensity contours serve as tools to test global properties of the density field. Furthermore, we identify coherent objects at various threshold levels and calculate their partial Minkowski function…
▽ More
We apply Minkowski functionals and various derived measures to decipher the morphological properties of large-scale structure seen in simulations of gravitational evolution. Minkowski functionals of isodensity contours serve as tools to test global properties of the density field. Furthermore, we identify coherent objects at various threshold levels and calculate their partial Minkowski functionals. We propose a set of two derived dimensionless quantities, planarity and filamentarity, which reduce the morphological information in a simple and intuitive way. Several simulations of the gravitational evolution of initial power-law spectra provide a framework for systematic tests of our method.
△ Less
Submitted 14 July, 1999; v1 submitted 27 April, 1999;
originally announced April 1999.
-
Multiple-streaming and the Probability Distribution of Density in Redshift Space
Authors:
Lam Hui,
Lev Kofman,
Sergei F. Shandarin
Abstract:
We examine several aspects of redshift distortions by expressing the redshift-space density in terms of the eigenvalues and orientation of the local Lagrangian deformation tensor. We explore the importance of multiple-streaming using the Zel'dovich approximation (ZA), and compute the average number of streams in real and redshift-space. It is found that multiple-streaming can be significant in r…
▽ More
We examine several aspects of redshift distortions by expressing the redshift-space density in terms of the eigenvalues and orientation of the local Lagrangian deformation tensor. We explore the importance of multiple-streaming using the Zel'dovich approximation (ZA), and compute the average number of streams in real and redshift-space. It is found that multiple-streaming can be significant in redshift-space but negligible in real-space, even at moderate values of the linear fluctuation amplitude ($σ< 1$). Moreover, unlike their real-space counter-parts, redshift-space multiple-streams can flow past each other with minimal interactions. Such nonlinear redshift-space effects, which operate even when the real-space density field is quite linear, could suppress the classic compression of redshift-structures predicted by linear theory (Kaiser 1987). We also compute using the ZA the probability distribution function (PDF) of density, as well as $S_3$, in real and redshift-space, and compare it with the PDF measured from N-body simulations. The role of caustics in defining the character of the high density tail is examined. It is found that (non-Lagrangian) smoothing, due to both finite resolution or discreteness and small-scale velocity dispersions, is very effective in erasing caustic structures, unless the initial power spectrum is sufficiently truncated.
△ Less
Submitted 9 January, 1999;
originally announced January 1999.
-
Comparison of Analytical Mass Functions with Numerical Simulations
Authors:
Jounghun Lee,
Sergei F. Shandarin
Abstract:
We present numerical testing results of our mass function derived in our previous paper, and compare the testing results with those of the popular Press-Schechter (PS) mass function. Two fiducial models are considered for the test: the scale-free power-law spectra P(k) \propto k^{n} with spectral indices n=-1, 0 and the standard cold dark matter (SCDM) model with Omega =1 and h=0.5. For the powe…
▽ More
We present numerical testing results of our mass function derived in our previous paper, and compare the testing results with those of the popular Press-Schechter (PS) mass function. Two fiducial models are considered for the test: the scale-free power-law spectra P(k) \propto k^{n} with spectral indices n=-1, 0 and the standard cold dark matter (SCDM) model with Omega =1 and h=0.5. For the power-law models, we use numerical data averaged over several different output times: ten output times from two N-body realizations for the n=-1 power-law model; four outputs from one realization for the n=0 model. While for the SCDM model, we consider four outputs separately at redshifts z = 0, 0.43, 1.14 and 1.86 from one large N-body simulations. The comparison results show that our mass function fits the numerical data in a much improved way over the PS one. Thus, we expect that our mass function can be a viable alternative of the PS mass function in applications to various areas.
△ Less
Submitted 23 March, 1999; v1 submitted 31 October, 1998;
originally announced November 1998.
-
Minkowski Functionals and Cluster Analysis for CMB Maps
Authors:
D. Novikov,
Hume A. Feldman,
Sergei F. Shandarin
Abstract:
We suggest novel statistics for the CMB maps that are sensitive to non-Gaussian features. These statistics are natural generalizations of the geometrical and topological methods that have been already used in cosmology such as the cumulative distribution function and genus. We compute the distribution functions of the Partial Minkowski Functionals for the excursion set above or bellow a constant…
▽ More
We suggest novel statistics for the CMB maps that are sensitive to non-Gaussian features. These statistics are natural generalizations of the geometrical and topological methods that have been already used in cosmology such as the cumulative distribution function and genus. We compute the distribution functions of the Partial Minkowski Functionals for the excursion set above or bellow a constant temperature threshold. Minkowski Functionals are additive and are translationally and rotationally invariant. Thus, they can be used for patchy and/or incomplete coverage. The technique is highly efficient computationally (it requires only O(N) operations, where N is the number of pixels per one threshold level). Further, it allows to split large data sets into smaller subsets. The full advantage of these statistics can be obtained only on very large data sets. We apply it to the 4-year DMR COBE data corrected for the Galaxy contamination as an illustration of the technique.
△ Less
Submitted 11 June, 1999; v1 submitted 18 September, 1998;
originally announced September 1998.
-
Morphology of Clusters and Superclusters in N-body simulations of Cosmological Gravitational Clustering
Authors:
B. S. Sathyaprakash,
Varun Sahni,
Sergei F. Shandarin
Abstract:
We analyse shapes of overdense regions (clusters and superclusters) in controlled N-body simulations of gravitational clustering with power law initial spectra P(k) \propto k^n, n = -3, -2, -1, 0. At values of the density just above the percolation transition the number of distinct (isolated) clusters peaks and we use this `natural threshold' to study the shapes and multiplicity function of clus…
▽ More
We analyse shapes of overdense regions (clusters and superclusters) in controlled N-body simulations of gravitational clustering with power law initial spectra P(k) \propto k^n, n = -3, -2, -1, 0. At values of the density just above the percolation transition the number of distinct (isolated) clusters peaks and we use this `natural threshold' to study the shapes and multiplicity function of clusters and superclusters. We find that the extent of both filamentarity and pancakeness increases as the simulation evolves, the former being appreciably larger than the latter at virtually all epochs and for all spectra considered by us. Our results also show that high density regions within very massive clusters/superclusters are likely to be noticeably filamentary or pancake/ribbon-like when compared to the less dense regions within these objects. We make a detailed study of two, moment-based `shape statistics' proposed, respectively, by Babul & Starkman (BS) and Luo & Vishniac (LV) and find that both LV and BS have certain built-in limitations: LV is biased towards oblate structures and tends to overemphasise this property; neither BS nor LV correctly describe the shape of strongly curved or topologically nontrivial objects. For instance, a thin filamentary torus and a ribbon are both described by BS and LV as being pancakes ! By contrast Shapefinders, a new shape diagnostic not constructed from density moments but from Minskowski functionals, does not suffer from these limitations and appears to faithfully reproduce the shapes of both simple and topologically complex objects.
△ Less
Submitted 21 May, 1998;
originally announced May 1998.
-
Filaments and Pancakes in the IRAS 1.2 Jy Redshift Catalogue
Authors:
B. S. Sathyaprakash,
Varun Sahni,
Sergei F. Shandarin,
Karl B. Fisher
Abstract:
We explore shapes of clusters and superclusters in the IRAS 1.2 Jy redshift survey with three reconstructions spanning the range β= 0.1, 0.5, 1.0, where β= Ω^{0.6}/b, b is the bias factor and Ωthe present value of the dimensionless matter density. Comparing our results to Gaussian randomized reconstructions of the IRAS catalogue, we find structures having both planar and filamentary properties.…
▽ More
We explore shapes of clusters and superclusters in the IRAS 1.2 Jy redshift survey with three reconstructions spanning the range β= 0.1, 0.5, 1.0, where β= Ω^{0.6}/b, b is the bias factor and Ωthe present value of the dimensionless matter density. Comparing our results to Gaussian randomized reconstructions of the IRAS catalogue, we find structures having both planar and filamentary properties. For β= 0.5, 1.0 the largest structures in the survey have a distinct tendency to be filament-like} in general agreement with the results of N-body simulations.
△ Less
Submitted 20 May, 1998;
originally announced May 1998.
-
The Large Scale Biasing and the Primordial Gravitational Potential
Authors:
Jounghun Lee,
Sergei F. Shandarin
Abstract:
We modified the Press-Schechter (PS) formalism and then analytically derived a constrained mass distribution function n(M|varphi) for the regions having some specified value of the primordial gravitational potential, varphi. The resulting modified PS theory predicts that gravitationally bound clumps with masses corresponding to rich clusters are significantly biased toward the regions of negativ…
▽ More
We modified the Press-Schechter (PS) formalism and then analytically derived a constrained mass distribution function n(M|varphi) for the regions having some specified value of the primordial gravitational potential, varphi. The resulting modified PS theory predicts that gravitationally bound clumps with masses corresponding to rich clusters are significantly biased toward the regions of negative primordial potential - the troughs of the potential. The prediction is quantitative, depending on the mass and the depth of the troughs, which can be tested in large N-body simulations. As an illustration of the magnitude of the effect we calculate the constrained mass function for the CDM model with Gamma = Omega h = 0.25 normalized to sigma_{8} = 1. In particular, we show that the probability of finding a clump of mass 10^{14} - 10^{15}h^{-1}M_{odot} in the region of negative initial potential is 1.3 - 3 times greater (depending on the mass) than that in the region of positiveinitial potential. The scale of the potential fluctuations R_{varphi}=sqrt{3} sigma_{varphi}/sigma_{varphi'} is shown to be approx 120 h^{-1} Mpc for the spectrum in question. The rms mass density contrast on this scale is only about sigma_{delta}(R_{varphi}) approx 0.03. Assuming that the modified PS theory is statistically correct, we conclude that clusters are significantly biased (b ge 10, b is a bias factor defined by Delta n_{cl}/ n_{cl} = b Delta rho_m/rho_m) toward the regions having negative initial potential.
△ Less
Submitted 25 July, 1998; v1 submitted 18 March, 1998;
originally announced March 1998.
-
Shapefinders: A new shape diagnostic for large scale structure
Authors:
Varun Sahni,
B. S. Sathyaprakash,
Sergei F. Shandarin
Abstract:
We construct a set of shape-finders which determine shapes of compact surfaces (iso-density surfaces in galaxy surveys or N-body simulations) without fitting them to ellipsoidal configurations as done earlier. The new indicators arise from simple, geometrical considerations and are derived from fundamental properties of a surface such as its volume, surface area, integrated mean curvature and co…
▽ More
We construct a set of shape-finders which determine shapes of compact surfaces (iso-density surfaces in galaxy surveys or N-body simulations) without fitting them to ellipsoidal configurations as done earlier. The new indicators arise from simple, geometrical considerations and are derived from fundamental properties of a surface such as its volume, surface area, integrated mean curvature and connectivity characterized by the Genus. These `Shapefinders' could be used to diagnose the presence of filaments, pancakes and ribbons in large scale structure. Their lower-dimensional generalization may be useful for the study of two-dimensional distributions such as temperature maps of the Cosmic Microwave Background.
△ Less
Submitted 7 January, 1998;
originally announced January 1998.
-
Dark Matter, Discreteness and Collision Error in Cosmological N-Body Simulations
Authors:
Randall J. Splinter,
Adrian L. Melott,
Sergei F. Shandarin
Abstract:
We report on a series of tests of agreement between three types of N-body simulations: PM, P$^3$M, and Tree codes. We find good agreement in both the individual and the statistical properties only on scales larger than the mean interparticle separation. As a result, we question most numerical results at and below below galaxy scales, either concerning primordial dark matter or baryonic matter co…
▽ More
We report on a series of tests of agreement between three types of N-body simulations: PM, P$^3$M, and Tree codes. We find good agreement in both the individual and the statistical properties only on scales larger than the mean interparticle separation. As a result, we question most numerical results at and below below galaxy scales, either concerning primordial dark matter or baryonic matter coupled to it by gravitation. Mass resolution (specifically, sufficient to have the mean interparticle separation comparable to the smoothing scales) appears be necessary for results to converge.
△ Less
Submitted 10 November, 1997;
originally announced November 1997.
-
The Cosmological Mass Distribution Function in the Zel'dovich Approximation
Authors:
Jounghun Lee,
Sergei F. Shandarin
Abstract:
An analytic approximation to the mass function for gravitationally bound objects is presented. We base on the Zel'dovich approximation to extend the Press-Schechter formalism to a nonspherical dynamical model. A simple extrapolation of that approximation suggests that the gravitational collapse along all three directions which eventually leads to the formation of real virialized objects - clumps…
▽ More
An analytic approximation to the mass function for gravitationally bound objects is presented. We base on the Zel'dovich approximation to extend the Press-Schechter formalism to a nonspherical dynamical model. A simple extrapolation of that approximation suggests that the gravitational collapse along all three directions which eventually leads to the formation of real virialized objects - clumps occur in the regions where the lowest eigenvalue of the deformation tensor,lambda_{3}, is positive. We derive the conditional probability of lambda_{3}>0 as a function of the linearly extrapolated density contrast, delta, and the conditional probability distribution of delta provided that lambda_{3}>0. These two conditional probability distributions show that the most probable density of the bound regions (lambda_{3}>0) is roughly 1.5 at the characteristic mass scale, and that the probability of lambda_{3}>0 is almost unity in the highly overdense regions (delta>3*sigma). Finally an analytic mass function of clumps is derived with a help of one simple ansatz which is employed to treat the multistream regions beyond the validity of the Zel'dovich approximation. The resulting mass function is renormalized by a factor of 12.5, which we justify with a sharp k-space filter by means of the modified Jedamzik analysis. Our mass function is shown to be different from the Press-Schechter one, having a lower peak and predicting more small-mass objects.
△ Less
Submitted 8 December, 1997; v1 submitted 19 September, 1997;
originally announced September 1997.
-
Fundamental Discreteness Limitations of Cosmological N-Body Clustering Simulations
Authors:
Randall J. Splinter,
Adrian L. Melott,
Sergei F. Shandarin,
Yasushi Suto
Abstract:
We explore some of the effects that discreteness and two-body scattering may have on N-body simulations with ``realistic'' cosmological initial conditions. We use an identical subset of particles from the initial conditions for a $128^3$ Particle-Mesh (PM) calculation as the initial conditions for a variety P$^3$M and Tree code runs. We investigate the effect of mass resolution (the mean interpa…
▽ More
We explore some of the effects that discreteness and two-body scattering may have on N-body simulations with ``realistic'' cosmological initial conditions. We use an identical subset of particles from the initial conditions for a $128^3$ Particle-Mesh (PM) calculation as the initial conditions for a variety P$^3$M and Tree code runs. We investigate the effect of mass resolution (the mean interparticle separation) since most ``high resolution'' codes only have high resolution in gravitational force. The phase-insensitive two--point statistics, such as the power spectrum (autocorrelation) are somewhat affected by these variations, but phase-sensitive statistics show greater differences. Results converge at the mean interparticle separation scale of the lowest mass-resolution code. As more particles are added, but the force resolution is held constant, the P$^3$M and the Tree runs agree more and more strongly with each other and with the PM run which had the same initial conditions. This shows high particle density is necessary for correct time evolution, since many different results cannot all be correct. However, they do not so converge to a PM run which continued the fluctuations to small scales. Our results show that ignoring them is a major source of error on comoving scales of the missing wavelengths. This can be resolved by putting in a high particle density. Since the codes never agree well on scales below the mean comoving interparticle separation, we find little justification for quantitative predictions on this scale. Some measures vary by 50%, but others can be off by a factor of three or more. Our results suggest possible problems with the density of galaxy halos, formation of early generation objects such as QSO absorber clouds, etc.
△ Less
Submitted 5 November, 1997; v1 submitted 10 June, 1997;
originally announced June 1997.
-
Detection of Network Structure in the Las Campanas Redshift Survey
Authors:
Sergei F. Shandarin,
Capp Yess
Abstract:
We employ a percolation technique developed for pointwise distributions to analyze two-dimensional projections of the three northern and three southern slices in the Las Campanas Redshift Survey. One of the goals of this paper is to compare the visual impressions of the structure within distributions with objective statistical analysis. We track the growth of the largest cluster as an indicator…
▽ More
We employ a percolation technique developed for pointwise distributions to analyze two-dimensional projections of the three northern and three southern slices in the Las Campanas Redshift Survey. One of the goals of this paper is to compare the visual impressions of the structure within distributions with objective statistical analysis. We track the growth of the largest cluster as an indicator of the network structure. We restrict our analysis to volume limited subsamples in the regions from 200 to 400 Mpc/h where the number density of galaxies is the highest. As a major result, we report a measurement of an unambiguous signal, with high signal-to-noise ratio (at least at the level of a few $σ$), indicating significant connectivity of the galaxy distribution which in two dimensions is indicative of a filamentary distribution. This is in general agreement with the visual impression and typical for the standard theory of the large-scale structure formation based on gravitational instability of initially Gaussian density fluctuations.
△ Less
Submitted 21 May, 1997; v1 submitted 20 May, 1997;
originally announced May 1997.