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.

Rotation of core-collapse supernovae (SNe) affects the neutrino emission characteristics. By comparing the neutrino properties of three three-dimensional SN simulations of a 15 M_sun progenitor (one non-rotating model and two models rotating at different velocities), we investigate how the neutrino emission varies with the flow dynamics in the SN core depending on the different degrees of rotation. The large-amplitude sinusoidal modulations due to the standing accretion-shock instability (SASI) are weaker in both the rotating models than in the non-rotating case. The SN progenitor rotation reduces the radial velocities and radial component of the kinetic energy associated with convection interior to the proto-neutron star. This effect seems to disfavor the growth of the hemispheric neutrino-emission asymmetries associated with the lepton-emission self-sustained asymmetry (LESA). An investigation of the multipole expansion of the neutrino luminosity and the electron neutrino lepton number flux shows a dominant quadrupolar mode in rotating SN models. Our findings highlight the power of using neutrinos as probes of SN hydrodynamics.

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.

Measuring the rotation of core-collapse supernovae (SN) and of their progenitor stars is extremely challenging. Here it is demonstrated that neutrinos may potentially be employed as stellar gyroscopes, if phases of activity by the standing accretion-shock instability (SASI) affect the neutrino emission prior to the onset of the SN explosion. This is shown by comparing the neutrino emission properties of self-consistent, three-dimensional (3D) SN simulations of a 15 M_sun progenitor without rotation as well as slow and fast rotation compatible with observational constraints. The explosion of the fast rotating model gives rise to long-lasting, massive polar accretion downflows with stochastic time-variability, detectable e.g. by the IceCube Neutrino Observatory for any observer direction. While spectrograms of the neutrino event rate of non-rotating SNe feature a well-known sharp peak due to SASI for observers located in the proximity of the SASI plane, the corresponding spectrograms of rotating models show activity over a wide range of frequencies, most notably above 200 Hz for rapid rotation. In addition, the Fourier power spectra of the event rate for rotating models exhibit a SASI peak with lower power than in non-rotating models. The spectra for the rotating models also show secondary peaks at higher frequencies with greater relative heights compared to the main SASI peak than for non-rotating cases. These rotational imprints will be detectable for SNe at 10 kpc or closer.

We revisit the diffuse supernova neutrino background in light of recent systematic studies of stellar core collapse that reveal the quantitative impacts of the progenitor conditions on the collapse process. In general, the dependence of the progenitor on the core-collapse neutrino emission is not monotonic in progenitor initial mass, but we show that it can, at first order, be characterized by the core compactness. For the first time, we incorporate the detailed variations in the neutrino emission over the entire mass range $8$-$100 {\rm M}_\odot$, based on (i) a long-term simulation of the core collapse of a $8.8 {\rm M}_\odot$ O-Ne-Mg core progenitor, (ii) over 100 simulations of iron core collapse to neutron stars, and (iii) half a dozen simulations of core collapse to black holes (the "failed channel"). The fraction of massive stars that undergo the failed channel remains uncertain, but in view of recent simulations which reveal high compactness to be conducive to collapse to black holes, we characterize the failed fraction by considering a threshold compactness above which massive stars collapse to black holes and below which the final remnant is a neutron star. We predict that future detections of the diffuse supernova neutrino background may have the power to reveal this threshold compactness, if its value is relatively small as suggested by interpretations of several recent astronomical observations.

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.

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.

The hard X-ray detector (HXD) onboard {\it Suzaku} measured soft $\gamma$-rays from the Type Ia supernova SN2014J at $77\pm2$ days after the explosion. Although the confidence level of the signal is about 90\% (i.e., $2 \sigma$), the $3 \sigma$ upper limit has been derived at $< 2.2 \times10^{-4}$ ph s$^{-1}$ cm$^{-2}$ in the 170 -- 250 keV band as the first independent measurement of soft $\gamma$-rays with an instrument other than {\it INTEGRAL}. For this analysis, we have examined the reproducibility of the NXB model of HXD/GSO using blank sky data. We find that the residual count rate in the 90 -- 500 keV band is distributed around an average of 0.19\% with a standard deviation of 0.42\% relative to the NXB rate. The averaged residual signals are consistent with that expected from the cosmic X-ray background. The flux of SN2014J derived from {\it Suzaku} measurements taken in one snapshot at $t=77\pm2$ days after the explosion is consistent with the {\it INTEGRAL} values averaged over the period between $t=$50 and 100 days and also with explosion models of single or double degenerate scenarios. Being sensitive to the total ejecta mass surrounding the radioactive material, the ratio between continuum and line flux in the soft gamma-ray regime might distinguish different progenitor models. The {\it Suzaku} data have been examined with this relation at $t=77\pm2$ days, but could not distinguish models between single and double degenerate-progenitors. We disfavor explosion models with larger $^{56}$Ni masses than 1 $M_\odot$, from our $1 \sigma$ error on the 170-250 keV X-ray flux of $(1.2\pm0.7) \times10^{-4}$ ph s$^{-1}$ cm$^{-2}$.

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 show that the X-ray line flux of the Mn Kalpha line at 5.9 keV from the decay of 55Fe is a promising diagnostic to distinguish between Type Ia supernova (SN Ia) explosion models. Using radiation transport calculations, we compute the line flux for two 3D explosion models: a near-Chandrasekhar mass delayed detonation and a violent merger of two white dwarfs. Both models are based on solar metallicity zero-age main sequence progenitors. Due to explosive nuclear burning at higher density, the delayed-detonation model synthesises 3.5 times more radioactive 55Fe than the merger model. As a result, we find that the peak Mn Kalpha line flux of the delayed-detonation model exceeds that of the merger model by a factor of 4.5. Since in both models the 5.9 keV X-ray flux peaks five to six years after the explosion, a single measurement of the X-ray line emission at this time can place a constraint on the explosion physics that is complementary to those derived from earlier phase optical spectra or light curves. We perform detector simulations of current and future X-ray telescopes to investigate the possibilities of detecting the X-ray line at 5.9 keV. For the delayed-detonation scenario, a line detection is feasible with Chandra up to 3 Mpc for an exposure time of 10^6 s. We find that it should be possible with currently existing X-ray instruments (with exposure times 5x10^5 s) to detect both of our models at sufficiently high S/N to distinguish between them for hypothetical events within the Local Group. The prospects for detection will be better with future missions. For example, the proposed Athena/X-IFU instrument could detect our delayed-detonation model out to a distance of 5 Mpc. This would make it possible to study future events occurring during its operational life at distances comparable to those of the recent supernovae SN 2011fe (6.4 Mpc) and SN 2014J (3.5 Mpc).

Besides the fact that the gamma-ray emission due to radioactive decays is responsible for powering the light curves of Type Ia supernovae (SNe Ia), gamma rays themselves are of particular interest as a diagnostic tool because they provide a direct way to obtain deeper insights into the nucleosynthesis and the kinematics of these explosion events. Focusing on two of the most broadly discussed SN Ia progenitor scenarios - a delayed detonation in a Chandrasekhar-mass white dwarf (WD) and a violent merger of two WDs - we use three-dimensional explosion models and perform radiative transfer simulations to obtain synthetic gamma-ray spectra. Both chosen models produce the same mass of 56Ni and have similar optical properties that are in reasonable agreement with the recently observed supernova SN 2011fe. In contrast to the optical regime, the gamma-ray emission of our two chosen models proves to be rather different. The almost direct connection of the emission of gamma rays to fundamental physical processes occuring in SNe Ia permits additional constraints concerning several explosion model properties that are not easily accessible within other wavelength ranges. Proposed future MeV missions such as GRIPS will resolve all spectral details only for nearby SNe Ia, but hardness ratio and light curve measurements still allow for a distinction of the two different models at 10 and 16 Mpc for an exposure time of 10^6 s, respectively. The possibility to detect the strongest line features up to the Virgo distance will offer the opportunity to build up a first sample of SN Ia detections in the gamma-ray energy range and underlines the importance of future space observatories for MeV gamma rays.

As the largest mass concentrations in the local Universe, nearby clusters of galaxies and their central galaxies are prime targets in searching for indirect signatures of dark matter annihilation (DMA). We seek to constrain the dark matter annihilation emission component from multi-frequency observations of the central galaxy of the Virgo cluster. The annihilation emission component is modeled by the prompt and inverse-Compton gamma rays from the hadronization of annihilation products from generic weakly interacting dark matter particles. This component is fitted to the excess of the observed data above the spectral energy distribution (SED) of the jet in M87, described with a best-fit synchrotron-self-Compton (SSC) spectrum. While this result is not sufficiently significant to claim a detection, we emphasize that a dark matter "double hump signature" can be used to unambiguously discriminate the dark matter emission component from the variable jet-related emission of M87 in future, more extended observation campaigns.

Clusters of galaxies, such as the Virgo cluster, host enormous quantities of dark matter, making them prime targets for efforts in indirect dark matter detection via potential radiative signatures from annihilation of dark matter particles and subsequent radiative losses of annihilation products. However, a careful study of ubiquitous astrophysical backgrounds is mandatory to single out potential evidence for dark matter annihilation. Here, we construct a multiwavelength spectral energy distribution for the central radio galaxy in the Virgo cluster, M87, using a state-of-the-art numerical Synchrotron Self Compton approach. Fitting recent Chandra, Fermi-LAT and Cherenkov observations, we probe different dark matter annihilation scenarios including a full treatment of the inverse Compton losses from electrons and positrons produced in the annihilation. It is shown that such a template can substantially improve upon existing dark matter detection limits.

The supernova remnant Cassiopeia A is a prime candidate for accelerating cosmic ray protons and ions. Gamma rays have been observed at GeV and TeV energies, which indicates hadronic interactions, but they could also be caused by inverse-Compton scattering of low-energy photons by accelerated electrons. We seek to predict the flux of nuclear de-excitation lines from Cas A through lower-energy cosmic rays and to compare it with COMPTEL measurements. Assuming a hadronic origin of the high-energy emission, we extrapolate the cosmic ray spectrum down to energies of 10 MeV, taking into account an equilibrium power-law momentum spectrum with a constant slope. We then calculate the nuclear line spectrum of Cassiopeia A, considering the most prominent chemical elements in the MeV band and their abundances as determined by X-ray spectroscopy. We show that the predicted line spectrum is close to the level of the COMPTEL sensitivity and agrees with conservative upper limits.

At high energies, the hadronic origin of gamma rays from supernova remnants is still debated. Assuming the observed gamma-rays from the Wolf-Rayet supernova remnant Cas A are due to accelerated protons and ions, we predict the nuclear de-excitation line emission arising from interactions with the heavy elements in the supernova ejecta. This illustrative example highlights the importance of MeV gamma ray observations of the hadronic fingerprint of cosmic rays. The lines could be observed in the MeV band with a future space mission such as GRIPS which has been proposed to ESA as an all-sky survey mission with gamma-ray, X-ray and near-infrared telescopes. In MeV gamma rays, its sensitivity will improve by a factor of 40 compared with previous missions.