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 ⊙ 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.

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 disk model. From our numerical experiments we find 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. Over the entire region of parameter space explored, turbulence driven by the MRI has many characteristic ratios such as that of the Maxwell stress to the magnetic pressure. We also find that the amplitudes of the spatial fluctuations in density and the time variability in the stress are characterized by the ratio of magnetic pressure to gas pressure in the nonlinear regime. 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 nonsteady magnetic reconnection.

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.Comment: to appear in the Astrophysical Journal; 25 pages, 10 figure

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 Ω. If the disk is threaded by a uniform vertical field oriented in the same sense as Ω, the axisymmetric evolution of the MRI is an exponentially growing two-channel flow without saturation. When the field is oppositely directed to Ω, 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.