skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Inferring Cosmological Parameters on SDSS via Domain-generalized Neural Networks and Light-cone Simulations
Abstract We present a proof-of-concept simulation-based inference on Ωmandσ8from the Sloan Digital Sky Survey (SDSS) Baryon Oscillation Spectroscopic Survey (BOSS) LOWZ Northern Galactic Cap (NGC) catalog using neural networks and domain generalization techniques without the need of summary statistics. Using rapid light-cone simulations L-picola, mock galaxy catalogs are produced that fully incorporate the observational effects. The collection of galaxies is fed as input to a point cloud-based network,Minkowski-PointNet. We also add relatively more accurate Gadgetmocks to obtain robust and generalizable neural networks. By explicitly learning the representations that reduce the discrepancies between the two different data sets via the semantic alignment loss term, we show that the latent space configuration aligns into a single plane in which the two cosmological parameters form clear axes. Consequently, during inference, the SDSS BOSS LOWZ NGC catalog maps onto the plane, demonstrating effective generalization and improving prediction accuracy compared to non-generalized models. Results from the ensemble of 25 independently trained machines find Ωm= 0.339 ± 0.056 andσ8= 0.801 ± 0.061, inferred only from the distribution of galaxies in the light-cone slices without relying on any indirect summary statistics. A single machine that best adapts to the Gadgetmocks yields a tighter prediction of Ωm= 0.282 ± 0.014 andσ8= 0.786 ± 0.036. We emphasize that adaptation across multiple domains can enhance the robustness of the neural networks in observational data.  more » « less
Award ID(s):
2406729
PAR ID:
10612803
Author(s) / Creator(s):
; ; ; ; ; ; ; ;
Publisher / Repository:
The Astrophysical Journal
Date Published:
Journal Name:
The Astrophysical Journal
Volume:
975
Issue:
1
ISSN:
0004-637X
Page Range / eLocation ID:
38
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract We analyze clustering measurements of BOSS galaxies using a simulation-based emulator of two-point statistics. We focus on the monopole and quadrupole of the redshift-space correlation function, and the projected correlation function, at scales of 0.1 ∼ 60h−1Mpc. Although our simulations are based onwCDM with general relativity (GR), we include a scaling parameter of the halo velocity field,γf, defined as the amplitude of the halo velocity field relative to the GR prediction. We divide the BOSS data into three redshift bins. After marginalizing over other cosmological parameters, galaxy bias parameters, and the velocity scaling parameter, we findfσ8(z= 0.25) = 0.413 ± 0.031,fσ8(z= 0.4) = 0.470 ± 0.026, andfσ8(z= 0.55) = 0.396 ± 0.022. Compared with Planck observations using a flat Lambda cold dark matter model, our results are lower by 1.9σ, 0.3σ, and 3.4σ, respectively. These results are consistent with other recent simulation-based results at nonlinear scales, including weak lensing measurements of BOSS LOWZ galaxies, two-point clustering of eBOSS LRGs, and an independent clustering analysis of BOSS LOWZ. All these results are generally consistent with a combination of γ f 1 / 2 σ 8 0.75 . We note, however, that the BOSS data is well fit assuming GR, i.e.,γf= 1. We cannot rule out an unknown systematic error in the galaxy bias model at nonlinear scales, but near-future data and modeling will enhance our understanding of the galaxy–halo connection, and provide a strong test of new physics beyond the standard model. 
    more » « less
  2. Abstract Recent works have discovered a relatively tight correlation between Ωmand the properties of individual simulated galaxies. Because of this, it has been shown that constraints on Ωmcan be placed using the properties of individual galaxies while accounting for uncertainties in astrophysical processes such as feedback from supernovae and active galactic nuclei. In this work, we quantify whether using the properties of multiple galaxies simultaneously can tighten those constraints. For this, we train neural networks to perform likelihood-free inference on the value of two cosmological parameters (Ωmandσ8) and four astrophysical parameters using the properties of several galaxies from thousands of hydrodynamic simulations of the CAMELS project. We find that using properties of more than one galaxy increases the precision of the Ωminference. Furthermore, using multiple galaxies enables the inference of other parameters that were poorly constrained with one single galaxy. We show that the same subset of galaxy properties are responsible for the constraints on Ωmfrom one and multiple galaxies. Finally, we quantify the robustness of the model and find that without identifying the model range of validity, the model does not perform well when tested on galaxies from other galaxy formation models. 
    more » « less
  3. Abstract The cluster mass–richness relation (MRR) is an observationally efficient and potentially powerful cosmological tool for constraining the matter density Ωmand the amplitude of fluctuationsσ8using the cluster abundance technique. We derive the MRR relation usingGalWCat19, a publicly available galaxy cluster catalog we created from the Sloan Digital Sky Survey-DR13 spectroscopic data set. In the MRR, cluster mass scales with richness as log M 200 = α + β log N 200 . We find that the MRR we derive is consistent with both the IllustrisTNG and mini-Uchuu cosmological numerical simulations, with a slope ofβ≈ 1. We use the MRR we derived to estimate cluster masses from theGalWCat19catalog, which we then use to set constraints on Ωmandσ8. Utilizing the all-member MRR, we obtain constraints of Ωm= 0.31 0.03 + 0.04 andσ8= 0.82 0.04 + 0.05 , and utilizing the red member MRR only, we obtain Ωm= 0.31 0.03 + 0.04 andσ8= 0.81 0.04 + 0.05 . Our constraints on Ωmandσ8are consistent and very competitive with the Planck 2018 results. 
    more » « less
  4. Abstract Observed evolution of the total mass distribution with redshift is crucial to testing galaxy evolution theories. To measure the total mass distribution, strong gravitational lenses complement the resolved dynamical observations that are currently limited toz≲ 0.5. Here we present the lens models for a pilot sample of seven galaxy-scale lenses from theASTRO3DGalaxy Evolution with Lenses (AGEL) survey. TheAGELlenses, modeled using HST/WFC3-F140W images with Gravitational Lens Efficient Explorer (GLEE) software, have deflector redshifts in the range 0.3 <zdefl< 0.9. Assuming a power-law density profile with slopeγ, we measure the total density profile for the deflector galaxies via lens modeling. We also measure the stellar velocity dispersions (σobs) for four lenses and obtainσobsfromSDSS-BOSSfor the remaining lenses to test our lens models by comparing observed and model-predicted velocity dispersions. For the sevenAGELlenses, we measure an average density profile slope of −1.95 ± 0.09 and aγ–zrelation that does not evolve with redshift atz< 1. Although our result is consistent with some observations and simulations, it differs from other studies atz< 1 that suggest theγ–zrelation evolves with redshift. The apparent conflicts among observations and simulations may be due to a combination of (1) systematics in the lensing and dynamical modeling; (2) challenges in comparing observations with simulations; and (3) assuming a simple power law for the total mass distribution. By providing more lenses atzdefl> 0.5, theAGELsurvey will provide stronger constraints on whether the mass profiles evolve with redshift as predicted by current theoretical models. 
    more » « less
  5. Abstract The rapid advancement of large-scale cosmological simulations has opened new avenues for cosmological and astrophysical research. However, the increasing diversity among cosmological simulation models presents a challenge to therobustness. In this work, we develop the Model-Insensitive ESTimator (Miest), a machine that canrobustlyestimate the cosmological parameters, Ωmandσ8, from neural hydrogen maps of simulation models in the Cosmology and Astrophysics with MachinE Learning Simulations project—IllustrisTNG,SIMBA, Astrid, and SWIFT-Eagle. An estimator is consideredrobustif it possesses a consistent predictive power across all simulations, including those used during the training phase. We train our machine using multiple simulation models and ensure that it only extracts common features between the models while disregarding the model-specific features. This allows us to develop a novel model that is capable of accurately estimating parameters across a range of simulation models, without being biased toward any particular model. Upon the investigation of the latent space—a set of summary statistics, we find that the implementation ofrobustnessleads to the blending of latent variables across different models, demonstrating the removal of model-specific features. In comparison to a standard machine lackingrobustness, the average performance of Mieston the unseen simulations during the training phase has been improved by ∼17% for Ωmand 38% forσ8. By using a machine learning approach that can extractrobust, yet physical features, we hope to improve our understanding of galaxy formation and evolution in a (subgrid) model-insensitive manner, and ultimately, gain insight into the underlying physical processes responsible forrobustness. 
    more » « less