-
Supernova 1987A: 3D Mixing and light curves for explosion models based on binary-merger progenitors
Authors:
V. P. Utrobin,
A. Wongwathanarat,
H. -Th. Janka,
E. Mueller,
T. Ertl,
A. Menon,
A. Heger
Abstract:
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 explosion…
▽ More
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.
△ Less
Submitted 21 April, 2021; v1 submitted 18 February, 2021;
originally announced February 2021.
-
Stellar Collapse Diversity and the Diffuse Supernova Neutrino Background
Authors:
Daniel Kresse,
Thomas Ertl,
Hans-Thomas Janka
Abstract:
The diffuse cosmic supernova neutrino background (DSNB) is observational target of the gadolinium-loaded Super-Kamiokande (SK) detector and the forthcoming JUNO and Hyper-Kamiokande detectors. Current predictions are hampered by our still incomplete understanding of the supernova (SN) explosion mechanism and of the neutron star (NS) equation of state and maximum mass. In our comprehensive study we…
▽ More
The diffuse cosmic supernova neutrino background (DSNB) is observational target of the gadolinium-loaded Super-Kamiokande (SK) detector and the forthcoming JUNO and Hyper-Kamiokande detectors. Current predictions are hampered by our still incomplete understanding of the supernova (SN) explosion mechanism and of the neutron star (NS) equation of state and maximum mass. In our comprehensive study we revisit this problem on grounds of the landscapes of successful and failed SN explosions obtained by Sukhbold et al. and Ertl et al. with parametrized one-dimensional neutrino engines for large sets of single-star and helium-star progenitors, with the latter serving as proxy of binary evolution effects. Besides considering engines of different strengths, leading to different fractions of failed SNe with black-hole (BH) formation, we also vary the NS mass limit, the spectral shape of the neutrino emission, and include contributions from poorly understood alternative NS-formation channels such as accretion-induced or merger-induced collapse events. Since the neutrino signals of our large model sets are approximate, we calibrate the associated degrees of freedom by using state-of-the-art simulations of proto-neutron star cooling. Our predictions are higher than other recent ones because of a large fraction of failed SNe with long delay to BH formation. Our best-guess model predicts a DSNB electron-antineutrino-flux of 28.8^{+24.6}_{-10.9} cm^{-2}s^{-1} with 6.0^{+5.1}_{-2.1} cm^{-2}s^{-1} in the favorable measurement interval of [10,30] MeV, and 1.3^{+1.1}_{-0.4} cm^{-2}s^{-1} with electron-antineutrino energies > 17.3 MeV, which is roughly a factor of two below the current SK limit. The uncertainty range is dominated by the still insufficiently constrained cosmic rate of stellar core-collapse events.
△ Less
Submitted 18 December, 2020; v1 submitted 9 October, 2020;
originally announced October 2020.
-
How much H and He is "hidden" in SNe Ib/c? -- II. Intermediate-mass objects: a 22 M$_{\odot}$ progenitor case study
Authors:
Jacob Teffs,
Thomas Ertl,
Paolo Mazzali,
Stephan Hachinger,
H. -Thomas Janka
Abstract:
Stripped envelope supernovae are a sub-class of core collapse supernovae showing several stages of H/He shell stripping that determines the class: H-free/He-poor SNe are classified as Type Ic, H-poor/He-rich are Type Ib, and H/He-rich are Type IIb. Stripping H/He with only stellar wind requires significantly higher mass loss rates than observed while binary-involved mass transfer may usually not s…
▽ More
Stripped envelope supernovae are a sub-class of core collapse supernovae showing several stages of H/He shell stripping that determines the class: H-free/He-poor SNe are classified as Type Ic, H-poor/He-rich are Type Ib, and H/He-rich are Type IIb. Stripping H/He with only stellar wind requires significantly higher mass loss rates than observed while binary-involved mass transfer may usually not strip enough to produce H/He free SNe. Type Ib/c SNe are sometimes found to include weak H/He transient lines as a product of a trace amount of H/He left over from stripping processes. The extent and mass of the H/He required to produce these lines is not well known. In this work, a 22 M$_{\odot}$ progenitor model is stripped of the H/He shells in five steps prior to collapse and then exploded at four explosion energies. Requiring both optical and NIR He I lines for helium identification does not allow much He mass to be hidden in SE--SNE. Increasing the mass of He above the CO core delays the visibility of O I 7774 in early spectra. Our SN Ib-like models are capable of reproducing the spectral evolution of a set of observed SNe with reasonable estimated $E_\mathrm{k}$ accuracy. Our SN\,IIb-like models can partially reproduce low energy observed SN IIb, but we find no observed comparison for the SN IIb-like models with high $E_\mathrm{k}$.
△ Less
Submitted 18 August, 2020;
originally announced August 2020.
-
Three-dimensional Models of Core-collapse Supernovae From Low-mass Progenitors With Implications for Crab
Authors:
G. Stockinger,
H. -Th. Janka,
D. Kresse,
T. Melson,
T. Ertl,
M. Gabler,
A. Gessner,
A. Wongwathanarat,
A. Tolstov,
S. -C. Leung,
K. Nomoto,
A. Heger
Abstract:
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…
▽ More
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.
△ Less
Submitted 10 June, 2020; v1 submitted 5 May, 2020;
originally announced May 2020.
-
Type Ic Supernova of a 22 $\text{M}_{\odot}$ Progenitor
Authors:
Jacob Teffs,
Thomas Ertl,
Paolo Mazzali,
Stephan Hachinger,
Thomas Janka
Abstract:
Type Ic supernovae (SNe Ic) are a sub-class of core-collapse supernovae that exhibit no helium or hydrogen lines in their spectra. Their progenitors are thought to be bare carbon-oxygen cores formed during the evolution of massive stars that are stripped of their hydrogen and helium envelopes sometime before collapse. SNe Ic present a range of luminosities and spectral properties, from luminous GR…
▽ More
Type Ic supernovae (SNe Ic) are a sub-class of core-collapse supernovae that exhibit no helium or hydrogen lines in their spectra. Their progenitors are thought to be bare carbon-oxygen cores formed during the evolution of massive stars that are stripped of their hydrogen and helium envelopes sometime before collapse. SNe Ic present a range of luminosities and spectral properties, from luminous GRB-SNe with broad-lined spectra to less luminous events with narrow-line spectra. Modelling SNe\,Ic reveals a wide range of both kinetic energies, ejecta masses, and $^{56}\text{Ni}$ masses. To explore this diversity and how it comes about, light curves and spectra are computed from the ejecta following the explosion of an initially 22 $\text{M}_{\odot}$ progenitor that was artificially stripped of its hydrogen and helium shells, producing a bare CO core of $\sim$ 5 $\text{M}_{\odot}$, resulting in an ejected mass of $\sim$ 4 $\text{M}_{\odot}$, which is an average value for SNe Ic. Four different explosion energies are used that cover a range of observed SNe. Finally, $^{56}\text{Ni}$ and other elements are artificially mixed in the ejecta using two approximations to determine how element distribution affects light curves and spectra. The combination of different explosion energy and degree of mixing produces spectra that roughly replicate the distribution of near-peak spectroscopic features of SNe Ic. High explosion energies combined with extensive mixing can produce red, broad-lined spectra, while minimal mixing and a lower explosion energy produce bluer, narrow-lined spectra.
△ Less
Submitted 29 January, 2020; v1 submitted 20 January, 2020;
originally announced January 2020.
-
The Explosion of Helium Stars Evolved With Mass Loss
Authors:
Thomas Ertl,
Stan E. Woosley,
Tuguldur Sukhbold,
H. -Thomas Janka
Abstract:
Light curves, explosion energies, and remnant masses are calculated for a grid of supernovae resulting from massive helium stars that have been evolved including mass loss. These presupernova stars should approximate the results of binary evolution for stars in interacting systems that lose their envelopes close to the time of helium core ignition. Initial helium star masses are in the range 2.5 t…
▽ More
Light curves, explosion energies, and remnant masses are calculated for a grid of supernovae resulting from massive helium stars that have been evolved including mass loss. These presupernova stars should approximate the results of binary evolution for stars in interacting systems that lose their envelopes close to the time of helium core ignition. Initial helium star masses are in the range 2.5 to 40\,\Msun, which correspond to main sequence masses of about 13 to 90\,\Msun. Common Type Ib and Ic supernovae result from stars whose final masses are approximately 2.5 to 5.6\,\Msun. For heavier stars, a large fraction of collapses lead to black holes, though there is an island of explodability for presupernova masses near 10\,\Msun. The median neutron star mass in binaries is 1.35--1.38\,\Msun \ and the median black hole mass is between 9 and 11\,\Msun. Even though black holes less massive than 5 \Msun\ are rare, they are predicted down to the maximum neutron star mass. There is no empty ``gap'', only a less populated mass range. For standard assumptions regarding the explosions and nucleosynthesis, the models predict light curves that are fainter than the brighter common Type Ib and Ic supernovae. Even with a very liberal, but physically plausible increase in $^{56}$Ni production, the highest energy models are fainter, at peak, than 10$^{42.6}$\,erg\,s$^{-1}$, and very few approach that limit. The median peak luminosity ranges from 10$^{42.0}$ to 10$^{42.3}$\,erg\,s$^{-1}$. Possible alternatives to the standard neutrino-powered and radioactive-illuminated models are explored. Magnetars are a promising alternative. Several other unusual varieties of Type I supernovae at both high and low mass are explored.
△ Less
Submitted 17 December, 2019; v1 submitted 3 October, 2019;
originally announced October 2019.
-
Three-dimensional mixing and light curves: constraints on the progenitor of supernova 1987A
Authors:
Victor Utrobin,
Annop Wongwathanarat,
H. -Thomas Janka,
Ewald Mueller,
T. Ertl,
Stan Woosley
Abstract:
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 pe…
▽ More
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)
△ Less
Submitted 19 March, 2019; v1 submitted 24 December, 2018;
originally announced December 2018.
-
Emission line models for the lowest-mass core collapse supernovae. I: Case study of a 9 $M_\odot$ one-dimensional neutrino-driven explosion
Authors:
A. Jerkstrand,
T. Ertl,
H. -T. Janka,
E. Müller,
T. Sukhbold,
S. E. Woosley
Abstract:
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 C…
▽ More
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.
△ Less
Submitted 12 October, 2017;
originally announced October 2017.
-
Erratum: Progenitor-explosion connection and remnant birth masses for neutrino-driven supernovae of iron-core progenitors (2012, ApJ, 757, 69)
Authors:
Thomas Ertl,
Marcella Ugliano,
H. -Thomas Janka,
Andreas Marek,
Almudena Arcones
Abstract:
An erroneous interpretation of the hydrodynamical results led to an incorrect determination of the fallback masses in Ugliano et al. (2012), which also (on a smaller level) affects the neutron star masses provided in that paper. This problem was already addressed and corrected in the follow-up works by Ertl et al. (2015) and Sukhbold et al. (2015). Therefore, the reader is advised to use the new d…
▽ More
An erroneous interpretation of the hydrodynamical results led to an incorrect determination of the fallback masses in Ugliano et al. (2012), which also (on a smaller level) affects the neutron star masses provided in that paper. This problem was already addressed and corrected in the follow-up works by Ertl et al. (2015) and Sukhbold et al. (2015). Therefore, the reader is advised to use the new data of the latter two publications. In the remaining text of this Erratum we present the differences of the old and new fallback results in detail and explain the origin of the mistake in the original analysis by Ugliano et al. (2012).
△ Less
Submitted 17 February, 2016;
originally announced February 2016.
-
Core-Collapse Supernovae from 9 to 120 Solar Masses Based on Neutrino-powered Explosions
Authors:
Tuguldur Sukhbold,
T. Ertl,
S. E. Woosley,
Justin M. Brown,
H. -T. Janka
Abstract:
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…
▽ More
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.
△ Less
Submitted 13 February, 2016; v1 submitted 15 October, 2015;
originally announced October 2015.
-
A two-parameter criterion for classifying the explodability of massive stars by the neutrino-driven mechanism
Authors:
T. Ertl,
H. -Th. Janka,
S. E. Woosley,
T. Sukhbold,
M. Ugliano
Abstract:
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 param…
▽ More
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.
△ Less
Submitted 8 January, 2016; v1 submitted 25 March, 2015;
originally announced March 2015.