Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userGlobal simulations of protoplanetary disks with net magnetic flux: I. Non-ideal MHD case
81
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
The planet-forming region of protoplanetary disks is cold, dense, and therefore weakly ionized. For this reason, magnetohydrodynamic (MHD) turbulence is thought to be mostly absent, and another mechanism has to be found to explain gas accretion. It has been proposed that magnetized winds, launched from the ionized disk surface, could drive accretion in the presence of a large-scale magnetic field. The efficiency and the impact of these surface winds on the disk structure is still highly uncertain. We present the first global simulations of a weakly ionized disk that exhibits large-scale magnetized winds. We also study the impact of self-organization, which was previously demonstrated only in non-stratified models. We perform numerical simulations of stratified disks with the PLUTO code. We compute the ionization fraction dynamically, and account for all three non-ideal MHD effects: ohmic and ambipolar diffusions, and the Hall drift. Simplified heating and cooling due to non-thermal radiation is also taken into account in the disk atmosphere. We find that disks can be accreting or not, depending on the configuration of the large-scale magnetic field. Magnetothermal winds, driven both by magnetic acceleration and heating of the atmosphere, are obtained in the accreting case. In some cases, these winds are asymmetric, ejecting predominantly on one side of the disk. The wind mass loss rate depends primarily on the average ratio of magnetic to thermal pressure in the disk midplane. The non-accreting case is characterized by a meridional circulation, with accretion layers at the disk surface and decretion in the midplane. Finally, we observe self-organization, resulting in axisymmetric rings of density and associated pressure bumps. The underlying mechanism and its impact on observable structures are discussed.
rate research
Read More
It has recently been established that the evolution of protoplanetary disks is primarily driven by magnetized disk winds, requiring large-scale magnetic flux threading the disks. The size of such disks is expected to shrink in time, as opposed to the conventional scenario of viscous expansion. We present the first global 2D non-ideal magnetohydrodynamic (MHD) simulations of protoplanetary disks that are truncated in the outer radius, aiming to understand the interaction of the disk with the interstellar environment, as well as global evolution of the disk and magnetic flux. We find that as the system relaxes, poloidal magnetic field threading the disk beyond the truncation radius collapses towards the midplane, leading to rapid reconnection. This process removes a substantial amount of magnetic flux from the system, and forms closed poloidal magnetic flux loops encircling the outer disk in quasi-steady-state. These magnetic flux loops can drive expansion beyond truncation radius, corresponding to substantial mass loss through magnetized disk outflow beyond truncation radius analogous to a combination of viscous spreading and external photoevaporation. The magnetic flux loops gradually shrink over time whose rates depend on level of disk magnetization and external environments, which eventually governs the long-term disk evolution.
We study dust transport in turbulent protoplanetary disks using three-dimensional global unstratified magnetohydrodynamic (MHD) simulations including Lagrangian dust particles. The turbulence is driven by the magnetorotational instability (MRI) with either ideal or non-ideal MHD that includes ambipolar diffusion (AD). In ideal MHD simulations, the surface density evolution (except for dust that drifts fastest), turbulent diffusion, and vertical scale height of dust can all be reproduced by simple one-dimensoinal and/or analytical models. However, in AD dominated simulations which simulate protoplanetary disks beyond 10s of AU, the vertical scale height of dust is larger than previously predicted. To understand this anomaly in more detail, we carry out both unstratified and stratified local shearing box simulations with Lagrangian particles, and find that turbulence in AD dominated disks has very different properties (e.g., temporal autocorrelation functions and power spectra) than turbulence in ideal MHD disks, which leads to quite different particle diffusion efficiency. For example, MRI turbulence with AD has a longer correlation time for the vertical velocity, which causes significant vertical particle diffusion and large dust scale height. In ideal MHD the Schmidt numbers ($Sc$) for radial and vertical turbulent diffusion are $Sc_{r}sim 1$ and $Sc_{z}gtrsim 3$, but in the AD dominated regime both $Sc_{r}$ and $Sc_{z}$ are $lesssim 1$. Particle concentration in pressure bumps induced by MRI turbulence has also been studied. Since non-ideal MHD effects dominate most regions in protoplanetary disks, our study suggests that modeling dust transport in turbulence driven by MRI with non-ideal MHD effects is important for understanding dust transport in realistic protoplanetary disks.
The structure and evolution of protoplanetary disks (PPDs) are largely governed by disk angular momentum transport, mediated by magnetic fields. In the most observable outer disk, PPD gas dynamics is primarily controlled by ambipolar diffusion as the dominant non-ideal magnetohydrodynamic (MHD) effect. In this work, we study the gas dynamics in outer PPDs by conducting a set of global 3D non-ideal MHD simulations with ambipolar diffusion and net poloidal magnetic flux, using the Athena++ MHD code, with resolution comparable to local simulations. Our simulations demonstrate the co-existence of magnetized disk wind and turbulence driven by the magneto-rotational instability (MRI). While MHD wind dominates disk angular momentum transport, the MRI turbulence also contributes significantly. We observe that magnetic flux spontaneously concentrate into axisymmetric flux sheets, leading to radial variations in turbulence levels, stresses, and accretion rates. Annular substructures arise as a natural consequence of magnetic flux concentration. The flux concentration phenomena show diverse properties with different levels of disk magnetization and ambipolar diffusion. The disk generally loses magnetic flux over time, though flux sheets could prevent the leak of magnetic flux in some cases. Our results demonstrate the ubiquity of disk annular substructures in weakly MRI turbulent outer PPDs, and imply a stochastic nature of disk evolution.
Outflows driven by large-scale magnetic fields likely play an important role in the evolution and dispersal of protoplanetary disks, and in setting the conditions for planet formation. We extend our 2-D axisymmetric non-ideal MHD model of these outflows by incorporating radiative transfer and simplified thermochemistry, with the twin aims of exploring how heating influences wind launching, and illustrating how such models can be tested through observations of diagnostic spectral lines. Our model disks launch magnetocentrifugal outflows primarily through magnetic tension forces, so the mass-loss rate increases only moderately when thermochemical effects are switched on. For typical field strengths, thermochemical and irradiation heating are more important than magnetic dissipation. We furthermore find that the entrained vertical magnetic flux diffuses out of the disk on secular timescales as a result of non-ideal MHD. Through post-processing line radiative transfer, we demonstrate that spectral line intensities and moment-1 maps of atomic oxygen, the HCN molecule, and other species show potentially observable differences between a model with a magnetically driven outflow and one with a weaker, photoevaporative outflow. In particular, the line shapes and velocity asymmetries in the moment-1 maps could enable the identification of outflows emanating from the disk surface.
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.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
