Explosions from Rotating Very Massive Star Collapses to Black Holes: Effects of Nuclear Burning

Abstract

We investigate the collapse of rotating very massive and supermassive stellar cores using numerical relativity simulations including an alpha-chain nuclear reaction network and neutrino cooling. Our main survey focuses on newly constructed models with initial core masses of $2 \times {10}^{3}-5\times 10^4M_\odot$. The collapse is triggered either by pair instability in lower-mass cores or by general-relativistic instability in higher-mass cores. We find that higher-mass cores undergo a nearly homologous collapse, whereas lower-mass cores show a more runaway-like collapse because neutrino cooling becomes more efficient at their higher densities and temperatures. Consequently, the black hole formed in lower-mass models initially contains a smaller fraction of the core mass, and disk formation occurs while a larger amount of matter remains outside the black hole. The lower compactness of pair-unstable cores also allows larger dimensionless angular momentum, favoring the formation of rapidly rotating black holes and massive disks. The disk bounce drives mass ejection with ejecta masses of order $10-10^3M_\odot$ and kinetic energies of order $10^{53}-10^{55}\,\mathrm{erg}$. Significant $^{56}$Ni production in the disk-bounce ejecta occurs only in the lowest-mass models. For selected models, we further follow the viscous evolution of the disk and find that viscosity enhances the ejecta mass and kinetic energy. In models with $\lesssim10^4M_\odot$, the viscosity-driven ejecta can originate from disk matter that has reached nuclear statistical equilibrium and can therefore become rich in $^{56}$Ni. These results suggest that rotating very massive star collapses can produce massive, energetic ejecta and, for sufficiently low core masses, substantial iron-group elements.