No Arabic abstract
Using the new state-of-the-art core-collapse supernova (CCSN) code F{sc{ornax}}, we have simulated the three-dimensional dynamical evolution of the cores of 9-, 10-, 11-, 12-, and 13-M$_{odot}$ stars from the onset of collapse. Stars from 8-M$_{odot}$ to 13-M$_{odot}$ constitute roughly 50% of all massive stars, so the explosive potential for this mass range is important to the overall theory of CCSNe. We find that the 9-, 10-, 11-, and 12-M$_{odot}$ models explode in 3D easily, but that the 13-M$_{odot}$ model does not. From these findings, and the fact that slightly more massive progenitors seem to explode citep{vartanyan2019}, we suggest that there is a gap in explodability near 12-M$_{odot}$ to 14-M$_{odot}$ for non-rotating progenitor stars. Factors conducive to explosion are turbulence behind the stalled shock, energy transfer due to neutrino-matter absorption and neutrino-matter scattering, many-body corrections to the neutrino-nucleon scattering rate, and the presence of a sharp silicon-oxygen interface in the progenitor. Our 3D exploding models frequently have a dipolar structure, with the two asymmetrical exploding lobes separated by a pinched waist where matter temporarily continues to accrete. This process maintains the driving neutrino luminosty, while partially shunting matter out of the way of the expanding lobes, thereby modestly facilitating explosion. The morphology of all 3D explosions is characterized by multiple bubble structures with a range of low-order harmonic modes. Though much remains to be done in CCSN theory, these and other results in the literature suggest that, at least for these lower-mass progenitors, supernova theory is converging on a credible solution.
We present results from an ab initio three-dimensional, multi-physics core collapse supernova simulation for the case of a 15 M progenitor. Our simulation includes multi-frequency neutrino transport with state-of-the-art neutrino interactions in the ray-by-ray approximation, and approximate general relativity. Our model exhibits a neutrino-driven explosion. The shock radius begins an outward trajectory at approximately 275 ms after bounce, giving the first indication of a developing explosion in the model. The onset of this shock expansion is delayed relative to our two-dimensional counterpart model, which begins at approximately 200 ms after core bounce. At a time of 441 ms after bounce, the angle-averaged shock radius in our three-dimensional model has reached 751 km. Further quantitative analysis of the outcomes in this model must await further development of the post-bounce dynamics and a simulation that will extend well beyond 1 s after stellar core bounce, based on the results for the same progenitor in the context of our two-dimensional, counterpart model. This more complete analysis will determine whether or not the explosion is robust and whether or not observables such as the explosion energy, 56Ni mass, etc. are in agreement with observations. Nonetheless, the onset of explosion in our ab initio three-dimensional multi-physics model with multi-frequency neutrino transport and general relativity is encouraging.
Most massive stars are born in binary or higher-order multiple systems and exchange mass with a companion during their lives. In particular, the progenitors of a large fraction of compact object mergers, and Galactic neutron stars (NSs) and black holes (BHs) have been stripped off their envelopes by a binary companion. Here, we study the evolution of single and stripped binary stars up to core collapse with the stellar evolution code MESA and their final fates with a parametric supernova (SN) model. We find that stripped binary stars can have systematically different pre-SN structures compared to genuine single stars and thus also different SN outcomes. The bases of these differences are already established by the end of core helium burning and are preserved up to core collapse. We find a non-monotonic pattern of NS and BH formation as a function of CO core mass that is different in single and stripped binary stars. In terms of initial masses, single stars of >35 Msun all form BHs, while this transition is only at 70 Msun in stripped stars. On average, stripped stars give rise to lower NS and BH masses, higher explosion energies, higher kick velocities and higher nickel yields. Within a simplified population synthesis model, we show that our results lead to a significant reduction of the rates of BH-NS and BH-BH mergers with respect to typical assumptions made on NS and BH formation. Therefore, we predict lower detection rates of such merger events by, e.g., advanced LIGO than is often considered. We further show how features in the NS-BH mass distribution of single and stripped stars relate to the chirp-mass distribution of compact object mergers. Further implications of our findings are discussed with respect to the missing red-supergiant problem, a possible mass gap between NSs and BHs, X-ray binaries and observationally inferred nickel masses from Type Ib/c and IIP Sne. [abridged]
Fallback in core-collapse supernovae is considered a major ingredient for explaining abundance anomalies in metal-poor stars and the natal kicks and spins of black holes (BHs). We present a first 3D simulation of BH formation and fallback in an aborted neutrino-driven explosion of a $40$ solar mass zero-metallicity progenitor from collapse to shock breakout. We follow the phase up to BH formation using the relativistic CoCoNuT-FMT code. For the subsequent evolution to shock breakout we apply the moving-mesh code Arepo to core-collapse supernovae for the first time. Our simulation shows that despite early BH formation, neutrino-heated bubbles can survive for tens of seconds before being accreted, leaving them sufficient time to transfer part of their energy to sustain the shock wave as is propagates through the envelope. Although the initial net energy ($sim 2$ Bethe) of the neutrino-heated ejecta barely equals the binding energy of the envelope, $11,mathrm{M}_odot$ of hydrogen are still expelled with an energy of $0.23$ Bethe. We find no significant mixing and only a modest BH kick and spin, but speculate that stronger effects could occur for slightly more energetic explosions or progenitors with less tightly bound envelopes.
We present a set of nonlocal thermodynamic equilibrium steady-state calculations of radiative transfer for one-year old type II supernovae (SNe) starting from state-of-the-art explosion models computed with detailed nucleosynthesis. This grid covers single-star progenitors with initial masses between 9 and 29$M_{odot}$, all evolved with KEPLER at solar metallicity and ignoring rotation. The [OI]$lambdalambda$$6300,6364$ line flux generally grows with progenitor mass, and H$alpha$ exhibits an equally strong and opposite trend. The [CaII]$lambdalambda$$7291,,7323$ strength increases at low $^{56}$Ni mass, low explosion energy, or with clumping. This CaII doublet, which forms primarily in the explosively-produced Si/S zones, depends little on the progenitor mass, but may strengthen if Ca$^+$ dominates in the H-rich emitting zones or if Ca is abundant in the O-rich zones. Indeed, Si-O shell merging prior to core collapse may boost the CaII doublet at the expense of the OI doublet, and may thus mimic the metal line strengths of a lower mass progenitor. We find that the $^{56}$Ni bubble effect has a weak impact, probably because it is too weak to induce much of an ionization shift in the various emitting zones. Our simulations compare favorably to observed SNe II, including SN2008bk (e.g., 9$M_{odot}$ model), SN2012aw (12$M_{odot}$ model), SN1987A (15$M_{odot}$ model), or SN2015bs (25$M_{odot}$ model with no Si-O shell merging). SNe II with narrow lines and a low $^{56}$Ni mass are well matched by the weak explosion of 9$-$11$M_{odot}$ progenitors. The nebular-phase spectra of standard SNe II can be explained with progenitors in the mass range 12$-$15$M_{odot}$, with one notable exception for SN2015bs. In the intermediate mass range, these mass estimates may increase by a few $M_{odot}$ with allowance for clumping of the O-rich material or CO molecular cooling.
We previously reported on a recursive method to generate the expansion of the lattice Green function of the $d$-dimensional face-centred cubic lattice (fcc). The method was used to generate many coefficients for d=7 and the corresponding linear differential equation has been obtained. In this paper, we show the strength and the limit of the method by producing the series and the corresponding linear differential equations for d=8, 9, 10, 11, 12. The differential Galois groups of these linear differential equations are shown to be symplectic for d=8, 10, 12 and orthogonal for d= 9, 11. The recursion relation naturally provides a 2-dimensional array $ T_d(n,j)$ where only the coefficients $ t_d(n,0)$ correspond to the coefficients of the lattice Green function of the d-dimensional fcc. The coefficients $ t_d(n,j)$ are associated to D-finite bivariate series annihilated by linear partial differential equations that we analyze.