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 more » 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.

« less
Award ID(s):
Publication Date:
Journal Name:
The Astrophysical Journal
Page Range or eLocation-ID:
Article No. 73
DOI PREFIX: 10.3847
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    The tidal disruption of stars by supermassive black holes (SMBHs) probes relativistic gravity. In the coming decade, the number of observed tidal disruption events (TDEs) will grow by several orders of magnitude, allowing statistical inferences of the properties of the SMBH and stellar populations. Here we analyze the probability distribution functions of the pericenter distances of stars that encounter an SMBH in the Schwarzschild geometry, where the results are completely analytic, and the Kerr metric. From this analysis we calculate the number of observable TDEs, defined to be those that come within the tidal radiusrtbut outside the direct capture radius (which is, in general, larger than the horizon radius). We find that relativistic effects result in a steep decline in the number of stars that have pericenter distancesrp≲ 10rg, whererg=GM/c2, and that for maximally spinning SMBHs the distribution function ofrpat such distances scales asfrprp4/3, or in terms ofβrt/rpscales asfββ−10/3. We find that spin has little effect on the TDE fraction until the very-high-mass end, where instead of being identically zero the rate is small (≲1% of the expected rate in the absence of relativistic effects). Effectively independent of spin, if the progenitorsmore »of TDEs reflect the predominantly low-mass stellar population and thus have masses ≲1M, we expect a substantial reduction in the rate of TDEs above 107M.

    « less
  2. Abstract

    We present a stellar dynamical mass measurement of a newly detected supermassive black hole (SMBH) at the center of the fast-rotating, massive elliptical galaxy NGC 2693 as part of the MASSIVE survey. We combine high signal-to-noise ratio integral field spectroscopy (IFS) from the Gemini Multi-Object Spectrograph with wide-field data from the Mitchell Spectrograph at McDonald Observatory to extract and model stellar kinematics of NGC 2693 from the central ∼150 pc out to ∼2.5 effective radii. Observations from Hubble Space Telescope WFC3 are used to determine the stellar light distribution. We perform fully triaxial Schwarzschild orbit modeling using the latest TriOS code and a Bayesian search in 6D galaxy model parameter space to determine NGC 2693's SMBH mass (MBH), stellar mass-to-light ratio, dark matter content, and intrinsic shape. We findMBH=1.7±0.4×109Mand a triaxial intrinsic shape with axis ratiosp=b/a= 0.902 ± 0.009 andq=c/a=0.7210.010+0.011, triaxiality parameterT= 0.39 ± 0.04. In comparison, the best-fit orbit model in the axisymmetric limit and (cylindrical) Jeans anisotropic model of NGC 2693 preferMBH=2.4±0.6×109MandMBH=2.9±0.3×109M, respectively. Neither model can account for the non-axisymmetric stellar velocity features present inmore »the IFS data.

    « less
  3. Abstract

    The existence of primordial black holes (PBHs), which may form from the collapse of matter overdensities shortly after the Big Bang, is still under debate. Among the potential signatures of PBHs are gravitational waves (GWs) emitted from binary black hole (BBH) mergers at redshiftsz≳ 30, where the formation of astrophysical black holes is unlikely. Future ground-based GW detectors, the Cosmic Explorer and Einstein Telescope, will be able to observe equal-mass BBH mergers with total mass of(10100)Mat such distances. In this work, we investigate whether the redshift measurement of a single BBH source can be precise enough to establish its primordial origin. We simulate BBHs of different masses, mass ratios and orbital orientations. We show that for BBHs with total masses between 20Mand 40Mmerging atz≥ 40, one can inferz> 30 at up to 97% credibility, with a network of one Einstein Telescope, one 40 km Cosmic Explorer in the US, and one 20 km Cosmic Explorer in Australia. This number reduces to 94% with a smaller network made of one Einstein Telescope and one 40 km Cosmic Explorer in the US. We also analyze how the measurement depends on the Bayesian priors usedmore »in the analysis and verify that priors that strongly favor the wrong model yield smaller Bayesian evidences.

    « less
  4. Abstract

    We present a Keck/MOSFIRE rest-optical composite spectrum of 16 typical gravitationally lensed star-forming dwarf galaxies at 1.7 ≲z≲ 2.6 (zmean= 2.30), all chosen independent of emission-line strength. These galaxies have a median stellar mass oflog(M*/M)med=8.290.43+0.51and a median star formation rate ofSFRHαmed=2.251.26+2.15Myr1. We measure the faint electron-temperature-sensitive [Oiii]λ4363 emission line at 2.5σ(4.1σ) significance when considering a bootstrapped (statistical-only) uncertainty spectrum. This yields a direct-method oxygen abundance of12+log(O/H)direct=7.880.22+0.25(0.150.06+0.12Z). We investigate the applicability at highzof locally calibrated oxygen-based strong-line metallicity relations, finding that the local reference calibrations of Bian et al. best reproduce (≲0.12 dex) our composite metallicity at fixed strong-line ratio. At fixedM*, our composite is well represented by thez∼ 2.3 direct-method stellar mass—gas-phase metallicity relation (MZR) of Sanders et al. When comparing to predicted MZRs from the IllustrisTNG and FIRE simulations, having recalculated our stellar masses with more realistic nonparametric star formation histories(log(M*/M)med=8.920.22+0.31), we find excellent agreement with the FIRE MZR. Our composite is consistent with no metallicity evolution, atmore »fixedM*and SFR, of the locally defined fundamental metallicity relation. We measure the doublet ratio [Oii]λ3729/[Oii]λ3726 = 1.56 ± 0.32 (1.51 ± 0.12) and a corresponding electron density ofne=10+215cm3(ne=10+74cm3) when considering the bootstrapped (statistical-only) error spectrum. This result suggests that lower-mass galaxies have lower densities than higher-mass galaxies atz∼ 2.

    « less
  5. Abstract

    We present the results of a new reverberation mapping campaign for the broad-line active galactic nucleus (AGN) in the edge-on spiral IC 4329A. Monitoring of the optical continuum withV-band photometry and broad emission-line flux variability with moderate-resolution spectroscopy allowed emission-line light curves to be measured for Hβ, Hγ, and Heiiλ4686. We find a time delay of16.32.3+2.6days for Hβ, a similar time delay of16.02.6+4.8days for Hγ, and an unresolved time delay of0.63.9+3.9days for Heii. The time delay for Hβis consistent with the predicted value from the relationship between AGN luminosity and broad-line region radius, after correction for the ∼2.4 mag of intrinsic extinction at 5100 Å. Combining the measured time delay for Hβwith the broad emission-line width and an adopted value of 〈f〉 = 4.8, we find a central supermassive black hole mass ofMBH=6.81.1+1.2×107M. Velocity-resolved time delays were measured across the broad Hβemission-line profile and may be consistent with an “M”-like shape. Modeling of the full reverberation response of Hβwas able to provide only modest constraints on some parameters, but does exhibit agreement with the black hole mass and average time delay. The models also suggest that themore »AGN structure is misaligned by a large amount from the edge-on galaxy disk. This is consistent with expectations from the unified model of AGNs, in which broad emission lines are expected to be visible only for AGNs that are viewed at relatively face-on inclinations.

    « less