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 ∼ 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 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 numerical results on two-(2D) and three-dimensional (3D) hydrodynamic core-collapse simulations of an 11.2M ⊙ 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.

We present numerical results on three-dimensional (3D) hydrodynamic core-collapse simulations of an 11.2M ⊙ 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. In agreement with previous study, our 1D model does not produce explosions for the 11.2 M ⊙ star, while the neutrino-driven revival of the stalled bounce shock is obtained both in the 2D and 3D models. The standing accretion-shock instability (SASI) is observed in the 3D models, in which the dominant mode of the SASI is bipolar (ℓ = 2) with its saturation amplitudes being slightly smaller than 2D. By performing a tracer-particle analysis, we show that the maximum residency time of material in the gain region becomes longer in 3D due to non-axisymmetric flow motions than in 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 in 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 are needed, which could be hopefully practicable by utilizing forthcoming Petaflops-class supercomputers.

We present results from fully relativistic three-dimensional core-collapse supernova simulations of a non-rotating  M 15 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 SASIinduced temporary perturbed mass accretion strikes 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.

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 (SASI) 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 that include rotation. We point out that models that form a north–south symmetric bipolar explosion can lead to larger explosion energies than the corresponding unipolar explosions can.

We investigated r-process nucleosynthesis in magneto-rotational supernovae, based on a new explosion mechanism induced by the magneto-rotational instability (MRI). A series of axisymmetric magnetohydrodynamical simulations with detailed microphysics including neutrino heating is performed, numerically resolving the MRI. Neutrino-heating dominated explosions, enhanced by magnetic fields, showed mildly neutronrich ejecta producing nuclei up toà 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 intermediate r-processes exhibit a variety of abundance distributions, compatible with several abundance patterns in r-processenhanced metal-poor stars. The amount of Eu ejectain 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 ( Ni 56 ) 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.

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 M 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 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 ∼ 100 ms after the stall of the prompt shock. We find that the shock-revival occurs when the magnetic pressure becomes strong enough, due to the field wrapping, 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 ∼ 10 15 G behind the jets. We point out that the time difference before the shock-revival has a strong correlation with the explosion 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.