In order to investigate the impact of collective neutrino oscillations (CNO) on the neutrino signal from a nearby supernova, we perform 3-flavor neutrino oscillation simulations employing the multiangle effect. The background hydrodynamic model is based on the neutrino hydrodynamic simulation of a 8.8 \Msun progenitor star. We find that CNO commences after some 100 ms post bounce. Before this, CNO is suppressed by matter-induced decoherence. In the inverted mass hierarchy, the spectrum of $\bar{\nu}_e$ becomes softer after the onset of CNO. To evaluate the detectability of this modification, we define a hardness ratio between the number of high energy neutrino events and low energy neutrino events adopting a fixed critical energy. We show that Hyper-Kamiokande (HK) can distinguish the effect of CNO for supernova distances out to $\sim 10$ kpc. On the other hand, for the normal mass hierarchy, the spectrum of $\nu_e$ becomes softer after the onset of CNO, and we show that DUNE can distinguish this feature for supernova distances out to $\sim 10$ kpc. More work is necessary to optimize the best value of critical energy for maximum sensitivity. We also show that if the spectrum of $\bar{\nu}_e$ in HK becomes softer due to CNO, the spectrum of $\nu_e$ in DUNE becomes harder, and vice versa. This synergistic observations in $\bar{\nu}_e$ and $\nu_e$, by HK and DUNE respectively, will be an intriguing opportunity to test the occurrence of CNO.

Core-collapse supernovae are dramatic explosions marking the catastrophic end of massive stars. The only means to get direct information about the supernova engine is from observations of neutrinos emitted by the forming neutron star, and through gravitational waves which are produced when the hydrodynamic flow or the neutrino flux is not perfectly spherically symmetric. The multidimensionality of the supernova engine, which breaks the sphericity of the central core such as convection, rotation, magnetic fields, and hydrodynamic instabilities of the supernova shock, is attracting great attention as the most important ingredient to understand the long-veiled explosion mechanism. Based on our recent work, we summarize properties of gravitational waves, neutrinos, and explosive nucleosynthesis obtained in a series of our multidimensional hydrodynamic simulations and discuss how the mystery of the central engines can be unraveled by deciphering these multimessengers produced under the thick veils of massive stars.

We study electron-neutrino and electron-antineutrino signals from a supernova with strong magnetic field detected by a 100 kton liquid Ar detector. The change of neutrino flavors by resonant spin-flavor conversions, matter effects, and neutrino self-interactions are taken into account. Different neutrino signals, characterized by neutronization burst event and the total event numbers of electron-neutrinos and electron-antineutrinos, are expected with different neutrino oscillation parameters and neutrino magnetic moment. Observations of supernova neutrino signals by a 100 kton liquid Ar detector would constrain oscillation parameters as well as neutrino magnetic moment in either normal and inverted mass hierarchies.

We study how the influence of the shock wave appears in neutrino oscillations and the neutrino spectrum using density profile of adiabatic explosion model of a core-collapse supernova which is calculated in an implicit Lagrangian code for general relativistic spherical hydrodynamics. We calculate expected event rates of neutrino detection at SK and SNO for various theta_{13} values and both normal and inverted hierarchies. The predicted event rates of bar{nu}_e and nu_e depend on the mixing angle theta_{13} for the inverted and normal hierarchies, respectively, and the influence of the shock appears for about 2 - 8 s when sin^2 2 theta_{13} is larger than 10^{-3}. These neutrino signals for the shock propagation is decreased by < 30 % for bar{nu}_e in inverted (SK) or by < 15 % for nu_e in normal hierarchy (SNO) compared with the case without shock. The obtained ratio of the total event for high-energy neutrinos (20 MeV < E_{nu} < 60 MeV) to low-energy neutrinos (5 MeV < E_{nu} < 20 MeV) is consistent with the previous studies in schematic semi-analytic or other hydrodynamic models of the shock propagation. The time dependence of the calculated ratio of the event rates of high-energy to low-energy neutrinos is a very useful observable which is sensitive to theta_{13} and hierarchies. Namely, time-dependent ratio shows clearer signal of the shock propagation that exhibits remarkable decrease by at most factor \sim 2 for bar{nu}_e in inverted (SK), whereas it exhibits smaller change by \sim 10 % for nu_e in normal hierarchy (SNO). Observing time-dependent high-energy to low-energy ratio of the neutrino events thus would provide a piece of very useful information to constrain theta_{13} and mass hierarchy, and eventually help understanding the propagation how the shock wave propagates inside the star.

Detailed dependence of resonant spin-flavor (RSF) conversion of supernova neutrinos on electron mole fraction Ye is investigated. Supernova explosion forms a hot-bubble and neutrino-driven wind region of which electron mole fraction exceeds 0.5 in several seconds after the core collapse. When a higher resonance of the RSF conversion is located in the innermost region, flavor change of the neutrinos strongly depends on the sign of 1-2Ye. At an adiabatic high RSF resonance the flavor conversion of bar{nu}_e -> nu_{mu,tau} occurs in Ye < 0.5 and normal mass hierarchy or in Ye > 0.5 and inverted mass hierarchy. In other cases of Ye values and mass hierarchies, the conversion of nu_e -> bar{nu}_{mu,tau} occurs. The final bar{nu}_e spectrum is evaluated in the cases of Ye < 0.5 and Ye > 0.5 taking account of the RSF conversion. Based on the obtained result, time variation of the event number ratios of low bar{nu}_e energy to high bar{nu}_e energy is discussed. In normal mass hierarchy, an enhancement of the event ratio should be seen in the period when the electron fraction in the innermost region exceeds 0.5. In inverted mass hierarchy, on the other hand, a dip of the event ratio should be observed. Therefore, the time variation of the event number ratio is useful to investigate the effect of the RSF conversion.

We investigate neutrino oscillations from core-collapse supernovae that produce magnetohydrodynamic (MHD) explosions. By calculating numerically the flavor conversion of neutrinos in the highly non-spherical envelope, we study how the explosion anisotropy has impacts on the emergent neutrino spectra through the Mikheyev-Smirnov-Wolfenstein effect. In the case of the inverted mass hierarchy with a relatively large theta_(13), we show that survival probabilities of electron type neutrinos and antineutrinos seen from the rotational axis of the MHD supernovae (i.e., polar direction), can be significantly different from those along the equatorial direction. The event numbers of electron type antineutrinos observed from the polar direction are predicted to show steepest decrease, reflecting the passage of the magneto-driven shock to the so-called high-resonance regions. Furthermore we point out that such a shock effect, depending on the original neutrino spectra, appears also for the low-resonance regions, which leads to a noticeable decrease in the electron type neutrino signals. This reflects a unique nature of the magnetic explosion featuring a very early shock-arrival to the resonance regions, which is in sharp contrast to the neutrino-driven delayed supernova models. Our results suggest that the two features in the electron type antineutrinos and neutrinos signals, if visible to the Super-Kamiokande for a Galactic supernova, could mark an observational signature of the magnetically driven explosions, presumably linked to the formation of magnetars and/or long-duration gamma-ray bursts.