skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: What You Don’t Know Can Hurt You: Use and Abuse of Astrophysical Models in Gravitational-wave Population Analyses
Abstract One of the goals of gravitational-wave astrophysics is to infer the number and properties of the formation channels of binary black holes (BBHs); to do so, one must be able to connect various models with the data. We explore benefits and potential issues with analyses using models informed by population synthesis. We consider five possible formation channels of BBHs, as in Zevin et al. (2021b). First, we confirm with the GWTC-3 catalog what Zevin et al. (2021b) found in the GWTC-2 catalog, i.e., that the data are not consistent with the totality of observed BBHs forming in any single channel. Next, using simulated detections, we show that the uncertainties in the estimation of the branching ratios can shrink by up to a factor of ∼1.7 as the catalog size increases from 50 to 250, within the expected number of BBH detections in LIGO–Virgo–KAGRA's fourth observing run. Finally, we show that this type of analysis is prone to significant biases. By simulating universes where all sources originate from a single channel, we show that the influence of the Bayesian prior can make it challenging to conclude that one channel produces all signals. Furthermore, by simulating universes where all five channels contribute but only a subset of channels are used in the analysis, we show that biases in the branching ratios can be as large as ∼50% with 250 detections. This suggests that caution should be used when interpreting the results of analyses based on strongly modeled astrophysical subpopulations.  more » « less
Award ID(s):
2045740
PAR ID:
10465322
Author(s) / Creator(s):
; ;
Publisher / Repository:
DOI PREFIX: 10.3847
Date Published:
Journal Name:
The Astrophysical Journal
Volume:
955
Issue:
2
ISSN:
0004-637X
Format(s):
Medium: X Size: Article No. 127
Size(s):
Article No. 127
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract The relative spin orientations of black holes (BHs) in binaries encode their evolutionary history: BHs assembled dynamically should have isotropically distributed spins, while spins of BHs originating in the field should be aligned with the orbital angular momentum. In this article, we introduce a simple population model for these dynamical and field binaries that uses spin orientations as an anchor to disentangle these two evolutionary channels. We then analyze binary BH mergers in the Third Gravitational-Wave Transient Catalog (GWTC-3) and ask whether BHs from the isotropic-spin population possess different distributions of mass ratios, spin magnitudes, or redshifts from the preferentially aligned-spin population. We find no compelling evidence that binary BHs in GWTC-3 have different source-property distributions depending on their spin alignment, but we do find that the dynamical and field channels cannot both have mass-ratio distributions that strongly favor equal masses. We give an example of how this can be used to provide insights into the various processes that drive these BHs to merge. We also find that the current detections are insufficient in extracting differences in spin magnitude or redshift distributions of isotropic and aligned-spin populations. 
    more » « less
  2. Abstract We present a search for outer solar system objects in the 6 yr of data from the Dark Energy Survey (DES). The DES covered a contiguous 5000 deg 2 of the southern sky with ≈80,000 3 deg 2 exposures in the grizY filters between 2013 and 2019. This search yielded 812 trans-Neptunian objects (TNOs), one Centaur and one Oort cloud comet, 458 reported here for the first time. We present methodology that builds upon our previous search on the first 4 yr of data. All images were reprocessed with an optimized detection pipeline that leads to an average completeness gain of 0.47 mag per exposure, as well as improved transient catalog production and algorithms for linkage of detections into orbits. All objects were verified by visual inspection and by the “sub-threshold significance,” the signal-to-noise ratio in the stack of images in which its presence is indicated by the orbit, but no detection was reported. This yields a pure catalog complete to r ≈ 23.8 mag and distances 29 < d < 2500 au. The TNOs have minimum (median) of 7 (12) nights’ detections and arcs of 1.1 (4.2) yr, and will have grizY magnitudes available in a further publication. We present software for simulating our observational biases for comparisons of models to our detections. Initial inferences demonstrating the catalog’s statistical power are: the data are inconsistent with the CFEPS-L7 model for the classical Kuiper Belt; the 16 “extreme” TNOs ( a > 150 au, q > 30 au) are consistent with the null hypothesis of azimuthal isotropy; and nonresonant TNOs with q > 38 au, a > 50 au show a significant tendency to be sunward of major mean-motion resonances. 
    more » « less
  3. 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
  4. This article presents tactile drafting techniques developed in collaboration with blind educators and students that have the potential to increase BLV students’ access to drafting and engineering graphic curriculum in K-12 and higher education. This work builds on previous work funded by the National Science Foundation (Goodridge et al., 2019; Ashby et al., 2018; Lopez et al., 2020; Goodridge et al., 2021a; Goodridge et al., 2021b) and it is the authors’ hope that some of the practices included herein will allow BLV youth to further develop technological and engineering literacy in related technology and engineering graphics courses. 
    more » « less
  5. During cosmic timescales, supermassive binary black holes (SMBBHs) form by galaxy merg- ers, where each galaxy hosts a supermassive black hole (SMBH) at its center. By studying SMBBHs, we can gain insights into galaxy evolution and black hole growth. However, the typical separation between black holes in SMBBHs is usually below 1 pc, making them dif- ficult to resolve using direct imaging or photometry. To be able to distinguish binary black holes (BBHs) from typical AGN powered by single black holes (SBHs), we conducted this research project to develop a new diagnostic method to identify a unique feature of SMBBHs, which can be used to distinguish AGN powered by BBHs from those powered by SBHs. The basic idea of this method is that BBHs have different configurations compared with those of SBHs, such as the circumbinary disk enveloping the whole binary system, the mini- disks around each black hole, and the streams between the circumbinary disk and minidisks. It is these different configurations of accretion disks of black holes that lead to the difference in the spectral energy distributions (SEDs), which in turn will differently photoionize the broad line region (BLR) and produce different strengths of broad emission lines. Thus, it is these differences in strengths of broad emission lines that we expect to see between two different configurations of black holes, either binary or single black holes. By identifying these differences in line strengths, we can distinguish between binary and single black holes. In order to achieve this goal of distinguishing BBHs from SBHs, we organized our project in the steps below: (1) obtained BBH SEDs from Gutiérrez et al., 2022, (2) generated SBH SEDs using XSPEC modeling code OPTXAGNF, (3) produced single-cloud models, which represent the BLR, input BBH and SBH SEDs into the models, and simulated the single0cloud response with a photoionization code Cloudy, (4) built cloud-ensemble models, which represent a more realistic BLR, input BBH and SBH SEDs into the models, and simulated the emission-line response within those clouds using a broad emission line mapping code BELMAC. In step (3), we have explored the differences in line ratios between BBH and SBH for a representative single-cloud photoionization model. The emission lines we used here are: Si IV λ1400Å, C III] λ1909Å, C IV λ1549Å, Mg II λ2798Å, and Lyα λ1216Å. It turned i out that differences do exist, but they are too small to be identified in observational data. Furthermore, we have investigated the line equivalent widths predicted by SBH and BBH models respectively. By doing so, we found some apparent differences between BBHs and SBHs in some specific emission lines: Lyα, CIV, and Hα. However, these differences vanish at the highest mass of black holes (109M⊙). In step (4), we continued the investigation of the equivalent width between SBH and BBH BLR cloud-ensemble models and found that some emission lines show the difference between BBHs and SBHs, such as CIV in the case of BH mass 107M⊙, U ∝ r−2, and log n/1 cm−3 = 10.5, 11.0. For the highest mass in the case of U ∝ r−2, the results are similar to the one in single-cloud models: no difference is shown between BBHs and SBHs across all emission lines. Most importantly, in this step, we found that for CIV λ1549 Å in the case of a black hole mass of 107M⊙, s = −2, and log n/1 cm−3 = 10.5, only the SBH EW falls inside the range of the observed range in SDSS DR7 Quasar Catalog while the BBH EW falls outside the range and becomes an outlier. This is what we want to find to distinguish BBHs from SBHs in the observational data. 
    more » « less