No Arabic abstract
Global evolution and dispersal of protoplanetary disks (PPDs) is governed by disk angular momentum transport and mass-loss processes. Recent numerical studies suggest that angular momentum transport in the inner region of PPDs is largely driven by magnetized disk wind, yet the wind mass-loss rate remains unconstrained. On the other hand, disk mass loss has conventionally been attributed to photoevaporation, where external heating on the disk surface drives a thermal wind. We unify the two scenarios by developing a 1D model of magnetized disk winds with a simple treatment of thermodynamics as a proxy for external heating. The wind properties largely depend on 1) the magnetic field strength at the wind base, characterized by the poloidal Alfven speed $v_{Ap}$, 2) the sound speed $c_s$ near the wind base, and 3) how rapidly poloidal field lines diverge (achieve $R^{-2}$ scaling). When $v_{Ap}gg c_s$, corotation is enforced near the wind base, resulting in centrifugal acceleration. Otherwise, the wind is accelerated mainly by the pressure of the toroidal magnetic field. In both cases, the dominant role played by magnetic forces likely yields wind outflow rates that well exceed purely hydrodynamical mechanisms. For typical PPD accretion-rate and wind-launching conditions, we expect $v_{Ap}$ to be comparable to $c_s$ at the wind base. The resulting wind is heavily loaded, with total wind mass loss rate likely reaching a considerable fraction of wind-driven accretion rate. Implications for modeling global disk evolution and planet formation are also discussed.
We present a model for the dispersal of protoplanetary disks by winds from either the central star or the inner disk. These winds obliquely strike the flaring disk surface and strip away disk material by entraining it in an outward radial-moving flow at the wind-disk interface which lies several disk scale heights above the mid-plane. The disk dispersal time depends on the entrainment velocity at which disk material flows into this turbulent shear layer interface. If the entrainment efficiency is ~10% of the local sound speed, a likely upper limit, the dispersal time at 1 AU is ~6 Myr for a disk with a surface density of 10^3 g cm^{-2}, a solar mass central star, and a wind with an outflow rate 10^{-8} Msun/yr and terminal velocity 200 km/s. When compared to photoevaporation and viscous evolution, wind stripping can be a dominant mechanism only for the combination of low accretion rates (< 10^{-8} Msun/yr) and wind outflow rates approaching these accretion rates. This case is unusual since generally outflow rates are < 0.1 of of accretion rates.
The gas dynamics of protoplanetary disks (PPDs) is largely controlled by non-ideal magnetohydrodynamic (MHD) effects including Ohmic resistivity, the Hall effect and ambipolar diffusion. Among these the role of the Hall effect is the least explored and most poorly understood. We have included all three non-ideal MHD effects in a self-consistent manner to investigate the role of the Hall effect on PPD gas dynamics using local shearing-box simulations. In this first paper, we focus on the inner region of PPDs, where previous studies excluding the Hall effect have revealed that the inner disk up to ~10 AU is largely laminar, with accretion driven by a magnetocentrifugal wind. We confirm this basic picture and show that the Hall effect introduces modest modifications to the wind solutions, depending on the polarity of the large-scale poloidal magnetic field B_0 threading the disk. When B_0.Omega>0, the horizontal magnetic field is strongly amplified toward the disk interior, leading to a stronger disk wind (by ~50% or less in terms of the wind-driven accretion rate). The enhanced horizontal field also leads to much stronger large-scale Maxwell stress (magnetic braking) that contributes to a considerable fraction of the wind-driven accretion rate. When B_0.Omega<0, the horizontal magnetic field is reduced, leading to a weaker disk wind (by ~20%) and negligible magnetic braking. Moreover, we find that when B_0.Omega>0, the laminar region extends farther to ~15 AU before the magneto-rotational instability sets in, while for B_0.Omega<0, the laminar region extends only to ~3-5 AU for a typical PPD accretion rates. Scaling relations for the wind properties, especially the wind-driven accretion rate, are provided for aligned and anti-aligned field geometries. Issues with the symmetry of the wind solutions and grain abundance are also discussed.
It has recently been shown that the inner region of protoplanetary disks (PPDs) is governed by wind-driven accretion, and the resulting accretion flow showing complex vertical profiles. Such complex flow structures are further enhanced due to the Hall effect, especially when the background magnetic field is aligned with disk rotation. We investigate how such flow structures impact global dust transport via Monte-Carlo simulations, focusing on two scenarios. In the first scenario, the toroidal magnetic field is maximized in the miplane, leading to accretion and decretion flows above and below. In the second scenario, the toroidal field changes sign across the midplane, leading to an accretion flow at the disk midplane, with decretion flows above and below. We find that in both cases, the contribution from additional gas flows can still be accurately incorporated into the advection-diffusion framework for vertically-integrated dust transport, with enhanced dust radial diffusion up to an effective $alpha^{rm eff}sim10^{-2}$ for strongly coupled dust, even when background turbulence is weak $alpha<10^{-4}$. Dust radial drift is also modestly enhanced in the second scenario. We provide a general analytical theory that accurately reproduces our simulation results, thus establishing a framework to model global dust transport that realistically incorporates vertical gas flow structures. We also note that the theory is equally applicable to the transport of chemical species.
Aims: We explore the long-term evolution of young protoplanetary disks with different approaches to computing the thermal structure determined by various cooling and heating processes in the disk and its surroundings. Methods: Numerical hydrodynamics simulations in the thin-disk limit were complemented with three thermal evolution schemes: a simplified $beta$-cooling approach with and without irradiation, in which the rate of disk cooling is proportional to the local dynamical time, a fiducial model with equal dust and gas temperatures calculated taking viscous heating, irradiation, and radiative cooling into account, and also a more sophisticated approach allowing decoupled dust and gas temperatures. Results: We found that the gas temperature may significantly exceed that of dust in the outer regions of young disks thanks to additional compressional heating caused by the infalling envelope material in the early stages of disk evolution and slow collisional exchange of energy between gas and dust in low-density disk regions. The outer envelope however shows an inverse trend with the gas temperatures dropping below that of dust. The global disk evolution is only weakly sensitive to temperature decoupling. Nevertheless, separate dust and gas temperatures may affect the chemical composition, dust evolution, and disk mass estimates. Constant-$beta$ models without stellar and background irradiation fail to reproduce the disk evolution with more sophisticated thermal schemes because of intrinsically variable nature of the $beta$-parameter. Constant-$beta$ models with irradiation can better match the dynamical and thermal evolution, but the agreement is still incomplete. Conclusions: Models allowing separate dust and gas temperatures are needed when emphasis is placed on the chemical or dust evolution in protoplanetary disks, particularly in sub-solar metallicity environments.
The global evolution of protoplanetary disks (PPDs) has recently been shown to be largely controlled by the amount of poloidal magnetic flux threading the disk, which is further controlled by the poorly understood process of magnetic flux transport. In weakly ionized gas as in PPDs, magnetic flux is largely frozen in the electron fluid, except when resistivity is large. When the disk is largely laminar, we show that the relative drift between the electrons and ions (the Hall-drift), and the ions and neutral fluids (ambipolar-drift) can play a dominant role on the transport of magnetic flux. Using two-dimensional simulations that incorporate the Hall effect and ambipolar diffusion (AD) with prescribed diffusivities, we show that when large-scale poloidal field is aligned with disk rotation, the Hall effect rapidly drags magnetic flux inward at the midplane region, while it slowly pushes flux outward above/below the midplane. This leads to a highly radially elongated field configuration as a global manifestation of the Hall-shear instability. This field configuration further promotes rapid outward flux transport by AD at the midplane, leading to instability saturation. In quasi-steady state, magnetic flux is transported outward at approximately the same rate at all heights, and the rate is comparable to the Hall-free case. For anti-aligned field polarity, the Hall effect consistently transports magnetic flux outward, leading to a largely vertical field configuration in the midplane region. The field lines in the upper layer first bend radially inward and then outward to launch a disk wind. Overall, the net rate of outward flux transport is about twice faster than the aligned case. In addition, the rate of flux transport increases with increasing disk magnetization. The absolute rate of transport is sensitive to disk microphysics which remains to be explored in future studies.