Gravitational-wave (GW) detections of binary black hole (BH) mergers have begun to sample the cosmic BH mass distribution. The evolution of single stellar cores predicts a gap in the BH mass distribution due to pair-instability supernovae (PISNe). Determining the upper and lower edges of the BH mass gap can be useful for interpreting GW detections of merging BHs. We use
Long-duration gamma-ray bursts (lGRBs) originate in relativistic collimated outflows—jets—that drill their way out of collapsing massive stars. Accurately modeling this process requires realistic stellar profiles for the jets to propagate through and break out of. Most previous studies have used simple power laws or pre-collapse models for massive stars. However, the relevant stellar profile for lGRB models is in fact that of a star after its core has collapsed to form a compact object. To self-consistently compute such a stellar profile, we use the open-source code GR1D to simulate the core-collapse process for a suite of low-metallicity rotating massive stellar progenitors that have undergone chemically homogeneous evolution. Our models span a range of zero-age main-sequence (ZAMS) masses:
- Award ID(s):
- 1815304
- NSF-PAR ID:
- 10397623
- Publisher / Repository:
- DOI PREFIX: 10.3847
- Date Published:
- Journal Name:
- The Astrophysical Journal Letters
- Volume:
- 944
- Issue:
- 2
- ISSN:
- 2041-8205
- Format(s):
- Medium: X Size: Article No. L38
- Size(s):
- ["Article No. L38"]
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
Abstract MESA to evolve single, nonrotating, massive helium cores with a metallicity ofZ = 10−5, until they either collapse to form a BH or explode as a PISN, without leaving a compact remnant. We calculate the boundaries of the lower BH mass gap for S-factors in the range S(300 keV) = (77,203) keV b, corresponding to the ±3σ uncertainty in our high-resolution tabulated12C(α ,γ )16O reaction rate probability distribution function. We extensively test temporal and spatial resolutions for resolving the theoretical peak of the BH mass spectrum across the BH mass gap. We explore the convergence with respect to convective mixing and nuclear burning, finding that significant time resolution is needed to achieve convergence. We also test adopting a minimum diffusion coefficient to help lower-resolution models reach convergence. We establish a new lower edge of the upper mass gap asM lower≃M ⊙from the ±3σ uncertainty in the12C(α ,γ )16O rate. We explore the effect of a larger 3α rate on the lower edge of the upper mass gap, findingM lower≃M ⊙. We compare our results with BHs reported in the Gravitational-Wave Transient Catalog. -
more » « less
Abstract We present extensive multifrequency Karl G. Jansky Very Large Array (VLA) and Very Long Baseline Array (VLBA) observations of the radio-bright supernova (SN) IIb SN 2004C that span ∼40–2793 days post-explosion. We interpret the temporal evolution of the radio spectral energy distribution in the context of synchrotron self-absorbed emission from the explosion’s forward shock as it expands in the circumstellar medium (CSM) previously sculpted by the mass-loss history of the stellar progenitor. VLBA observations and modeling of the VLA data point to a blastwave with average velocity ∼0.06
c that carries an energy of ≈1049erg. Our modeling further reveals a flat CSM density profileρ CSM∝R −0.03±0.22up to a break radiusR br≈ (1.96 ± 0.10) × 1016cm, with a steep density gradient followingρ CSM∝R −2.3±0.5at larger radii. We infer that the flat part of the density profile corresponds to a CSM shell with mass ∼0.021M ☉, and that the progenitor’s effective mass-loss rate varied with time over the range (50–500) × 10−5M ☉yr−1for an adopted wind velocityv w = 1000 km s−1and shock microphysical parametersϵ e = 0.1,ϵ B = 0.01. These results add to the mounting observational evidence for departures from the traditional single-wind mass-loss scenarios in evolved, massive stars in the centuries leading up to core collapse. Potentially viable scenarios include mass loss powered by gravity waves and/or interaction with a binary companion. -
more » « less
Abstract We compare the core-collapse evolution of a pair of 15.8
M ☉stars with significantly different internal structures, a consequence of the bimodal variability exhibited by massive stars during their late evolutionary stages. The 15.78 and 15.79M ☉progenitors have core masses (masses interior to an entropy of 4k Bbaryon−1) of 1.47 and 1.78M ☉and compactness parametersξ 1.75of 0.302 and 0.604, respectively. The core-collapse simulations are carried out in 2D to nearly 3 s postbounce and show substantial differences in the times of shock revival and explosion energies. The 15.78M ☉model begins exploding promptly at 120 ms postbounce when a strong density decrement at the Si–Si/O shell interface, not present in the 15.79M ☉progenitor, encounters the stalled shock. The 15.79M ☉model takes 100 ms longer to explode but ultimately produces a more powerful explosion. Both the larger mass accretion rate and the more massive core of the 15.79M ☉model during the first 0.8 s postbounce time result in largerν e /M ☉model resulted in the ejection of twice as much56Ni. Most of the ejecta in both models are moderately proton rich, though counterintuitively the highest electron fraction (Y e = 0.61) ejecta in either model are in the less-energetic 15.78M ☉model, while the lowest electron fraction (Y e = 0.45) ejecta in either model are in the 15.79M ☉model. -
more » « less
Abstract The fate of stars in the zero-age main-sequence (ZAMS) range ≈8–12
M ⊙is unclear. They could evolve to form white dwarfs or explode as electron-capture supernovae (SNe) or iron core-collapse SNe (CCSNe). Even though the initial mass function indicates that this mass range should account for over 40% of all CCSN progenitors, few have been observationally confirmed, likely due to the faintness and rapid evolution of some of these transients. In this paper, we present a sample of nine Ca-rich/O-poor Type IIb SNe detected by the Zwicky Transient Facility with progenitors likely in this mass range. These sources have a [Caii ]λ λ 7291, 7324/[Oi ]λ λ 6300, 6364 flux ratio of ≳2 in their nebular spectra. Comparing the measured [Oi ] luminosity (≲1039erg s−1) and derived oxygen mass (≈0.01M ⊙) with theoretical models, we infer that the progenitor ZAMS mass for these explosions is less than 12M ⊙. The ejecta properties (M ej≲ 1M ⊙andE kin∼ 1050erg) are also consistent. The low ejecta mass of these sources indicates a class of strongly-stripped SNe that is a transition between the regular stripped-envelope SNe and ultra-stripped SNe. The progenitor could be stripped by a main-sequence companion and result in the formation of a neutron star−main sequence binary. Such binaries have been suggested to be progenitors of neutron star−white dwarf systems that could merge within a Hubble time and be detectable with LISA. -
more » « less
Abstract Many of the short-lived radioactive nuclei that were present in the early solar system can be produced in massive stars. In the first paper in this series, we focused on the production of26Al in massive binaries. In our second paper, we considered rotating single stars; two more short-lived radioactive nuclei,36Cl and41Ca; and the comparison to the early solar system data. In this work, we update our previous conclusions by further considering the impact of binary interactions. We used the MESA stellar evolution code with an extended nuclear network to compute massive (10–80
M ⊙), binary stars at various initial periods and solar metallicity (Z = 0.014), up to the onset of core collapse. The early solar system abundances of26Al and41Ca can be matched self-consistently by models with initial masses ≥25M ⊙, while models with initial primary masses ≥35M ⊙can also match36Cl. Almost none of the models provide positive net yields for19F, while for22Ne the net yields are positive from 30M ⊙and higher. This leads to an increase by a factor of approximately 4 in the amount of22Ne produced by a stellar population of binary stars, relative to single stars. In addition, besides the impact on the stellar yields, our 10M ⊙primary star undergoing Case A mass transfer ends its life as a white dwarf instead of as a core-collapse supernova. This demonstrates that binary interactions can also strongly impact the evolution of stars close to the supernova boundary.
