Flexible Proposal Distributions for Gravitational-Wave Population Inference with Gaussian Mixture Models
Year: 2026
Type: preprint
Abstract: Population analyses of gravitational-wave (GW) sources use hierarchical Bayesian inference to constrain astrophysical models while accounting for selection effects. The likelihood is usually estimated by Monte Carlo integration over finite sets of event posterior samples and simulated injections. As catalogs grow, more samples are needed to control the Monte Carlo variance, but the available event samples are fixed by parameter estimation. Sharp population features can also lie in regions poorly sampled by the event posteriors. We fit Gaussian mixture models (GMMs) to the event posteriors and evaluate the likelihood by importance sampling with two proposals: an ensemble of event posteriors and a mixture of that ensemble with a fiducial population. Together with a GMM representation of detected injections, this allows additional samples to be drawn for both event and selection integrals. Increasing the sample count reduces Monte Carlo variance but leaves uncertainty from fitting the densities to finite training sets. In simulated catalogs, the GMM-based estimators show no detectable additional bias relative to the posterior-sample estimator; their calibration is comparable within the limitations of the tests. An application to GWTC-3 with a different mass model gives consistent results across the estimators. The method separates the numerical sampling budget from the original sample count, while requiring explicit checks of density-fit uncertainty and approximation error.
Source: arXiv (Cornell University)
Institutions: University of Wisconsin–Milwaukee, Carnegie Mellon University
Language: English
Subfield: Astronomy and Astrophysics
Field: Physics and Astronomy
Domain: Physical Sciences
Open Access status: green
OpenAlex ID: W7223345644