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 carried out 2D-axisymmetric MHD simulations of core-collapse supernovae for rapidly-rotating magnetized progenitors. By changing both the strength of the magnetic field and the spatial resolution, the evolution of the magnetorotational instability (MRI) and its impacts upon the dynamics are investigated. We found that the MRI greatly amplifies the seed magnetic fields in the regime where not the Alfv\'en mode but the buoyant mode plays a primary role in the exponential growth phase. The MRI indeed has a powerful impact on the supernova dynamics. It makes the shock expansion faster and the explosion more energetic, with some models being accompanied by the collimated-jet formations. These effects, however, are not made by the magnetic pressure except for the collimated-jet formations. The angular momentum transfer induced by the MRI causes the expansion of the heating region, by which the accreting matter gain an additional time to be heated by neutrinos. The MRI also drifts low-$Y_p$ matter from the deep inside of the core to the heating region, which makes the net neutrino heating rate larger by the reduction of the cooling due to the electron capture. These two effects enhance the efficiency of the neutrino heating, which is found to be the key to boost the explosion. Indeed we found that our models explode far more weakly when the net neutrino heating is switched off. The contribution of the neutrino heating to the explosion energy could reach 60\% even in the case of strongest magnetic field in the current simulations.

The majority of massive stars are formed in binary systems. It is hence reasonable to expect that most core-collapse supernovae (CCSNe) take place in binaries and the existence of a companion star may leave some imprints in observed features. Having this in mind, we have conducted two-dimensional hydrodynamical simulations of the collisions of CCSNe ejecta with the companion star in an almost-equal-mass ($\sim 10M_\odot$) binary to find out possible consequences of such events. In particular we pay attention to the amount of mass removed and its dependence on the binary separation. In contrast to the previous surmise, we find that the companion mass is stripped not by momentum transfer but by shock heating. Up to $25\%$ of the original mass can be removed for the closest separations and the removed mass decreases as $M_{ub} \propto a^{-4.3}$ with the binary separation $a$. By performing some experimental computations with artificially-modified densities of incident ejecta, we show that if the velocity of ejecta is fixed, the density of incident ejecta is the single important parameter that actually determines the removed mass as $M_{ub} \propto \rho_{ej} ^{1.4} $. On the other hand, another set of simulations with modified velocities of incident ejecta demonstrate that the strength of the forward shock, which heats up the stellar material and causes the mass loss of the companion star, is actually the key parameter for the removed mass.

We investigated the impacts of magnetorotational instability (MRI) on the dynamics of weakly magnetized, rapidly rotating core-collapse by conducting high resolution MHD simulations in axisymmetry with simplified neutrino transfer. We found that an initially sub-magnetar class magnetic field is drastically amplified by MRI and substantially affects the dynamics thereafter. Although the magnetic pressure is not strong enough to eject matter, the amplified magnetic field efficiently transfers angular momentum from higher to lower latitudes, which causes the expansion of the heating region at low latitudes due to the extra centrifugal force. This then enhance the efficiency of neutrino heating and eventually leads to neutrino-driven explosion. This is a new scenario of core-collapse supernovae that has never been demonstrated by numerical simulations so far.

We performed the first numerical simulations of magnetorotational instability from a sub-magnetar-class seed magnetic field in core collapse supernovae. As a result of axisymmetric ideal MHD simulations, we found that the magnetic field is greatly amplified to magnetar-class strength. In saturation phase, a substantial part of the core is dominated by turbulence, and the magnetic field possesses dominant large scale components, comparable to the size of the proto-neutron star. A pattern of coherent chanel flows, which generally appears during exponential growth phase in previous local simulations, is not observed in our global simulations. While the approximate convergence in the exponential growth rate is attained by increasing spatial resolution, that of the saturation magnetic field is not achieved due to still large numerical diffusion. Although the effect of magnetic field on the dynamics is found to be mild, a simulation with a high-enough resolution might result in a larger impact.

We studied roles of a turbulent resistivity in the core-collapse of a strongly magnetized massive star, carrying out 2D-resistive-MHD simulations. The three cases with different initial strengths of magnetic field and rotation are investigated; 1. strongly magnetized rotating core; 2.moderately magnetized rotating core; 3. very strongly magnetized non-rotating core. In each case, both an ideal-MHD model and resistive-MHD models are computed. As a result of computations, each model shows a matter eruption helped by a magnetic acceleration (and also by a centrifugal acceleration in the rotating cases). We found that a resistivity attenuates the explosion in case~1 and 2, while it enhances the explosion in case~3. We also found that in the rotating cases, main mechanisms for the amplification of a magnetic field in the post-bounce phase are an outward advection of magnetic field and a winding of poloidal magnetic field-lines by differential rotation, which are somewhat dampened down with the presence of a resistivity. Although the magnetorotational instability seems to occur in the rotating models, it will play only a minor role in a magnetic field amplification. Another impact of resistivity is that on the aspect ratio. In the rotating cases, a large aspect ratio of the ejected matters, $> 2.5$, attained in a ideal-MHD model is reduced to some extent in a resistive model. These results indicate that a resistivity possibly plays an important role in the dynamics of strongly magnetized supernovae.

A series of numerical simulations on magnetorotational core-collapse supernovae are carried out. Dipole-like configurations which are offset northward are assumed for the initially strong magnetic fields together with rapid differential rotations. Aims of our study are to investigate effects of the offset magnetic field on magnetar kicks and on supernova dynamics. Note that we study a regime where the proto-neutron star formed after collapse has a large magnetic field strength approaching that of a ``magnetar'', a highly magnetized slowly rotating neutron star. As a result, equatorially-asymmetric explosions occur with a formation of the bipolar jets. Resultant magnetar's kick velocities are $\sim 300-1000$ km s$^{-1}$. We find that the acceleration is mainly due to the magnetic pressure while the somewhat weaker magnetic tension works toward the opposite direction, which is due to stronger magnetic field in the northern hemisphere. Noted that observations of magnetar's proper motions are very scarce, our results supply a prediction for future observations. Namely, magnetars possibly have large kick velocities, several hundred km s$^{-1}$, as ordinary neutron stars do, and in an extreme case they could have those up to 1000 km s$^{-1}$.

We perform two-dimensional numerical simulations on the core-collapse of a massive star with strong magnetic fields and differential rotations using a numerical code ZEUS-2D. Changing field configurations and laws of differential rotation parametrically, we compute 14 models and investigate effects of these parameters on the dynamics. In our models, we do not solve the neutrino transport and instead employ a phenomenological parametric EOS that takes into account the neutrino emissions. As a result of the calculations, we find that the field configuration plays a significant role in the dynamics of the core if the initial magnetic field is large enough. Models with initially concentrated fields produce more energetic explosions and more prolate shock waves than the uniform field. Quadrapole-like fields produce remarkably collimated and fast jet, which might be important for gamma-ray bursts(GRB). The Lorentz forces exerted in the region where the plasma-beta is less than unity are responsible for these dynamics. The pure toroidal field, on the other hand, does not lead to any explosion or matter ejection. This suggests the presupernova models of Heger et al.(2003), in which toroidal fields are predominant, is disadvantageous for the magnetorotation-induced supernova considered here. Models with initially weak magnetic fields do not lead to explosion or matter ejection, either. In these models magnetic fields play no role as they do not grow on the timescale considered in this paper so that the magnetic pressure could be comparable to the matter pressure. This is because the exponential field growth as expected in MRI is not seen in our models. The magnetic field is amplified mainly by field-compression and field-wrapping in our simulations.