Search

Foreground removal presents a significant obstacle in both current and forthcoming intensity mapping surveys. While numerous techniques have been developed that show promise in simulated datasets, their efficacy often diminishes when applied to real-world data. A primary issue is the frequency-dependent variations in the instrumental response. In this paper, we propose a novel approach utilizing the internal cross-correlation among different frequencies to calibrate the beam's frequency fluctuations. Using a simulated dataset that incorporates frequency-dependent random fluctuations into the beam model, we illustrate that our method can achieve considerable improvements over traditional techniques. Our results represent a step forward in enhancing the precision and reliability of foreground removal in intensity mapping surveys., Comment: Accepted for publication in ApJS

We propose a lightweight deep convolutional neural network (lCNN) to estimate cosmological parameters from simulated three-dimensional dark matter (DM) halo distributions and associated statistics. The training dataset comprises 2000 realizations of a cubic box with a side length of 1000 $h^{-1}{\rm Mpc}$, and interpolated over a cubic grid of $300^3$ voxels, with each simulation produced using $512^3$ DM particles and $512^3$ neutrinos. Under the flat $\Lambda$CDM model, simulations vary standard six cosmological parameters including $\Omega_m$, $\Omega_b$, $h$, $n_s$, $\sigma_8$, $w$, along with the neutrino mass sum, $M_\nu$. We find that: 1) within the framework of lCNN, extracting large-scale structure information is more efficient from the halo density field compared to relying on the statistical quantities including the power spectrum, the two-point correlation function, and the coefficients from wavelet scattering transform; 2) combining the halo density field with its Fourier transformed counterpart enhances predictions, while augmenting the training dataset with measured statistics further improves performance; 3) achieving high accuracy in inferring $\Omega_m$, $h$, and $\sigma_8$ by the neural network model, while being inefficient in predicting $\Omega_b$, { $n_s$}, $M_\nu$ and $w$; 4) { compared to the simple fully connected network trained with three statistical quantities, our CNN yields statistically reduced errors, showing improvements of approximately 23\% for $\Omega_m$, 11\% for $h$, 8\% for $n_s$, and 21\% for $\sigma_8$. Additionally, in comparison with the likelihood-based analysis on $P(k)$ data, our CNN provides much tighter constraints on parameters, especially on $\Omega_m$ and $\sigma_8$.} Our study emphasizes this lCNN-based novel approach in extracting large-scale structure information and estimating cosmological parameters., Comment: v2: matches the version published in PRD,17 pages,12 figures

This paper introduces ASTRA (Algorithm for Stochastic Topological RAnking), a novel method designed for the classification of galaxies into principal cosmic web structures -- voids, sheets, filaments, and knots -- especially tailored for large spectroscopic surveys. ASTRA utilizes observed data and a corresponding random catalog, generating probabilities for data and random points to belong to each cosmic web type. Extensive evaluation includes testing on datasets from N-body simulations (dark matter only and hydrodynamical) and observations, incorporating visual inspections and quantitative analyses of mass and volume filling fractions. ASTRA demonstrates proficiency in constructing void size function catalogs and is benchmarked against established cosmic web classification methods. Spatial correlations among cosmic web types are explored using 2-point correlation functions. ASTRA's approach to the cosmic web holds promise for enhancing the retrieval of non-standard clustering statistics, thereby advancing the estimation of cosmological parameters. With simplicity, speed, and applicability across various datasets, ASTRA is a valuable tool in the quantitative study of the large-scale structure of the Universe., Comment: 12 pages, 9 figures. Comments welcome

The $\beta$-skeleton approach can be conveniently utilized to construct the cosmic web based on the spatial geometry distribution of galaxies, particularly in sparse samples. This method plays a key role in establishing the three-dimensional structure of the Universe and serves as a tool for quantitatively characterizing the nature of the cosmic web. This study is the first application of $\beta$-skeleton information as weights in mark weighted correlation functions (MCFs), presenting a novel statistical measure. We have applied the $\beta$-skeleton approach to the CMASS NGC galaxy samples from SDSS BOSS DR12 in the redshift interval $0.45 \leq z \leq 0.55$. Additionally, we applied this approach to three COLA cosmological simulations with different settings ($\Omega_m=0.25, \Omega_m=0.31, \Omega_m=0.4$) for comparison. We measured three MCFs, each weighted by i) the number of neighboring galaxies around each galaxy, ii) the average distance of each galaxy from its surrounding neighbors, and iii) the reciprocal of the average distance of each galaxy from its surrounding neighbors. By comparing measurements and calculating corresponding $\chi^2$ statistics, we observe high sensitivity to the cosmological parameter $\Omega_m$ through a joint analysis of the two-point correlation and three MCFs., Comment: 16 pages,12 figues

The spectroscopic survey of China's Space Survey Telescope (CSST) is expected to obtain a huge number of slitless spectra, including more than one hundred million galaxy spectra and millions of active galactic nuclei (AGN) spectra. By making use of these spectra, we can measure the Baryon Acoustic Oscillation (BAO) signals over large redshift ranges with excellent precisions. In this work, we predict the CSST measurements of the post-reconstruction galaxy power spectra at $0

We investigate the feasibility of using COmoving Lagrangian Acceleration (COLA) technique to efficiently generate galaxy mock catalogues that can accurately reproduce the statistical properties of observed galaxies. Our proposed scheme combines the subhalo abundance matching (SHAM) procedure with COLA simulations, utilizing only three free parameters: the scatter magnitude ($\sigma_{\rm scat}$) in SHAM, the initial redshift ($z_{\rm init}$) of the COLA simulation, and the time stride ($da$) used by COLA. In this proof-of-concept study, we focus on a subset of BOSS CMASS NGC galaxies within the redshift range $z\in [0.45, 0.55]$. We perform $\mathtt{GADGET}$ simulation and low-resolution COLA simulations with various combinations of $(z_{\rm init}, da)$, each using $1024^{3}$ particles in an $800~h^{-1}{\rm Mpc}$ box. By minimizing the difference between COLA mock and CMASS NGC galaxies for the monopole of the two-point correlation function (2PCF), we obtain the optimal $\sigma_{\rm scat}$. We have found that by setting $z_{\rm init}=29$ and $da=1/30$, we achieve a good agreement between COLA mock and CMASS NGC galaxies within the range of 4 to $20~h^{-1}{\rm Mpc}$, with a computational cost two orders of magnitude lower than that of the N-body code. Moreover, a detailed verification is performed by comparing various statistical properties, such as anisotropic 2PCF, three-point clustering, and power spectrum multipoles, which shows similar performance between GADGET mock and COLA mock catalogues with the CMASS NGC galaxies. Furthermore, we assess the robustness of the COLA mock catalogues across different cosmological models, demonstrating consistent results in the resulting 2PCFs. Our findings suggest that COLA simulations are a promising tool for efficiently generating mock catalogues for emulators and machine learning analyses in exploring the large-scale structure of the Universe., Comment: 24 pages, 14 figures, 4 tables

Environmental and instrumental conditions can cause anomalies in astronomical images, which can potentially bias all kinds of measurements if not excluded. Detection of the anomalous images is usually done by human eyes, which is slow and sometimes not accurate. This is an important issue in weak lensing studies, particularly in the era of large scale galaxy surveys, in which image qualities are crucial for the success of galaxy shape measurements. In this work we present two automatic methods for detecting anomalous images in astronomical datasets. The anomalous features can be divided into two types: one is associated with the source images, and the other appears on the background. Our first method, called the Entropy Method, utilizes the randomness of the orientation distribution of the source shapes and the background gradients to quantify the likelihood of an exposure being anomalous. Our second method involves training a neural network (autoencoder) to detect anomalies. We evaluate the effectiveness of the Entropy Method on the CFHTLenS and DECaLS DR3 data. In CFHTLenS, with 1171 exposures, the Entropy Method outperforms human inspection by detecting 12 of the 13 anomalous exposures found during human inspection and uncovering 10 new ones. In DECaLS DR3, with 17112 exposures, the Entropy method detects a significant number of anomalous exposures while keeping a low false positive rate. We find that although the neural network performs relatively well in detecting source anomalies, its current performance is not as good as the Entropy Method.

The accuracy of Bayesian inference can be negatively affected by the use of inaccurate forward models. In the case of gravitational-wave inference, accurate but computationally expensive waveform models are sometimes substituted with faster but approximate ones. The model error introduced by this substitution can be mitigated in various ways, one of which is by interpolating and marginalizing over the error using Gaussian process regression. However, the use of Gaussian process regression is limited by the curse of dimensionality, which makes it less effective for analyzing higher-dimensional parameter spaces and longer signal durations. In this work, to address this limitation, we focus on gravitational-wave signals from extreme-mass-ratio inspirals as an example, and propose several significant improvements to the base method: an improved prescription for constructing the training set, GPU-accelerated training algorithms, and a new likelihood that better adapts the base method to the presence of detector noise. Our results suggest that the new method is more viable for the analysis of realistic gravitational-wave data.

The peculiar velocities of dark matter halos are crucial to study many issues in cosmology and galaxy evolution. In this study, by using the state-of-the-art deep learning technique, a UNet-based neural network, we propose to reconstruct the peculiar velocity field from the redshift-space distribution of dark matter halos. Through a point-to-point comparison and examination of various statistical properties, we demonstrate that, the reconstructed velocity field is in good agreement with the ground truth. The power spectra of various velocity field components, including velocity magnitude, divergence and vorticity, can be successfully recovered when $k\lesssim 1.1$ $h/\rm Mpc$ (the Nyquist frequency of the simulations) at about 80% accuracy. This approach is very promising and presents an alternative method to correct the redshift-space distortions using the measured 3D spatial information of halos. Additionally, for the reconstruction of the momentum field of halos, UNet achieves similar good results. Hence the applications in various aspects of cosmology are very broad, such as correcting redshift errors and improving measurements in the structure of the cosmic web, the kinetic Sunyaev-Zel'dovich effect, BAO reconstruction, etc., Comment: 15 pages, 13 figures, 4 tables

The star formation history (SFH) of galaxies is critical for understanding galaxy evolution. Hydrodynamical simulations enable us to precisely reconstruct the SFH of galaxies and establish a link to the underlying physical processes. In this work, we present a model to describe individual galaxies' SFHs from three simulations: TheThreeHundred, Illustris-1 and TNG100-1. This model divides the galaxy SFH into two distinct components: the "main sequence" and the "variation". The "main sequence" part is generated by tracing the history of the $SFR-M_*$ main sequence of galaxies across time. The "variation" part consists of the scatter around the main sequence, which is reproduced by fractional Brownian motions. We find that: 1) The evolution of the main sequence varies between simulations; 2) fractional Brownian motions can reproduce many features of SFHs, however, discrepancies still exist; 3) The variations and mass-loss rate are crucial for reconstructing the SFHs of the simulations. This model provides a fair description of the SFHs in simulations. On the other hand, by correlating the fractional Brownian motion model to simulation data, we provide a 'standard' against which to compare simulations., Comment: 22 pages, 18 figures, 2 tables

As one of Stage IV space-based telescopes, China Space Station Telescope (CSST) can perform photometric and spectroscopic surveys simultaneously to efficiently explore the Universe in extreme precision. In this work, we investigate several powerful CSST cosmological probes, including cosmic shear, galaxy-galaxy lensing, photometric and spectroscopic galaxy clustering, and number counts of galaxy clusters, and study the capability of these probes by forecasting the results of joint constraints on the cosmological parameters. By referring to real observational results, we generate mock data and estimate the measured errors based on CSST observational and instrumental designs. To study the systematical effects on the results, we also consider a number of systematics in CSST photometric and spectroscopic surveys, such as the intrinsic alignment, shear calibration uncertainties, photometric redshift uncertainties, galaxy bias, non-linear effects, instrumental effects, etc. The Fisher matrix method is used to derive the constraint results from individual or joint surveys on the cosmological and systematical parameters. We find that the joint constraints by including all these CSST cosmological probes can significantly improve the results from current observations by one order of magnitude at least, which gives $\Omega_m$ and $\sigma_8$ $<$1% accuracy, and $w_0$ and $w_a$ $<$5% and 20% accuracies, respectively. This indicates that the CSST photometric and spectroscopic multi-probe surveys could provide powerful tools to explore the Universe and greatly improve the studies of relevant cosmological problems., Comment: 17 pages, 12 figures, 3 tables. Accepted for publication in MNRAS

The mark weighted correlation function (MCF) $W(s,\mu)$ is a computationally efficient statistical measure which can probe clustering information beyond that of the conventional 2-point statistics. In this work, we extend the traditional mark weighted statistics by using powers of the density field gradient $|\nabla \rho/\rho|^\alpha$ as the weight, and use the angular dependence of the scale-averaged MCFs to constrain cosmological parameters. The analysis shows that the gradient based weighting scheme is statistically more powerful than the density based weighting scheme, while combining the two schemes together is more powerful than separately using either of them. Utilising the density weighted or the gradient weighted MCFs with $\alpha=0.5,\ 1$, we can strengthen the constraint on $\Omega_m$ by factors of 2 or 4, respectively, compared with the standard 2-point correlation function, while simultaneously using the MCFs of the two weighting schemes together can be $1.25$ times more statistically powerful than using the gradient weighting scheme alone. The mark weighted statistics may play an important role in cosmological analysis of future large-scale surveys. Many issues, including the possibility of using other types of weights, the influence of the bias on this statistics, as well as the usage of MCFs in the tomographic Alcock-Paczynski method, are worth further investigations., Comment: 11 pages, 8 figures. Accepted by MNRAS

We investigate impacts of massive neutrinos on the cosmic velocity fields, employing high-resolution cosmological $N$-body simulations provided by the information-optimized CUBE code, where cosmic neutrinos are evolved using collisionless hydrodynamics and their perturbations can be accurately resolved. In this study we focus, for the first time, on the analysis of massive-neutrino induced suppression effects in various cosmic velocity field components of velocity magnitude, divergence, vorticity and dispersion. By varying the neutrino mass sum $M_\nu$ from 0 -- 0.4 eV, the simulations show that, the power spectra of vorticity -- exclusively sourced by non-linear structure formation that is affected by massive neutrinos significantly -- is very sensitive to the mass sum, which potentially provide novel signatures in detecting massive neutrinos. Furthermore, using the chi-square statistic, we quantitatively test the sensitivity of the density and velocity power spectra to the neutrino mass sum. Indeed, we find that, the vorticity spectrum has the highest sensitivity, and the null hypothesis of massless neutrinos is incompatible with both vorticity and divergence spectra from $M_\nu=0.1$ eV at high significance ($p$-value $= 0.03$ and $0.07$, respectively). These results demonstrate clearly the importance of peculiar velocity field measurements, in particular of vorticity and divergence components, in determination of neutrino mass and mass hierarchy., Comment: 13 pages, 9 figures. Monthly Notices of the Royal Astronomical Society, 2022

Precise cosmic web classification of observed galaxies in massive spectroscopic surveys can be either highly uncertain or computationally expensive. As an alternative, we explore a fast Machine Learning-based approach to infer the underlying dark matter tidal cosmic web environment of a galaxy distribution from its $\beta$-skeleton graph. We develop and test our methodology using the cosmological magnetohydrodynamic simulation Illustris-TNG at $z=0$. We explore three different tree-based machine-learning algorithms to find that a random forest classifier can best use graph-based features to classify a galaxy as belonging to a peak, filament or sheet as defined by the T-Web classification algorithm. The best match between the galaxies and the dark matter T-Web corresponds to a density field smoothed over scales of $2$ Mpc, a threshold over the eigenvalues of the dimensionless tidal tensor of $\lambda_{\rm{th}}=0.0$ and galaxy number densities around $8\times 10^{-3}$ Mpc$^{-3}$. This methodology results on a weighted F1 score of 0.728 and a global accuracy of 74\%. More extensive tests that take into account lightcone effects and redshift space distortions (RSD) are left for future work. We make one of our highest ranking random forest models available on a public repository for future reference and reuse., Comment: ApJ accepted, 13 pages, 9 figures

We develop a deep learning technique to infer the non-linear velocity field from the dark matter density field. The deep learning architecture we use is an "U-net" style convolutional neural network, which consists of 15 convolution layers and 2 deconvolution layers. This setup maps the 3-dimensional density field of $32^3$-voxels to the 3-dimensional velocity or momentum fields of $20^3$-voxels. Through the analysis of the dark matter simulation with a resolution of $2 {h^{-1}}{\rm Mpc}$, we find that the network can predict the the non-linearity, complexity and vorticity of the velocity and momentum fields, as well as the power spectra of their value, divergence and vorticity and its prediction accuracy reaches the range of $k\simeq1.4$ $h{\rm Mpc}^{-1}$ with a relative error ranging from 1% to $\lesssim$10%. A simple comparison shows that neural networks may have an overwhelming advantage over perturbation theory in the reconstruction of velocity or momentum fields., Comment: 10 pages, 6 figures, 4 tables, accepted for publication in ApJ

Space-borne gravitational wave detectors like TianQin are expected to detect GW signals emitted by the mergers of massive black hole binaries. Luminosity distance information can be obtained from GW observations, and one can perform cosmological inference if redshift information can also be extracted, which would be straightforward if an electromagnetic counterpart exists. In this paper, we concentrate on the conservative scenario where the EM counterparts are not available, and comprehensively study if cosmological parameters can be inferred through a statistical approach, utilizing the non-uniform distribution of galaxies as well as the black hole mass-host galaxy bulge luminosity relationship. By adopting different massive black hole binary merger models, and assuming different detector configurations, we conclude that the statistical inference of cosmological parameters is indeed possible. TianQin is expected to constrain the Hubble constant to a relative error of about 4%-7%, depending on the underlying model. The multidetector network of TianQin and LISA can significantly improve the precision of cosmological parameters. In the most favorable model, it is possible to achieve a level of 1.7% with a network of TianQin and LISA. We find that without EM counterparts, constraints on all other parameters need a larger number of events or more precise sky localization of GW sources, which can be achieved by the multidetector network or under a favorable model for massive black hole mergers. However, in the optimistic case, where EM counterparts are available, one can obtain useful constraints on all cosmological parameters in the Lambda-CDM cosmology, regardless of the population model. Moreover, we can also constrain the equation of state of the dark energy without the EM counterparts, and it is even possible to study the evolution of EoS of the DE when the EM counterparts are observed., Comment: 34 pages, 17 figures, update to match published version

We study the role of gravitational waves (GW) in the heat death of the universe. Due to the GW emission, in a very long period, dynamical systems in the universe suffer from persistent mechanical energy dissipation, evolving to a state of universal rest and death. With N-body simulations, we adopt a simple yet representative scheme to calculate the energy loss due to the GW emission. For current dark matter systems with mass $\sim10^{12}-10^{15} M_\odot$, we estimate their GW emission timescale as $\sim10^{19}-10^{25}$ years. This timescale is significantly larger than any baryon processes in the universe, but still $\sim10^{80}$ times shorter than that of the Hawking radiation. We stress that our analysis could be invalid due to many unknowns such as the dynamical chaos, the quadrupole momentum of halos, the angular momentum loss, the dynamic friction, the central black hole accretion, the dark matter decays or annihilations, the property of dark energy and the future evolution of the universe., Comment: Accepted for publication in MNRAS. 6 pages, 4 figures

We explore the information theory entropy of a graph as a scalar to quantify the cosmic web. We find entropy values in the range between 1.5 and 3.2 bits. We argue that this entropy can be used as a discrete analogue of scalars used to quantify the connectivity in continuous density fields. After showing that the entropy clearly distinguishes between clustered and random points, we use simulations to gauge the influence of survey geometry, cosmic variance, redshift space distortions, redshift evolution, cosmological parameters and spatial number density. Cosmic variance shows the least important influence while changes from the survey geometry, redshift space distortions, cosmological parameters and redshift evolution produce larger changes on the order of $10^{-2}$ bits. The largest influence on the graph entropy comes from changes in the number density of clustered points. As the number density decreases, and the cosmic web is less pronounced, the entropy can diminish up to 0.2 bits. The graph entropy is simple to compute and can be applied both to simulations and observational data from large galaxy redshift surveys; it is a new statistic that can be used in a complementary way to other kinds of topological or clustering measurements., Comment: 5 pages, 4 figures, Accepted for publication in MNRAS

We used the mark weighted correlation functions (MCFs), $W(s)$, to study the large scale structure of the Universe. We studied five types of MCFs with the weighting scheme $\rho^\alpha$, where $\rho$ is the local density, and $\alpha$ is taken as $-1,\ -0.5,\ 0,\ 0.5$, and 1. We found that different MCFs have very different amplitudes and scale-dependence. Some of the MCFs exhibit distinctive peaks and valleys that do not exist in the standard correlation functions. Their locations are robust against the redshifts and the background geometry, however it is unlikely that they can be used as ``standard rulers'' to probe the cosmic expansion history. Nonetheless we find that these features may be used to probe parameters related with the structure formation history, such as the values of $\sigma_8$ and the galaxy bias. Finally, after conducting a comprehensive analysis using the full shapes of the $W(s)$s and $W_{\Delta s}(\mu)$s, we found that, combining different types of MCFs can significantly improve the cosmological parameter constraints. Compared with using only the standard correlation function, the combinations of MCFs with $\alpha=0,\ 0.5,\ 1$ and $\alpha=0,\ -1,\ -0.5,\ 0.5,\ 1$ can improve the constraints on $\Omega_m$ and $w$ by $\approx30\%$ and $50\%$, respectively. We find highly significant evidence that MCFs can improve cosmological parameter constraints., Comment: 15pages, 17figures, APJ accepted