Mixing of 56Ni, whose nuclear decay energy is a major luminosity source in stripped-envelope supernovae, is known to affect the observational properties of stripped-envelope supernovae such as light-curve and color evolution. Here we systematically investigate the effect of 56Ni mixing on the photospheric velocity evolution in stripped-envelope supernovae. We show that 56Ni mixing significantly affects the early photospheric velocity evolution. The photospheric velocity, which is often used to constrain the ejecta mass and explosion energy, significantly varies by just changing the degree of 56Ni mixing. In addition, the models with a small degree of 56Ni mixing show a flattening in the early photospheric velocity evolution, while the fully mixed models show a monotonic decrease. The velocity flattening appears in both helium and carbon+oxygen progenitor explosions with a variety of ejecta mass, explosion energy, and 56Ni mass. Some stripped-envelope supernovae with early photospheric velocity information do show such a flattening. We find that Type Ib SN 2007Y, which has early photospheric velocity information, has a signature of a moderate degree of 56Ni mixing in the photospheric velocity evolution and about a half of the ejecta is mixed in it. The immediate spectroscopic follow-up observations of stripped-envelope supernovae shortly after the explosion providing the early photospheric evolution give an important clue to constrain 56Ni mixing in the ejecta.

A number of astrophyical transients originating from stellar explosions are powered by the collision of the ejected material with the circumstellar medium, which efficiently produces thermal radiation via shock dissipation. We investigate how such interaction-powered transients are distributed in the peak bolometric luminosity vs the rise time phase space. Taking the advantage of less time-consuming one-dimensional simulations with spherical symmetry, we calculated more than 500 models with different circumstellar mass and radius, ejecta mass and energy, and chemical compositions. The peak bolometric luminosity, the total radiated energy, and the rise time of the interaction-powered emission are measured for each simulated light curve. We consider how these characteristic quantities are determined as a function of the model parameters and discuss possible implications for the observed populations of (potential) interaction-powered transients, such as type IIn supernovae and fast blue optical transients.

We argue that the neutrino halo, a population of neutrinos that have undergone direction-changing scattering in the stellar envelope of a core-collapse supernova (CCSNe), is sensitive to neutrino emission history through time of flight. We show that the constant time approximation, commonly used in calculating the neutrino halo, does not capture the spatiotemporal evolution of the halo neutrino population and that correcting for time of flight can produce conditions which may trigger fast neutrino flavor conversion. We also find that there exists a window of time early in all CCSNe where the neutrino halo population is sufficiently small that it may be negligible. This suggests that collective neutrino oscillation calculations which neglect the Halo may be well founded at sufficiently early times.

We extend the multi-angle computational framework and investigate the time evolution of the neutrino halo on collective neutrino oscillation in the core collapse of an iron core progenitor. We find that in the case of the $9.6\, \rm M_\odot$ progenitor adopted in this work, there are windows of time when the effects of neutrino halo and collective neutrino oscillation are not simultaneously large. Inside the shock, the impact of the inward-scattered halo neutrino cannot in general be neglected compared to the outward-propagating neutrino flux. However, during early epochs, collective neutrino oscillation is effectively shut down by multi-angle matter suppression. During the intermediate epoch, collective neutrino oscillation is not suppressed, but its onset radius is beyond the still relatively small explosion shock front where the halo is prominent. We also find in the case of the $9.6\, \rm M_\odot$ progenitor the halo neutrinos induce a delay in the onset of collective neutrino oscillations. This causes novel flavor conversions which sharpen collective neutrino oscillation spectral features. We predict that the inclusion of neutrino halo effects makes neutrino signals that are more clearly distinct from thermal emission that when halo neutrinos are omitted.

We present results of three-dimensional (3D), radiation-magnetohydrodynamics (MHD) simulations of core-collapse supernovae in full general relativity (GR) with spectral neutrino transport. In order to study the effects of progenitor's rotation and magnetic fields, we compute three models, where the precollapse rotation rate and magnetic fields are included parametrically to a 20 M$_{\odot}$ star. While we find no shock revival in our two non-magnetized models during our simulation times ($\sim500$ ms after bounce), the magnetorotationally (MR) driven shock expansion immediately initiates after bounce in our rapidly rotating and strongly magnetized model. We show that the expansion of the MR-driven flows toward the polar directions is predominantly driven by the magnetic pressure, whereas the shock expansion toward the equatorial direction is supported by neutrino heating. Our detailed analysis indicates that the growth of the so-called kink instability may hinder the collimation of jets, resulting in the formation of broader outflows. Furthermore we find a dipole emission of lepton number, only in the MR explosion model, whose asymmetry is consistent with the explosion morphology. Although it is similar to the lepton-number emission self-sustained asymmetry (LESA), our analysis shows that the dipole emission occurs not from the protoneutron star convection zone but from above the neutrino sphere indicating that it is not associated with the LESA. We also report several unique neutrino signatures, which are significantly dependent on both the time and the viewing angle, if observed, possibly providing a rich information regarding the onset of the MR-driven explosion.

We examine the eigenfrequencies of gravitational waves from the protoneutron stars (PNSs) provided via the core-collapse supernovae, focusing on how the frequencies depend on the dimension of numerical simulations especially for the early phase after core bounce. As expected, we find that the time evolution of gravitational wave frequencies depends strongly on the dimension of numerical simulation as well as the equation of state for high density matter. Even so, we find that the fundamental frequency as a function of PNS average density and the ratio of the specific eigenfrequencies to the fundamental frequency as a function of PNS properties are independent of the dimension of the numerical simulations, where the dependence of the equation of state is also very weak in this early postbounce phase. Thus, one can safely discuss the gravitational wave frequencies of the PNSs as a function of the PNS average density or compactness even with the frequencies obtained from the one dimensional simulations. We also provide phenomenological relation between the compactness and average density as well as the relation among the $f$-, $p$-, and $g$- mode frequencies.

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.

A dense neutrino gas, such as the one anticipated in the supernova environment, can experience fast neutrino flavor conversions on scales much shorter than those expected in vacuum probably provided that the angular distributions of $\nu_e$ and $\bar\nu_e$ cross each other. We perform a detailed investigation of the neutrino angular distributions obtained by solving the Boltzmann equations for fixed matter profiles of some representative snapshots during the post-bounce phase of core-collapse supernovae in multidimensional calculations of an $11.2\mathrm{M}_{\odot}$ and a $27\mathrm{M}_{\odot}$ progenitor models. Although the $11.2\mathrm{M}_{\odot}$ model features $\nu_e - \bar\nu_e$ angular crossings and the associated fast modes at different time snapshots, the $27\mathrm{M}_{\odot}$ model does not show any crossings within the decoupling region. We show that this can be understood by studying the multipole components of the neutrino distributions. In fact, there is a higher chance for the occurrence of $\nu_e - \bar\nu_e$ angular crossings for the zones where the multipole components of the neutrino distributions are strong enough. We also show that there can exist more than one crossings between the angular distributions of $\nu_e$ and $\bar{\nu}_e$. In addition, apart from the crossings within the neutrino decoupling region, there is a class of $\nu_e-\bar\nu_e$ angular crossings which appears very deep inside the proto-neutron star.

We present results from a full general relativistic three-dimensional hydrodynamics simulation of rapidly rotating core-collapse of a 70 M$_{\odot}$ star with three-flavor spectral neutrino transport. We find a strong gravitational wave (GW) emission that originates from the growth of the one- and two-armed spiral waves extending from the nascent proto-neutron star (PNS). The GW spectrogram shows several unique features that are produced by the non-axisymmetric instabilities. After bounce, the spectrogram first shows a transient quasi-periodic time modulation at $\sim$ 450 Hz. In the second active phase, it again shows the quasi-periodic modulation but with the peak frequency increasing with time, which continues until the final simulation time. From our detailed analysis, such features can be well explained by a combination of the so-called low-$T/|W|$ instability and the PNS core contraction.

We report results from a convection dynamo simulation of proto-neutron star (PNS), with a nuclear equation of state (EOS) and the initial hydrodynamic profile taken from a neutrino radiation-hydrodynamics simulation of a massive stellar core-collapse. A moderately-rotating PNS with the spin period of $170$ ms in the lepton-driven convection stage is focused. We find that large-scale flow and thermodynamic fields with north-south asymmetry develop in the turbulent flow, as a consequence of the convection in the central part of the PNS, which we call as a "deep core convection". Intriguingly, even with such a moderate rotation, large-scale, $10^{15}$ G, magnetic field with dipole symmetry is spontaneously built up in the PNS. The turbulent electro-motive force arising from rotationally-constrained core convection is shown to play a key role in the large-scale dynamo. The large-scale structures organized in the PNS may impact the explosion dynamics of supernovae and subsequent evolution to the neutron stars.

The dense neutrino medium in a core-collapse supernova or a neutron-star merger event can experience fast flavor conversions on time/distance scales that are much smaller than those of vacuum oscillations. It is believed that fast neutrino flavor transformation occurs in the region where the angular distributions of $\nu_e$ and $\bar\nu_e$ cross each other. We present the first study of this crossing phenomenon and the fast neutrino flavor conversions in multidimensional (multi-D) supernova models. We examine the neutrino distributions obtained by solving the Boltzmann transport equation for several fixed profiles which are representative snapshots taken from separate 2D and 3D supernova simulations with an $11.2 M_\odot$ progenitor model. Our research shows that the spherically asymmetric patterns of the $\nu_e$ and $\bar\nu_e$ fluxes in multi-D models can assist the appearance of the crossing between the $\nu_e$ and $\bar\nu_e$ angular distributions. In the models that we have studied, there exist unstable neutrino oscillation modes in and beyond the neutrino decoupling region which have amplitude growth rates as large as an $e$-fold per nanosecond. This finding can have important consequences for the explosion mechanism, nucleosynthesis, and neutrino signals of core-collapse supernovae.

We perform a series of two-dimensional radiation-hydrodynamic simulations of the collision between supernova ejecta and circumstellar media (CSM). The hydrodynamic interaction of a fast flow and the surrounding media efficiently dissipates the kinetic energy of the fast flow and considered as a dominant energy source for a specific class of core-collapse supernovae. Despite some observational evidence for aspherical ejecta and/or CSM structure, multi-dimensional effects in the ejecta-CSM interaction are relatively unexplored. Our numerical simulations equipped with an adaptive mesh refinement technique successfully reproduce hydrodynamic instabilities developing around the ejecta-CSM interface. We also investigate effects of disk-like CSM on the dynamical evolution of supernova ejecta and bolometric light curves. We find that emission powered by ejecta-disk interaction exhibits significant viewing angle dependence. For a line of sight close to the symmetry axis, the observer directly sees the supernova ejecta, leading to a short brightening timescale. For an observer seeing the emission through the CSM disk, thermal photons diffuse throughout the CSM and thus the light curve is severely smeared out.

Core-collapse supernovae (CCSNe) are the final stage of massive stars, marking the birth of neutron stars (NSs). The aspherical mass ejection drives a natal kick of the forming NS. In this work we study the properties of the NS kick based on our long-term hydrodynamics CCSN simulations. We perform two-dimensional (2D) simulations for ten progenitors from a 10.8 to 20 $\, M_{\odot}$ star covering a wide range of the progenitor's compactness parameter, and two three-dimensional (3D) simulations for an 11.2 $\, M_{\odot}$ star. Our 2D models present a variety of explosion energies between $\sim 1.3 \times 10^{50}$ erg and $\sim 1.2 \times 10^{51}$ erg, and NS kick velocities between $\sim 100$ km s$^{-1}$ and $\sim 1500$ km s$^{-1}$. For the 2D exploding models, we find that the kick velocities tend to become higher with the progenitor's compactness. This is because the high progenitor compactness results in high neutrino luminosity from the proto-neutron star (PNS), leading to more energetic explosions. Since high-compactness progenitors produce massive PNSs, we point out that the NS masses and the kick velocities can be correlated, which is moderately supported by observation. Comparing 2D and 3D models of the 11.2 $\, M_{\odot}$ star, the diagnostic explosion energy in 3D is, as previously identified, higher than that in 2D, whereas the 3D model results in a smaller asymmetry in the ejecta distribution and a smaller kick velocity than in 2D. Our results confirm the importance of self-consistent CCSN modeling covering a long-term postbounce evolution in 3D for a quantitative prediction of the NS kicks.

We perform two- (2D) and three-dimensional (3D) hydrodynamics simulations of convective oxygen shell-burning that takes place deep inside a massive progenitor star of a core-collapse supernova. Using one dimensional (1D) stellar evolution code, we first calculate the evolution of massive stars with an initial mass of 9-40 $M_\odot$. Four different overshoot parameters are applied, and CO core mass trend similar to previous works is obtained in the 1D models. Selecting eleven 1D models that have a silicon and oxygen coexisting layer, we perform 2D hydrodynamics simulations of the evolution $\sim$100 s until the onset of core-collapse. We find that convection with large-scale eddies and the turbulent Mach number $\sim$0.1 is obtained in the models having a Si/O layer with a scale of 10$^8$ cm, whereas most models that have an extended O/Si layer up to a few $\times 10^9$ cm exhibit lower turbulent velocity. Our results indicate that the supernova progenitors that possess a thick Si/O layer could provide a preferable condition for perturbation-aided explosions. We perform 3D simulation of a 25 $M_\odot$ model, which exhibits large-scale convection in the 2D models. The 3D model develops large ($\ell = 2$) convection similar to the 2D model, however, the turbulent velocity is lower. By estimating the neutrino emission properties of the 3D model, we point out that a time modulation of the event rates, if observed in KamLAND and Hyper-Kamiokande, would provide an important information about structural changes in the presupernova convective layer.

To obtain the eigenfrequencies of a protoneutron star (PNS) in the postbounce phase of core-collapse supernovae (CCSNe), we perform a linear perturbation analysis of the angle-averaged PNS profiles using results from a general relativistic CCSN simulation of a $15 M_{\odot}$ star. In this work, we investigate how the choice of the outer boundary condition could affect the PNS oscillation modes in the linear analysis. By changing the density at the outer boundary of the PNS surface in a parametric manner, we show that the eigenfrequencies strongly depend on the surface density. By comparing with the gravitational wave (GW) signatures obtained in the hydrodynamics simulation, the so-called surface $g$-mode of the PNS can be well ascribed to the fundamental oscillations of the PNS. The frequency of the fundamental oscillations can be fitted by a function of the mass and radius of the PNS similar to the case of cold neutron stars. In the case that the position of the outer boundary is chosen to cover not only the PNS but also the surrounding postshock region, we obtain the eigenfrequencies close to the modulation frequencies of the standing accretion-shock instability (SASI). However, we point out that these oscillation modes are unlikely to have the same physical origin of the SASI modes seen in the hydrodynamics simulation. We discuss possible limitations of applying the angle-averaged, linear perturbation analysis to extract the full ingredients of the CCSN GW signatures.

Macroscopic evolution of relativistic charged matter with chirality imbalance is described by the chiral magnetohydrodynamics (chiral MHD). One such astrophysical system is high-density lepton matter in core-collapse supernovae where the chirality imbalance of leptons is generated by the parity-violating weak processes. After developing the chiral MHD equations for this system, we perform numerical simulations for the real-time evolutions of magnetic and flow fields, and study the properties of the chiral MHD turbulence. In particular, we observe the inverse cascade of the magnetic energy and the fluid kinetic energy. Our results suggest that the chiral effects that have been neglected so far can reverse the turbulent cascade direction from direct to inverse cascade, which would impact the magnetohydrodynamics evolution in the supernova core toward explosion.

Recent core-collapse supernova (CCSN) simulations have predicted several distinct features in gravitational-wave (GW) spectrograms, including a ramp-up signature due to the g-mode oscillation of the proto-neutron star (PNS) and an excess in the low-frequency domain (100-300 Hz) potentially induced by the standing accretion shock instability (SASI). These predictions motivated us to perform a sophisticated time-frequency analysis (TFA) of the GW signals, aimed at preparation for future observations. By reanalyzing a gravitational waveform obtained in a three-dimensional general-relativistic CCSN simulation, we show that both the spectrogram with an adequate window and the quadratic TFA separate the multimodal GW signatures much more clearly compared with the previous analysis. We find that the observed low-frequency excess during the SASI active phase is divided into two components, a stronger one at 130 Hz and an overtone at 260 Hz, both of which evolve quasi-statically during the simulation time. We also identify a new mode whose frequency varies from 700 to 600 Hz. Furthermore, we develop the quadratic TFA for the Stokes I, Q, U, and V parameters as a new tool to investigate the GW circular polarization. We demonstrate that the polarization states that randomly change with time after bounce are associated with the PNS g-mode oscillation, whereas a slowly changing polarization state in the low-frequency domain is connected to the PNS core oscillation. This study demonstrates the capability of the sophisticated TFA for diagnosing the polarized CCSN GWs in order to explore their complex nature.

We investigate a progenitor mass distribution of core-collapse supernova remnants (CCSNRs) in our Galaxy and the Large and Small Magellanic Clouds, for the first time. We count the number of CCSNRs in three mass ranges divided by the zero-age main-sequence mass, $M_{\rm ZAMS}$; A: $M_{\rm ZAMS} < 15\ {\rm M}_\odot$, B: $15\ {\rm M}_\odot < M_{\rm ZAMS} < 22.5\ {\rm M}_\odot$, C: $M_{\rm ZAMS} > 22.5\ {\rm M}_\odot$. Simple compilation of progenitor masses in the literature yields a progenitor mass distribution of $f_{\rm A}: f_{\rm B}: f_{\rm C} =0.24:0.28:0.48$, where $f$ is the number fraction of the progenitors. The distribution is inconsistent with any standard initial mass functions. We notice, however, that previous mass estimates are subject to large systematic uncertainties because most of the relative abundances (X/Si) are not really good probe for the progenitor masses. Instead, we propose to rely only on the Fe/Si ratio which is sensitive to the CO core mass ($M_{\rm COcore}$) and $M_{\rm ZAMS}$. Comparing Fe/Si ratios in SNRs in the literature with the newest theoretical model, we estimate 33 $M_{\rm COcore}$ and $M_{\rm ZAMS}$, leading to a revised progenitor mass distribution of $f_{\rm A}: f_{\rm B}: f_{\rm C} = 0.47: 0.32 : 0.21$. This is consistent with the standard Salpeter initial mass function. However, the relation between $M_{\rm COcore}$ and $M_{\rm ZAMS}$ could be affected by binary evolution, which is not taken into account in this study and should be considered in the future work to derive a better progenitor mass distribution estimate.

We present a comparison between several simulation codes designed to study the core-collapse supernova mechanism. We pay close attention to controlling the initial conditions and input physics in order to ensure a meaningful and informative comparison. Our goal is three-fold. First, we aim to demonstrate the current level of agreement between various groups studying the core-collapse supernova central engine. Second, we desire to form a strong basis for future simulation codes and methods to compare to. Lastly, we want this work to be a stepping stone for future work exploring more complex simulations of core-collapse supernovae, i.e., simulations in multiple dimensions and simulations with modern neutrino and nuclear physics. We compare the early (first ~500ms after core bounce) spherically-symmetric evolution of a 20 solar mass progenitor star from six different core-collapse supernovae codes: 3DnSNe-IDSA, AGILE-BOLTZTRAN, FLASH, F{\sc{ornax}}, GR1D, and PROMETHEUS-VERTEX. Given the diversity of neutrino transport and hydrodynamic methods employed, we find excellent agreement in many critical quantities, including the shock radius evolution and the amount of neutrino heating. Our results provide an excellent starting point from which to extend this comparison to higher dimensions and compare the development of hydrodynamic instabilities that are crucial to the supernova explosion mechanism, such as turbulence and convection.

It has been suggested that whether a star explodes or not, and what kind of explosion properties it shows, is strongly dependent on the progenitor's core structure. We present the results from 101 axisymmetric core-collapse supernova simulations performed with progenitors spanning initial masses in the range from 10.8 to 75 solar masses, and focus on their connections to the compactness of the progenitor's core. Our simulations confirm a correlation between the neutrinos emitted during the accretion phase and the progenitor's compactness. We suggest that the ratio of observed neutrino events during the first hundreds of milliseconds can be used to infer the progenitor's inner mass density structure.

We study the final fate of a very massive star by performing full general relativistic (GR), three-dimensional (3D) simulation with three-flavor multi-energy neutrino transport. Utilizing a 70 solar mass zero metallicity progenitor, we self-consistently follow the radiation-hydrodynamics from the onset of gravitational core-collapse until the second collapse of the proto-neutron star (PNS), leading to black hole (BH) formation. Our results show that the BH formation occurs at a post-bounce time of ~300 ms for the 70 Msun star. This is significantly earlier than those in the literature where lower mass progenitors were employed. At a few ~10 ms before BH formation, we find that the stalled bounce shock is revived by intense neutrino heating from the very hot PNS, which is aided by violent convection behind the shock. In the context of 3D-GR core-collapse modeling with multi-energy neutrino transport, our numerical results present the first evidence to validate a fallback BH formation scenario of the 70 Msun star.

The birth properties of neutron stars yield important information on the still debated physical processes that trigger the explosion and on intrinsic neutron-star physics. These properties include the high space velocities of young neutron stars with average values of several 100 km/s, whose underlying "kick" mechanism is not finally clarified. There are two competing possibilities that could accelerate NSs during their birth: anisotropic ejection of either stellar debris or neutrinos. We here present new evidence from X-ray measurements that chemical elements between silicon and calcium in six young gaseous supernova remnants are preferentially expelled opposite to the direction of neutron star motion. There is no correlation between the kick velocities and magnetic field strengths of these neutron stars. Our results support a hydrodynamic origin of neutron-star kicks connected to asymmetric explosive mass ejection, and they conflict with neutron-star acceleration scenarios that invoke anisotropic neutrino emission caused by particle and nuclear physics in combination with very strong neutron-star magnetic fields.

We present an analysis of the circular polarization of gravitational-waves (GWs) using results from three-dimensional (3D), general relativistic (GR) core-collapse simulations of a non-rotating $15 M_{\odot}$ star. For the signal detection, we perform a coherent network analysis taking into account the four interferometers of LIGO Hanford, LIGO Livingston, VIRGO, and KAGRA. We focus on the Stokes $V$ parameter, which directly characterizes the asymmetry of the GW circular polarization. We find that the amplitude of the GW polarization becomes bigger for our 3D-GR model that exhibits strong activity of the standing accretion shock instability (SASI). Our results suggest that the SASI-induced accretion flows to the proto-neutron star (PNS) lead to a characteristic, low-frequency modulation (100 $\sim$ 200 Hz) in both the waveform and the GW circular polarization. During the vigorous SASI phase, we observe that the GW polarization switches from the right- to left-handed mode, which is clearly visible in the spectrogram. By estimating the signal-to-noise ratio of the GW polarization, we demonstrate that the detection horizon of the circular polarization extends by more than a factor of several times farther comparing to that of the GW amplitude. Our results suggest that the GW circular polarization, if detected, could provide a new probe into the pre-explosion hydrodynamics such as the SASI activity and the $g$-mode oscillation of the PNS.

Accurate description of neutrino opacities is central both to the core-collapse supernova (CCSN) phenomenon and to the validity of the explosion mechanism itself. In this work, we study in a systematic fashion the role of a variety of well-selected neutrino opacities in CCSN simulations where multi-energy, three-flavor neutrino transport is solved by the isotropic diffusion source approximation (IDSA) scheme. To verify our code, we first present results from one-dimensional (1D) simulations following core-collapse, bounce, and up to ~ 250 ms postbounce of a 15 $M_{\odot}$ star using a standard set of neutrino opacities by Bruenn (1985). Detailed comparison with published results supports the reliability of our three-flavor IDSA scheme using the standard opacity set. We then investigate in 1D simulations how the individual opacity update leads to the difference from the base-line run with the standard opacity set. By making a detailed comparison with previous work, we check the validity of our implementation of each update in a step-by-step manner. Individual neutrino opacities with the largest impact on the overall evolution in 1D simulations are selected for a systematic comparison in our two-dimensional (2D) simulations. Special emphasis is devoted to the criterion of explodability in the 2D models. We discuss the implications of these results as well as the limitations and requirements for future towards more elaborate CCSN modeling.

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 analysis of neutrino and GW signals based on three-dimensional (3D) core-collapse supernova simulations of a rapidly rotating 27 $M_{\odot}$ star. We find a new neutrino signature that is produced by a lighthouse effect where the spinning of strong neutrino emission regions around the rotational axis leads to quasi-periodic modulation in the neutrino signal. Depending on the observer's viewing angle, the time modulation will be clearly detectable in IceCube and the future Hyper-Kamiokande. The GW emission is also anisotropic where the GW signal is emitted, as previously identified, most strongly toward the equator at rotating core-collapse and bounce, and the non-axisymmetric instabilities in the postbounce phase lead to stronger GW emission toward the spin axis. We show that these GW signals can be a target of LIGO-class detectors for a Galactic event. The origin of the postbounce GW emission naturally explains why the peak GW frequency is about twice of the neutrino modulation frequency. We point out that the simultaneous detection of the rotation-induced neutrino and GW signatures could provide a smoking-gun signature of a rapidly rotating proto-neutron star at the birth.

We present results from general-relativistic (GR) three-dimensional (3D) core-collapse simulations with approximate neutrino transport for three non-rotating progenitors (11.2, 15, and 40 Msun) using different nuclear equations of state (EOSs). We find that the combination of progenitor's higher compactness at bounce and the use of softer EOS leads to stronger activity of the standing accretion shock instability (SASI). We confirm previous predications that the SASI produces characteristic time modulations both in neutrino and gravitational-wave (GW) signals. By performing a correlation analysis of the SASI-modulated neutrino and GW signals, we find that the correlation becomes highest when we take into account the time-delay effect due to the advection of material from the neutrino sphere to the proto-neutron star core surface. Our results suggest that the correlation of the neutrino and GW signals, if detected, would provide a new signature of the vigorous SASI activity in the supernova core, which can be hardly seen if neutrino-convection dominates over the SASI.

We suggest the future detection of neutrinos from a Galactic core-collapse supernova can be used to infer the progenitor's inner mass density structure. We present the results from 20 axisymmetric core-collapse supernova simulations performed with progenitors spanning initial masses in the range 11--30Msun, and focus on their connections to the progenitor compactness. The compactness is a measure of the mass density profile of the progenitor core and recent investigations have suggested its salient connections to the outcomes of core collapse. Our simulations confirm a correlation between the neutrinos emitted during the accretion phase and the progenitor's compactness, and that the ratio of observed neutrino events during the first hundreds of milliseconds provides a promising handle on the progenitor's inner structure. Neutrino flavor mixing during the accretion phase remains a large source of uncertainty.

Core-collapse supernovae are the first polluters of heavy elements in the galactic history. As such, it is important to study the nuclear compositions of their ejecta, and understand their dependence on the progenitor structure (e.g., mass, compactness, metallicity). Here, we present a detailed nucleosynthesis study based on two long-term, two-dimensional core-collapse supernova simulations of a 11.2 M$_{\odot}$ and a 17.0 M$_{\odot}$ star. We find that in both models nuclei well beyond the iron group (up to $Z \approx 44$) can be produced, and discuss in detail also the nucleosynthesis of the p-nuclei $^{92,94}$Mo and $^{96,98}$Ru. While we observe the production of $^{92}$Mo and $^{94}$Mo in slightly neutron-rich conditions in both simulations, $^{96,98}$Ru can only be produced efficiently via the $\nu$p-process. Furthermore, the production of Ru in the $\nu$p-process heavily depends on the presence of very proton-rich material in the ejecta. This disentanglement of production mechanisms has interesting consequences when comparing to the abundance ratios between these isotopes in the solar system and in presolar grains.

The gravitational-wave (GW) asteroseismology is a powerful technique for extracting interior information of compact objects. In this work, we focus on spacetime modes, the so-called $w$-modes, of GWs emitted from a proto-neutron star (PNS) in the postbounce phase of core-collapse supernovae. Using results from recent three-dimensional supernova models, we study how to infer the properties of the PNS based on a quasi-normal mode analysis in the context of the GW asteroseismology. We find that the $w_1$-mode frequency multiplied by the PNS radius is expressed as a linear function with respect to the ratio of the PNS mass to the PNS radius. This relation is insensitive to the nuclear equation of state (EOS) employed in this work. Combining with another universal relation of the $f$-mode oscillations, we point out that the time dependent mass-radius relation of the PNS can be obtained by observing both the $f$- and $w_1$-mode GWs simultaneously. Our results suggest that the simultaneous detection of the two modes could provide a new probe into finite-temperature nuclear EOS that predominantly determines the PNS evolution.

We study the effects of collective neutrino oscillations on $\nu p$ process nucleosynthesis in proton-rich neutrino-driven winds by including both the multi-angle $3\times3$ flavor mixing and the nucleosynthesis network calculation. The number flux of energetic electron antineutrinos is raised by collective neutrino oscillations in a $1$D supernova model for $40 M_{\odot}$ progenitor. When the gas temperature decreases down to $\sim2-3\times10^{9}$ K, the increased flux of electron antineutrinos promotes $\nu p$ process more actively, resulting in the enhancement of $p$-nuclei. In the early phase of neutrino-driven wind, blowing at $0.6$ s after core bounce, oscillation effects are prominent in inverted mass hierarchy and $p$-nuclei are synthesized up to $^{106}\mathrm{Cd}$ and $^{108}\mathrm{Cd}$. On the other hand, in the later wind trajectory at $1.1$ s after core bounce, abundances of $p$-nuclei are increased remarkably by $\sim10-10^{4}$ times in normal mass hierarchy and even reaching heavier $p$-nuclei such as $^{124}\mathrm{Xe}$, $^{126}\mathrm{Xe}$ and $^{130}\mathrm{Ba}$. The averaged overproduction factor of $p$-nuclei is dominated by the later wind trajectories. Our results demonstrate that collective neutrino oscillations can strongly influence $\nu p$ process, which indicates that they should be included in the network calculations in order to obtain precise abundances of $p$-nuclei. The conclusions of this paper depend on the difference of initial neutrino parameters between electron and non-electron antineutrino flavors which is large in our case. Further systematic studies on input neutrino physics and wind trajectories are necessary to draw a robust conclusion. However, this finding would help understand the origin of solar-system isotopic abundances of $p$-nuclei such as $^{92,94}\mathrm{Mo}$ and $^{96,98}\mathrm{Ru}$.

The next Galactic supernova is expected to bring great opportunities for the direct detection of gravitational waves (GW), full flavor neutrinos, and multi-wavelength photons. To maximize the science return from such a rare event, it is essential to have established classes of possible situations and preparations for appropriate observations. To this end, we use a long-term numerical simulation of the core-collapse supernova (CCSN) of a 17 solar-mass red supergiant progenitor to self-consistently model the multi-messenger signals expected in GW, neutrino, and electromagnetic messengers. This supernova model takes into account the formation and evolution of a protoneutron star, neutrino-matter interaction, and neutrino transport, all within a two-dimensional shock hydrodynamics simulation. With this, we separately discuss three situations: (i) a CCSN at the Galactic Center, (ii) an extremely nearby CCSN within hundreds of parsecs, and (iii) a CCSN in nearby galaxies within several Mpc. These distance regimes necessitate different strategies for synergistic observations. In a Galactic CCSN, neutrinos provide strategic timing and pointing information. We explore how these in turn deliver an improvement in the sensitivity of GW analyses and help to guarantee observations of early electromagnetic signals. To facilitate the detection of multi-messenger signals of CCSNe in extremely nearby and extragalactic distances, we compile a list of nearby red supergiant candidates and a list of nearby galaxies with their expected CCSN rates. By exploring the sequential multi-messenger signals of a nearby CCSN, we discuss preparations for maximizing successful studies of such an unprecedented stirring event.

We investigated r-process nucleosynthesis in magneto-rotational supernovae, based on a new explosion mechanism induced by the magneto-rotational instability. A series of axisymmetric magneto-hydrodynamical simulations with detailed microphysics including neutrino heating is performed, numerically resolving the magneto-rotational instability. Neutrino-heating dominated explosions, enhanced by magnetic fields, showed mildly neutron-rich ejecta producing nuclei up to $A \sim 130$ (i.e. the weak r-process), while explosion models with stronger magnetic fields reproduce a solar-like r-process pattern. More commonly seen abundance patterns in our models are in between the weak and regular r-process, producing lighter and intermediate mass nuclei. These {\it intermediate r-processes} exhibit a variety of abundance distributions, compatible with several abundance patterns in r-process-enhanced metal-poor stars. The amount of Eu ejecta $\sim 10^{-5} M_\odot$ in magnetically-driven jets agrees with predicted values in the chemical evolution of early galaxies. In contrast, neutrino-heating dominated explosions have a significant amount of Fe ($^{56}{\rm Ni}$) and Zn, comparable to regular supernovae and hypernovae, respectively. These results indicate magneto-rotational supernovae can produce a wide range of heavy nuclei from iron-group to r-process elements, depending on the explosion dynamics.

We present results from fully relativistic three-dimensional core-collapse supernova (CCSN) simulations of a non-rotating 15 Msun star using three different nuclear equations of state (EoSs). From our simulations covering up to ~350 ms after bounce, we show that the development of the standing accretion shock instability (SASI) differs significantly depending on the stiffness of nuclear EoS. Generally, the SASI activity occurs more vigorously in models with softer EoS. By evaluating the gravitational-wave (GW) emission, we find a new GW signature on top of the previously identified one, in which the typical GW frequency increases with time due to an accumulating accretion to the proto-neutron star (PNS). The newly observed quasi-periodic signal appears in the frequency range from ~100 to 200 Hz and persists for ~150 ms before neutrino-driven convection dominates over the SASI. By analyzing the cycle frequency of the SASI sloshing and spiral modes as well as the mass accretion rate to the emission region, we show that the SASI frequency is correlated with the GW frequency. This is because the SASI-induced temporary perturbed mass accretion strike the PNS surface, leading to the quasi-periodic GW emission. Our results show that the GW signal, which could be a smoking-gun signature of the SASI, is within the detection limits of LIGO, advanced Virgo, and KAGRA for Galactic events.

We examine the time evolution of the frequencies of the gravitational wave after the bounce within the framework of relativistic linear perturbation theory using the results of one dimensional numerical simulations of core-collapse supernovae. Protoneutron star models are constructed in such a way that the mass and radius of protoneutron star become equivalent to the results obtained from the numerical simulations. Then, we find that the frequencies of gravitational waves radiating from protoneutron stars strongly depend on the mass and radius of protoneutron stars, but almost independently of the profiles of electron fraction and entropy per baryon inside the star. Additionally, we find that the frequencies of gravitational waves can be characterized by the square root of the average density of protoneutron star irrespectively the progenitor models, which are completely different from the empirical formula for cold neutron stars. The dependence of the spectra on the mass and radius is different from that of the $g$-mode: the oscillations around the surface of protoneutron stars due to the convection and the standing accretion-shock instability. Careful observations of the these modes of gravitational waves can determine the evolution of the mass and radius of protoneutron stars after core-bounce. Furthermore, the expected frequencies of gravitational waves are around a few hundred hertz in the early stage after bounce, which must be a good candidate for the ground-based gravitational wave detectors.

We report results from a series of three-dimensional (3D) rotational core-collapse simulations for $11.2$ and 27 $M_{/odot}$ stars employing neutrino transport scheme by the isotropic diffusion source approximation. By changing the initial strength of rotation systematically, we find a rotation-assisted explosion for the 27$M_{/odot}$ progenitor, which fails in the absence of rotation. The unique feature was not captured in previous two-dimensional (2D) self-consistent rotating models because the growing non-axisymmetric instabilities play a key role. In the rapidly rotating case, strong spiral flows generated by the so-called low $T/|W|$ instability enhance the energy transport from the proto-neutron star (PNS) to the gain region, which makes the shock expansion more energetic. The explosion occurs more strongly in the direction perpendicular to the rotational axis, which is different from previous 2D predictions.

We present a new multi-dimensional radiation-hydrodynamics code for massive stellar core-collapse in full general relativity (GR). Employing an M1 analytical closure scheme, we solve spectral neutrino transport of the radiation energy and momentum based on a truncated moment formalism. Regarding neutrino opacities, we take into account a baseline set in state-of-the-art simulations, in which inelastic neutrinoelectron scattering, thermal neutrino production via pair annihilation and nucleonnucleon bremsstrahlung are included. While the Einstein field equations and the spatial advection terms in the radiation-hydrodynamics equations are evolved explicitly, the source terms due to neutrino-matter interactions and energy shift in the radiation moment equations are integrated implicitly by an iteration method. To verify our code, we first perform a series of standard radiation tests with analytical solutions that include the check of gravitational redshift and Doppler shift. A good agreement in these tests supports the reliability of the GR multi-energy neutrino transport scheme. We then conduct several test simulations of core-collapse, bounce, and shock-stall of a 15Msun star in the Cartesian coordinates and make a detailed comparison with published results. Our code performs quite well to reproduce the results of full-Boltzmann neutrino transport especially before bounce. In the postbounce phase, our code basically performs well, however, there are several differences that are most likely to come from the insufficient spatial resolution in our current 3D-GR models. For clarifying the resolution dependence and extending the code comparison in the late postbounce phase, we discuss that next-generation Exaflops-class supercomputers are at least needed.

By performing neutrino-radiation hydrodynamic simulations in spherical symmetry (1D) and axial symmetry (2D) with different progenitor models by Woosley & Heger (2007) from 12 $M_{\odot}$ to 100 $M_{\odot}$, we find that all 1D runs fail to produce an explosion and several 2D runs succeed. The difference in the shock evolutions for different progenitors can be interpreted by the difference in their mass accretion histories, which are in turn determined by the density structures of progenitors. The mass accretion history has two phases in the majority of the models: the earlier phase in which the mass accretion rate is high and rapidly decreasing and the later phase with a low and almost constant accretion rate. They are separated by the so-called turning point, the origin of which is a change of the accreting layer. We argue that shock revival will most likely occur around the turning point and hence that its location in the $\dot M$-$L_\nu$ plane will be a good measure for the possibility of shock revival: if the turning point lies above the critical curve and the system stays there for a long time, shock revival will obtain. In addition, we develop a phenomenological model to approximately evaluate the trajectories in the $\dot M$-$L_\nu$ plane, which, after calibrating free parameters by a small number of 1D simulations, reproduces the location of the turning point reasonably well by using the initial density structure of progenitor alone. We suggest the application of the phenomenological model to a large collection of progenitors in order to infer without simulations which ones are more likely to explode.

Using predictions from three-dimensional (3D) hydrodynamics simulations of core-collapse supernovae (CCSNe), we present a coherent network analysis to detection, reconstruction, and the source localization of the gravitational-wave (GW) signals. We use the {\tt RIDGE} pipeline for the analysis, in which the network of LIGO Hanford, LIGO Livingston, VIRGO, and KAGRA is considered. By combining with a GW spectrogram analysis, we show that several important hydrodynamics features in the original waveforms persist in the waveforms of the reconstructed signals. The characteristic excess in the spectrograms originates not only from rotating core-collapse, bounce and the subsequent ring down of the proto-neutron star (PNS) as previously identified, but also from the formation of magnetohydrodynamics jets and non-axisymmetric instabilities in the vicinity of the PNS. Regarding the GW signals emitted near at the rotating core bounce, the horizon distance extends up to $\sim$ 18 kpc for the most rapidly rotating 3D model in this work. Following the rotating core bounce, the dominant source of the GW emission shifts to the non-axisymmetric instabilities. The horizon distances extend maximally up to $\sim$ 40 kpc seen from the spin axis. With an increasing number of 3D models trending towards explosion recently, our results suggest that in addition to the best studied GW signals due to rotating core-collapse and bounce, the time is ripe to consider how we can do science from GWs of CCSNe much more seriously than before. Particularly the quasi-periodic signals due to the non-axisymmetric instabilities and the detectability should deserve further investigation to elucidate the inner-working of the rapidly rotating CCSNe.

The r-process nucleosynthesis in core-collapse supernovae (CC-SNe) is studied, with a focus on the explosion scenario induced by rotation and strong magnetic fields. Nucleosynthesis calculations are conducted based on magneto-hydrodynamical explosion models with a wide range of parameters for initial rotation and magnetic fields. The explosion models are classified in two different types: i.e., prompt-magnetic-jet and delayed-magnetic-jet, for which the magnetic fields of proto-neutron stars (PNSs) during collapse and the core-bounce are strong and comparatively moderate, respectively. Following the hydrodynamical trajectories of each explosion model, we confirmed that r-processes successfully occur in the prompt-magnetic-jets, which produce heavy nuclei including actinides. On the other hand, the r-process in the delayed-magnetic-jet is suppressed, which synthesizes only nuclei up to the second peak ($A \sim 130$). Thus, the r-process in the delayed-magnetic-jets could explain only "weak r-process" patterns observed in metal-poor stars rather than the "main r-process", represented by the solar abundances. Our results imply that core-collapse supernovae are possible astronomical sources of heavy r-process elements if their magnetic fields are strong enough, while weaker magnetic explosions may produce "weak r-process" patterns ($A \lesssim 130$). We show the potential importance and necessity of magneto-rotational supernovae for explaining the galactic chemical evolution, as well as abundances of r-process enhanced metal-poor stars. We also examine the effects of the remaining uncertainties in the nature of PNSs due to weak interactions that determine the final neutron-richness of ejecta. Additionally, we briefly discuss radioactive isotope yields in primary jets (e.g., $^{56}$Ni), with relation to several optical observation of SNe and relevant high-energy astronomical phenomena.

We present an overview of two-dimensional (2D) core-collapse supernova simulations employing neutrino transport scheme by the isotropic diffusion source approximation. We study 101 solar-metallicity, 247 ultra metal-poor, and 30 zero-metal progenitors covering zero-age main sequence mass from $10.8 M_{\odot} $ to $75.0 M_{\odot} $. Using the 378 progenitors in total, we systematically investigate how the differences in the structures of these multiple progenitors impact the hydrodynamics evolution. By following a long-term evolution over 1.0 s after bounce, most of the computed models exhibit neutrino-driven revival of the stalled bounce shock at $\sim$ 200 - 800 ms postbounce, leading to the possibility of explosion. Pushing the boundaries of expectations in previous one-dimensional (1D) studies, our results confirm that the compactness parameter $\xi$ that characterizes the structure of the progenitors is also a key in 2D to diagnose the properties of neutrino-driven explosions. Models with high $\xi$ undergo high ram pressure from the accreting matter onto the stalled shock, which affects the subsequent evolution of the shock expansion and the mass of the protoneutron star under the influence of neutrino-driven convection and the standing accretion-shock instability. We show that the accretion luminosity becomes higher for models with high $\xi$, which makes the growth rate of the diagnostic explosion energy higher and the synthesized nickel mass bigger. We find that these explosion characteristics tend to show a monotonic increase as a function of the compactness parameter $\xi$.

Magnetorotational instability (MRI) in a convectively-stable layer around the neutrinosphere is simulated by a three-dimensional model of supernova core. To resolve MRI-unstable modes, a thin layer approximation considering only the radial global stratification is adopted. Our intriguing finding is that the convectively-stable layer around the neutrinosphere becomes fully-turbulent due to the MRI and its nonlinear penetration into the strongly-stratified MRI-stable region. The intensity of the MRI-driven turbulence increases with magnetic flux threading the core, but is limited by a free energy stored in the differential rotation. The turbulent neutrinosphere is a natural consequence of rotating core-collapse and could exert a positive impact on the supernova mechanism.

We study the multi-dimensional properties of neutrino transfer inside supernova cores by solving the Boltzmann equations for neutrino distribution functions in genuinely six dimensional (6D) phase space. Adopting representative snapshots of the post-bounce core from other supernova simulations in three dimensions, we solve the temporal evolutions to stationary states of neutrino distribution functions by our Boltzmann solver. Taking advantage of the multi-angle and multi-energy feature realized by the S$_n$ method in our code, we reveal the genuine characteristics of spatially three dimensional (3D) neutrino transfer such as non-radial fluxes and non-diagonal Eddington tensors. In addition, we assess the ray-by-ray approximation, turning off the lateral-transport terms in our code. We demonstrate that the ray-by-ray approximation tends to propagate fluctuations in thermodynamical states around the neutrino-sphere along each radial ray and overestimate the variations between the neutrino distributions on different radial rays. We find that the difference in the densities and fluxes of neutrinos between the ray-by-ray approximation and the full Boltzmann transport becomes ~20%, which is also the case for the local heating rate, whereas the volume-integrated heating rate in the Boltzmann transport is found to be only slightly larger (~2%) than the counterpart in the ray-by-ray approximation due to cancellation among different rays. These results suggest that we had better assess carefully the possible influences of various approximations in the neutrino transfer employed in the current simulations on supernova dynamics. Detailed information on the angle and energy moments of neutrino distribution functions will be profitable for the future development of numerical methods in neutrino-radiation hydrodynamics.

Mapping supernovae to their progenitors is fundamental to understanding the collapse of massive stars. We investigate the red supergiant problem, which concerns why red supergiants with masses $\sim16$-$30 M_\odot$ have not been identified as progenitors of Type IIP supernovae, and the supernova rate problem, which concerns why the observed cosmic supernova rate is smaller than the observed cosmic star formation rate. We find key physics to solving these in the compactness parameter, which characterizes the density structure of the progenitor. If massive stars with compactness above $\xi_{2.5} \sim 0.2$ fail to produce canonical supernovae, (i) stars in the mass range $16$-$30 M_\odot$ populate an island of stars that have high $\xi_{2.5}$ and do not produce canonical supernovae, and (ii) the fraction of such stars is consistent with the missing fraction of supernovae relative to star formation. We support this scenario with a series of two- and three-dimensional radiation hydrodynamics core-collapse simulations. Using more than 300 progenitors covering initial masses $10.8$-$75 M_\odot$ and three initial metallicities, we show that high compactness is conducive to failed explosions. We then argue that a critical compactness of $\sim 0.2$ as the divide between successful and failed explosions is consistent with state-of-the-art three-dimensional core-collapse simulations. Our study implies that numerical simulations of core collapse need not produce robust explosions in a significant fraction of compact massive star initial conditions.

We perform a series of simplified numerical experiments to explore how rotation impacts on the three-dimensional (3D) hydrodynamics of core-collapse supernovae. For the sake of our systematic study, we employ a light-bulb scheme to trigger explosions and a three-flavor neutrino leakage scheme to treat deleptonization effects and neutrino losses from proto-neutron star interior. Using a 15 solar mass progenitor, we compute thirty models in 3D with a wide variety of initial angular momentum and light-bulb neutrino luminosity. We find that the rotation can help onset of neutrino-driven explosions for the models in which the initial angular momentum is matched to that obtained in recent stellar evolutionary calculations (0.3-3 rad/s at the center). For the models with larger initial angular momentum, the shock surface deforms to be more oblate due to larger centrifugal force. This makes not only a gain region more concentrated around the equatorial plane, but also the mass in the gain region bigger. As a result, buoyant bubbles tend to be coherently formed and rise in the equatorial region, which pushes the revived shock ever larger radii until a global explosion is triggered. We find that these are the main reasons that the preferred direction of explosion in 3D rotating models is often perpendicular to the spin axis, which is in sharp contrast to the polar explosions around the axis that was obtained in previous 2D simulations.

We study properties of gravitational waves (GWs) from rotating core-collapse of a 15M_odot star by performing three-dimensional general-relativistic hydrodynamic simulations with an approximate neutrino transport. By parametrically changing the precollapse angular momentum, we focus on the effects of rotation on the GW signatures in the early postbounce evolution. Regarding three-flavor neutrino transport, we solve the energy-averaged set of radiation energy and momentum based on the Thorne's momentum formalism. In addition to the gravitational quadrupole radiation from matter motions, we take into account GWs from anisotropic neutrino emission. With these computations, our results present several supporting evidences for the previous anticipation that non-axisymmetric instabilities play an essential role in determining the postbounce GW signatures. During prompt convection, we find that the waveforms show narrow-band and highly quasi-periodic signals which persist until the end of simulations. We point out that such feature reflects the growth of the one-armed spiral modes. The typical frequency of the quasi-periodic waveforms can be well explained by the propagating acoustic waves between the stalled shock and the rotating proto-neutron star surface, which suggests the appearance of the standing-accretion-shock instability. Although the GW signals exhibit strong variability between the two polarizations and different observer directions, they are within the detection limits of next generation detectors such as by KAGRA and Advanced LIGO, if the source with sufficient angular momentum is located in our Galaxy.

We revisit potential impacts of nuclear burning on the onset of the neutrino-driven explosions of core-collapse supernovae. By changing the neutrino luminosity and its decay time to obtain parametric explosions in one-(1D) and two-dimensional (2D) models with or without a 13-isotope alpha network, we study how the inclusion of nuclear burning could affect the postbounce dynamics for four progenitor models; three for 15.0 Msun stars, one for an 11.2 Msun star. We find that the energy supply due to nuclear burning of infalling material behind the shock can energize the shock expansion especially for models that produce only marginal explosions in the absence of nuclear burning. These models are energized by nuclear energy deposition when the shock front passes through the silicon-rich layer and/or later it touches the oxygen-rich layer. Depending on the neutrino luminosity and its decay time, a diagnostic energy of explosion increases up to a few times 10^50 erg for models with nuclear burning compared to the corresponding models without. We point out that these features are most remarkable for the Limongi-Chieffi progenitor in both 1D and 2D, because the progenitor model possesses a massive oxygen layer with its inner-edge radius being smallest among the employed progenitors, so that the shock can touch the rich fuel on a shorter timescale after bounce. The energy difference is generally smaller (~0.1-0.2 times 10^51 erg) in 2D than in 1D (at most ~0.6 times 10^51 erg). This is because neutrino-driven convection and the shock instability in 2D models enhance the neutrino heating efficiency, which makes the contribution of nuclear burning relatively smaller compared to 1D models. Considering uncertainties in progenitor models, our results indicate that nuclear burning should remain as one of the important ingredients to foster the onset of neutrino-driven explosions.

We present numerical results on two- (2D) and three-dimensional (3D) hydrodynamic core-collapse simulations of an 11.2$M_\odot$ star. By changing numerical resolutions and seed perturbations systematically, we study how the postbounce dynamics is different in 2D and 3D. The calculations were performed with an energy-dependent treatment of the neutrino transport based on the isotropic diffusion source approximation scheme, which we have updated to achieve a very high computational efficiency. All the computed models in this work including nine 3D models and fifteen 2D models exhibit the revival of the stalled bounce shock, leading to the possibility of explosion. All of them are driven by the neutrino-heating mechanism, which is fostered by neutrino-driven convection and the standing-accretion-shock instability (SASI). Reflecting the stochastic nature of multi-dimensional (multi-D) neutrino-driven explosions, the blast morphology changes from models to models. However, we find that the final fate of the multi-D models whether an explosion is obtained or not, is little affected by the explosion stochasticity. In agreement with some previous studies, higher numerical resolutions lead to slower onset of the shock revival in both 3D and 2D. Based on the self-consistent supernova models leading to the possibility of explosions, our results systematically show that the revived shock expands more energetically in 2D than in 3D.

By implementing widely-used equations of state (EOS) from Lattimer & Swesty (LS) and H. Shen et al. (SHEN) in core-collapse supernova simulations, we explore possible impacts of these EOS on the post-bounce dynamics prior to the onset of neutrino-driven explosions. Our spherically symmetric (1D) and axially symmetric (2D) models are based on neutrino radiation hydrodynamics including spectral transport, which is solved by the isotropic diffusion source approximation. We confirm that in 1D simulations neutrino-driven explosions cannot be obtained for any of the employed EOS. Impacts of the EOS on the post-bounce hydrodynamics are more clearly visible in 2D simulations. In 2D models of a 15 M_sun progenitor using the LS EOS, the stalled bounce shock expands to increasingly larger radii, which is not the case using the SHEN EOS. Keeping in mind that the omission of the energy drain by heavy-lepton neutrinos in the present scheme could facilitate explosions, we find that 2D models of an 11.2 M_sun progenitor produce neutrino-driven explosions for all the EOS under investigation. Models using the LS EOS are slightly more energetic compared to those with the SHEN EOS. The more efficient neutrino heating in the LS models coincides with a higher electron antineutrino luminosity and a larger mass that is enclosed within the gain region. The models based on the LS EOS also show a more vigorous and aspherical downflow of accreting matter to the surface of the protoneutron star (PNS). The accretion pattern is essential for the production and strength of outgoing pressure waves, that can push in turn the shock to larger radii and provide more favorable conditions for the explosion. [abbreviated]

Bearing in mind the application to core-collapse supernovae, we study nonlinear properties of the magneto-rotational instability (MRI) by means of three- dimensional simulations in the framework of a local shearing box approximation. By changing systematically the shear rates that symbolize the degree of differential rotation in nascent proto-neutron stars (PNSs), we derive a scaling relation between the turbulent stress sustained by the MRI and the shear- vorticity ratio. Our parametric survey shows a power-law scaling between the turbulent stress ($<< w_{\rm tot}>>$) and the shear- vorticity ratio ($g_q$) as $<<w_{\rm tot}>> \propto g_q^{\delta}$ with its index $\delta \sim 0.5$. The MRI-amplified magnetic energy has a similar scaling relative to the turbulent stress, while the Maxwell stress has slightly smaller power-law index ($\sim 0.36$). By modeling the effect of viscous heating rates due to the MRI turbulence, we show that the stronger magnetic fields or the larger shear rates initially imposed lead to the higher dissipation rates. For a rapidly rotating PNS with the spin period in milliseconds and with strong magnetic fields of $10^{15}$ G, the energy dissipation rate is estimated to exceed $10^{51} {\rm erg\ sec^{-1}}$. Our results suggest that the conventional magnetohydrodynamic (MHD) mechanism of core-collapse supernovae is likely to be affected by the MRI-driven turbulence, which we speculate, on one hand, could harm the MHD-driven explosions due to the dissipation of the shear rotational energy at the PNS surface, on the other hand the energy deposition there might be potentially favorable for the working of the neutrino-heating mechanism.

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 present results from the first generation of multi-dimensional hydrodynamic core-collapse simulations in full general relativity (GR) that include an approximate treatment of neutrino transport. Using a M1 closure scheme with an analytic variable Eddington factor, we solve the energy-independent set of radiation energy and momentum based on the Thorne's momentum formalism. To simplify the source terms of the transport equations, a methodology of multiflavour neutrino leakage scheme is partly employed. Our newly developed code is designed to evolve the Einstein field equation together with the GR radiation hydrodynamic equations. We follow the dynamics starting from the onset of gravitational core-collapse of a 15 $M_{\odot}$ star, through bounce, up to about 100 ms postbounce in this study to study how the spacial multi-dimensionality and GR would affect the dynamics in the early postbounce phase. Our 3D results support the anticipation in previous 1D results that the neutrino luminosity and average neutrino energy of any neutrino flavor in the postbounce phase increase when switching from SR to GR hydrodynamics. This is because the deeper gravitational well of GR produces more compact core structures, and thus hotter neutrino spheres at smaller radii. By analyzing the residency timescale to the neutrino-heating timescale in the gain region, we show that the criterion to initiate neutrino-driven explosions can be most easily satisfied in 3D models, irrespective of SR or GR hydrodynamics. Our results suggest that the combination of GR and 3D hydrodynamics provides the most favorable condition to drive a robust neutrino-driven explosion.

By performing two-dimensional special relativistic (SR) magnetohydrodynamic simulations, we study possible signatures of gravitational waves (GWs) in the context of the collapsar model for long-duration gamma-ray bursts. In our SR simulations, the central black hole is treated as an absorbing boundary. By doing so, we focus on the GWs generated by asphericities in neutrino emission and matter motions in the vicinity of the hyperaccreting disks. We compute nine models by adding initial angular momenta and magnetic fields parametrically to a precollapse core of a $35 M_{\odot}$ progenitor star. As for the microphysics, a realistic equation of state is employed and the neutrino cooling is taken into account via a multiflavor neutrino leakage scheme. To accurately estimate GWs produced by anisotropic neutrino emission, we perform a ray-tracing analysis in general relativity by a post-processing procedure. By employing a stress formula that includes contributions both from magnetic fields and special relativistic corrections, we study also the effects of magnetic fields on the gravitational waveforms. We find that the GW amplitudes from anisotropic neutrino emission show a monotonic increase with time, whose amplitudes are much larger than those from matter motions of the accreting material. We show that the increasing trend of the neutrino GWs stems from the excess of neutrino emission in the direction near parallel to the spin axis illuminated from the hyperaccreting disks. We point out that a recently proposed future space-based interferometer like Fabry-Perot type DECIGO would permit the detection of these GW signals within $\approx$ 100 Mpc.

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 present numerical results on three-dimensional (3D) hydrodynamic core-collapse simulations of an $11.2 M_{\odot}$ star. By comparing one-(1D) and two-dimensional(2D) results with those of 3D, we study how the increasing spacial multi-dimensionality affects the postbounce supernova dynamics. The calculations were performed with an energy-dependent treatment of the neutrino transport that is solved by the isotropic diffusion source approximation scheme. By performing a tracer-particle analysis, we show that the maximum residency time of material in the gain region is shown to be longer for 3D due to non-axisymmetric flow motions than 2D, which is one of advantageous aspects of 3D models to obtain neutrino-driven explosions. Our results show that convective matter motions below the gain radius become much more violent in 3D than 2D, making the neutrino luminosity larger for 3D. Nevertheless the emitted neutrino energies are made smaller due to the enhanced cooling. Our results indicate whether these advantages for driving 3D explosions could or could not overwhelm the disadvantages is sensitive to the employed numerical resolutions. An encouraging finding is that the shock expansion tends to become more energetic for models with finer resolutions. To draw a robust conclusion, 3D simulations with much more higher numerical resolutions and also with more advanced treatment of neutrino transport as well as of gravity is needed, which could be hopefully practicable by utilizing forthcoming Petaflops-class supercomputers.

By performing a series of two-dimensional, special relativistic magnetohydrodynamic (MHD) simulations, we study signatures of gravitational waves (GWs) in the magnetohydrodynamically-driven core-collapse supernovae. In order to extract the gravitational waveforms, we present a stress formula including contributions both from magnetic fields and special relativistic corrections. By changing the precollapse magnetic fields and initial angular momentum distributions parametrically, we compute twelve models. As for the microphysics, a realistic equation of state is employed and the neutrino cooling is taken into account via a multiflavor neutrino leakage scheme. With these computations, we find that the total GW amplitudes show a monotonic increase after bounce for models with a strong precollapse magnetic field ($10^{12}$G) also with a rapid rotation imposed. We show that this trend stems both from the kinetic contribution of MHD outflows with large radial velocities and also from the magnetic contribution dominated by the toroidal magnetic fields that predominantly trigger MHD explosions. For models with weaker initial magnetic fields, the total GW amplitudes after bounce stay almost zero, because the contribution from the magnetic fields cancels with the one from the hydrodynamic counterpart. These features can be clearly understood with a careful analysis on the explosion dynamics. We point out that the GW signals with the increasing trend, possibly visible to the next-generation detectors for a Galactic supernova, would be associated with MHD explosions with the explosion energies exceeding $10^{51}$ erg.

By performing a series of one- and two-dimensional (1-, 2D) hydrodynamic simulations with spectral neutrino transport, we study possible impacts of collective neutrino oscillations on the dynamics of core-collapse supernovae. To model the spectral swapping which is one of the possible outcome of the collective neutrino oscillations, we parametrize the onset time when the spectral swap begins, the radius where the spectral swap occurs, and the threshold energy above which the spectral interchange between heavy-lepton neutrinos and electron/anti-electron neutrinos takes place, respectively. By doing so, we systematically study how the neutrino heating enhanced by the spectral swapping could affect the shock evolution as well as the matter ejection. We also investigate the progenitor dependence using a suite of progenitor models (13, 15, 20, and 25 $M_\odot$). We find that there is a critical heating rate induced by the spectral swapping to trigger explosions, which significantly differs between the progenitors. The critical heating rate is generally smaller for 2D than 1D due to the multidimensionality that enhances the neutrino heating efficiency. For the progenitors employed in this paper, the final remnant masses are estimated to range in 1.1-1.5$M_\odot$. For our 2D model of the $15M_\odot$ progenitor, we find a set of the oscillation parameters that could account for strong supernova explosions ($\sim 10^{51}$ erg), simultaneously leaving behind the remnant mass close to $\sim 1.4 M_\odot$.

By performing axisymmetric hydrodynamic simulations of core-collapse supernovae with spectral neutrino transport based on the isotropic diffusion source approximation scheme, we support the assumption that the neutrino-heating mechanism aided by the standing accretion shock instability and convection can initiate an explosion of a 13 $M_{\odot}$ star. Our results show that bipolar explosions are more likely to be associated with models which include rotation. We point out that models, which form a north-south symmetric bipolar explosion, can lead to larger explosion energies than for the corresponding unipolar explosions.

Bearing in mind the application to the collapsar models of gamma-ray bursts (GRBs), we develop a numerical scheme and code for estimating the deposition of energy and momentum due to the neutrino pair annihilation ($\nu + {\bar \nu} \rightarrow e^{-} + e^{+}$) in the vicinity of accretion tori around a Kerr black hole. Our code is designed to solve the general relativistic neutrino transfer by a ray-tracing method. To solve the collisional Boltzmann equation in curved spacetime, we numerically integrate the so-called rendering equation along the null geodesics. For the neutrino opacity, the charged-current $\beta$-processes are taken into account, which are dominant in the vicinity of the accretion tori. The numerical accuracy of the developed code is certificated by several tests, in which we show comparisons with the corresponding analytic solutions. Based on the hydrodynamical data in our collapsar simulation, we estimate the annihilation rates in a post-processing manner. Increasing the Kerr parameter from 0 to 1, it is found that the general relativistic effect can increase the local energy deposition rate by about one order of magnitude, and the net energy deposition rate by several tens of percents. After the accretion disk settles into a stationary state (typically later than $\sim 9$ s from the onset of gravitational collapse), we point out that the neutrino-heating timescale in the vicinity of the polar funnel region can be shorter than the dynamical timescale. Our results suggest the neutrino pair annihilation has a potential importance equal to the conventional magnetohydrodynamic mechanism for igniting the GRB fireballs.

We develop a numerical scheme and code for estimating the energy and momentum transfer via neutrino pair annihilation ($\nu + {\bar \nu} \to e^{-}+ e^{+}$), bearing in mind the application to the collapsar models of gamma-ray bursts (GRBs). To calculate the neutrino flux illuminated from the accretion disk, we perform a ray-tracing calculation in the framework of special relativity. The numerical accuracy of the developed code is certificated by several tests, in which we show comparisons with the corresponding analytical solutions. Using hydrodynamical data in our collapsar simulation, we estimate the annihilation rates in a post-processing manner. We show that the neutrino energy deposition and momentum transfers are strongest near the inner edge of the accretion disk. The beaming effects of special relativity are found to change the annihilation rates by several factors in the polar funnel region. After the accretion disk settles into a stationary state (typically later than $\sim 9$ s from the onset of gravitational collapse), we find that the neutrino-heating timescale in the vicinity of the polar funnel ($\lesssim 80$ km) can become shorter than the hydrodynamical timescale, indicating that the neutrino-heated outflows can be launched there. We point out that the momentum transfer can play as important role as the energy deposition for the efficient acceleration of neutrino-driven outflows. Our results suggest that the neutrino pair annihilation has a potential importance equal to the conventional magnetohydrodynamic mechanism for igniting the GRB fireballs.

We present our numerical results of two-dimensional hydrodynamic (HD) simulations and magnetohydrodynamic (MHD) simulations of the collapse of rotating massive stars in light of the collapsar model of gamma-ray bursts (GRBs). Pushed by recent evolution calculations of GRB progenitors, we focus on lower angular momentum of the central core than the ones taken mostly in previous studies. By performing special relativistic simulations including both realistic equation of state and neutrino cooling, we follow a long-term evolution of the slowly rotating collapsars up to 10 s, accompanied by the formation of jets and accretion disks. We find such outfows can be launched both by MHD process and neutrino process. We investigeate the properties of these jets whether it can become GRBs or remains primary weak outflow.

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.

We present our numerical results of two-dimensional magnetohydrodynamic (MHD) simulations of the collapse of rotating massive stars in light of the collapsar model of gamma-ray bursts (GRBs). Pushed by recent evolution calculations of GRB progenitors, we focus on lower angular momentum of the central core than the ones taken mostly in previous studies. By performing special relativistic simulations including both realistic equation of state and neutrino coolings, we follow a unprecedentedly long-term evolution of the slowly rotating collapsars up to $\sim$ 10 s, accompanied by the formation of jets and accretion disks. Our results show that for the GRB progenitors to function as collapsars, there is a critical initial angular momentum, below which matter is quickly swallowed to the central objects, no accretion disks and no MHD outflows are formed. When larger than the criteria, we find the launch of the MHD jets in the following two ways. For models with stronger initial magnetic fields, the magnetic pressure amplified inside the accretion disk can drive the MHD outflows, which makes the strong magnetic explosions like a 'magnetic tower' (type II). For the models with weaker initial magnetic fields, the magnetic tower stalls first and the subsequent MHD outflows are produced by some eventual inflows of the accreting material from the equator to the polar regions (type I). Regardless of type I or II, the jets can attain only mildly relativistic speeds with the explosion energy less than $10^{49} \erg$, which could possibly be related to the X-ray flashes. Due to high opacity for neutrinos inside the disk, we find that the luminosities of $\nu_e$ and $\bar{\nu}_e$ become almost comparable, which is advantageous for making the energy deposition rate larger.

We perform a series of two-dimensional magnetohydrodynamic core-collapse simulations of rapidly rotating and strongly magnetized massive stars. To study the properties of magnetic explosions for a longer time stretch of postbounce evolution, we develop a new code under the framework of special relativity including a realistic equation of state with a multiflavor neutrino leakage scheme. Our results show the generation of the magnetically-dominated jets in the two ways. One is launched just after the core-bounce in a prompt way and another is launched at $ \sim 100 $ ms after the stall of the prompt shock. We find that the shock-revival occurs when the magnetic pressure becomes strong, due to the field wrapping, enough to overwhelm the ram pressure of the accreting matter. The critical toroidal magnetic fields for the magnetic shock-revival are found to be universal of $\sim 10^{15}\mathrm{G}$ behind the jets. We point out that the time difference before the shock-revival has a strong correlation with the explosions energies. Our results suggest that the magnetically dominated jets are accompanied by the formation of the magnetars. Since the jets are mildly relativistic, we speculate that they might be the origin of some observed X-ray flashes.

We study the effects of rotation on the neutrino emission from Population III (Pop III) stars by performing a series of two-dimensional rotational collapse simulations of Pop III stellar cores. Our results show that rotation enhances the neutrino luminosities and the average energies of emitted neutrinos. This is because the thermalized inner core, which is the dominant neutrino source from Pop III stars, can be enlarged, due to rotational flattening, enough to extend the inner core outside the neutrinospheres. This is in sharp contrast to the case of spherical collapse, in which the case of inner core shrinks deeper inside the neutrinospheres before black hole formation, which hinders the efficient neutrino emission. In the case of rotational core-collapse, the emitted neutrino energies are found to become larger in the vicinity near the pole than the ones near the equatorial plane. These factors make the emergent neutrino spectrum broader and harder than the spherical collapse case. By computing the overall neutrino signals produced by the ensemble of individual rotating Pop III stars, we find that the amplitudes of the relic neutrinos, depending on their star formation rates, can dominate over the contributions from ordinary core-collapse supernovae below a few MeV. A detection of this signal could be an important tool to probe star formation history in the early universe.

We estimate the gravitational wave (GW) background from Population III (Pop III) stars using the results from our hydrodynamic simulations (Suwa et al. 2007). We calculate gravitational waveforms from matter motions and anisotropic neutrino emissions for single Pop III stars. We find that the GW amplitudes from matter motions are dominant until just after bounce, but those from neutrinos dominate later on at frequencies below $\sim 10$ Hz in the GW spectrum. Computing the overall signal produced by the ensemble of such Pop III stars, we find that the resultant density parameter of the GW background peaks at the amplitude of $\Omega_\mathrm{GW}\sim 10^{-10}$ in the frequency interval $\sim 1-10$ Hz. We show that such signals, depending on the formation rate of Pop III stars, can be within the detection limits of future planned interferometers such as DECIGO and BBO in the frequency interval of $\sim 0.1-1$ Hz. Our results suggest that the detection of the GW background from Pop III stars can be an important tool to supply the information about the star formation history in the early universe.

Nuclear ``pasta'', nonspherical nuclei in dense matter, is predicted to occur in collapsing supernova cores. We show how pasta phases affect the neutrino transport cross section via weak neutral current using several nuclear models. This is the first calculation of the neutrino opacity of the phases with rod-like and slab-like nuclei taking account of finite temperature effects, which are well described by the quantum molecular dynamics. We also show that pasta phases can occupy 10--20% of the mass of supernova cores in the later stage of the collapse.

We perform a series of two-dimensional magnetorotational core-collapse simulations of Pop III stars. Changing the initial distributions of rotation and magnetic fields prior to collapse in a parametric manner, we compute 19 models. By so doing, we systematically investigate how rotation and magnetic fields affect the collapse dynamics and explore how the properties of the black-hole formations and neutrino emissions could be affected. As for the microphysics, we employ a realistic equation of state and approximate the neutrino transfer by a multiflavour leakage scheme. With these computations, we find that the jet-like explosions are obtained by the magnetodriven shock waves if the initial magnetic field is as large as $10^{12}$ G. We point out that the black-hole masses at the formation decrease with the initial field strength, on the other hand, increase with the initial rotation rates. As for the neutrino properties, we point out that the degree of the differential rotation plays an important role to determine which species of the neutrino luminosity is more dominant than the others. Furthermore, we find that the stronger magnetic fields make the peak neutrino luminosities smaller, because the magnetic pressure acts to halt the collapse in the central regions, leading to the suppression of the releasable gravitational binding energies.

Two-dimensional magnetohydrodynamic simulations are performed using the ZEUS-2D code to investigate the dynamics of a collapsar that generates a GRB jet, taking account of realistic equation of state, neutrino cooling and heating processes, magnetic fields, and gravitational force from the central black hole and self gravity. It is found that neutrino heating processes are not so efficient to launch a jet in this study. It is also found that a jet is launched mainly by B_\phi fields that are amplified by the winding-up effect. However, since the ratio of total energy relative to the rest mass energy in the jet is not so high as several hundred, we conclude that the jets seen in this study are not be a GRB jet. This result suggests that general relativistic effects, which are not included in this study, will be important to generate a GRB jet. Also, the accretion disk with magnetic fields may still play an important role to launch a GRB jet, although a simulation for much longer physical time (\sim 10-100 s) is required to confirm this effect. It is shown that considerable amount of 56Ni is synthesized in the accretion disk. Thus there will be a possibility for the accretion disk to supply sufficient amount of 56Ni required to explain the luminosity of a hypernova. Also, it is shown that neutron-rich matter due to electron captures with high entropy per baryon is ejected along the polar axis. Moreover, it is found that the electron fraction becomes larger than 0.5 around the polar axis near the black hole by \nu_e capture at the region. Thus there will be a possibility that r-process and r/p-process nucleosynthesis occur at these regions. Finally, much neutrons will be ejected from the jet, which suggests that signals from the neutron decays may be observed as the delayed bump of afterglow or gamma-rays.

We perform a series of two-dimensional magnetohydrodynamic simulations of the rotational core-collapse of a magnetized massive star. We employ a realistic equation of state and take into account the neutrino cooling by the so-called leakage scheme. In this study we systematically investigate how the strong magnetic field and the rapid rotation affect the propagation of the shock waves. Our results show that in the case of the strong initial poloidal magnetic field, the toroidal magnetic field amplified by the differential rotation, becomes strong enough to generate a tightly collimated shock wave along the rotational axis. On the other hand, in the case of the weak initial magnetic field, although the differential rotation amplifies toroidal magnetic field over the long rotational period, the launched shock wave is weak and the shape of it becomes wider. The former case is expected to be accompanied by the formation of the so-called magnetar. Our models with rapid rotation and strong magnetic field can create a nozzle formed by the collimated shock wave. This might be the analogous situation of the collapsar that is plausible for the central engine of the Gamma-Ray Bursts.