No Arabic abstract
By performing local three-dimensional MHD simulations of stratified accretion disks, we investigate disk winds driven by MHD turbulence. Initially given weak vertical magnetic fields are effectively amplified by magnetorotational instability and winding due to differential rotation. Large scale channel flows develop most effectively at 1.5 - 2 times the scale heights where the magnetic pressure is comparable to but slightly smaller than the gas pressure. The breakup of these channel flows drives structured disk winds by transporting the Poynting flux to the gas. These features are universally observed in the simulations of various initial fields. This disk wind process should play an essential role in the dynamical evaporation of proto-planetary disks. The breakup of channel flows also excites the momentum fluxes associated with Alfvenic and (magneto-)sonic waves toward the mid-plane, which possibly contribute to the sedimentation of small dust grains in protoplanetary disks.
We investigate the roles of magnetically driven disk wind (MDW) and thermally driven photoevaporative wind (PEW) in the long-time evolution of protoplanetary disks. We start simulations from the early phase in which the disk mass is $0.118,{mathrm{M}_{odot}}$ around a $1,{mathrm{M}_{odot}}$ star and track the evolution until the disk is completely dispersed. We incorporate the mass loss by PEW and the mass loss and magnetic braking (wind torque) by MDW, in addition to the viscous accretion, viscous heating, and stellar irradiation. We find that MDW and PEW respectively have different roles: magnetically driven wind ejects materials from an inner disk in the early phase, whereas photoevaporation has a dominant role in the late phase in the outer ($gtrsim1,$au) disk. The disk lifetime, which depends on the combination of MDW, PEW, and viscous accretion, shows a large variation of $sim1$-$20,$Myr; the gas is dispersed mainly by the MDW and the PEW in the cases with a low viscosity and the lifetime is sensitive to the mass-loss rate and torque of the MDW, whereas the lifetime is insensitive to these parameters when the viscosity is high. Even in disks with very weak turbulence, the cooperation of MDW and PEW enables the disk dispersal within a few Myr.
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/ uOmega $, where $v_A$ is the Alfv{e}n velocity, $ u$ 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/etaOmega$ 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} = u/eta$ and $eta$ is the magnetic diffusivity.
The intermittent small-scale structure of turbulence governs energy dissipation in many astrophysical plasmas and is often believed to have universal properties for sufficiently large systems. In this work, we argue that small-scale turbulence in accretion disks is universal in the sense that it is insensitive to the magnetorotational instability (MRI) and background shear, and therefore indistinguishable from standard homogeneous magnetohydrodynamic (MHD) turbulence at small scales. We investigate the intermittency of current density, vorticity, and energy dissipation in numerical simulations of incompressible MHD turbulence driven by the MRI in a shearing box. We find that the simulations exhibit a similar degree of intermittency as in standard MHD turbulence. We perform a statistical analysis of intermittent dissipative structures and find that energy dissipation is concentrated in thin sheet-like structures that span a wide range of scales up to the box size. We show that these structures exhibit strikingly similar statistical properties to those in standard MHD turbulence. Additionally, the structures are oriented in the toroidal direction with a characteristic tilt of approximately 17.5 degrees, implying an effective guide field in that direction.
We perform 3D vertically-stratified local shearing-box ideal MHD simulations of the magnetorotational instability (MRI) that include a net vertical magnetic flux, which is characterized by beta_0 (ratio of gas pressure to magnetic pressure of the net vertical field at midplane). We have considered beta_0=10^2, 10^3 and 10^4 and in the first two cases the most unstable linear MRI modes are well resolved in the simulations. We find that the behavior of the MRI turbulence strongly depends on beta_0: The radial transport of angular momentum increases with net vertical flux, achieving alpha=0.08 for beta_0=10^4 and alpha>1.0 for beta_0=100, where alpha is the Shakura-Sunyaev parameter. A critical value lies at beta_0=10^3: For beta_0>10^3, the disk consists of a gas pressure dominated midplane and a magnetically dominated corona. The turbulent strength increases with net flux, and angular momentum transport is dominated by turbulent fluctuations. The magnetic dynamo that leads to cyclic flips of large-scale fields still exists, but becomes more sporadic as net flux increases. For beta_0<10^3, the entire disk becomes magnetic dominated. The turbulent strength saturates, and the magnetic dynamo is quenched. Stronger large-scale fields are generated with increasing net flux, which dominates angular momentum transport. A strong outflow is launched from the disk by the magnetocentrifugal mechanism, and the mass flux increases linearly with net vertical flux and shows sign of saturation at beta_0=10^2. However, the outflow is unlikely to be directly connected to a global wind: for beta_0>10^3, the large-scale field has no permanent bending direction due to dynamo activities, while for beta_0<10^3, the outflows from the top and bottom sides of the disk bend towards opposite directions, inconsistent with a physical disk wind geometry. Global simulations are needed to address the fate of the outflow.
We review the properties of dust in protoplanetary disks around optically visible pre-main sequence stars obtained with a variety of observational techniques, from measurements of scattered light at visual and infrared wavelengths to mid-infrared spectroscopy and millimeter interferometry. A general result is that grains in disks are on average much larger than in the diffuse interstellar medium (ISM). In many disks, there is evidence that a large mass of dust is in grains with millimeter and centimeter sizes, more similar to sand and pebbles than to grains. Smaller grains (with micron-sizes) exist closer to the disk surface, which also contains much smaller particles, e.g., polycyclic aromatic hydrocarbons. There is some evidence of a vertical stratification, with smaller grains closer to the surface. Another difference with ISM is the higher fraction of crystalline relative to amorphous silicates found in disk surfaces. There is a large scatter in dust properties among different sources, but no evidence of correlation with the stellar properties, for samples that include objects from intermediate to solar mass stars and brown dwarfs. There is also no apparent correlation with the age of the central object, over a range roughly between 1 and 10 Myr. This suggests a scenario where significant grain processing may occur very early in the disk evolution, possibly when it is accreting matter from the parental molecular core. Further evolution may occur, but not necessarily rapidly, since we have evidence that large amounts of grains, from micron to centimeter size, can survive for periods as long as 10 Myr.