SIDE-real: Supernova Ia Dust Extinction with truncated marginal neural ratio estimation applied to real data

Abstract

We present the first fully simulation-based hierarchical analysis of the light curves of a population of low-redshift type Ia supernovae (SNae Ia). Our hardware-accelerated forward model, released in the Python package slicsim, includes stochastic variations of each SN's spectral flux distribution (based on the pre-trained BayeSN model), extinction from dust in the host and in the Milky Way, redshift, and realistic instrumental noise. By utilising truncated marginal neural ratio estimation (TMNRE), a neural network-enabled simulation-based inference technique, we implicitly marginalise over 4000 latent variables (for a set of $\approx 100$ SNae Ia) to efficiently infer SN Ia absolute magnitudes and host-galaxy dust properties at the population level while also constraining the parameters of individual objects. Amortisation of the inference procedure allows us to obtain coverage guarantees for our results through Bayesian validation and frequentist calibration. Furthermore, we show a detailed comparison to full likelihood-based inference, implemented through Hamiltonian Monte Carlo, on simulated data and then apply TMNRE to the light curves of 86 SNae Ia from the Carnegie Supernova Project, deriving marginal posteriors in excellent agreement with previous work. Given its ability to accommodate arbitrarily complex extensions to the forward model -- e.g. different populations based on host properties, redshift evolution, complicated photometric redshift estimates, selection effects, and non-Ia contamination -- without significant modifications to the inference procedure, TMNRE has the potential to become the tool of choice for cosmological parameter inference from future, large SN Ia samples.

16 pages, 10 figures; Published in MNRAS; (c): 2024 The Author(s). Published by Oxford University Press on behalf of Royal Astronomical Society

References in corpus (44)

• Improving neural networks by preventing co-adaptation of feature detectors
• A Comprehensive Measurement of the Local Value of the Hubble Constant with 1 km/s/Mpc Uncertainty from the Hubble Space Telescope and the SH0ES Team
• In the Realm of the Hubble tension $-$ a Review of Solutions
• The Pantheon+ Analysis: Cosmological Constraints
• PolyChord: nested sampling for cosmology
• Cosmological parameters from the comparison of peculiar velocities with predictions from the 2M++ density field
• Correcting Type Ia Supernova Distances for Selection Biases and Contamination in Photometrically Identified Samples
• Composable Effects for Flexible and Accelerated Probabilistic Programming in NumPyro
• Fast likelihood-free cosmology with neural density estimators and active learning
• The Carnegie Supernova Project I: Third Photometry Data Release of Low-Redshift Type Ia Supernovae and Other White Dwarf Explosions
• The Pantheon+ Analysis: SuperCal-Fragilistic Cross Calibration, Retrained SALT2 Light Curve Model, and Calibration Systematic Uncertainty
• SALT3: An Improved Type Ia Supernova Model for Measuring Cosmic Distances
• Type Ia Supernova Light Curve Inference: Hierarchical Bayesian Analysis in the Near Infrared
• Nuisance hardened data compression for fast likelihood-free inference
• Bayesian Estimation Applied to Multiple Species: Towards cosmology with a million supernovae
• Improved Treatment of Host-Galaxy Correlations in Cosmological Analyses With Type Ia Supernovae
• Think Global, Act Local: The Influence of Environment Age and Host Mass on Type Ia Supernova Light Curves
• Union Through UNITY: Cosmology with 2,000 SNe Using a Unified Bayesian Framework
• Testing the Consistency of Dust Laws in SN Ia Host Galaxies: A BayeSN Examination of Foundation DR1
• A Revised SALT2 Surface for Fitting Type Ia Supernova Light Curves
• SNIa-Cosmology Analysis Results from Simulated LSST Images: from Difference Imaging to Constraints on Dark Energy
• SICRET: Supernova Ia Cosmology with truncated marginal neural Ratio EsTimation
• Estimating the warm dark matter mass from strong lensing images with truncated marginal neural ratio estimation
• Constraining the X-ray heating and reionization using 21-cm power spectra with Marginal Neural Ratio Estimation
• Dark energy by natural evolution: Constraining dark energy using Approximate Bayesian Computation
• Constraining the SN Ia Host Galaxy Dust Law Distribution and Mass Step: Hierarchical BayeSN Analysis of Optical and Near-Infrared Light Curves
• Likelihood-free Inference with Mixture Density Network
• One never walks alone: the effect of the perturber population on subhalo measurements in strong gravitational lenses
• The Dark Energy Survey Supernova Program: Cosmological Analysis and Systematic Uncertainties
• Using Rest-Frame Optical and NIR Data from the RAISIN Survey to Explore the Redshift Evolution of Dust Laws in SN Ia Host Galaxies
• A Probabilistic Autoencoder for Type Ia Supernovae Spectral Time Series
• Towards Reliable Simulation-Based Inference with Balanced Neural Ratio Estimation
• Test of Artificial Neural Networks in Likelihood-free Cosmological Constraints: A Comparison of IMNN and DAE
• What to do when things get crowded? Scalable joint analysis of overlapping gravitational wave signals
• Amortized Bayesian Inference for Supernovae in the Era of the Vera Rubin Observatory Using Normalizing Flows
• The Pantheon+ Analysis: Forward-Modeling the Dust and Intrinsic Colour Distributions of Type Ia Supernovae, and Quantifying their Impact on Cosmological Inferences
• Detection is truncation: studying source populations with truncated marginal neural ratio estimation
• Stratified Learning: A General-Purpose Statistical Method for Improved Learning under Covariate Shift
• SimSIMS: Simulation-based Supernova Ia Model Selection with thousands of latent variables
• Bayesian Simulation-based Inference for Cosmological Initial Conditions
• Fishnets: Information-Optimal, Scalable Aggregation for Sets and Graphs
• The Impact of Dust on Cepheid and Type Ia Supernova Distances
• Scalable hierarchical BayeSN inference: Investigating dependence of SN Ia host galaxy dust properties on stellar mass and redshift
• Analytic auto-differentiable $\Lambda$CDM cosmography

Cited by in corpus (7)

• Simulation-Based Inference of the sky-averaged 21-cm signal from CD-EoR with REACH
• Fast likelihood-free inference in the LSS Stage IV era
• STAR NRE: Solving supernova selection effects with set-based truncated auto-regressive neural ratio estimation
• CIGaRS I: Combined simulation-based inference from type Ia supernovae and host photometry
• FlowSN: Neural Simulation-Based Inference under Realistic Selection Effects applied to Supernova Cosmology
• Calibrating Bayesian Tension Statistics using Neural Ratio Estimation
• Diagnosing Systematic Effects Using the Inferred Initial Power Spectrum