Search

Six binary-merger progenitors of Supernova 1987A (SN 1987A) with properties close to those of the blue supergiant Sanduleak -69 202 are exploded by neutrino heating and evolved until long after shock breakout in three dimensions (3D), and continued for light-curve calculations in spherical symmetry. Our results confirm previous findings for single-star progenitors: (1) 3D neutrino-driven explosions with SN 1987A-like energies synthesize Ni-56 masses consistent with the radioactive light-curve tail; (2) hydrodynamic models mix hydrogen inward to minimum velocities below 40 km/s compatible with spectral observations of SN 1987A; and (3) for given explosion energy the efficiency of outward radioactive Ni-56 mixing depends mainly on high growth factors of Rayleigh-Taylor instabilities at the (C+O)/He and He/H composition interfaces and a weak interaction of fast plumes with the reverse shock occurring below the He/H interface. All binary-merger models possess presupernova radii matching the photometric radius of Sanduleak -69 202 and a structure of the outer layers allowing them to reproduce the observed initial luminosity peak in the first about 7 days. Models that mix about 0.5 Msun of hydrogen into the He-shell and exhibit strong outward mixing of Ni-56 with maximum velocities exceeding the 3000 km/s observed for the bulk of ejected Ni-56 have light-curve shapes in good agreement with the dome of the SN 1987A light curve. A comparative analysis of the best representatives of our 3D neutrino-driven explosion models of SN 1987A based on single-star and binary-merger progenitors reveals that only one binary model fulfills all observational constraints, except one., Comment: 31 pages, 14 figures, 6 tables. Accepted for publication in ApJ

We present 3D full-sphere supernova simulations of non-rotating low-mass (~9 Msun) progenitors, covering the entire evolution from core collapse through bounce and shock revival, through shock breakout from the stellar surface, until fallback is completed several days later. We obtain low-energy explosions [~(0.5-1.0)x 10^{50} erg] of iron-core progenitors at the low-mass end of the core-collapse supernova (LMCCSN) domain and compare to a super-AGB (sAGB) progenitor with an oxygen-neon-magnesium core that collapses and explodes as electron-capture supernova (ECSN). The onset of the explosion in the LMCCSN models is modelled self-consistently using the Vertex-Prometheus code, whereas the ECSN explosion is modelled using parametric neutrino transport in the Prometheus-HOTB code, choosing different explosion energies in the range of previous self-consistent models. The sAGB and LMCCSN progenitors that share structural similarities have almost spherical explosions with little metal mixing into the hydrogen envelope. A LMCCSN with less 2nd dredge-up results in a highly asymmetric explosion. It shows efficient mixing and dramatic shock deceleration in the extended hydrogen envelope. Both properties allow fast nickel plumes to catch up with the shock, leading to extreme shock deformation and aspherical shock breakout. Fallback masses of <~5x10^{-3} Msun have no significant effects on the neutron star (NS) masses and kicks. The anisotropic fallback carries considerable angular momentum, however, and determines the spin of the newly-born NS. The LMCCSNe model with less 2nd dredge-up results in a hydrodynamic and neutrino-induced NS kick of >40 km/s and a NS spin period of ~30 ms, both not largely different from those of the Crab pulsar at birth., Comment: 47 pages, 27 figures, 6 tables; minor revisions, accepted by MNRAS

With the same method as used previously, we investigate neutrino-driven explosions of a larger sample of blue supergiant models. The larger sample includes three new presupernova stars. The results are compared with light-curve observations of the peculiar type IIP SN 1987A. The explosions were modeled in 3D with the neutrino-hydrodynamics code PROMETHEUS-HOTB, and light-curve calculations were performed in spherical symmetry with the radiation-hydrodynamics code CRAB. Our results confirm the basic findings of the previous work: 3D neutrino-driven explosions with SN 1987A-like energies synthesize an amount of Ni-56 that is consistent with the radioactive tail of the light curve. Moreover, the models mix hydrogen inward to minimum velocities below 400 km/s as required by spectral observations. Hydrodynamic simulations with the new progenitor models, which possess smaller radii than the older ones, show much better agreement between calculated and observed light curves in the initial luminosity peak and during the first 20 days. A set of explosions with similar energies demonstrated that a high growth factor of Rayleigh-Taylor instabilities at the (C+O)/He composition interface combined with a weak interaction of fast Rayleigh-Taylor plumes, where the reverse shock occurs below the He/H interface, provides a sufficient condition for efficient outward mixing of Ni-56 into the hydrogen envelope. This condition is realized to the required extent only in one of the older stellar models, which yielded a maximum velocity of around 3000 km/s for the bulk of ejected Ni-56, but failed to reproduce the helium-core mass of 6 Msun inferred from the absolute luminosity of the presupernova star. We conclude that none of the single-star progenitor models proposed for SN 1987A to date satisfies all constraints set by observations. (Abridged), Comment: 16 pages, 9 figures, 3 tables; accepted for publication in Astron. Astrophys

A large fraction of core-collapse supernovae (CCSNe), 30-50%, are expected to originate from the low-mass end of progenitors with $M_{\rm ZAMS}~= 8-12~M_\odot$. However, degeneracy effects make stellar evolution modelling of such stars challenging, and few predictions for their supernova light curves and spectra have been presented. Here we calculate synthetic nebular spectra of a 9 $M_\odot$ Fe CCSN model exploded with the neutrino mechanism. The model predicts emission lines with FWHM$\sim$1000 km/s, including signatures from each deep layer in the metal core. We compare this model to observations of the three subluminous IIP SNe with published nebular spectra; SN 1997D, SN 2005cs, and SN 2008bk. The prediction of both line profiles and luminosities are in good agreement with SN 1997D and SN 2008bk. The close fit of a model with no tuning parameters provides strong evidence for an association of these objects with low-mass Fe CCSNe. For SN 2005cs, the interpretation is less clear, as the observational coverage ended before key diagnostic lines from the core had emerged. We perform a parameterised study of the amount of explosively made stable nickel, and find that none of these three SNe show the high $^{58}$Ni/$^{56}$Ni ratio predicted by current models of electron capture SNe (ECSNe) and ECSN-like explosions. Combined with clear detection of lines from O and He shell material, these SNe rather originate from Fe core progenitors. We argue that the outcome of self-consistent explosion simulations of low-mass stars, which gives fits to many key observables, strongly suggests that the class of subluminous Type IIP SNe is the observational counterpart of the lowest mass CCSNe., Comment: Resubmitted to MNRAS after referee comments

Nucleosynthesis, light curves, explosion energies, and remnant masses are calculated for a grid of supernovae resulting from massive stars with solar metallicity and masses from 9.0 to 120 solar masses. The full evolution is followed using an adaptive reaction network of up to 2000 nuclei. A novel aspect of the survey is the use of a one-dimensional neutrino transport model for the explosion. This explosion model has been calibrated to give the observed energy for SN 1987A, using several standard progenitors, and for the Crab supernova using a 9.6 solar mass progenitor. As a result of using a calibrated central engine, the final kinetic energy of the supernova is variable and sensitive to the structure of the presupernova star. Many progenitors with extended core structures do not explode, but become black holes, and the masses of exploding stars do not form a simply connected set. The resulting nucleosynthesis agrees reasonably well with the sun provided that a reasonable contribution from Type Ia supernovae is also allowed, but with a deficiency of light s-process isotopes. The resulting neutron star IMF has a mean gravitational mass near 1.4 solar masses. The average black hole mass is about 9 solar masses if only the helium core implodes, and 14 solar masses if the entire presupernova star collapses. Only ~10% of supernovae come from stars over 20 solar masses and some of these are Type Ib or Ic. Some useful systematics of Type IIp light curves are explored., Comment: 49 pages, 39 figures and 9 tables. Accepted to The Astrophysical Journal 26/1/2016

Thus far, judging the fate of a massive star (either a neutron star (NS) or a black hole) solely by its structure prior to core collapse has been ambiguous. Our work and previous attempts find a non-monotonic variation of successful and failed supernovae with zero-age main-sequence mass, for which no single structural parameter can serve as a good predictive measure. However, we identify two parameters computed from the pre-collapse structure of the progenitor, which in combination allow for a clear separation of exploding and non-exploding cases with only few exceptions (~1-2.5%) in our set of 621 investigated stellar models. One parameter is M4, defining the normalized enclosed mass for a dimensionless entropy per nucleon of s=4, and the other is mu4 = d(m/M_sun)/d(r/1000 km) at s=4, being the normalized mass-derivative at this location. The two parameters mu4 and M4*mu4 can be directly linked to the mass-infall rate, Mdot, of the collapsing star and the electron-type neutrino luminosity of the accreting proto-NS, L_nue ~ M_ns*Mdot, which play a crucial role in the "critical luminosity" concept for the theoretical description of neutrino-driven explosions as runaway phenomenon of the stalled accretion shock. All models were evolved employing the approach of Ugliano et al. for simulating neutrino-driven explosions in spherical symmetry. The neutrino emission of the accretion layer is approximated by a gray transport solver, while the uncertain neutrino emission of the 1.1 M_sun proto-NS core is parametrized by an analytic model. The free parameters connected to the core-boundary prescription are calibrated to reproduce the observables of Supernova 1987A for five different progenitor models., Comment: 23 pages, 12 figures; accepted by ApJ; revised version considerably enlarged (Fig. 7 and Sect.3.6 added)