-
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.
-
Constraints on neutrino natal kicks from black-hole binary VFTS 243
Authors:
Alejandro Vigna-Gómez,
Reinhold Willcox,
Irene Tamborra,
Ilya Mandel,
Mathieu Renzo,
Tom Wagg,
Hans-Thomas Janka,
Daniel Kresse,
Julia Bodensteiner,
Tomer Shenar,
Thomas M. Tauris
Abstract:
The recently reported observation of VFTS 243 is the first example of a massive black-hole binary system with negligible binary interaction following black-hole formation. The black-hole mass ($\approx 10\ M_{\odot}$) and near-circular orbit ($e\approx 0.02$) of VFTS 243 suggest that the progenitor star experienced complete collapse, with energy-momentum being lost predominantly through neutrinos.…
▽ More
The recently reported observation of VFTS 243 is the first example of a massive black-hole binary system with negligible binary interaction following black-hole formation. The black-hole mass ($\approx 10\ M_{\odot}$) and near-circular orbit ($e\approx 0.02$) of VFTS 243 suggest that the progenitor star experienced complete collapse, with energy-momentum being lost predominantly through neutrinos. VFTS 243 enables us to constrain the natal kick and neutrino-emission asymmetry during black-hole formation. At 68% C.L., the natal kick velocity (mass decrement) is $\lesssim 10$ km/s ($\lesssim 1.0\ M_{\odot}$), with a full probability distribution that peaks when $\approx 0.3\ M_{\odot}$ were ejected, presumably in neutrinos, and the black hole experienced a natal kick of $4$ km/s. The neutrino-emission asymmetry is $\lesssim 4$%, with best fit values of $\sim$0-0.2%. Such a small neutrino natal kick accompanying black-hole formation is in agreement with theoretical predictions.
△ Less
Submitted 2 April, 2024; v1 submitted 2 October, 2023;
originally announced October 2023.
-
Production of 44Ti and iron-group nuclei in the ejecta of 3D neutrino-driven supernovae
Authors:
Andre Sieverding,
Daniel Kresse,
Hans-Thomas Janka
Abstract:
The radioactive isotopes of 44Ti and 56Ni are important products of explosive nucleosynthesis, which play a key role for supernova (SN) diagnostics and were detected in several nearby young SN remnants. However, most SN models based on non-rotating single stars predict yields of 44Ti that are much lower than the values inferred from observations. We present, for the first time, the nucleosynthesis…
▽ More
The radioactive isotopes of 44Ti and 56Ni are important products of explosive nucleosynthesis, which play a key role for supernova (SN) diagnostics and were detected in several nearby young SN remnants. However, most SN models based on non-rotating single stars predict yields of 44Ti that are much lower than the values inferred from observations. We present, for the first time, the nucleosynthesis yields from a self-consistent three-dimensional (3D) SN simulation of an approximately 19 Msun progenitor star that reaches an explosion energy comparable to that of SN 1987A and that covers the evolution of the neutrino-driven explosion until more than 7 seconds after core bounce. We find a significant enhancement of the Ti/Fe yield compared to recent spherically symmetric (1D) models and demonstrate that the long-time evolution is crucial to understand the efficient production of 44Ti due to the non-monotonic temperature and density histories of ejected mass elements. Additionally, we identify characteristic signatures of the nucleosynthesis in proton-rich ejecta, in particular high yields of 45Sc and 64Zn.
△ Less
Submitted 17 October, 2023; v1 submitted 18 August, 2023;
originally announced August 2023.
-
Low-luminosity type IIP supernovae: SN 2005cs and SN 2020cxd as very low-energy iron core-collapse explosions
Authors:
Alexandra Kozyreva,
Hans-Thomas Janka,
Daniel Kresse,
Stefan Taubenberger,
Petr Baklanov
Abstract:
SN 2020cxd is a representative of the family of low-energy, underluminous Type IIP supernovae (SNe), whose observations and analysis were recently reported by Yang et al. (2021). Here we re-evaluate the observational data for the diagnostic SN properties by employing the hydrodynamic explosion model of a 9 MSun red supergiant progenitor with an iron core and a pre-collapse mass of 8.75 Msun. The e…
▽ More
SN 2020cxd is a representative of the family of low-energy, underluminous Type IIP supernovae (SNe), whose observations and analysis were recently reported by Yang et al. (2021). Here we re-evaluate the observational data for the diagnostic SN properties by employing the hydrodynamic explosion model of a 9 MSun red supergiant progenitor with an iron core and a pre-collapse mass of 8.75 Msun. The explosion of the star was obtained by the neutrino-driven mechanism in a fully self-consistent simulation in three dimensions (3D). Multi-band light curves and photospheric velocities for the plateau phase are computed with the one-dimensional radiation-hydrodynamics code STELLA, applied to the spherically averaged 3D explosion model as well as spherisized radial profiles in different directions of the 3D model. We find that the overall evolution of the bolometric light curve, duration of the plateau phase, and basic properties of the multi-band emission can be well reproduced by our SN model with its explosion energy of only 0.7x10^50 erg and an ejecta mass of 7.4 Msun. These values are considerably lower than the previously reported numbers, but they are compatible with those needed to explain the fundamental observational properties of the prototype low-luminosity SN 2005cs. Because of the good compatibility of our photospheric velocities with line velocities determined for SN 2005cs, we conclude that the line velocities of SN 2020cxd are probably overestimated by up to a factor of about 3. The evolution of the line velocities of SN 2005cs compared to photospheric velocities in different explosion directions might point to intrinsic asymmetries in the SN ejecta.
△ Less
Submitted 31 May, 2022; v1 submitted 1 March, 2022;
originally announced March 2022.
-
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.
-
Stellar Collapse Diversity and the Diffuse Supernova Neutrino Background
Authors:
Daniel Kresse,
Thomas Ertl,
Hans-Thomas Janka
Abstract:
The diffuse cosmic supernova neutrino background (DSNB) is observational target of the gadolinium-loaded Super-Kamiokande (SK) detector and the forthcoming JUNO and Hyper-Kamiokande detectors. Current predictions are hampered by our still incomplete understanding of the supernova (SN) explosion mechanism and of the neutron star (NS) equation of state and maximum mass. In our comprehensive study we…
▽ More
The diffuse cosmic supernova neutrino background (DSNB) is observational target of the gadolinium-loaded Super-Kamiokande (SK) detector and the forthcoming JUNO and Hyper-Kamiokande detectors. Current predictions are hampered by our still incomplete understanding of the supernova (SN) explosion mechanism and of the neutron star (NS) equation of state and maximum mass. In our comprehensive study we revisit this problem on grounds of the landscapes of successful and failed SN explosions obtained by Sukhbold et al. and Ertl et al. with parametrized one-dimensional neutrino engines for large sets of single-star and helium-star progenitors, with the latter serving as proxy of binary evolution effects. Besides considering engines of different strengths, leading to different fractions of failed SNe with black-hole (BH) formation, we also vary the NS mass limit, the spectral shape of the neutrino emission, and include contributions from poorly understood alternative NS-formation channels such as accretion-induced or merger-induced collapse events. Since the neutrino signals of our large model sets are approximate, we calibrate the associated degrees of freedom by using state-of-the-art simulations of proto-neutron star cooling. Our predictions are higher than other recent ones because of a large fraction of failed SNe with long delay to BH formation. Our best-guess model predicts a DSNB electron-antineutrino-flux of 28.8^{+24.6}_{-10.9} cm^{-2}s^{-1} with 6.0^{+5.1}_{-2.1} cm^{-2}s^{-1} in the favorable measurement interval of [10,30] MeV, and 1.3^{+1.1}_{-0.4} cm^{-2}s^{-1} with electron-antineutrino energies > 17.3 MeV, which is roughly a factor of two below the current SK limit. The uncertainty range is dominated by the still insufficiently constrained cosmic rate of stellar core-collapse events.
△ Less
Submitted 18 December, 2020; v1 submitted 9 October, 2020;
originally announced October 2020.
-
Detecting the Diffuse Supernova Neutrino Background in the future Water-based Liquid Scintillator Detector Theia
Authors:
Julia Sawatzki,
Michael Wurm,
Daniel Kresse
Abstract:
A large-scale neutrino observatory based on Water-based Liquid Scintillator (WbLS) will be excellently suited for a measurement of the Diffuse Supernova Neutrino Background (DSNB). The WbLS technique offers high signal efficiency and effective suppression of the otherwise overwhelming background from neutral-current interactions of atmospheric neutrinos. To illustrate this, we investigate the DSNB…
▽ More
A large-scale neutrino observatory based on Water-based Liquid Scintillator (WbLS) will be excellently suited for a measurement of the Diffuse Supernova Neutrino Background (DSNB). The WbLS technique offers high signal efficiency and effective suppression of the otherwise overwhelming background from neutral-current interactions of atmospheric neutrinos. To illustrate this, we investigate the DSNB sensitivity for two configurations of the future Theia detector by developing the expected signal and background rejection efficiencies along a full analysis chain. Based on a statistical analysis of the remaining signal and background rates, we find that a rather moderate exposure of 190kt$\cdot$yrs will be sufficient to claim a ($5σ$) discovery of the faint DSNB signal for standard model assumptions. We conclude that, in comparison with other experimental techniques, WbLS offers the highest signal efficiency of more than 80\% and best signal significance over background.
△ Less
Submitted 27 January, 2021; v1 submitted 29 July, 2020;
originally announced July 2020.
-
Three-dimensional Models of Core-collapse Supernovae From Low-mass Progenitors With Implications for Crab
Authors:
G. Stockinger,
H. -Th. Janka,
D. Kresse,
T. Melson,
T. Ertl,
M. Gabler,
A. Gessner,
A. Wongwathanarat,
A. Tolstov,
S. -C. Leung,
K. Nomoto,
A. Heger
Abstract:
We present 3D full-sphere supernova simulations of non-rotating low-mass (~9 Msun) progenitors, covering the entire evolution from core collapse through bounce and shock revival, through shock breakout from the stellar surface, until fallback is completed several days later. We obtain low-energy explosions [~(0.5-1.0)x 10^{50} erg] of iron-core progenitors at the low-mass end of the core-collapse…
▽ More
We present 3D full-sphere supernova simulations of non-rotating low-mass (~9 Msun) progenitors, covering the entire evolution from core collapse through bounce and shock revival, through shock breakout from the stellar surface, until fallback is completed several days later. We obtain low-energy explosions [~(0.5-1.0)x 10^{50} erg] of iron-core progenitors at the low-mass end of the core-collapse supernova (LMCCSN) domain and compare to a super-AGB (sAGB) progenitor with an oxygen-neon-magnesium core that collapses and explodes as electron-capture supernova (ECSN). The onset of the explosion in the LMCCSN models is modelled self-consistently using the Vertex-Prometheus code, whereas the ECSN explosion is modelled using parametric neutrino transport in the Prometheus-HOTB code, choosing different explosion energies in the range of previous self-consistent models. The sAGB and LMCCSN progenitors that share structural similarities have almost spherical explosions with little metal mixing into the hydrogen envelope. A LMCCSN with less 2nd dredge-up results in a highly asymmetric explosion. It shows efficient mixing and dramatic shock deceleration in the extended hydrogen envelope. Both properties allow fast nickel plumes to catch up with the shock, leading to extreme shock deformation and aspherical shock breakout. Fallback masses of <~5x10^{-3} Msun have no significant effects on the neutron star (NS) masses and kicks. The anisotropic fallback carries considerable angular momentum, however, and determines the spin of the newly-born NS. The LMCCSNe model with less 2nd dredge-up results in a hydrodynamic and neutrino-induced NS kick of >40 km/s and a NS spin period of ~30 ms, both not largely different from those of the Crab pulsar at birth.
△ Less
Submitted 10 June, 2020; v1 submitted 5 May, 2020;
originally announced May 2020.
-
Neutrino emission characteristics of black hole formation in three-dimensional simulations of stellar collapse
Authors:
Laurie Walk,
Irene Tamborra,
Hans-Thomas Janka,
Alexander Summa,
Daniel Kresse
Abstract:
Neutrinos are unique probes of core-collapse supernova dynamics, especially in the case of black hole (BH) forming stellar collapses, where the electromagnetic emission may be faint or absent. By investigating two 3D hydrodynamical simulations of BH-forming stellar collapses of mass 40 and 75 M_sun, we identify the physical processes preceding BH formation through neutrinos, and forecast the neutr…
▽ More
Neutrinos are unique probes of core-collapse supernova dynamics, especially in the case of black hole (BH) forming stellar collapses, where the electromagnetic emission may be faint or absent. By investigating two 3D hydrodynamical simulations of BH-forming stellar collapses of mass 40 and 75 M_sun, we identify the physical processes preceding BH formation through neutrinos, and forecast the neutrino signal expected in the existing IceCube and Super-Kamiokande detectors, as well as in the future generation DUNE facility. Prior to the abrupt termination of the neutrino signal corresponding to BH formation, both models develop episodes of strong and long-lasting activity by the spiral standing accretion shock instability (SASI). We find that the spiral SASI peak in the Fourier power spectrum of the neutrino event rate will be distinguishable at 3 sigma above the detector noise for distances up to O(30) kpc in the most optimistic scenario, with IceCube having the highest sensitivity. Interestingly, given the long duration of the spiral SASI episodes, the spectrograms of the expected neutrino event rate carry clear signs of the evolution of the blue spiral SASI frequency as a function of time, as the shock radius and post-shock fluid velocity evolve. Due to the high accretion luminosity and its large-amplitude SASI-induced modulations, any contribution from asymmetric (dipolar or quadrupolar) neutrino emission associated with the lepton emission self-sustained asymmetry (LESA) is far subdominant in the neutrino signal.
△ Less
Submitted 25 May, 2020; v1 submitted 28 October, 2019;
originally announced October 2019.
-
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.