Showing 1–50 of 187 results for author: Janka, H -

Shockingly Bright Warm Carbon Monoxide Molecular Features in the Supernova Remnant Cassiopeia A Revealed by JWST

Authors: J. Rho, S. -H. Park, R. Arendt, M. Matsuura, D. Milisavljevic, T. Temim, I. De Looze, W. P. Blair, A. Rest, O. Fox, A. P. Ravi, B. -C. Koo, M. Barlow, A. Burrows, R. Chevalier, G. Clayton, R. Fesen, C. Fransson, C. Fryer, H. L. Gomez, H. -T. Janka, F. Kirchschlarger, J. M. Laming, S. Orlando, D. Patnaude , et al. (14 additional authors not shown)

Abstract: We present JWST NIRCam (F356W and F444W filters) and MIRI (F770W) images and NIRSpec- IFU spectroscopy of the young supernova remnant Cassiopeia A (Cas A). We obtained the data as part of a JWST survey of Cas A. The NIRCam and MIRI images map the spatial distributions of synchrotron radiation, Ar-rich ejecta, and CO on both large and small scales, revealing remarkably complex structures. The CO em… ▽ More We present JWST NIRCam (F356W and F444W filters) and MIRI (F770W) images and NIRSpec- IFU spectroscopy of the young supernova remnant Cassiopeia A (Cas A). We obtained the data as part of a JWST survey of Cas A. The NIRCam and MIRI images map the spatial distributions of synchrotron radiation, Ar-rich ejecta, and CO on both large and small scales, revealing remarkably complex structures. The CO emission is stronger at the outer layers than the Ar ejecta, which indicates the reformation of CO molecules behind the reverse shock. NIRSpec-IFU spectra (3 - 5.5 microns) were obtained toward two representative knots in the NE and S fields. Both regions are dominated by the bright fundamental rovibrational band of CO in the two R and P branches, with strong [Ar VI] and relatively weaker, variable strength ejecta lines of [Si IX], [Ca IV], [Ca V] and [Mg IV]. The NIRSpec-IFU data resolve individual ejecta knots and filaments spatially and in velocity space. The fundamental CO band in the JWST spectra reveals unique shapes of CO, showing a few tens of sinusoidal patterns of rovibrational lines with pseudo-continuum underneath, which is attributed to the high-velocity widths of CO lines. The CO also shows high J lines at different vibrational transitions. Our results with LTE modeling of CO emission indicate a temperature of 1080 K and provide unique insight into the correlations between dust, molecules, and highly ionized ejecta in supernovae, and have strong ramifications for modeling dust formation that is led by CO cooling in the early Universe. △ Less

Submitted 5 June, 2024; originally announced June 2024.

Comments: accepted for the ApJ letter (17 pages and 10 figures)

Simple fits for the neutrino luminosities from protoneutron star cooling

Authors: Giuseppe Lucente, Malte Heinlein, H. -Thomas Janka, Alessandro Mirizzi

Abstract: We propose a simple fit function, $L_{ν_i}(t) = C\, t^{-α}\, e^{-(t/τ)^{n}}$, to parametrize the luminosities of neutrinos and antineutrinos of all flavors during the protoneutron star (PNS) cooling phase at post-bounce times $t \gtrsim 1$ s. This fit is based on results from a set of neutrino-hydrodynamics simulations of core-collapse supernovae in spherical symmetry. The simulations were perform… ▽ More We propose a simple fit function, $L_{ν_i}(t) = C\, t^{-α}\, e^{-(t/τ)^{n}}$, to parametrize the luminosities of neutrinos and antineutrinos of all flavors during the protoneutron star (PNS) cooling phase at post-bounce times $t \gtrsim 1$ s. This fit is based on results from a set of neutrino-hydrodynamics simulations of core-collapse supernovae in spherical symmetry. The simulations were performed with an energy-dependent transport for six neutrino species and took into account the effects of convection and muons in the dense and hot PNS interior. We provide values of the fit parameters $C$, $α$, $τ$, and $n$ for different neutron star masses and equations of state as well as correlations between these fit parameters. Our functional description is useful for analytic supernova modeling, for characterizing the neutrino light curves in large underground neutrino detectors, and as a tool to extract information from measured signals on the mass and equation of state of the PNS and on secondary signal components on top of the PNS's neutrino emission. △ Less

Submitted 8 May, 2024; v1 submitted 1 May, 2024; originally announced May 2024.

Comments: 30 pages, 16 figures, 11 tables. Minor changes, additional references included

Interplay Between Neutrino Kicks and Hydrodynamic Kicks of Neutron Stars and Black Holes

Authors: H. -Thomas Janka, Daniel Kresse

Abstract: Neutron stars (NSs) are observed with high space velocities and elliptical orbits in binaries. The magnitude of these effects points to natal kicks that originate from asymmetries during the supernova (SN) explosions. Using a growing set of long-time 3D SN simulations with the Prometheus-Vertex code, we explore the interplay of NS kicks that are induced by asymmetric neutrino emission and by asymm… ▽ More Neutron stars (NSs) are observed with high space velocities and elliptical orbits in binaries. The magnitude of these effects points to natal kicks that originate from asymmetries during the supernova (SN) explosions. Using a growing set of long-time 3D SN simulations with the Prometheus-Vertex code, we explore the interplay of NS kicks that are induced by asymmetric neutrino emission and by asymmetric mass ejection. Anisotropic neutrino emission can arise from a large-amplitude dipolar convection asymmetry inside the proto-NS (PNS) termed LESA (Lepton-number Emission Self-sustained Asymmetry) and from aspherical accretion downflows around the PNS, which can lead to anisotropic neutrino emission (absorption) with a neutrino-induced NS kick roughly opposite to (aligned with) the kick by asymmetric mass ejection. In massive progenitors hydrodynamic kicks can reach up to more than 1300 km/s, whereas our calculated neutrino kicks reach (55-140) km/s (estimated upper bounds of (170-265) km/s) and only about (10-50) km/s, if LESA is the main cause of asymmetric neutrino emission. Therefore hydrodynamic NS kicks dominate in explosions of high-mass progenitors, whereas LESA-induced neutrino kicks dominate for NSs born in low-energy SNe of the lowest-mass progenitors, when these explode nearly spherically. Our models suggest that the Crab pulsar with its velocity of about 160 km/s, if born in the low-energy explosion of a low-mass, single-star progenitor, should have received a hydrodynamic kick in a considerably asymmetric explosion. Black holes, if formed by the collapse of short-lived PNSs and solely kicked by anisotropic neutrino emission, obtain velocities of only some km/s. △ Less

Submitted 3 June, 2024; v1 submitted 24 January, 2024; originally announced January 2024.

Comments: 68 pages, 29 figures, 3 tables; submitted to Astrophysics and Space Science (Ap&SS), 2022 Astronomy Prize Awardees Collection; significantly extended version because of relevant comments by referee and readers

Dynamics and Equation of State Dependencies of Relevance for Nucleosynthesis in Supernovae and Neutron Star Mergers

Authors: H. -Thomas Janka, Andreas Bauswein

Abstract: Neutron stars (NSs) and black holes (BHs) are born when the final collapse of the stellar core terminates the lives of stars more massive than about 9 Msun. This can trigger the powerful ejection of a large fraction of the star's material in a core-collapse supernova (CCSN), whose extreme luminosity is energized by the decay of radioactive isotopes such as 56Ni and 56Co. When evolving in close bin… ▽ More Neutron stars (NSs) and black holes (BHs) are born when the final collapse of the stellar core terminates the lives of stars more massive than about 9 Msun. This can trigger the powerful ejection of a large fraction of the star's material in a core-collapse supernova (CCSN), whose extreme luminosity is energized by the decay of radioactive isotopes such as 56Ni and 56Co. When evolving in close binary systems, the compact relics of such infernal catastrophes spiral towards each other on orbits gradually decaying by gravitational-wave emission. Ultimately, the violent collision of the two components forms a more massive, rapidly spinning remnant, again accompanied by the ejection of considerable amounts of matter. These merger events can be observed by high-energy bursts of gamma rays with afterglows and electromagnetic transients called kilonovae, which radiate the energy released in radioactive decays of freshly assembled rapid neutron-capture elements. By means of their mass ejection and the nuclear and neutrino reactions taking place in the ejecta, both CCSNe and compact object mergers (COMs) are prominent sites of heavy-element nucleosynthesis and play a central role in the cosmic cycle of matter and the chemical enrichment history of galaxies. The nuclear equation of state (EoS) of NS matter, from neutron-rich to proton-dominated conditions and with temperatures ranging from about zero to ~100 MeV, is a crucial ingredient in these astrophysical phenomena. It determines their dynamical processes, their remnant properties even at the level of deciding between NS or BH, and the properties of the associated emission of neutrinos, whose interactions govern the thermodynamic conditions and the neutron-to-proton ratio for nucleosynthesis reactions in the innermost ejecta. This chapter discusses corresponding EoS dependent effects of relevance in CCSNe as well as COMs. (slightly abridged) △ Less

Submitted 20 February, 2023; v1 submitted 14 December, 2022; originally announced December 2022.

Comments: 96 pages, 23 figures; figure citations added; accepted for the Handbook of Nuclear Physics, eds. Isao Tanihata, Hiroshi Toki, Toshitaka Kajino

Parameterisations of thermal bomb explosions for core-collapse supernovae and 56Ni production

Authors: Liliya Imasheva, H. -Thomas Janka, Achim Weiss

Abstract: Thermal bombs are a widely used method to artificially trigger explosions of core-collapse supernovae (CCSNe) to determine their nucleosynthesis or ejecta and remnant properties. Recently, their use in spherically symmetric (1D) hydrodynamic simulations led to the result that {56,57}Ni and 44Ti are massively underproduced compared to observational estimates for Supernova 1987A, if the explosions a… ▽ More Thermal bombs are a widely used method to artificially trigger explosions of core-collapse supernovae (CCSNe) to determine their nucleosynthesis or ejecta and remnant properties. Recently, their use in spherically symmetric (1D) hydrodynamic simulations led to the result that {56,57}Ni and 44Ti are massively underproduced compared to observational estimates for Supernova 1987A, if the explosions are slow, i.e., if the explosion mechanism of CCSNe releases the explosion energy on long timescales. It was concluded that rapid explosions are required to match observed abundances, i.e., the explosion mechanism must provide the CCSN energy nearly instantaneously on timescales of some ten to order 100 ms. This result, if valid, would disfavor the neutrino-heating mechanism, which releases the CCSN energy on timescales of seconds. Here, we demonstrate by 1D hydrodynamic simulations and nucleosynthetic post-processing that these conclusions are a consequence of disregarding the initial collapse of the stellar core in the thermal-bomb modelling before the bomb releases the explosion energy. We demonstrate that the anti-correlation of 56Ni yield and energy-injection timescale vanishes when the initial collapse is included and that it can even be reversed, i.e., more 56Ni is made by slower explosions, when the collapse proceeds to small radii similar to those where neutrino heating takes place in CCSNe. We also show that the 56Ni production in thermal-bomb explosions is sensitive to the chosen mass cut and that a fixed mass layer or fixed volume for the energy deposition cause only secondary differences. Moreover, we propose a most appropriate setup for thermal bombs. △ Less

Submitted 7 November, 2022; v1 submitted 22 September, 2022; originally announced September 2022.

Comments: 24 pages, 12 figures; version accepted by MNRAS after minor revision

Impact of systematic nuclear uncertainties on composition and decay heat of dynamical and disk ejecta in compact binary mergers

Authors: I. Kullmann, S. Goriely, O. Just, A. Bauswein, H. -T. Janka

Abstract: Theoretically predicted yields of elements created by the rapid neutron capture (r-) process carry potentially large uncertainties associated with incomplete knowledge of nuclear properties and approximative hydrodynamical modelling of the matter ejection processes. We present an in-depth study of the nuclear uncertainties by varying theoretical nuclear input models that describe the experimentall… ▽ More Theoretically predicted yields of elements created by the rapid neutron capture (r-) process carry potentially large uncertainties associated with incomplete knowledge of nuclear properties and approximative hydrodynamical modelling of the matter ejection processes. We present an in-depth study of the nuclear uncertainties by varying theoretical nuclear input models that describe the experimentally unknown neutron-rich nuclei. This includes two frameworks for calculating the radiative neutron capture rates and 14 different models for nuclear masses, $β$-decay rates and fission properties. Our r-process nuclear network calculations are based on detailed hydrodynamical simulations of dynamically ejected material from NS-NS or NS-BH binary mergers plus the secular ejecta from BH-torus systems. The impact of nuclear uncertainties on the r-process abundance distribution and the early radioactive heating rate is found to be modest (within a factor of $\sim20$ for individual $A>90$ abundances and a factor of 2 for the heating rate). However, the impact on the late-time heating rate is more significant and depends strongly on the contribution from fission. We witness significantly larger sensitivity to the nuclear physics input if only a single trajectory is used compared to considering ensembles of $\sim$200-300 trajectories, and the quantitative effects of the nuclear uncertainties strongly depend on the adopted conditions for the individual trajectory. We use the predicted Th/U ratio to estimate the cosmochronometric age of six metal-poor stars and find the impact of the nuclear uncertainties to be up to 2 Gyr. △ Less

Submitted 11 May, 2023; v1 submitted 15 July, 2022; originally announced July 2022.

Comments: 28 pages, 24 figures, accepted to MNRAS

A new scenario for magnetar formation: Tayler-Spruit dynamo in a proto-neutron star spun up by fallback

Authors: P. Barrère, J. Guilet, A. Reboul-Salze, R. Raynaud, H. -T. Janka

Abstract: Magnetars are isolated young neutron stars characterized by the most intense magnetic fields known in the universe. The origin of their magnetic field is still a challenging question. In situ magnetic field amplification by dynamo action is a promising process to generate ultra-strong magnetic fields in fast-rotating progenitors. However, it is unclear whether the fraction of progenitors harboring… ▽ More Magnetars are isolated young neutron stars characterized by the most intense magnetic fields known in the universe. The origin of their magnetic field is still a challenging question. In situ magnetic field amplification by dynamo action is a promising process to generate ultra-strong magnetic fields in fast-rotating progenitors. However, it is unclear whether the fraction of progenitors harboring fast core rotation is sufficient to explain the entire magnetar population. To address this point, we propose a new scenario for magnetar formation, in which a slow-rotating proto-neutron star is spun up by the supernova fallback. We argue that this can trigger the development of the Tayler-Spruit dynamo while other dynamo processes are disfavored. Using previous works done on this dynamo and simulations to characterize the fallback, we derive equations modelling the coupled evolution of the proto-neutron star rotation and magnetic field. Their time integration for different fallback masses is successfully compared with analytical estimates of the amplification timescales and saturation value of the magnetic field. We find that the magnetic field is amplified within $20$ to $40$s after the core bounce, and that the radial magnetic field saturates at intensities $10^{14}-10^{15}$G, therefore spanning the full range of magnetar's dipolar magnetic fields. We also compare predictions of two proposed saturation mechanisms showing that magnetar-like magnetic fields can be generated for neutron star spun up to rotation periods $\lesssim8$ms and $\lesssim28$ms, corresponding to fallback masses $\gtrsim4\times10^{-2}{\rm M}_{\odot}$ and $\gtrsim10^{-2}{\rm M}_{\odot}$. Thus, our results suggest that magnetars can be formed from slow-rotating progenitors for fallback masses compatible with recent supernova simulations and leading to plausible initial rotation periods of the proto-neutron star. △ Less

Submitted 19 October, 2022; v1 submitted 2 June, 2022; originally announced June 2022.

Comments: 14 pages, 5 figures, accepted for publication in A&A

Journal ref: A&A 668, A79 (2022)

Low-Energy Physics in Neutrino LArTPCs

Authors: D. Caratelli, W. Foreman, A. Friedland, S. Gardiner, I. Gil-Botella, G. Karagiorgi, M. Kirby, G. Lehmann Miotto, B. R. Littlejohn, M. Mooney, J. Reichenbacher, A. Sousa, K. Scholberg, J. Yu, T. Yang, S. Andringa, J. Asaadi, T. J. C. Bezerra, F. Capozzi, F. Cavanna, E. Church, A. Himmel, T. Junk, J. Klein, I. Lepetic , et al. (264 additional authors not shown)

Abstract: In this white paper, we outline some of the scientific opportunities and challenges related to detection and reconstruction of low-energy (less than 100 MeV) signatures in liquid argon time-projection chamber (LArTPC) detectors. Key takeaways are summarized as follows. 1) LArTPCs have unique sensitivity to a range of physics and astrophysics signatures via detection of event features at and below… ▽ More In this white paper, we outline some of the scientific opportunities and challenges related to detection and reconstruction of low-energy (less than 100 MeV) signatures in liquid argon time-projection chamber (LArTPC) detectors. Key takeaways are summarized as follows. 1) LArTPCs have unique sensitivity to a range of physics and astrophysics signatures via detection of event features at and below the few tens of MeV range. 2) Low-energy signatures are an integral part of GeV-scale accelerator neutrino interaction final states, and their reconstruction can enhance the oscillation physics sensitivities of LArTPC experiments. 3) BSM signals from accelerator and natural sources also generate diverse signatures in the low-energy range, and reconstruction of these signatures can increase the breadth of BSM scenarios accessible in LArTPC-based searches. 4) Neutrino interaction cross sections and other nuclear physics processes in argon relevant to sub-hundred-MeV LArTPC signatures are poorly understood. Improved theory and experimental measurements are needed. Pion decay-at-rest sources and charged particle and neutron test beams are ideal facilities for experimentally improving this understanding. 5) There are specific calibration needs in the low-energy range, as well as specific needs for control and understanding of radiological and cosmogenic backgrounds. 6) Novel ideas for future LArTPC technology that enhance low-energy capabilities should be explored. These include novel charge enhancement and readout systems, enhanced photon detection, low radioactivity argon, and xenon doping. 7) Low-energy signatures, whether steady-state or part of a supernova burst or larger GeV-scale event topology, have specific triggering, DAQ and reconstruction requirements that must be addressed outside the scope of conventional GeV-scale data collection and analysis pathways. △ Less

Submitted 1 March, 2022; originally announced March 2022.

Comments: Contribution to Snowmass 2021

Evidence for past interaction with an asymmetric circumstellar shell in the young SNR Cassiopeia A

Authors: S. Orlando, A. Wongwathanarat, H. -T. Janka, M. Miceli, S. Nagataki, M. Ono, F. Bocchino, J. Vink, D. Milisavljevic, D. J. Patnaude, G. Peres

Abstract: Observations of the SNR Cassiopeia A (Cas A) show asymmetries in the reverse shock that cannot be explained by models describing a remnant expanding through a spherically symmetric wind of the progenitor star. We investigate whether a past interaction of Cas A with an asymmetric circumstellar shell can account for the observed asymmetries. We performed 3D MHD simulations that describe the remnant… ▽ More Observations of the SNR Cassiopeia A (Cas A) show asymmetries in the reverse shock that cannot be explained by models describing a remnant expanding through a spherically symmetric wind of the progenitor star. We investigate whether a past interaction of Cas A with an asymmetric circumstellar shell can account for the observed asymmetries. We performed 3D MHD simulations that describe the remnant evolution from the SN to its interaction with a circumstellar shell. The initial conditions are provided by a 3D neutrino-driven SN model whose morphology resembles Cas A. We explored the parameter space of the shell, searching for a set of parameters able to produce reverse shock asymmetries at the age of 350 years analogous to those observed in Cas A. The interaction of the remnant with the shell can produce asymmetries resembling those observed in the reverse shock if the shell was asymmetric with the densest portion in the nearside to the northwest (NW). The reverse shock shows the following asymmetries at the age of Cas A: i) it moves inward in the observer frame in the NW region, while it moves outward in other regions; ii) the geometric center of the reverse shock is offset to the NW from the geometric center of the forward shock; iii) the reverse shock in the NW region has enhanced nonthermal emission because, there, the ejecta enter the reverse shock with a higher velocity (between 4000 and 7000 km/s) than in other regions (below 2000 km/s). The asymmetries observed in the reverse shock of Cas A can be interpreted as signatures of the interaction of the remnant with an asymmetric circumstellar shell that occurred between 180 and 240 years after the SN event. We suggest that the shell was, most likely, the result of a massive eruption from the progenitor star that occurred between $10^4$ and $10^5$ years prior to core-collapse. We estimate a total mass of the shell of the order 2 Msun. △ Less

Submitted 12 August, 2022; v1 submitted 3 February, 2022; originally announced February 2022.

Comments: 20 pages, 11 figures, 2 Tables, 1 interactive graphics. Accepted for publication on A&A, high quality figures will be available on the published paper

Journal ref: A&A 666, A2 (2022)

Nebular phase properties of supernova Ibc from He-star explosions

Authors: L. Dessart, D. J. Hillier, T. Sukhbold, S. E. Woosley, H. -T. Janka

Abstract: Following our recent work on Type II supernovae (SNe), we present a set of 1D nonlocal thermodynamic equilibrium radiative transfer calculations for nebular-phase Type Ibc SNe starting from state-of-the-art explosion models with detailed nucleosynthesis. Our grid of progenitor models is derived from He stars that were subsequently evolved under the influence of wind mass loss. These He stars, whic… ▽ More Following our recent work on Type II supernovae (SNe), we present a set of 1D nonlocal thermodynamic equilibrium radiative transfer calculations for nebular-phase Type Ibc SNe starting from state-of-the-art explosion models with detailed nucleosynthesis. Our grid of progenitor models is derived from He stars that were subsequently evolved under the influence of wind mass loss. These He stars, which most likely form through binary mass exchange, synthesize less oxygen than their single-star counterparts with the same zero-age main sequence (ZAMS) mass. This reduction is greater in He-star models evolved with an enhanced mass loss rate. We obtain a wide range of spectral properties at 200d. In models from He stars with an initial mass >6Msun, the [OI] 6300, 6364 is of comparable or greater strength than [CaII] 7291,7323 -- the strength of [OI] 6300, 6364 increases with He-star initial mass. In contrast, models from lower mass He stars exhibit a weak [OI] 6300, 6364, strong [CaII] 7291,7323, but also strong NII lines and FeII emission below 5500A. The ejecta density, modulated by the ejecta mass, the explosion energy, and clumping, has a critical impact on the gas ionization, line cooling, and the spectral properties. FeII dominates the emission below 5500A and is stronger at earlier nebular epochs. It ebbs as the SN ages, while the fractional flux in [OI] 6300, 6364 and [CaII] 7291,7323 increases, with a similar rate, as the ejecta recombine. Although the results depend on the adopted wind mass loss rate and pre-SN mass, we find that He stars of 6-8Msun initially (ZAMS mass of 23-28Msun) match adequately the properties of standard SNe Ibc. Our results for less massive He stars are more perplexing, since the predicted spectra are not seen in nature. They may be missed by current surveys or associated with Type Ibn SNe in which interaction dominates over decay power. [Abridged] △ Less

Submitted 25 September, 2021; originally announced September 2021.

Comments: Accepted for publication in A&A

Journal ref: A&A 656, A61 (2021)

Dynamical ejecta of neutron star mergers with nucleonic weak processes I: Nucleosynthesis

Authors: I. Kullmann, S. Goriely, O. Just, R. Ardevol-Pulpillo, A. Bauswein, H. -T. Janka

Abstract: We present a coherent study of the impact of neutrino interactions on the r-process element nucleosynthesis and the heating rate produced by the radioactive elements synthesised in the dynamical ejecta of neutron star-neutron star (NS-NS) mergers. We have studied the material ejected from four NS-NS merger systems based on hydrodynamical simulations which handle neutrino effects in an elaborate wa… ▽ More We present a coherent study of the impact of neutrino interactions on the r-process element nucleosynthesis and the heating rate produced by the radioactive elements synthesised in the dynamical ejecta of neutron star-neutron star (NS-NS) mergers. We have studied the material ejected from four NS-NS merger systems based on hydrodynamical simulations which handle neutrino effects in an elaborate way by including neutrino equilibration with matter in optically thick regions and re-absorption in optically thin regions. We find that the neutron richness of the dynamical ejecta is significantly affected by the neutrinos emitted by the post-merger remnant, in particular when compared to a case neglecting all neutrino interactions. Our nucleosynthesis results show that a solar-like distribution of r- process elements with mass numbers $A \gtrsim 90$ is produced, including a significant enrichment in Sr and a reduced production of actinides compared to simulations without inclusion of the nucleonic weak processes. The composition of the dynamically ejected matter as well as the corresponding rate of radioactive decay heating are found to be rather independent of the system mass asymmetry and the adopted equation of state. This approximate degeneracy in abundance pattern and heating rates can be favourable for extracting the ejecta properties from kilonova observations, at least if the dynamical component dominates the overall ejecta. Part II of this work will study the light curve produced by the dynamical ejecta of our four NS merger models. △ Less

Submitted 20 November, 2021; v1 submitted 6 September, 2021; originally announced September 2021.

Comments: 18 pages, 21 figures, 2 Tables, accepted by MNRAS 2021-11-19

The explosion of 9$-$29$M_\odot$ stars as Type II supernovae : results from radiative-transfer modeling at one year after explosion

Authors: Luc Dessart, D. John Hillier, Tuguldur Sukhbold, Stan Woosley, H. -T. Janka

Abstract: 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]… ▽ More 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]$λλ$$6300,6364$ line flux generally grows with progenitor mass, and H$α$ exhibits an equally strong and opposite trend. The [CaII]$λλ$$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. △ Less

Submitted 27 May, 2021; originally announced May 2021.

Comments: Accepted for publication in A&A

Journal ref: A&A 652, A64 (2021)

Supernova Fallback as Origin of Neutron Star Spins and Spin-kick Alignment

Authors: H. -Thomas Janka, Annop Wongwathanarat, Michael Kramer

Abstract: Natal kicks and spins are characteristic properties of neutron stars (NSs) and black holes (BHs). Both offer valuable clues to dynamical processes during stellar core collapse and explosion. Moreover, they influence the evolution of stellar multiple systems and the gravitational-wave signals from their inspiral and merger. Observational evidence of possibly generic spin-kick alignment has been int… ▽ More Natal kicks and spins are characteristic properties of neutron stars (NSs) and black holes (BHs). Both offer valuable clues to dynamical processes during stellar core collapse and explosion. Moreover, they influence the evolution of stellar multiple systems and the gravitational-wave signals from their inspiral and merger. Observational evidence of possibly generic spin-kick alignment has been interpreted as indication that NS spins are either induced with the NS kicks or inherited from progenitor rotation, which thus might play a dynamically important role during stellar collapse. Current three-dimensional supernova simulations suggest that NS kicks are transferred in the first seconds of the explosion, mainly by anisotropic mass ejection and, on a secondary level, anisotropic neutrino emission. In contrast, the NS spins are only determined minutes to hours later by angular momentum associated with fallback of matter that does not become gravitationally unbound in the supernova. Here, we propose a novel scenario to explain spin-kick alignment as a consequence of tangential vortex flows in the fallback matter that is accreted mostly from the direction of the NS's motion. For this effect the initial NS kick is crucial, because it produces a growing offset of the NS away from the explosion center, thus promoting onesided accretion. In this new scenario conclusions based on traditional concepts are reversed. For example, pre-kick NS spins are not required, and rapid progenitor-core rotation can hamper spin-kick alignment. We also discuss implications for natal BH kicks and the possibility of tossing the BH's spin axis during its formation. △ Less

Submitted 5 December, 2021; v1 submitted 15 April, 2021; originally announced April 2021.

Comments: 38 pages, 9 figures, 1 table; includes new Sect.2 summarizing observations and new Fig. 8; considerably extended because of comments by referee and readers; accepted by ApJ

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.

Comments: 31 pages, 14 figures, 6 tables. Accepted for publication in ApJ

Fission fragment distributions and their impact on the r-process nucleosynthesis in neutron star mergers

Authors: J. -F. Lemaître, S. Goriely, A. Bauswein, H. -T. Janka

Abstract: Neutron star (NS) merger ejecta offer a viable site for the production of heavy r-process elements with nuclear mass numbers A >140. The crucial role of fission recycling is responsible for the robustness of this site against many astrophysical uncertainties. Here, we introduce new improvements to our scission-point model, called SPY, to derive the fission fragment distribution for all neutron-ric… ▽ More Neutron star (NS) merger ejecta offer a viable site for the production of heavy r-process elements with nuclear mass numbers A >140. The crucial role of fission recycling is responsible for the robustness of this site against many astrophysical uncertainties. Here, we introduce new improvements to our scission-point model, called SPY, to derive the fission fragment distribution for all neutron-rich fissioning nuclei of relevance in r-process calculations. These improvements include a phenomenological modification of the scission distance and a smoothing procedure of the distribution. Such corrections lead to a much better agreement with experimental fission yields. Those yields are also used to estimate the number of neutrons emitted by the excited fragments on the basis of different neutron evaporation models. Our new fission yields are extensively compared to those predicted by the so-called GEF model. The impact of fission on the r-process nucleosynthesis in binary neutron mergers is also reanalyzed. Two scenarios are considered, the first one with low initial electron fraction is subject to intense fission recycling, in contrast to the second one which includes weak interactions on nucleons. The various regions of the nuclear chart responsible for fission recycling during the neutron irradiation as well as after freeze-out are discussed. The contribution fission processes may have to the final abundance distribution is also studied in detail in the light of newly defined quantitative indicators describing the fission recycling, the fission seeds and the fission progenitors. In particular, those allow us to estimate the contribution of fission to the final abundance distribution stemming from specific heavy nuclei. Calculations obtained with SPY and GEF fission fragment distributions are compared for both r-process scenarios. △ Less

Submitted 6 February, 2021; originally announced February 2021.

Linking core-collapse supernova explosions to supernova remnants through 3D MHD modeling

Authors: S. Orlando, A. Wongwathanarat, H. -T. Janka, M. Miceli, M. Ono, S. Nagataki, F. Bocchino, G. Peres

Abstract: The structure and morphology of supernova remnants (SNRs) reflect the properties of the parent supernovae (SNe) and the characteristics of the inhomogeneous environments through which the remnants expand. Linking the morphology of SNRs to anisotropies developed in their parent SNe can be essential to obtain key information on many aspects of the explosion processes associated with SNe. Nowadays, o… ▽ More The structure and morphology of supernova remnants (SNRs) reflect the properties of the parent supernovae (SNe) and the characteristics of the inhomogeneous environments through which the remnants expand. Linking the morphology of SNRs to anisotropies developed in their parent SNe can be essential to obtain key information on many aspects of the explosion processes associated with SNe. Nowadays, our capability to study the SN-SNR connection has been largely improved thanks to multi-dimensional models describing the long-term evolution from the SN to the SNR as well as to observational data of growing quality and quantity across the electromagnetic spectrum which allow to constrain the models. Here we used the numerical resources obtained in the framework of the "Accordo Quadro INAF-CINECA (2017)" together with a CINECA ISCRA Award N.HP10BARP6Y to describe the full evolution of a SNR from the core-collapse to the full-fledged SNR at the age of 2000 years. Our simulations were compared with observations of SNR Cassiopeia A (Cas A) at the age of $\sim 350$~years. Thanks to these simulations we were able to link the physical, chemical and morphological properties of a SNR to the physical processes governing the complex phases of the SN explosion. △ Less

Submitted 25 December, 2020; originally announced December 2020.

Comments: 4 pages, 2 figures, report on the numerical resources obtained in the framework of the "Accordo Quadro INAF-CINECA (2017)". arXiv admin note: text overlap with arXiv:2009.01789

Fast neutrino flavor conversions in one-dimensional core-collapse supernova models with and without muon creation

Authors: Francesco Capozzi, Sajad Abbar, Robert Bollig, H. -Thomas Janka

Abstract: In very dense environments, neutrinos can undergo fast flavor conversions on scales as short as a few centimeters provided that the angular distribution of the neutrino lepton number crosses zero. This work presents the first attempt to establish whether the non-negligible abundance of muons and their interactions with neutrinos in the core of supernovae can affect the occurrence of such crossings… ▽ More In very dense environments, neutrinos can undergo fast flavor conversions on scales as short as a few centimeters provided that the angular distribution of the neutrino lepton number crosses zero. This work presents the first attempt to establish whether the non-negligible abundance of muons and their interactions with neutrinos in the core of supernovae can affect the occurrence of such crossings. For this purpose we employ state-of-the-art one-dimensional core-collapse supernova simulations, considering models that include muon-neutrino interactions as well as models without these reactions. Although a consistent treatment of muons in the equation of state and neutrino transport does not seem to modify significantly the conditions for the occurrence of fast modes, it allows for the existence of an interesting phenomenon, namely fast instabilities in the $μ-τ$ sector. We also show that crossings below the supernova shock are a relatively generic feature of the one-dimensional simulations under investigation, which contrasts with the previous reports in the literature. Our results highlight the importance of multi-dimensional simulations with muon creation, where our results must be tested in the future. △ Less

Submitted 11 March, 2021; v1 submitted 15 December, 2020; originally announced December 2020.

Comments: 10 pages, 5 figures in the main text. An appendix has been added with 3 figures showing a comparison between results obtained with different assumptions. The main results are unchanged, but their interpretation has been slightly modified and sometimes extended to improve clarity. Matches the version accepted on PRD

Report number: MPP-2020-226

Journal ref: Phys. Rev. D 103, 063013 (2021)

On the characteristics of fast neutrino flavor instabilities in three-dimensional core-collapse supernova models

Authors: Sajad Abbar, Francesco Capozzi, Robert Glas, H. -Thomas Janka, Irene Tamborra

Abstract: We assess the occurrence of fast neutrino flavor instabilities in two three-dimensional state-of-the-art core-collapse supernova simulations performed using a two-moment three-species neutrino transport scheme: one with an exploding 9$\mathrm{M_{\odot}}$ and one with a non-exploding 20$\mathrm{M_{\odot}}$ model. Apart from confirming the presence of fast instabilities occurring within the neutrino… ▽ More We assess the occurrence of fast neutrino flavor instabilities in two three-dimensional state-of-the-art core-collapse supernova simulations performed using a two-moment three-species neutrino transport scheme: one with an exploding 9$\mathrm{M_{\odot}}$ and one with a non-exploding 20$\mathrm{M_{\odot}}$ model. Apart from confirming the presence of fast instabilities occurring within the neutrino decoupling and the supernova pre-shock regions, we detect flavor instabilities in the post-shock region for the exploding model. These instabilities are likely to be scattering-induced. In addition, the failure in achieving a successful explosion in the heavier supernova model seems to seriously hinder the occurrence of fast instabilities in the post-shock region. This is a consequence of the large matter densities behind the stalled or retreating shock, which implies high neutrino scattering rates and thus more isotropic distributions of neutrinos and antineutrinos. Our findings suggest that the supernova model properties and the fate of the explosion can remarkably affect the occurrence of fast instabilities. Hence, a larger set of realistic hydrodynamical simulations of the stellar collapse is needed in order to make reliable predictions on the flavor conversion physics. △ Less

Submitted 11 March, 2021; v1 submitted 11 December, 2020; originally announced December 2020.

Comments: 12 pages, 7 figures; Minor changes to match the published version in PRD

Report number: MPP-2020-224

Journal ref: Phys. Rev. D 103, 063033 (2021)

Self-consistent 3D Supernova Models From -7 Minutes to +7 Seconds: a 1-bethe Explosion of a ~19 Solar-mass Progenitor

Authors: R. Bollig, N. Yadav, D. Kresse, H. -Th. Janka, B. Mueller, A. Heger

Abstract: To date, modern three-dimensional (3D) supernova (SN) simulations have not demonstrated that explosion energies of 10^{51} erg (=1 bethe = 1B) or more are possible for neutrino-driven SNe of non/slow-rotating M < 20 solar-mass progenitors. We present the first such model, considering a non-rotating, solar-metallicity 18.88 solar-mass progenitor, whose final 7 minutes of convective oxygen-shell bur… ▽ More To date, modern three-dimensional (3D) supernova (SN) simulations have not demonstrated that explosion energies of 10^{51} erg (=1 bethe = 1B) or more are possible for neutrino-driven SNe of non/slow-rotating M < 20 solar-mass progenitors. We present the first such model, considering a non-rotating, solar-metallicity 18.88 solar-mass progenitor, whose final 7 minutes of convective oxygen-shell burning were simulated in 3D and showed a violent oxygen-neon shell merger prior to collapse. A large set of 3D SN-models was computed with the Prometheus-Vertex code, whose improved convergence of the two-moment equations with Boltzmann closure allows now to fully exploit the implicit neutrino-transport treatment. Nuclear burning is treated with a 23-species network. We vary the angular grid resolution and consider different nuclear equations of state and muon formation in the proto-neutron star (PNS), which requires six-species transport with coupling of all neutrino flavors across all energy-momentum groups. Elaborate neutrino transport was applied until ~2 seconds after bounce. In one case the simulation was continued to >7 seconds with an approximate treatment of neutrino effects that allows for seamless continuation without transients. A spherically symmetric neutrino-driven wind does not develop. Instead, accretion downflows to the PNS and outflows of neutrino-heated matter establish a monotonic rise of the explosion energy until ~7 seconds post bounce, when the outgoing shock reaches about 50,000 km and enters the He-layer. The converged value of the explosion energy at infinity (with overburden subtracted) is roughly 1B and the ejected 56Ni mass up to 0.087 solar masses, both within a few 10 percent of the SN 1987A values. The final NS mass and kick are about 1.65 solar masses and over 450 km/s, respectively. △ Less

Submitted 15 April, 2021; v1 submitted 20 October, 2020; originally announced October 2020.

Comments: 29 pages, 17 figures, 1 table, 3 appendices; considerably extended due to comments by referee and readers; accepted by ApJ

The fully developed remnant of a neutrino-driven supernova: Evolution of ejecta structure and asymmetries in SNR Cassiopeia A

Authors: S. Orlando, A. Wongwathanarat, H. -T. Janka, M. Miceli, M. Ono, S. Nagataki, F. Bocchino, G. Peres

Abstract: Abridged. We aim at exploring to which extent the remnant keeps memory of the asymmetries that develop stochastically in the neutrino-heating layer due to hydrodynamic instabilities (e.g., convective overturn and the standing accretion shock instability) during the first second after core bounce. We coupled a 3D HD model of a neutrino-driven SN explosion with 3D MHD/HD simulations of the remnant f… ▽ More Abridged. We aim at exploring to which extent the remnant keeps memory of the asymmetries that develop stochastically in the neutrino-heating layer due to hydrodynamic instabilities (e.g., convective overturn and the standing accretion shock instability) during the first second after core bounce. We coupled a 3D HD model of a neutrino-driven SN explosion with 3D MHD/HD simulations of the remnant formation. The simulations cover 2000 years of expansion and include all physical processes relevant to describe the complexities in the SN evolution and the subsequent interaction of the stellar debris with the wind of the progenitor star. The interaction of large-scale asymmetries left from the earliest phases of the explosion with the reverse shock produces, at the age of $\approx 350$~years, an ejecta structure and a remnant morphology which are remarkably similar to those observed in Cas A. Small-scale structures in the large-scale Fe-rich plumes created during the initial stages of the SN, combined with HD instabilities that develop after the passage of the reverse shock, naturally produce a pattern of ring- and crown-like structures of shocked ejecta. The consequence is a spatial inversion of the ejecta layers with Si-rich ejecta being physically interior to Fe-rich ejecta. The full-fledged remnant shows voids and cavities in the innermost unshocked ejecta resulting from the expansion of Fe-rich plumes and their inflation due to the decay of radioactive species. The asymmetric distributions of $^{44}$Ti and $^{56}$Fe and their abundance ratio are both compatible with those inferred from high-energy observations of Chandra and NuSTAR. The main asymmetries observed in the ejecta distribution of Cas A can be explained by the interaction of the reverse shock with the large-scale asymmetries that developed from stochastic processes that originate during the first seconds of the SN blast. △ Less

Submitted 1 December, 2020; v1 submitted 3 September, 2020; originally announced September 2020.

Comments: 33 pages, 29 figures, 1 Table, 4 movies, 3 interactive graphics, accepted for publication on A&A, high quality figures will be available on the published paper

Journal ref: A&A 645, A66 (2021)

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.

Comments: 19 pages, 15 figures. Accepted by MNRAS, awaiting publication

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.

Comments: 47 pages, 27 figures, 6 tables; minor revisions, accepted by MNRAS

Properties of gamma-ray decay lines in 3D core-collapse supernova models, with application to SN 1987A and Cas A

Authors: A. Jerkstrand, A. Wongwathanarat, H. -T. Janka, M. Gabler, D. Alp, R. Diehl, K. Maeda, J. Larsson, C. Fransson, A. Menon, A. Heger

Abstract: Comparison of theoretical line profiles to observations provides important tests for supernova explosion models. We study the shapes of radioactive decay lines predicted by current 3D core-collapse explosion simulations, and compare these to observations of SN 1987A and Cas A. Both the widths and shifts of decay lines vary by several thousand kilometers per second depending on viewing angle. The l… ▽ More Comparison of theoretical line profiles to observations provides important tests for supernova explosion models. We study the shapes of radioactive decay lines predicted by current 3D core-collapse explosion simulations, and compare these to observations of SN 1987A and Cas A. Both the widths and shifts of decay lines vary by several thousand kilometers per second depending on viewing angle. The line profiles can be complex with multiple peaks. By combining observational constraints from 56Co decay lines, 44Ti decay lines, and Fe IR lines, we delineate a picture of the morphology of the explosive burning ashes in SN 1987A. For M_ZAMS=15-20 Msun progenitors exploding with ~1.5 *10^51 erg, ejecta structures suitable to reproduce the observations involve a bulk asymmetry of the 56Ni of at least ~400 km/s and a bulk velocity of at least ~1500 km/s. By adding constraints to reproduce the UVOIR bolometric light curve of SN 1987A up to 600d, an ejecta mass around 14 Msun is favoured. We also investigate whether observed decay lines can constrain the neutron star (NS) kick velocity. The model grid provides a constraint V_NS > V_redshift, and applying this to SN 1987A gives a NS kick of at least 500 km/s. For Cas A, our single model provides a satisfactory fit to the NuSTAR observations and reinforces the result that current neutrino-driven core-collapse SN models can achieve enough bulk asymmetry in the explosive burning material. Finally, we investigate the internal gamma-ray field and energy deposition, and compare the 3D models to 1D approximations. △ Less

Submitted 8 April, 2020; v1 submitted 11 March, 2020; originally announced March 2020.

Comments: 29 pages, published in MNRAS

The Birth Function for Black Holes and Neutron Stars in Close Binaries

Authors: S. Woosley, Tuguldur Sukhbold, H. -T. Janka

Abstract: The mass function for black holes and neutron stars at birth is explored for mass-losing helium stars. These should resemble, more closely than similar studies of single hydrogen-rich stars, the results of evolution in close binary systems. The effects of varying the mass-loss rate and metallicity are calculated using a simple semi-analytic approach to stellar evolution that is tuned to reproduce… ▽ More The mass function for black holes and neutron stars at birth is explored for mass-losing helium stars. These should resemble, more closely than similar studies of single hydrogen-rich stars, the results of evolution in close binary systems. The effects of varying the mass-loss rate and metallicity are calculated using a simple semi-analytic approach to stellar evolution that is tuned to reproduce detailed numerical calculations. Though the total fraction of black holes made in stellar collapse events varies considerably with metallicity, mass-loss rate, and mass cutoff, from 5$\%$ to 30$\%$, the shapes of their birth functions are very similar for all reasonable variations in these quantities. Median neutron star masses are in the range 1.32 - 1.37 $M_\odot$ regardless of metallicity. The median black hole mass for solar metallicity is typically 8 to 9 $M_\odot$ if only initial helium cores below 40 $M_\odot$ (ZAMS mass less than 80 $M_\odot$) are counted, and 9 - 13 $M_\odot$, in most cases, if helium cores with initial masses up to 150 $M_\odot$ (ZAMS mass less than 300 $M_\odot$) contribute. As long as the mass-loss rate as a function of mass exhibits no strong non-linearities, the black hole birth function from 15 to 35 $M_\odot$ has a slope that depends mostly on the initial mass function for main sequence stars. These findings imply the possibility of constraining the initial mass function and the properties of mass loss in close binaries using ongoing measurements of gravitational wave radiation. The expected rotation rates of the black holes are briefly discussed. △ Less

Submitted 5 April, 2020; v1 submitted 28 January, 2020; originally announced January 2020.

Comments: submitted to ApJ

Fast Neutrino Flavor Instability in the Neutron-star Convection Layer of Three-dimensional Supernova Models

Authors: Robert Glas, H. -Thomas Janka, Francesco Capozzi, Manibrata Sen, Basudeb Dasgupta, Alessandro Mirizzi, Guenter Sigl

Abstract: Neutrinos from a supernova (SN) might undergo fast flavor conversions near the collapsed stellar core. We perform a detailed study of this intriguing possibility, analyzing time-dependent state-of-the-art 3D SN models of 9 and 20 Msun. Both models were computed with multi-D three-flavor neutrino transport based on a two-moment solver, and both exhibit the presence of the lepton-number emission sel… ▽ More Neutrinos from a supernova (SN) might undergo fast flavor conversions near the collapsed stellar core. We perform a detailed study of this intriguing possibility, analyzing time-dependent state-of-the-art 3D SN models of 9 and 20 Msun. Both models were computed with multi-D three-flavor neutrino transport based on a two-moment solver, and both exhibit the presence of the lepton-number emission self-sustained asymmetry (LESA). The transport solution does not provide the angular distributions of the neutrino fluxes, which are crucial to track the fast flavor instability. To overcome this limitation, we use a recently proposed approach based on the angular moments of the energy-integrated electron lepton-number distribution. With this method we find the possibility of fast neutrino flavor instability at radii <~20 km, which is well interior to the neutrinosphere. Our results confirm recent observations in a 2D SN model and in 2D/3D models with fixed matter background, which were computed with Boltzmann neutrino transport. However, the flavor unstable locations are not isolated points as discussed previously, but thin skins surrounding volumes where electron antineutrinos are more abundant than electron neutrinos. These volumes grow with time and appear first in the convective layer of the proto-neutron star (PNS), where a decreasing electron fraction (Ye) and high temperatures favor the occurrence of regions with negative neutrino chemical potential. Since Ye remains higher in the LESA dipole direction, where convective lepton-number transport out from the nonconvective PNS core slows down the deleptonization, flavor unstable conditions become more widespread in the opposite hemisphere. This interesting phenomenon deserves further investigation, since its impact on SN modeling and possible consequences for SN dynamics and neutrino observations are presently unclear. (abridged) △ Less

Submitted 17 January, 2020; v1 submitted 30 November, 2019; originally announced December 2019.

Comments: 21 pages, 13 figures; more detailed explanations added due to referee comments; accepted by PRD

Report number: MPP-2019-239

Journal ref: Phys. Rev. D 101, 063001 (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.

Comments: 50 pages, 36 figures, 10 tables; submitted to ApJ

J0453+1559: a neutron star-white dwarf binary from a thermonuclear electron-capture supernova?

Authors: Thomas M. Tauris, H. -Thomas Janka

Abstract: The compact binary radio pulsar system J0453+1559 (Martinez et al. 2015) consists of a recycled pulsar as primary component of $1.559(5)\,M_\odot$ and an unseen companion star of $1.174(4)\,M_\odot$. Because of the relatively large orbital eccentricity of $e = 0.1125$, it was argued that the companion is a neutron star, making it the neutron star with the lowest accurately determined mass to date.… ▽ More The compact binary radio pulsar system J0453+1559 (Martinez et al. 2015) consists of a recycled pulsar as primary component of $1.559(5)\,M_\odot$ and an unseen companion star of $1.174(4)\,M_\odot$. Because of the relatively large orbital eccentricity of $e = 0.1125$, it was argued that the companion is a neutron star, making it the neutron star with the lowest accurately determined mass to date. However, a direct observational determination of the nature of the companion is currently not feasible. Moreover, state-of-the-art stellar evolution and supernova modeling are contradictive concerning the possibility to produce such a low-mass neutron star remnant. Here we challenge the neutron star interpretation by reasoning that the lower-mass component could instead be a white dwarf born in a thermonuclear electron-capture supernova (tECSN) event, in which oxygen-neon deflagration in the degenerate stellar core of an ultra-stripped progenitor ejects several $0.1\,M_\odot$ of matter and leaves a bound ONeFe white dwarf as the second-formed compact remnant. We determine the ejecta mass and remnant kick needed in this scenario to explain the properties of PSR~J0453+1559 by a neutron star--white dwarf system. More work on tECSNe is needed to assess the viability of this scenario. △ Less

Submitted 12 November, 2019; v1 submitted 26 September, 2019; originally announced September 2019.

Comments: ApJL, accepted for publication (6 pages, 4 figures)

Resolution Study for Three-dimensional Supernova Simulations with the Prometheus-Vertex Code

Authors: Tobias Melson, Daniel Kresse, H. -Thomas Janka

Abstract: We present a carefully designed, systematic study of the angular resolution dependence of simulations with the Prometheus-Vertex neutrino-hydrodynamics code. Employing a simplified neutrino heating-cooling scheme in the Prometheus hydrodynamics module allows us to sample the angular resolution between 4 degrees and 0.5 degrees. With a newly-implemented static mesh refinement (SMR) technique on the… ▽ More We present a carefully designed, systematic study of the angular resolution dependence of simulations with the Prometheus-Vertex neutrino-hydrodynamics code. Employing a simplified neutrino heating-cooling scheme in the Prometheus hydrodynamics module allows us to sample the angular resolution between 4 degrees and 0.5 degrees. With a newly-implemented static mesh refinement (SMR) technique on the Yin-Yang grid, the angular coordinates can be refined in concentric shells, compensating for the diverging structure of the spherical grid. In contrast to previous studies with Prometheus and other codes, we find that higher angular resolution and therefore lower numerical viscosity provides more favorable explosion conditions and faster shock expansion. We discuss the possible reasons for the discrepant results. The overall dynamics seem to converge at a resolution of about 1 degree. Applying the SMR setup to marginally exploding progenitors is disadvantageous for the shock expansion, however, because kinetic energy of downflows is dissipated to internal energy at resolution interfaces, leading to a loss of turbulent pressure support and a steeper temperature gradient. We also present a way to estimate the numerical viscosity on grounds of the measured turbulent kinetic-energy spectrum, leading to smaller values that are better compatible with the flow behavior witnessed in our simulations than results following calculations in previous literature. Interestingly, the numerical Reynolds numbers in the turbulent, neutrino-heated postshock layer (some 10 to several 100) are in the ballpark of expected neutrino-drag effects on the relevant length scales in the turbulent postshock layer. We provide a formal derivation and quantitative assessment of the neutrino drag terms in an appendix. △ Less

Submitted 3 February, 2020; v1 submitted 2 April, 2019; originally announced April 2019.

Comments: 37 pages, 14 figures, 4 tables; revised version with neutrino drag discussion extended for numerical evaluation; accepted by ApJ

NADA-FLD: A General Relativistic, Multi-dimensional Neutrino-hydrodynamics Code Employing Flux-limited Diffusion

Authors: Ninoy Rahman, Oliver Just, H. -Thomas Janka

Abstract: We present the new code NADA-FLD to solve multi-dimensional neutrino-hydrodynamics in full general relativity (GR) in spherical polar coordinates. The energy-dependent neutrino transport assumes the flux-limited diffusion (FLD) approximation and evolves the neutrino energy densities measured in the frame comoving with the fluid. Operator splitting is used to avoid multi-dimensional coupling of gri… ▽ More We present the new code NADA-FLD to solve multi-dimensional neutrino-hydrodynamics in full general relativity (GR) in spherical polar coordinates. The energy-dependent neutrino transport assumes the flux-limited diffusion (FLD) approximation and evolves the neutrino energy densities measured in the frame comoving with the fluid. Operator splitting is used to avoid multi-dimensional coupling of grid cells in implicit integration steps involving matrix inversions. Terms describing lateral diffusion and advection are integrated explicitly using the Allen-Cheng or the Runge-Kutta-Legendre method, which remain stable even in the optically thin regime. We discuss several toy-model problems in one and two dimensions to test the basic functionality and individual components of the transport scheme. We also perform fully dynamic core-collapse supernova (CCSN) simulations in spherical symmetry. For a Newtonian model we find good agreement with the M1 code ALCAR, and for a GR model we reproduce the main effects of GR in CCSNe already found by previous works. △ Less

Submitted 2 October, 2019; v1 submitted 29 January, 2019; originally announced January 2019.

Comments: 28 pages, 11 figures; accepted by MNRAS; supplemented by RKL method, resolution tests, and appendix with detailed derivations of FLD flux in response to referee comments

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.

Comments: 16 pages, 9 figures, 3 tables; accepted for publication in Astron. Astrophys

Journal ref: A&A 624, A116 (2019)

Gravitational waves from three-dimensional core-collapse supernova models: The impact of moderate progenitor rotation

Authors: H. Andresen, E. Müller, H. -Th. Janka, A. Summa, K. Gill, M. Zanolin

Abstract: We present predictions for the gravitational-wave (GW) emission of three-dimensional supernova (SN) simulations performed for a 15 solar-mass progenitor with the Prometheus-Vertex code using energy-dependent, three-flavor neutrino transport. The progenitor adopted from stellar evolution calculations including magnetic fields had a fairly low specific angular momentum (j_Fe <~ 10^{15} cm^2/s) in th… ▽ More We present predictions for the gravitational-wave (GW) emission of three-dimensional supernova (SN) simulations performed for a 15 solar-mass progenitor with the Prometheus-Vertex code using energy-dependent, three-flavor neutrino transport. The progenitor adopted from stellar evolution calculations including magnetic fields had a fairly low specific angular momentum (j_Fe <~ 10^{15} cm^2/s) in the iron core (central angular velocity ~0.2 rad/s), which we compared to simulations without rotation and with artificially enhanced rotation (j_Fe <~ 2*10^{16} cm^2/s; central angular velocity ~0.5 rad/s). Our results confirm that the time-domain GW signals of SNe are stochastic, but possess deterministic components with characteristic patterns at low frequencies (<~200 Hz), caused by mass motions due to the standing accretion shock instability (SASI), and at high frequencies, associated with gravity-mode oscillations in the surface layer of the proto-neutron star (PNS). Non-radial mass motions in the post-shock layer as well as PNS convection are important triggers of GW emission, whose amplitude scales with the power of the hydrodynamic flows. There is no monotonic increase of the GW amplitude with rotation, but a clear correlation with the strength of SASI activity. Our slowly rotating model is a fainter GW emitter than the non-rotating model because of weaker SASI activity and damped convection in the post-shock layer and PNS. In contrast, the faster rotating model exhibits a powerful SASI spiral mode during its transition to explosion, producing the highest GW amplitudes with a distinctive drift of the low-frequency emission peak from ~80-100 Hz to ~40-50 Hz. This migration signifies shock expansion, whereas non-exploding models are discriminated by the opposite trend. △ Less

Submitted 13 May, 2019; v1 submitted 17 October, 2018; originally announced October 2018.

Comments: Added new figure, figure 9. Updated figure 9, now figure 10. Modified the discussion of the proto-neutron star convection. Added a figure showing the average rotation rate as a function of radius. Added a section discussing where the low-frequency gravitational waves are generated, this information is visualized in figure 9. We also made some minor changes to the text and selected plots

Effects of LESA in Three-Dimensional Supernova Simulations with Multi-Dimensional and Ray-by-Ray-plus Neutrino Transport

Authors: Robert Glas, H. -Thomas Janka, Tobias Melson, Georg Stockinger, Oliver Just

Abstract: A set of eight self-consistent, time-dependent supernova (SN) simulations in three spatial dimensions (3D) for 9 solar-mass and 20 solar-mass progenitors is evaluated for the presence of dipolar asymmetries of the electron lepton-number emission as discovered by Tamborra et al. and termed lepton-number emission self-sustained asymmetry (LESA). The simulations were performed with the Aenus-Alcar ne… ▽ More A set of eight self-consistent, time-dependent supernova (SN) simulations in three spatial dimensions (3D) for 9 solar-mass and 20 solar-mass progenitors is evaluated for the presence of dipolar asymmetries of the electron lepton-number emission as discovered by Tamborra et al. and termed lepton-number emission self-sustained asymmetry (LESA). The simulations were performed with the Aenus-Alcar neutrino/hydrodynamics code, which treats the energy- and velocity-dependent transport of neutrinos of all flavors by a two-moment scheme with algebraic M1 closure. For each of the progenitors, results with fully multi-dimensional (FMD) neutrino transport and with ray-by-ray-plus (RbR+) approximation are considered for two different grid resolutions. While the 9 solar-mass models develop explosions, the 20 solar-mass progenitor does not explode with the employed version of simplified neutrino opacities. In all 3D models we observe the growth of substantial dipole amplitudes of the lepton-number (electron neutrino minus antineutrino) flux with stable or slowly time-evolving direction and overall properties fully consistent with the LESA phenomenon. Models with RbR+ transport develop LESA dipoles somewhat faster and with temporarily higher amplitudes, but the FMD calculations exhibit cleaner hemispheric asymmetries with a far more dominant dipole. In contrast, the RbR+ results display much wider multipole spectra of the neutrino-emission anisotropies with significant power also in the quadrupole and higher-order modes. Our results disprove speculations that LESA is a numerical artifact of RbR+ transport. We also discuss LESA as consequence of a dipolar convection flow inside of the nascent neutron star and establish, tentatively, a connection to Chandrasekhar's linear theory of thermal instability in spherical shells. △ Less

Submitted 7 June, 2019; v1 submitted 26 September, 2018; originally announced September 2018.

Comments: 20 pages, 9 figures; revised version accepted by ApJ; new Figs. 6,7, and new panels in Fig.8 added; Sects. 4,5,6 considerably extended in reply to referee questions

Three-Dimensional Core-Collapse Supernova Simulations with Multi-Dimensional Neutrino Transport Compared to the Ray-by-Ray-plus Approximation

Authors: Robert Glas, Oliver Just, H. -Thomas Janka, Martin Obergaulinger

Abstract: Self-consistent, time-dependent supernova (SN) simulations in three spatial dimensions (3D) are conducted with the Aenus-Alcar code, comparing, for the first time, calculations with fully multi-dimensional (FMD) neutrino transport and the ray-by-ray-plus (RbR+) approximation, both based on a two-moment solver with algebraic M1 closure. We find good agreement between 3D results with FMD and RbR+ tr… ▽ More Self-consistent, time-dependent supernova (SN) simulations in three spatial dimensions (3D) are conducted with the Aenus-Alcar code, comparing, for the first time, calculations with fully multi-dimensional (FMD) neutrino transport and the ray-by-ray-plus (RbR+) approximation, both based on a two-moment solver with algebraic M1 closure. We find good agreement between 3D results with FMD and RbR+ transport for both tested grid resolutions in the cases of a 20 solar-mass progenitor, which does not explode with the employed simplified set of neutrino opacities, and of an exploding 9 solar-mass model. This is in stark contrast to corresponding axisymmetric (2D) simulations, which confirm previous claims that the RbR+ approximation can foster explosions in 2D in particular in models with powerful axial sloshing of the stalled shock due to the standing accretion shock instability (SASI). However, while local and instantaneous variations of neutrino fluxes and heating rates can still be considerably higher with RbR+ transport in 3D, the time-averaged quantities are very similar to FMD results because of the absence of a fixed, artificial symmetry axis that channels the flow. Therefore, except for stochastic fluctuations, the neutrino signals and the post-bounce evolution of 3D simulations with FMD and RbR+ transport are also very similar, in particular for our calculations with the better grid resolution. Higher spatial resolution has clearly a more important impact than the differences by the two transport treatments. Our results back up the use of the RbR+ approximation for neutrino transport in 3D SN modeling. △ Less

Submitted 5 February, 2019; v1 submitted 26 September, 2018; originally announced September 2018.

Comments: 25 pages, 16 figures; referee comments included, new appendix added; accepted by ApJ

Improved Leakage-Equilibration-Absorption Scheme (ILEAS) for Neutrino Physics in Compact Object Mergers

Authors: Ricard Ardevol-Pulpillo, H. -Thomas Janka, Oliver Just, Andreas Bauswein

Abstract: We present a new, computationally efficient, energy-integrated approximation for neutrino effects in hot and dense astrophysical environments such as supernova cores and compact binary mergers and their remnants. Our new method, termed ILEAS for Improved Leakage-Equilibration-Absorption Scheme, improves the lepton-number and energy losses of traditional leakage descriptions by a novel prescription… ▽ More We present a new, computationally efficient, energy-integrated approximation for neutrino effects in hot and dense astrophysical environments such as supernova cores and compact binary mergers and their remnants. Our new method, termed ILEAS for Improved Leakage-Equilibration-Absorption Scheme, improves the lepton-number and energy losses of traditional leakage descriptions by a novel prescription of the diffusion time-scale based on a detailed energy integral of the flux-limited diffusion equation. The leakage module is supplemented by a neutrino-equilibration treatment that ensures the proper evolution of the total lepton number and medium plus neutrino energies as well as neutrino-pressure effects in the neutrino-trapping domain. Moreover, we employ a simple and straightforwardly applicable ray-tracing algorithm for including re-absorption of escaping neutrinos especially in the decoupling layer and during the transition to semi-transparent conditions. ILEAS is implemented on a three-dimensional (3D) Cartesian grid with a minimum of free and potentially case-dependent parameters and exploits the basic physics constraints that should be fulfilled in the neutrino-opaque and free-streaming limits. We discuss a suite of tests for stationary and time-dependent proto-neutron star models and post-merger black-hole-torus configurations, for which 3D ILEAS results are demonstrated to agree with energy-dependent 1D and 2D two-moment (M1) neutrino transport on the level of 10--15 percent in basic neutrino properties. This also holds for the radial profiles of the neutrino luminosities and of the electron fraction. Even neutrino absorption maps around torus-like neutrino sources are qualitatively similar without any fine-tuning, confirming that ILEAS can satisfactorily reproduce local losses and re-absorption of neutrinos as found in sophisticated transport calculations. △ Less

Submitted 28 February, 2019; v1 submitted 31 July, 2018; originally announced August 2018.

Comments: 39 pages, 19 figures; revised version incl. referee requests; new Sect.3.4 with first NS-NS merger models added; accepted by MNRAS

Parallelized Solution Method of the Three-dimensional Gravitational Potential on the Yin-Yang Grid

Authors: Marius Almanstoetter, Tobias Melson, H. -Thomas Janka, Ewald Mueller

Abstract: We present a new method for solving the three-dimensional gravitational potential of a density field on the Yin-Yang grid. Our algorithm is based on a multipole decomposition and completely symmetric with respect to the two Yin-Yang grid patches. It is particularly efficient on distributed-memory machines with a large number of compute tasks, because the amount of data being explicitly communicate… ▽ More We present a new method for solving the three-dimensional gravitational potential of a density field on the Yin-Yang grid. Our algorithm is based on a multipole decomposition and completely symmetric with respect to the two Yin-Yang grid patches. It is particularly efficient on distributed-memory machines with a large number of compute tasks, because the amount of data being explicitly communicated is minimized. All operations are performed on the original grid without the need for interpolating data onto an auxiliary spherical mesh. △ Less

Submitted 12 July, 2018; v1 submitted 12 June, 2018; originally announced June 2018.

Comments: 8 pages, 4 figures; two minor additions after refereeing; accepted by ApJ

Journal ref: The Astrophysical Journal, Volume 863, Issue 2, article id. 142, 7 pp. (2018)

X-ray Absorption in Young Core-Collapse Supernova Remnants

Authors: Dennis Alp, Josefin Larsson, Claes Fransson, M. Gabler, Annop Wongwathanarat, H. -Thomas Janka

Abstract: The material expelled by core-collapse supernova (SN) explosions absorbs X-rays from the central regions. We use SN models based on three-dimensional neutrino-driven explosions to estimate optical depths to the center of the explosion, compare different progenitor models, and investigate the effects of explosion asymmetries. The optical depths below 2 keV for progenitors with a remaining hydrogen… ▽ More The material expelled by core-collapse supernova (SN) explosions absorbs X-rays from the central regions. We use SN models based on three-dimensional neutrino-driven explosions to estimate optical depths to the center of the explosion, compare different progenitor models, and investigate the effects of explosion asymmetries. The optical depths below 2 keV for progenitors with a remaining hydrogen envelope are expected to be high during the first century after the explosion due to photoabsorption. A typical optical depth is $100 t_4^{-2} E^{-2}$, where $t_4$ is the time since the explosion in units of 10 000 days (${\sim}$27 years) and $E$ the energy in units of keV. Compton scattering dominates above 50 keV, but the scattering depth is lower and reaches unity already at ${\sim}$1000 days at 1 MeV. The optical depths are approximately an order of magnitude lower for hydrogen-stripped progenitors. The metallicity of the SN ejecta is much higher than in the interstellar medium, which enhances photoabsorption and makes absorption edges stronger. These results are applicable to young SN remnants in general, but we explore the effects on observations of SN 1987A and the compact object in Cas A in detail. For SN 1987A, the absorption is high and the X-ray upper limits of ${\sim}$100 Lsun on a compact object are approximately an order of magnitude less constraining than previous estimates using other absorption models. The details are presented in an accompanying paper. For the central compact object in Cas A, we find no significant effects of our more detailed absorption model on the inferred surface temperature. △ Less

Submitted 19 June, 2018; v1 submitted 11 May, 2018; originally announced May 2018.

Comments: 21 pages, 7 figures. Accepted for publication in ApJ. Updated to match accepted version; added Section 2.5 on asymmetries and discussion on homologous expansion in preamble of Section 4

The 30-Year Search for the Compact Object in SN 1987A

Authors: Dennis Alp, Josefin Larsson, Claes Fransson, Remy Indebetouw, Anders Jerkstrand, Antero Ahola, David Burrows, Peter Challis, Phil Cigan, Aleksandar Cikota, Robert P. Kirshner, Jacco Th. van Loon, Seppo Mattila, C. -Y. Ng, Sangwook Park, Jason Spyromilio, S. E. Woosley, Maarten Baes, Patrice Bouchet, Roger A. Chevalier, Kari A. Frank, Bryan M. Gaensler, Haley L. Gomez, H. -Thomas Janka, Bruno Leibundgut , et al. (10 additional authors not shown)

Abstract: Despite more than 30 years of searches, the compact object in Supernova (SN) 1987A has not yet been detected. We present new limits on the compact object in SN 1987A using millimeter, near-infrared, optical, ultraviolet, and X-ray observations from ALMA, VLT, HST, and Chandra. The limits are approximately 0.1 mJy ($0.1\times 10^{-26}$ erg s$^{-1}$ cm$^{-2}$ Hz$^{-1}$) at 213 GHz, 1 Lsun (… ▽ More Despite more than 30 years of searches, the compact object in Supernova (SN) 1987A has not yet been detected. We present new limits on the compact object in SN 1987A using millimeter, near-infrared, optical, ultraviolet, and X-ray observations from ALMA, VLT, HST, and Chandra. The limits are approximately 0.1 mJy ($0.1\times 10^{-26}$ erg s$^{-1}$ cm$^{-2}$ Hz$^{-1}$) at 213 GHz, 1 Lsun ($6\times 10^{-29}$ erg s$^{-1}$ cm$^{-2}$ Hz$^{-1}$) in optical if our line-of-sight is free of ejecta dust, and $10^{36}$ erg s$^{-1}$ ($2\times 10^{-30}$ erg s$^{-1}$ cm$^{-2}$ Hz$^{-1}$) in 2-10 keV X-rays. Our X-ray limits are an order of magnitude less constraining than previous limits because we use a more realistic ejecta absorption model based on three-dimensional neutrino-driven SN explosion models (presented in an accompanying article). The allowed bolometric luminosity of the compact object is 22 Lsun if our line-of-sight is free of ejecta dust, or 138 Lsun if dust-obscured. Depending on assumptions, these values limit the effective temperature of a neutron star to <4-8 MK and do not exclude models, which typically are in the range 3-4 MK. For the simplest accretion model, the accretion rate for an efficiency $η$ is limited to $< 10^{-11} η^{-1}$ Msun yr$^{-1}$, which excludes most predictions. For pulsar activity modeled by a rotating magnetic dipole in vacuum, the limit on the magnetic field strength ($B$) for a given spin period ($P$) is $B < 10^{14} P^2$ G s$^{-2}$. By combining information about radiation reprocessing and geometry, it is likely that the compact object is a dust-obscured thermally-emitting neutron star, which may appear as a region of higher-temperature ejecta dust emission. △ Less

Submitted 30 July, 2018; v1 submitted 11 May, 2018; originally announced May 2018.

Comments: 35 pages, 10 figures. Accepted for publication in ApJ. Updated to match accepted version; slightly shortened and clarified the discussion

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.

Comments: Resubmitted to MNRAS after referee comments

Simulations of Electron Capture and Low-Mass Iron Core Supernovae

Authors: B. Müller, S. Wanajo, H. -Th. Janka, A. Heger, D. Gay, S. A. Sim

Abstract: The evolutionary pathways of core-collapse supernova progenitors at the low-mass end of the spectrum are beset with major uncertainties. In recent years, a variety of evolutionary channels has been discovered in addition to the classical electron capture supernova channel of super-AGB stars. The few available progenitor models at the low-mass end have been studied with great success in supernova s… ▽ More The evolutionary pathways of core-collapse supernova progenitors at the low-mass end of the spectrum are beset with major uncertainties. In recent years, a variety of evolutionary channels has been discovered in addition to the classical electron capture supernova channel of super-AGB stars. The few available progenitor models at the low-mass end have been studied with great success in supernova simulations as the peculiar density structure makes for robust neutrino-driven explosions in this mass range. Detailed nucleosynthesis calculations have been conducted both for models of electron capture supernovae and low-mass iron core supernovae and revealed an interesting production of the lighter trans-iron elements (such as Zn, Sr, Y, Zr) as well as rare isotopes like Ca-48 and Fe-60. We stress the need to explore the low-mass end of the supernova spectrum further and link various observables to understand the diversity of explosions in this regime. △ Less

Submitted 7 October, 2017; originally announced October 2017.

Comments: 7 page, 3 figures, proceedings of the conference "The AGB-Supernova Mass Transition", to appear in Memorie della Società Astronomica Italiana

Rotation-supported Neutrino-driven Supernova Explosions in Three Dimensions and the Critical Luminosity Condition

Authors: A. Summa, H. -Th. Janka, T. Melson, A. Marek

Abstract: We present the first self-consistent, three-dimensional (3D) core-collapse supernova simulations performed with the Prometheus-Vertex code for a rotating progenitor star. Besides using the angular momentum of the 15 solar-mass model as obtained in the stellar evolution calculation with an angular frequency of about 0.001 rad/s (spin period of more than 6000 s) at the Si/Si-O interface, we also com… ▽ More We present the first self-consistent, three-dimensional (3D) core-collapse supernova simulations performed with the Prometheus-Vertex code for a rotating progenitor star. Besides using the angular momentum of the 15 solar-mass model as obtained in the stellar evolution calculation with an angular frequency of about 0.001 rad/s (spin period of more than 6000 s) at the Si/Si-O interface, we also computed 2D and 3D cases with no rotation and with a ~300 times shorter rotation period and different angular resolutions. In 2D, only the nonrotating and slowly rotating models explode, while rapid rotation prevents an explosion within 500 ms after bounce because of lower radiated neutrino luminosities and mean energies and thus reduced neutrino heating. In contrast, only the fast rotating model develops an explosion in 3D when the Si/Si-O interface collapses through the shock. The explosion becomes possible by the support of a powerful SASI spiral mode, which compensates for the reduced neutrino heating and pushes strong shock expansion in the equatorial plane. Fast rotation in 3D leads to a "two-dimensionalization" of the turbulent energy spectrum (yielding roughly a -3 instead of a -5/3 power-law slope at intermediate wavelengths) with enhanced kinetic energy on the largest spatial scales. We also introduce a generalization of the "universal critical luminosity condition" of Summa et al. (2016) to account for the effects of rotation, and demonstrate its viability for a set of more than 40 core-collapse simulations including 9 and 20 solar-mass progenitors as well as black-hole forming cases of 40 and 75 solar-mass stars to be discussed in forthcoming papers. △ Less

Submitted 22 November, 2017; v1 submitted 14 August, 2017; originally announced August 2017.

Comments: 24 pages, 19 figures; refereed version with additional section on resolution dependence; accepted by ApJ

Formation of Double Neutron Star Systems

Authors: T. M. Tauris, M. Kramer, P. C. C. Freire, N. Wex, H. -T. Janka, N. Langer, Ph. Podsiadlowski, E. Bozzo, S. Chaty, M. U. Kruckow, E. P. J. van den Heuvel, J. Antoniadis, R. P. Breton, D. J. Champion

Abstract: Double neutron star (DNS) systems represent extreme physical objects and the endpoint of an exotic journey of stellar evolution and binary interactions. Large numbers of DNS systems and their mergers are anticipated to be discovered using the Square-Kilometre-Array searching for radio pulsars and high-frequency gravitational wave detectors (LIGO/VIRGO), respectively. Here we discuss all key proper… ▽ More Double neutron star (DNS) systems represent extreme physical objects and the endpoint of an exotic journey of stellar evolution and binary interactions. Large numbers of DNS systems and their mergers are anticipated to be discovered using the Square-Kilometre-Array searching for radio pulsars and high-frequency gravitational wave detectors (LIGO/VIRGO), respectively. Here we discuss all key properties of DNS systems, as well as selection effects, and combine the latest observational data with new theoretical progress on various physical processes with the aim of advancing our knowledge on their formation. We examine key interactions of their progenitor systems and evaluate their accretion history during the high-mass X-ray binary stage, the common envelope phase and the subsequent Case BB mass transfer, and argue that the first-formed NSs have accreted at most $\sim 0.02\;M_{\odot}$. We investigate DNS masses, spins and velocities, and in particular correlations between spin period, orbital period and eccentricity. Numerous Monte Carlo simulations of the second supernova (SN) events are performed to extrapolate pre-SN stellar properties and probe the explosions. All known close-orbit DNS systems are consistent with ultra-stripped exploding stars. Although their resulting NS kicks are often small, we demonstrate a large spread in kick magnitudes which may, in general, depend on the past interaction history of the exploding star and thus correlate with the NS mass. We analyze and discuss NS kick directions based on our SN simulations. Finally, we discuss the terminal evolution of close-orbit DNS systems until they merge and possibly produce a short $γ$-ray burst. △ Less

Submitted 18 September, 2017; v1 submitted 28 June, 2017; originally announced June 2017.

Comments: 62 pages (emulateapj), 40 figures, 7 tables, published in ApJ. Minor typos fixed, two figures updated and text on ULXs updated. Published version is available here: https://doi.org/10.3847/1538-4357/aa7e89

Journal ref: The Astrophysical Journal, 846, 170 (2017)

Very Deep Inside the SN 1987A Core Ejecta: Molecular Structures Seen in 3D

Authors: F. J. Abellán, R. Indebetouw, J. M. Marcaide, M. Gabler, C. Fransson, J. Spyromilio, D. N. Burrows, R. Chevalier, P. Cigan, B. M. Gaensler, H. L. Gomez, H. -Th. Janka, R. Kirshner, J. Larsson, P. Lundqvist, M. Matsuura, R. McCray, C. -Y. Ng, S. Park, P. Roche, L. Staveley-Smith, J. Th. Van Loon, J. C. Wheeler, S. E. Woosley

Abstract: Most massive stars end their lives in core-collapse supernova explosions and enrich the interstellar medium with explosively nucleosynthesized elements. Following core collapse, the explosion is subject to instabilities as the shock propagates outwards through the progenitor star. Observations of the composition and structure of the innermost regions of a core-collapse supernova provide a direct p… ▽ More Most massive stars end their lives in core-collapse supernova explosions and enrich the interstellar medium with explosively nucleosynthesized elements. Following core collapse, the explosion is subject to instabilities as the shock propagates outwards through the progenitor star. Observations of the composition and structure of the innermost regions of a core-collapse supernova provide a direct probe of the instabilities and nucleosynthetic products. SN 1987A in the Large Magellanic Cloud (LMC) is one of very few supernovae for which the inner ejecta can be spatially resolved but are not yet strongly affected by interaction with the surroundings. Our observations of SN 1987A with the Atacama Large Millimeter/submillimeter Array (ALMA) are of the highest resolution to date and reveal the detailed morphology of cold molecular gas in the innermost regions of the remnant. The 3D distributions of carbon and silicon monoxide (CO and SiO) emission differ, but both have a central deficit, or torus-like distribution, possibly a result of radioactive heating during the first weeks ("nickel heating"). The size scales of the clumpy distribution are compared quantitatively to models, demonstrating how progenitor and explosion physics can be constrained. △ Less

Submitted 14 June, 2017; originally announced June 2017.

Comments: ApJL accepted, 9 pages, 5 figures

Muon Creation in Supernova Matter Facilitates Neutrino-driven Explosions

Authors: R. Bollig, H. -Th. Janka, A. Lohs, G. Martinez-Pinedo, C. J. Horowitz, T. Melson

Abstract: Muons can be created in nascent neutron stars (NSs) due to the high electron chemical potentials and the high temperatures. Because of their relatively lower abundance compared to electrons, their role has so far been ignored in numerical simulations of stellar core collapse and NS formation. However, the appearance of muons softens the NS equation of state, triggers faster NS contraction and thus… ▽ More Muons can be created in nascent neutron stars (NSs) due to the high electron chemical potentials and the high temperatures. Because of their relatively lower abundance compared to electrons, their role has so far been ignored in numerical simulations of stellar core collapse and NS formation. However, the appearance of muons softens the NS equation of state, triggers faster NS contraction and thus leads to higher luminosities and mean energies of the emitted neutrinos. This strengthens the postshock heating by neutrinos and can facilitate explosions by the neutrino-driven mechanism. △ Less

Submitted 17 November, 2017; v1 submitted 14 June, 2017; originally announced June 2017.

Comments: 6 pages, 4 figures; one more model added upon referee request; accepted by PRL

Journal ref: Phys. Rev. Lett. 119, 242702 (2017)

Spatial distribution of radionuclides in 3D models of SN 1987A and Cas A

Authors: H. -Thomas Janka, Michael Gabler, Annop Wongwathanarat

Abstract: Fostered by the possibilities of multi-dimensional computational modeling, in particular the advent of three-dimensional (3D) simulations, our understanding of the neutrino-driven explosion mechanism of core-collapse supernovae (SNe) has experienced remarkable progress over the past decade. First self-consistent, first-principle models have shown successful explosions in 3D, and even failed cases… ▽ More Fostered by the possibilities of multi-dimensional computational modeling, in particular the advent of three-dimensional (3D) simulations, our understanding of the neutrino-driven explosion mechanism of core-collapse supernovae (SNe) has experienced remarkable progress over the past decade. First self-consistent, first-principle models have shown successful explosions in 3D, and even failed cases may be cured by moderate changes of the microphysics inside the neutron star (NS), better grid resolution, or more detailed progenitor conditions at the onset of core collapse, in particular large-scale perturbations in the convective Si and O burning shells. 3D simulations have also achieved to follow neutrino-driven explosions continuously from the initiation of the blast wave, through the shock breakout from the progenitor surface, into the radioactively powered evolution of the SN, and towards the free expansion phase of the emerging remnant. Here we present results from such simulations, which form the basis for direct comparisons with observations of SNe and SN remnants in order to derive constraints on the still disputed explosion mechanism. It is shown that predictions based on hydrodynamic instabilities and mixing processes associated with neutrino-driven explosions yield good agreement with measured NS kicks, light-curve properties of SN 1987A, and asymmetries of iron and 44Ti distributions observed in SN 1987A and Cassiopeia A. △ Less

Submitted 2 May, 2017; originally announced May 2017.

Comments: 9 pages, 6 figures; submitted to: "SN 1987A, 30 years later", Proceedings IAU Symposium No. 331, 2017; eds. M. Renaud et al

Supernova Simulations from a 3D Progenitor Model -- Impact of Perturbations and Evolution of Explosion Properties

Authors: B. Müller, T. Melson, A. Heger, H. -Th. Janka

Abstract: We study the impact of large-scale perturbations from convective shell burning on the core-collapse supernova explosion mechanism using three-dimensional (3D) multi-group neutrino hydrodynamics simulations of an 18 solar mass progenitor. Seed asphericities in the O shell, obtained from a recent 3D model of O shell burning, help trigger a neutrino-driven explosion 330ms after bounce whereas the sho… ▽ More We study the impact of large-scale perturbations from convective shell burning on the core-collapse supernova explosion mechanism using three-dimensional (3D) multi-group neutrino hydrodynamics simulations of an 18 solar mass progenitor. Seed asphericities in the O shell, obtained from a recent 3D model of O shell burning, help trigger a neutrino-driven explosion 330ms after bounce whereas the shock is not revived in a model based on a spherically symmetric progenitor for at least another 300ms. We tentatively infer a reduction of the critical luminosity for shock revival by ~20% due to pre-collapse perturbations. This indicates that convective seed perturbations play an important role in the explosion mechanism in some progenitors. We follow the evolution of the 18 solar mass model into the explosion phase for more than 2s and find that the cycle of accretion and mass ejection is still ongoing at this stage. With a preliminary value of 0.77 Bethe for the diagnostic explosion energy, a baryonic neutron star mass of 1.85 solar masses, a neutron star kick of ~600km/s and a neutron star spin period of ~20ms at the end of the simulation, the explosion and remnant properties are slightly atypical, but still lie comfortably within the observed distribution. Although more refined simulations and a larger survey of progenitors are still called for, this suggests that a solution to the problem of shock revival and explosion energies in the ballpark of observations are within reach for neutrino-driven explosions in 3D. △ Less

Submitted 28 July, 2017; v1 submitted 1 May, 2017; originally announced May 2017.

Comments: 23 pages, 22 figures, accepted for publication in MNRAS

Light curve analysis of ordinary type IIP supernovae based on neutrino-driven explosion simulations in three dimensions

Authors: V. P. Utrobin, A. Wongwathanarat, H. -Th. Janka, E. Mueller

Abstract: Type II-plateau supernovae (SNe IIP) are the most numerous subclass of core-collapse SNe originating from massive stars. In the framework of the neutrino-driven explosion mechanism, we study the SN outburst properties for a red supergiant progenitor model and compare the corresponding light curves with observations of the ordinary Type IIP SN 1999em. Three-dimensional (3D) simulations of (parametr… ▽ More Type II-plateau supernovae (SNe IIP) are the most numerous subclass of core-collapse SNe originating from massive stars. In the framework of the neutrino-driven explosion mechanism, we study the SN outburst properties for a red supergiant progenitor model and compare the corresponding light curves with observations of the ordinary Type IIP SN 1999em. Three-dimensional (3D) simulations of (parametrically triggered) neutrino-driven explosions are performed with the (explicit, finite-volume, Eulerian, multifluid hydrodynamics) code PROMETHEUS, using a presupernova model of a 15 Msun star as initial data. At approaching homologous expansion, the hydrodynamical and composition variables of the 3D models are mapped to a spherically symmetric configuration, and the simulations are continued with the (implicit, Lagrangian radiation-hydrodynamics) code CRAB to follow the blast-wave evolution during the SN outburst. Our 3D neutrino-driven explosion model with an explosion energy of about 0.5x10^51 erg produces Ni-56 in rough agreement with the amount deduced from fitting the radioactively powered light-curve tail of SN 1999em. The considered presupernova model, 3D explosion simulations, and light-curve calculations can explain the basic observational features of SN 1999em, except for those connected to the presupernova structure of the outer stellar layers. Our 3D simulations show that the distribution of Ni-rich matter in velocity space is asymmetric with a strong dipole component that is consistent with the observations of SN 1999em. The monotonic luminosity decline from the plateau to the radioactive tail in ordinary SNe IIP is a manifestation of the intense turbulent mixing at the He/H composition interface. △ Less

Submitted 8 August, 2017; v1 submitted 12 April, 2017; originally announced April 2017.

Comments: 16 pages, 13 figures, 2 tables; added figure, discussions, and references; accepted for publication in ApJ

Neutrino-driven Explosions

Authors: H. -Th. Janka

Abstract: The question why and how core-collapse supernovae (SNe) explode is one of the central and most long-standing riddles of stellar astrophysics. A solution is crucial for deciphering the SN phenomenon, for predicting observable signals such as light curves and spectra, nucleosynthesis, neutrinos, and gravitational waves, for defining the role of SNe in the evolution of galaxies, and for explaining th… ▽ More The question why and how core-collapse supernovae (SNe) explode is one of the central and most long-standing riddles of stellar astrophysics. A solution is crucial for deciphering the SN phenomenon, for predicting observable signals such as light curves and spectra, nucleosynthesis, neutrinos, and gravitational waves, for defining the role of SNe in the evolution of galaxies, and for explaining the birth conditions and properties of neutron stars (NSs) and stellar-mass black holes. Since the formation of such compact remnants releases over hundred times more energy in neutrinos than the SN in the explosion, neutrinos can be the decisive agents for powering the SN outburst. According to the standard paradigm of the neutrino-driven mechanism, the energy transfer by the intense neutrino flux to the medium behind the stagnating core-bounce shock, assisted by violent hydrodynamic mass motions (sometimes subsumed by the term "turbulence"), revives the outward shock motion and thus initiates the SN blast. Because of the weak coupling of neutrinos in the region of this energy deposition, detailed, multidimensional hydrodynamic models including neutrino transport and a wide variety of physics are needed to assess the viability of the mechanism. Owing to advanced numerical codes and increasing supercomputer power, considerable progress has been achieved in our understanding of the physical processes that have to act in concert for the success of neutrino-driven explosions. First studies begin to reveal observational implications and avenues to test the theoretical picture by data from individual SNe and SN remnants but also from population-integrated observables. While models will be further refined, a real breakthrough is expected through the next Galactic core-collapse SN, when neutrinos and gravitational waves can be used to probe the conditions deep inside the dying star. (abridged) △ Less

Submitted 28 February, 2017; originally announced February 2017.

Comments: Author version of chapter for 'Handbook of Supernovae,' edited by A. Alsabti and P. Murdin, Springer. 54 pages, 13 figures

Neutrino Emission from Supernovae

Authors: H. -Th. Janka

Abstract: Supernovae are the most powerful cosmic sources of MeV neutrinos. These elementary particles play a crucial role when the evolution of a massive star is terminated by the collapse of its core to a neutron star or a black hole and the star explodes as supernova. The release of electron neutrinos, which are abundantly produced by electron captures, accelerates the catastrophic infall and causes a gr… ▽ More Supernovae are the most powerful cosmic sources of MeV neutrinos. These elementary particles play a crucial role when the evolution of a massive star is terminated by the collapse of its core to a neutron star or a black hole and the star explodes as supernova. The release of electron neutrinos, which are abundantly produced by electron captures, accelerates the catastrophic infall and causes a gradual neutronization of the stellar plasma by converting protons to neutrons as dominant constituents of neutron star matter. The emission of neutrinos and antineutrinos of all flavors carries away the gravitational binding energy of the compact remnant and drives its evolution from the hot initial to the cold final state. The absorption of electron neutrinos and antineutrinos in the surroundings of the newly formed neutron star can power the supernova explosion and determines the conditions in the innermost supernova ejecta, making them an interesting site for the nucleosynthesis of iron-group elements and trans-iron nuclei. In this Chapter the basic neutrino physics in supernova cores and nascent neutron stars will be discussed. This includes the most relevant neutrino production, absorption, and scattering processes, elementary aspects of neutrino transport in dense environments, the characteristic neutrino emission phases with their typical signal features, and the perspectives connected to a measurement of the neutrino signal from a future galactic supernova. △ Less

Submitted 28 February, 2017; originally announced February 2017.

Comments: Author version of chapter for 'Handbook of Supernovae,' edited by A. Alsabti and P. Murdin, Springer. 30 pages, 9 figures

Neutron star kicks by the gravitational tug-boat mechanism in asymmetric supernova explosions: progenitor and explosion dependence

Authors: H. -Th. Janka

Abstract: Asymmetric mass ejection in the early phase of supernova (SN) explosions can impart a kick velocity to the new-born neutron star (NS). For neutrino-driven explosions the NS acceleration was shown to be mainly caused by the gravitational attraction of the anisotropically expelled inner ejecta, while hydrodynamic forces contribute on a subdominant level, and asymmetric neutrino emission plays only a… ▽ More Asymmetric mass ejection in the early phase of supernova (SN) explosions can impart a kick velocity to the new-born neutron star (NS). For neutrino-driven explosions the NS acceleration was shown to be mainly caused by the gravitational attraction of the anisotropically expelled inner ejecta, while hydrodynamic forces contribute on a subdominant level, and asymmetric neutrino emission plays only a secondary role. Two- and three-dimensional hydrodynamic simulations demonstrated that this gravitational tug-boat mechanism can explain the observed space velocities of young NSs up to more than 1000 km/s. Here, we discuss how the NS kick depends on the energy, ejecta mass, and asymmetry of the SN explosion, and which role the compactness of the pre-collapse stellar core plays for the momentum transfer to the NS. We also provide simple analytic expressions for the NS velocity in terms of these quantities. Referring to results of hydrodynamic simulations in the literature, we argue why within the discussed scenario of NS acceleration, electron-capture SNe, low-mass Fe-core SNe and ultra-stripped SNe can be expected to have considerably lower intrinsic NS kicks than core-collapse SNe of massive stellar cores. Our basic arguments remain valid also if progenitor stars possess large-scale asymmetries in their convective silicon and oxygen burning layers. Possible scenarios for spin-kick alignment are sketched. Much of our discussion remains on a conceptual and qualitative level and more work is needed on the numerical modeling side to determine the dependences of involved parameters, whose prescriptions will be needed for recipes that can be used to better describe NS kicks in binary evolution and population synthesis studies. △ Less

Submitted 16 February, 2017; v1 submitted 18 November, 2016; originally announced November 2016.

Comments: 11 pages, 1 figure; revised and accepted by ApJ; section on possibilities for spin-kick alignment added

Production and Distribution of $^{44}$Ti and $^{56}$Ni in a Three-dimensional Supernova Model Resembling Cassiopeia A

Authors: A. Wongwathanarat, H. -Th. Janka, E. Mueller, E. Pllumbi, S. Wanajo

Abstract: The spatial and velocity distributions of nuclear species synthesized in the innermost regions of core-collapse supernovae can yield important clues about explosion asymmetries and the operation of the still disputed explosion mechanism. Recent observations of radioactive $^{44}$Ti with high-energy satellite telescopes (Nuclear Spectroscopic Telescope Array [NuSTAR], INTEGRAL) have measured gamma-… ▽ More The spatial and velocity distributions of nuclear species synthesized in the innermost regions of core-collapse supernovae can yield important clues about explosion asymmetries and the operation of the still disputed explosion mechanism. Recent observations of radioactive $^{44}$Ti with high-energy satellite telescopes (Nuclear Spectroscopic Telescope Array [NuSTAR], INTEGRAL) have measured gamma-ray line details, which provide direct evidence of large-scale explosion asymmetries in SN 1987A and in Cassiopeia A (Cas A) even by mapping of the spatial brightness distribution (NuSTAR). Here we discuss a 3D simulation of a neutrino-driven explosion, using a parameterized neutrino engine, whose $^{44}$Ti distribution is mostly concentrated in one hemisphere pointing opposite to the neutron star (NS) kick velocity. Both exhibit intriguing resemblance to the observed morphology of the Cas A remnant, although neither the progenitor nor the explosion was fine-tuned for a perfect match. Our results demonstrate that the asymmetries observed in this remnant can, in principle, be accounted for by a neutrino-driven explosion, and that the high $^{44}$Ti abundance in Cas A may be explained without invoking rapid rotation or a jet-driven explosion, because neutrino-driven explosions generically eject large amounts of high-entropy matter. The recoil acceleration of the NS is connected to mass ejection asymmetries and is opposite to the direction of the stronger explosion, fully compatible with the gravitational tugboat mechanism. Our results also imply that Cas A and SN 1987A could possess similarly "one-sided" Ti and Fe asymmetries, with the difference that Cas A is viewed from a direction with large inclination angle to the NS motion, whereas the NS in SN 1987A should have a dominant velocity component pointing toward us. △ Less

Submitted 12 June, 2017; v1 submitted 14 October, 2016; originally announced October 2016.

Comments: 23 pages, 12 figures; revised version; discussion of new 3D analysis of 44Ti in Cas A and 3 new figures added; accepted by ApJ

Journal ref: ApJ 842, 13 (2017)