No Arabic abstract
Nonspherical mass motions are a generic feature of core-collapse supernovae, and hydrodynamic instabilities play a crucial role for the explosion mechanism. First successful neutrino-driven explosions could be obtained with self-consistent, first-principle simulations in three spatial dimensions (3D). But 3D models tend to be less prone to explosion than corresponding axisymmetric (2D) ones. This has been explained by 3D turbulence leading to energy cascading from large to small spatial scales, inversely to the 2D case, thus disfavoring the growth of buoyant plumes on the largest scales. Unless the inertia to explode simply reflects a lack of sufficient resolution in relevant regions, it suggests that some important aspect may still be missing for robust and sufficiently energetic neutrino-powered explosions. Such deficits could be associated with progenitor properties like rotation, magnetic fields or pre-collapse perturbations, or with microphysics that could lead to an enhancement of neutrino heating behind the shock. 3D simulations have also revealed new phenomena that are not present in 2D, for example spiral modes of the standing accretion shock instability (SASI) and a stunning dipolar lepton-emission self-sustained asymmetry (LESA). Both impose time- and direction-dependent variations on the detectable neutrino signal. The understanding of these effects and of their consequences is still in its infancy.
We develop a numerical code to calculate the neutrino transfer with multi-energy and multi-angle in three dimensions (3D) for the study of core-collapse supernovae. The numerical code solves the Boltzmann equations for neutrino distributions by the discrete-ordinate (S_n) method with a fully implicit differencing for time advance. The Boltzmann equations are formulated in the inertial frame with collision terms being evaluated to the zeroth order of v/c. A basic set of neutrino reactions for three neutrino species is implemented together with a realistic equation of state of dense matter. The pair process is included approximately in order to keep the system linear. We present numerical results for a set of test problems to demonstrate the ability of the code. The numerical treatments of advection and collision terms are validated first in the diffusion and free streaming limits. Then we compute steady neutrino distributions for a background extracted from a spherically symmetric, general relativistic simulation of 15Msun star and compare them with the results in the latter computation. We also demonstrate multi-D capabilities of the 3D code solving neutrino transfers for artificially deformed supernova cores in 2D and 3D. Formal solutions along neutrino paths are utilized as exact solutions. We plan to apply this code to the 3D neutrino-radiation hydrodynamics simulations of supernovae. This is the first article in a series of reports on the development.
The importance of detecting neutrinos from a Milky Way core-collapse supernova is well known. An under-studied phase is proto-neutron star cooling. For SN 1987A, this seemingly began at about 2 s, and is thus probed by only 6 of the 19 events (and only the $bar{ u}_e$ flavor) in the Kamiokande-II and IMB detectors. With the higher statistics expected for present and near-future detectors, it should be possible to measure detailed neutrino signals out to very late times. We present the first comprehensive study of neutrino detection during the proto-neutron star cooling phase, considering a variety of outcomes, using all flavors, and employing detailed detector physics. For our nominal model, the event yields (at 10 kpc) after 10 s -- the approximate duration of the SN 1987A signal -- far exceed the entire SN 1987A yield, with $simeq$250 $bar{ u}_e$ events (to 50 s) in Super-Kamiokande, $simeq$110 $ u_e$ events (to 40 s) in DUNE, and $simeq$10 $ u_mu, u_tau, bar{ u}_mu, bar{ u}_tau$ events (to 20 s) in JUNO. These data would allow unprecedented probes of the proto-neutron star, including the onset of neutrino transparency and hence its transition to a neutron star. If a black hole forms, even at very late times, this can be clearly identified. But will the detectors fulfill their potential for this perhaps once-ever opportunity for an all-flavor, high-statistics detection of a core collapse? Maybe. Further work is urgently needed, especially for DUNE to thoroughly investigate and improve its MeV capabilities.
Knowledge of the progenitors of core-collapse supernovae is a fundamental component in understanding the explosions. The recent progress in finding such stars is reviewed. The minimum initial mass that can produce a supernova has converged to 8 +/- 1 solar masses, from direct detections of red supergiant progenitors of II-P SNe and the most massive white dwarf progenitors, although this value is model dependent. It appears that most type Ibc supernovae arise from moderate mass interacting binaries. The highly energetic, broad-lined Ic supernovae are likely produced by massive, Wolf-Rayet progenitors. There is some evidence to suggest that the majority of massive stars above ~20 solar masses may collapse quietly to black-holes and that the explosions remain undetected. The recent discovery of a class of ultra-bright type II supernovae and the direct detection of some progenitor stars bearing luminous blue variable characteristics suggests some very massive stars do produce highly energetic explosions. The physical mechanism is open to debate and these SNe pose a challenge to stellar evolutionary theory.
We investigate the post-explosion phase in core-collapse supernovae with 2D hydrodynamical simulations and a simple neutrino treatment. The latter allows us to perform 46 simulations and follow the evolution of the 32 successful explosions during several seconds. We present a broad study based on three progenitors (11.2 $M_odot$, 15 $M_odot$, and 27 $M_odot$), different neutrino-heating efficiencies, and various rotation rates. We show that the first seconds after shock revival determine the final explosion energy, remnant mass, and properties of ejected matter. Our results suggest that a continued mass accretion increases the explosion energy even at late times. We link the late-time mass accretion to initial conditions such as rotation strength and shock deformation at explosion time. Only some of our simulations develop a neutrino-driven wind that survives for several seconds. This indicates that neutrino-driven winds are not a standard feature expected after every successful explosion. Even if our neutrino treatment is simple, we estimate the nucleosynthesis of the exploding models for the 15 $M_odot$ progenitor after correcting the neutrino energies and luminosities to get a more realistic electron fraction.
We present self-consistent, axisymmetric core-collapse supernova simulations performed with the Prometheus-Vertex code for 18 pre-supernova models in the range of 11-28 solar masses, including progenitors recently investigated by other groups. All models develop explosions, but depending on the progenitor structure, they can be divided into two classes. With a steep density decline at the Si/Si-O interface, the arrival of this interface at the shock front leads to a sudden drop of the mass-accretion rate, triggering a rapid approach to explosion. With a more gradually decreasing accretion rate, it takes longer for the neutrino heating to overcome the accretion ram pressure and explosions set in later. Early explosions are facilitated by high mass-accretion rates after bounce and correspondingly high neutrino luminosities combined with a pronounced drop of the accretion rate and ram pressure at the Si/Si-O interface. Because of rapidly shrinking neutron star radii and receding shock fronts after the passage through their maxima, our models exhibit short advection time scales, which favor the efficient growth of the standing accretion-shock instability. The latter plays a supportive role at least for the initiation of the re-expansion of the stalled shock before runaway. Taking into account the effects of turbulent pressure in the gain layer, we derive a generalized condition for the critical neutrino luminosity that captures the explosion behavior of all models very well. We validate the robustness of our findings by testing the influence of stochasticity, numerical resolution, and approximations in some aspects of the microphysics.