Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userHall-effect Controlled Gas Dynamics in Protoplanetary Disks: II. Full 3D Simulations toward the Outer Disk
134
0
0.0
(
0
)
Authors
Xue-Ning Bai
Ask ChatGPT about the research
No Arabic abstract
We perform 3D stratified shearing-box MHD simulations on the gas dynamics of protoplanetary disks threaded by net vertical magnetic field Bz. All three non-ideal MHD effects, Ohmic resistivity, the Hall effect and ambipolar diffusion are included in a self-consistent manner based on equilibrium chemistry. We focus on regions toward outer disk radii, from 5-60AU, where Ohmic resistivity tends to become negligible, ambipolar diffusion dominates over an extended region across disk height, and the Hall effect largely controls the dynamics near the disk midplane. We find that around R=5AU, the system launches a laminar/weakly turbulent magnetocentrifugal wind when the net vertical field Bz is not too weak, as expected. Moreover, the wind is able to achieve and maintain a configuration with reflection symmetry at disk midplane. The case with anti-aligned field polarity (Omega. Bz<0) is more susceptible to the MRI when Bz drops, leading to an outflow oscillating in radial directions and very inefficient angular momentum transport. At the outer disk around and beyond R=30AU, the system shows vigorous MRI turbulence in the surface layer due to far-UV ionization, which efficiently drives disk accretion. The Hall effect affects the stability of the midplane region to the MRI, leading to strong/weak Maxwell stress for aligned/anti-aligned field polarities. Nevertheless, the midplane region is only very weakly turbulent. Overall, the basic picture is analogous to the conventional layered accretion scenario applied to the outer disk. In addition, we find that the vertical magnetic flux is strongly concentrated into thin, azimuthally extended shells in most of our simulations beyond 15AU. This is a generic phenomenon unrelated to the Hall effect, and leads to enhanced zonal flow. Observational and theoretical implications, as well as future prospects are briefly discussed.
rate research
Read More
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.
Theoretical models of the ionization state in protoplanetary disks suggest the existence of large areas with low ionization and weak coupling between the gas and magnetic fields. In this regime hydrodynamical instabilities may become important. In this work we investigate the gas and dust structure and dynamics for a typical T Tauri system under the influence of the vertical shear instability (VSI). We use global 3D radiation hydrodynamics simulations covering all $360^circ$ of azimuth with embedded particles of 0.1 and 1mm size, evolved for 400 orbits. Stellar irradiation heating is included with opacities for 0.1- to 10-$mu$m-sized dust. Saturated VSI turbulence produces a stress-to-pressure ratio of $alpha simeq 10^{-4}$. The value of $alpha$ is lowest within 30~au of the star, where thermal relaxation is slower relative to the orbital period and approaches the rate below which VSI is cut off. The rise in $alpha$ from 20 to 30~au causes a dip in the surface density near 35~au, leading to Rossby wave instability and the generation of a stationary, long-lived vortex spanning about 4~au in radius and 40~au in azimuth. Our results confirm previous findings that mm size grains are strongly vertically mixed by the VSI. The scale height aspect ratio for 1mm grains is determined to be 0.037, much higher than the value $H/r=0.007$ obtained from millimeter-wave observations of the HL~Tau system. The measured aspect ratio is better fit by non-ideal MHD models. In our VSI turbulence model, the mm grains drift radially inwards and many are trapped and concentrated inside the vortex. The turbulence induces a velocity dispersion of $sim 12$~m/s for the mm grains, indicating that grain-grain collisions could lead to fragmentation.
We present ALMA Band 6 observations of a complete sample of protoplanetary disks in the young (1-3 Myr) Lupus star-forming region, covering the 1.33 mm continuum and the 12CO, 13CO, and C18O J=2-1 lines. The spatial resolution is 0.25 arcsec with a medium 3-sigma continuum sensitivity of 0.30 mJy, corresponding to M_dust ~ 0.2 M_earth. We apply Keplerian masking to enhance the signal-to-noise ratios of our 12CO zero-moment maps, enabling measurements of gas disk radii for 22 Lupus disks; we find that gas disks are universally larger than mm dust disks by a factor of two on average, likely due to a combination of the optically thick gas emission as well as the growth and inward drift of the dust. Using the gas disk radii, we calculate the dimensionless viscosity parameter, alpha_visc, finding a broad distribution and no correlations with other disk or stellar parameters, suggesting that viscous processes have not yet established quasi-steady states in Lupus disks. By combining our 1.33 mm continuum fluxes with our previous 890 micron continuum observations, we also calculate the mm spectral index, alpha_mm, for 70 Lupus disks; we find an anti-correlation between alpha_mm and mm flux for low-mass disks (M_dust < 5), followed by a flattening as disks approach alpha_mm = 2, which could indicate faster grain growth in higher-mass disks, but may also reflect their larger optically thick components. In sum, this work demonstrates the continuous stream of new insights into disk evolution and planet formation that can be gleaned from unbiased ALMA disk surveys.
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.
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.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
