Accelerated Bayesian SED Modeling using Amortized Neural Posterior Estimation

Abstract

State-of-the-art spectral energy distribution (SED) analyses use a Bayesian framework to infer the physical properties of galaxies from observed photometry or spectra. They require sampling from a high-dimensional space of SED model parameters and take $>10-100$ CPU hours per galaxy, which renders them practically infeasible for analyzing the $billions$ of galaxies that will be observed by upcoming galaxy surveys ($e.g.$ DESI, PFS, Rubin, Webb, and Roman). In this work, we present an alternative scalable approach to rigorous Bayesian inference using Amortized Neural Posterior Estimation (ANPE). ANPE is a simulation-based inference method that employs neural networks to estimate the posterior probability distribution over the full range of observations. Once trained, it requires no additional model evaluations to estimate the posterior. We present, and publicly release, $\mathrm{SED}flow$, an ANPE method to produce posteriors of the recent Hahn et al. (2022) SED model from optical photometry. $\mathrm{SED}flow$ takes $\sim 1$ $secondpergalaxy$ to obtain the posterior distributions of 12 model parameters, all of which are in excellent agreement with traditional Markov Chain Monte Carlo sampling results. We also apply $\mathrm{SED}flow$ to 33,884 galaxies in the NASA-Sloan Atlas and publicly release their posteriors: see https://changhoonhahn.github.io/SEDflow.

arXiv comment21 pages, 5 figures; submitted to ApJ; code available at https://changhoonhahn.github.io/SEDflow

References in corpus (61)

• Maps of Dust IR Emission for Use in Estimation of Reddening and CMBR Foregrounds
• Stellar population synthesis at the resolution of 2003
• The Sloan Digital Sky Survey: Technical Summary
• Galactic Stellar and Substellar Initial Mass Function
• Modules for Experiments in Stellar Astrophysics (MESA)
• Modules for Experiments in Stellar Astrophysics (MESA): Giant Planets, Oscillations, Rotation, and Massive Stars
• LSST: from Science Drivers to Reference Design and Anticipated Data Products
• Modules for Experiments in Stellar Astrophysics (MESA): Binaries, Pulsations, and Explosions
• MESA Isochrones and Stellar Tracks (MIST). I: Solar-Scaled Models
• Introducing the Illustris Project: Simulating the coevolution of dark and visible matter in the Universe
• The Cosmic Evolution Survey (COSMOS) -- Overview
• The James Webb Space Telescope
• MESA Isochrones and Stellar Tracks (MIST) 0: Methods for the construction of stellar isochrones
• Evolutionary population synthesis: models, analysis of the ingredients and application to high-z galaxies
• 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
• MILES: A Medium resolution INT Library of Empirical Spectra
• The Eighth Data Release of the Sloan Digital Sky Survey: First Data from SDSS-III
• The propagation of uncertainties in stellar population synthesis modeling III: model calibration, comparison, and evaluation
• Introducing the Illustris Project: the evolution of galaxy populations across cosmic time
• The frontier of simulation-based inference
• 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
• Extragalactic Science, Cosmology and Galactic Archaeology with the Subaru Prime Focus Spectrograph (PFS)
• Stellar Population Inference with Prospector
• A standard stellar library for evolutionary synthesis: II. The M dwarf extension
• Improved background subtraction for the Sloan Digital Sky Survey images
• How to Measure Galaxy Star Formation Histories II: Nonparametric Models
• Science Objectives and Early Results of the DEEP2 Redshift Survey
• The Illustris Simulation: Public Data Release
• Deriving Physical Properties from Broadband Photometry with Prospector: Description of the Model and a Demonstration of its Accuracy Using 129 Galaxies in the Local Universe
• A Substantial Population of Low Mass Stars in Luminous Elliptical Galaxies
• Modelling and interpreting spectral energy distributions of galaxies with BEAGLE
• The initial mass function of early-type galaxies
• The Dust Attenuation Law in Distant Galaxies: Evidence for Variation with Spectral Type
• Fitting the integrated Spectral Energy Distributions of Galaxies
• How to measure galaxy star-formation histories I: Parametric models
• Real-time gravitational-wave science with neural posterior estimation
• 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
• Fast likelihood-free cosmology with neural density estimators and active learning
• Fast, Slow, Early, Late: Quenching Massive Galaxies at z~0.8
• Mining gold from implicit models to improve likelihood-free inference
• Likelihood-free inference with neural compression of DES SV weak lensing map statistics
• 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
• Approximate Bayesian Computation for Astronomical Model Analysis: A Case Study in Galaxy Demographics and Morphological Transformation at High Redshift
• Likelihood-Free Cosmological Inference with Type Ia Supernovae: Approximate Bayesian Computation for a Complete Treatment of Uncertainty
• SED fitting with Markov Chain Monte Carlo: Methodology and Application to z=3.1 Lyman Alpha Emitting Galaxies
• Dimensionality Reduction of SDSS Spectra with Variational Autoencoders
• SPECULATOR: Emulating stellar population synthesis for fast and accurate galaxy spectra and photometry
• Chemical tagging can work: Identification of stellar phase-space structures purely by chemical-abundance similarity
• Survey of gravitationally-lensed objects in HSC imaging (SuGOHI). III. Statistical strong lensing constraints on the stellar IMF of CMASS galaxies
• Approximate Bayesian Computation in Large Scale Structure: constraining the galaxy-halo connection
• Gravitational wave population inference with deep flow-based generative network
• Likelihood Non-Gaussianity in Large-Scale Structure Analyses
• Recovering stellar population parameters via different population models and stellar libraries
• Real-Time Likelihood-Free Inference of Roman Binary Microlensing Events with Amortized Neural Posterior Estimation
• Accelerating Approximate Bayesian Computation with Quantile Regression: Application to Cosmological Redshift Distributions
• The Influence of Galaxy Environment on the Stellar Initial Mass Function of Early-Type Galaxies
• Accurate X-ray Timing in the Presence of Systematic Biases With Simulation-Based Inference

Cited by in corpus (4)

• $\mathrm{S}\mathrm{IM}\mathrm{BIG}$: Mock Challenge for a Forward Modeling Approach to Galaxy Clustering
• Hierarchical Bayesian inference of photometric redshifts with stellar population synthesis models
• Stochastic Modelling of Star Formation Histories III. Constraints from Physically-Motivated Gaussian Processes
• DSPS: Differentiable Stellar Population Synthesis