We have investigated the effects of a smooth transition layer at the contact discontinuity on the growth of the Richtmyer-Meshkov instability (RMI) by hydrodynamic numerical simulations and derived an empirical condition for the suppression of the instability. The transition layer has little influence on the RMI when the thickness $L$ is narrower than the wavelength of an interface modulation $\lambda$. However, if the transition layer becomes broader than $\lambda$, the perturbed velocity associated with the RMI is reduced considerably. The suppression condition is interpreted as the cases that the shock transit time through the transition layer is longer than the sound crossing time of the modulation wavelength. The fluctuation kinetic energy decreases as $L^{-p}$ with $p = 2.5$, which indicates that the growth velocity of the RMI decreases in proportion to $L^{-p/2}$ by the presence of the transition layer. This feature is found to be quite universal and appeared in a wide range of shock-interface interactions.

Thermal fusion plasmas initiated by standing whistler waves are investigated numerically by two- and one-dimensional Particle-in-Cell simulations. When a standing whistler wave collapses due to the wave breaking of ion plasma waves, the energy of the electromagnetic waves transfers directly to the ion kinetic energy. Here we find that the ion heating by the standing whistler wave is operational even in multi-dimensional simulations of multi-ion species targets, such as deuterium-tritium (DT) ices and solid ammonia borane (H$_6$BN). The energy conversion efficiency to ions becomes as high as 15% of the injected laser energy, which depends significantly on the target thickness and laser pulse duration. The ion temperature could reach a few tens of keV or much higher if appropriate laser-plasma conditions are selected. DT fusion plasmas generated by this method must be useful as efficient neutron sources. Our numerical simulations suggest that the neutron generation efficiency exceeds 10$^9$ n/J per steradian, which is beyond the current achievements of the state-of-the-art laser experiments. The standing whistler wave heating would expand the experimental possibility for an alternative ignition design of magnetically confined laser fusion, and also for more difficult fusion reactions including the aneutronic proton-boron reaction.

Efficient energy transfer from electromagnetic waves to ions has been demanded to control laboratory plasmas for various applications and could be useful to understand the nature of space and astrophysical plasmas. However, there exists a severe unsolved problem that most of the wave energy is converted quickly to electrons, but not to ions. Here, an energy conversion process to ions in overdense plasmas associated with whistler waves is investigated by numerical simulations and theoretical model. Whistler waves propagating along a magnetic field in space and laboratories often form the standing waves by the collision of counter-propagating waves or through the reflection. We find that ions in the standing whistler waves acquire a large amount of energy directly from the waves in a short timescale comparable to the wave oscillation period. Thermalized ion temperature increases in proportion to the square of the wave amplitude and becomes much higher than the electron temperature in a wide range of wave-plasma conditions. This efficient ion-heating mechanism applies to various plasma phenomena in space physics and fusion energy sciences.

We developed an experimental platform for studying magnetic reconnection in an external magnetic field with simultaneous measurements of plasma imaging, flow velocity, and magnetic-field variation. Here, we investigate the stagnation and acceleration in counter-streaming plasmas generated by high-power laser beams. A plasma flow perpendicular to the initial flow directions is measured with laser Thomson scattering. The flow is, interestingly, accelerated toward the high-density region, which is opposite to the direction of the acceleration by pressure gradients. This acceleration is possibly interpreted by the interaction of two magnetic field loops initially generated by Biermann battery effect, resulting in a magnetic reconnection forming a single field loop and additional acceleration by a magnetic tension force.

A preliminary numerical experiment is conducted for laboratory experiments on the generation of magnetized collisionless shocks with high-power lasers by using one-dimensional particle-in-cell simulation. The present study deals with the interaction between a moving Aluminum plasma and a Nitrogen plasma at rest. In the numerical experiment, the Nitrogen plasma is unmagnetized or magnetized by a weak external magnetic field. Since the previous study suggested the generation of spontaneous magnetic field in the piston (Aluminum) plasma due to the Biermann battery, the effect of the magnetic field is of interest. Sharp jumps of electron density and magnetic field are observed around the interface between the two plasmas as long as one of the two plasmas is magnetized, which indicates the formation of tangential electron-magneto-hydro-dynamic discontinuity. When the Aluminum plasma is magnetized, strong compression of both density and magnetic field takes place in the pure Aluminum plasma during the gyration of Nitrogen ions in the Aluminum plasma region. The formation of a shock downstream is indicated from the shock jump condition. The result suggests that the spontaneous magnetic field in the piston (Aluminum) plasma plays an essential role in the formation of a perpendicular collisionless shock.

Water, ethanol, and ammonia are the key components of the mantles of Uranus and Neptune. To improve structure and evolution models and give an explanation of the magnetic fields and luminosities of the icy giants, those components need to be characterised at planetary conditions (some Mbar and a few $10^3$ K). Those conditions are typical of the Warm Dense Matter regime, which exhibits a rich phase diagram, with the coexistence of many states of matter and a large variety of chemical processes. H$_2$O, C:H:O, and C:H:N:O mixtures have been compressed up to 2.8 Mbar along the principal Hugoniot using laser-driven decaying shocks. The experiments were performed at the GEKKO XII and LULI 2000 laser facilities using standard optical diagnostics (Doppler velocimetry and pyrometry) to characterise equation of state and optical reflectivity of the shocked states. The results show that H$_2$O and the C:H:N:O mixture share the same equation of state with a density scaling, while the reflectivity behaves differently by what concerns both the onset pressures and the saturation values. The reflectivity measurement at two frequencies allows to estimate the conductivity and the complex refractive index using a Drude model.

A new target design is presented to model high-energy radiative accretion shocks in polars. In this paper, we present the experimental results obtained on the GEKKO XII laser facility for the POLAR project. The experimental results are compared with 2D FCI2 simulations to characterize the dynamics and the structure of plasma flow before and after the collision. The good agreement between simulations and experimental data confirm the formation of a reverse shock where cooling losses start modifying the post-shock region. With the multi-material structure of the target, a hydrodynamic collimation is exhibited and a radiative structure coupled with the reverse shock is highlighted in both experimental data and simulations. The flexibility on the laser energy produced on GEKKO XII, allowed us to produce high-velocity flows and study new and interesting radiation hydrodynamic regimes between those obtained on the LULI2000 and Orion laser facilities.

The interaction of dense plasmas with an intense laser under a strong external magnetic field has been investigated. When the cyclotron frequency for the ambient magnetic field is higher than the laser frequency, the laser's electromagnetic field is converted to the whistler mode that propagates along the field line. Because of the nature of the whistler wave, the laser light penetrates into dense plasmas with no cutoff density, and produces superthermal electrons through cyclotron resonance. It is found that the cyclotron resonance absorption occurs effectively under the broadened conditions, or a wider range of the external field, which is caused by the presence of relativistic electrons accelerated by the laser field. The upper limit of the ambient field for the resonance increases in proportion to the square root of the relativistic laser intensity. The propagation of a large-amplitude whistler wave could raise the possibility for plasma heating and particle acceleration deep inside dense plasmas.

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.

We report an experimental study of the phase diagrams of periclase (MgO), enstatite (MgSiO3) and forsterite (Mg2SiO4) at high pressures. We investigated with laser driven decaying shocks the pressure/temperature curves of MgO, MgSiO3 and Mg2SiO4 between 0.2-1.2 TPa, 0.12-0.5 TPa and 0.2-0.85 TPa respectively. A melting signature has been observed in MgO at 0.47 TPa and 9860 K, while no phase changes were observed neither in MgSiO3 nor in Mg2SiO4. An increasing of reflectivity of MgO, MgSiO3 and Mg2SiO4 liquids have been detected at 0.55 TPa -12 760 K, 0.15 TPa - 7540 K, 0.2 TPa - 5800 K, respectively. In contrast to SiO2, melting and metallization of these compounds do not coincide implying the presence of poor electrically conducting liquids close to the melting lines. This has important implications for the generation of dynamos in Super-earths mantles.

Weakly Interacting Massive Particles (WIMPs) may constitute a large fraction of the matter in the Universe. There are excess events in the data of DAMA/LIBRA, CoGeNT, CRESST-II, and recently CDMS-Si, which could be consistent with WIMP masses of approximately 10 GeV/c2. However, for MDM > 10 GeV/c2 null results of the CDMS-Ge, XENON, and LUX detectors may be in tension with the potential detections for certain dark matter scenarios and assuming a certain light response. We propose the use of a new class of biological dark matter (DM) detectors to further examine this light dark matter hypothesis, taking advantage of new signatures with low atomic number targets, Two types of biological DM detectors are discussed here: DNA-based detectors and enzymatic reactions (ER) based detectors. In the case of DNA-based detectors, we discuss a new implementation. In the case of ER detectors, there are four crucial phases of the detection process: a) change of state due to energy deposited by a particle; b) amplification due to the release of energy derived from the action of an enzyme on its substrate; c) sustainable but non-explosive enzymatic reaction; d) self-termination due to the denaturation of the enzyme, when the temperature is raised. This paper provides information of how to design as well as optimize these four processes.

A large set of numerical simulations of magnetohydrodynamic (MHD) turbulence induced by the magnetorotational instability (MRI) is presented. Revisiting the previous survey conducted by Sano et al. (2004), we investigate the gas pressure dependence of the saturation level. In ideal MHD simulations, the gas pressure dependence is found to be very sensitive to the choice of a numerical scheme. This is because the numerical magnetic Prandtl number varies according to the scheme as well as the pressure, which considerably affects the results. The saturation level is more sensitive to the numerical magnetic Prandtl number than the pressure. In MHD simulations with explicit viscosity and resistivity, the saturation level increases with the physical magnetic Prandtl number, and it is almost independent of the gas pressure when the magnetic Prandtl number is constant. This is indicative of the incompressible turbulence saturated by the secondary tearing instability.

Weakly Interacting Massive Particles (WIMPs) may constitute most of the matter in the Universe. The ability to detect the directionality of recoil nuclei will considerably facilitate detection of WIMPs. In this paper we propose a novel type of dark matter detector: detectors made of DNA or RNA could provide nanometer resolution for tracking, an energy threshold of 0.5 keV, and can operate at room temperature. When a WIMP from the Galactic Halo elastically scatters off of a nucleus in the detector, the recoiling nucleus then traverses hundreds of strings of single stranded nucleic acids (ssNA) with known base sequences and severs ssNA strands along its trajectory. The location of the break can be identified by amplifying and identifying the segments of cut ssNA using techniques well known to biologists. Thus the path of the recoiling nucleus can be tracked to nanometer accuracy. In one such detector concept, the transducers are nanometer-thick Au-foils of 1m x 1m, and the direction of recoiling nuclei is measured by "NA Tracking Chamber" consisting of ordered array of ssNA strands. Polymerase Chain Reaction (PCR) and ssNA sequencing are used to read-out the detector. The proposed detector is smaller and cheaper than other alternatives: 1 kg of gold and 0.1 to 4 kg of ssNA (depending on length and strand density), packed into 0.01m$^3$, can be used to study 10 GeV WIMPs. A variety of other detector target elements could be used in this detector to optimize for different WIMP masses and to identify WIMP properties. By leveraging advances in molecular biology, we aim to achieve about 1,000-fold better spatial resolution than in conventional WIMP detectors at reasonable cost.

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.

We present the results of local, vertically stratified, radiation MHD shearing box simulations of MRI turbulence appropriate for the hydrogen ionizing regime of dwarf nova and soft X-ray transient outbursts. We incorporate the frequency-integrated opacities and equation of state for this regime, but neglect non-ideal MHD effects and surface irradiation, and do not impose net vertical magnetic flux. We find two stable thermal equilibrium tracks in the effective temperature versus surface mass density plane, in qualitative agreement with the S-curve picture of the standard disk instability model. We find that the large opacity at temperatures near $10^4$K, a corollary of the hydrogen ionization transition, triggers strong, intermittent thermal convection on the upper stable branch. This convection strengthens the magnetic turbulent dynamo and greatly enhances the time-averaged value of the stress to thermal pressure ratio $\alpha$, possibly by generating vertical magnetic field that may seed the axisymmetric magnetorotational instability, and by increasing cooling so that the pressure does not rise in proportion to the turbulent dissipation. These enhanced stress to pressure ratios may alleviate the order of magnitude discrepancy between the $\alpha$-values observationally inferred in the outburst state and those that have been measured from previous local numerical simulations of magnetorotational turbulence that lack net vertical magnetic flux.

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.

We examine the conditions under which the disks of gas and dust orbiting young gas giant planets are sufficiently conducting to experience turbulence driven by the magneto-rotational instability. By modeling the ionization and conductivity in the disk around proto-Jupiter, we find that turbulence is possible if the X-rays emitted near the Sun reach the planet's vicinity and either (1) the gas surface densities are in the range of the minimum-mass models constructed by augmenting Jupiter's satellites to Solar composition, while dust is depleted from the disk atmosphere, or (2) the surface densities are much less, and in the range of gas-starved models fed with material from the Solar nebula, but not so low that ambipolar diffusion decouples the neutral gas from the plasma. The results lend support to both minimum-mass and gas-starved models of the protojovian disk: (1) The dusty minimum-mass models have internal conductivities low enough to prevent angular momentum transfer by magnetic forces, as required for the material to remain in place while the satellites form. (2) The gas-starved models have magnetically-active surface layers and a decoupled interior "dead zone". Similar active layers in the Solar nebula yield accretion stresses in the range assumed in constructing the circumjovian gas-starved models. Our results also point to aspects of both classes of models that can be further developed. Non-turbulent minimum-mass models will lose dust from their atmospheres by settling, enabling gas to accrete through a thin surface layer. For the gas-starved models it is crucial to learn whether enough stellar X-ray and ultraviolet photons reach the circumjovian disk. Additionally the stress-to-pressure ratio ought to increase with distance from the planet, likely leading to episodic accretion outbursts.

The critical strength of a magnetic field required for the suppression of the Richtmyer-Meshkov instability (RMI) is investigated numerically by using a two-dimensional single-mode analysis. For the cases of MHD parallel shocks, the RMI can be stabilized as a result of the extraction of vorticity from the interface. A useful formula describing a critical condition for MHD RMI has been introduced, and which is successfully confirmed by the direct numerical simulations. The critical field strength is found to be largely depending on the Mach number of the incident shock. If the shock is strong enough, even low-$\beta$ plasmas can be subject to the growth of the RMI.

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 amplification of a magnetic field due to the Richtmyer-Meshkov instability (RMI) is investigated by two-dimensional MHD simulations. Single-mode analysis is adopted to reveal definite relation between the nonlinear evolution of RMI and the field enhancement. It is found that an ambient magnetic field is stretched by fluid motions associated with the RMI, and the strength is amplified significantly by more than two orders of magnitude. The saturation level of the field is determined by a balance between the amplified magnetic pressure and the thermal pressure after shock passage. This effective amplification can be achieved in a wide range of the conditions for the RMI such as the Mach number of an incident shock and the density ratio at a contact discontinuity. The results suggest that the RMI could be a robust mechanism of the amplification of interstellar magnetic fields and cause the origin of localized strong fields observed at the shock of supernova remnants.

The principal Hugoniot for liquid hydrogen was obtained up to 55 GPa under laser-driven shock loading. Pressure and density of compressed hydrogen were determined by impedance-matching to a quartz standard. The shock temperature was independently measured from the brightness of the shock front. Hugoniot data of hydrogen provide a good benchmark to modern theories of condensed matter. The initial number density of liquid hydrogen is lower than that for liquid deuterium, and this results in shock compressed hydrogen having a higher compression and higher temperature than deuterium at the same shock pressure.

We apply ionization balance and MHD calculations to investigate whether magnetic activity moderated by recombination on dust can account for the mass accretion rates and the mid-infrared spectra and variability of protostellar disks. The MHD calculations use the stratified shearing-box approach and include grain settling and the feedback from the changing dust abundance on the resistivity of the gas. The two-decade spread in accretion rates among T Tauri stars is too large to result solely from variety in the grain size and stellar X-ray luminosity, but can be produced by varying these together with the disk magnetic flux. The diversity in the silicate bands can come from the coupling of grain settling to the distribution of the magneto-rotational turbulence, through three effects: (1) Recombination on grains yields a magnetically inactive dead zone extending above two scale heights, while turbulence in the magnetically active disk atmosphere overshoots the dead zone boundary by only about one scale height. (2) Grains deep in the dead zone oscillate vertically in waves driven by the turbulent layer above, but on average settle at the laminar rates, so the interior of the dead zone is a particle sink and the disk atmosphere becomes dust-depleted. (3) With sufficient depletion, the dead zone is thinner and mixing dredges grains off the midplane. The MHD results also show that the magnetic activity intermittently lifts clouds of dust into the atmosphere. The photosphere height changes by up to one-third over a few orbits, while the extinction along lines of sight grazing the disk surface varies by factors of two over times down to 0.1 orbit. We suggest that the changing shadows cast by the dust clouds on the outer disk are a cause of the daily to monthly mid-infrared variability in some young stars. (Abridged.)

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.

Planets form inside protostellar disks in a dead zone where the electrical resistivity of the gas is too high for magnetic forces to drive turbulence. We show that much of the dead zone nevertheless is active and flows toward the star while smooth, large-scale magnetic fields transfer the orbital angular momentum radially outward. Stellar X-ray and radionuclide ionization sustain a weak coupling of the dead zone gas to the magnetic fields, despite the rapid recombination of free charges on dust grains. Net radial magnetic fields are generated in the magneto-rotational turbulence in the electrically conducting top and bottom surface layers of the disk, and reach the midplane by Ohmic diffusion. A toroidal component to the fields is produced near the midplane by the orbital shear. The process is similar to the magnetization of the Solar tachocline. The result is a laminar, magnetically-driven accretion flow in the region where the planets form.

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.

We investigate the conditions for the presence of a magnetically inactive dead zone in protostellar disks, using 3-D shearing-box MHD calculations including vertical stratification, Ohmic resistivity and time-dependent ionization chemistry. Activity driven by the magnetorotational instability fills the whole thickness of the disk at 5 AU, provided cosmic ray ionization is present, small grains are absent and the gas-phase metal abundance is sufficiently high. At 1 AU the larger column density of 1700 g/cm^2 means the midplane is shielded from ionizing particles and remains magnetorotationally stable even under the most favorable conditions considered. Nevertheless the dead zone is effectively eliminated. Turbulence mixes free charges into the interior as they recombine, leading to a slight coupling of the midplane gas to the magnetic fields. Weak, large-scale radial fields diffuse to the midplane where they are sheared out to produce stronger azimuthal fields. The resulting midplane accretion stresses are just a few times less than in the surface layers on average.

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.

We analyze the ionization state of the magnetohydrodynamically turbulent protoplanetary disks and propose a new mechanism of sustaining ionization. First, we show that in the quasi-steady state of turbulence driven by magnetorotational instability in a typical protoplanetary disk with dust grains, the amount of energy dissipation should be sufficient for providing the ionization energy that is required for activating magnetorotational instability. Second, we show that in the disk with dust grains the energetic electrons that compose electric currents in weakly ionized gas can provide collisional ionization, depending on the actual saturation state of magnetorotational turbulence. On the other hand, we show that in the protoplanetary disks with the reduced effect of dust grains, the turbulent motion can homogenize the ionization degree, leading to the activation of magnetorotational instability even in the absence of other ionization processes. The results in this Letter indicate that most of the regions in protoplanetary disks remain magnetically active, and we thus require a change in the theoretical modeling of planet formation.

The saturation level of the magnetorotational instability (MRI) is investigated using three-dimensional MHD simulations. The shearing box approximation is adopted and the vertical component of gravity is ignored, so that the evolution of the MRI is followed in a small local part of the disk. We focus on the dependence of the saturation level of the stress on the gas pressure, which is a key assumption in the standard alpha disk model. From our numerical experiments it is found that there is a weak power-law relation between the saturation level of the Maxwell stress and the gas pressure in the nonlinear regime; the higher the gas pressure, the larger the stress. Although the power-law index depends slightly on the initial field geometry, the relationship between stress and gas pressure is independent of the initial field strength, and is unaffected by Ohmic dissipation if the magnetic Reynolds number is at least 10. The relationship is the same in adiabatic calculations, where pressure increases over time, and nearly-isothermal calculations, where pressure varies little with time. Our numerical results are qualitatively consistent with an idea that the saturation level of the MRI is determined by a balance between the growth of the MRI and the dissipation of the field through reconnection. The quantitative interpretation of the pressure-stress relation, however, may require advances in the theoretical understanding of non-steady magnetic reconnection.

We examine the small-scale dynamics of black hole accretion disks in which radiation pressure exceeds gas pressure. Local patches of disk are modeled by numerically integrating the equations of radiation MHD in the flux-limited diffusion approximation. The shearing-box approximation is used, and the vertical component of gravity is neglected. Magneto-rotational instability (MRI) leads to turbulence in which accretion stresses are due primarily to magnetic torques. When radiation is locked to gas over the length and time scales of fluctuations in the turbulence, the accretion stress, density contrast, and dissipation differ little from those in the corresponding calculations with radiation replaced by extra gas pressure. However, when radiation diffuses each orbit a distance that is comparable to the RMS vertical wavelength of the MRI, radiation pressure is less effective in resisting squeezing. Large density fluctuations occur, and radiation damping of compressive motions converts PdV work into photon energy. The accretion stress in calculations having a net vertical magnetic field is found to be independent of opacity over the range explored, and approximately proportional to the square of the net field. In calculations with zero net magnetic flux, the accretion stress depends on the portion of the total pressure that is effective in resisting compression. The stress is lower when radiation diffuses rapidly with respect to the gas. We show that radiation-supported Shakura-Sunyaev disks accreting via internal magnetic stresses are likely in their interiors to have radiation marginally coupled to turbulent gas motions.

The evolution of the magnetorotational instability (MRI) during the transition from outburst to quiescence in a dwarf nova disk is investigated using three-dimensional MHD simulations. The shearing box approximation is adopted for the analysis, so that the efficiency of angular momentum transport is studied in a small local patch of the disk: this is usually referred as to a one-zone model. To take account of the low ionization fraction of the disk, the induction equation includes both ohmic dissipation and the Hall effect. We induce a transition from outburst to quiescence by an instantaneous decrease of the temperature. The evolution of the MRI during the transition is found to be very sensitive to the temperature of the quiescent disk. As long as the temperature is higher than a critical value of about 2000 K, MHD turbulence and angular momentum transport is sustained by the MRI. However, MHD turbulence dies away within an orbital time if the temperature falls below this critical value. In this case, the stress drops off by more than 2 orders of magnitude, and is dominated by the Reynolds stress associated with the remnant motions from the outburst. The critical temperature depends slightly on the distance from the central star and the local density of the disk.

The nonlinear evolution of the magnetorotational instability (MRI) in weakly ionized accretion disks, including the effect of the Hall term and ohmic dissipation, is investigated using local three-dimensional MHD simulations and various initial magnetic field geometries. When the magnetic Reynolds number, Re_M \equiv v_A^2 / \eta \Omega (where v_A is the Alfven speed, \eta the magnetic diffusivity, and \Omega the angular frequency), is initially larger than a critical value Re_{M, crit}, the MRI evolves into MHD turbulence in which angular momentum is transported efficiently by the Maxwell stress. If Re_M < Re_{M, crit}, however, ohmic dissipation suppresses the MRI, and the stress is reduced by several orders of magnitude. The critical value is in the range of 1 - 30 depending on the initial field configuration. The Hall effect does not modify the critical magnetic Reynolds number by much, but enhances the saturation level of the Maxwell stress by a factor of a few. We show that the saturation level of the MRI is characterized by v_{Az}^2 / \eta \Omega, where v_{Az} is the Alfven speed in the nonlinear regime along the vertical component of the field. The condition for turbulence and significant transport is given by v_{Az}^2 / \eta \Omega \gtrsim 1, and this critical value is independent of the strength and geometry of the magnetic field or the size of the Hall term. If the magnetic field strength in an accretion disk can be estimated observationally, and the magnetic Reynolds number v_A^2 / \eta \Omega is larger than about 30, this would imply the MRI is operating in the disk.

The effect of the Hall term on the evolution of the magnetorotational instability (MRI) in weakly ionized accretion disks is investigated using local axisymmetric simulations. First, we show that the Hall term has important effects on the MRI when the temperature and density in the disk is below a few thousand K and between 10^13 and 10^18 cm^{-3} respectively. Such conditions can occur in the quiescent phase of dwarf nova disks, or in the inner part (inside 10 - 100 AU) of protoplanetary disks. When the Hall term is important, the properties of the MRI are dependent on the direction of the magnetic field with respect to the angular velocity vector \Omega. If the disk is threaded by a uniform vertical field oriented in the same sense as \Omega, the axisymmetric evolution of the MRI is an exponentially growing two-channel flow without saturation. When the field is oppositely directed to \Omega, however, small scale fluctuations prevent the nonlinear growth of the channel flow and the MRI evolves into MHD turbulence. These results are anticipated from the characteristics of the linear dispersion relation. In axisymmetry on a field with zero-net flux, the evolution of the MRI is independent of the size of the Hall term relative to the inductive term. The evolution in this case is determined mostly by the effect of ohmic dissipation.

We perform numerical simulations of magneto-rotational instability in a local patch of accretion disk in which radiation pressure exceeds gas pressure. Such conditions may occur in the central regions of disks surrounding compact objects in active galactic nuclei and Galactic X-ray sources. We assume axisymmetry, and neglect vertical stratification. The growth rates of the instability on initially uniform magnetic fields are consistent with the linear analysis of Blaes & Socrates 2001. As is the case when radiation effects are neglected, the non-linear development of the instability leads to transitory turbulence when the initial magnetic field has no net vertical flux. During the turbulent phase, angular momentum is transported outwards. The Maxwell stress is a few times the Reynolds stress, and their sum is about four times the mean pressure in the vertical component of the magnetic field. For magnetic pressure exceeding gas pressure, turbulent fluctuations in the field produce density contrasts about equal to the ratio of magnetic to gas pressure. These are many times larger than in the corresponding gas pressure dominated situation, and may have profound implications for the steady-state vertical structure of radiation-dominated disks. Diffusion of radiation from compressed regions damps turbulent motions, converting kinetic energy into photon energy.

The nonlinear evolution and the saturation mechanism of the magnetorotational instability (MRI) are investigated using three-dimensional resistive MHD simulations. A local shearing box is used for our numerical analysis and the simulations begin with a purely vertical magnetic field. We find that the magnetic stress in the nonlinear stage of the MRI is strongly fluctuating. The time evolution shows the quasi-periodic recurrence of spike-shape variations typically for a few orbits which correspond to the rapid amplification of the magnetic field by the nonlinear growth of a two-channel solution followed by the decay through magnetic reconnections. The increase rate of the total energy in the shearing box system is analytically related to the volume-averaged torque in the system. We find that at the saturated state this energy gain of the system is balanced with the increase of the thermal energy mostly due to the joule heating. The spike-shape time evolution is a general feature of the nonlinear evolution of the MRI in the disks threaded by vertical fields and can be seen if the effective magnetic Reynolds number is larger than about unity.

We investigate where in protoplanetary disks magnetorotational instability operates, which can cause angular momentum transport in the disks. We investigate the spatial distribution of various charged particles and the unstable regions for a variety of models for protoplanetary disks taking into account the recombination of ions and electrons at grain surfaces, which is an important process in most parts of the disks. We find that for all the models there is an inner region which is magnetorotationally stable due to ohmic dissipation. This must make the accretion onto the central star non-steady. For the model of the minimum-mass solar nebula, the critical radius, inside of which the disk is stable, is about 20 AU, and the mass accretion rate just outside the critical radius is 10^{-7} - 10^{-6} M_{\odot} yr^{-1}. The stable region is smaller in a disk of lower column density. Dust grains in protoplanetary disks may grow by mutual sticking and may sediment toward the midplane of the disks. We find that the stable region shrinks as the grain size increases or the sedimentation proceeds. Therefore in the late evolutionary stages, protoplanetary disks can be magnetorotationally unstable even in the inner regions.