No Arabic abstract
Most massive stars, the progenitors of core-collapse supernovae, are in close binary systems and may interact with their companion through mass transfer or merging. We undertake a population synthesis study to compute the delay-time distribution of core-collapse supernovae, that is, the supernova rate versus time following a starburst, taking into account binary interactions. We test the systematic robustness of our results by running various simulations to account for the uncertainties in our standard assumptions. We find that a significant fraction, $15^{+9}_{-8}$%, of core-collapse supernovae are `late, that is, they occur 50-200 Myrs after birth, when all massive single stars have already exploded. These late events originate predominantly from binary systems with at least one, or, in most cases, with both stars initially being of intermediate mass ($4-8M_{odot}$). The main evolutionary channels that contribute often involve either the merging of the initially more massive primary star with its companion or the engulfment of the remaining core of the primary by the expanding secondary that has accreted mass at an earlier evolutionary stage. Also, the total number of core-collapse supernovae increases by $14^{+15}_{-14}$% because of binarity for the same initial stellar mass. The high rate implies that we should have already observed such late core-collapse supernovae, but have not recognized them as such. We argue that $phi$ Persei is a likely progenitor and that eccentric neutron star - white dwarf systems are likely descendants. Late events can help explain the discrepancy in the delay-time distributions derived from supernova remnants in the Magellanic Clouds and extragalactic type Ia events, lowering the contribution of prompt Ia events. We discuss ways to test these predictions and speculate on the implications for supernova feedback in simulations of galaxy evolution.
Ground-based optical spectra and Hubble Space Telescope images of ten core-collapse supernovae (CCSNe) obtained several years to decades after outburst are analyzed with the aim of understanding the general properties of their late-time emissions. New observations of SN 1957D, 1970G, 1980K, and 1993J are included as part of the study. Blueshifted line emissions in oxygen and/or hydrogen with conspicuous line substructure are a common and long-lasting phenomenon in the late-time spectra. Followed through multiple epochs, changes in the relative strengths and velocity widths of the emission lines are consistent with expectations for emissions produced by interaction between SN ejecta and the progenitor stars circumstellar material. The most distinct trend is an increase in the strength of [O III]/([O I]+[O II]) with age, and a decline in Halpha/([O I]+[O II]) which is broadly consistent with the view that the reverse shock has passed through the H envelope of the ejecta in many of these objects. We also present a spatially integrated spectrum of the young Galactic supernova remnant Cassiopeia A (Cas A). Similarities observed between the emission line profiles of the 330 yr old Cas A remnant and decades old CCSNe suggest that observed emission line asymmetry in evolved CCSN spectra may be associated with dust in the ejecta, and that minor peak substructure typically interpreted as clumps or blobs of ejecta may instead be linked with large-scale rings of SN debris.
The importance of detecting neutrinos from a Milky Way core-collapse supernova is well known. An under-studied phase is proto-neutron star cooling. For SN 1987A, this seemingly began at about 2 s, and is thus probed by only 6 of the 19 events (and only the $bar{ u}_e$ flavor) in the Kamiokande-II and IMB detectors. With the higher statistics expected for present and near-future detectors, it should be possible to measure detailed neutrino signals out to very late times. We present the first comprehensive study of neutrino detection during the proto-neutron star cooling phase, considering a variety of outcomes, using all flavors, and employing detailed detector physics. For our nominal model, the event yields (at 10 kpc) after 10 s -- the approximate duration of the SN 1987A signal -- far exceed the entire SN 1987A yield, with $simeq$250 $bar{ u}_e$ events (to 50 s) in Super-Kamiokande, $simeq$110 $ u_e$ events (to 40 s) in DUNE, and $simeq$10 $ u_mu, u_tau, bar{ u}_mu, bar{ u}_tau$ events (to 20 s) in JUNO. These data would allow unprecedented probes of the proto-neutron star, including the onset of neutrino transparency and hence its transition to a neutron star. If a black hole forms, even at very late times, this can be clearly identified. But will the detectors fulfill their potential for this perhaps once-ever opportunity for an all-flavor, high-statistics detection of a core collapse? Maybe. Further work is urgently needed, especially for DUNE to thoroughly investigate and improve its MeV capabilities.
Hydrogen-rich supernovae, known as Type II (SNe II), are the most common class of explosions observed following the collapse of the core of massive stars. We use analytical estimates and population synthesis simulations to assess the fraction of SNe II progenitors that are expected to have exchanged mass with a companion prior to explosion. We estimate that 1/3 to 1/2 of SN II progenitors have a history of mass exchange with a binary companion before exploding. The dominant binary channels leading to SN II progenitors involve the merger of binary stars. Mergers are expected to produce a diversity of SN II progenitor characteristics, depending on the evolutionary timing and properties of the merger. Alternatively, SN II progenitors from interacting binaries may have accreted mass from their companion, and subsequently been ejected from the binary system after their companion exploded. We show that the overall fraction of SN II progenitors that are predicted to have experienced binary interaction is robust against the main physical uncertainties in our models. However, the relative importance of different binary evolutionary channels is affected by changing physical assumptions. We further discuss ways in which binarity might contribute to the observed diversity of SNe II by considering potential observational signatures arising from each binary channel. For supernovae which have a substantial H-rich envelope at explosion (i.e., excluding Type IIb SNe), a surviving non-compact companion would typically indicate that the supernova progenitor star was in a wide, non-interacting binary. We argue that a significant fraction of even Type II-P SNe are expected to have gained mass from a companion prior to explosion.
We infer the progenitor mass distribution for 22 historic core-collapse supernovae (CCSNe) using a Bayesian hierarchical model. For this inference, we use the local star formation histories to estimate the age for each supernova (SN). These star formation histories often show multiple bursts of star formation; our model assumes that one burst is associated with the SN progenitor and the others are random bursts of star formation. The primary inference is the progenitor age distribution. Due to the limited number of historic SNe and highly uncertain star formation at young ages, we restrict our inference to the slope of the age distribution and the maximum age for CCSNe. Using single-star evolutionary models, we transform the progenitor age distribution into a progenitor mass distribution. Under these assumptions, the minimum mass for CCSNe is ${M_textrm{min}}~=~8.60^{+0.37}_{-0.41} M_odot$ and the slope of the progenitor mass distribution is $alpha = -2.61^{+1.05}_{-1.18}$. The power-law slope for the progenitor mass distribution is consistent with the standard Salpeter initial mass function ($alpha = -2.35$). These values are consistent with previous estimates using precursor imaging and the age-dating technique, further confirming that using stellar populations around SN and supernova remnants is a reliable way to infer the progenitor masses.
Convective instabilities in the advanced stages of nuclear shell burning can play an important role in neutrino-driven supernova explosions. In our previous work, we studied the interaction of vorticity and entropy waves with the supernova shock using a linear perturbations theory. In this paper, we extend our work by studying the effect of acoustic waves. As the acoustic waves cross the shock, the perturbed shock induces a field of entropy and vorticity waves in the post-shock flow. We find that, even when the upstream flow is assumed to be dominated by sonic perturbations, the shock-generated vorticity waves contain most of the turbulent kinetic energy in the post-shock region, while the entropy waves produced behind the shock are responsible for most of the density perturbations. The entropy perturbations are expected to become buoyant as a response to the gravity force and then generate additional turbulence in the post-shock region. This leads to a modest reduction of the critical neutrino luminosity necessary for producing an explosion, which we estimate to be less than $ sim 5 % $.