skip to main content


Title: Ain’t No Mountain High Enough: Semiparametric Modeling of LIGO–Virgo’s Binary Black Hole Mass Distribution
Abstract

We introduce a semiparametric model for the primary mass distribution of binary black holes (BBHs) observed with gravitational waves (GWs) that applies a cubic-spline perturbation to a power law. We apply this model to the 46 BBHs included in the second gravitational-wave transient catalog (GWTC-2). The spline perturbation model recovers a consistent primary mass distribution with previous results, corroborating the existence of a peak at 35M(>97% credibility) found with the Powerlaw+Peakmodel. The peak could be the result of pulsational pair-instability supernovae. The spline perturbation model finds potential signs of additional features in the primary mass distribution at lower masses similar to those previously reported by Tiwari and Fairhurst. However, with fluctuations due to small-number statistics, the simpler Powerlaw+Peakand BrokenPowerlawmodels are both still perfectly consistent with observations. Our semiparametric approach serves as a way to bridge the gap between parametric and nonparametric models to more accurately measure the BBH mass distribution. With larger catalogs we will be able to use this model to resolve possible additional features that could be used to perform cosmological measurements and will build on our understanding of BBH formation, stellar evolution, and nuclear astrophysics.

 
more » « less
Award ID(s):
2146528
NSF-PAR ID:
10361793
Author(s) / Creator(s):
; ; ;
Publisher / Repository:
DOI PREFIX: 10.3847
Date Published:
Journal Name:
The Astrophysical Journal
Volume:
924
Issue:
2
ISSN:
0004-637X
Format(s):
Medium: X Size: Article No. 101
Size(s):
["Article No. 101"]
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Several features in the mass spectrum of merging binary black holes (BBHs) have been identified using data from the Third Gravitational Wave Transient Catalog (GWTC-3). These features are of particular interest as they may encode the uncertain mechanism of BBH formation. We assess if the features are statistically significant or the result of Poisson noise due to the finite number of observed events. We simulate catalogs of BBHs whose underlying distribution does not have the features of interest, apply the analysis previously performed on GWTC-3, and determine how often such features are spuriously found. We find that one of the features found in GWTC-3, the peak at ∼35M, cannot be explained by Poisson noise alone: peaks as significant occur in 1.7% of catalogs generated from a featureless population. This peak is therefore likely to be of astrophysical origin. The data is suggestive of an additional significant peak at ∼10M, though the exact location of this feature is not resolvable with current observations. Additional structure beyond a power law, such as the purported dip at ∼14M, can be explained by Poisson noise. We also provide a publicly available package,GWMockCat, that creates simulated catalogs of BBH events with correlated measurement uncertainty and selection effects according to user-specified underlying distributions and detector sensitivities.

     
    more » « less
  2. Abstract

    We introduce the first complete nonparametric model for the astrophysical distribution of the binary black hole (BBH) population. Constructed from basis splines, we use these models to conduct the most comprehensive data-driven investigation of the BBH population to date, simultaneously fitting nonparametric models for the BBH mass ratio, spin magnitude and misalignment, and redshift distributions. With GWTC-3, we report the same features previously recovered with similarly flexible models of the mass distribution, most notably the peaks in merger rates at primary masses of ∼10Mand ∼35M. Our model reports a suppressed merger rate at low primary masses and a mass-ratio distribution consistent with a power law. We infer a distribution for primary spin misalignments that peaks away from alignment, supporting conclusions of recent work. We find broad agreement with the previous inferences of the spin magnitude distribution: the majority of BBH spins are small (a< 0.5), the distribution peaks ata∼ 0.2, and there is mild support for a nonspinning subpopulation, which may be resolved with larger catalogs. With a modulated power law describing the BBH merger rate’s evolution in redshift, we see hints of the rate evolution either flattening or decreasing atz∼ 0.2–0.5, but the full distribution remains entirely consistent with a monotonically increasing power law. We conclude with a discussion of the astrophysical context of our new findings and how nonparametric methods in gravitational-wave population inference are uniquely poised to complement to the parametric approach as we enter the data-rich era of gravitational-wave astronomy.

     
    more » « less
  3. ABSTRACT

    The black hole (BH) masses measured from gravitational wave observations appear to cluster around specific mass values. Consequently, the primary (and chirp) mass distribution of binary black holes (BBHs) inferred using these measurements shows four emerging peaks. These peaks are approximately located at a primary (chirp) mass value of 10 $\, \mathrm{M}_\odot$ (8$\, \mathrm{M}_\odot$), 20 $\, \mathrm{M}_\odot$ (14 $\, \mathrm{M}_\odot$), 35 $\, \mathrm{M}_\odot$ (28 $\, \mathrm{M}_\odot$), and 63 $\, \mathrm{M}_\odot$ (49 $\, \mathrm{M}_\odot$). Although the presence of the first and third peaks has been attributed to BBH formation in star clusters or due to the evolution of stellar binaries in isolation, the second peak has received relatively less attention because it lacks significance in the primary mass distribution. In this article, we report that confidence in the second peak depends on the mass parameter we choose to model the population on. Unlike primary mass, this peak is significant when modelled on the chirp mass. We discuss the disparity as a consequence of mass asymmetry in the observations that cluster at the second peak. Finally, we report this asymmetry as part of a potential trend in the mass ratio distribution manifested as a function of the chirp mass, but not as a function of primary mass, when we include the observation GW190814 in our modelling. The chirp mass is not a parameter of astrophysical relevance. Features present in the chirp mass, but not in the primary mass, are relatively difficult to explain and expected to garner significant interest.

     
    more » « less
  4. Abstract

    The population-level distributions of the masses, spins, and redshifts of binary black holes (BBHs) observed using gravitational waves can shed light on how these systems form and evolve. Because of the complex astrophysical processes shaping the inferred BBH population, models allowing for correlations among these parameters will be necessary to fully characterize these sources. We hierarchically analyze the BBH population detected by LIGO and Virgo with a model allowing for correlations between the effective aligned spin and the primary mass and redshift. We find that the width of the effective spin distribution grows with redshift at 98.6% credibility. We determine this trend to be robust under the application of several alternative models and additionally verify that such a correlation is unlikely to be spuriously introduced using a simulated population. We discuss the possibility that this correlation could be due to a change in the natal black hole spin distribution with redshift.

     
    more » « less
  5. Abstract

    In their most recent observing run, the LIGO-Virgo-KAGRA Collaboration observed gravitational waves from compact binary mergers with highly asymmetric mass ratios, including both binary black holes (BBHs) and neutron star-black holes (NSBHs). It appears that NSBHs with mass ratiosq≃ 0.2 are more common than equally asymmetric BBHs, but the reason for this remains unclear. We use the binary population synthesis codecosmicto investigate the evolutionary pathways leading to the formation and merger of asymmetric compact binaries. We find that within the context of isolated binary stellar evolution, most asymmetric mergers start off as asymmetric stellar binaries. Because of the initial asymmetry, these systems tend to first undergo a dynamically unstable mass transfer phase. However, after the first star collapses into a compact object, the mass ratio is close to unity and the second phase of mass transfer is usually stable. According to our simulations, this stable mass transfer fails to shrink the orbit enough on its own for the system to merge. Instead, the natal kick received by the second-born compact object during its collapse is key in determining how many of these systems can merge. For the most asymmetric systems with mass ratios ofq≤ 0.1, the merging systems in our models receive an average kick magnitude of 255 km s−1during the second collapse, while the average kick for non-merging systems is 59 km s−1. Because lower mass compact objects, like neutron stars, are expected to receive larger natal kicks than higher mass BHs, this may explain why asymmetric NSBH systems merge more frequently than asymmetric BBH systems.

     
    more » « less