Abstract We propose a Bayesian inference framework to predict the merger history of LIGO-Virgo binary black holes (BHs), whose binary components may have undergone hierarchical mergers in the past. The framework relies on numerical relativity predictions for the mass, spin, and kick velocity of the remnant BHs. This proposed framework computes the masses, spins, and kicks imparted to the remnant of the parent binaries, given the initial masses and spin magnitudes of the binary constituents. We validate our approach by performing an “injection study” based on a constructed sequence of hierarchically formed binaries. Noise is added to the final binary in the sequence, and the parameters of the “parent” and “grandparent” binaries in the merger chain are then reconstructed. This method is then applied to three GWTC-3 events: GW190521, GW200220_061928, and GW190426_190642. These events were selected because at least one of the binary companions lies in the putative pair-instability supernova mass gap, in which stellar processes alone cannot produce BHs. Hierarchical mergers offer a natural explanation for the formation of BHs in the pair-instability mass gap. We use the backward evolution framework to predict the parameters of the parents of the primary companion of these three binaries. For instance, the parent binary of GW190521 has masses and within the 90% credible interval. Astrophysical environments with escape speeds ≥100 km s−1are preferred sites to host these events. Our approach can be readily applied to future high-mass gravitational wave events to predict their formation history under the hierarchical merger assumption.
more »
« less
Two of a Kind: Comparing Big and Small Black Holes in Binaries with Gravitational Waves
Abstract When modeling the population of merging binary black holes, analyses have generally focused on characterizing the distribution of primary (i.e., more massive) black holes in the binary, while using simplistic prescriptions for the distribution of secondary masses. However, the secondary mass distribution and its relationship to the primary mass distribution provide a fundamental observational constraint on the formation history of coalescing binary black holes. If both black holes experience similar stellar evolutionary processes prior to collapse, as might be expected in dynamical formation channels, the primary and secondary mass distributions would show similar features. If they follow distinct evolutionary pathways (for example, due to binary interactions that break symmetry between the initially more massive and less massive stars), their mass distributions may differ. We present the first analysis of the binary black hole population that explicitly fits for the secondary mass distribution. We find that the data is consistent with a ∼30M⊙peak existing only in the distribution of secondary rather than primary masses. This would have major implications for our understanding of the formation of these binaries. Alternatively, the data is consistent with the peak existing in both component mass distributions, a possibility not included in most previous studies. In either case, the peak is observed at , which is shifted lower than the value obtained in previous analyses of the marginal primary mass distribution, placing this feature in further tension with expectations from a pulsational pair-instability supernova pileup.
more »
« less
- Award ID(s):
- 2110507
- PAR ID:
- 10489869
- Publisher / Repository:
- DOI PREFIX: 10.3847
- Date Published:
- Journal Name:
- The Astrophysical Journal
- Volume:
- 962
- Issue:
- 1
- ISSN:
- 0004-637X
- Format(s):
- Medium: X Size: Article No. 69
- Size(s):
- Article No. 69
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract We search for features in the mass distribution of detected compact binary coalescences which signify the transition between neutron stars (NSs) and black holes (BHs). We analyze all gravitational-wave (GW) detections by the LIGO Scientific Collaboration, the Virgo Collaboration, and the KAGRA Collaboration (LVK) made through the end of the first half of the third observing run, and find clear evidence for two different populations of compact objects based solely on GW data. We confidently (99.3%) find a steepening relative to a single power law describing NSs and low-mass BHs below , which is consistent with many predictions for the maximum NS mass. We find suggestions of the purported lower mass gap between the most massive NSs and the least massive BHs, but are unable to conclusively resolve it with current data. If it exists, we find the lower mass gap’s edges to lie at and . We reexamine events that have been deemed “exceptional” by the LVK collaborations in the context of these features. We analyze GW190814 self-consistently in the context of the full population of compact binaries, finding support for its secondary to be either a NS or a lower mass gap object, consistent with previous claims. Our models are the first to accommodate this event, which is an outlier with respect to the binary BH population. We find that GW200105 and GW200115 probe the edges of, and may have components within, the lower mass gap. As future data improve global population models, the classification of these events will also improve.more » « less
-
Abstract Gravitational-wave (GW) detections of merging neutron star–black hole (NSBH) systems probe astrophysical neutron star (NS) and black hole (BH) mass distributions, especially at the transition between NS and BH masses. Of particular interest are the maximum NS mass, minimum BH mass, and potential mass gap between them. While previous GW population analyses assumed all NSs obey the same maximum mass, if rapidly spinning NSs exist, they can extend to larger maximum masses than nonspinning NSs. In fact, several authors have proposed that the ∼2.6M⊙object in the event GW190814—either the most massive NS or least massive BH observed to date—is a rapidly spinning NS. We therefore infer the NSBH mass distribution jointly with the NS spin distribution, modeling the NS maximum mass as a function of spin. Using four LIGO–Virgo NSBH events including GW190814, if we assume that the NS spin distribution is uniformly distributed up to the maximum (breakup) spin, we infer the maximum nonspinning NS mass is (90% credibility), while assuming only nonspinning NSs, the NS maximum mass must be >2.53M⊙(90% credibility). The data support the mass gap’s existence, with a minimum BH mass at . With future observations, under simplified assumptions, 150 NSBH events may constrain the maximum nonspinning NS mass to ±0.02M⊙, and we may even measure the relation between the NS spin and maximum mass entirely from GW data. If rapidly rotating NSs exist, their spins and masses must be modeled simultaneously to avoid biasing the NS maximum mass.more » « less
-
Abstract In the pursuit of understanding the population of stellar remnants within the Milky Way, we analyze the sample of ∼950 microlensing events observed by the Spitzer Space Telescope between 2014 and 2019. In this study we focus on a subsample of nine microlensing events, selected based on their long timescales, small microlensing parallaxes, and joint observations by the Gaia mission, to increase the probability that the chosen lenses are massive and the mass is measurable. Among the selected events we identify lensing black holes and neutron star candidates, with potential confirmation through forthcoming release of the Gaia time-series astrometry in 2026. Utilizing Bayesian analysis and Galactic models, along with the Gaia Data Release 3 proper-motion data, four good candidates for dark remnants were identified: OGLE-2016-BLG-0293, OGLE-2018-BLG-0483, OGLE-2018-BLG-0662, and OGLE-2015-BLG-0149, with lens masses of , , and , respectively. Notably, the first two candidates are expected to exhibit astrometric microlensing signals detectable by Gaia, offering the prospect of validating the lens masses. The methodologies developed in this work will be applied to the full Spitzer microlensing sample, populating and analyzing the timescale (tE) versus parallax (πE) diagram to derive constraints on the population of lenses in general and massive remnants in particular.more » « less
-
Abstract The scale ofα-element yields is difficult to predict from theory because of uncertainties in massive star evolution, supernova physics, and black hole formation, and it is difficult to constrain empirically because the impact of higher yields can be compensated by greater metal loss in galactic winds. We use a recent measurement of the mean iron yield of core collapse supernovae (CCSN) by Rodriguez et al., , to infer the scale ofα-element yields by assuming that the plateau of [α/Fe] abundance ratios observed in low-metallicity stars represents the yield ratio of CCSN. For a plateau at [α/Fe]cc= 0.45, we find that the population-averaged yields of O and Mg are about equal to the solar abundance of these elements, , where is the mass of element X produced by massive stars per unit mass of star formation. The inferred O and Fe yields agree with predictions of the Sukhbold et al. CCSN models assuming their Z9.6+N20 neutrino-driven engine, a scenario in which many progenitors withM< 40M⊙implode to black holes rather than exploding. The yields are lower than assumed in many models of the galaxy mass–metallicity relation, reducing the level of outflows needed to match observed abundances. Our one-zone chemical evolution models with evolve to solar metallicity at late times. By further requiring that models reach [α/Fe] ≈ 0 at late times, we infer a Hubble-time integrated Type Ia supernova rate of , compatible with estimates from supernova surveys.more » « less
An official website of the United States government
