Simulation-based inference of galaxy properties from JWST pixels
arXiv:2506.04336 · doi:10.1051/0004-6361/202555810
Abstract
We present an efficient Bayesian SED-fitting framework tailored to multiwavelength pixel photometry from the JWST Advanced Deep Extragalactic Survey (JADES). Our method employs simulation-based inference to enable rapid posterior sampling across galaxy pixels, leveraging the unprecedented spatial resolution, wavelength coverage, and depth provided by the survey. It is trained on synthetic photometry generated from MILES stellar population models, incorporating both parametric and non-parametric SFHs, realistic noise, and JADES-like filter sensitivity thresholds. We validate this amortised inference approach on mock datasets, achieving robust and well-calibrated posterior distributions, with an score of 0.99 for stellar mass. Applying our pipeline to real observations, we derive spatially resolved maps of stellar population properties down to (averaged over F277W, F356W, F444W) for 1083 JADES galaxies and ~2 million pixels with spectroscopic redshifts. These maps enable the identification of dusty or starburst regions and offer insights into mass growth and the structural assembly. We assess the outshining phenomenon by comparing pixel-based and integrated stellar mass estimates, finding limited impact only in low-mass galaxies () but systematic differences of ~0.20 dex linked to SFH priors. With an average posterior sampling speed of seconds per pixel and a total inference time of ~1 CPU-day for the full dataset, our model offers a scalable solution for extracting high-fidelity stellar population properties from HST+JWST datasets, opening the way for statistical studies at sub-galactic scales.
13 pages, 7 figures. Accepted for publication in Astronomy & Astrophysics (A&A) after minor revisions and an extended discussion. Software available at https://github.com/patriglesias/SBIPIX.git
References in corpus (66)
- Stellar population synthesis at the resolution of 2003
- Galactic Stellar and Substellar Initial Mass Function
- MESA Isochrones and Stellar Tracks (MIST). I: Solar-Scaled Models
- dynesty: A Dynamic Nested Sampling Package for Estimating Bayesian Posteriors and Evidences
- CANDELS: The Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey
- Infrared Emission from Interstellar Dust. IV. The Silicate-Graphite-PAH Model in the Post-Spitzer Era
- CANDELS: The Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey - The Hubble Space Telescope Observations, Imaging Data Products and Mosaics
- 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
- The Great Observatories Origins Deep Survey: Initial Results From Optical and Near-Infrared Imaging
- Adaptive spatial binning of integral-field spectroscopic data using Voronoi tessellations
- CIGALE: a python Code Investigating GALaxy Emission
- The propagation of uncertainties in stellar population synthesis modeling III: model calibration, comparison, and evaluation
- Modeling the Panchromatic Spectral Energy Distributions of Galaxies
- The Hubble Ultra Deep Field
- Inferring the star-formation histories of massive quiescent galaxies with BAGPIPES: Evidence for multiple quenching mechanisms
- Stellar Population Inference with Prospector
- The stellar mass assembly of galaxies from z=0 to z=4. Analysis of a sample selected in the rest-frame near-infrared with Spitzer
- How to Measure Galaxy Star Formation Histories II: Nonparametric Models
- Resolved stellar mass maps of galaxies. I: method and implications for global mass estimates
- The Abundance of Star-Forming Galaxies in the Redshift Range 8.5 to 12: New Results from the 2012 Hubble Ultra Deep Field Campaign
- Preparing Red-Green-Blue (RGB) Images from CCD Data
- NIRCam Performance on JWST In Flight
- Smooth(er) Stellar Mass Maps in CANDELS: Constraints on the Longevity of Clumps in High-redshift Star-forming Galaxies
- The MUSE Hubble Ultra Deep Field Survey: I. Survey description, data reduction and source detection
- The HST eXtreme Deep Field XDF: Combining all ACS and WFC3/IR Data on the HUDF Region into the Deepest Field Ever
- How to measure galaxy star-formation histories I: Parametric models
- Recovering galaxy star formation and metallicity histories from spectra using VESPA
- Multi-Wavelength View of Kiloparsec-Scale Clumps in Star-Forming Galaxies at z~2
- The Art of Measuring Physical Parameters in Galaxies: A Critical Assessment of Spectral Energy Distribution Fitting Techniques
- JADES NIRSpec Initial Data Release for the Hubble Ultra Deep Field: Redshifts and Line Fluxes of Distant Galaxies from the Deepest JWST Cycle 1 NIRSpec Multi-Object Spectroscopy
- How Well Can We Measure the Stellar Mass of a Galaxy: The Impact of the Assumed Star Formation History Model in SED Fitting
- Fast, Slow, Early, Late: Quenching Massive Galaxies at z~0.8
- The 2012 Hubble Ultra Deep Field (UDF12): Observational Overview
- Non-parametric Star Formation History Reconstruction with Gaussian Processes I: Counting Major Episodes of Star Formation
- CEERS Key Paper V: A triality on the nature of HST-dark galaxies
- Relative merits of different types of rest-frame optical observations to constrain galaxy physical parameters
- Stochastic modelling of star-formation histories II: star-formation variability from molecular clouds and gas inflow
- Reconstruction of Galaxy Star Formation Histories through SED Fitting: The Dense Basis Approach
- The MUSE Hubble Ultra Deep Field surveys: Data release II
- Stochastic modeling of star-formation histories I: the scatter of the star-forming main sequence
- A Detailed Study of Photometric Redshifts for GOODS-South Galaxies
- Spatially Resolved Properties of High Redshift Galaxies in the SMACS0723 JWST ERO Field
- Spatially-unresolved SED fitting can underestimate galaxy masses: a solution to the missing mass problem
- On biases in the predictions of stellar population synthesis models
- Missing Stellar Mass in SED Fitting: Spatially Unresolved Photometry can Underestimate Galaxy Masses
- Stars in the Hubble Ultra Deep Field
- Understanding the Scatter in the Spatially-resolved Star Formation Main Sequence of Local Massive Spiral Galaxies
- Introducing piXedfit -- a Spectral Energy Distribution Fitting Code Designed for Resolved Sources
- Circumstellar Dust Around AGB Stars and Implications for Infrared Emission from Galaxies
- Stellar populations in local star-forming galaxies. I.-Data and modelling procedure
- Panchromatic SED modelling of spatially-resolved galaxies
- Outshining in the Spatially Resolved Analysis of a Strongly-Lensed Galaxy at z=6.072 with JWST NIRCam
- Star formation histories from multi-band photometry: A new approach
- The DESI PRObabilistic Value-Added Bright Galaxy Survey (PROVABGS) Mock Challenge
- Stochastic Modelling of Star Formation Histories III. Constraints from Physically-Motivated Gaussian Processes
- The Structures of Distant Galaxies V: The Evolution of Galaxy Structure in Stellar Mass at z < 1
- PopSED: Population-Level Inference for Galaxy Properties from Broadband Photometry with Neural Density Estimation
- Deriving the star formation histories of galaxies from spectra with simulation-based inference
- More is better: Strong constraints on the stellar properties of LEGA-C z ~ 1 galaxies with Prospector
- JWST PRIMER: A lack of outshining in four normal z =4-6 galaxies from the ALMA-CRISTAL Survey
- Solution to the conflict between the resolved and unresolved galaxy stellar mass estimation from the perspective of JWST
- Population Models for Star Formation Timescales in Early Galaxies: The First Step Towards Solving Outshining in Star Formation History Inference
- The impact of medium-width bands on the selection, and subsequent luminosity function measurements, of high-z galaxies
- JADES Ultra-red Flattened Objects: Morphologies and Spatial Gradients in Color and Stellar Populations
- Behind the Spotlight: A systematic assessment of outshining using NIRCam medium-bands in the JADES Origins Field
- Unveiling Multiple Physical Processes on a Cluster Galaxy at z=0.3 Using JWST