No abstract

We present a method called Significant Cosmic Holes in Universe (SCHU) for identifying cosmic voids and loops of filaments in cosmological datasets and assigning their statistical significance using techniques from topological data analysis. In particular, persistent homology is used to find different dimensional holes. For dark matter halo catalogs and galaxy surveys, the 0-, 1-, and 2dimensional holes can be identified with connected components (i.e. clusters), loops of filaments, and voids. The procedure overlays dark matter halos/galaxies on a three-dimensional grid, and a distance-to-measure (DTM) function is calculated at each point of the grid. A nested set of simplicial complexes (a filtration) is generated over the lower-level sets of the DTM across increasing threshold values. The filtered simplicial complex can then be used to summarize the birth and death times of the different dimension homology group generators (i.e., the holes). Persistent homology summary diagrams, called persistence diagrams, are produced from the dimension, birth times, and death times of each homology group generator. Using the persistence diagrams and bootstrap sampling, we explain how p-values can be assigned to each homology group generator. The homology group generators on a persistence diagram are not, in general, uniquely located back in the original dataset volume so we propose a method for finding a representation of the homology group generators. This method provides a novel, statistically rigorous approach for locating informative generators in cosmological datasets, which may be useful for providing complementary cosmological constraints on the effects of, for example, the sum of the neutrino masses. The method is tested on a Voronoi foam simulation, and then subsequently applied to a subset of the SDSS galaxy survey and a cosmological simulation. Lastly, we calculate Betti functions for two of the MassiveNuS simulations and discuss implications for using the persistent homology of the density field to help break degeneracy in the cosmological parameters.

Based on the photometric redshift catalog of Zou et al., we apply a fast clustering algorithm to identify 540,432 galaxy clusters at z ≲ 1 in the DESI legacy imaging surveys, which cover a sky area of about 20,000 deg2. Monte Carlo simulations indicate that the false-detection rate of our detecting method is about 3.1%. The total masses of galaxy clusters are derived using a calibrated richness–mass relation that is based on the observations of X-ray emission and the Sunyaev and Zel’dovich effect. The median redshift and mass of our detected clusters are about 0.53 and 1.23 × 1014 M ⊙, respectively. Comparing with previous clusters identified using the data of the Sloan Digital Sky Survey; we can recognize most of them, especially those with high richness. Our catalog will be used for further statistical studies on galaxy clusters and environmental effects on galaxy evolution, etc.

Continuum normalization of echelle spectra is an important data analysis step that is difficult to automate. Polynomial fitting requires a reasonably high order model to follow the steep slope of the blaze function. However, in the presence of deep spectral lines, a high order polynomial fit can result in ripples in the normalized continuum that increase errors in spectral analysis. Here, we present two algorithms for flattening the spectrum continuum. The Alpha-shape Fitting to Spectrum algorithm (AFS) is completely data-driven, using an alpha shape to obtain an initial estimate of the blaze function. The Alpha-shape and Lab Source Fitting to Spectrum algorithm (ALSFS) incorporates a continuum constraint from a lab source reference spectrum for the blaze function estimation. These algorithms are tested on a simulated spectrum, where we demonstrate improved normalization compared to polynomial regression for continuum