No Arabic abstract
Photoevaporation and magnetically driven winds are two independent mechanisms to remove mass from protoplanetary disks. In addition to accretion, the effect of these two principles acting concurrently could be significant and the transition between those two has not been extensively studied and quantified in the literature yet. In order to contribute to the understanding of disk winds, we present the phenomena emerging in the framework of two-dimensional axisymmetric, non-ideal magnetohydrodynamic simulations including EUV-/ X-ray driven photoevaporation. Of particular interest are the examination of the transition region between photoevaporation and magnetically driven wind, the possibility of emerging magneto-centrifugal wind effects, as well as the morphology of the wind itself depending on the strength of the magnetic field. We use the PLUTO code in a 2.5D axisymmetric configuration with additional treatment of EUV-/ X-ray heating and dynamic ohmic diffusion based on a semi-analytical chemical model. We identify the transition between both outflow types to occur for values of the initial plasma beta $beta geq 10^7$, while magnetically driven winds generally outperform photoevaporation for stronger fields. In our simulations we observe irregular and asymmetric outflows for stronger magnetic fields. In the weak field regime the photoevaporation rates are slightly lowered by perturbations of the gas density in the inner regions of the disk. Overall, our results predict a wind with a lever arm smaller than 1.5, consistent with a hot magneto-thermal wind. Stronger accretion flows are present for values of $beta < 10^7$.
We discovered a new growth mode of dust grains to km-sized bodies in protoplanetary disks that evolve by viscous accretion and magnetically driven disk winds (MDWs). We solved an approximate coagulation equation of dust grains with time-evolving disks that consist of both gas and solid components by a one-dimensional model. With the grain growth, all solid particles initially drift inward toward the central star by the gas drag force. However, the radial profile of gas pressure, $P$, is modified by the MDW that disperses the gas in an inside-out manner. Consequently, a local concentration of solid particles is created by the converging radial flux of drifting dust grains at the location with the convex upward profile of $P$. When the dimensionless stopping time, ${rm St}$, there exceeds unity, the solid particles spontaneously reach the growth dominated state because of the positive feedback between the suppressed radial drift and the enhanced accumulation of dust particles that drift from the outer part. Once the solid particles are in the drift limited state, the above-mentioned condition of ${rm St} gtrsim 1$ for the dust growth is equivalent with begin{equation} Sigma_{rm d}/Sigma_{rm g}gtrsim eta, onumber end{equation} where $Sigma_{rm d}/Sigma_{rm g}$ is the dust-to-gas surface-density ratio and $eta$ is dimensionless radial pressure-gradient force. As a consequence of the successful growth of dust grains, a ring-like structure containing planetesimal-sized bodies is formed at the inner part of the protoplanetary disks. Such a ring-shaped concentration of planetesimals is expected to play a vital role in the subsequent planet formation.
We present a novel mechanism for the outward transport of crystalline dust particles: the outward radial drift of pebbles. The dust ring structure is frequently observed in protoplanetary disks. One of the plausible mechanisms of the formation of dust rings is the accumulation of pebbles around the pressure maximum, which is formed by the mass loss due to magnetically driven disk winds. In evolving protoplanetary disks due to magnetically driven disk winds, dust particles can migrate outwardly from the crystallization front to the pressure maximum by radial drift. We found that the outward radial drift process can transport crystalline dust particles efficiently when the radial drift timescale is shorter than the advection timescale. Our model predicts that the crystallinity of silicate dust particles could be as high as 100% inside the dust ring position.
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.
Recent multi-wavelength observations suggest that inner parts of protoplanetary disks (PPDs) have shorter lifetimes for heavier host stars. Since PPDs around high-mass stars are irradiated by strong ultra-violet radiation, photoevaporation may provide an explanation for the observed trend. We perform radiation hydrodynamics simulations of photoevaporation of PPDs for a wide range of host star mass of $M_* =0.5$-$7.0 M_{odot}$. We derive disk mass-loss rate $dot{M}$, which has strong stellar dependence as $dot{M} approx 7.30times10^{-9}(M_{*}/M_{odot})^{2}M_{odot}rm{yr}^{-1}$. The absolute value of $dot{M}$ scales with the adopted far-ultraviolet and X-ray luminosities. We derive the surface mass-loss rates and provide polynomial function fits to them. We also develop a semi-analytic model that well reproduces the derived mass-loss rates. The estimated inner disk lifetime decreases as the host star mass increases, in agreement with the observational trend. We thus argue that photoevaporation is a major physical mechanism for PPD dispersal for a wide range of the stellar mass and can account for the observed stellar mass dependence of the inner disk lifetime.
The gas dynamics of weakly ionized protoplanetary disks (PPDs) is largely governed by the coupling between gas and magnetic fields, described by three non-ideal magnetohydrodynamical (MHD) effects (Ohmic, Hall, ambipolar). Previous local simulations incorporating these processes have revealed that the inner regions of PPDs are largely laminar accompanied by wind-driven accretion. We conduct 2D axisymmetric, fully global MHD simulations of these regions ($sim1-20$ AU), taking into account all non-ideal MHD effects, with tabulated diffusion coefficients and approximate treatment of external ionization and heating. With net vertical field aligned with disk rotation, the Hall-shear instability strongly amplifies horizontal magnetic field, making the overall dynamics dependent on initial field configuration. Following disk formation, the disk likely relaxes into an inner zone characterized by asymmetric field configuration across the midplane that smoothly transitions to a more symmetric outer zone. Angular momentum transport is driven by both MHD winds and laminar Maxwell stress, with both accretion and decretion flows present at different heights, and modestly asymmetric winds from the two disk sides. With anti-aligned field polarity, weakly magnetized disks settle into an asymmetric field configuration with supersonic accretion flow concentrated at one side of disk surface, and highly asymmetric winds between the two disk sides. In all cases, the wind is magneto-thermal in nature characterized by mass loss rate exceeding the accretion rate. More strongly magnetized disks give more symmetric field configuration and flow structures. Deeper far-UV penetration leads to stronger and less stable outflows. Implications for observations and planet formation are also discussed.