No Arabic abstract
We investigate correlated gravitational wave and neutrino signals from rotating core-collapse supernovae with simulations. Using an improved mode identification procedure based on mode function matching, we show that a linear quadrupolar mode of the core produces a dual imprint on gravitational waves and neutrinos in the early post-bounce phase of the supernova. The angular harmonics of the neutrino emission are consistent with the mode energy around the neutrinospheres, which points to a mechanism for the imprint on neutrinos. Thus, neutrinos carry information about the mode amplitude in the outer region of the core, whereas gravitational waves probe deeper in. We also find that the best-fit mode function has a frequency bounded above by $sim 420$ Hz, and yet the modes frequency in our simulations is $sim 15%$ higher, due to the use of Newtonian hydrodynamics and a widely used pseudo-Newtonian gravity approximation. This overestimation is particularly important for the analysis of gravitational wave detectability and asteroseismology, pointing to limitations of pseudo-Newtonian approaches for these purposes, possibly even resulting in excitation of incorrect modes. In addition, mode frequency matching (as opposed to mode function matching) could be resulting in mode misidentification in recent work. Lastly, we evaluate the prospects of a multimessenger detection of the mode using current technology. The detection of the imprint on neutrinos is most challenging, with a maximum detection distance of $sim!1$ kpc using the IceCube Neutrino Observatory. The maximum distance for detecting the complementary gravitational wave imprint is $sim!5$ kpc using Advanced LIGO at design sensitivity.
We study the three-dimensional (3D) hydrodynamics of the post-core-bounce phase of the collapse of a 27-solar-mass star and pay special attention to the development of the standing accretion shock instability (SASI) and neutrino-driven convection. To this end, we perform 3D general-relativistic simulations with a 3-species neutrino leakage scheme. The leakage scheme captures the essential aspects of neutrino cooling, heating, and lepton number exchange as predicted by radiation-hydrodynamics simulations. The 27-solar-mass progenitor was studied in 2D by B. Mueller et al. (ApJ 761:72, 2012), who observed strong growth of the SASI while neutrino-driven convection was suppressed. In our 3D simulations, neutrino-driven convection grows from numerical perturbations imposed by our Cartesian grid. It becomes the dominant instability and leads to large-scale non-oscillatory deformations of the shock front. These will result in strongly aspherical explosions without the need for large-scale SASI shock oscillations. Low-l-mode SASI oscillations are present in our models, but saturate at small amplitudes that decrease with increasing neutrino heating and vigor of convection. Our results, in agreement with simpler 3D Newtonian simulations, suggest that once neutrino-driven convection is started, it is likely to become the dominant instability in 3D. Whether it is the primary instability after bounce will ultimately depend on the physical seed perturbations present in the cores of massive stars. The gravitational wave signal, which we extract and analyze for the first time from 3D general-relativistic models, will serve as an observational probe of the postbounce dynamics and, in combination with neutrinos, may allow us to determine the primary hydrodynamic instability.
We present a first study of the progenitor star dependence of the three-dimensional (3D) neutrino mechanism of core-collapse supernovae. We employ full 3D general-relativistic multi-group neutrino radiation-hydrodynamics and simulate the post-bounce evolutions of progenitors with zero-age main sequence masses of $12$, $15$, $20$, $27$, and $40,M_odot$. All progenitors, with the exception of the $12,M_odot$ star, experience shock runaway by the end of their simulations. In most cases, a strongly asymmetric explosion will result. We find three qualitatively distinct evolutions that suggest a complex dependence of explosion dynamics on progenitor density structure, neutrino heating, and 3D flow. (1) Progenitors with massive cores, shallow density profiles, and high post-core-bounce accretion rates experience very strong neutrino heating and neutrino-driven turbulent convection, leading to early shock runaway. Accretion continues at a high rate, likely leading to black hole formation. (2) Intermediate progenitors experience neutrino-driven, turbulence-aided explosions triggered by the arrival of density discontinuities at the shock. These occur typically at the silicon/silicon-oxygen shell boundary. (3) Progenitors with small cores and density profiles without strong discontinuities experience shock recession and develop the 3D standing-accretion shock instability (SASI). Shock runaway ensues late, once declining accretion rate, SASI, and neutrino-driven convection create favorable conditions. These differences in explosion times and dynamics result in a non-monotonic relationship between progenitor and compact remnant mass.
The proposal that core collapse supernovae are neutrino driven is still the subject of active investigation more than fifty years after the seminal paper by Colgate and White. The modern version of this paradigm, which we owe to Wilson, proposes that the supernova shock wave is powered by neutrino heating, mediated by the absorption of electron-flavor neutrinos and antineutrinos emanating from the proto-neutron star surface, or neutrinosphere. Neutrino weak interactions with the stellar core fluid, the theory of which is still evolving, are flavor and energy dependent. The associated neutrino mean free paths extend over many orders of magnitude and are never always small relative to the stellar core radius. Thus, neutrinos are never always fluid like. Instead, a kinetic description of them in terms of distribution functions that determine the number density of neutrinos in the six-dimensional phase space of position, direction, and energy, for both neutrinos and antineutrinos of each flavor, or in terms of angular moments of these neutrino distributions that instead provide neutrino number densities in the four-dimensional phase-space subspace of position and energy, is needed. In turn, the computational challenge is twofold: (i) to map the kinetic equations governing the evolution of these distributions or moments onto discrete representations that are stable, accurate, and, perhaps most important, respect physical laws such as conservation of lepton number and energy and the Fermi--Dirac nature of neutrinos and (ii) to develop efficient, supercomputer-architecture-aware solution methods for the resultant nonlinear algebraic equations. In this review, we present the current state of the art in attempts to meet this challenge.
The next galactic core-collapse supernova (CCSN) has already exploded, and its electromagnetic (EM) waves, neutrinos, and gravitational waves (GWs) may arrive at any moment. We present an extensive study on the potential sensitivity of prospective detection scenarios for GWs from CCSNe within 5Mpc, using realistic noise at the predicted sensitivity of the Advanced LIGO and Advanced Virgo detectors for 2015, 2017, and 2019. We quantify the detectability of GWs from CCSNe within the Milky Way and Large Magellanic Cloud, for which there will be an observed neutrino burst. We also consider extreme GW emission scenarios for more distant CCSNe with an associated EM signature. We find that a three detector network at design sensitivity will be able to detect neutrino-driven CCSN explosions out to ~5.5 kpc, while rapidly rotating core collapse will be detectable out to the Large Magellanic Cloud at 50kpc. Of the phenomenological models for extreme GW emission scenarios considered in this study, such as long-lived bar-mode instabilities and disk fragmentation instabilities, all models considered will be detectable out to M31 at 0.77 Mpc, while the most extreme models will be detectable out to M82 at 3.52 Mpc and beyond.
Hydrogen-rich core collapse supernovae, known as Type II supernovae, are the most common type of stellar explosion realized in nature. They are defined by the presence of prominent hydrogen lines in their spectra. Type II supernovae are observed only in star-forming galaxies, and several events have been directly linked to massive star progenitors. Five main subclasses are identified: Type IIP (displaying a plateau in their light curve), Type IIL (displaying a light curve decline), Type IIn (displaying narrow emission lines), Type IIb (displaying increasingly strong He features with time) and 87A-likes (displaying long-rising light curves similar to that of SN 1987A). Type IIP supernovae have been robustly established as the explosions of red supergiants, while the progenitors of Type IILs remain elusive. Type IIns are likely linked to luminous blue variables, Type IIb progenitors may be interacting binary systems and the prototype of the 87A-like class was observed to be the explosion of a blue supergiant. The diversity in progenitor mass, metallicity, binarity and rotation is likely responsible for the diversity in observed explosion types, but the connection between progenitor parameters and supernova properties is not yet entirely understood theoretically nor fully mapped observationally. New observational methods for constraining this connection are currently being implemented, including the analyses of large samples of events, making use of very early data (obtained hours to days from explosion) and statistical studies of host-galaxy properties.