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.

A cocoon is a by-product of a propagating jet that results from shock heating at the jet head. Herein, considering simultaneous cocoon formation, we study the stability of relativistic jets propagating through the uniform ambient medium. Using a simple analytic argument, we demonstrate that independent from the jet launching condition, the effective inertia of the jet is larger than that of the cocoon when the fully relativistic jet oscillates radially owing to the pressure mismatch between jet and cocoon. In such situations, it is expected that the onset condition for the oscillation-induced Rayleigh-Taylor instability is satisfied at the jet interface, resulting in the destabilization of the relativistic jet during its propagation. We have quantitatively verified and confirmed our prior expectation by performing relativistic hydrodynamic simulations in three dimensions. The possible occurrences of the Richtmyer-Meshkov instability, oscillation-induced centrifugal instability, and Kelvin-Helmholtz instability are also discussed.

We investigate the stability of a neutrino-dominated accretion flow (NDAF), which is expected to be formed in the gravitational collapse of a massive star or the merger of a neutron star binary, based on the variable-$\alpha$ prescription. Recent magnetohydrodynamic (MHD) simulations shows that the viscosity parameter $\alpha$ is proportional to the power of the magnetic Prandtl number ${\rm Pm}=\nu/\eta$, where $\nu$ and $\eta$ are the kinematic viscosity and electric resistivity of the fluid, respectively. In the inner region of a hyperaccretion flow, the viscosity and resistivity are carried by mildly, relativistically degenerated electrons. We fit the dependence of the magnetic Prandtl number on density and temperature by a simple analytic form, and derive the condition for an NDAF to be dynamically unstable. As demonstrations we perform simple one-dimensional simulations of NDAFs with the variable-$\alpha$ and show that the mass accretion becomes highly time-variable in the unstable branch. This mechanism may account for the rapid variability observed in the prompt emission of gamma-ray bursts (GRBs). The mass ejection from a hyperaccretion flow due to viscous heating, which makes a kilonova/macronova emission in the merger of a neutron star binary, is also briefly discussed.

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.

Context: Convectively-driven flows play a crucial role in the dynamo processes that are responsible for producing magnetic activity in stars and planets. It is still not fully understood why many astrophysical magnetic fields have a significant large-scale component. Aims: Our aim is to investigate the dynamo properties of compressible convection in a rapidly rotating Cartesian domain, focusing upon a parameter regime in which the underlying hydrodynamic flow is known to be unstable to a large-scale vortex instability. Methods: The governing equations of three-dimensional nonlinear magnetohydrodynamics (MHD) are solved numerically. Different numerical schemes are compared and we propose a possible benchmark case for other similar codes. Results: In keeping with previous related studies, we find that convection in this parameter regime can drive a large-scale dynamo. The components of the mean horizontal magnetic field oscillate, leading to a continuous overall rotation of the mean field. Whilst the large-scale vortex instability dominates the early evolution of the system, it is suppressed by the magnetic field and makes a negligible contribution to the mean electromotive force that is responsible for driving the large-scale dynamo. The cycle period of the dynamo is comparable to the ohmic decay time, with longer cycles for dynamos in convective systems that are closer to onset. In these particular simulations, large-scale dynamo action is found only when vertical magnetic field boundary conditions are adopted at the upper and lower boundaries. Strongly modulated large-scale dynamos are found at higher Rayleigh numbers, with periods of reduced activity ("grand minima"-like events) occurring during transient phases in which the large-scale vortex temporarily re-establishes itself, before being suppressed again by the magnetic field.

We report the first successful simulation of spontaneous formation of surface magnetic structures from a large-scale dynamo by strongly-stratified thermal convection in Cartesian geometry. The large-scale dynamo observed in our strongly-stratified model has physical properties similar to those in earlier weakly-stratified convective dynamo simulations, indicating that the $\alpha^2$-type mechanism is responsible for it. Additionally to the large-scale dynamo, we find that large-scale structures of the vertical magnetic field are spontaneously formed in the convection zone surface only for the case of strongly-stratified atmosphere. The organization of the vertical magnetic field proceeds in the upper convection zone within tens of convective turn-over time and band-like bipolar structures are recurrently-appeared in the dynamo-saturated stage. We examine possibilities of several candidates as the origin of the surface magnetic structure formation, and then suggest the existence of an as-yet-unknown mechanism for the self-organization of the large-scale magnetic structure, which should be inherent in the strongly-stratified convective atmosphere.

Effects of magnetic field on stellar differential rotation are studied by comparing magnetohydrodynamic (MHD) models and their hydrodynamic (HD) counterparts in the broad range of rotation rate and in varying initial rotation profile. Fully-compressible MHD simulations of rotating penetrative convection are performed in a full-spherical shell geometry. Critical conditions for the transition of the differential rotation between faster equator (solar-type) and slower equator (anti-solar type) are explored with focusing on the "Rossby number (${\rm Ro}$)" and the "convective Rossby number (${\rm Ro}_{\rm conv}$)". It is confirmed that the transition is more gradual and the critical value for it is higher in the MHD model than the HD model in the view of the ${\rm Ro}_{\rm conv}$-dependence. The rotation profile shows, as observed in earlier studies, the bistability near the transition in the HD model, while it disappears when allowing the growth of magnetic fields except for the model with taking anti-solar type solution as the initial condition. We find that the transition occurs at ${\rm Ro} \simeq 1$ both in the MHD and HD models independently of the hysteresis. Not only the critical value, the sharpness of the transition is also similar between the two models in the view of the ${\rm Ro}$-dependence. The influences of the dynamo-generated magnetic field and/or the hysteresis on convective motion are reflected in the ${\rm Ro}$. This would be the reason why the transition is unified in the view of the ${\rm Ro}$-dependence. We finally discuss the ${\rm Ro}$-dependence of magnetic dynamo activities with emphasis on its possible relation to the kinetic helicity profile.

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.

The mechanism of large-scale dynamos in rigidly rotating stratified convection is explored by direct numerical simulations (DNS) in Cartesian geometry. A mean-field dynamo model is also constructed using turbulent velocity profiles consistently extracted from the corresponding DNS results. By quantitative comparison between the DNS and our mean-field model, it is demonstrated that the oscillatory $\alpha^2$ dynamo wave, excited and sustained in the convection zone, is responsible for large-scale magnetic activities such as cyclic polarity reversal and spatiotemporal migration. The results provide strong evidence that a nonuniformity of the $\alpha$-effect, which is a natural outcome of rotating stratified convection, can be an important prerequisite for large-scale stellar dynamos, even without the $\Omega$-effect.

Convective dynamo simulations are performed in local Cartesian geometry. We report the first successful simulation of a large-scale oscillatory dynamo in rigidly rotating convection without stably stratified layers. A key requirement for exciting the large-scale dynamo is a sufficiently long integration time comparable to the ohmic diffusion time. By comparing two models with and without stably stratified layers, their effect on the large-scale dynamo is also studied. The spatiotemporal evolution of the large-scale magnetic field is similar in both models. However, it is intriguing that the magnetic cycle is much shorter in the model without the stable layer than with the stable layer. This suggests that the stable layer impedes the cyclic variations of the large-scale magnetic field.

Spherical solar dynamo simulations are performed. Self-consistent, fully compressible magnetohydrodynamic system with a stably stratified layer below the convective envelope is numerically solved with a newly developed simulation code based on the Yin-Yang grid. The effects of penetrative convection are studied by comparing two models with and without the stable layer. The differential rotation profile in both models is reasonably solar-like with equatorial acceleration. When considering the penetrative convection, a tachocline-like shear layer is developed and maintained beneath the convection zone without assuming any forcing. While turbulent magnetic field becomes predominant in the region where the convective motion is vigorous, mean-field component is preferentially organized in the region where the convective motion is less vigorous. Especially in the stable layer, the strong large-scale field with a dipole symmetry is spontaneously built up. The polarity reversal of the mean-field component takes place globally and synchronously throughout the system regardless the presence of the stable layer. Our results suggest that the stably stratified layer is a key component for organizing the large-scale strong magnetic field, but is not essential for the polarity reversal.

We study the stability of a non-rotating single-component jet using two-dimensional special relativistic hydrodynamic simulations. By assuming translational invariance along the jet axis, we exclude the destabilization effect by Kelvin-Helmholtz mode. The nonlinear evolution of the transverse structure of the jet with a normal jet velocity is highlighted. An intriguing finding in our study is that Rayleigh-Taylor and Richtmeier-Meshkov type instabilities can destroy cylindrical jet configuration as a result of spontaneously induced radial oscillating motion. This is powered by in-situ energy conversion between the thermal and bulk kinetic energies. The effective inertia ratio of the jet to the surrounding medium $\eta$ determines a threshold for the onset of instabilities. The condition $\eta < 1$ should be satisfied for the transverse structure of the jet being persisted.

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.

The effect of rarefaction acceleration on the propagation dynamics and structure of relativistically hot jets is studied through relativistic hydrodynamic simulations. We emphasize the nonlinear interaction of rarefaction waves excited at the interface between a cylindrical jet and the surrounding medium. From simplified one-dimensional models with radial jet structure, we find that a decrease in the relativistic pressure due to the interacting rarefaction waves in the central zone of the jet transiently yields a more powerful boost of the bulk jet than that expected from single rarefaction acceleration. This leads to a cyclic in-situ energy conversion between thermal and bulk kinetic energies which induces radial oscillating motion of the jet. The oscillation timescale is characterized by the initial pressure ratio of the jet to the ambient medium, and follows a simple scaling relation $\tau_{\rm oscillation} \propto (P_{\rm jet,0}/P_{\rm amb,0})^{1/2}$. It is confirmed from extended two-dimensional simulations that this radial oscillating motion in the one-dimensional system manifests as modulation of the structure of the jet in a more realistic situation where a relativistically hot jet propagates through an ambient medium. It is found that when the ambient medium has a power law pressure distribution, the size of the reconfinement region along the propagation direction of the jet in the modulation structure $\lambda$ evolves according to a self-similar relation $\lambda \propto t^{\alpha/2}$ where $\alpha$ is the power-law index of the pressure distribution.

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.

Relativistic Sweet-Parker type magnetic reconnection is investigated by relativistic resistive magnetohydrodynamic (RRMHD) simulations. As an initial setting, we assume anti-parallel magnetic fields and a spatially uniform resistivity. A perturbation imposed on the magnetic fields triggers magnetic reconnection around a current sheet, and the plasma inflows into the reconnection region. The inflows are then heated due to ohmic dissipation in the diffusion region, and finally become relativistically hot outflows. The outflows are not accelerated to ultra-relativistic speeds (i.e., Lorentz factor ~ 1), even when the magnetic energy dominates the thermal and rest mass energies in the inflow region. Most of the magnetic energy in the inflow region is converted into the thermal energy of the outflow during the reconnection process. The energy conversion from magnetic to thermal energy in the diffusion region results in an increase in the plasma inertia. This prevents the outflows from being accelerated to ultra-relativistic speeds. We find that the reconnection rate R obeys the scaling relation R S^{-0.5}, where S is the Lundquist number. This feature is the same as that of non-relativistic reconnection. Our results are consistent with the theoretical predictions of Lyubarsky (2005) for Sweet-Parker type magnetic reconnection.

The nonlinear dynamics of outflows driven by magnetic explosion on the surface of a compact star is investigated through special relativistic magnetohydrodynamic simulations. We adopt, as the initial equilibrium state, a spherical stellar object embedded in hydrostatic plasma which has a density $\rho(r) \propto r^{- \alpha}$ and is threaded by a dipole magnetic field. The injection of magnetic energy at the surface of compact star breaks the equilibrium and triggers a two-component outflow. At the early evolutionary stage, the magnetic pressure increases rapidly around the stellar surface, initiating a magnetically driven outflow. A strong forward shock driven outflow is then excited. The expansion velocity of the magnetically driven outflow is characterized by the Alfv\'en velocity on the stellar surface, and follows a simple scaling relation $v_{\rm mag} \propto {v_{\rm A}}^{1/2}$. When the initial density profile declines steeply with radius, the strong shock is accelerated self-similarly to relativistic velocity ahead of the magnetically driven component. We find that it evolves according to a self-similar relation $\Gamma_{\rm sh} \propto r_{\rm sh}$, where $\Gamma_{\rm sh}$ is the Lorentz factor of the plasma measured at the shock surface $r_{\rm sh}$. Purely hydrodynamic process would be responsible for the acceleration mechanism of the shock driven outflow. Our two-component outflow model, which is the natural outcome of the magnetic explosion, can provide a better understanding of the magnetic active phenomena on various magnetized compact stars.

Based on the characteristics of the magnetorotational instability (MRI) and the MRI-driven turbulence, we construct a steady model for a geometrically thin disk using "non-standard" $\alpha$-prescription. The efficiency of the angular momentum transport depends on the magnetic Prandtl number, $Pm = \nu/\eta$, where $\nu$ and $\eta$ are the microscopic viscous and magnetic diffusivities. In our disk model, Shakura-Sunyaev's $\alpha$-parameter has a power-law dependence on the magnetic Prandtl number, that is $\alpha \propto Pm^\delta$ where $\delta$ is the constant power-law index. Adopting Spitzer's microscopic diffusivities, the magnetic Prandtl number becomes a decreasing function of the disk radius when $\delta > 0$. The transport efficiency of the angular momentum and the viscous heating rate are thus smaller in the outer part of the disk, while these are impacted by the size of index $\delta$. We find that the disk becomes more unstable to the gravitational instability for a larger value of index $\delta$. The most remarkable feature of our disk model is that the thermal and secular instabilities can grow in its middle part even if the radiation pressure is negligibly small in the condition $\delta > 2/3$. In the realistic disk systems, it would be difficult to maintain the steady mass accretion state unless the $Pm$-dependence of MRI-driven turbulence is relatively weak.

The possible role of magneto-rotational instability (MRI) and its driven MHD turbulence in the solar interior is studied on the basis of the linear and nonlinear theories coupling with physical parameters, providing solar rotation profile inverted from the helioseismic observation and a standard model for the internal structure of the sun. We find that the MRI venue is confined to the higher-latitude tachocline and lower-latitude near-surface shear layer. It is especially interesting that the MRI-active region around the tachocline closely overlaps with the area indicating a steep entropy rise which is required from the thermal wind balance in the sun. This suggests that the MRI-driven turbulence plays a crucial role in maintaining the thermal wind balance in the sun via the exceptional turbulent heating and equatorward angular momentum transports. The warm pole existing around the tachocline might be a natural outcome of the turbulent activities energized by the MRI.

We present a theoretical model describing magnetar giant flares on the basis of solar flare/coronal mass ejection theory. In our model, a preflare activity plays a crucial role in driving evaporating flows and supplying baryonic matters into the magnetosphere. The loaded baryonic matter, that is called "prominence", is then gradually uplifted via crustal cracking with maintaining a quasi-force-free equilibrium of the magnetosphere. Finally the prominence is erupted by the magnetic pressure force due to the loss of equilibrium triggered by the explosive magnetic reconnection. The giant flare should be induced as a final outcome of the prominence eruption accompanied by large-scale field reconfigurations. An essential difference between solar and magnetar flares is the control process of their evolutionary dynamics. The flaring activity on magnetars is mainly controlled by the radiative process unlike the solar flare governed by the electron conduction. It is highly suggestive that our model is accountable for the physical properties of the extraordinary giant flare observed on 2004 December 27 from SGR1806-20, including the source of baryonic matters loaded in the expanding ejecta observed after the giant burst.

Axisymmetric magnetorotational instability (MRI) in viscous accretion disks is investigated by linear analysis and two-dimensional nonlinear simulations. The linear growth of the viscous MRI is characterized by the Reynolds number defined as $R_{\rm MRI} \equiv v_A^2/\nu\Omega $, where $v_A$ is the Alfv{\'e}n velocity, $\nu$ is the kinematic viscosity, and $\Omega$ is the angular velocity of the disk. Although the linear growth rate is suppressed considerably as the Reynolds number decreases, the nonlinear behavior is found to be almost independent of $R_{\rm MRI}$. At the nonlinear evolutionary stage, a two-channel flow continues growing and the Maxwell stress increases until the end of calculations even though the Reynolds number is much smaller than unity. A large portion of the injected energy to the system is converted to the magnetic energy. The gain rate of the thermal energy, on the other hand, is found to be much larger than the viscous heating rate. Nonlinear behavior of the MRI in the viscous regime and its difference from that in the highly resistive regime can be explained schematically by using the characteristics of the linear dispersion relation. Applying our results to the case with both the viscosity and resistivity, it is anticipated that the critical value of the Lundquist number $S_{\rm MRI} \equiv v_A^2/\eta\Omega$ for active turbulence depends on the magnetic Prandtl number $S_{{\rm MRI},c} \propto Pm^{1/2}$ in the regime of $Pm \gg 1$ and remains constant when $Pm \ll 1$, where $Pm \equiv S_{\rm MRI}/R_{\rm MRI} = \nu/\eta$ and $\eta$ is the magnetic diffusivity.

The central engine of gamma-ray bursts (GRBs) is believed to be a hot and dense disk with hyperaccretion onto a few solar-mass black hole. We investigate where the magnetorotational instability (MRI) actively operates in the hyperaccretion disk, which can cause angular momentum transport in the disk. The inner region of hyperaccretion disks can be neutrino opaque, and the energy- and momentum-transport by neutrinos could affect the growth of the MRI significantly. Assuming reasonable disk models and a weak magnetic field $B \lesssim 10^{14} \rm{G}$, it is found that the MRI is strongly suppressed by the neutrino viscosity in the inner region of hyperaccretion disks. On the other hand, the MRI can drive active MHD turbulence in the outer neutrino-transparent region regardless of the field strength. This suggests that the baryonic matter is accumulated into the inner dead zone where the MRI grows inactively and the angular momentum transport is inefficient. When the dead zone gains a large amount of mass and becomes gravitationally unstable, intense mass accretion onto the central black hole would occur episodically through the gravitational torque. This process can be a physical mechanism of the short-term variability in the prompt emission of GRBs. Finally, the origin of flaring activities in the X-ray afterglow is predicted in the context of our episodic accretion scenario.

Neutrino radiation takes a major role in the momentum, heat, and lepton transports in proto-neutron stars (PNSs). These diffusive processes affect the growth of magnetorotational instability (MRI) in PNSs. We perform a local linear analysis for the axisymmetric and nonaxisymmetric MRI including the effects of neutrino transports and ohmic dissipation. We find that the MRI can grow even in the multi-diffusive situations that are realized in neutrino loaded PNSs. When the toroidal magnetic component dominates over the poloidal one, nonaxisymmetric MRI modes grow much faster than axisymmetric modes. These results suggest the importance of the nonaxisymmetric MRI in PNSs. Thus the understandings of three-dimensional nonlinear evolutions of the MRI are necessary to reveal the explosion mechanism of core-collapse supernovae.

We investigate the stability of differentially rotating proto-neutron stars (PNSs) with a toroidal magnetic field. Stability criteria for nonaxisymmetric MHD instabilities are derived using a local linear analysis. PNSs are expected to have much stronger radial shear in the rotation velocity compared to normal stars. We find that nonaxisymmetric magnetorotational instability (NMRI) with a large azimuthal wavenumber $m$ is dominant over the kink mode ($m=1$) in differentially rotating PNSs. The growth rate of the NMRI is of the order of the angular velocity $\Omega$ which is faster than that of the kink-type instability by several orders of magnitude. The stability criteria are analogous to those of the axisymmetric magnetorotational instability with a poloidal field, although the effects of leptonic gradients are considered in our analysis. The NMRI can grow even in convectively stable layers if the wavevectors of unstable modes are parallel to the restoring force by the Brunt-V\"ais\"al\"a oscillation. The nonlinear evolution of NMRI could amplify the magnetic fields and drive MHD turbulence in PNSs, which may lead to enhancement of the neutrino luminosity.