Based on our recent three-dimensional core-collapse supernova (CCSN) simulations including both exploding and non-exploding models, we study the detailed neutrino signals in representative terrestrial neutrino observatories, Super-Kamiokande (Hyper-Kamiokande), DUNE, JUNO, and IceCube. We find that the physical origin of difference in the neutrino signals between 1D and 3D is mainly proto-neutron-star (PNS) convection. We study the temporal and angular variations of the neutrino signals and discuss the detectability of the time variations driven by the spiral Standing Accretion Shock Instability (spiral SASI) when it emerges for non-exploding models. In addition, we determine that there can be a large angular asymmetry in the event rate ($\gtrsim 50 \%$), but that the time-integrated signal has a relatively modest asymmetry ($\lesssim 20 \%$). Both features are associated with the lepton-number emission self-sustained asymmetry (LESA) and the spiral SASI. Moreover, our analysis suggests that there is an interesting correlation between the total neutrino energy (TONE) and the cumulative number of neutrino events in each detector, a correlation that can facilitate data analyses of real observations. We demonstrate the retrieval of neutrino energy spectra for all flavors of neutrino by applying a novel spectrum reconstruction technique to the data from multiple detectors. We find that this new method is capable of estimating the TONE within the error of $\sim$20\% if the distance to the CCSN is $\lesssim 6$ kpc.

We examine the sensitivity of neutrino emissions to stellar evolution models for a 15$M_\odot$ progenitor, paying particular attention to a phase prior to the collapse. We demonstrate that the number luminosities in both electron-type neutrinos ($\nu_e$) and their anti-partners ($\bar{\nu}_e$) differ by more than an order of magnitude by changing spatial resolutions and nuclear network sizes on stellar evolution models. We also develop a phenomenological model to capture the essential trend of the diversity, in which neutrino luminosities are expressed as a function of central density, temperature and electron fraction. In the analysis, we show that neutrino luminosity can be well characterized by these central quantities. This analysis also reveals that the most influential quantity to the time evolution of $\nu_e$ luminosity is matter density, while it is temperature for $\bar{\nu}_e$. These qualitative trends will be useful and applicable to constrain the physical state of progenitors at the final stages of stellar evolution from future neutrino observations, although more detailed systematic studies including various mass progenitors are required to assess the applicability.

We report on the core-collapse supernova simulation we conducted for a 11.2M progenitor model in three-dimensional space up to 20 ms after bounce, using a radiation hydrodynamics code with full Boltzmann neutrino transport. We solve the six-dimensional Boltzmann equations for three neutrino species and the three-dimensional compressible Euler equations with Furusawa and Togashi's nuclear equation of state. We focus on the prompt convection at ~10 ms after bounce and investigate how neutrinos are transported in the convective matter. We apply a new analysis based on the eigenvalues and eigenvectors of the Eddington tensor and make a comparison between the Boltzmann transport results and the M1 closure approximation in the transition regime between the optically thick and thin limits. We visualize the eigenvalues and eigenvectors using an ellipsoid, in which each principal axis is parallel to one of the eigenvectors and has a length proportional to the corresponding eigenvalue. This approach enables us to understand the difference between the Eddington tensor derived directly from the Boltzmann simulation and the one given by the M1 prescription from a new perspective. We find that the longest principal axis of the ellipsoid is almost always nearly parallel to the energy flux in the M1 closure approximation whereas in the Boltzmann simulation it becomes perpendicular in some transition regions, where the mean free path is ~0.1 times the radius. In three spatial dimensions, the convective motions make it difficult to predict where this happens and possibly improves the closure relation there.

Using the Boltzmann-radiation-hydrodynamics code which solves the Boltzmann equation for the neutrino transport, we present the results of the simulations with the nuclear equations of state (EOSs) of Lattimer-Swesty (LS) and Furusawa-Shen (FS). We extend the simulation time of the LS model and conducted thorough investigations, though the previous paper reported some of the results briefly. Only the LS model shows the shock revival. This seems to originate from the nuclear composition: the different nuclear composition results in the different energy loss by the photodissociation and hence the different strength of the prompt convection and the later neutrino-driven convection. The proto-neutron star seen in the FS model is more compact than that in the LS model since the existence of multi-nuclear species softens the EOS. For the behavior of neutrinos, we examined the flux and the Eddington tensor of neutrinos. In the optically thick region, the diffusion of neutrinos and the dragging by the matter motion determine the flux. In the optically thin region, the free-streaming determines it. The Eddington tensor is compared with that obtained from the M1-closure relation. The M1-closure scheme overestimates the contribution from the velocity-dependent terms in the semi-transparent region.

Recent theoretical work indicates that the neutrino radiation in core-collapse supernovae may be susceptible to flavor instabilities that set in far behind the shock, grow extremely rapidly, and have the potential to profoundly affect supernova dynamics and composition. Here we analyze the nonlinear collective oscillations that are prefigured by these instabilities. We demonstrate that a zero-crossing in $n_{\nu_e} - n_{\bar{\nu}_e}$ as a function of propagation angle is not sufficient to generate instability. Our analysis accounts for this fact and allows us to formulate complementary criteria. Using Fornax simulation data, we show that fast collective oscillations qualitatively depend on how forward-peaked the neutrino angular distributions are.

The small energy exchange via nucleon recoils in neutrino-nucleon scattering is now supposed to be one of the important factors for successful explosion of core-collapse supernovae (CCSNe) as they can change neutrino spectra through accumulation of a large number of scatterings. In finite-difference methods employed for neutrino transport in CCSN simulations, we normally can not afford to deploy a large enough number of energy bins needed to resolve this small energy exchange and sub-grid techniques are employed one way or another. In this paper we study quantitatively with the Monte Carlo (MC) method how well such a treatment performs. We first investigate the effects of nucleon recoils on the neutrino spectra and confirm that the average energy is reduced by $\sim$15% for heavy-lepton neutrinos and by much smaller quantities for other types of neutrinos in a typical post-bounce situation. It is also observed that the nucleon scattering dominates the electron scattering in the thermalization of neutrino spectra in all flavors. We then study possible artifacts that the coarse energy grid may produce in the finite-difference methods. In order to mimic the latter calculation, we re-distribute MC particles in each energy bin after a certain interval in a couple of ways and study how the results are affected and depend on the energy-resolution. We also discuss possible implications of our results for the finite-difference methods.

This paper presents the first systematic study of proto-neutron star (PNS) convection in three dimensions (3D) based on our latest numerical Fornax models of core-collapse supernova (CCSN). We confirm that PNS convection commonly occurs, and then quantify the basic physical characteristics of the convection. By virtue of the large number of long-term models, the diversity of PNS convective behavior emerges. We find that the vigor of PNS convection is not correlated with CCSN dynamics at large radii, but rather with the mass of PNS $-$ heavier masses are associated with stronger PNS convection. We find that PNS convection boosts the luminosities of $\nu_{\mu}$, $\nu_{\tau}$, $\bar{\nu}_{\mu}$, and $\bar{\nu}_{\tau}$ neutrinos, while the impact on other species is complex due to a competition of factors. Finally, we assess the consequent impact on CCSN dynamics and the potential for PNS convection to generate pulsar magnetic fields.

Neutrinos are believed to have a key role in the explosion mechanism of core-collapse supernovae as they carry most of the energy released by the gravitational collapse of a massive star. If their flavor is converted fast inside the neutrino sphere, the supernova explosion may be influenced. This paper is reporting the results of the extended work of our previous paper. We perform a thorough survey of the ELN crossing in one of our self-consistent, realistic Boltzmann simulations in two spatial dimensions under axisymmetry for the existence of the crossings between $\nu_e$ and $\bar\nu_e$ angular distributions, or the electron lepton number (ELN) crossing. We report for the first time the positive detections deep inside the core of the massive star in the vicinity of neutrino sphere at $r$ $\approx$ 16 - 21 km. We find that low values of the electron fraction $Y_e$ produced by convective motions together with the appearance of light elements are critically important to give rise to the ELN crossing by enhancing the chemical potential difference between proton and neutron, and hence by mitigating the Fermi-degeneracy of $\nu_e$. Since the region of positive detection are sustained and, in fact, expanding with time, it may have an impact on the explosion of core-collapse supernovae, observational neutrino astronomy and nucleosynthesis of heavy nuclei.

We have conducted nineteen state-of-the-art 3D core-collapse supernova simulations spanning a broad range of progenitor masses. This is the largest collection of sophisticated 3D supernova simulations ever performed. We have found that while the majority of these models explode, not all do, and that even models in the middle of the available progenitor mass range may be less explodable. This does not mean that those models for which we did not witness explosion would not explode in Nature, but that they are less prone to explosion than others. One consequence is that the "compactness" measure is not a metric for explodability. We find that lower-mass massive star progenitors likely experience lower-energy explosions, while the higher-mass massive stars likely experience higher-energy explosions. Moreover, most 3D explosions have a dominant dipole morphology, have a pinched, wasp-waist structure, and experience simultaneous accretion and explosion. We reproduce the general range of residual neutron-star masses inferred for the galactic neutron-star population. The most massive progenitor models, however, in particular vis \`a vis explosion energy, need to be continued for longer physical times to asymptote to their final states. We find that while the majority of the inner ejecta have Y$_e = 0.5$, there is a substantial proton-rich tail. This result has important implications for the nucleosynthetic yields as a function of progenitor. Finally, we find that the non-exploding models eventually evolve into compact inner configurations that experience a quasi-periodic spiral SASI mode. We otherwise see little evidence of the SASI in the exploding models.

We present a linear stability analysis of the fast-pairwise neutrino flavor conversion based on a result of our latest axisymmetric core-collapse supernova (CCSN) simulation with full Boltzmann neutrino transport. In the CCSN simulation, coherent asymmetric neutrino emissions of electron-type neutrinos ($\nu_{\rm e}$) and their anti-particles ($\bar{\nu}_{\rm e}$), in which the asymmetry of $\nu_{\rm e}$ and $\bar{\nu}_{\rm e}$ is anti-correlated with each other, occur at almost the same time as the onset of aspherical shock expansion. We find that the asymmetric neutrino emissions play a crucial role on occurrences of fast flavor conversions. The linear analysis shows that unstable modes appear in both pre- and post-shock flows; for the latter they appear only in the hemisphere of higher $\bar{\nu}_{\rm e}$ emissions (the same hemisphere with stronger shock expansion). We analyze in depth the characteristics of electron-lepton-number (ELN) crossing by closely inspecting the angular distributions of neutrinos in momentum space. The ELN crossing happens in various ways, and the property depends on the radius: in the vicinity of neutron star, $\bar{\nu}_{\rm e}$ ($\nu_{\rm e}$) dominates over $\nu_{\rm e}$ ($\bar{\nu}_{\rm e}$) in the forward (backward) direction: at the larger radius the ELN crossing occurs in the opposite way. We also find that the non-radial ELN crossing occurs at the boundary between no ELN crossing and the radial one, which is an effect of genuine multi-D transport. Our findings indicate that the collective neutrino oscillation may occur more commonly in CCSNe and suggest that the CCSN community needs to accommodate these oscillations self-consistently in the modelling of CCSNe.

We make a strong case that the fast neutrino-flavor conversion, one of the collective flavor oscillation modes, commonly occurs in core-collapse supernovae (CCSNe). It is confirmed in the numerical data obtained in realistic simulations of CCSNe but the argument is much more generic and applicable universally: the coherent neutrino-nucleus scattering makes the electron lepton number (ELN) change signs at some inward direction and trigger the flavor conversion in the outward direction in the pre-shock region. Although the ELN crossing is tiny and that is why it has eluded recognition so far, it is still large enough to induce the flavor conversion. Our findings will have an important observational consequences for CCSNe neutrinos.

Using our new state-of-the-art core-collapse supernova (CCSN) code Fornax, we explore the dependence upon spatial resolution of the outcome and character of three-dimensional (3D) supernova simulations. For the same 19-M$_{\odot}$ progenitor star, energy and radial binning, neutrino microphysics, and nuclear equation of state, changing only the number of angular bins in the $\theta$ and $\phi$ directions, we witness that our lowest resolution 3D simulation does not explode. However, when jumping progressively up in resolution by factors of two in each angular direction on our spherical-polar grid, models then explode, and explode slightly more vigorously with increasing resolution. This suggests that there can be a qualitative dependence of the outcome of 3D CCSN simulations upon spatial resolution. The critical aspect of higher spatial resolution is the adequate capturing of the physics of neutrino-driven turbulence, in particular its Reynolds stress. The greater numerical viscosity of lower-resolution simulations results in greater drag on the turbulent eddies that embody turbulent stress, and, hence, in a diminution of their vigor. Turbulent stress not only pushes the temporarily stalled shock further out, but bootstraps a concomitant increase in the deposited neutrino power. Both effects together lie at the core of the resolution dependence we observe.

We investigate the impact of asymmetric neutrino emissions on the explosive nucleosynthesis in neutrino-driven core-collapse supernovae (CCSNe). We find that the asymmetric emissions tend to yield larger amounts of proton-rich ejecta (electron fraction, $Y_e > 0.51$) in the hemisphere of the higher $\nu_{\rm e}$ emissions, meanwhile neutron-rich matter ($Y_e < 0.49$) are ejected in the opposite hemisphere of the higher ${\bar \nu}_{\rm e}$ emissions. For larger asymmetric cases with $\ge 30\%$, the neutron-rich ejecta is abundantly produced, in which there are too much elements heavier than Zn compared to the solar abundances. This may place an upper limit of the asymmetric neutrino emissions in CCSNe. The characteristic features are also observed in elemental distribution; (1) abundances lighter than Ca are insensitive to the asymmetric neutrino emissions: (2) the production of Zn and Ge is larger in the neutron-rich ejecta even for smaller asymmetric cases with $\le 10\%$. We discuss these observational consequences, which may account for the (anti-)correlations among asymmetries of heavy elements and neutron star kicks in supernova remnants (SNRs). Future SNR observations of the direct measurement for the mass and spatial distributions of $\alpha$ elements, Fe, Zn and Ge will provide us the information on the asymmetric degree of neutrino emissions.

In this {\it Letter}, we present the result of an axisymmetric CCSN simulation conducted with appropriate treatments of neutrino transport and proper motions of proto-neutron star (PNS), in which a remarkable PNS acceleration is observed in association with asymmetric neutrino emissions by the time of $300$ms after bounce. We find that these asymmetric neutrino emissions play important roles in the acceleration of PNS in this phase. The correlation between the PNS proper motion and the asymmetric ejecta is similar to that in the neutron star (NS) kick of hydrodynamic origin. Both electron-type neutrinos ($\nu_{\rm e}$) and their anti-particles ($\bar{\nu}_{\rm e}$) have $\sim 10\%$ level of asymmetry between the northern and southern hemispheres, while other heavy-leptonic neutrinos ($\nu_x$) have much smaller asymmetry of $\sim 1 \%$. The emissions of $\bar{\nu}_{\rm e}$ and $\nu_x$ are higher in the hemisphere of stronger shock expansion whereas the $\nu_{\rm e}$ emission is enhanced in the opposite hemisphere: in total the neutrinos carry some linear momentum to the hemisphere of the stronger shock expansion. The asymmetry is attributed to the non-spherical distribution of electron-fraction ($Y_e$) in the envelope of PNS. Although it is similar to LESA (lepton-emission self-sustained asymmetry), the $Y_e$ asymmetry seems to be associated with the PNS motion: the latter triggers lateral circular motions in the envelope of PNS by breaking symmetry of the matter distribution there, which are then sustained by a combination of convection, lateral neutrino diffusion and matter-pressure gradient. Our findings may have an influence on the current theories on the NS kick mechanism although long-term simulations are required to assess their impact on the later evolution.

We present a new method for neutrino-matter coupling in multi-dimensional radiation-hydrodynamic simulations of core-collapse supernova (CCSN) with the full Boltzmann neutrino transport. The development is motivated by the fact that the accurate conservation of momentum is required for reliable numerical modelings of CCSN dynamics including a recoil of proto-neutron star (PNS). The new method is built on a hybrid approach, in which we use the energy-momentum tensor of neutrinos to compute the momentum feedback from neutrino to matter in the optically thick region while we employ the collision term in the optically thin region. In this method we utilize a general relativistic description of radiation-hydrodynamics with angular moments, which allows us to evaluate the momentum feedback from neutrino to matter without inconsistency with our Boltzmann solver. We demonstrate that the new method substantially improves the accuracy of linear momentum conservation in our CCSN simulations under reasonable angular resolutions in momentum space, alleviating the difficulty in giving the diffusion limit precisely with the discrete ordinate (Sn) method. It is the first-ever demonstration that the PNS kick can be handled directly and properly in multi-dimensional radiation-hydrodynamic simulations with the full Boltzmann neutrino transport.

We study the gravitational wave signal from eight new 3D core-collapse supernova simulations. We show that the signal is dominated by $f$- and $g$-mode oscillations of the protoneutron star and its frequency evolution encodes the contraction rate of the latter, which, in turn, is known to depend on the star's mass, on the equation of state, and on transport properties in warm nuclear matter. A lower-frequency component of the signal, associated with the standing accretion shock instability, is found in only one of our models. Finally, we show that the energy radiated in gravitational waves is proportional to the amount of turbulent energy accreted by the protoneutron star.

Neutrinos are densely populated deep inside the core of massive stars after their gravitational collapse to produce supernova explosions and form compact stars such as neutron stars (NS) and black holes (BH). It has been considered that they may change their flavor identities through so-called fast-pairwise conversions induced by mutual forward scatterings. If that is really the case, the dynamics of supernova explosion will be influenced, since the conversion may occur near the neutrino sphere, from which neutrinos are effectively emitted. In this paper, we conduct a pilot study of such possibilities based on the results of fully self-consistent, realistic simulations of a core-collapse supernova explosion in two spatial dimensions under axisymmetry. As we solved the Boltzmann equations for neutrino transfer in the simulation not as a post-process but in real time, the angular distributions of neutrinos in momentum space for all points in the core at all times are available, a distinct feature of our simulations. We employ some of these distributions extracted at a few selected points and times from the numerical data and apply linear analysis to assess the possibility of the conversion. We focus on the vicinity of the neutrino sphere, where different species of neutrinos move in different directions and have different angular distributions as a result. This is a pilot study for a more thorough survey that will follow soon. We find no positive sign of conversion unfortunately at least for the spatial points and times we studied in this particular model. We hence investigate rather in detail the condition for the conversion by modifying the neutrino distributions rather arbitrarily by hand.

With the Boltzmann-radiation-hydrodynamics code, which we have developed to solve numerically the Boltzmann equations for neutrino transfer, the Newtonian hydrodynamics equations, and the Newtonian self-gravity simultaneously and consistently, we simulate the collapse of a rotating core of the progenitor with a zero-age-main-sequence mass of $11.2\,M_\odot$ and a shelluler rotation of $1\,{\rm rad\,s^{-1}}$ at the center. We pay particular attention in this paper to the neutrino distribution in phase space, which is affected by the rotation. By solving the Boltzmann equations directly, we can assess the rotation-induced distortion of the angular distribution in momentum space, which gives rise to the rotational component of the neutrino flux. We compare the Eddington tensors calculated both from the raw data and from the M1-closure approximation. We demonstrate that the Eddington tensor is determined by complicated interplays of the fluid velocity and the neutrino interactions and that the M1-closure, which assumes that the Eddington factor is determined by the flux factor, fails to fully capture this aspect, especially in the vicinity of the shock. We find that the error in the Eddington factor reaches $\sim 20\%$ in our simulation. This is due not to the resolution but to the different dependence of the Eddington and flux factors on the angular profile of the neutrino distribution function, and hence modification to the closure relation is needed.

We perform an extensive study of the influence of nuclear weak interactions on core-collapse supernovae (CCSNe), paying particular attention to consistency between nuclear abundances in the equation of state (EOS) and nuclear weak interactions. We compute properties of uniform matter based on the variational method. For inhomogeneous nuclear matter, we take a full ensemble of nuclei into account with various finite-density and thermal effects and directly use the nuclear abundances to compute nuclear weak interaction rates. To quantify the impact of a consistent treatment of nuclear abundances on CCSN dynamics, we carry out spherically symmetric CCSN simulations with full Boltzmann neutrino transport, systematically changing the treatment of weak interactions, EOSs, and progenitor models. We find that the inconsistent treatment of nuclear abundances between the EOS and weak interaction rates weakens the EOS dependence of both the dynamics and neutrino signals. We also test the validity of two artificial prescriptions for weak interactions of light nuclei and find that both prescriptions affect the dynamics. Furthermore, there are differences in neutrino luminosities by ~10% and in average neutrino energies by 0.25-1 MeV from those of the fiducial model. We also find that the neutronization burst neutrino signal depends on the progenitor more strongly than on the EOS, preventing a detection of this signal from constraining the EOS.

The progenitor stars of core-collapse supernovae (CCSNe) are asymmetrically fluctuating due to turbulent convections in the late stages of their lives. The progenitor asymmetry at the pre-supernova stage has recently caught the attention as a new ingredient to facilitate shock revival in the delayed neutrino-heating mechanism. In this paper, we investigate the importance of the progenitor asymmetries to shock revival with a semi-analytical approach. Free parameters were chosen such that the time evolution of shock radii and mass accretion rates are compatible with the results of detailed numerical simulations of CCSNe in spherical symmetry. We first estimate the amplitude of asymmetries required for the shock revival by the impulsive change of pre-shock flows in the context of neutrino heating mechanism, and then convert the amplitude to the corresponding amplitude in the pre-supernova phase by taking into account the growth of asymmetries during infall. We apply our model to various types of progenitors and find that the requisite amplitude of pre-supernova asymmetry is roughly three times larger than the prediction by current stellar evolution models unless other additional physical ingredients such as multi-dimensional fluid instabilities and turbulent convections in post-shock flows aid shock revival. We thus conclude that progenitor asymmetries can not trigger the shock revival by the impulsive way but rather play a supplementary role in reality.

We propose a new class of method for solving nonlinear systems of equations, which, among other things,has four nice features: (i) it is inspired by the mathematical property of damped oscillators, (ii) it can be regarded as a simple extention to the Newton-Raphson(NR) method, (iii) it has the same local convergence as the NR method does, (iv) it has a significantly wider convergence region or the global convergence than that of the NR method. In this article, we present the evidence of these properties, applying our new method to some examples and comparing it with the NR method.

We present the first results of our spatially axisymmetric core-collapse supernova simulations with full Boltzmann neutrino transport, which amount to a time-dependent 5-dimensional (2 in space and 3 in momentum space) problem in fact. Special relativistic effects are fully taken into account with a two-energy-grid technique. We performed two simulations for a progenitor of 11.2M, employing different nuclear equations-of-state (EOS's): Lattimer and Swesty's EOS with the incompressibility of K = 220MeV (LS EOS) and Furusawa's EOS based on the relativistic mean field theory with the TM1 parameter set (FS EOS). In the LS EOS the shock wave reaches ~700km at 300ms after bounce and is still expanding whereas in the FS EOS it stalled at ~200km and has started to recede by the same time. This seems to be due to more vigorous turbulent motions in the former during the entire post-bounce phase, which leads to higher neutrino-heating efficiency in the neutrino-driven convection. We also look into the neutrino distributions in momentum space, which is the advantage of the Boltzmann transport over other approximate methods. We find non-axisymmetric angular distributions with respect to the local radial direction, which also generate off-diagonal components of the Eddington tensor. We find that the r {\theta}-component reaches ~10% of the dominant rr-component and, more importantly, it dictates the evolution of lateral neutrino fluxes, dominating over the {\theta}{\theta}-component, in the semi-transparent region. These data will be useful to further test and possibly improve the prescriptions used in the approximate methods.

This paper is a sequel to our previous one (Kato et al.2015), which calculated the luminosities and spectra of electron-type anti-neutrinos ($\bar{\nu}_e$'s) from the progenitors of core-collapse supernovae. Expecting that a capability to detect electron-type neutrinos ($\nu_e$'s) will increase dramatically with the emergence of liquid-argon detectors such as DUNE, we broaden the scope in this study to include all-flavors of neutrinos emitted from the pre-bounce phase. We pick up three progenitor models of an electron capture supernova (ECSN) and iron-core collapse supernovae (FeCCSNe). We find that the number luminosities reach $\sim10^{57} \mathrm{s^{-1}}$ and $\sim10^{53} \mathrm{s^{-1}}$ at maximum for $\nu_e$ and $\bar{\nu}_e$, respectively. We also estimate the numbers of detection events at terrestrial neutrino detectors including DUNE, taking flavor oscillations into account and assuming the distance to the progenitors to be 200 pc. It is demonstrated that $\bar{\nu}_e$'s from the ECSN-progenitor will be undetected at almost all detectors, whereas we will be able to observe $\gtrsim$15900 $\nu_e$'s at DUNE for the inverted mass hierarchy. From the FeCCSN-progenitors, the number of $\bar{\nu}_e$ events will be largest for JUNO, 200-900 $\bar{\nu}_e$'s, depending on the mass hierarchy whereas the number of $\nu_e$ events at DUNE is $\gtrsim$2100 for the inverted mass hierarchy. These results imply that the detection of $\bar{\nu}_e$'s is useful to distinguish FeCCSN- from ECSN-progenitors, while $\nu_e$'s will provide us with detailed information on the collapse phase regardless of the type and mass of progenitor.

The mechanism driving core-collapse supernovae is sensitive to the interplay between matter and neutrino radiation. However, neutrino radiation transport is very difficult to simulate, and several radiation transport methods of varying levels of approximation are available. We carefully compare for the first time in multiple spatial dimensions the discrete ordinates (DO) code of Nagakura, Yamada, and Sumiyoshi and the Monte Carlo (MC) code Sedonu, under the assumptions of a static fluid background, flat spacetime, elastic scattering, and full special relativity. We find remarkably good agreement in all spectral, angular, and fluid interaction quantities, lending confidence to both methods. The DO method excels in determining the heating and cooling rates in the optically thick region. The MC method predicts sharper angular features due to the effectively infinite angular resolution, but struggles to drive down noise in quantities where subtractive cancellation is prevalent, such as the net gain in the protoneutron star and off-diagonal components of the Eddington tensor. We also find that errors in the angular moments of the distribution functions induced by neglecting velocity dependence are sub-dominant to those from limited momentum-space resolution. We briefly compare directly computed second angular moments to those predicted by popular algebraic two-moment closures, and find that the errors from the approximate closures are comparable to the difference between the DO and MC methods. Included in this work is an improved Sedonu code, which now implements a fully special relativistic, time-independent version of the grid-agnostic Monte Carlo random walk approximation.

We developed a three-dimensional radiative transfer code for an ultra-relativistic background flow-field by using the Monte Carlo (MC) method in the context of gamma-ray burst (GRB) emission. For obtaining reliable simulation results in the coupled computation of MC radiation transport with relativistic hydrodynamics which can reproduce GRB emission, we validated radiative transfer computation in the ultra-relativistic regime and assessed the appropriate simulation conditions. The radiative transfer code was validated through two test calculations: (1) computing in different inertial frames and (2) computing in flow-fields with discontinuous and smeared shock fronts. The simulation results of the angular distribution and spectrum were compared among three different inertial frames and in good agreement with each other. If the time duration for updating the flow-field was sufficiently small to resolve a mean free path of a photon into ten steps, the results were thoroughly converged. The spectrum computed in the flow-field with a discontinuous shock front obeyed a power-law in frequency whose index was positive in the range from 1 to 10 MeV. The number of photons in the high-energy side decreased with the smeared shock front because the photons were less scattered immediately behind the shock wave due to the small electron number density. The large optical depth near the shock front was needed for obtaining high-energy photons through bulk Compton scattering. Even one-dimensional structure of the shock wave could affect the results of radiation transport computation.

We have constructed a nuclear equation of state (EOS) that includes a full nuclear ensemble for use in core-collapse supernova simulations. It is based on the EOS for uniform nuclear matter that two of the authors derived recently, applying a variational method to realistic two- and there-body nuclear forces. We utilize an extended liquid drop model of heavy nuclei and a quantum-theoretical mass evaluation for light nuclei. In addition to realistic nuclear forces, the inclusion of in-medium effects on the full ensemble of nuclei makes the new EOS one of the most realistic EOS's for supernova simulations. We make comparisons with the FYSS EOS, which is based on the same formulation for the nuclear ensemble but adopts the relativistic mean field (RMF) theory with the TM1 parameter set for uniform nuclear matter. The new EOS is softer than the FYSS EOS around and above nuclear saturation densities. We find that neutron-rich nuclei with small mass numbers are more abundant in the new EOS than in the FYSS EOS because of the larger saturation densities and smaller symmetry energy of nuclei in the former. We apply the two EOS's to 1D supernova simulations and find that the new EOS gives lower electron fractions and higher temperatures in the collapse phase owing to the smaller symmetry energy. As a result, the inner core has smaller masses for the new EOS. It is more compact, on the other hand, due to the softness of the new EOS and bounces at higher densities. The ensuing outward propagations of the shock wave in the outer core are very similar in the two simulations, which may be an artifact, though, caused by the use of the same tabulated electron capture rates for heavy nuclei ignoring differences in the nuclear composition between the two EOS's in these computations.

We investigate the criterion for the acoustic mechanism to work successfully in core-collapse supernovae. The acoustic mechanism is an alternative to the neutrino-heating mechanism. It was proposed by Burrows et al., who claimed that acoustic waves emitted by $g$-mode oscillations in proto-neutron stars (PNS) energize a stalled shock wave and eventually induce an explosion. Previous works mainly studied to which extent the $g$-modes are excited in the PNS. In this paper, on the other hand, we investigate how strong the acoustic wave needs to be if it were to revive a stalled shock wave. By adding the acoustic power as a new axis, we draw a critical surface, an extension of the critical curve commonly employed in the context of neutrino heating. We perform both 1D and 2D parametrized simulations, in which we inject acoustic waves from the inner boundary. In order to quantify the power of acoustic waves, we use the extended Myers theory to take neutrino reactions into proper account. We find for the 1D simulations that rather large acoustic powers are required to relaunch the shock wave, since the additional heating provided by the secondary shocks developed from acoustic waves is partially canceled by the neutrino cooling that is also enhanced. In 2D, the required acoustic powers are consistent with those of Burrows et al. Our results seem to imply, however, that it is the sum of neutrino heating and acoustic powers that matters for shock revival.

We investigate the influences of the nuclear composition on the weak interaction rates of heavy nuclei during the core collapse of massive stars. The nuclear abundances in nuclear statistical equilibrium (NSE) are calculated by some equation of state (EOS) models including in-medium effects on nuclear masses. We systematically examine the sensitivities of electron capture and neutrino-nucleus scattering on heavy nuclei to the nuclear shell effects and the single nucleus approximation. We find that the washout of shell effects at high temperatures brings significant change to weak rates by smoothing the nuclear abundance distribution: the electron capture rate decreases by $\sim$20$\%$ in the early phase and increases by $\sim$40$\%$ in the late phase at most, while the cross section for neutrino-nucleus scattering is reduced by $\sim$15$\%$. This is because the open-shell nuclei become abundant instead of those with closed neutron shells as the shell effects disappear. We also find that the single-nucleus description based on the average values leads to underestimations of weak rates. Electron captures and neutrino coherent scattering on heavy nuclei are reduced by $\sim$80$\%$ in the early phase and by $\sim$5$\%$ in the late phase, respectively. These results indicate that NSE like EOS accounting for shell washout is indispensable for the reliable estimation of weak interaction rates in simulations of core-collapse supernovae.

We present a newly developed moving-mesh technique for the multi-dimensional Boltzmann-Hydro code for the simulation of core-collapse supernovae (CCSNe). What makes this technique different from others is the fact that it treats not only hydrodynamics but also neutrino transfer in the language of the 3+1 formalism of general relativity (GR), making use of the shift vector to specify the time evolution of the coordinate system. This means that the transport part of our code is essentially general relativistic although in this paper it is applied only to the moving curvilinear coordinates in the flat Minknowski spacetime, since the gravity part is still Newtonian. The numerical aspect of the implementation is also described in detail. Employing the axisymmetric two-dimensional version of the code, we conduct two test computations: oscillations and runaways of proto-neutron star (PNS). We show that our new method works fine, tracking the motions of PNS correctly. We believe that this is a major advancement toward the realistic simulation of CCSNe.

A new computationally efficient method has been introduced to treat self-gravity in mesh based hydrodynamical simulations. It is applied simply by slightly modifying the Poisson equation into an inhomogeneous wave equation. This roughly corresponds to the weak field limit of the Einstein equations in general relativity, and as long as the gravitation propagation speed is taken to be larger than the hydrodynamical characteristic speed, the results agree with solutions for the Poisson equation. The solutions almost perfectly agree if the domain is taken large enough, or appropriate boundary conditions are given. Our new method can not only significantly reduce the computational time compared with existent methods, but is also fully compatible with massive parallel computation, nested grids and adaptive mesh refinement techniques, all of which can accelerate the progress in computational astrophysics and cosmology.

We examine the bright radio synchrotron counterparts of low-luminosity gamma-ray bursts (llGRBs) and relativistic supernovae (SNe) and find that they can be powered by spherical hypernova (HN) explosions. Our results imply that radio-bright HNe are driven by relativistic jets that are choked deep inside the progenitor stars or quasi-spherical magnetized winds from fast-rotating magnetars. We also consider the optical synchrotron counterparts of radio-bright HNe and show that they can be observed as precursors several days before the SN peak with an r-band absolute magnitude of M_r ~ -14 mag. While previous studies suggested that additional trans-relativistic components are required to power the bright radio emission, we find that they overestimated the energy budget of the trans-relativistic component by overlooking some factors related to the minimum energy of non-thermal electrons. If an additional trans-relativistic component exists, then a much brighter optical precursor with M_r ~ -20 mag can be expected. Thus, the scenarios of radio-bright HNe can be distinguished by using optical precursors, which can be detectable from < 100 Mpc by current SN surveys like the Kiso SN Survey, Palomar Transient Factory, and Panoramic Survey Telescope & Rapid Response System.

We propose a novel numerical method for solving multi-dimensional, special relativistic Boltzmann equations for neutrinos coupled to hydrodynamics equations. It is meant to be applied to simulations of core-collapse supernovae. We handle special relativity in a non-conventional way, taking account of all orders of v/c. Consistent treatment of advection and collision terms in the Boltzmann equations is the source of difficulties, which we overcome by employing two different energy grids: Lagrangian remapped and laboratory fixed grids. We conduct a series of basic tests and perform a one-dimensional simulation of core-collapse, bounce and shock-stall for a 15M_{sun} progenitor model with a minimum but essential set of microphysics. We demonstrate in the latter simulation that our new code is capable of handling all phases in core-collapse supernova. For comparison, a non-relativistic simulation is also conducted with the same code, and we show that they produce qualitatively wrong results in neutrino transfer. Finally, we discuss a possible incorporation of general relativistic effects in our method.

The effect of rotation on the explosion of core-collapse supernovae is investigated systematically in three-dimensional simulations. In order to obtain the critical conditions for explosion as a function of mass accretion rate, neutrino luminosity, and specific angular momentum, rigidly rotating matter was injected from the outer boundary with an angular momentum, which is increased every 500 ms. It is found that there is a critical value of the specific angular momentum, above which the standing shock wave revives, for a given combination of mass accretion rate and neutrino luminosity, i.e. an explosion can occur by rotation even if the neutrino luminosity is lower than the critical value for a given mass accretion rate in non-rotational models. The coupling of rotation and hydrodynamical instabilities plays an important role to characterize the dynamics of shock revival for the range of specific angular momentum that are supposed to be realistic. Contrary to expectations from past studies, the most rapidly expanding direction of the shock wave is not aligned with the rotation axis. Being perpendicular to the rotation axis on average, it can be oriented in various directions. Its dispersion is small when the spiral mode of the standing accretion shock instability (SASI) governs the dynamics, while it is large when neutrino-driven convection is dominant. As a result of the comparison between 2D and 3D rotational models, it is found that m=!0 modes of neutrino-driven convection or SASI are important for shock revival around the critical surface.

In this study, we conduct three-dimensional hydrodynamic simulations systematically to investigate the flow patterns behind the accretion shock waves that are commonly formed in the post-bounce phase of core-collapse supernovae. Adding small perturbations to spherically symmetric, steady, shocked accretion flows, we compute the subsequent evolutions to find what flow pattern emerges as a consequence of hydrodynamical instabilities such as convection and standing accretion shock instability (SASI) for different neutrino luminosities and mass accretion rates. Depending on these two controlling parameters, various flow patterns are indeed realized. We classify them into three basic patterns and two intermediate ones; the former includes sloshing motion (SL), spiral motion (SP) and multiple buoyant bubble formation (BB); the latter consists of spiral motion with buoyant-bubble formation (SPB) and spiral motion with pulsationally changing rotational velocities (SPP). Although the post-shock flow is highly chaotic, there is a clear trend in the pattern realization. The sloshing and spiral motions tend to be dominant for high accretion rates and low neutrino luminosities, and multiple buoyant bubbles prevail for low accretion rates and high neutrino luminosities. It is interesting that the dominant pattern is not always identical between the semi-nonlinear and nonlinear phases near the critical luminosity; the intermediate cases are realized in the latter case. Running several simulations with different random perturbations, we confirm that the realization of flow pattern is robust in most cases.

The observations of jet breaks in the afterglows of short gamma-ray bursts (SGRBs) indicate that the jet has a small opening angle of < 10{\deg}. The collimation mechanism of the jet is a longstanding theoretical problem. We numerically analyze the jet propagation in the material ejected by double neutron star merger, and demonstrate that if the ejecta mass is > 10^{-2} M_{sun}, the jet is well confined by the cocoon and emerges from the ejecta with the required collimation angle. Our results also suggest that there are some populations of choked (failed) SGRBs or low-luminous new types of event. By constructing a model for SGRB 130603B, which is associated with the first kilonova/macronova can- didate, we infer that the equation-of-state of neutron stars would be soft enough to provide sufficient ejecta to collimate the jet, if this event was associated with a double neutron star merger.

We perform numerical experiments to investigate the influence of inelastic neutrino reactions with light nuclei on the standing accretion shock instability (SASI). The time evolution of shock waves is calculated with a simple light-bulb approximation for the neutrino transport and a multi-nuclei equation of state. The neutrino absorptions and inelastic interactions with deuterons, tritons, helions and alpha particles are taken into account in the hydrodynamical simulations in addition to the ordinary charged-current interactions with nucleons. Axial symmetry is assumed but no equatorial symmetry is imposed. We show that the heating rates of deuterons reach as high as 10% of those of nucleons around the bottom of the gain region. On the other hand, alpha particles are heated near the shock wave, which is important when the shock wave expands and the density and temperature of matter become low. It is also found that the models with heating by light nuclei have different evolutions from those without it in the non-linear phase of SASI. This results is because matter in the gain region has a varying density and temperature and there appear sub-regions that are locally rich in deuterons and alpha particles. Although the light nuclei are never dominant heating sources and they work favorably for shock revival in some cases and unfavorably in other cases, they are non-negligible and warrant further investigation.

We perform some experimental simulations in spherical symmetry and axisymmetry to understand the post-shock-revival evolution of core-collapse supernovae. Assuming that the stalled shock wave is relaunched by neutrino heating and employing the so-called light bulb approximation, we induce shock revival by raising the neutrino luminosity by hand up to the critical value, which is also de- termined by dynamical simulations. A 15M_{sun} progenitor model is employed. We incorporate nuclear network calculations with a consistent equation of state in the simulations to account for the energy release by nuclear reactions and their feedback to hydrodynamics. Varying the shock-relaunch time rather arbitrarily, we investigate the ensuing long-term evolutions systematically, paying particular attention to the explosion energy and nucleosynthetic yields as a function of this relaunch time, or equivalently the accretion rate at shock revival. We study in detail how the diagnostic explosion energy approaches the asymptotic value and which physical processes contribute to the explosion energy in what proportions as well as their dependence on the relaunch time and the dimension of dynamics. We find that the contribution of nuclear reactions to the explosion energy is comparable to or greater than that of neutrino heating. In particular, recombinations are dominant over burnings in the contributions of nuclear reactions. Interestingly 1D models studied in this paper cannot produce the appropriate explosion energy and nickel mass simultaneously, overproducing nickels, whereas this problem is resolved in 2D models if the shock is relaunched at 300-400ms after bounce.

We develop a new semi-dynamical method to study shock revival by neutrino heating in core- collapse supernovae. Our new approach is an extension of the previous studies that employ spherically symmetric, steady, shocked accretion flows together with the light bulb approximation. The latter has been widely used in the supernova community for the phenomenological investigation of the criteria for successful supernova explosions. In the present approach, on the other hand, we get rid of the steady-state condition and take into account shock wave motions instead. We have in mind the scenario that not the critical luminosity but the critical fluctuation generated by hydrodynamical instabilities such as SASI and neutrino-driven convection in the post-shock region determines the onset of shock revival. After confirming that the new approach indeed captures the dynamics of revived shock wave qualitatively, we then apply the method to various initial conditions and find that there is a critical fluctuation for shock revival, which can be well fit by the following formula: f_crit ~ 0.8 * (M_in/1.4M_sun) * {1- (rsh/10^8cm)}, in which fcrit denotes the critical pressure fluctuation normalized by the unperturbed post-shock value. Min and rsh stand for the mass of the central compact object and the shock radius, respectively. The critical fluctuation decreases with the shock radius, whereas it increases with the mass of the central object. We discuss the possible implications of our results for 3D effects on shock revival, which is currently controversial in the supernova community.

We numerically investigate the jet propagation through a rotating collapsing Wolf-Rayet star with detailed central engine physics constructed based on the neutrino-driven collapsar model. The collapsing star determines the evolution of mass accretion rate, black hole mass and spin, all of which are important ingredients for determining the jet luminosity. We reveal that neutrino-driven jets in rapidly spinning Wolf-Rayet stars are capable of breaking out from the stellar envelope, while those propagating in slower rotating progenitors fail to jet breakout due to insufficient kinetic power. For progenitor models with successful jet breakouts, the kinetic energy accumulated in the cocoon could be as large as ~1051erg and might significantly contribute to the luminosity of the afterglow emission or to the kinetic energy of the accompanying supernova if nickel production takes place. We further analyze the post breakout phase using a simple analytical prescription and conclude that the relativistic jet component could produce events with an isotropic-luminosity Lp(iso) ~1052erg/s and isotropic-energy Ej(iso) ~1054erg. Our findings support the idea of rapidly rotating Wolf-Rayet stars as plausible progenitors of GRBs, while slowly rotational ones could be responsible for low luminosity or failed GRBs.

This is a status report on our endeavor to reveal the mechanism of core-collapse supernovae (CCSNe) by large-scale numerical simulations. Multi-dimensionality of the supernova engine, general relativistic magnetohydrodynamics, energy and lepton number transport by neutrinos emitted from the forming neutron star as well as nuclear interactions there, are all believed to play crucial roles in repelling infalling matter and producing energetic explosions. These ingredients are nonlinearly coupled with one another in the dynamics of core-collapse, bounce, and shock expansion. Serious quantitative studies of CCSNe hence make extensive numerical computations mandatory. Since neutrinos are neither in thermal nor in chemical equilibrium in general, their distributions in the phase space should be computed. This is a six dimensional (6D) neutrino transport problem and quite a challenge even for those with an access to the most advanced numerical resources such as the "K computer". To tackle this problem, we have embarked on multi-front efforts. In particular we report in this paper our recent progresses in the treatments of multi-dimensional (multi-D) radiation-hydrodynamics. We are currently proceeding on two different paths to the ultimate goal; in one approach we employ an approximate but highly efficient scheme for neutrino transport and treat 3D hydrodynamics and/or general relativity rigorously; some neutrino-driven explosions will be presented and comparisons will be made between 2D and 3D models quantitatively; in the second approach, on the other hand, exact but so far Newtonian Boltzmann equations are solved in two and three spatial dimensions; we will show some demonstrative test simulations. We will also address the perspectives of exa-scale computations on the next generation supercomputers.

We investigate the propagation of accretion-powered jets in various types of massive stars such as Wolf-Rayet stars, light Population III (Pop III) stars, and massive Pop III stars, all of which are the progenitor candidates of Gamma-Ray Bursts (GRBs). We perform two dimensional axisymmetric simulations of relativistic hydrodynamics taking into account both the envelope collapse and the jet propagation (i.e., the negative feedback of the jet on the accretion). Based on our hydrodynamic simulations, we show for the first time that the accretion-powered jet can potentially break out relativistically from the outer layers of Pop III progenitors. In our simulations, the accretion rate is estimated by the mass flux going through the inner boundary, and the jet is injected with a fixed accretion-to-jet conversion efficiency $\eta$. By varying the efficiency $\eta$ and opening angle $\theta_{op}$ for more than 40 models, we find that the jet can make a relativistic breakout from all types of progenitors for GRBs if a simple condition $\eta \gtrsim 10^{-4} (\theta_{op}/8^{\circ})^2$ is satisfied, which is consistent with analytical estimates. Otherwise no explosion or some failed spherical explosions occur.

We investigate by two-dimensional axisymmetric relativistic hydrodynamical simulations (1) jet propagations through an envelope of a rapidly rotating and collapsing massive star, which is supposed to be a progenitor of long duration gamma ray bursts (GRBs), (2) breakouts and subsequent expansions into stellar winds and (3) accompanying photospheric emissions. We find that if the envelope rotates uniformly almost at the mass shedding limit, its outer part stops contracting eventually when the centrifugal force becomes large enough. Then another shock wave is formed, propagates outwards and breaks out of the envelope into the stellar wind. Which breaks out earlier, the jet or the centrifugal bounce-induced shock, depends on the timing of jet injection. If the shock breakout occurs earlier owing to a later injection, the jet propagation and subsequent photospheric emissions are affected substantially. We pay particular attention to observational consequences of the difference in the timing of jet injection. We calculate optical depths to find the location of photospheres, extracting densities and temperatures at appropriate retarded times from the hydrodynamical data. We show that the luminosity and observed temperature of the photospheric emissions are both much lower than those reported in previous studies. Although luminosities are still high enough for GRBs, the observed temperature are lower than the energy at the spectral peak expected by the Yonetoku-relation. This may imply that energy exchanges between photons and matter are terminated deeper inside or some non-thermal processes are operating to boost photon energies.

We present a new formula to numerically construct configurations in rotational equilibrium, which consist of multiple layers. Each layer rotates uniformly or differentially according to cylindrical rotation-laws that are different from layer to layer. Assuming a different barotropic equation of state (EOS) for each layer, we solve the Bernoulli equation in each layer separately and combine the solutions by imposing continuity of the pressure at each boundary of the layers. It is confirmed that a single continuous barotropic EOS is incompatible with the junction condition. Identifying appropriate variables to be solved, we construct a convergent iteration scheme. For demonstration, we obtain two-layered configurations, each layer of which rotates rapidly with either an "$\Omega$-constant law" or a "$j$-constant law" or a "$v$-constant law". Other rotation laws and/or a larger number of layers can be treated similarly. We hope that this formula will be useful in studying the stellar evolution in multi-dimension with the non-spherical configuration induced by rotation being fully taken into account.

Dissipations, e.g. heat flow and bulk and shear viscosities, cause the transport of energy, momentum and angular momentum, which is the essence of accretion of matters onto celestial objects. Dissipations are usually described by the Fourier and Navier-Stokes laws ("classic laws" of dissipations). However the classic laws result in an infinitely fast propagation of dissipations. In relativistic formulation, the classic laws of dissipations violate the causality, and hence no relativistic theory of accretion flow onto celestial object is formulated. In this short report, we summarize the causal dissipative hydrodynamics, so-called "Extended Irreversible Thermodynamics" (EIT), with a supplemental comment of which the original works of EIT are not aware, and then show two theorems about relativistic dissipative flows around a Schwarzschild black hole. By these theorems, a significant property of EIT in contrast with classic laws of dissipations is clarified, and a dissipative instability of an exact solution of relativistic perfect fluid flow is also obtained.

We show for the first time the direct time variable radio images in the context of shocked accretion flows around a black hole under the general relativistic treatment of both hydrodynamics and radiation transfer. Time variability around a black hole can be induced by the non-axisymmetric standing accretion shock instability (namely Black Hole SASI). Since the spiral arm shock waves generate the density and temperature waves at the post shock region, they cause time variability in the black hole vicinity. Based on our dynamical simulations, we discuss a possibility of detection for the time variable radio images of M87 by the future space telescope VSOP2/ASTRO-G satellite. The most luminous part of the images is predicted to be near 15 Schwarzschild radii for some snapshots. We show that our results are consistent with existing observational data such as timeaveraged radio spectra, VLBA images and variability timescale for M87. We also discuss observations of M87 with millimeter and sub-millimeter interferometers.

We report two different types of backflow from jets by performing 2D special relativistic hydrodynamical simulations. One is anti-parallel and quasi-straight to the main jet (quasi-straight backflow), and the other is bent path of the backflow (bent backflow). We find that the former appears when the head advance speed is comparable to or higher than the local sound speed at the hotspot while the latter appears when the head advance speed is slower than the sound speed bat the hotspot. Bent backflow collides with the unshocked jet and laterally squeezes the jet. At the same time, a pair of new oblique shocks are formed at the tip of the jet and new bent fast backflows are generated via these oblique shocks. The hysteresis of backflow collisions is thus imprinted in the jet as a node and anti-node structure. This process also promotes broadening of the jet cross sectional area and it also causes a decrease in the head advance velocity. This hydrodynamic process may be tested by observations of compact young jets.

This paper is a sequel to our previous work for accretion onto a Schwarzschild black hole and the so-called standing accretion shock instability (SASI), in this paper we investigate non-axisymmetric perturbations for a Kerr black hole. The linear and non-linear phases for the shock evolution are analyzed in detail by both 2D general relativistic hydrodynamical simulations and linear analysis. Since the structure of steady axisymmetric accretion flows with a standing shock wave is very sensitive to the inner transonic flow, their properties such as Mach numbers, which are important for the stability, depend on the Kerr parameter very much. Although the essential features of the instability do not differ from the previous results for the Schwarzschild black hole, the frame dragging effects specific to the Kerr black hole is also evident. Interestingly, the oscillation periods of the fundamental unstable modes are dependent only on the shock radius irrespective of the injection parameters.

We study the stability of standing shock waves in advection-dominated accretion flows into a Schwarzschild black hole by 2D general relativistic hydrodynamic simulations as well as linear analysis in the equatorial plane. We demonstrate that the accretion shock is stable against axisymmetric perturbations but becomes unstable to non-axisymmetric perturbations. The results of dynamical simulations show good agreement with linear analysis on the stability, oscillation and growing time scales. The comparison of different wave-travel times with the growth time scales of the instability suggests that the instability is likely to be of the Papaloizou-Pringle type, induced by the repeated propagations of acoustic waves. However, the wavelengths of perturbations are too long to clearly define the reflection point. By analyzing the non-linear phase in the dynamical simulations, it is shown that quadratic mode couplings precede the non-linear saturation. It is also found that not only short-term random fluctuations by turbulent motions but also quasi periodic oscillations take place on longer time scales in the non-linear phase. We give some possible implications of the instability for quasi periodic oscillations (QPOs) and the central engine for gamma ray bursts (GRBs).