PopSED: Population-Level Inference for Galaxy Properties from Broadband Photometry with Neural Density Estimation

Abstract

We present PopSED, a framework for the population-level inference of galaxy properties from photometric data. Unlike the traditional approach of first analyzing individual galaxies and then combining the results to determine the physical properties of the entire galaxy population, we directly make the population distribution the inference objective. We train normalizing flows to approximate the population distribution by minimizing the Wasserstein distance between the synthetic photometry of the galaxy population and the observed data. We validate our method using mock observations and apply it to galaxies from the GAMA survey. PopSED reliably recovers the redshift and stellar mass distribution of $10^5$ galaxies using broadband photometry within $<1$ GPU hr, being $10^{5-6}$ times faster than the traditional spectral energy distribution modeling method. From the population posterior, we also recover the star-forming main sequence for GAMA galaxies at $z<0.1$. With the unprecedented number of galaxies in upcoming surveys, our method offers an efficient tool for studying galaxy evolution and deriving redshift distributions for cosmological analyses.

arXiv comment13 pages, 4 figures. AJ accepted

References in corpus (67)

• Array Programming with NumPy
• Astropy: A Community Python Package for Astronomy
• Planck 2015 results. XIII. Cosmological parameters
• Stellar population synthesis at the resolution of 2003
• Galactic Stellar and Substellar Initial Mass Function
• The Astropy Project: Building an inclusive, open-science project and status of the v2.0 core package
• The Dust Content and Opacity of Actively Star-Forming Galaxies
• The Seventh Data Release of the Sloan Digital Sky Survey
• The Astropy Project: Sustaining and Growing a Community-oriented Open-source Project and the Latest Major Release (v5.0) of the Core Package
• The Origin of the Mass--Metallicity Relation: Insights from 53,000 Star-Forming Galaxies in the SDSS
• LSST: from Science Drivers to Reference Design and Anticipated Data Products
• The physical properties of star forming galaxies in the low redshift universe
• MESA Isochrones and Stellar Tracks (MIST). I: Solar-Scaled Models
• The 2.5 m Telescope of the Sloan Digital Sky Survey
• Introducing the Illustris Project: Simulating the coevolution of dark and visible matter in the Universe
• MESA Isochrones and Stellar Tracks (MIST) 0: Methods for the construction of stellar isochrones
• The propagation of uncertainties in stellar population synthesis modeling I: The relevance of uncertain aspects of stellar evolution and the IMF to the derived physical properties of galaxies
• A Simple Model for the Absorption of Starlight by Dust in Galaxies
• Overview of the SDSS-IV MaNGA Survey: Mapping Nearby Galaxies at Apache Point Observatory
• MILES: A Medium resolution INT Library of Empirical Spectra
• A Highly Consistent Framework for the Evolution of the Star-Forming "Main Sequence" from z~0-6
• The propagation of uncertainties in stellar population synthesis modeling III: model calibration, comparison, and evaluation
• Galaxy and Mass Assembly (GAMA): survey diagnostics and core data release
• Introducing the Illustris Project: the evolution of galaxy populations across cosmic time
• Analysis of galaxy SEDs from far-UV to far-IR with CIGALE: Studying a SINGS test sample
• Modeling the Panchromatic Spectral Energy Distributions of Galaxies
• Inferring the star-formation histories of massive quiescent galaxies with BAGPIPES: Evidence for multiple quenching mechanisms
• A standard stellar library for evolutionary synthesis: I. Calibration of theoretical spectra
• Stellar Population Inference with Prospector
• Galaxy And Mass Assembly: Stellar Mass Estimates
• A standard stellar library for evolutionary synthesis: II. The M dwarf extension
• A population of red candidate massive galaxies ~600 Myr after the Big Bang
• How to Measure Galaxy Star Formation Histories II: Nonparametric Models
• The Illustris Simulation: Public Data Release
• A Modified Magnitude System that Produces Well-Behaved Magnitudes, Colors, and Errors Even for Low Signal-to-Noise Ratio Measurements
• An Objective Definition for the Main Sequence of Star-Forming Galaxies
• Weak lensing for precision cosmology
• The Mass-Metallicity and the Fundamental Metallicity Relation revisited on a fully Te-based abundance scale for galaxies
• Photometric Response Functions of the SDSS Imager
• How to measure galaxy star-formation histories I: Parametric models
• DESI Bright Galaxy Survey: Final Target Selection, Design, and Validation
• Galaxy And Mass Assembly (GAMA): the G02 field, Herschel-ATLAS target selection and Data Release 3
• A standard stellar library for evolutionary synthesis. III. Metallicity calibration
• An Older, More Quiescent Universe from Panchromatic SED Fitting of the 3D-HST Survey
• Fast likelihood-free cosmology with neural density estimators and active learning
• Galaxy And Mass Assembly (GAMA): Panchromatic Data Release (far-UV --- far-IR) and the low-z energy budget
• Extreme deconvolution: Inferring complete distribution functions from noisy, heterogeneous and incomplete observations
• Parameters of 220 million stars from Gaia BP/RP spectra
• Galaxy And Mass Assembly (GAMA): The galaxy stellar mass function to $z=0.1$ from the r-band selected equatorial regions
• Deep Extragalactic VIsible Legacy Survey (DEVILS): SED Fitting in the D10-COSMOS Field and the Evolution of the Stellar Mass Function and SFR-$M_{\star }$ relation
• Photometric Redshifts for Next-Generation Surveys
• The Euclid mission design
• SDSS IV MaNGA - An archaeological view of the Cosmic Star Formation History
• SPECULATOR: Emulating stellar population synthesis for fast and accurate galaxy spectra and photometry
• The Star Forming Main Sequence of Dwarf Low Surface Brightness Galaxies
• Hierarchical Bayesian inference of galaxy redshift distributions from photometric surveys
• Translation and Rotation Equivariant Normalizing Flow (TRENF) for Optimal Cosmological Analysis
• Gravitational wave population inference with deep flow-based generative network
• SBI++: Flexible, Ultra-fast Likelihood-free Inference Customized for Astronomical Applications
• Forward modeling of galaxy populations for cosmological redshift distribution inference
• Real-Time Likelihood-Free Inference of Roman Binary Microlensing Events with Amortized Neural Posterior Estimation
• The DESI PRObabilistic Value-Added Bright Galaxy Survey (PROVABGS) Mock Challenge
• Autoencoding Galaxy Spectra II: Redshift Invariance and Outlier Detection
• Quasar Factor Analysis -- An Unsupervised and Probabilistic Quasar Continuum Prediction Algorithm with Latent Factor Analysis
• Deep Potential: Recovering the gravitational potential from a snapshot of phase space
• DIGS: Deep Inference of Galaxy Spectra with Neural Posterior Estimation
• Wasserstein distance as a new tool for discriminating cosmologies through the topology of large scale structure

Cited by in corpus (12)

• Quantifying the Effects of Known Unknowns on Inferred High-redshift Galaxy Properties: Burstiness, the IMF, and Nebular Physics
• pop-cosmos: A comprehensive picture of the galaxy population from COSMOS data
• pop-cosmos: Scaleable inference of galaxy properties and redshifts with a data-driven population model
• Simulation-based inference of galaxy properties from JWST pixels
• Modelling the selection of galaxy groups with end to end simulations
• pop-cosmos: Insights from generative modeling of a deep, infrared-selected galaxy population
• Data-Space Validation of High-Dimensional Models by Comparing Sample Quantiles
• galsbi: A Python package for the GalSBI galaxy population model
• Flexible Simulation Based Inference for Galaxy Photometric Fitting with Synthesizer
• GalSBI: Phenomenological galaxy population model for cosmology using simulation-based inference
• GalSBI-SPS: a stellar population synthesis-based galaxy population model for cosmology and galaxy evolution applications
• Can AI Dream of Unseen Galaxies? Conditional Diffusion Model for Galaxy Morphology Augmentation