STROOPWAFEL: Simulating rare outcomes from astrophysical populations, with application to gravitational-wave sources
arXiv:1905.00910 · doi:10.1093/mnras/stz2558
Abstract
Gravitational-wave observations of double compact object (DCO) mergers are providing new insights into the physics of massive stars and the evolution of binary systems. Making the most of expected near-future observations for understanding stellar physics will rely on comparisons with binary population synthesis models. However, the vast majority of simulated binaries never produce DCOs, which makes calculating such populations computationally inefficient. We present an importance sampling algorithm, STROOPWAFEL, that improves the computational efficiency of population studies of rare events, by focusing the simulation around regions of the initial parameter space found to produce outputs of interest. We implement the algorithm in the binary population synthesis code COMPAS, and compare the efficiency of our implementation to the standard method of Monte Carlo sampling from the birth probability distributions. STROOPWAFEL finds 25-200 times more DCO mergers than the standard sampling method with the same simulation size, and so speeds up simulations by up to two orders of magnitude. Finding more DCO mergers automatically maps the parameter space with far higher resolution than when using the traditional sampling. This increase in efficiency also leads to a decrease of a factor 3-10 in statistical sampling uncertainty for the predictions from the simulations. This is particularly notable for the distribution functions of observable quantities such as the black hole and neutron star chirp mass distribution, including in the tails of the distribution functions where predictions using standard sampling can be dominated by sampling noise.
arXiv commentAccepted. Data and scripts to reproduce main results is publicly available. The code for the STROOPWAFEL algorithm will be made publicly available. Early inquiries can be addressed to the lead author
References in corpus (40)
- Observation of Gravitational Waves from a Binary Black Hole Merger
- The NumPy array: a structure for efficient numerical computation
- GW170817: Observation of Gravitational Waves from a Binary Neutron Star Inspiral
- On the variation of the Initial Mass Function
- GWTC-1: A Gravitational-Wave Transient Catalog of Compact Binary Mergers Observed by LIGO and Virgo during the First and Second Observing Runs
- Binary interaction dominates the evolution of massive stars
- Evolution of binary stars and the effect of tides on binary populations
- Comprehensive analytic formulae for stellar evolution as a function of mass and metallicity
- Mind your Ps and Qs: the Interrelation between Period (P) and Mass-ratio (Q) Distributions of Binary Stars
- Compact Remnant Mass Function: Dependence on the Explosion Mechanism and Metallicity
- A Comprehensive Study of Binary Compact Objects as Gravitational Wave Sources: Evolutionary Channels, Rates, and Physical Properties
- The first gravitational-wave source from the isolated evolution of two 40-100 Msun stars
- Pulsational Pair-Instability Supernovae
- Double Compact Objects I: The Significance Of The Common Envelope On Merger Rates
- A new route towards merging massive black holes
- Merging binary black holes formed through chemically homogeneous evolution in short-period stellar binaries
- The chemically homogeneous evolutionary channel for binary black hole mergers: rates and properties of gravitational-wave events detectable by advanced LIGO
- Progenitors of gravitational wave mergers: Binary evolution with the stellar grid-based code ComBinE
- Formation of the first three gravitational-wave observations through isolated binary evolution
- On the formation history of Galactic double neutron stars
- The effect of the metallicity-specific star formation history on double compact object mergers
- Spin-Orbit Misalignment in Close Binaries with Two Compact Objects
- Remnant baryon mass outside of the black hole after a neutron star-black hole merger
- A Highly Spinning and Aligned Binary Black Hole Merger in the Advanced LIGO First Observing Run
- Merging stellar-mass binary black holes
- Orbital Characteristics of Binary Systems after Asymmetric Supernova Explosions
- Adaptive Importance Sampling in General Mixture Classes
- Where are LIGO's Big Black Holes?
- The impact of pair-instability mass loss on the binary black hole mass distribution
- Formation of Low-Mass X-Ray Binaries. II. Common Envelope Evolution of Primordial Binaries with Extreme Mass Ratios
- Accuracy of inference on the physics of binary evolution from gravitational-wave observations
- Mining Gravitational-wave Catalogs To Understand Binary Stellar Evolution: A New Hierarchical Bayesian Framework
- Impact of inter-correlated initial binary parameters on double black hole and neutron star mergers
- r-process enrichment of ultra-faint dwarf galaxies by fast merging double neutron stars
- Dirichlet Process Gaussian-mixture model: An application to localizing coalescing binary neutron stars with gravitational-wave observations
- Adaptive Importance Sampling for Estimation in Structured Domains
- dart_board: Binary Population Synthesis with Markov Chain Monte Carlo
- Exploring the Parameter Space of Compact Binary Population Synthesis
- Optimal mixture weights in multiple importance sampling
- Astrophysics of stellar black holes
Cited by in corpus (35)
- COSMIC Variance in Binary Population Synthesis
- The effect of the metallicity-specific star formation history on double compact object mergers
- Hierarchical mergers of stellar-mass black holes and their gravitational-wave signatures
- The origin of spin in binary black holes: Predicting the distributions of the main observables of Advanced LIGO
- The impact of pair-instability mass loss on the binary black hole mass distribution
- Impact of Massive Binary Star and Cosmic Evolution on Gravitational Wave Observations I: Black Hole-Neutron Star Mergers
- The redshift evolution of the binary black hole merger rate: a weighty matter
- Polluting the pair-instability mass gap for binary black holes through super-Eddington accretion in isolated binaries
- Kilonova Luminosity Function Constraints based on Zwicky Transient Facility Searches for 13 Neutron Star Mergers
- Probing multiple populations of compact binaries with third-generation gravitational-wave detectors
- No peaks without valleys: The stable mass transfer channel for gravitational-wave sources in light of the neutron star-black hole mass gap
- Multiple Stellar Evolution: a population synthesis algorithm to model the stellar, binary, and dynamical evolution of multiple-star systems
- Modelling Double Neutron Stars: Radio and Gravitational Waves
- Gravitational wave sources in our Galactic backyard: Predictions for BHBH, BHNS and NSNS binaries detectable with LISA
- Signatures of mass ratio reversal in gravitational waves from merging binary black holes
- On the single-event-based identification of primordial black hole mergers at cosmological distances
- Constraints on the contributions to the observed binary black hole population from individual evolutionary pathways in isolated binary evolution
- Measuring properties of primordial black hole mergers at cosmological distances: effect of higher order modes in gravitational waves
- Exploring the impact of IMF and binary parameter stochasticity with a binary population synthesis code
- Population synthesis of accreting white dwarfs: Rates and evolutionary pathways of H and He novae
- Formation of GW230529 from Isolated Binary Evolution
- Active Learning for Computationally Efficient Distribution of Binary Evolution Simulations
- Surrogate forward models for population inference on compact binary mergers
- Recovering Injected Astrophysics from the LISA Galactic Double White Dwarf Binaries
- Formation of the First Two Black Hole-Neutron Star Mergers (GW200115 and GW200105) from Isolated Binary Evolution
- The role of triple evolution in the formation of LISA double white dwarfs
- Where to find over-massive brown dwarfs: new benchmark systems for binary evolution
- Cosmology with Binary Neutron Stars: Does Mass-Redshift Correlation Matter?
- Luminous Fast Blue Optical Transients as "Failed" Gravitational-wave Sources: Helium CoreBlack Hole Mergers Following Delayed Dynamical Instability
- Single-event likelihood of star cluster properties with LIGO-Virgo-Kagra binary black hole observations
- Aluminium-26 production in low- and intermediate-mass binary systems
- Using Binary Population Synthesis to Calculate the Yields of Low- and Intermediate-Mass Binary Populations at Low Metallicity
- Evolutionary tracks of binary neutron star progenitors across cosmic times
- Detailed properties of gravitational-wave mergers from flyby perturbations of wide binary black holes in the field
- Impact of Massive Binary Star and Cosmic Evolution on Gravitational Wave Observations II: Double Compact Object Rates and Properties