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 investigate explosive nucleosynthesis in a non-rotating 15$M_\odot$ star with solar metallicity that explodes by a neutrino-heating supernova (SN) mechanism aided by both standing accretion shock instability (SASI) and convection. To trigger explosions in our two-dimensional hydrodynamic simulations, we approximate the neutrino transport with a simple light-bulb scheme and systematically change the neutrino fluxes emitted from the protoneutron star. By a post-processing calculation, we evaluate abundances and masses of the SN ejecta for nuclei with the mass number $\le 70$ employing a large nuclear reaction network. Aspherical abundance distributions, which are observed in nearby core-collapse SN remnants, are obtained for the non-rotating spherically-symmetric progenitor, due to the growth of low-mode SASI. Abundance pattern of the supernova ejecta is similar to that of the solar system for models whose masses ranges $(0.4-0.5) \Ms$ of the ejecta from the inner region ($\le 10,000\km$) of the precollapse core. For the models, the explosion energies and the \nuc{Ni}{56} masses are $ \simeq 10^{51} \rm erg$ and $(0.05-0.06) \Ms$, respectively; their estimated baryonic masses of the neutron star are comparable to the ones observed in neutron-star binaries. These findings may have little uncertainty because most of the ejecta is composed by matter that is heated via the shock wave and has relatively definite abundances. The abundance ratios for Ne, Mg, Si and Fe observed in Cygnus loop are well reproduced with the SN ejecta from an inner region of the $15\Ms$ progenitor.

By performing three-dimensional (3D) simulations that demonstrate the neutrino-driven core-collapse supernovae aided by the standing accretion shock instability (SASI), we study how the spiral modes of the SASI can have impacts on the properties of the gravitational-wave (GW) emission. To see the effects of rotation in the non-linear postbounce phase, we give a uniform rotation on the flow advecting from the outer boundary of the iron core, whose specific angular momentum is assumed to agree with recent stellar evolution models. We compute fifteen 3D models in which the initial angular momentum as well as the input neutrino luminosities from the protoneutron star are changed in a systematic manner. By performing a ray-tracing analysis, we accurately estimate the GW amplitudes generated by anisotropic neutrino emission. Our results show that the gravitational waveforms from neutrinos in models that include rotation exhibit a common feature otherwise they vary much more stochastically in the absence of rotation. The breaking of the stochasticity stems from the excess of the neutrino emission parallel to the spin axis. This is because the compression of matter is more enhanced in the vicinity of the equatorial plane due to the growth of the spiral SASI modes, leading to the formation of spiral flows circulating around the spin axis with higher temperatures. We point out that a recently proposed future space interferometers like Fabry-Perot type DECIGO would permit detection of these signals for a Galactic supernova.

We propose a ray-tracing method to estimate gravitational waves (GWs) generated by anisotropic neutrino emission in supernova cores. To calculate the gravitational waveforms, we derive analytic formulae in a useful form, which are applicable also for three-dimensional computations. Pushed by evidence of slow rotation prior to core-collapse, we focus on asphericities in neutrino emission and matter motions outside the protoneutron star. Based on the two-dimensional (2D) models, which mimic SASI-aided neutrino heating explosions, we compute the neutrino anisotropies via the ray-tracing method in a post-processing manner and calculate the resulting waveforms. With these computations, it is found that the waveforms exhibit more variety in contrast to the ones previously estimated by the ray-by-ray analysis (e.g., Kotake et al. (2007)). In addition to a positively growing feature, which was predicted to determine the total wave amplitudes predominantly, the waveforms are shown to exhibit large negative growth for some epochs during the growth of SASI. These features are found to stem from the excess of neutrino emission in lateral directions, which can be precisely captured by the ray-tracing calculation. Due to the negative contributions and the neutrino absorptions appropriately taken into account by the ray-tracing method, the wave amplitudes become more than one-order-of magnitude smaller than the previous estimation, thus making their detections very hard for a galactic source.On the other hand, it is pointed out that the GW spectrum from matter motions have its peak near $\sim 100$ Hz, which could be characteristic for the SASI-induced supernova explosions.(abridged)

We study properties of gravitational waves based on the three-dimensional simulations, which demonstrate the neutrino-driven explosions aided by the standing accretion shock instability (SASI). Pushed by evidence supporting slow rotation prior to core-collapse, we focus on the asphericities in neutrino emissions and matter motions outside the protoneutron star. By performing a ray-tracing calculation in 3D, we estimate accurately the gravitational waveforms from anisotropic neutrino emissions. In contrast to the previous work assuming axisymmetry, we find that the gravitational waveforms vary much more stochastically because the explosion anisotropies depend sensitively on the growth of the SASI which develops chaotically in all directions. Our results show that the gravitational-wave spectrum has its peak near $\sim 100$ Hz, reflecting the SASI-induced matter overturns of $\sim O(10)$ ms. We point out that the detection of such signals, possibly visible to the LIGO-class detectors for a Galactic supernova, could be an important probe into the long-veiled explosion mechanism.

We studied the effects of rotation on standing accretion shock instability (SASI) by performing three-dimensional hydrodynamics simulations. Taking into account a realistic equation of state and neutrino heating/cooling, we prepared a spherically symmmetric and steady accretion flow through a standing shock wave onto a proto-neutron star (PNS). When the SASI entered the nonlinear phase, we imposed uniform rotation on the flow advecting from the outer boundary of the iron core, whose specific angular momentum was assumed to agree with recent stellar evolution models. Using spherical harmonics in space and Fourier decompositions in time, we performed mode analysis of the nonspherical deformed shock wave to observe rotational effects on the SASI in the nonlinear phase. We found that rotation imposed on the axisymmetric SASI did not make any spiral modes and hardly affected sloshing modes, except for steady l=2, m=0 modes. In contrast, rotation imposed on the non-axisymmetric flow increased the amplitude of spiral modes so that some spiral flows accreting on the PNS were more clearly formed inside the shock wave than without rotation. The amplitudes of spiral modes increased significantly with rotation in the progressive direction.

We have studied non-axisymmetric standing accretion shock instability, or SASI, by 3D hydrodynamical simulations. This is an extention of our previous study on axisymmetric SASI. We have prepared a spherically symmetric and steady accretion flow through a standing shock wave onto a proto-neutron star, taking into account a realistic equation of state and neutrino heating and cooling. This unperturbed model is supposed to represent approximately the typical post-bounce phase of core-collapse supernovae. We then have added a small perturbation (~1%) to the radial velocity and computed the ensuing evolutions. Not only axisymmetric but non-axisymmetric perturbations have been also imposed. We have applied mode analysis to the non-spherical deformation of the shock surface, using the spherical harmonics. We have found that (1) the growth rates of SASI are degenerate with respect to the azimuthal index m of the spherical harmonics Y_l^m, just as expected for a spherically symmetric background, (2) nonlinear mode couplings produce only m=0 modes for the axisymmetric perturbations, whereas m=!0 modes are also generated in the non-axisymmetric cases according to the selection rule for the quadratic couplings, (3) the nonlinear saturation level of each mode is lower in general for 3D than for 2D because a larger number of modes are contributing to turbulence in 3D, (4) low l modes are dominant in the nonlinear phase, (5) the equi-partition is nearly established among different m modes in the nonlinear phase, (6) the spectra with respect to l obey power laws with a slope slightly steeper for 3D, and (7) although these features are common to the models with and without a shock revival at the end of simulation, the dominance of low l modes is more remarkable in the models with a shock revival.

The so-called "acoustic revival mechanism" of core-collapse supernova proposed recently by the Arizona group is an interesting new possibility. Aiming to understand the elementary processes involved in the mechanism, we have calculated the eigen frequencies and eigen functions for the g-mode oscillations of a non-rotating proto neutron star. The possible excitation of these modes by the standing accretion shock instability, or SASI, is discussed based on these eigen functions. We have formulated the forced oscillations of $g$-modes by the external pressure perturbations exerted on the proto neutron star surface. The driving pressure fluctuations have been adopted from our previous computations of the axisymmetric SASI in the non-linear regime. We have paid particular attention to low l modes, since these are the modes that are dominant in SASI and that the Arizona group claimed played an important role in their acoustic revival scenario. Here l is the index of the spherical harmonic functions, $Y_l^m$. Although the frequency spectrum of the non-linear SASI is broadened substantially by non-linear couplings, the typical frequency is still much smaller than those of g-modes, the fact leading to a severe impedance mismatch. As a result, the excitations of various $g$-modes are rather inefficient and the energy of the saturated g-modes is $\sim 10^{50}$erg or smaller, with the g_2-mode being the largest in our model. Here the g_2-mode has two radial nodes and is confined to the interior of the convection region. The energy transfer rate from the g-modes to out-going sound waves is estimated from the growth of the g-modes and found to be $\sim 10^{51}$erg/s in the model studied in this paper.

We present the results of numerical experiments, in which we study how the asphericities induced by the growth of the standing accretion shock instability (SASI) produce the gravitational waveforms in the postbounce phase of core-collapse supernovae. To obtain the neutrino-driven explosions, we parameterize the neutrino fluxes emitted from the central protoneutron star and approximate the neutrino transfer by a light-bulb scheme. We find that the waveforms due to the anisotropic neutrino emissions show the monotonic increase with time, whose amplitudes are up to two order-of-magnitudes larger than the ones from the convective matter motions outside the protoneutron stars. We point out that the amplitudes begin to become larger when the growth of the SASI enters the nonlinear phase, in which the deformation of the shocks and the neutrino anisotropy become large. From the spectrum analysis of the waveforms, we find that the amplitudes from the neutrinos are dominant over the ones from the matter motions at the frequency below $\sim 100$ Hz, which are suggested to be within the detection limits of the detectors in the next generation such as LCGT and the advanced LIGO for a supernova at 10 kpc. As a contribution to the gravitational wave background, we show that the amplitudes from this source could be larger at the frequency above $\sim$ 1 Hz than the primordial gravitational wave backgrounds, but unfortunately, invisible to the proposed space-based detectors.

We present the results of numerical experiments, in which we have investigated the influence of the inelastic neutrino-helium interactions on the standing accretion shock instability supposed to occur in the post-bounce supernova core. The axisymmetric hydrodynamical simulations of accretion flows through the standing accretion shock wave onto the protoneutron star show that the interactions are relatively minor and the linear growth of the shock instability is hardly affected. The extra heating given by the inelastic reactions becomes important for the shock revival after the instability enters the non-linear regime, but only when the neutrino luminosity is very close to the critical value, at which the shock would be revived without the interactions. We have also studied the dependence of the results on the initial amplitudes of perturbation and the temperatures of mu and tau neutrinos.

We report a current status of our radiation-magnetohydrodynamic code for the study of core-collapse supernovae. In this contribution, we discuss the accuracy of our newly developed numerical code by presenting the test problem in a static background model. We also present the application to the spherically symmetric core-collapse simulations. Since close comparison with the previously published models is made, we are now applying it for the study of magnetorotational core-collapse supernovae.

We have numerically studied the instability of the spherically symmetric standing accretion shock wave against non-spherical perturbations. We have in mind the application to the collapse-driven supernovae in the post bounce phase, where the prompt shock wave generated by core bounce is commonly stalled. We take an experimental stand point in this paper. Using spherically symmetric, completely steady, shocked accretion flows as unperturbed states, we have clearly observed both the linear growth and the subsequent nonlinear saturation of the instability. In so doing, we have employed a realistic equation of state together with heating and cooling via neutrino reactions with nucleons. We have done a mode analysis based on the spherical harmonics decomposition and found that the modes with l=1, 2 are dominant not only in the linear regime, but also after the nonlinear couplings generate various modes and the saturation occurs. Varying the neutrino luminosity, we have constructed the unperturbed states both with and without a negative entropy-gradient. We have found that in both cases the growth of the instability is similar, suggesting the convection does not play a dominant role, which also appears to be supported by the recent linear analysis of the convection in accretion flows by Foglizzo et al. The real part of the eigen frequency seems to be mainly determined by the advection time rather than by the sound-crossing time. Whatever the cause may be, the instability is favorable for the shock revival.

Two-dimensional numerical simulations of an accretion flow in a close binary system are performed by solving the Euler equations with radiative transfer. In the present study, the specific heat ratio is assumed to be constant while radiative cooling effect is included as a non-adiabatic process. The cooling effect of the disc is considered by discharging energy in the vertical directions from the top and bottom surfaces of the disc. We use the flux-limited diffusion approximation to calculate the radiative heat flux values. Our calculations show that a disc structure appears and the spiral shocks are formed on the disc. These features are similar to that observed in the case of an adiabatic gas with a lower specific heat ratio, $\gamma=1.01$. It is found that when radiative cooling effect is accounted for, the mass of the disc becomes larger than that assuming $\gamma=5/3$, and smaller than that assuming $\gamma=1.01$. It is concluded that employing an adiabatic gas with a lower specific heat ratio is almost a valid assumption for simulating accretion disc with radiative cooling effect.