Black Hole Supernovae Outcomes Across a Wide Progenitor Range

TL;DR

This study investigates the outcomes of black hole supernovae across a broad range of progenitor star structures using long-term simulations, revealing their occurrence in various progenitors and detailed remnant properties.

Contribution

It provides the first extensive set of long-term simulations showing BHSNe arise across a wide progenitor mass and structure spectrum, with detailed analysis of remnant masses and explosion energies.

Findings

18 BHSN outcomes found across the progenitor mass range.

BH formation occurs between 0.7 s and 4.4 s after bounce.

Final explosion energies range from 2×10^{49} to 3×10^{51} erg.

Abstract

Black hole supernovae (BHSNe), the term we use for core-collapse events in which black hole (BH) formation occurs after shock revival but before the explosion is complete, have emerged as a natural outcome of multidimensional simulations as these calculations have been extended to seconds after bounce. Yet they remain one of the least studied outcomes of core collapse. Here, we assess whether they are confined to the most compact and massive progenitors, whose birth rates are low, or whether they arise systematically across a wider range of progenitor structures. We perform 23 long-term axisymmetric core-collapse simulations of progenitors spanning 19.51-60$\,M_\odot$ and compactnesses $0.31 \lesssim \xi_{2.5} \lesssim 0.63$. We find 18 BHSN outcomes across nearly the full ZAMS mass range considered, corresponding to progenitors with $0.40 \lesssim \xi_{2.5} \lesssim 0.63$. BH formation occurs between $\sim0.7$ s and $\sim4.4$ s after bounce. After BH formation, we continue the evolution with an excision treatment to at least 5000 s. The final explosion energies span $\sim2\times10^{49}$-$3\times10^{51}$ erg, while the final BH gravitational masses span $\sim3$-$26\,M_\odot$. We find a clear remnant-mass trend with CO-core mass, but show that the CO core alone is not an adequate proxy for the final BH mass, especially for progenitors at the low- and high-mass ends of the CO-core distribution. Except for the highest CO-core mass models, no single spherical mass coordinate cleanly separates ejecta from remnant material. Finally, a 2D axisymmetric and a 3D model are compared as we discuss differences between the two geometries.