Abstract We perform three-dimensional shearing box hydrodynamical simulations to explore the outcome of gravitational instability in the outer regions of neutrino-cooled disks such as those formed from the collapse of rotating massive stars (“collapsars”). We employ a physical equation of state and optically thin neutrino cooling and assume an electron fraction set by the balance ofe±pair-capture reactions. Disks in a marginally stable initial state (Toomre parameterQ ≈ 1) undergo runaway cooling and fragmentation when the dimensionless cooling timescale obeysτcool ≡ tcoolΩ ≲ 10, where Ω is the orbital frequency; these conditions correspond to accretion rates ≳M⊙s−1on the upper end of those achieved by collapsar progenitor stars. Fragmentation leads to the formation of neutron-rich clumps (electron fractionYe ≲ 0.1) spanning a range of masses ∼0.01–1M⊙around the local Jeans value. Most clumps exceed the local Chandrasekhar mass and hence will continue to collapse to nuclear densities, forming neutron stars (NSs) with subsolar masses otherwise challenging to create through ordinary stellar core collapse. Even cool disks dominated byαparticles (Ye ≃ 0.5) can fragment and collapse into neutron-rich clumps capable of forming subsolar NSs. Although our simulations cannot follow this process directly, if the disk-formed NSs subsequently pair into binaries, the GW chirps from their rapid mergers are potentially detectable by ground-based observatories. The temporal coincidence of such a hierarchical NS merger chain with the collapsar gamma-ray burst and supernova would offer a uniquely spectacular multimessenger “symphony.”
more »
« less
Fragmentation in Gravitationally Unstable Collapsar Disks and Subsolar Neutron Star Mergers
Abstract Although stable neutron stars (NSs) can in principle exist down to massesMns≈ 0.1M⊙, standard models of stellar core-collapse predict a robust lower limitMns≳ 1.2M⊙, roughly commensurate with the Chandrasekhar massMChof the progenitor’s iron core (electron fractionYe≈ 0.5). However, this limit may be circumvented in sufficiently dense neutron-rich environments (Ye< 0.5) for which is reduced to ≲1M⊙. Such physical conditions could arise in the black hole accretion disks formed from the collapse of rapidly rotating stars (“collapsars”), as a result of gravitational instabilities and cooling-induced fragmentation, similar to models for planet formation in protostellar disks. We confirm that the conditions to form subsolar-mass NS (ssNS) may be marginally satisfied in the outer regions of massive neutrino-cooled collapsar disks. If the disk fragments into multiple ssNSs, their subsequent coalescence offers a channel for precipitating subsolar mass LIGO/Virgo gravitational-wave mergers that does not implicate primordial black holes. The model makes several additional predictions: (1) ∼Hz frequency Doppler modulation of the ssNS-merger gravitational-wave signals due to the binary’s orbital motion in the disk; (2) at least one additional gravitational-wave event (coincident within ≲hours), from the coalescence of the ssNS-merger remnant(s) with the central black hole; (3) an associated gamma-ray burst and supernova counterpart, the latter boosted in energy and enriched withr-process elements from the NS merger(s) embedded within the exploding stellar envelope (“kilonovae inside a supernova”).
more »
« less
- Award ID(s):
- 2002577
- PAR ID:
- 10538533
- Publisher / Repository:
- ApJL
- Date Published:
- Journal Name:
- The Astrophysical Journal Letters
- Volume:
- 971
- Issue:
- 2
- ISSN:
- 2041-8205
- Page Range / eLocation ID:
- L34
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract We report the observation of a coalescing compact binary with component masses 2.5–4.5M⊙and 1.2–2.0M⊙(all measurements quoted at the 90% credible level). The gravitational-wave signal GW230529_181500 was observed during the fourth observing run of the LIGO–Virgo–KAGRA detector network on 2023 May 29 by the LIGO Livingston observatory. The primary component of the source has a mass less than 5M⊙at 99% credibility. We cannot definitively determine from gravitational-wave data alone whether either component of the source is a neutron star or a black hole. However, given existing estimates of the maximum neutron star mass, we find the most probable interpretation of the source to be the coalescence of a neutron star with a black hole that has a mass between the most massive neutron stars and the least massive black holes observed in the Galaxy. We provisionally estimate a merger rate density of for compact binary coalescences with properties similar to the source of GW230529_181500; assuming that the source is a neutron star–black hole merger, GW230529_181500-like sources may make up the majority of neutron star–black hole coalescences. The discovery of this system implies an increase in the expected rate of neutron star–black hole mergers with electromagnetic counterparts and provides further evidence for compact objects existing within the purported lower mass gap.more » « less
-
Abstract On 2025 August 18, the LIGO–Virgo–KAGRA collaboration reported a subthreshold gravitational-wave candidate detection consistent with a subsolar-mass neutron star merger, denoted S250818k. An optical transient, AT2025ulz, was discovered within the localization region. AT2025ulz initially appeared to meet the expected behavior of kilonova emission, the telltale signature of a binary neutron star merger. The transient subsequently rebrightened after ∼5 days and was classified as a Type IIb supernova. In this work, we analyze the observations of its host galaxy obtained by the Dark Energy Spectroscopic Instrument (DESI). From the DESI spectrum, we obtain a secure redshift ofz = 0.084840 ± 0.000006. If S250818k has an astrophysical origin, this places the transient within 2σof the gravitational-wave distance and results in an integral overlap between the gravitational-wave alert and the transient location of . Our analysis of the host galaxy’s spectral energy distribution reveals a star-forming, dusty galaxy with stellar mass ∼1010M⊙, broadly consistent with the population of both short gamma-ray bursts and core-collapse supernova host galaxies. We also present our follow-up of DESI-selected candidate host galaxies using the Fraunhofer Telescope at the Wendelstein Observatory, and show the promise of DESI for associating or rejecting candidate electromagnetic counterparts to gravitational-wave alerts. These results emphasize the value of DESI’s extensive spectroscopic dataset in rapidly characterizing host galaxies, enabling spectroscopic host subtraction, and guiding targeted follow-up.more » « less
-
Abstract On 2023 November 23, the two LIGO observatories both detected GW231123, a gravitational-wave signal consistent with the merger of two black holes with masses and (90% credible intervals), at a luminosity distance of 0.7–4.1 Gpc, a redshift of , and with a network signal-to-noise ratio of ∼20.7. Both black holes exhibit high spins— and , respectively. A massive black hole remnant is supported by an independent ringdown analysis. Some properties of GW231123 are subject to large systematic uncertainties, as indicated by differences in the inferred parameters between signal models. The primary black hole lies within or above the theorized mass gap where black holes between 60–130M⊙should be rare, due to pair-instability mechanisms, while the secondary spans the gap. The observation of GW231123 therefore suggests the formation of black holes from channels beyond standard stellar collapse and that intermediate-mass black holes of mass ∼200M⊙form through gravitational-wave-driven mergers.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

