A large DESI sample reveals thousands of infalling cold gas absorbers at low redshift, with velocity distributions indicating multiple accretion pathways including radial inflows and satellite accretion.
hub Canonical reference
Importance Nested Sampling and the MultiNest Algorithm
Canonical reference. 80% of citing Pith papers cite this work as background.
abstract
Bayesian inference involves two main computational challenges. First, in estimating the parameters of some model for the data, the posterior distribution may well be highly multi-modal: a regime in which the convergence to stationarity of traditional Markov Chain Monte Carlo (MCMC) techniques becomes incredibly slow. Second, in selecting between a set of competing models the necessary estimation of the Bayesian evidence for each is, by definition, a (possibly high-dimensional) integration over the entire parameter space; again this can be a daunting computational task, although new Monte Carlo (MC) integration algorithms offer solutions of ever increasing efficiency. Nested sampling (NS) is one such contemporary MC strategy targeted at calculation of the Bayesian evidence, but which also enables posterior inference as a by-product, thereby allowing simultaneous parameter estimation and model selection. The widely-used MultiNest algorithm presents a particularly efficient implementation of the NS technique for multi-modal posteriors. In this paper we discuss importance nested sampling (INS), an alternative summation of the MultiNest draws, which can calculate the Bayesian evidence at up to an order of magnitude higher accuracy than `vanilla' NS with no change in the way MultiNest explores the parameter space. This is accomplished by treating as a (pseudo-)importance sample the totality of points collected by MultiNest, including those previously discarded under the constrained likelihood sampling of the NS algorithm. We apply this technique to several challenging test problems and compare the accuracy of Bayesian evidences obtained with INS against those from vanilla NS.
hub tools
citation-role summary
citation-polarity summary
representative citing papers
Direct imaging reveals enhanced metallicity and aligned orbit in 29 Cyg b, supporting protoplanetary disk formation near the deuterium burning limit.
A differentiable forward model and likelihood enable probabilistic inference over many spatial morphologies for the Galactic Center gamma-ray Excess using variational methods on GPUs.
SPICE-RACS DR2 delivers the largest single Faraday rotation measure catalog from a radio survey, with 250,000-340,000 RMs across most of the sky at median uncertainty of 2 rad m^{-2}.
Simulation-based inference with a Gaussian process emulator trained on ~1300 POSSIS simulations enables rapid, robust kilonova parameter estimation that avoids MCMC biases from likelihood misspecification.
GJ 1132 b is estimated to have received at least 50 times the cumulative XUV flux of modern Earth with over 95% probability across models, supporting its classification as an atmosphere-free world.
Shape correlations in CEνNS allow likelihood and CNN analyses to discriminate sterile neutrinos from NSI and approximately localize sterile parameters in favorable regions.
dynesty is an open-source Python package for dynamic nested sampling that improves efficiency in Bayesian posterior and evidence estimation compared to MCMC on certain problems.
BAGPIPES fitting of 9289 massive quiescent galaxies shows most SFHs rise gradually then quench in 1-2 Gyr, with faster quenching at z>1 and slower at z<1, interpreted as multiple AGN feedback and gas-supply mechanisms.
The minimal majoron framework permits simultaneous majoron dark matter and thermal leptogenesis in a constrained cosmological window set by freeze-in production, warm dark matter bounds, and indirect detection limits.
Joint NICER+IXPE pulse-profile modeling of SRGA J144459.2-604207 favors large neutron-star mass and radius with two independent hotspots but shows strong sensitivity to joint-analysis methodology.
Massive quiescent galaxies at cosmic noon are compact and bulge-dominated with inside-out quenching, where inner regions formed stars ~0.5 Gyr earlier and quenched faster than outskirts.
Augmented SVGD with birth-death sampling tracks posterior sample changes as prior shape parameters vary in non-conjugate hierarchical models, enabling mode discovery during transitions from unimodal to multimodal posteriors.
Galactic synchrotron emissions above 20 MHz can set tighter upper limits on the abundance of primordial black holes with masses above 10^16 grams than previous cosmic-ray electron data.
Symbolic emulators approximate key Lambda CDM functions to 0.001-0.05% accuracy across relevant redshifts and Omega_m values, enabling faster 3x2pt inference with consistent results.
No evidence for deviations from general relativity is found in LIGO-Virgo binary black hole events, with improved constraints on waveform parameters, graviton mass, and ringdown properties.
For ultra-light PBHs with extended mass distributions, new CMB-derived 95% exclusion bounds on f_PBH are obtained by jointly varying ΛCDM parameters, yielding f_PBH < 1.6 × 10^{-5} for a uniform distribution between 10^{15} and 10^{17} g.
Bilby introduces a user-friendly Python library for accurate Bayesian inference on gravitational-wave signals from compact binaries and other sources, including hierarchical population modeling.
Bayesian hierarchical modeling of photometric redshifts in KiDS+VIKING-450 raises S8 to 0.756 ± 0.039 and reduces Planck tension to 1.9σ.
New measurement of the Hubble parameter H(z=0.12) = 71.33 ± 4.20 km s^{-1} Mpc^{-1} obtained from cosmology-independent stellar ages of passive galaxies in DESI Data Release 1.
No evidence for physics beyond general relativity is found in the analysis of 15 GW events from GWTC-3, with consistency in residuals, PN parameters, and remnant properties.
BILBY is validated on simulated compact binary signals and reproduces the eleven GWTC-1 results with configuration and output files provided for reproduction.
citing papers explorer
-
Peering down the barrel with DESI DR2: 10 000+ inflows at $z$ < 0.6 reveal how galaxies accrete cold gas
A large DESI sample reveals thousands of infalling cold gas absorbers at low redshift, with velocity distributions indicating multiple accretion pathways including radial inflows and satellite accretion.
-
Direct Images of CO2 Absorption in the Atmosphere of a Super-Jupiter: Enhanced Metallicity Suggestive of Formation in a Disk
Direct imaging reveals enhanced metallicity and aligned orbit in 29 Cyg b, supporting protoplanetary disk formation near the deuterium burning limit.
-
High-dimensional inference for the $\gamma$-ray sky with differentiable programming
A differentiable forward model and likelihood enable probabilistic inference over many spatial morphologies for the Galactic Center gamma-ray Excess using variational methods on GPUs.
-
The Rapid ASKAP Continuum Survey VII: Spectra and Polarisation In Cutouts of Extragalactic Sources (SPICE-RACS) Second Data Release -- Unveiling the Magnetised Sky
SPICE-RACS DR2 delivers the largest single Faraday rotation measure catalog from a radio survey, with 250,000-340,000 RMs across most of the sky at median uncertainty of 2 rad m^{-2}.
-
Rapid and robust simulation-based inference for kilonovae
Simulation-based inference with a Gaussian process emulator trained on ~1300 POSSIS simulations enables rapid, robust kilonova parameter estimation that avoids MCMC biases from likelihood misspecification.
-
The Range of Cumulative XUV Flux on GJ 1132 b
GJ 1132 b is estimated to have received at least 50 times the cumulative XUV flux of modern Earth with over 95% probability across models, supporting its classification as an atmosphere-free world.
-
Analytical and Machine Learning Methods for Model Discernment at CE$\nu$NS Experiments
Shape correlations in CEνNS allow likelihood and CNN analyses to discriminate sterile neutrinos from NSI and approximately localize sterile parameters in favorable regions.
-
dynesty: A Dynamic Nested Sampling Package for Estimating Bayesian Posteriors and Evidences
dynesty is an open-source Python package for dynamic nested sampling that improves efficiency in Bayesian posterior and evidence estimation compared to MCMC on certain problems.
-
Inferring the star-formation histories of massive quiescent galaxies with BAGPIPES: Evidence for multiple quenching mechanisms
BAGPIPES fitting of 9289 massive quiescent galaxies shows most SFHs rise gradually then quench in 1-2 Gyr, with faster quenching at z>1 and slower at z<1, interpreted as multiple AGN feedback and gas-supply mechanisms.
-
The Majoron Cosmological Window: Dark Matter and Thermal Leptogenesis
The minimal majoron framework permits simultaneous majoron dark matter and thermal leptogenesis in a constrained cosmological window set by freeze-in production, warm dark matter bounds, and indirect detection limits.
-
Pulse profile modelling of the 2024 outburst of the accreting millisecond pulsar SRGA J144459.2-604207
Joint NICER+IXPE pulse-profile modeling of SRGA J144459.2-604207 favors large neutron-star mass and radius with two independent hotspots but shows strong sensitivity to joint-analysis methodology.
-
Morphological and Star Formation Properties of Cosmic Noon Massive Quiescent Galaxies
Massive quiescent galaxies at cosmic noon are compact and bulge-dominated with inside-out quenching, where inner regions formed stars ~0.5 Gyr earlier and quenched faster than outskirts.
-
Sample continuation in Bayesian hierarchical model via variational inference
Augmented SVGD with birth-death sampling tracks posterior sample changes as prior shape parameters vary in non-conjugate hierarchical models, enabling mode discovery during transitions from unimodal to multimodal posteriors.
-
Constraints on Primordial Black Holes from Galactic Diffuse Synchrotron Emissions
Galactic synchrotron emissions above 20 MHz can set tighter upper limits on the abundance of primordial black holes with masses above 10^16 grams than previous cosmic-ray electron data.
-
Symbolic Emulators for Cosmology: Accelerating Cosmological Analyses Without Sacrificing Precision
Symbolic emulators approximate key Lambda CDM functions to 0.001-0.05% accuracy across relevant redshifts and Omega_m values, enabling faster 3x2pt inference with consistent results.
-
Tests of General Relativity with Binary Black Holes from the second LIGO-Virgo Gravitational-Wave Transient Catalog
No evidence for deviations from general relativity is found in LIGO-Virgo binary black hole events, with improved constraints on waveform parameters, graviton mass, and ringdown properties.
-
CMB constraints on ultra-light primordial black holes with extended mass distributions
For ultra-light PBHs with extended mass distributions, new CMB-derived 95% exclusion bounds on f_PBH are obtained by jointly varying ΛCDM parameters, yielding f_PBH < 1.6 × 10^{-5} for a uniform distribution between 10^{15} and 10^{17} g.
-
Bilby: A user-friendly Bayesian inference library for gravitational-wave astronomy
Bilby introduces a user-friendly Python library for accurate Bayesian inference on gravitational-wave signals from compact binaries and other sources, including hierarchical population modeling.
-
KiDS+VIKING-450 cosmology with Bayesian hierarchical model redshift distributions
Bayesian hierarchical modeling of photometric redshifts in KiDS+VIKING-450 raises S8 to 0.756 ± 0.039 and reduces Planck tension to 1.9σ.
-
New $H(z)$ measurement at Redshift = 0.12 with DESI Data Release 1
New measurement of the Hubble parameter H(z=0.12) = 71.33 ± 4.20 km s^{-1} Mpc^{-1} obtained from cosmology-independent stellar ages of passive galaxies in DESI Data Release 1.
-
Tests of General Relativity with GWTC-3
No evidence for physics beyond general relativity is found in the analysis of 15 GW events from GWTC-3, with consistency in residuals, PN parameters, and remnant properties.
-
Bayesian inference for compact binary coalescences with BILBY: Validation and application to the first LIGO--Virgo gravitational-wave transient catalogue
BILBY is validated on simulated compact binary signals and reproduces the eleven GWTC-1 results with configuration and output files provided for reproduction.