skip to main content

Title: Inferring the Neutron Star Maximum Mass and Lower Mass Gap in Neutron Star–Black Hole Systems with Spin

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.6Mobject 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 is2.70.4+0.5M(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 at5.41.0+0.7M. 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
Award ID(s):
Author(s) / Creator(s):
Publisher / Repository:
DOI PREFIX: 10.3847
Date Published:
Journal Name:
The Astrophysical Journal
Medium: X Size: Article No. 73
Article No. 73
Sponsoring Org:
National Science Foundation
More Like this
  1. 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 below2.40.5+0.5M, 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 at2.20.5+0.7Mand6.01.4+2.4M. 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
  2. Abstract

    We measure the correlation between black hole massMBHand host stellar massM*for a sample of 38 broad-line quasars at 0.2 ≲z≲ 0.8 (median redshiftzmed= 0.5). The black hole masses are derived from a dedicated reverberation mapping program for distant quasars, and the stellar masses are derived from two-band optical+IR Hubble Space Telescope imaging. Most of these quasars are well centered within ≲1 kpc from the host galaxy centroid, with only a few cases in merging/disturbed systems showing larger spatial offsets. Our sample spans two orders of magnitude in stellar mass (∼109–1011M) and black hole mass (∼107–109M) and reveals a significant correlation between the two quantities. We find a best-fit intrinsic (i.e., selection effects corrected)MBHM*,hostrelation oflog(MBH/M)=7.010.33+0.23+1.740.64+0.64log(M*,host/1010M), with an intrinsic scatter of0.470.17+0.24dex. Decomposing our quasar hosts into bulges and disks, there is a similarMBHM*,bulgerelation with slightly larger scatter, likely caused by systematic uncertainties in the bulge–disk decomposition. TheMBHM*,hostrelation atzmed= 0.5 is similar to that in local quiescent galaxies, with negligible evolution over the redshift range probed by our sample. With direct black hole masses from reverberation mapping and the large dynamical range of the sample, selection biases do not appear to affect our conclusions significantly. Our results, along with other samples in the literature, suggest that the locally measured black hole mass–host stellar mass relation is already in place atz∼ 1.

    more » « less
  3. Abstract

    We conduct a systematic tidal disruption event (TDE) demographics analysis using the largest sample of optically selected TDEs. A flux-limited, spectroscopically complete sample of 33 TDEs is constructed using the Zwicky Transient Facility over 3 yr (from 2018 October to 2021 September). We infer the black hole (BH) mass (MBH) with host galaxy scaling relations, showing that the sampleMBHranges from 105.1Mto 108.2M. We developed a survey efficiency corrected maximum volume method to infer the rates. The rest-frameg-band luminosity function can be well described by a broken power law ofϕ(Lg)Lg/Lbk0.3+Lg/Lbk2.61, withLbk= 1043.1erg s−1. In the BH mass regime of 105.3≲ (MBH/M) ≲ 107.3, the TDE mass function followsϕ(MBH)MBH0.25, which favors a flat local BH mass function (dnBH/dlogMBHconstant). We confirm the significant rate suppression at the high-mass end (MBH≳ 107.5M), which is consistent with theoretical predictions considering direct capture of hydrogen-burning stars by the event horizon. At a host galaxy mass ofMgal∼ 1010M, the average optical TDE rate is ≈3.2 × 10−5galaxy−1yr−1. We constrain the optical TDE rate to be [3.7, 7.4, and 1.6] × 10−5galaxy−1yr−1in galaxies with red, green, and blue colors.

    more » « less
  4. Abstract

    We present optical follow-up of IGR J16194-2810, a hard X-ray source discovered by the INTEGRAL mission. The optical counterpart is a ∼500Lred giant at a distance of 2.1 kpc. We measured 17 radial velocities (RVs) of the giant over a period of 271 days. Fitting these RVs with a Keplerian model, we find an orbital period ofPorb= 192.73 ± 0.01 days and a companion mass functionf(M2) = 0.365 ± 0.003M. We detect ellipsoidal variability with the same period in optical light curves from the ASAS-SN survey. Joint fitting of the RVs, light curves, and the broadband spectral energy distribution allows us to robustly constrain the masses of both components. We find a giant mass ofM=0.990.03+0.02Mand a companion mass ofM2=1.230.03+0.05M, implying that the companion is a neutron star (NS). We recover a 4.06 hr period in the system’s TESS light curve, which we tentatively associate with the NS spin period. The giant does not yet fill its Roche lobe, suggesting that current mass transfer is primarily via winds. Modules for Experiments in Stellar Astrophysics evolutionary models predict that the giant will overflow its Roche lobe in 5–10 Myr, eventually forming a recycled pulsar + white dwarf binary with a ∼900 days period. IGR J16194-2810 provides a window on the future evolution of wide NS + main sequence binaries recently discovered via Gaia astrometry. As with those systems, the binary’s formation history is uncertain. Before the formation of the NS, it likely survived a common envelope episode with a donor-to-accretor mass ratio ≳10 and emerged in a wide orbit. The NS likely formed with a weak kick (vkick≲ 50 km s−1), as stronger kicks would have disrupted the orbit.

    more » « less
  5. Abstract

    The cluster mass–richness relation (MRR) is an observationally efficient and potentially powerful cosmological tool for constraining the matter density Ωmand the amplitude of fluctuationsσ8using the cluster abundance technique. We derive the MRR relation usingGalWCat19, a publicly available galaxy cluster catalog we created from the Sloan Digital Sky Survey-DR13 spectroscopic data set. In the MRR, cluster mass scales with richness aslogM200=α+βlogN200. We find that the MRR we derive is consistent with both the IllustrisTNG and mini-Uchuu cosmological numerical simulations, with a slope ofβ≈ 1. We use the MRR we derived to estimate cluster masses from theGalWCat19catalog, which we then use to set constraints on Ωmandσ8. Utilizing the all-member MRR, we obtain constraints of Ωm=0.310.03+0.04andσ8=0.820.04+0.05, and utilizing the red member MRR only, we obtain Ωm=0.310.03+0.04andσ8=0.810.04+0.05. Our constraints on Ωmandσ8are consistent and very competitive with the Planck 2018 results.

    more » « less